Exon-level expression analyses identify MYCN and NTRK1 as major determinants of alternative exon usage and robustly predict primary neuroblastoma outcome.

BACKGROUND: Using mRNA expression-derived signatures as predictors of individual patient outcome has been a goal ever since the introduction of microarrays. Here, we addressed whether analyses of tumour mRNA at the exon level can improve on the predictive power and classification accuracy of gene-based expression profiles using neuroblastoma as a model. METHODS: In a patient cohort comprising 113 primary neuroblastoma specimens expression profiling using exon-level analyses was performed to define predictive signatures using various machine-learning techniques. Alternative transcript use was calculated from relative exon expression. Validation of alternative transcripts was achieved using qPCR- and cell-based approaches. RESULTS: Both predictors derived from the gene or the exon levels resulted in prediction accuracies >80% for both event-free and overall survival and proved as independent prognostic markers in multivariate analyses. Alternative transcript use was most prominently linked to the amplification status of the MYCN oncogene, expression of the TrkA/NTRK1 neurotrophin receptor and survival. CONCLUSION: As exon level-based prediction yields comparable, but not significantly better, prediction accuracy than gene expression-based predictors, gene-based assays seem to be sufficiently precise for predicting outcome of neuroblastoma patients. However, exon-level analyses provide added knowledge by identifying alternative transcript use, which should deepen the understanding of neuroblastoma biology.

Correlation between porcine PPARGC1A mRNA expression and its downstream target genes in backfat and longissimus dorsi muscle.

Knowledge of in vivo relationship between the coactivator PPARGC1A and its target genes is very limited, especially in the pig. In this study, a real-time PCR experiment was performed on longissimus dorsi muscle (MLD) and backfat with 10 presumed PPARGC1A downstream target genes, involved in energy and fat metabolism, to identify possible relationships with PPARGC1A mRNA expression in vivo in the pig (n = 20). Except for UCP3 and LPL, a very significant difference in expression was found between MLD and backfat for all genes (P < 0.01). Hierarchical cluster analysis and the significant pairing of mRNA expression data between sampling locations suggested a genetic regulation of the expression of several target genes. A positive correlation with PPARGC1A was found for CPT1B, GLUT4, PDK4, and TFAM (P < 0.0001). A negative correlation was found for UCP2, FABP4, LEP (P < 0.0001), and TNF (P = 0.0071). No significant correlation was detected for UCP3 and LPL. This study provides evidence for a clear difference in mRNA expression of crucial genes in fat and energy metabolism between 2 important tissues. Our data suggest a clear impact of PPARGC1A on energy and lipid metabolism in vivo in the pig, through several of these downstream target genes.

Emerging patterns of cryptic chromosomal imbalance in patients with idiopathic mental retardation and multiple congenital anomalies: a new series of 140 patients and review of published reports.

BACKGROUND: Chromosomal abnormalities are a major cause of mental retardation and multiple congenital anomalies (MCA/MR). Screening for these chromosomal imbalances has mainly been done by standard karyotyping. Previous array CGH studies on selected patients with chromosomal phenotypes and normal karyotypes suggested an incidence of 10-15% of previously unnoticed de novo chromosomal imbalances. OBJECTIVE: To report array CGH screening of a series of 140 patients (the largest published so far) with idiopathic MCA/MR but normal karyotype. RESULTS: Submicroscopic chromosomal imbalances were detected in 28 of the 140 patients (20%) and included 18 deletions, seven duplications, and three unbalanced translocations. Seventeen of 24 imbalances were confirmed de novo and 19 were assumed to be causal. Excluding subtelomeric imbalances, our study identified 11 clinically relevant interstitial submicroscopic imbalances (8%). Taking this and previously reported studies into consideration, array CGH screening with a resolution of at least 1 Mb has been undertaken on 432 patients with MCA/MR. Most imbalances are non-recurrent and spread across the genome. In at least 8.8% of these patients (38 of 432) de novo intrachromosomal alterations have been identified. CONCLUSIONS: Array CGH should be considered an essential aspect of the genetic analysis of patients with MCA/MR. In addition, in the present study three patients were mosaic for a structural chromosome rearrangement. One of these patients had monosomy 7 in as few as 8% of the cells, showing that array CGH allows detection of low grade mosaicisims.

[New insights into the genetic basis of neuroblastoma]. / Nieuwe inzichten in de genetische basis van neuroblastoom.

Neuroblastoma (NB) is, next to acute lymphoblastic leukaemia, brain tumours and lymphoma the most frequent paediatric tumour (8-10%). Our research group aims to contribute to the unravelling of the genetic basis of NB. Insight into the genes and signalling pathways involved in tumour formation and development can represent an essential step towards the development of more efficient molecular targeted therapies. A first part of our research work was devoted to the analysis of genomic alterations in NB. By means of a new highly sensitive method for detecting gains and losses of chromosomal segments, we recognised three major prognostic relevant genomic subtypes of NB. In addition smaller subgroups with deviating genomic patterns were detected. In addition, this work yielded important information regarding delineation of critical regions of gain and loss in NB which should facilitate further selection of candidate oncogenes or tumour suppressor genes. A second important part of our work focussed on the gene expression profiling of NB precursor cells. We were able as the first to isolate these cells and determine their transcriptome, a finding of fundamental importance for future expression studies in NB. Another study focussed on the identification of MYCN transcriptional target genes. Gene expression analyses of model systems developed in our lab and of a large panel of cell lines and tumours allowed us to subtract a list of candidate genes which are now under further study. Finally, we initiated research towards the understanding of the role of methylation in NB oncogenesis. From this, we were able to create a list of potentially relevant methylated genes in NB. From the above it is clear that our team has made important contributions to the understanding of the complex biology and clinical behaviour of NB. Also, a broad technically innovative research platform has been developed which will allow us to dissect NB genetics with greater speed and accuracy.

Genome wide measurement of DNA copy number changes in neuroblastoma: dissecting amplicons and mapping losses, gains and breakpoints.

In the past few years high throughput methods for assessment of DNA copy number alterations have witnessed rapid progress. Both 'in house' developed BAC, cDNA, oligonucleotide and commercial arrays are now available and widely applied in the study of the human genome, particularly in the context of disease. Cancer cells are known to exhibit DNA losses, gains and amplifications affecting tumor suppressor genes and proto-oncogenes. Moreover, these patterns of genomic imbalances may be associated with particular tumor types or subtypes and may have prognostic value. Here we summarize recent array CGH findings in neuroblastoma, a pediatric tumor of the sympathetic nervous system. A total of 176 primary tumors and 53 cell lines have been analyzed on different platforms. Through these studies the genomic content and boundaries of deletions, gains and amplifications were characterized with unprecedented accuracy. Furthermore, in conjunction with cytogenetic findings, array CGH allows the mapping of breakpoints of unbalanced translocations at a very high resolution.

A new recurrent inversion, inv(7)(p15q34), leads to transcriptional activation of HOXA10 and HOXA11 in a subset of T-cell acute lymphoblastic leukemias.

Chromosomal translocations with breakpoints in T-cell receptor (TCR) genes are recurrent in T-cell malignancies. These translocations involve the TCRalphadelta gene (14q11), the TCRbeta gene (7q34) and to a lesser extent the TCRgamma gene at chromosomal band 7p14 and juxtapose T-cell oncogenes next to TCR regulatory sequences leading to deregulated expression of those oncogenes. Here, we describe a new recurrent chromosomal inversion of chromosome 7, inv(7)(p15q34), in a subset of patients with T-cell acute lymphoblastic leukemia characterized by CD2 negative and CD4 positive, CD8 negative blasts. This rearrangement juxtaposes the distal part of the HOXA gene cluster on 7p15 to the TCRbeta locus on 7q34. Real time quantitative PCR analysis for all HOXA genes revealed high levels of HOXA10 and HOXA11 expression in all inv(7) positive cases. This is the first report of a recurrent chromosome rearrangement targeting the HOXA gene cluster in T-cell malignancies resulting in deregulated HOXA gene expression (particularly HOXA10 and HOXA11) and is in keeping with a previous report suggesting HOXA deregulation in MLL-rearranged T- and B cell lymphoblastic leukemia as the key factor in leukaemic transformation. Finally, our observation also supports the previous suggested role of HOXA10 and HOXA11 in normal thymocyte development.

Quantification of splice variants using real-time PCR.

A reliable and robust method for measuring the expression of alternatively spliced transcripts is an important step in investigating the significance of each variant. So far, accurate quantification of splice variants has been laborious and difficult due to the intrinsic limitations of conventional methods. The many advantages of real-time PCR have made this technique attractive to study its application in quantification of splice isoforms. We use skipping of exon 37 in the NF1 gene as a model to compare and evaluate the different strategies for quantitating splice variants using real-time PCR. An overview of three different possibilities for detecting alternative transcripts is given. We propose the use of a boundary-spanning primer to quantify isoforms that differ greatly in abundance. We describe here a novel method for creating a reliable standard curve using one plasmid containing both alternative transcripts. In addition, we validate the use of an absolute standard curve based on a dilution series of fluorometrically quantified PCR products.

Characteristic pattern of chromosomal gains and losses in Merkel cell carcinoma detected by comparative genomic hybridization.

Merkel cell carcinoma or small cell carcinoma of the skin is a rare skin cancer seen in increasing numbers in Queensland, Australia. In its clinical course and histopathology, it resembles small cell lung carcinoma (SCLC). Little is known of the genetic basis of this disease except for a number of cytogenetic studies and three loss of heterozygosity studies. Therefore, comparative genomic hybridization was performed to determine the characteristic DNA gains and losses that occur in this tumor. Comparative genomic hybridization analysis of 34 specimens from 24 patients revealed a pattern of gains and losses that closely resembles that seen in SCLC. Overall frequent loss was seen for chromosomes 3p (46%), 5q (21%), 8p (21%), 10 (33%), 11q (17%), 13q (33%), and 17p (25%). Significant gains were seen for chromosomes 1 (63%), 3q (33%), 5p (38%), 8q (38%), 19 (63%), and X (41%), with smaller numbers having gains for chromosomes 6, 7, 20, and 21. In contrast to SCLC, amplification in Merkel cell carcinoma is a rare event.

RT-qPCR gene expression analysis in zebrafish: Preanalytical precautions and use of expressed repetitive elements for normalization.

Gene expression analysis is increasingly important in many fields of biological research. Understanding patterns of expressed genes is assumed to provide insight into complex regulatory networks and can lead to the identification of genes relevant to specific biological processes, including disease. Among different techniques, reverse transcription quantitative polymerase chain reaction (RT-qPCR) is currently regarded as the gold standard for targeted quantification of RNA gene expression, especially because of its high sensitivity, specificity, accuracy, and precision, and also because of its practical simplicity and processing speed. However, different critical factors can influence the outcome of RT-qPCR studies, including isolation of RNA, reverse transcription to cDNA, and data analysis. These factors need to be addressed in order to obtain biologically meaningful results. In this chapter, we describe how RT-qPCR can be used in a reliable way to successfully study differential gene expression in zebrafish. Hereby, we especially focus on how expressed repetitive elements can be employed as reference targets in zebrafish RT-qPCR studies and how they can further improve the quality of the data.

Quantitative and Functional Requirements for Bioluminescent Cancer Models.

BACKGROUND: Bioluminescent cancer models are widely used but detailed quantification of the luciferase signal and functional comparison with a non-transfected control cell line are generally lacking. In the present study, we provide quantitative and functional tests for luciferase-transfected cells. MATERIALS AND METHODS: We quantified the luciferase expression in BLM and HCT8/E11 transfected cancer cells, and examined the effect of long-term luciferin exposure. The present study also investigated functional differences between parental and transfected cancer cells. RESULTS: Our results showed that quantification of different single-cell-derived populations are superior with droplet digital polymerase chain reaction. Quantification of luciferase protein level and luciferase bioluminescent activity is only useful when there is a significant difference in copy number. Continuous exposure of cell cultures to luciferin leads to inhibitory effects on mitochondrial activity, cell growth and bioluminescence. These inhibitory effects correlate with luciferase copy number. Cell culture and mouse xenograft assays showed no significant functional differences between luciferase-transfected and parental cells. CONCLUSION: Luciferase-transfected cells should be validated by quantitative and functional assays before starting large-scale experiments.

MicroRNA-193b-3p acts as a tumor suppressor by targeting the MYB oncogene in T-cell acute lymphoblastic leukemia.

The MYB oncogene is a leucine zipper transcription factor essential for normal and malignant hematopoiesis. In T-cell acute lymphoblastic leukemia (T-ALL), elevated MYB levels can arise directly through T-cell receptor-mediated MYB translocations, genomic MYB duplications or enhanced TAL1 complex binding at the MYB locus or indirectly through the TAL1/miR-223/FBXW7 regulatory axis. In this study, we used an unbiased MYB 3'untranslated region-microRNA (miRNA) library screen and identified 33 putative MYB-targeting miRNAs. Subsequently, transcriptome data from two independent T-ALL cohorts and different subsets of normal T-cells were used to select miRNAs with relevance in the context of normal and malignant T-cell transformation. Hereby, miR-193b-3p was identified as a novel bona fide tumor-suppressor miRNA that targets MYB during malignant T-cell transformation thereby offering an entry point for efficient MYB targeting-oriented therapies for human T-ALL.

Modulation of the peroxisomal gene expression pattern by dehydroepiandrosterone and vitamin D: therapeutic implications.

Peroxisomes are ubiquitous organelles required for several metabolic functions. Their dysfunction is responsible for a group of human inherited disorders. In the search for endogenous factors regulating the peroxisomal compartment in normal liver, we treated female rats with dehydroepiandrosterone (DHEA) and 25-hydroxycholecalciferol for 1 and 6 days. Relative transcription levels of 39 selected genes were evaluated by real-time quantitative RT-PCR analysis. Catalase (peroxisomal marker)-specific activity was assayed in total liver homogenate and peroxisomes were visualized by catalase localization. DHEA induced peroxisome proliferation and raised catalase specific activity. Expression levels of 16 (of which 11 were peroxisomal) genes were altered. Pex 11, acyl-CoA oxidase,l - andd -multifunctional enzyme, thiolase 1, phytanoyl-CoA hydroxylase, 70 kDa peroxisomal membrane protein and very long chain acyl-CoA synthetase were upregulated, three others were downregulated. Vitamin D caused downregulation of six genes. Administration of vitamin D to peroxisomal disorder patients may be contraindicated. The adrenocortical hormone DHEA is a potential natural regulator of the peroxisomal compartment. Its therapeutic use in X-linked adrenoleukodystrophy, some other beta-oxidation defects and classical Refsum should be considered.

DNA methylation silences miR-132 in prostate cancer.

Silencing of microRNAs (miRNAs) by promoter CpG island methylation may be an important mechanism in prostate carcinogenesis. To screen for epigenetically silenced miRNAs in prostate cancer (PCa), we treated prostate normal epithelial and carcinoma cells with 5-aza-2'-deoxycytidine (AZA) and subsequently examined expression changes of 650 miRNAs by megaplex stemloop reverse transcription-quantitative PCR. After applying a selection strategy, we analyzed the methylation status of CpG islands upstream to a subset of miRNAs by methylation-specific PCR. The CpG islands of miR-18b, miR-132, miR-34b/c, miR-148a, miR-450a and miR-542-3p showed methylation patterns congruent with their expression modulations in response to AZA. Methylation analysis of these CpG islands in a panel of 50 human prostate carcinoma specimens and 24 normal controls revealed miR-132 to be methylated in 42% of human cancer cases in a manner positively correlated to total Gleason score and tumor stage. Expression analysis of miR-132 in our tissue panel confirmed its downregulation in methylated tumors. Re-expression of miR-132 in PC3 cells induced cell detachment followed by cell death (anoikis). Two pro-survival proteins-heparin-binding epidermal growth factor and TALIN2-were confirmed as direct targets of miR-132. The results of this study point to miR-132 as a methylation-silenced miRNA with an antimetastatic role in PCa controlling cellular adhesion.

MYCN and ALKF1174L are sufficient to drive neuroblastoma development from neural crest progenitor cells.

Neuroblastoma is an embryonal tumor with a heterogeneous clinical course. The tumor is presumed to be derived from the neural crest, but the cells of origin remain to be determined. To date, few recurrent genetic changes contributing to neuroblastoma formation, such as amplification of the MYCN oncogene and activating mutations of the ALK oncogene, have been identified. The possibility to model neuroblastoma in mice allows investigation of the cell of origin hypothesis in further detail. Here we present the evidence that murine neural crest progenitor cells can give rise to neuroblastoma upon transformation with MYCN or ALK(F1174L). For this purpose we used JoMa1, a multipotent neural crest progenitor cell line, which is kept in a viable and undifferentiated state by a tamoxifen-activated c-Myc transgene (c-MycER(T)). Expression of MYCN or ALK(F1174L), one of the oncogenic ALK variants identified in primary neuroblastomas, enabled these cells to grow independently of c-MycER(T) activity in vitro and caused formation of neuroblastoma-like tumors in vivo in contrast to parental JoMa1 cells and JoMa1 cells-expressing TrkA or GFP. Tumorigenicity was enhanced upon serial transplantation of tumor-derived cells, and tumor cells remained susceptible to the MYC-inhibitor, NBT-272, indicating that cell growth depended on functional MYCN. Our findings support neural crest progenitor cells as the precursor cells of neuroblastoma, and indicate that neuroblastomas arise as their malignant progeny.

Modulation of neuroblastoma disease pathogenesis by an extensive network of epigenetically regulated microRNAs.

MicroRNAs (miRNAs) contribute to the pathogenesis of many forms of cancer, including the pediatric cancer neuroblastoma, but the underlying mechanisms leading to altered miRNA expression are often unknown. Here, a novel integrated approach for analyzing DNA methylation coupled with miRNA and mRNA expression data sets identified 67 epigenetically regulated miRNA in neuroblastoma. A large proportion (42%) of these miRNAs was associated with poor patient survival when underexpressed in tumors. Moreover, we demonstrate that this panel of epigenetically silenced miRNAs targets a large set of genes that are overexpressed in tumors from patients with poor survival in a highly redundant manner. The genes targeted by the epigenetically regulated miRNAs are enriched for a number of biological processes, including regulation of cell differentiation. Functional studies involving ectopic overexpression of several of the epigenetically silenced miRNAs had a negative impact on neuroblastoma cell viability, providing further support to the concept that inactivation of these miRNAs is important for neuroblastoma disease pathogenesis. One locus, miR-340, induced either differentiation or apoptosis in a cell context dependent manner, indicating a tumor suppressive function for this miRNA. Intriguingly, it was determined that miR-340 is upregulated by demethylation of an upstream genomic region that occurs during the process of neuroblastoma cell differentiation induced by all-trans retinoic acid (ATRA). Further biological studies of miR-340 revealed that it directly represses the SOX2 transcription factor by targeting of its 3'-untranslated region, explaining the mechanism by which SOX2 is downregulated by ATRA. Although SOX2 contributes to the maintenance of stem cells in an undifferentiated state, we demonstrate that miR-340-mediated downregulation of SOX2 is not required for ATRA induced differentiation to occur. In summary, our results exemplify the dynamic nature of the miRNA epigenome and identify a remarkable network of miRNA/mRNA interactions that significantly contribute to neuroblastoma disease pathogenesis.

Lysine-specific demethylase 1 restricts hematopoietic progenitor proliferation and is essential for terminal differentiation.

Lysine (K)-specific demethylase 1A (LSD1/KDM1A) has been identified as a potential therapeutic target in solid cancers and more recently in acute myeloid leukemia. However, the potential side effects of a LSD1-inhibitory therapy remain elusive. Here, we show, with a newly established conditional in vivo knockdown model, that LSD1 represents a central regulator of hematopoietic stem and progenitor cells. LSD1 knockdown (LSD1-kd) expanded progenitor numbers by enhancing their proliferative behavior. LSD1-kd led to an extensive expansion of granulomonocytic, erythroid and megakaryocytic progenitors. In contrast, terminal granulopoiesis, erythropoiesis and platelet production were severely inhibited. The only exception was monopoiesis, which was promoted by LSD1 deficiency. Importantly, we showed that peripheral blood granulocytopenia, monocytosis, anemia and thrombocytopenia were reversible after LSD1-kd termination. Extramedullary splenic hematopoiesis contributed to the phenotypic reversion, and progenitor populations remained expanded. LSD1-kd was associated with the upregulation of key hematopoietic genes, including Gfi1b, Hoxa9 and Meis1, which are known regulators of the HSC/progenitor compartment. We also demonstrated that LSD1-kd abrogated Gfi1b-negative autoregulation by crossing LSD1-kd with Gfi1b:GFP mice. Taken together, our findings distinguish LSD1 as a critical regulator of hematopoiesis and point to severe, but reversible, side effects of a LSD1-targeted therapy.

Oncogenic activation of FOXR1 by 11q23 intrachromosomal deletion-fusions in neuroblastoma.

Neuroblastoma tumors frequently show loss of heterozygosity of chromosome 11q with a shortest region of overlap in the 11q23 region. These deletions are thought to cause inactivation of tumor suppressor genes leading to haploinsufficiency. Alternatively, micro-deletions could lead to gene fusion products that are tumor driving. To identify such events we analyzed a series of neuroblastomas by comparative genomic hybridization and single-nucleotide polymorphism arrays and integrated these data with Affymetrix mRNA profiling data with the bioinformatic tool R2 (http://r2.amc.nl). We identified three neuroblastoma samples with small interstitial deletions at 11q23, upstream of the forkhead-box R1 transcription factor (FOXR1). Genes at the proximal side of the deletion were fused to FOXR1, resulting in fusion transcripts of MLL-FOXR1 and PAFAH1B2-FOXR1. FOXR1 expression has only been detected in early embryogenesis. Affymetrix microarray analysis showed high FOXR1 mRNA expression exclusively in the neuroblastomas with micro-deletions and rare cases of other tumor types, including osteosarcoma cell line HOS. RNAi silencing of FOXR1 strongly inhibited proliferation of HOS cells and triggered apoptosis. Expression profiling of these cells and reporter assays suggested that FOXR1 is a negative regulator of fork-head box factor-mediated transcription. The neural crest stem cell line JoMa1 proliferates in culture conditional to activity of a MYC-ER transgene. Over-expression of the wild-type FOXR1 could functionally replace MYC and drive proliferation of JoMa1. We conclude that FOXR1 is recurrently activated in neuroblastoma by intrachromosomal deletion/fusion events, resulting in overexpression of fusion transcripts. Forkhead-box transcription factors have not been previously implicated in neuroblastoma pathogenesis. Furthermore, this is the first identification of intrachromosomal fusion genes in neuroblastoma.

Hsa-mir-145 is the top EWS-FLI1-repressed microRNA involved in a positive feedback loop in Ewing's sarcoma.

EWS-FLI1 is a chromosome translocation-derived chimeric transcription factor that has a central and rate-limiting role in the pathogenesis of Ewing's sarcoma. Although the EWS-FLI1 transcriptomic signature has been extensively characterized on the mRNA level, information on its impact on non-coding RNA expression is lacking. We have performed a genome-wide analysis of microRNAs affected by RNAi-mediated silencing of EWS-FLI1 in Ewing's sarcoma cell lines, and differentially expressed between primary Ewing's sarcoma and mesenchymal progenitor cells. Here, we report on the identification of hsa-mir-145 as the top EWS-FLI1-repressed microRNA. Upon knockdown of EWS-FLI1, hsa-mir-145 expression dramatically increases in all Ewing's sarcoma cell lines tested. Vice versa, ectopic expression of the microRNA in Ewing's sarcoma cell lines strongly reduced EWS-FLI1 protein, whereas transfection of an anti-mir to hsa-mir-145 increased the EWS-FLI1 levels. Reporter gene assays revealed that this modulation of EWS-FLI1 protein was mediated by the microRNA targeting the FLI1 3'-untranslated region. Mutual regulations of EWS-FLI1 and hsa-mir-145 were mirrored by an inverse correlation between their expression levels in four of the Ewing's sarcoma cell lines tested. Consistent with the role of EWS-FLI1 in Ewing's sarcoma growth regulation, forced hsa-mir-145 expression halted Ewing's sarcoma cell line growth. These results identify feedback regulation between EWS-FLI1 and hsa-mir-145 as an important component of the EWS-FLI1-mediated Ewing's sarcomagenesis that may open a new avenue to future microRNA-mediated therapy of this devastating malignant disease.