Expressional alterations in functional ultra-conserved non-coding rnas in response to all-transretinoic acid - induced differentiation in neuroblastoma cells
© Watters et al.; licensee BioMed Central Ltd. 2013
Received: 13 August 2012
Accepted: 20 March 2013
Published: 8 April 2013
Ultra-conserved regions (UCRs) are segments of the genome (≥ 200 bp) that exhibit 100% DNA sequence conservation between human, mouse and rat. Transcribed UCRs (T-UCRs) have been shown to be differentially expressed in cancers versus normal tissue, indicating a possible role in carcinogenesis. All-trans-retinoic acid (ATRA) causes some neuroblastoma (NB) cell lines to undergo differentiation and leads to a significant decrease in the oncogenic transcription factor MYCN. Here, we examine the impact of ATRA treatment on T-UCR expression and investigate the biological significance of these changes.
We designed a custom tiling microarray to profile the expression of 481 T-UCRs in sense and anti-sense orientation (962 potential transcripts) in untreated and ATRA-treated neuroblastoma cell lines (SH-SY5Y, SK-N-BE, LAN-5). Following identification of significantly differentially expressed T-UCRs, we carried out siRNA knockdown and gene expression microarray analysis to investigate putative functional roles for selected T-UCRs.
Following ATRA-induced differentiation, 32 T-UCRs were differentially expressed (16 up-regulated, 16 down-regulated) across all three cell lines. Further insight into the possible role of T-UC.300A, an independent transcript whose expression is down-regulated following ATRA was achieved by siRNA knockdown, resulting in the decreased viability and invasiveness of ATRA-responsive cell lines. Gene expression microarray analysis following knockdown of T-UC.300A revealed a number of genes whose expression was altered by changing T-UC.300A levels and that might play a role in the increased proliferation and invasion of NB cells prior to ATRA-treatment.
Our results indicate that significant numbers of T-UCRs have altered expression levels in response to ATRA. While the precise roles that T-UCRs might play in cancer or in normal development are largely unknown and an important area for future study, our findings strongly indicate that the function of non-coding RNA T-UC.300A is connected with proliferation, invasion and the inhibition of differentiation of neuroblastoma cell lines prior to ATRA treatment.
KeywordsATRA neuroblastoma Transcribed ultra-conserved regions Differentiation
Neuroblastoma (NB) is a highly heterogenous childhood cancer that arises from precursor cells of the sympathetic nervous system . Clinical behaviour of these tumors, ranging from spontaneous regression to rapid progression and death due to disease, is highly correlated with a number of genomic alterations involving ploidy, MYCN amplification (MNA), and large-scale genomic imbalances such as loss of chromosome 1p, 3p, 11q and gain of 17q. MNA and loss of heterozygosity on chromosome 11q are particularly associated with aggressive disease course and represent independent genetic subtypes of NB. Each genetic subtype of NB, such as MNA or 11q-, also has significant differences in the expression patterns of large sets of protein coding genes [2–6], and in the expression profiles of non-coding RNAs such as microRNAs [7–11].
Ultra-conserved regions (UCRs) are by definition DNA segments that are at least 200 bp in length and that are 100% conserved between human, rat and mouse genomes. Four hundred and eighty-one such regions have been identified . UCRs are comprised of three basic types – intragenic (39%), intronic (43%) and exonic (15%), which also includes ‘partly exonic’ and ‘exon containing’. Approximately 3% of UCRs are not easily classified, due to their juxtaposition with alternative splice variants of host genes and the resulting variable annotation. Calin et al.  carried out the first analysis of transcribed UCRs (T-UCRs) in cancer, demonstrating that approximately 9% of the 962 possible T-UCRs (sense + anti-sense) were aberrantly transcribed in either carcinomas or leukemias relative to normal tissue. Most significantly, the authors further demonstrated that siRNA-mediated down-regulation of one T-UCR (T-UC.73A) significantly increased apoptosis in a colorectal cancer cell line. Two recent studies have demonstrated that analysis of UCR expression signatures can also be applied to the evaluation of NB tumors [14, 15]. Differential UCR expression profiles were shown to be associated with outcome in short-term versus long-term survivors with high-risk, stage 4 NB . In addition, Mestdagh et al., found an expression signature of up-regulated T-UCRs in MNA compared to non-MNA tumors .
The synthetic retinoic acid, 13-cis-retinoic acid, is an established component of the treatment given to children with high-risk NB to reduce minimal residual disease [16, 17] and exposure of a number of NB cell lines, such as SK-N-BE, to ATRA induces neural cell differentiation along with down-regulation of MYCN . Here, we identify T-UCRs that are responsive to the retinoid, all-trans-retinoic acid (ATRA), across three ATRA-sensitive cell lines and investigate the functional role of the deregulated transcript T-UC.300A.
Previous studies analyzing UCR expression in NB have used qPCR, involving reverse transcription with random primers, which is unable to distinguish between transcripts originating from the sense or the anti-sense genomic strand. Our approach involved the construction of tiling arrays for 962 UCR regions, allowing for the detection of both sense and anti-sense transcripts and for expression of host genes. We identified and validated a number of T-UCRs that are differentially expressed following ATRA-induced differentiation. Further insight into a functional role for T-UC.300A was achieved by siRNA knockdown resulting in the decreased viability and invasiveness of ATRA-responsive cell lines. As T-UC.300A is down-regulated following ATRA treatment, our findings strongly indicate that its function is connected with the increased proliferation and invasion of NB cells prior to ATRA-treatment.
The NB cell lines SK-N-BE and SH-SY5Y were obtained from the American Type Culture Collection (ATCC). Cell culture medium consisted of Ham’s F12 and EMEM (50:50) supplemented with fetal bovine serum (10%), non-essential amino acids (0.5%), L-glutamine (0.5%) and penicillin/streptomycin (1%). The LAN-5 cell line was obtained from the Children’s Oncology Group Repository. LAN-5 culture medium consisted of RPMI supplemented with penicillin/streptomycin (1%). All cell culture reagents were obtained from GIBCO.
All-trans retinoic acid (ATRA) was administered daily (5 μM final conc.) to cells over a period of 7 days. Treated and untreated cells were fixed using 4% paraformaldhyde (Sigma) and permeabilised in 0.5% Triton X-100. Cells were probed using the neuronal marker βIII Tubulin (Abcam), and subsequently were incubated with the flourescein-conjugated goat anti-rabbit Alexa Flour 488 antibody (Invitrogen). Cells were washed in PBS and then counterstained using DAPI. The Nikon TE2000s Fluorescence microscope was used to examine the cells and photographs were taken with the Hamamatsu (Orca 285) CCD Camera.
The UCR custom tiling microarray was developed using Roche NimbleGens 4 × 72 K array, composed of four identical subarrays, tiling 962 ultra-conserved sequences (481 sequences in both sense and anti-sense orientation). Oligo lengths ranged from 50mer to 72mer, in order to maintain a similar Tm across all probes. The genomic coordinates of the 481 UCRs (Build HG17, May 2004) were obtained from http://users.soe.ucsc.edu/~jill/ultra.html and converted to Build HG18, using UCSC’s Batch Coordinate Conversion tool (http://genome.ucsc.edu/cgi-bin/hgLiftOver). In addition to the ultra-conserved sequences themselves, tiling coverage also spanned 2500 bases upstream and 500 bases downstream of the location of each UCR on both sense and anti-sense strands. Probes were designed in two separate containers (sense and anti-sense) on each subarray to facilitate independent data analysis. Array design is available in ArrayExpress– Accession number A-MEXP-1899.
For gene expression, the Homo Sapiens 4 × 72 K gene expression array from Roche NimbleGen was used.
The QIAGEN RNeasy Mini Kit (Cat. No. 74101) was used to extract RNA from cells (untreated at Day 0 and ATRA-treated at Day 7). DNase treatment was included to ensure complete removal of any genomic DNA that could affect results by also hybridising to the microarrays. RNA integrity was confirmed with the Agilent RNA Nano 6000 kit (Cat. No. 5067–1511) and an Agilent Bioanalyzer. Only RNA with an RNA Integrity Number (RIN) of >8 was used for microarray analysis. For tiling microarrays, double-stranded cDNA was synthesised from 3 ug total RNA using the ExpressArt TRinucleotide mRNA Amplification Micro Kit (AmpTec). In vitro transcription using the ExpressArt AminoAllyl Add-on Module (AmpTec) generated aminoallyl modified-anti-sense RNA (aRNA), which was subsequently incubated with NHS-Cy3 (Amersham). Following purification, 4 μg Cy3-aRNA was hybridised to microarrays.
For gene expression microarray analysis following siRNA knockdown of T-UC.300A or ATRA treatment, sample preparation was carried out as previously described .
Reverse transcription for transcribed UCRs was carried out on 1ug total RNA with gene-specific primers and the SuperScript III First-Strand Synthesis System for RT-PCR (Invitrogen) in a total reaction volume of 20 ul (T-UC.324: 5′ CCCCATCCCATATGACACTC 3′; T-UC.300A: 5′ AAAAGTGGAAATCAATTTTGAAGG 3′). Real-time PCR was carried out using custom designed TaqMan assays for T-UC.324 and T-UC.300A (Applied Biosystems).
Total protein was isolated from cells using a radioimmunoprecipitation assay (RIPA) lysis buffer (Sigma). Cell pellets were washed with PBS and solubilized in RIPA for 30 mins. Protein concentration was measured using the BCA assay from Pierce. Proteins were fractionized on 6% or 10% polyacrylamide gels, and blotted onto nitrocellulose membrane. MYCN protein and the neuronal marker β – III Tubulin were detected by Western Blot using the mouse monoclonal antibody SC-53993 (Santa Cruz) and the rabbit polyclonal antibody AB-8191 (Abcam) respectively.
siRNAs against T-UC.324 and T-UC.300A were designed using Dharmacons siRNA Design Centre. siRNAs were as follows:
T-UC.324: TTACCTAACCAGTGATTAA (sense strand sequence)
T-UC.300A: ATTCATGGATGGAGATTGA (sense strand sequence)
Cells were transfected with T-UCR siRNAs (final concentration 50 nM) or negative control siRNA (Dharmacon Negative Control #1, final concentration 50 nM) using the transfection reagent Lipofectamine (Invitrogen). Media was changed after 24 hrs. RNA was extracted 120 hrs after transfection.
Acid phosphatase assay
Cells were transfected with siRNAs in 96-well plates using Lipofectamine, and plates were set up for timepoints 24-120 hrs. At each timepoint, the appropriate plate was washed twice with PBS. 10 mM p-nitrophenol phosphate in 0.1 M sodium acetate with 0.1% triton X-100 was added. Plates were incubated at 37°C for two hours and the reaction was stopped with 50 μL 1 M sodium hydroxide per well. Absorbance was measured at 405 nm using the Victor X3 Multi-Label Reader (Perkin Elmer).
Cell proliferation assay
The effect of siRNA knockdown of T-UC.300A on the rate of cell proliferation in SH-SY5Y cells was assayed using the Cell Proliferation ELISA, BrdU (colorimetric) kit (Roche, Cat. No. 11 647 229 001). Cells were transfected with siRNA against T-UC.300A or a scrambled control and were cultured in a tissue culture grade, flat bottom 96-well plate for 96 hours (n=4). BrdU labeling solution was added to each well (final concentration 10 μM BrdU) and the cells were re-incubated for an additional 2 hours. The labeling medium was then removed and the cells were fixed and denatured. Denatured DNA was incubated with an anti-BrdU monoclonal antibody conjugated with peroxidase for 90 minutes, followed by washing. 100 μl/well substrate solution was added and the plate was incubated at +15 to +25°C for 5 minutes. Stop solution (25 μl 1M H2SO4) was added to each well and mixed thoroughly. Absorbance of samples was measured using the Victor X3 Multi-Label Reader (Perkin Elmer) at 450 nm, reference wavelength 690 nm. Following subtraction of the absorbance value for the blank control from all readings, the change in proliferation rate induced by knockdown of T-UC.300A was determined.
Cells were transfected in a 6-well plate. 48 hours after transfection, cells were trypsinized and counted. 2.5×104 cells were added to BD BioCoat™ Growth Factor Reduced MATRIGEL™ Invasion Chambers (BD Biosciences) as per manufacturers’ instructions and incubated for a further 72 hrs. To determine the average number of invading cells, MATRIGEL inserts were stained with crystal violet and viewed under the microscope. The number of cells/field, in 5 random fields were counted at 200× magnification. Mean values of duplicate experiments were calculated and results subjected to t-test.
Data processing and bioinformatics analysis
Gene expression arrays
The mRNA expression data was analysed using NimbleGens NimbleScan software version 2.4, which applied quantile normalization to the data  and expression values were obtained using the Robust Multi-Chip Average algorithm as described by Irizarry et al. . Expressional alterations of 1.5 fold across both biological repeats were considered significant.
T-UCR tiling arrays
Microarrays were scanned using the GenePix 4000B scanner. Pair reports were generated using NimbleGen’s NimbleScan software version 2.4, and normalised by applying quantile normalisation. Median smoothing was first carried out on both the forward and reverse T-UCR tiling arrays with a window of 3 probes to remove outlying probe values. T-UCRs were considered expressed if their mean probe value was above the array background (median expression of the array). Only T-UCRs that were expressed in all 3 cell lines were considered. T-UCRs that were differentially expressed at least 1.5 fold in 2 out of 3 cell lines (n ≥ 4) were selected. For these selected T-UCRs, p-values were calculated across all 3 cell lines (n=6 biological samples) using Student’s t-test, and only T-UCRs with a p-value < 0.05 were considered significant.
Correlation of intragenic t-ucrs and host gene expression
Using the expression data from the ATRA-treated and untreated tiling arrays, intragenic T-UCRs transcribed in the same directions as their host genes were assessed for independent transcription relative to exonic expression, if exons were covered by the tiled region. Pearson’s correlation of T-UCR expression (mean probe value) and gene (mean expression of exonic probes) over experimental samples was used to assess T-UCR/Host gene expression correlation. Significance of this correlation was determined from a table of critical values for Pearson’s Correlation over various degrees of freedom.
Cluster analysis of intragenic t-ucrs expression
T-UCR expression was measured on the forward and reverse strand however this causes problems when calculating object similarity in cluster analyses i.e. every object will be represented twice in the expression dataset by two independent features sets, one on the forward strand and one on the reverse strand. The matrix was re-ordered by first re-orientating the arbitrary direction of T-UCR expression so that it was relative to the host gene with sense (S) referring to those T-UCRs transcribed with the host gene and anti-sense (A) referring to those T-UCRs transcribed in the opposite direction to the host gene. The S and A expressions for each T-UCR were then merged into a single feature set, thus the data matrix of 962 and F expression features was converted into a data matrix of 481 T-UCRS and 2F features (F S + F A ). Similarity of expression of intragenic T-UCRs was then determined using Spearman’s rank correlation. Hierarchical clustering was subsequently carried out using Ward’s metric and presented with accompanying heatmap visualization. Analysis was carried out using algorithms implemented in Java v1.6 and packages from R statistical programming language R version 2.8.0.
In order to examine the expression levels of all T-UCRs, a custom array was designed tiling across all UCRs on both plus and minus genomic strands. Thus, a potential 962 T-UCRs, transcribed from 481 genomic regions when taking into account forward and reverse transcripts, could be detected. Tiling also incorporated 2500 bases upstream and 500 bases downstream of each UCR. Discernment of the genomic strand from which each T-UCR originated was possible by synthesizing anti-sense RNA (aRNA), which is complementary to the original transcript, for Cy3 labelling and array hybridization.
Differential expression of t-ucrs following atra treatment
T-UCRs altered by ATRA
Orientation with host gene
Previously reported data
FRA<2MB (Calin et al., 2007)
HOX cluster (Calin et al., 2007)
FRA<2MB (Calin et al., 2007)
FRA<2MB (Calin et al., 2007)
Down-regulated in HCC (Calin et al., 2007)
No data available
No data available
Inferred role in differentiation for UCR49 (Mestdagh et al., 2010)
No data available
LOH (Calin et al., 2007)
FRA<2MB, HPV16<2.5MB (Calin et al., 2007)
LOH, HOX gene, Up-regulated in CRC (Calin et al., 2007)
No data available
AMPLIF (Calin et al., 2007)
No data available
Up-regulated in CRC (Calin et al., 2007) Down-regulated in CRC (Sana et al., 2012)
HOX gene (Calin et al., 2007)
Investigation of regulation through transcriptional interference
Independent intragenic T-UCRs, either sense or anti-sense, could be involved in the regulation of their host genes through different types of transcriptional interference - where transcriptional elongation has a direct, cis-acting suppressive effect on a second transcriptional process . Analysis of all intragenic T-UCRs and their host gene expression revealed seven T-UCRs whose expression was significantly anti-correlated with that of their host gene (Additional file 3: Table S1). None of the seven T-UCRs showed a consistent change in expression across cell lines following ATRA-treatment, however in individual samples T-UCR expression was significantly anti-correlated with host gene expression (p<0.05) which might suggest a form of regulation through transcriptional interference. The anti-correlation between T-UC.350 and T-UC.90 and their host genes can be seen in Figure 2B, C. Of particular interest is the T-UC.350/DACH1 anti-correlation. DACH1 is a tumor suppressor gene, whose expression is reduced in prostate and endometrial cancer correlating with increased tumor progression and invasion [24, 25]. The significant anti-correlation seen here between the two transcripts (p<0.05) suggests a possible cis regulation between the two in different cancer types.
siRNA knockdown and functional analysis of t-uc.300a
The prioritization of differentially expressed T-UCRs for functional assessment was based on the statistical significance of differential expression between untreated and ATRA-treated cells, in addition to other features of interest such as proximity to cancer-associated genomic regions or deregulation in other cancer types. The transcript of most interest from our list was T-UC.300A, which is down-regulated following ATRA treatment (Table 1). This UCR is located within a frequent loss of heterozygosity (LOH) region on chromosome arm 10q, anti-sense to the developmental gene PAX2, a gene whose under- or over-expression is associated with both pediatric and adult kidney pathology . In addition, T-UC.300A levels were previously shown as up-regulated in colorectal carcinoma tumors . Considering its down-regulation in ATRA-treated neuroblastoma cells, this suggests that T-UC.300A could play a role in both pediatric and adult cancers. T-UC.324 was selected as our second T-UCR for functional analysis. The most significantly up-regulated transcript across the three ATRA-treated cell lines (p=0.002) (Table 1, Figure 1), this ncRNA is expressed anti-sense to the Mpedd2 gene, which displays an anti-tumorigenic function in differentiated NB cells and is also highly expressed . Both T-UC.324 and T-UC.300A are independent transcripts, and not transcribed as part of a host gene.
Top 15 up- and down-regulated protein coding genes following siRNA knockdown of T-UC.300A in SH-SY5Y cells
Mean fold change T-UC.300A siRNA knockdown
Mean fold change after ATRA
Unsupervised analysis of T-UCR expression
In order to determine if T-UCRs had sub-groupings with characteristic expression behaviour, unsupervised analysis was performed on the intragenic T-UCR cell line expression profiles. The T-UCR data was pre-processed as described in the Methods section.
This study represents the first profiling and functional analysis of transcribed long ultra-conserved ncRNAs that are differentially expressed during ATRA induced differentiation of NB cell lines. Although their 100% homology between human, mouse and rat denotes possible critical roles in vertebrate development, there has been limited investigation of their roles in cancer development [13–15, 28, 29]. Here we sought to identify and investigate the functionality of T-UCRs differentially expressed in NB cell lines following ATRA-treatment. We have identified a role for T-UC.300A in the promotion of proliferation and invasion of NB cells prior to ATRA-induced differentiation. This intronic transcript is anti-sense to the developmental transcription factor PAX2, a gene whose function is essential during embryonic development and morphogenesis but whose overexpression also contributes to the pathogenesis of diseases such as renal cell carcinoma and medulloblastoma [30, 31]. There is an emerging role for ncRNAs in the regulation of HOX genes both in cis, through the process of transcription itself, and in trans through a direct role in gene regulation such as chromatin re-modelling [32, 33]. The repressive role of T-UC.300A was alleviated by siRNA knockdown arguing against a role for transcriptional interference and, indicating that this ncRNA most likely exerts its effect in trans. PAX2 expression was not detected in either ATRA-treated or untreated cells, also supporting a trans-acting role for T-UC.300A.
Gene expression microarray analysis showed that 150 genes were up-regulated following T-UC.300A knockdown, 46 of which were also up-regulated by ATRA treatment. The gene with the greatest increase in expression, INPP5D, encodes the SH2-containing inositol phosphatase 1 (SHIP1), a negative regulator of the PI3K/AKT pathway. Gene transfer of INPP5D has been shown to reduce proliferation in CD34+ cells from patients with acute myeloid leukemia (AML) . A number of the commonly differentially expressed genes are of particular interest in identifying a putative pathway for T-UC.300A function. At least two of the most highly up-regulated genes (CALCA, KLF2) have been reported as hypermethylated in a number of cancers, including neuroblastoma [35–38]. In addition, KLF2 is involved in the inhibition of proliferation and migration and is silenced in cancer by the Polycomb repressive complex 2 protein EZH2 [39, 40]. Furthermore KLF2 expression in acute promyelocytic leukemia (APL) can be restored by ATRA treatment . A recent study has shown that overexpression of the long ncRNA HOTAIR altered PRC2 occupancy in breast cancer cells . With this in mind, it is possible that the corresponding up-regulation of genes such as KLF2, CALCA and COL1A2 following knockdown of T-UC.300A could be indicative of a role for this T-UCR in chromatin re-programming possibly by targeting methylation–associated complexes such as PRC2. Therefore reduced T-UC.300A levels (either by siRNA or ATRA treatment) could co-ordinate the de-repression of these genes by reducing methylation at their promoters and explain the observed increase in gene expression.
Our results show that the expression of a number of T-UCRs is altered by ATRA-induced differentiation. However it is clear that a number of mechanisms are at work in their regulation, and consequently in the downstream regulation of T-UCR targets themselves.
This work was supported in part by grants from Cancer Research Ireland and the Children’s Medical and Research Foundation.
- Brodeur GM: Neuroblastoma: biological insights into a clinical enigma. Nat Rev Cancer. 2003, 3 (3): 203-216. 10.1038/nrc1014.View ArticlePubMedGoogle Scholar
- De Preter K: Human fetal neuroblast and neuroblastoma transcriptome analysis confirms neuroblast origin and highlights neuroblastoma candidate genes. Genome Biol. 2006, 7 (9): R84-10.1186/gb-2006-7-9-r84.View ArticlePubMedPubMed CentralGoogle Scholar
- McArdle L: Oligonucleotide microarray analysis of gene expression in neuroblastoma displaying loss of chromosome 11q. Carcinogenesis. 2004, 25 (9): 1599-1609. 10.1093/carcin/bgh173.View ArticlePubMedGoogle Scholar
- Wang Q: Integrative genomics identifies distinct molecular classes of neuroblastoma and shows that multiple genes are targeted by regional alterations in DNA copy number. Cancer Res. 2006, 66 (12): 6050-6062. 10.1158/0008-5472.CAN-05-4618.View ArticlePubMedGoogle Scholar
- Oberthuer A: Customized oligonucleotide microarray gene expression-based classification of neuroblastoma patients outperforms current clinical risk stratification. J Clin Oncol. 2006, 24 (31): 5070-5078. 10.1200/JCO.2006.06.1879.View ArticlePubMedGoogle Scholar
- Vermeulen J: Predicting outcomes for children with neuroblastoma using a multigene-expression signature: a retrospective SIOPEN/COG/GPOH study. Lancet Oncol. 2009, 10 (7): 663-671. 10.1016/S1470-2045(09)70154-8.View ArticlePubMedPubMed CentralGoogle Scholar
- Bray I: Widespread dysregulation of MiRNAs by MYCN amplification and chromosomal imbalances in neuroblastoma: association of miRNA expression with survival. PLoS One. 2009, 4 (11): e7850-10.1371/journal.pone.0007850.View ArticlePubMedPubMed CentralGoogle Scholar
- Buckley PG: Chromosomal and microRNA expression patterns reveal biologically distinct subgroups of 11q- neuroblastoma. Clin Cancer Res. 2010, 16 (11): 2971-2978. 10.1158/1078-0432.CCR-09-3215.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen Y, Stallings RL: Differential patterns of microRNA expression in neuroblastoma are correlated with prognosis, differentiation, and apoptosis. Cancer Res. 2007, 67 (3): 976-983. 10.1158/0008-5472.CAN-06-3667.View ArticlePubMedGoogle Scholar
- Mestdagh P: MYCN/c-MYC-induced microRNAs repress coding gene networks associated with poor outcome in MYCN/c-MYC-activated tumors. Oncogene. 2010, 29 (9): 1394-1404. 10.1038/onc.2009.429.View ArticlePubMedGoogle Scholar
- Schulte JH: Deep sequencing reveals differential expression of microRNAs in favorable versus unfavorable neuroblastoma. Nucleic Acids Res. 2010, 38 (17): 5919-5928. 10.1093/nar/gkq342.View ArticlePubMedPubMed CentralGoogle Scholar
- Bejerano G: Ultraconserved elements in the human genome. Science. 2004, 304 (5675): 1321-1325. 10.1126/science.1098119.View ArticlePubMedGoogle Scholar
- Calin GA: Ultraconserved regions encoding ncRNAs are altered in human leukemias and carcinomas. Cancer Cell. 2007, 12 (3): 215-229. 10.1016/j.ccr.2007.07.027.View ArticlePubMedGoogle Scholar
- Mestdagh P: An integrative genomics screen uncovers ncRNA T-UCR functions in neuroblastoma tumours. Oncogene. 2010, 29 (24): 3583-3592. 10.1038/onc.2010.106.View ArticlePubMedGoogle Scholar
- Scaruffi P: Transcribed-Ultra Conserved Region expression is associated with outcome in high-risk neuroblastoma. BMC Cancer. 2009, 9: 441-10.1186/1471-2407-9-441.View ArticlePubMedPubMed CentralGoogle Scholar
- Pahlman S: Retinoic acid-induced differentiation of cultured human neuroblastoma cells: a comparison with phorbolester-induced differentiation. Cell Differ. 1984, 14 (2): 135-144. 10.1016/0045-6039(84)90038-1.View ArticlePubMedGoogle Scholar
- Wagner LM, Danks MK: New therapeutic targets for the treatment of high-risk neuroblastoma. J Cell Biochem. 2009, 107 (1): 46-57. 10.1002/jcb.22094.View ArticlePubMedGoogle Scholar
- Thiele CJ, Reynolds CP, Israel MA: Decreased expression of N-myc precedes retinoic acid-induced morphological differentiation of human neuroblastoma. Nature. 1985, 313 (6001): 404-406. 10.1038/313404a0.View ArticlePubMedGoogle Scholar
- Das S: MicroRNA mediates DNA demethylation events triggered by retinoic acid during neuroblastoma cell differentiation. Cancer Res. 2010, 70 (20): 7874-7881. 10.1158/0008-5472.CAN-10-1534.View ArticlePubMedPubMed CentralGoogle Scholar
- Bolstad BM: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.View ArticlePubMedGoogle Scholar
- Irizarry RA: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249.View ArticlePubMedGoogle Scholar
- Kim DY: Up-regulated hoxC4 induces CD14 expression during the differentiation of acute promyelocytic leukemia cells. Leuk Lymphoma. 2005, 46 (7): 1061-1066. 10.1080/10428190500102589.View ArticlePubMedGoogle Scholar
- Callen BP: Transcriptional interference between convergent promoters caused by elongation over the promoter. Mol Cell. 2004, 14 (5): 647-656. 10.1016/j.molcel.2004.05.010.View ArticlePubMedGoogle Scholar
- Wu K: The cell fate determination factor dachshund inhibits androgen receptor signaling and prostate cancer cellular growth. Cancer Res. 2009, 69 (8): 3347-3355. 10.1158/0008-5472.CAN-08-3821.View ArticlePubMedPubMed CentralGoogle Scholar
- Nan F: Altered expression of DACH1 and cyclin D1 in endometrial cancer. Cancer Biol Ther. 2009, 8 (16): 1534-1539. 10.4161/cbt.8.16.8963.View ArticlePubMedGoogle Scholar
- Harshman LA, Brophy PD: PAX2 in human kidney malformations and disease. Pediatr Nephrol. 2011Google Scholar
- Liguori L: The metallophosphodiesterase Mpped2 impairs tumorigenesis in neuroblastoma. Cell Cycle. 2012, 11 (3): 569-81. 10.4161/cc.11.3.19063.View ArticlePubMedGoogle Scholar
- Braconi C: Expression and functional role of a transcribed noncoding RNA with an ultraconserved element in hepatocellular carcinoma. Proc Natl Acad Sci USA. 2011, 108 (2): 786-91. 10.1073/pnas.1011098108.View ArticlePubMedGoogle Scholar
- Sana J: Expression levels of transcribed ultraconserved regions uc.73 and uc.388 are altered in colorectal cancer. Oncology. 2012, 82 (2): 114-8. 10.1159/000336479.View ArticlePubMedGoogle Scholar
- Gruss P, Walther C: Pax in development. Cell. 1992, 69 (5): 719-22. 10.1016/0092-8674(92)90281-G.View ArticlePubMedGoogle Scholar
- Gnarra JR, Dressler GR: Expression of Pax-2 in human renal cell carcinoma and growth inhibition by anti-sense oligonucleotides. Cancer Res. 1995, 55 (18): 4092-8.PubMedGoogle Scholar
- Sessa L: Noncoding RNA synthesis and loss of Polycomb group repression accompanies the colinear activation of the human HOXA cluster. RNA. 2007, 13 (2): 223-39.View ArticlePubMedPubMed CentralGoogle Scholar
- Rinn JL: Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell. 2007, 129 (7): 1311-23. 10.1016/j.cell.2007.05.022.View ArticlePubMedPubMed CentralGoogle Scholar
- Metzner A: Reduced proliferation of CD34(+) cells from patients with acute myeloid leukemia after gene transfer of INPP5D. Gene Ther. 2009, 16 (4): 570-3. 10.1038/gt.2008.184.View ArticlePubMedGoogle Scholar
- Yao D: Quantitative assessment of gene methylation and their impact on clinical outcome in gastric cancer. Clin Chim Acta. 2012, 413 (7–8): 787-94.View ArticlePubMedGoogle Scholar
- Olk-Batz C: Aberrant DNA methylation characterizes juvenile myelomonocytic leukemia with poor outcome. Blood. 2011, 117 (18): 4871-80. 10.1182/blood-2010-08-298968.View ArticlePubMedGoogle Scholar
- Caren H: Identification of epigenetically regulated genes that predict patient outcome in neuroblastoma. BMC Cancer. 2011, 11: 66-10.1186/1471-2407-11-66.View ArticlePubMedPubMed CentralGoogle Scholar
- Moran A: Methylation profiling in non-small cell lung cancer: clinical implications. Int J Oncol. 2012, 40 (3): 739-46.PubMedGoogle Scholar
- Taniguchi H: Silencing of Kruppel-like factor 2 by the histone methyltransferase EZH2 in human cancer. Oncogene. 2012, 31 (15): 1988-94. 10.1038/onc.2011.387.View ArticlePubMedGoogle Scholar
- Black AR, Black JD, Azizkhan-Clifford J: Sp1 and kruppel-like factor family of transcription factors in cell growth regulation and cancer. J Cell Physiol. 2001, 188 (2): 143-60. 10.1002/jcp.1111.View ArticlePubMedGoogle Scholar
- Humbert M: Deregulated expression of Kruppel-like factors in acute myeloid leukemia. Leuk Res. 2011, 35 (7): 909-13. 10.1016/j.leukres.2011.03.010.View ArticlePubMedGoogle Scholar
- Gupta RA: Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010, 464 (7291): 1071-6. 10.1038/nature08975.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/13/184/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.