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 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.

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.

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.

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.

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.

MicroRNA miR-885-5p targets CDK2 and MCM5, activates p53 and inhibits proliferation and survival.

Several microRNA (miRNA) loci are found within genomic regions frequently deleted in primary neuroblastoma, including miR-885-5p at 3p25.3. In this study, we demonstrate that miR-885-5p is downregulated on loss of 3p25.3 region in neuroblastoma. Experimentally enforced miR-885-5p expression in neuroblastoma cell lines inhibits proliferation triggering cell cycle arrest, senescence and/or apoptosis. miR-885-5p leads to the accumulation of p53 protein and activates the p53 pathway, resulting in upregulation of p53 targets. Enforced miR-885-5p expression consistently leads to downregulation of cyclin-dependent kinase (CDK2) and mini-chromosome maintenance protein (MCM5). Both genes are targeted by miR-885-5p via predicted binding sites within the 3'-untranslated regions (UTRs) of CDK2 and MCM5. Transcript profiling after miR-885-5p introduction in neuroblastoma cells reveals alterations in expression of multiple genes, including several p53 target genes and a number of factors involved in p53 pathway activity. Taken together, these data provide evidence that miR-885-5p has a tumor suppressive role in neuroblastoma interfering with cell cycle progression and cell survival.

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.

An integrative genomics screen uncovers ncRNA T-UCR functions in neuroblastoma tumours.

Different classes of non-coding RNAs, including microRNAs, have recently been implicated in the process of tumourigenesis. In this study, we examined the expression and putative functions of a novel class of non-coding RNAs known as transcribed ultraconserved regions (T-UCRs) in neuroblastoma. Genome-wide expression profiling revealed correlations between specific T-UCR expression levels and important clinicogenetic parameters such as MYCN amplification status. A functional genomics approach based on the integration of multi-level transcriptome data was adapted to gain insights into T-UCR functions. Assignments of T-UCRs to cellular processes such as TP53 response, differentiation and proliferation were verified using various cellular model systems. For the first time, our results define a T-UCR expression landscape in neuroblastoma and suggest widespread T-UCR involvement in diverse cellular processes that are deregulated in the process of tumourigenesis.

Monoallelic but not biallelic loss of Dicer1 promotes tumorigenesis in vivo.

Human tumors are characterized by widespread reduction in microRNA (miRNA) expression, although it is unclear how such changes come about and whether they have an etiological role in the disease. Importantly, miRNA knockdown has been shown to enhance the tumorigenic potential of human lung adenocarcinoma cells. A defect in miRNA processing is one possible mechanism for global downregulation. To explore this possibility in more detail in vivo, we have manipulated Dicer1 gene dosage in a mouse model of retinoblastoma. We show that although monoallelic loss of Dicer1 does not affect normal retinal development, it dramatically accelerates tumor formation on a retinoblastoma-sensitized background. Importantly, these tumors retain one wild-type Dicer1 allele and exhibit only a partial decrease in miRNA processing. Accordingly, in silico analysis of human cancer genome data reveals frequent hemizygous, but not homozygous, deletions of DICER1. Strikingly, complete loss of Dicer1 function in mice did not accelerate retinoblastoma formation. miRNA profiling of these tumors identified members of the let-7 and miR-34 families as candidate tumor suppressors in retinoblastoma. We conclude that Dicer1 functions as a haploinsufficient tumor suppressor. This finding has implications for cancer etiology and cancer therapy.

MYCN/c-MYC-induced microRNAs repress coding gene networks associated with poor outcome in MYCN/c-MYC-activated tumors.

Increased activity of MYC protein-family members is a common feature in many cancers. Using neuroblastoma as a tumor model, we established a microRNA (miRNA) signature for activated MYCN/c-MYC signaling in two independent primary neuroblastoma tumor cohorts and provide evidence that c-MYC and MYCN have overlapping functions. On the basis of an integrated approach including miRNA and messenger RNA (mRNA) gene expression data we show that miRNA activation contributes to widespread mRNA repression, both in c-MYC- and MYCN-activated tumors. c-MYC/MYCN-induced miRNA activation was shown to be dependent on c-MYC/MYCN promoter binding as evidenced by chromatin immunoprecipitation. Finally, we show that pathways, repressed through c-MYC/MYCN miRNA activation, are highly correlated to tumor aggressiveness and are conserved across different tumor entities suggesting that c-MYC/MYCN activate a core set of miRNAs for cooperative repression of common transcriptional programs related to disease aggressiveness. Our results uncover a widespread correlation between miRNA activation and c-MYC/MYCN-mediated coding gene expression modulation and further substantiate the overlapping functions of c-MYC and MYCN in the process of tumorigenesis.

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.

Escape from p53-mediated tumor surveillance in neuroblastoma: switching off the p14(ARF)-MDM2-p53 axis.

A primary failsafe program against unrestrained proliferation and oncogenesis is provided by the p53 tumor suppressor protein, inactivation of which is considered as a hallmark of cancer. Intriguingly, mutations of the TP53 gene are rarely encountered in neuroblastoma tumors, suggesting that alternative p53-inactivating lesions account for escape from p53 control in this childhood malignancy. Several recent studies have shed light on the mechanisms by which neuroblastoma cells circumvent the p53-driven antitumor barrier. We review here these mechanisms for evasion of p53-mediated growth control and conclude that deregulation of the p14(ARF)-MDM2-p53 axis seems to be the principal mode of p53 inactivation in neuroblastoma, opening new perspectives for targeted therapeutic intervention.

[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.