Research resource: transcriptional response to glucocorticoids in childhood acute lymphoblastic leukemia

The XPO1 Inhibitor KPT-8602 Synergizes with Dexamethasone in Acute Lymphoblastic Leukemia

PURPOSE

KPT-8602 (Eltanexor) is a second-generation exportin-1 (XPO1) inhibitor with potent activity against acute lymphoblastic leukemia (ALL) in preclinical models and with minimal effects on normal cells. In this study, we evaluated whether KPT-8602 would synergize with dexamethasone, vincristine, or doxorubicin, three drugs currently used for the treatment of ALL.

EXPERIMENTAL DESIGN

First, we searched for the most synergistic combination of KPT-8602 with dexamethasone, vincristine, or doxorubicin in vitro in both B-ALL and T-ALL cell lines using proliferation and apoptosis as a readout. Next, we validated this synergistic effect by treatment of clinically relevant B- and T-ALL patient-derived xenograft models in vivo. Finally, we performed RNA-sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) to determine the mechanism of synergy.

RESULTS

KPT-8602 showed strong synergism with dexamethasone on human B-ALL and T-ALL cell lines as well as in vivo in three patient-derived ALL xenografts. Compared with single-drug treatment, the drug combination caused increased apoptosis and led to histone depletion. Mechanistically, integration of ChIP-seq and RNA-seq data revealed that addition of KPT-8602 to dexamethasone enhanced the activity of the glucocorticoid receptor (NR3C1) and led to increased inhibition of E2F-mediated transcription. We observed strong inhibition of E2F target genes related to cell cycle, DNA replication, and transcriptional regulation.

CONCLUSIONS

Our preclinical study demonstrates that KPT-8602 enhances the effects of dexamethasone to inhibit B-ALL and T-ALL cells via NR3C1- and E2F-mediated transcriptional complexes, allowing to achieve increased dexamethasone effects for patients.

3D-cultivation of NSCLC cell lines induce gene expression alterations of key cancer-associated pathways and mimic in-vivo conditions

This work evaluated gene expression differences between a hanging-drop 3D NSCLC model and 2D cell cultures and their in-vivo relevance by comparison to patient-derived data from The Cancer Genome Atlas.

Gene expression of 2D and 3D cultures for Colo699 and A549 were assessed using Affymetrix HuGene 1.0 ST gene chips. Biostatistical analyses tested for reproducibility, comparability and significant differences in gene expression profiles between cell lines, experiments and culture methods.

The analyses revealed a high interassay correlation within specific culture systems proving a high validity. 979 genes were altered in A549 and 1106 in Colo699 cells due to 3D cultivation. The overlap of changed genes between the cell lines was small (149), but the involved pathways in the reactome and GO- analyses showed a high overlap with DNA methylation, cell cycle, SIRT1, PKN1 pathway, DNA repair and oxidative stress as well known cancer-associated representatives. Additional specific GSEA-analyses revealed changes in immunologic and endothelial cell proliferation pathways, whereas hypoxic, EMT and angiogenic pathways were downregulated.

Gene enrichment analyses showed 3D-induced gene up-regulations in the cell lines 38 to be represented in in-vivo samples of NSCLC patients using data of The Cancer Genome Atlas.

Thus, our 3D NSCLC model might provide a tool for early drug development and investigation of microenvironment-associated mechanisms. However, this work also highlights the need for further individualization and model adaption to address remaining challenges.

Dexamethasone-induced cell death is restricted to specific molecular subgroups of multiple myeloma

Due to its cytotoxic effect in lymphoid cells, dexamethasone is widely used in the treatment of multiple myeloma (MM). However, only a subset of myeloma patients responds to high-dose dexamethasone. Despite the undeniable anti-myeloma benefits of dexamethasone, significant adverse effects have been reported. We re-evaluate the anti-tumor effect of dexamethasone according to the molecular heterogeneity of MM. We demonstrated that the pro-death effect of dexamethasone is related to the genetic heterogeneity of MM because sensitive cell lines were restricted to MAF and MMSET signature subgroups, whereas all CCND1 cell lines (n = 10) were resistant to dexamethasone. We demonstrated that the glucocorticoid receptor expression was an important limiting factor for dexamethasone-induced cell death and we found a correlation between glucocorticoid receptor levels and the induction of glucocorticoid-induced leucine zipper (GILZ) under dexamethasone treatment. By silencing GILZ, we next demonstrated that GILZ is necessary for Dex induced apoptosis while triggering an imbalance between anti- and pro-apoptotic Bcl-2 proteins. Finally, the heterogeneity of the dexamethasone response was further confirmed in vivo using myeloma xenograft models. Our findings suggested that the effect of dexamethasone should be re-evaluated within molecular subgroups of myeloma patients to improve its efficacy and reduce its adverse effects.

New insights into the genetics of glioblastoma multiforme by familial exome sequencing

Glioblastoma multiforme (GBM) is the most aggressive and malignant subtype of human brain tumors. While a family clustering of GBM has long been acknowledged, relevant hereditary factors still remained elusive. Exome sequencing of families offers the option to discover respective genetic factors.We sequenced blood samples of one of the rare affected families: while both parents were healthy, both children were diagnosed with GBM. We report 85 homozygous non-synonymous single nucleotide variations (SNVs) in both siblings that were heterozygous in the parents. Beyond known key players for GBM such as ERBB2, PMS2, or CHI3L1, we identified over 50 genes that have not been associated to GBM so far. We also discovered three accumulative effects potentially adding to the tumorigenesis in the siblings: a clustering of multiple variants in single genes (e.g., PTPRB, CROCC), the aggregation of affected genes on specific molecular pathways (e.g., Focal adhesion or ECM receptor interaction) and genomic proximity (e.g., chr22.q12.2, chr1.p36.33). We found a striking accumulation of SNVs in specific genes for the daughter, who developed not only a GBM at the age of 12 years but was subsequently diagnosed with a pilocytic astrocytoma, a common acute lymphatic leukemia and a diffuse pontine glioma.The reported variants underline the relevance of genetic predisposition and cancer development in this family and demonstrate that GBM has a complex and heterogeneous genetic background. Sequencing of other affected families will help to further narrow down the driving genetic causes for this disease.

Dintor: functional annotation of genomic and proteomic data

Background

During the last decade, a great number of extremely valuable large-scale genomics and proteomics datasets have become available to the research community. In addition, dropping costs for conducting high-throughput sequencing experiments and the option to outsource them considerably contribute to an increasing number of researchers becoming active in this field. Even though various computational approaches have been developed to analyze these data, it is still a laborious task involving prudent integration of many heterogeneous and frequently updated data sources, creating a barrier for interested scientists to accomplish their own analysis.

Results

We have implemented Dintor, a data integration framework that provides a set of over 30 tools to assist researchers in the exploration of genomics and proteomics datasets. Each of the tools solves a particular task and several tools can be combined into data processing pipelines. Dintor covers a wide range of frequently required functionalities, from gene identifier conversions and orthology mappings to functional annotation of proteins and genetic variants up to candidate gene prioritization and Gene Ontology-based gene set enrichment analysis. Since the tools operate on constantly changing datasets, we provide a mechanism to unambiguously link tools with different versions of archived datasets, which guarantees reproducible results for future tool invocations. We demonstrate a selection of Dintor’s capabilities by analyzing datasets from four representative publications. The open source software can be downloaded and installed on a local Unix machine. For reasons of data privacy it can be configured to retrieve local data only. In addition, the Dintor tools are available on our public Galaxy web service at http://dintor.eurac.edu.

Conclusions

Dintor is a computational annotation framework for the analysis of genomic and proteomic datasets, providing a rich set of tools that cover the most frequently encountered tasks. A major advantage is its capability to consistently handle multiple versions of tool-associated datasets, supporting the researcher in delivering reproducible results.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2279-5) contains supplementary material, which is available to authorized users.

Quantitative proteomic analysis reveals maturation as a mechanism underlying glucocorticoid resistance in B lineage ALL and re-sensitization by JNK inhibition

Glucocorticoid (GC) resistance is a continuing clinical problem in childhood acute lymphoblastic leukaemia (ALL) but the underlying mechanisms remain unclear. A proteomic approach was used to compare profiles of the B-lineage ALL GC-sensitive cell line, PreB 697, and its GC-resistant sub-line, R3F9, pre- and post-dexamethasone exposure. PAX5, a transcription factor critical to B-cell development was differentially regulated in the PreB 697 compared to the R3F9 cell line in response to GC. PAX5 basal protein expression was less in R3F9 compared to its GC-sensitive parent and confirmed to be lower in other GC-resistant sub-lines of Pre B 697 and was associated with a decreased expression of the PAX5 transcriptional target, CD19. Gene set enrichment analysis showed that increasing GC-resistance was associated with differentiation from preB-II to an immature B-lymphocyte stage. GC-resistant sub-lines were shown to have higher levels of phosphorylated JNK compared to the parent line and JNK inhibition caused re-sensitization to GC. Exploiting this maturation may be key to overcoming GC resistance and targeting signalling pathways linked to the maturation state, such as JNK, may be a novel approach.

Allergenic Can f 1 and its human homologue Lcn-1 direct dendritic cells to induce divergent immune responses

Why and when the immune system skews to Th2 mediated allergic immune responses is still poorly characterized. With two homologous lipocalins, the major respiratory dog allergen Can f 1 and the human endogenous, non-allergenic Lipocalin-1, we investigated their impact on human monocyte-derived dendritic cells (DC). The two lipocalins had differential effects on DC according to their allergenic potential. Compared to Lipocalin-1, Can f 1 persistently induced lower levels of the Th1 skewing maturation marker expression, tryptophan breakdown and interleukin (IL)-12 production in DC. As a consequence, T cells stimulated by DC treated with Can f 1 produced more of the Th2 signature cytokine IL-13 and lower levels of the Th1 signature cytokine interferon-γ than T cells stimulated by Lipocalin-1 treated DC. These data were partially verified by a second pair of homologous lipocalins, the cat allergen Fel d 4 and its putative human homologue major urinary protein. Our data indicate that the crosstalk of DC with lipocalins alone has the potential to direct the type of immune response to these particular antigens. A global gene expression analysis further supported these results and indicated significant differences in intracellular trafficking, sorting and antigen presentation pathways when comparing Can f 1 and Lipocalin-1 stimulated DC. With this study we contribute to a better understanding of the induction phase of a Th2 immune response.

miR-22 and miR-29a Are Members of the Androgen Receptor Cistrome Modulating LAMC1 and Mcl-1 in Prostate Cancer

The normal prostate as well as early stages and advanced prostate cancer (PCa) require a functional androgen receptor (AR) for growth and survival. The recent discovery of microRNAs (miRNAs) as novel effector molecules of AR disclosed the existence of an intricate network between AR, miRNAs and downstream target genes. In this study DUCaP cells, characterized by high content of wild-type AR and robust AR transcriptional activity, were chosen as the main experimental model. By integrative analysis of chromatin immunoprecipitation-sequencing (ChIP-seq) and microarray expression profiling data, miRNAs putatively bound and significantly regulated by AR were identified. A direct AR regulation of miR-22, miR-29a, and miR-17-92 cluster along with their host genes was confirmed. Interestingly, endogenous levels of miR-22 and miR-29a were found to be reduced in PCa cells expressing AR. In primary tumor samples, miR-22 and miR-29a were less abundant in the cancerous tissue compared with the benign counterpart. This specific expression pattern was associated with a differential DNA methylation of the genomic AR binding sites. The identification of laminin gamma 1 (LAMC1) and myeloid cell leukemia 1 (MCL1) as direct targets of miR-22 and miR-29a, respectively, suggested a tumor-suppressive role of these miRNAs. Indeed, transfection of miRNA mimics in PCa cells induced apoptosis and diminished cell migration and viability. Collectively, these data provide additional information regarding the complex regulatory machinery that guides miRNAs activity in PCa, highlighting an important contribution of miRNAs in the AR signaling.

The tumor-modulatory effects of Caspase-2 and Pidd1 do not require the scaffold protein Raidd

The receptor-interacting protein-associated ICH-1/CED-3 homologous protein with a death domain (RAIDD/CRADD) functions as a dual adaptor and is a constituent of different multi-protein complexes implicated in the regulation of inflammation and cell death. Within the PIDDosome complex, RAIDD connects the cell death-related protease, Caspase-2, with the p53-induced protein with a death domain 1 (PIDD1). As such, RAIDD has been implicated in DNA-damage-induced apoptosis as well as in tumorigenesis. As loss of Caspase-2 leads to an acceleration of tumor onset in the Eμ-Myc mouse lymphoma model, whereas loss of Pidd1 actually delays onset of this disease, we set out to interrogate the role of Raidd in cancer in more detail. Our data obtained analyzing Eμ-Myc/Raidd^−/− mice indicate that Raidd is unable to protect from c-Myc-driven lymphomagenesis. Similarly, we failed to observe a modulatory effect of Raidd deficiency on DNA-damage-driven cancer. The role of Caspase-2 as a tumor suppressor and that of Pidd1 as a tumor promoter can therefore be uncoupled from their ability to interact with the Raidd scaffold, pointing toward the existence of alternative signaling modules engaging these two proteins in this context.

The synthetic glucocorticoids prednisolone and dexamethasone regulate the same genes in acute lymphoblastic leukemia cells

Background

Glucocorticoids (GCs) cause apoptosis in malignant cells of lymphoid lineage by transcriptionally regulating a plethora of genes. As a result, GCs are included in almost all treatment protocols for lymphoid malignancies, particularly childhood acute lymphoblastic leukemia (chALL). The most commonly used synthetic GCs in the clinical setting are prednisolone and dexamethasone. While the latter has a higher activity and more effectively reduces the tumor load in patients, it is also accompanied by more serious adverse effects than the former. Whether this difference might be explained by regulation of different genes by the two GCs has never been addressed.

Results

Using a recently developed GC bioassay based on a GC-responsive reporter construct in human Jurkat T-ALL cells, we found ~7-fold higher biological activity with dexamethasone than prednisolone. Similarly, 1.0e-7 M dexamethasone and 7.0e-7 M prednisolone triggered similar cell death rates in CCRF-CEM-C7H2 T-chALL cells after 72 hours of treatment. Using microarray-based whole genome expression profiling and a variety of statistical and other approaches, we compared the transcriptional response of chALL cells to 6 hour exposure to both synthetic GCs at the above concentrations. Our experiments did not detect any gene whose regulation by dexamethasone differed significantly from that by prednisolone.

Conclusions

Our findings suggest that the reported differences in treatment efficacy and cytotoxicity of dexamethasone and prednisolone are not caused by inherent differences of the 2 drugs to regulate the expression of certain genes, but rather result either from applying them in biologically in-equivalent concentrations and/or from differences in their pharmacokinetics and - dynamics resulting in different bioactivities in tumor cells and normal tissues.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-662) contains supplementary material, which is available to authorized users.

SOCS2 correlates with malignancy and exerts growth-promoting effects in prostate cancer

Deregulation of cytokine and growth factor signaling due to an altered expression of endogenous regulators is well recognized in prostate cancer (PCa) and other cancers. Suppressor of cytokine signaling 2 (SOCS2) is a key regulator of the GH, IGF, and prolactin signaling pathways that have been implicated in carcinogenesis. In this study, we evaluated the expression patterns and functional significance of SOCS2 in PCa. Protein expression analysis employing tissue microarrays from two independent patient cohorts revealed a significantly enhanced expression in tumor tissue compared with benign tissue as well as association with Gleason score and disease progression. In vitro and in vivo assays uncovered the involvement of SOCS2 in the regulation of cell growth and apoptosis. Functionally, SOCS2 knockdown inhibited PCa cell proliferation and xenograft growth in a CAM assay. Decreased cell growth after SOCS2 downregulation was associated with cell-cycle arrest and apoptosis. In addition, we proved that SOCS2 expression is significantly elevated upon androgenic stimulation in androgen receptor (AR)-positive cell lines, providing a possible mechanistic explanation for high SOCS2 levels in PCa tissue. Consequently, SOCS2 expression correlated with AR expression in the malignant tissue of patients. On the whole, our study linked increased SOCS2 expression in PCa with a pro-proliferative role in vitro and in vivo.

Translational profiling in childhood acute lymphoblastic leukemia: no evidence for glucocorticoid regulation of mRNA translation

Background

Glucocorticoids (GCs) are natural stress induced steroid hormones causing cell cycle arrest and cell death in lymphoid tissues. Therefore they are the central component in the treatment of lymphoid malignancies, in particular childhood acute lymphoblastic leukemia (chALL). GCs act mainly via regulating gene transcription, which has been intensively studied by us and others. GC control of mRNA translation has also been reported but has never been assessed systematically. In this study we investigate the effect of GCs on mRNA translation on a genome-wide scale.

Results

Childhood T- (CCRF-CEM) and precursor B-ALL (NALM6) cells were exposed to GCs and subjected to “translational profiling”, a technique combining sucrose-gradient fractionation followed by Affymetrix Exon microarray analysis of mRNA from different fractions, to assess the translational efficiency of the expressed genes. Analysis of GC regulation in ribosome-bound fractions versus transcriptional regulation revealed no significant differences, i.e., GC did not entail a significant shift between ribosomal bound and unbound mRNAs.

Conclusions

In the present study we analyzed for the first time possible effects of GC on the translational efficiency of expressed genes in two chALL model systems employing whole genome polysome profiling. Our results did not reveal significant differences in translational efficiency of expressed genes thereby arguing against a potential widespread regulatory effect of GCs on translation at least in the investigated in vitro systems.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-844) contains supplementary material, which is available to authorized users.

EpCAM overexpression prolongs proliferative capacity of primary human breast epithelial cells and supports hyperplastic growth

Introduction

The Epithelial Cell Adhesion Molecule (EpCAM) has been shown to be strongly expressed in human breast cancer and cancer stem cells and its overexpression has been supposed to support tumor progression and metastasis. However, effects of EpCAM overexpression on normal breast epithelial cells have never been studied before. Therefore, we analyzed effects of transient adenoviral overexpression of EpCAM on proliferation, migration and differentiation of primary human mammary epithelial cells (HMECs).

Methods

HMECs were transfected by an adenoviral system for transient overexpression of EpCAM. Thereafter, changes in cell proliferation and migration were studied using a real time measurement system. Target gene expression was evaluated by transcriptome analysis in proliferating and polarized HMEC cultures. A Chicken Chorioallantoic Membrane (CAM) xenograft model was used to study effects on in vivo growth of HMECs.

Results

EpCAM overexpression in HMECs did not significantly alter gene expression profile of proliferating or growth arrested cells. Proliferating HMECs displayed predominantly glycosylated EpCAM isoforms and were inhibited in cell proliferation and migration by upregulation of p27^KIP1 and p53. HMECs with overexpression of EpCAM showed a down regulation of E-cadherin. Moreover, cells were more resistant to TGF-β1 induced growth arrest and maintained longer capacities to proliferate in vitro. EpCAM overexpressing HMECs xenografts in chicken embryos showed hyperplastic growth, lack of lumen formation and increased infiltrates of the chicken leukocytes.

Conclusions

EpCAM revealed oncogenic features in normal human breast cells by inducing resistance to TGF-β1-mediated growth arrest and supporting a cell phenotype with longer proliferative capacities in vitro. EpCAM overexpression resulted in hyperplastic growth in vivo. Thus, we suggest that EpCAM acts as a prosurvival factor counteracting terminal differentiation processes in normal mammary glands.

MicroRNAs and Glucocorticoid-Induced Apoptosis in Lymphoid Malignancies

The initial response of lymphoid malignancies to glucocorticoids (GCs) is a critical parameter predicting successful treatment. Although being known as a strong inducer of apoptosis in lymphoid cells for almost a century, the signaling pathways regulating the susceptibility of the cells to GCs are only partly revealed. There is still a need to develop clinical tests that can predict the outcome of GC therapy. In this paper, I discuss important parameters modulating the pro-apoptotic effects of GCs, with a specific emphasis on the microRNA world comprised of small players with big impacts. The journey through the multifaceted complexity of GC-induced apoptosis brings forth explanations for the differential treatment response and raises potential strategies for overcoming drug resistance.

Identification and characterization of alternative exon usage linked glioblastoma multiforme survival

Background

Alternative exon usage (AEU) is an important component of gene regulation. Exon expression platforms allow the detection of associations between AEU and phenotypes such as cancer. Numerous studies have identified associations between gene expression and the brain cancer glioblastoma multiforme (GBM). The few consistent gene expression biomarkers of GBM that have been reported may be due to the limited consideration of AEU and the analytical approaches used. The objectives of this study were to develop a model that accounts for the variations in expression present between the exons within a gene and to identify AEU biomarkers of GBM survival.

Methods

The expression of exons corresponding to 25,403 genes was related to the survival of 250 individuals diagnosed with GBM in a training data set. Genes exhibiting AEU in the training data set were confirmed in an independent validation data set of 78 patients. A hierarchical mixed model that allows the consideration of covariation between exons within a gene and of the effect of the epidemiological characteristics of the patients was developed to identify associations between exon expression and patient survival. This general model describes all three possible scenarios: multi-exon genes with and without AEU, and single-exon genes.

Results

AEU associated with GBM survival was identified on 2477 genes (P-value < 5.0E-04 or FDR-adjusted P-value < 0.05). G-protein coupled receptor 98 (Gpr98) and epidermal growth factor (Egf) were among the genes exhibiting AEU with 30 and 9 exons associated with GBM survival, respectively. Pathways enriched among the AEU genes included focal adhesion, ECM-receptor interaction, ABC transporters and pathways in cancer. In addition, 24 multi-exon genes without AEU and 8 single-exon genes were associated with GBM survival (FDR-adjusted P-value < 0.05).

Conclusions

The inferred patterns of AEU were consistent with in silico AS models. The hierarchical model used offered a flexible and simple way to interpret and identify associations between survival that accommodates multi-exon genes with or without AEU and single exon genes. Our results indicate that differential expression of AEU could be used as biomarker for GBM and potentially other cancers.