- Open Access
A large scale expression study associates uc.283-plus lncRNA with pluripotent stem cells and human glioma
Genome Medicine volume 6, Article number: 76 (2014)
There are 481 ultra-conserved regions (UCRs) longer than 200 bases in the genomes of human, mouse and rat. These DNA sequences are absolutely conserved and show 100% identity with no insertions or deletions. About half of these UCRs are reported as transcribed and many correspond to long non-coding RNAs (lncRNAs).
We used custom microarrays with 962 probes representing sense and antisense sequences for the 481 UCRs to examine their expression across 374 normal samples from 46 different tissues and 510 samples representing 10 different types of cancer. The expression in embryonic stem cells of selected UCRs was validated by real time PCR.
We identified tissue selective UCRs and studied UCRs in embryonic and induced pluripotent stem cells. Among the normal tissues, the uc.283 lncRNA was highly specific for pluripotent stem cells. Intriguingly, the uc.283-plus lncRNA was highly expressed in some solid cancers, particularly in one of the most untreatable types, glioma.
Our results suggest that uc.283-plus lncRNA might have a role in pluripotency of stem cells and in the biology of glioma.
Long non-coding RNAs (lncRNAs) are involved in many biological processes . Transcribed ultra-conserved regions (T-UCRs) are a large portion of the so-called ultra-conserved regions (UCRs). The 'ultra-conserved' term was originally proposed for genomic regions longer than 200 bp that are absolutely conserved (100% homology with no insertions or deletions) in human, mouse, and rat genomes . Many of these elements possess tissue-specific enhancer activity - and others have been shown to associate with splicing regulators. Evolutionary conservation has become a powerful tool to identify functionally important regions in the human genome . A high proportion of UCRs show extreme conservation within mammals only and nearly 47% of UCRs in human have been localized to exons of genes involved in RNA processing or in the regulation of transcription and development ,. The reasons for this extreme conservation remain a mystery, but it was proposed that UCRs play a role in the ontogeny and phylogeny of mammals and other vertebrates. This idea is supported by the identification of a distal enhancer and an ultra-conserved exon derived from a retroposon active more than 400 million years ago in lobe-finned fishes and terrestrial vertebrates, and maintained as active in the 'living fossil' coelacanth . Another study showed the concurrent presence of enhancer and transcript functions in non-exonic UCRs, and suggested that they may belong to non-coding RNAs (ncRNAs) . Recently, a positive correlation of expression with conservation and epigenetic marks was described in T-UCRs, among other ncRNAs . Despite still having largely unknown roles, T-UCRs are thus now thought to act as 'regulators' of other RNAs .
Recent studies suggested that UCRs could contribute to the development of malignancies ,. Genome-wide profiling revealed that UCRs have distinct signatures in human leukemias and carcinomas  and are frequently located at fragile sites and in cancer-associated genomic regions . Clinical findings also suggested that UCR signatures can have independent prognostic value in high-risk neuroblastoma patients  by providing additional prognostic value in conjunction with N-MYC activity/amplification . Additionally, SNPs within UCRs were associated with increased familial breast cancer risk .
We performed this large scale study to identify UCR activity in cancer, analyzing almost 900 human samples from a panel of 46 normal tissues and 10 solid cancers using a custom-made microarray platform.
UCR expression arrays
We studied the expression of UCRs using the Ohio State University Comprehensive Cancer Center (OSUCCC) custom microarray . The Gene Expression Omnibus describes the OSU-CCC 4.0 platform under accession number GPL14184. Briefly, a sense and an anti-sense 40-mer probe were designed for each of 481 UCRs. Each probe was printed in duplicate in two different slide locations, and therefore quadruplicate measures were available. Total RNA (2 μg) were used for labeling and hybridization. The microarrays were hybridized in 6X SSPE (0.9 M NaCl/60 MM NaH2PO4 · H2O/8 MM EDTA, pH 7.4)/30% formamide at 25°C for 18 h, washed in 0.75X TNT (Tris HCl/NaCl/Tween 20) at 37°C for 40 Minutes. Processed slides were scanned using a microarray scanner (Axon Molecular Devices, Sunnyvale, CA, USA), with the laser set to 635 nm, at a fixed PMT (photomultiplier tube) setting, and a scan resolution of 10 Mm. Microarray images were analyzed using GenePix Pro and post-processing was performed essentially as described earlier .
T-UCRs were retained when present in at least 20% of samples and when at least 20% of them had a fold change of more than 1.5 from the gene median. Absent calls were thresholded prior to normalization and statistical analysis. Normalization was performed by using quantiles . First, all samples were classified according to organ, tissue and cell type; the samples were then grouped in systems (Additional file 1) and cancer types (Additional file 2). To assess the specificity of UCR expression across groups we used information content (IC) . Differentially expressed RNAs were identified using t-tests over two-class experiments or F-tests over multiple classes (that is, various normal tissues) within the class comparison tool .
Computational methods and folding free energies
Computational methods were used to investigate uc.283-plus and predict secondary structure for its RNA sequence. To investigate secondary structure, we used the RNAfold web server in the freely available ViennaRNA package version 2.0 to compute the centroid secondary structure and the corresponding free energy changes for folding, the minimum free energy (MFE; kcal/mol). In order to detect putative microRNA (miRNA) target sites in the considered region, we used IntaRNA, freely available online ,. Calculation of accessibility is based on ensemble free energies. Ensemble free energies were calculated using a partition function approach assuming global folding of the ncRNA and local folding of the mRNA. For this purpose, RNAfold and RNAup are integrated into IntaRNA via the ViennaRNA library . Highly stable miRNA-target duplexes are represented as having a very low hybridization energy.
RNA expression and validation
Mouse embryonic stem cell (ESC) lines (V6.4; hybrid 129/C57Bl/6) were cultured at 37°C in 5% CO2 in specific medium and total RNA was extracted using TRIzol. RNA samples were treated with 1 U DNase I Ampl Grade (Invitrogen Life Technology, Carlsbad, CA, USA) to remove any contaminating genomic DNA.
Total RNA was reverse transcribed using random hexamers (Invitrogen Life Technologies, Carlsbad, CA, USA). T-UCR expression was quantified by quantitative PCR (qPCR) with SYBR green. All reactions were performed in triplicate and 18S rRNA was used as reference. The triplicate Ct values were averaged and normalized Ct (ΔCt) calculated. The living conditions of the mice were appropriate for their species and all mouse experiments were approved by the Institutional Animal Care and Use Committee (IACUC) and University Laboratory Animal Resources (ULAR) of The Ohio State University. The animals were euthanized as per the IACUC approved guidelines and protocols before harvesting the tissues. RNA purity was assessed by the ratio of absorbance at 260 and 280 nm (A260/280 nm) using a NanoDrop ND-1000 (NanoDrop Inc., Wilmington, DE, USA). All tissues were obtained under the guidelines of approved protocols from the Ohio State University Internal Review Board and informed consent was obtained from each subject. RNA profiles for the human tissue and cell line samples have been deposited at ArrayExpress (E-TABM-969 and E-TABM-970) and at the NCBI Gene Expression Omnibus (GSE16654). Cancers samples were deposited at ArrayExpress (E-TABM-971 for breast carcinoma, E-TABM-46 for colorectal adenocarcinoma, E-TABM-22 for lung cancer, E-TABM-343 for ovarian carcinoma, and E-TABM-49 for prostate adenocarcinoma), and at the Gene Expression Omnibus (GSE7828 for colorectal carcinoma, GSE20099 and GSE24839 for esophagus carcinoma, GSE53504 for glioma and GSE14936 for lung adenocarcinoma and squamous carcinoma).
Results and discussion
T-UCR expression in normal human tissues
We tested the expression of UCRs in 374 samples from 46 types of normal tissues, belonging to 16 histological groups . For UCR profiling we used the OSUMC microarray platform, previously validated in two large scale studies ,. This platform has probes for 481 putative T-UCRs in either genomic strand (designated `plus’ or `minus’; in some other studies defined `+’ or `+A’, respectively). A global analysis showed that only a portion of the UCRs (296 out of 962) was expressed in human tissues (P-value <0.001; Additional file 3); 48% of these were non-exonic, 26% putative exonic and 26% exonic, proportions that reflect previously published studies ,. Fifty-seven T-UCRs were transcribed bi-directionally (Additional file 4). Tissue selectivity was calculated using the IC . The most tissue-selective UCRs were represented by a group of 15 UCRs, mostly expressed in epidermis, with an IC value ranging from 3.23 to 1.70 (Figure 1; Additional file 5). These 15 T-UCRs were not located in the same genomic cluster. Another tissue-specific RNA was uc.450-plus, which is highly expressed in the central nervous system. This finding confirms previous data showing that uc.450 was expressed in the dorsal root ganglion and the neural tube . In particular, Visel et al.  hypothesized an enhancer function for hs385, which fully includes uc.450. Interestingly, the opposite probe, uc.450-minus, was specific for the respiratory system. Uc.174-minus, located on the opposite strand of the MATR3 exon, was expressed mainly in the respiratory system and epidermis. The placenta was characterized by the expression of uc.319-minus (intergenic), while uc.237-minus (intronic) was restricted to the gastrointestinal system. Uc.43-minus (intronic; our probe is homologous to Mus musculus CN668140 EST), uc.75-minus, located on the opposite strand of the last exon of ZEB2, and uc.42-plus (intronic) were expressed in both liver and the respiratory system. In addition, uc.417-minus was expressed in adipose tissue and in the gastrointestinal system. The only UCR selective for embryo was uc.283-plus, on which we focus for the rest of this study.
A single T-UCR is differentially expressed in pluripotent stem cells
Our extensive study of UCR tissue selectivity identified only one RNA specific for the embryo. To further investigate this finding we studied ESCs and induced pluripotent stem cells. Uc.283-plus, a 277 nucleotide-long sequence located at chr10:50,604,757-50,605,033, was sufficient to discriminate between adult tissues and pluripotent stem cells. Figure 2A shows the expression of uc.283-plus at distinct stages of differentiation. Uc.283-plus displayed high values in human ESCs and induced pluripotent stem cells, increased in trophoblasts at 7 and 14 days of embryoid body differentiation, and decreased in definitive endoderm to reach the lowest values in spontaneously differentiated monolayers.
To validate our microarray data that suggest uc.283-plus is pluripotent specific, we used qPCR on RNAs extracted from the mouse ESCs and RNAs from five mouse adult tissues (spleen, seminal vesicles, thymus, lung, and bone marrow). We used ESC-specific miR-302 as a positive ncRNA control ,. The qPCR (Figure 2B) confirmed strong and mouse ESC-selective expression for both miR-302 and uc.283-plus (Wilcoxon test P-value <0.01) (Additional file 6). No enhancer activity was associated in Vista Enhancer Browser  with uc.283. In a further quest to investigate the possible function of uc.283, we analyzed a 2,000 bp region surrounding this UCR using the UCSC Genome Browser. In the H1-hESC chromatin state segmentation (HMM track from ENCODE/Broad), this sequence might be an inactive or poised promoter . In addition, the edges of uc.283 overlap with the initial tracts of two open chromatin regions (OpenChrom_15681 and OpenChrom_15682). Data from the methyl 450 K bead array track (ENCODE/HAIB) showed that, in close proximity to the start of uc.283-plus, the CpG sequence was not methylated in H1-hESC and human umbilical vein endothelial cells, and partially or totally methylated in other cell types. These findings, together with the proximity to open chromatin regions, suggested that uc.283-plus could be located in a euchromatic region during embryonic development and in ESCs. Furthermore, uc.283-plus is not present in the catalogue of lncRNAs originated from divergent transcription at promoters of active protein-coding genes . In particular, there seems to be no connection between uc.283-plus and the adjacent CpG island, containing the DRGX (Dorsal root ganglia homeobox) promoter (on the opposite strand).
Among solid cancers uc.283-plus is mostly expressed in glioma
Cancer stem cells are a highly debated issue in oncology ,. Since uc.283 was associated with pluripotency, we assessed its RNA levels in about 500 tumors from several types of solid cancers. Surprisingly, uc.283-plus was over-expressed in prostate adenocarcinoma and glioma samples (Figure 3). Recently, lncRNAs were investigated in glioma to define grade and histological differentiation of the tumor . The highest levels of expression were found in glioma, considered one of the most aggressive cancers with high propensity for proliferation and tissue invasion. It is tempting to speculate that the high expression of uc.283-plus in glioma is correlated with a 'cancer stem cell phenotype', a well-studied event occurring in glioma . Interestingly, Lujambio et al.  identified an RNA in the uc.283 genomic region but transcribed from the opposite strand in various types of cancer cell lines. They also showed that uc.283-minus undergoes specific CpG island hypermethylation, suggesting that it could be regulated by epigenetic alteration. Recently, Hudson et al.  confirmed the up-regulation of uc.283-minus in a prostate cancer cell line treated with the DNA hypomethylating agent 5-azacytidine and with the histone deacetylase inhibitor trichostatin A. Hudson et al. produced a list of all the possible ucRNA-mRNA interactions based on sequence complementarity according to the thermodynamics of the loop-loop RNA interactions ,. In order to assess whether uc.283-plus has a functional role, we checked for possible interactions with the mRNAs listed but did not find any. Subsequently, we hypothesized that it has a 'sponge function' for recruitment of miRNAs or other class of ncRNAs. We examined the possible secondary structure of uc.283-plus (Figure S1A,B in Additional file 7)  and submitted the RNA sequence to the web-tool application RegRNA2.0 , analyzing predicted miRNA target sites with a score ≥150 and a free energy of -20 or less. We found that our transcript sequence could be a target of three miRNAs: hsa-miR-455-5p, has-miR-640 and has-miR-1909-3p. We verified these possible interactions using another algorithm, IntaRNA (Figure S1C in Additional file 7) . Interestingly, hsa-miR-1909-3p was discovered in human ESCs by deep sequencing of small RNA libraries  and targeted genes, such as DICER1, SOX2 and NOTCH1. Hsa-miR-455-5p and hsa-miR-640 were deregulated in several cancers but were not the most abundant isoforms and they have not been well characterized yet ,. Very recently, Liz et al.  showed that the long ncRNA uc.283-minus controlled pri-miRNA processing. This ncRNA-ncRNA interaction prevents pri-miRNA-195 cleavage by Drosha. Therefore, understanding the interactions of this kind of ncRNA is of particular importance to pinpoint their biological meaning.
We report here the genome-wide analysis of UCR tissue selectivity among 16 human histologic groups, corresponding to 46 different normal tissues. A fraction of UCRs were tissue-selective, while others were broadly expressed. Uc.283-plus was highly expressed in pluripotent ESCs and induced pluripotent stem cells. Uc.283-plus may play an important role in pluripotency and among solid cancers it is highly expressed in glioma. The understanding of the biological roles of UCRs, as those of the other lncRNAs, remains an open challenge. This study can be a starting point for the further characterization of UCR activities in normal and cancer tissues.
MG and SV designed the study and drafted the manuscript. JP and SW helped in acquisition of data. PD and SS carried out the RT-PCR. MG, MES and RZ helped with statistical analysis. MAT, GC, VC and CMC helped in analysis and interpretation of the data. RAJ, CD and RG helped to revise the manuscript. JP, MP and SV helped in writing the paper. All authors read and approved the final manuscript.
embryonic stem cell
long non-coding RNA
quantitative polymerase chain reaction
transcribed ultra-conserved region
Mercer TR, Dinger ME, Mattick JS: Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009, 10: 155-159. 10.1038/nrg2521.
Bejerano G, Pheasant M, Makunin I, Stephen S, Kent WJ, Mattick JS, Haussler D: Ultraconserved elements in the human genome. Science. 2004, 304: 1321-1325. 10.1126/science.1098119.
Bejerano G, Lowe CB, Ahituv N, King B, Siepel A, Salama SR, Rubin EM, Kent WJ, Haussler D: A distal enhancer and an ultraconserved exon are derived from a novel retroposon. Nature. 2006, 441: 87-90. 10.1038/nature04696.
Poulin F, Nobrega MA, Plajzer-Frick I, Holt A, Afzal V, Rubin EM, Pennacchio LA: In vivo characterization of a vertebrate ultraconserved enhancer. Genomics. 2005, 85: 774-781. 10.1016/j.ygeno.2005.03.003.
Pennacchio LA, Ahituv N, Moses AM, Prabhakar S, Nobrega MA, Shoukry M, Minovitsky S, Dubchak I, Holt A, Lewis KD, Plajzer-Frick I, Akiyama J, De Val S, Afzal V, Black BL, Couronne O, Eisen MB, Visel A, Rubin EM: In vivo enhancer analysis of human conserved non-coding sequences. Nature. 2006, 444: 499-502. 10.1038/nature05295.
Baira E, Greshock J, Coukos G, Zhang L: Ultraconserved elements: genomics, function and disease. RNA Biol. 2008, 5: 132-134. 10.4161/rna.5.3.6673.
Katzman S, Kern AD, Bejerano G, Fewell G, Fulton L, Wilson RK, Salama SR, Haussler D: Human genome ultraconserved elements are ultraselected. Science. 2007, 317: 915-10.1126/science.1142430.
Licastro D, Gennarino VA, Petrera F, Sanges R, Banfi S, Stupka E: Promiscuity of enhancer, coding and non-coding transcription functions in ultraconserved elements. BMC Genomics. 2010, 11: 151-10.1186/1471-2164-11-151.
Nielsen MM, Tehler D, Vang S, Sudzina F, Hedegaard J, Nordentoft I, Orntoft TF, Lund AH, Pedersen JS: Identification of expressed and conserved human noncoding RNAs. RNA. 2014, 20: 236-251. 10.1261/rna.038927.113.
Liz J, Portela A, Soler M, Gomez A, Ling H, Michlewski G, Calin GA, Guil S, Esteller M: Regulation of pri-miRNA processing by a long noncoding RNA transcribed from an ultraconserved region. Mol Cell. 2014, 55: 3-4. 10.1016/j.molcel.2014.05.005.
Fabbri M, Garzon R, Andreeff M, Kantarjian HM, Garcia-Manero G, Calin GA: MicroRNAs and noncoding RNAs in hematological malignancies: molecular, clinical and therapeutic implications. Leukemia. 2008, 22: 1095-1105. 10.1038/leu.2008.30.
Galasso M, Elena Sana M, Volinia S: Non-coding RNAs: a key to future personalized molecular therapy?. Genome Med. 2010, 2: 12-10.1186/gm133.
Calin GA, Liu CG, Ferracin M, Hyslop T, Spizzo R, Sevignani C, Fabbri M, Cimmino A, Lee EJ, Wojcik SE, Shimizu M, Tili E, Rossi S, Taccioli C, Pichiorri F, Liu X, Zupo S, Herlea V, Gramantieri L, Lanza G, Alder H, Rassenti L, Volinia S, Schmittgen TD, Kipps TJ, Negrini M, Croce CM: Ultraconserved regions encoding ncRNAs are altered in human leukemias and carcinomas. Cancer Cell. 2007, 12: 215-229. 10.1016/j.ccr.2007.07.027.
Rossi S, Sevignani C, Nnadi SC, Siracusa LD, Calin GA: Cancer-associated genomic regions (CAGRs) and noncoding RNAs: bioinformatics and therapeutic implications. Mamm Genome. 2008, 19: 526-540. 10.1007/s00335-008-9119-8.
Scaruffi P, Stigliani S, Moretti S, Coco S, De Vecchi C, Valdora F, Garaventa A, Bonassi S, Tonini GP: Transcribed-ultra conserved region expression is associated with outcome in high-risk neuroblastoma. BMC Cancer. 2009, 9: 441-10.1186/1471-2407-9-441.
Mestdagh P, Fredlund E, Pattyn F, Rihani A, Van Maerken T, Vermeulen J, Kumps C, Menten B, De Preter K, Schramm A, Schulte J, Noguera R, Schleiermacher G, Janoueix-Lerosey I, Laureys G, Powel R, Nittner D, Marine JC, Ringnér M, Speleman F, Vandesompele J: An integrative genomics screen uncovers ncRNA T-UCR functions in neuroblastoma tumours. Oncogene. 2010, 29: 3583-3592. 10.1038/onc.2010.106.
Yang R, Frank B, Hemminki K, Bartram CR, Wappenschmidt B, Sutter C, Kiechle M, Bugert P, Schmutzler RK, Arnold N, Weber BH, Niederacher D, Meindl A, Burwinkel B: SNPs in ultraconserved elements and familial breast cancer risk. Carcinogenesis. 2008, 29: 351-355. 10.1093/carcin/bgm290.
Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, Visone R, Iorio M, Roldo C, Ferracin M, Prueitt RL, Yanaihara N, Lanza G, Scarpa A, Vecchione A, Negrini M, Harris CC, Croce CM: A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci U S A. 2006, 103: 2257-2261. 10.1073/pnas.0510565103.
Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.
Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, Lin C, Socci ND, Hermida L, Fulci V, Chiaretti S, Foà R, Schliwka J, Fuchs U, Novosel A, Müller RU, Schermer B, Bissels U, Inman J, Phan Q, Chien M, Weir DB, Choksi R, De Vita G, Frezzetti D, Trompeter HI, et al: A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007, 129: 1401-1414. 10.1016/j.cell.2007.04.040.
Simon R, Lam A, Li MC, Ngan M, Menenzes S, Zhao Y: Analysis of gene expression data using BRB-ArrayTools. Cancer Inform. 2007, 3: 11-17. [http://linus.nci.nih.gov/BRB-ArrayTools.html],
Busch A, Richter AS, Backofen R: IntaRNA: efficient prediction of bacterial sRNA targets incorporating target site accessibility and seed regions. Bioinformatics. 2008, 24: 2849-2856. 10.1093/bioinformatics/btn544.
Richter AS, Schleberger C, Backofen R, Steglich C: Seed-based INTARNA prediction combined with GFP-reporter system identifies mRNA targets of the small RNA Yfr1. Bioinformatics. 2010, 26: 1-5. 10.1093/bioinformatics/btp609.
Gruber AR, Lorenz R, Bernhart SH, Neubock R, Hofacker IL: The Vienna RNA websuite. Nucleic Acids Res. 2008, 36: W70-W74. 10.1093/nar/gkn188.
Volinia S, Galasso M, Costinean S, Tagliavini L, Gamberoni G, Drusco A, Marchesini J, Mascellani N, Sana ME, Abu Jarour R, Desponts C, Teitell M, Baffa R, Aqeilan R, Iorio MV, Taccioli C, Garzon R, Di Leva G, Fabbri M, Catozzi M, Previati M, Ambs S, Palumbo T, Garofalo M, Veronese A, Bottoni A, Gasparini P, Harris CC, Visone R, Pekarsky Y, et al: Reprogramming of miRNA networks in cancer and leukemia. Genome Res. 2010, 20: 589-599. 10.1101/gr.098046.109.
Calin GA, Ferracin M, Cimmino A, Di Leva G, Shimizu M, Wojcik SE, Iorio MV, Visone R, Sever NI, Fabbri M, Iuliano R, Palumbo T, Pichiorri F, Roldo C, Garzon R, Sevignani C, Rassenti L, Alder H, Volinia S, Liu CG, Kipps TJ, Negrini M, Croce CM: A MicroRNA signature associated with prognosis and progression in chronic lymphocytic leukemia. N Engl J Med. 2005, 353: 1793-1801. 10.1056/NEJMoa050995.
Visel A, Minovitsky S, Dubchak I, Pennacchio LA: VISTA Enhancer Browser-a database of tissue-specific human enhancers. Nucleic Acids Res. 2007, 35: D88-D92. 10.1093/nar/gkl822.
Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, Zhang X, Wang L, Issner R, Coyne M, Ku M, Durham T, Kellis M, Bernstein BE: Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011, 473: 43-49. 10.1038/nature09906.
Sigova AA, Mullen AC, Molinie B, Gupta S, Orlando DA, Guenther MG, Almada AE, Lin C, Sharp PA, Giallourakis CC, Young RA: Divergent transcription of long noncoding RNA/mRNA gene pairs in embryonic stem cells. Proc Natl Acad Sci U S A. 2013, 110: 2876-2881. 10.1073/pnas.1221904110.
Magee JA, Piskounova E, Morrison SJ: Cancer stem cells: impact, heterogeneity, and uncertainty. Cancer Cell. 2012, 21: 283-296. 10.1016/j.ccr.2012.03.003.
Watkins S, Sontheimer H: Unique biology of gliomas: challenges and opportunities. Trends Neurosci. 2012, 35: 546-556. 10.1016/j.tins.2012.05.001.
Zhang X, Sun S, Pu JK, Tsang AC, Lee D, Man VO, Lui WM, Wong ST, Leung GK: Long non-coding RNA expression profiles predict clinical phenotypes in glioma. Neurobiol Dis. 2012, 48: 1-8. 10.1016/j.nbd.2012.06.004.
Gronych J, Pfister SM, Jones DT: Connect four with glioblastoma stem cell factors. Cell. 2014, 157: 525-527. 10.1016/j.cell.2014.04.001.
Lujambio A, Portela A, Liz J, Melo SA, Rossi S, Spizzo R, Croce CM, Calin GA, Esteller M: CpG island hypermethylation-associated silencing of non-coding RNAs transcribed from ultraconserved regions in human cancer. Oncogene. 2010, 29: 6390-6401. 10.1038/onc.2010.361.
Hudson RS, Yi M, Volfovsky N, Prueitt RL, Esposito D, Volinia S, Liu CG, Schetter AJ, Van Roosbroeck K, Stephens RM, Calin GA, Croce CM, Ambs S: Transcription signatures encoded by ultraconserved genomic regions in human prostate cancer. Mol Cancer. 2013, 12: 13-10.1186/1476-4598-12-13.
Lorenz R, Bernhart SH, Honer Zu Siederdissen C, Tafer H, Flamm C, Stadler PF, Hofacker IL: ViennaRNA Package 2.0. Algorithms for molecular biology. AMB. 2011, 6: 26-
Muckstein U, Tafer H, Hackermuller J, Bernhart SH, Stadler PF, Hofacker IL: Thermodynamics of RNA-RNA binding. Bioinformatics. 2006, 22: 1177-1182. 10.1093/bioinformatics/btl024.
Chang TH, Huang HY, Hsu JB, Weng SL, Horng JT, Huang HD: An enhanced computational platform for investigating the roles of regulatory RNA and for identifying functional RNA motifs. BMC Bioinformatics. 2013, 14: S4-
Bar M, Wyman SK, Fritz BR, Qi J, Garg KS, Parkin RK, Kroh EM, Bendoraite A, Mitchell PS, Nelson AM, Ruzzo WL, Ware C, Radich JP, Gentleman R, Ruohola-Baker H, Tewari M: MicroRNA discovery and profiling in human embryonic stem cells by deep sequencing of small RNA libraries. Stem Cells. 2008, 26: 2496-2505. 10.1634/stemcells.2008-0356.
Li Y, Mine T, Ioannides CG: Short GC-rich RNA similar to miR 1909 and 1915 folds in silico with the 5’-UTR and ORF of Notch and responders: potential for the elimination of cancer stem cells. Oncol Rep. 2010, 24: 1443-1453.
Saito K, Inagaki K, Kamimoto T, Ito Y, Sugita T, Nakajo S, Hirasawa A, Iwamaru A, Ishikura T, Hanaoka H, Okubo K, Onozaki T, Zama T: MicroRNA-196a is a putative diagnostic biomarker and therapeutic target for laryngeal cancer. PloS One. 2013, 8: e71480-10.1371/journal.pone.0071480.
Li X, Lu Y, Chen Y, Lu W, Xie X: MicroRNA profile of paclitaxel-resistant serous ovarian carcinoma based on formalin-fixed paraffin-embedded samples. BMC Cancer. 2013, 13: 216-10.1186/1471-2407-13-216.
Microarray analyses were performed using BRB-ArrayTools developed by Dr Richard Simon and the BRB-ArrayTools Development Team, GenePattern (Broad Institute), BioConductor and R.
SV was supported by AIRC (IG 13585), PRIN MIUR 2010 and CNR EPIGEN grants. MG was supported by Contributo UNIFE Fondi 5X1000 - ANNO 2011.
The authors declare that they have no competing interests.
Electronic supplementary material
Additional file 1: Table S1.: The number of microarray chip samples used is reported for each system type. (XLSX 11 KB)
Additional file 2: Table S2.: The number of microarray chip samples used is reported for each cancer type. (XLSX 9 KB)
Additional file 3: Table S3.: Class comparison analysis: expression levels of tissue-selective T-UCRs. (XLSX 885 KB)
Additional file 5: Table S5.: Tissue specific T-UCRs in normal tissues. Sorted by information content (IC). The genomic coordinates came from hg19. (XLSX 13 KB)
Additional file 6: Table S6.: Expression levels are calculated as 2(-ΔΔCT) using 18S rRNA as a reference. Average of at least three experiments, each PCR performed in triplicate. Standard errors (ERR.ST.) are reported. (XLSX 12 KB)
Additional file 7: Figure S1.: (A) The uc.283-plus folding structure for the minimum free energy prediction. The optimal secondary structure has a minimum free energy of -65.70 Kcal/mol. (B) The centroid secondary structure with a minimum free energy of -32.80 Kcal/mol. (C) Hsa-miR-640, hsa-miR-1909-3p and hsa-miR-455-4p target uc.283-plus in three different positions with different hybridization energy values. (TIFF 153 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Galasso, M., Dama, P., Previati, M. et al. A large scale expression study associates uc.283-plus lncRNA with pluripotent stem cells and human glioma. Genome Med 6, 76 (2014). https://doi.org/10.1186/s13073-014-0076-4
- Pluripotent Stem Cell
- Definitive Endoderm
- Open Chromatin Region
- Extreme Conservation
- Familial Breast Cancer Risk