A microRNA activity map of human mesenchymal tumors: connections to oncogenic pathways; an integrative transcriptomic study
- Elena Fountzilas†1,
- Andrew D Kelly†1,
- Antonio R Perez-Atayde4,
- Jeffrey Goldsmith2,
- Panagiotis A Konstantinopoulos1,
- Nancy Francoeur1,
- Mick Correll3,
- Renee Rubio3,
- Lan Hu3,
- Mark C Gebhardt5,
- John Quackenbush3 and
- Dimitrios Spentzos1Email author
© Fountzilas et al.; licensee BioMed Central Ltd. 2012
Received: 15 March 2012
Accepted: 6 July 2012
Published: 23 July 2012
MicroRNAs (miRNAs) are nucleic acid regulators of many human mRNAs, and are associated with many tumorigenic processes. miRNA expression levels have been used in profiling studies, but some evidence suggests that expression levels do not fully capture miRNA regulatory activity. In this study we integrate multiple gene expression datasets to determine miRNA activity patterns associated with cancer phenotypes and oncogenic pathways in mesenchymal tumors – a very heterogeneous class of malignancies.
Using a computational method, we identified differentially activated miRNAs between 77 normal tissue specimens and 135 sarcomas and we validated many of these findings with microarray interrogation of an independent, paraffin-based cohort of 18 tumors. We also showed that miRNA activity is imperfectly correlated with miRNA expression levels. Using next-generation miRNA sequencing we identified potential base sequence alterations which may explain differential activity. We then analyzed miRNA activity changes related to the RAS-pathway and found 21 miRNAs that switch from silenced to activated status in parallel with RAS activation. Importantly, nearly half of these 21 miRNAs were predicted to regulate integral parts of the miRNA processing machinery, and our gene expression analysis revealed significant reductions of these transcripts in RAS-active tumors. These results suggest an association between RAS signaling and miRNA processing in which miRNAs may attenuate their own biogenesis.
Our study represents the first gene expression-based investigation of miRNA regulatory activity in human sarcomas, and our findings indicate that miRNA activity patterns derived from integrated transcriptomic data are reproducible and biologically informative in cancer. We identified an association between RAS signaling and miRNA processing, and demonstrated sequence alterations as plausible causes for differential miRNA activity. Finally, our study highlights the value of systems level integrative miRNA/mRNA assessment with high-throughput genomic data, and the applicability of paraffin-tissue-derived RNA for validation of novel findings.
KeywordsMicroRNA Microarray RAS Mesenchymal tumors MicroRNA biogenesis
Early research on microRNAs (miRNAs) has demonstrated their critical function in a variety of neoplastic processes and has further highlighted the molecular complexity of cancer [1–5]. One of the most complex and heterogeneous cancer types is the group of malignant mesenchymal tumors (also known as sarcomas). There are few reliable biomarkers for sarcoma classification and the molecular underpinning of their heterogeneous behavior remains poorly understood [6, 7]. Early work has shown that miRNA expression levels can be used to distinguish between sarcoma subtypes . However, expression levels do not necessarily signify activity in terms of effects on their target mRNAs; there is evidence that miRNA activity can be increased irrespectively of miRNA expression levels [9, 10]. Given the increasing amount of gene expression data now available in the public domain, the concept of inferring miRNA activity using gene expression profiles as a surrogate has been proposed [9, 11]. This method combines miRNA target predictions based on sequence complementarity with concerted changes in the expression levels of corresponding target mRNAs . Thus, the output is an inferred level of miRNA regulatory activity. In this study we sought to identify miRNA activity patterns in sarcomas by integrating gene expression data from multiple sources and using a recently developed computational algorithm . On this basis, miRNAs were defined as either activated or silenced in tumors (not necessarily equivalent to over or under-expressed). We then validated potentially altered miRNAs by profiling an independent paraffin-derived sarcoma cohort and investigating their possible connection with oncogenic pathway activity. We also performed RNA-sequencing to identify possible miRNA sequence alterations and we propose a link between the RAS pathway and mature miRNA biogenesis.
Gene expression datasets
We used four public datasets, (oligonucleotide Affymetrix U133A), from Japan , Memorial Sloan Kettering Cancer Center (MSKCC) , UK  and Genomics Institute of the Novartis Research Foundation (GINRF) . Raw data were retrieved for a total of 77 normal tissue samples, including epithelial/adenoid (44), hematopoietic (1), neuroendocrine (6), gonadal (4), neural (9) and mesenchymal tissues (13), and 135 sarcoma samples (including 28 non-myxoid liposarcomas comprised of 6 well-differentiated, 3 pleomorphic, and 19 dedifferentiated, 30 round cell/myxoid liposarcomas, 16 fibrosarcomas, 30 synovial sarcomas, 20 leiomyosarcomas, and 11 osteosarcomas – available in only one dataset). The data were processed using the Robust Multi-Array Average (RMA) algorithm. Non-biological experimental variation (batch effect) between the datasets was corrected using a previously described algorithm . The compendium of these public datasets was used as a discovery set to identify candidate miRNAs with deregulated activity.
For comparison purposes we processed raw data in a similar manner from non-sarcoma datasets. Specifically, we used three publicly available ovarian cancer (Duke , Michigan , UPenn ) and three head and neck cancer datasets (UPenn , University of Medicine and Dentistry of New Jersey , UWisconsin ), all oligonoucleotide Affymetrix U133A or U133 2.0 plus.
Paraffin-based validation cohort
We used 18 formalin-fixed paraffin-embedded (FFPE) sarcoma samples from the pathology archive of Beth Israel Deaconess Medical Center (BIDMC) and Boston Children's Hospital (BCH). This work was done in accordance with a protocol for archival tissue collection and use which was approved by the Institutional Review Board (IRB) at both institutions. The requirement for a patient consent form was waived by the IRB at BIDMC. This cohort included 4 liposarcomas (all well-differentiated, non-myxoid), 3 leiomyosarcomas, 2 synovial sarcomas, and 9 osteosarcomas.
FFPE RNA isolation, whole genome and miRNA profiling
FFPE samples were cut into 1–3 mm cores. Total RNA was isolated using the Qiagen RNeasy FFPE protocol. Whole genome c-DNA-mediated annealing, selection, extension, and ligation (DASL) arrays, (Illumina, CA) containing probes for 24,000 annotated genes, were used for profiling. The DASL assay is a bead-based method for expression profiling of degraded RNA, such as that extracted from FFPE samples [23–27]. Similarly, miRNA expression profiling was performed using miRNA DASL assays, containing probes for 1146 miRNAs [28, 29]. Raw miRNA and mRNA DASL data have been deposited in NCBI’s Gene Expression Omnibus (GSE35851, and GSE35852) .
The expression profiling experiments were performed at the Molecular Genetics Core at BCH. Normalization was performed following manufacturer instructions (Genome StudioTM, Gene Expression Module v1.0 User Guide, Illumina). Background subtracted sample intensities were scaled by a factor equal to the ratio of average intensity of a virtual reference sample to the average intensity of a given sample.
Small RNA sequencing
Total RNA samples were prepared for smRNA sequencing using Illumina’s Small RNA v1.5 Sample Preparation Guide. Total RNA input ranged from 5-10 μg and first underwent 3′ and 5′ adaptor ligation followed by reverse transcription and 12 cycles of amplification on a Bio-Rad iCycler. cDNA constructs were then purified using a 6% Novex TBE PAGE gel on Invitrogen’s XCell SureLock Novex Mini-Cell System. Band sizes ranging from 80-100 bp were cut from the gel and purified. cDNA constructs were eluted from the gel and purified by ethanol precipitation according to Illumina’s protocol. Libraries were analyzed on Agilent’s 2100 Bioanalyzer with a High Sensitivity DNA Chip specific for next generation sequencing. Final libraries were immobilized onto a single read Illumina flowcell at a concentration of 12pM and underwent cluster amplification on Illumina’s Cluster Station using their DGE Small RNA Cluster Generation Kit. The amplified flowcell was then sequenced on Illumina’s GAIIx with 36 cycles of sequencing.
miRNA read mapping and quantification
The leading 21 bases were trimmed from the 36-bp reads based on the quality score and the length of mature miRNAs. The trimmed reads were mapped to miRNA precursor sequences in miRBase 16.0  to achieve more sensitive expression profiles using the software miRExpress . One base difference between the reads and the miRNA precursor sequences was allowed, which covered exact match, one gap, one base insertion, and one base difference. The number of reads mapped to a miRNA sequence was taken to represent the expression level of that miRNA.
miRNA activity algorithm
To assess miRNA activity patterns we used a recently described algorithm  designed to take a set of gene expression changes as a surrogate to determine relative miRNA activity across two conditions. The algorithm is based on the premise that expression changes of the target genes (miRanda target prediction algorithm) of a certain miRNA between two conditions reflect its activity. In brief, the expression changes are ranked in a decreasing order (expression change vector). Next, the expression change vector is screened for the distribution of genes with high binding affinity for a certain miRNA. Under the null hypothesis of no miRNA activity change, genes with high and low binding affinities will position randomly in the expression change vector. Thus, miRNA activity (or silencing) inference can be made if the distribution of gene targets for a specific miRNA is skewed on the expression vector. A positive activity score (AS) indicates the miRNA has inferred activation, while a negative activity score indicates miRNA silencing.
Estimation of false discovery rate
An estimated false discovery rate (FDR) was based on permutations of the gene expression data as previously described . In brief, for each miRNA (x) activity scores are calculated for the original data (AS(x)), and also for each of 1000 random permutations (k) of the gene labels in the mRNA expression data (NS(x, k)). NS(x, k) for all x and k is then used as the null distribution for FDR calculation for a given AS(x) = AS*. If AS* ≥ 0, the FDR estimate for miRNA x* is then defined as the ratio of the percentage of all (x, k) where NS(x, k) ≥ 0, and NS(x, k) ≥ AS*, divided by the percentage of miRNAs with AS(x) ≥ 0, where AS(x) ≥ AS*, and similarly if AS* < 0 .
Functional representational analysis
To explore biological themes in the miRNA activity patterns we used functional representational analysis, as previously described . For each biologic theme, an EASE (Expression Analysis Systematic Explorer) score is calculated based on the over-representation, or lack thereof, of genes belonging to that theme in the gene pattern discriminating two conditions. The EASE score is an adjusted Fisher’s test, further modified by the FDR method.
Predictions of RAS activation
We retrieved gene expression “read outs” of RAS activation previously validated by controlled RAS activation in vitro. These “read outs” were used to train Bayesian probit regression models of pathway activity . We applied these models to assign a probability of pathway activation in individual sarcoma samples in our study. Non-biological experimental variation between datasets was corrected using the batch effect adjustment algorithm as above. In order to afford high confidence for activity calls a probability of 0.8 was the minimum for predicted pathway activation.
Assessment of RAS-associated miRNA targets
The predicted mRNA targets of “RAS-switching” miRNAs were identified using the TargetScan and miRanda algorithms (both available online) [37, 38]. Relevant transcript levels between RAS-active and RAS-inactive tumors were compared using a 1-tailed t-test assuming heteroskedasticity.
miRNA activity in the different sarcoma histologies
miRNA activity patterns in sarcoma subtypes
VS ALL NORMAL TISSUE
VS MESENCHYMAL NORMAL TISSUE
COMMON in all histologies
Histology-specific miRNA deregulation patterns
OSTEOSARCOMA SPECIFIC miRNAs
Validation of miRNA activity patterns in a paraffin tissue cohort
Differentially active miRNAs may harbor a sequence alteration
Example of a potential miRNA sequence alteration
Counts (≤1 Bases Different)
Counts (0 Bases Different)
ACUGGACUU - GGGUCAGAAGGC
ACUGGACUU - GGGUCAGAAGGC
ACUGGACUU - GGGUCAGAAGGC
ACUGGACUU - GGGUCAGAAGGC
Differentially activated miRNAs with possible sequence alterations
VS ALL NORMAL TISSUE SAMPLES
VS MESENCHYMAL NORMAL TISSUE SAMPLES
Imperfect correlation between miRNA activity and miRNA expression levels
Sarcomas demonstrate partially different miRNA activity patterns compared to epithelial cancers
To investigate the degree to which miRNA activity patterns that we discovered are unique to sarcoma, we compared samples from three ovarian and three head and neck cancer datasets with the same normal tissue samples (from GINRF) that we had previously used for the sarcoma analysis, and identified differentially activated miRNAs. This analysis revealed that the majority of the histology specific miRNAs described above were unique to sarcoma and were not shared with the epithelial tumors (23/28 miRNAs were unique to sarcoma; Table 2).
However, we also found that the miRNAs which were commonly activated in sarcomas with respect to both all normal tissue and mesenchymal normal tissue highly overlapped (50 out of 53 miRNAs; Fisher’s exact test p < 2e-45) with the miRNAs which were commonly activated in both the HNC and Ovarian cancer tissue samples. Interestingly, the same was not true of commonly silenced miRNAs in the sarcoma subtypes. Of the 17 miRNAs commonly silenced in sarcomas with respect to both all normal and mesenchymal tissue, only 1 was also commonly silenced in both the HNC and Ovarian samples (Additional file 1: Table S4). Therefore, it appears that many activated miRNAs are common to epithelial cancers, and may represent a more general cancer phenomenon. There are, however, several silenced miRNAs which are common to all sarcoma histological subtypes which appear to be silenced only in sarcomas.
RAS pathway status is associated with miRNA activity and mature miRNA biogenesis
In order to further explore possible biological connections with important cancer pathways, we hypothesized that sarcoma phenotypes characterized by distinct activation of a known oncogenic pathway may demonstrate different miRNA activity patterns. In order to test this, we compared miRNA activity patterns between the sarcoma samples that demonstrated RAS pathway activation to those that did not. The pathway activation predictions were made based on published gene expression signatures of oncogenic pathway activation . There was some variation in the prevalence of RAS activity across histological subtypes. The fractions of “RAS active” samples were 19/28, 6/30, 9/16, 3/30, and 8/20 for non-myxoid liposarcoma, myxoid liposarcoma, fibrosarcoma, synovial sarcoma, and leiomyosarcoma respectively. Indeed, we found that both in aggregate (all subtypes taken together) and in a subtype-specific manner, samples separated by RAS activity status demonstrated different activity profiles. Specifically, we identified 42 miRNAs activated in the aggregate “RAS active” group and 30 miRNAs silenced in aggregate “RAS non-active” group (Additional file 1: Table S5, Table S6; p = 0.005 and FDR = 0.01). Among these miRNAs, 21 were present in both lists (Fisher’s exact test p < 6e-11), suggesting that these miRNAs may reverse their activity upon transition to RAS-active tumor status (“RAS-switching” miRNAs). Furthermore, it has been shown that miR-7 – one member of this list – is transcribed as a result of RAS signaling . When we examine RAS-associated miRNA activity changes by specific subtype, the results for non-myxoid liposarcomas and synovial sarcomas are largely overlapping with the aggregated analysis, however, it does not appear that there are significant RAS-switching miRNAs in the other histological subtypes (Additional file 1: Table S7).
Summary of predicted RAS-related miRNA targets
RAS + Tumors
Target Scan Predictions
EIF2C2, DICER1, TARBP2
EIF2C2, DICER1, TARBP2
Biological themes represented in distinct miRNA activity patterns
To identify other possible biological mechanisms that may be perturbed by miRNA activity changes we used predicted gene targets for each histology-specific miRNA to discover biological themes overrepresented in these target gene sets. We identified a number of biological themes that seem to be shared by the majority of the sarcoma subtypes. However, there were some unique themes in each histological subtype, for instance the extracellular matrix and inflammatory response pathways in synovial sarcoma. The full list of biological themes is presented in Additional file 1: Table S8 (EASE Score = 0.05, global FDR = 0).
miRNAs have been shown to play a critical role in many biological processes, including cell proliferation, cell cycle, differentiation and apoptosis [1–5]. Their primary function was initially thought to be the direct inhibition of translation, but they are now recognized to target mRNAs for degradation . It has been suggested that the effect of a miRNA on its target mRNA depends on the strength of their binding and the degree of sequence complementarity. Under this paradigm, perfect pairing leads to mRNA degradation, while imperfect pairing results in translation inhibition . Until recently, most miRNA studies have focused on expression levels, but clinical data on miRNA activity are lacking and it is unclear if miRNA expression levels are a good surrogate for activity.
Sarcomas – a uniquely complex group of mesenchymal tumors – are perfect candidates for exploring the regulatory role of miRNAs with the aims of better understanding their biology, and developing clinical biomarkers and therapeutic targets. To our knowledge, there is limited information on the role of miRNAs in sarcoma. Subramanian et al. used miRNA expression levels to characterize various sarcoma subtypes with distinct miRNA profiles, thereby supporting the possible importance of miRNAs in the biology of these tumors .
Our goal was to determine miRNA activity in some of the most common sarcoma subtypes with a recently developed algorithm which uses sarcoma gene expression data as a surrogate. We identified several miRNAs that appear specifically deregulated in each sarcoma subtype, using normal tissue as a comparator. Despite the technical challenges associated with confirming in silico findings, we validated the deregulated activity of many of these candidate miRNAs using a paraffin-based cohort. The majority of these miRNAs were shared in all histological subtypes, suggesting that they are perhaps related to a general neoplastic transformation. Another subset, however, appeared to be unique to each sarcoma subtype. In order to further corroborate the miRNA specificity for each sarcoma subtype, we performed miRNA activity analysis using ovarian and head and neck cancer datasets and the same normal tissue cohort as a comparator. This analysis demonstrated that the majority of the sarcoma subtype-specific miRNAs were also truly unique to sarcoma subtypes. At the same time, our findings support the notion that certain common miRNA activity changes in sarcomas may be related to a general cancer phenotype as nearly all of these miRNAs were also activated in both the ovarian and head and neck tumors.
We were then interested in uncovering potential etiologies for differential activity, one example being mature miRNA sequence alterations. Using RNA-sequencing on two of the FFPE specimens from our validation cohort we found that several miRNAs which we identified as differentially activated in all sarcomas relative to normal tissue harbor possible sequence alterations. Whether this is indicative of mutation or post-transcriptional processing is unclear because we did not perform genomic DNA sequencing, but nevertheless, an impact on miRNA activity could be explained by either phenomenon. We reason that a miRNA base deletion could conceivably lead to either increased or decreased activity because target complementarity may be either increased or decreased as a result. Another explanation for differential activity could be the presence of a chromosomal translocation. We identified the chromosomal locations of miRNAs identified as sarcoma subtype-specific in our study, and we found that miR-221, which was uniquely silenced in synovial sarcoma in our analysis, is located at Xp11.3, very near the common synovial chromosomal translocation t(X;18)(p11.2, q11.2) . Rigorously investigating all possible reasons for differential activity is beyond the scope of this study, but our findings regarding potential miRNA sequence alterations suggest that mutation, post-transcriptional modification, and/or chromosomal aberrations may play a prominent role.
To explore how differential miRNA activity may manifest characteristic phenotypic states in cancer, we evaluated the relationship between miRNA activation and RAS signaling. We categorized sarcoma samples as RAS-active versus RAS-inactive using previously validated expression “read outs” of RAS activity . The data demonstrated that, in aggregate, sarcomas with active RAS were characterized by different miRNA activity profiles compared to sarcomas without active RAS and, interestingly, a subset of miRNAs appeared to “switch” activity between the two pathway “classes.” We also examined the distributions of RAS status with respect to histological subtype and found considerable variability in the rates of RAS activation. This suggests that RAS pathway activity may be sarcoma-subtype-specific per se. Performing the activity analysis separately on each of the histological subtypes revealed that significant “RAS-switching” miRNAs were present in only the non-myxoid liposarcomas and the synovial sarcomas. Interestingly, one of these miRNAs, miR-7, has been shown to promote tumorigenesis via regulation by a mechanism in which RAS signaling increases miR-7 transcription . We propose that the increased expression of miR-7 in some RAS-active sarcomas also leads to increased miRNA activity as determined by our computational approach. A very interesting finding is that many of the identified RAS-switching miRNAs have predicted mRNA targets which encode proteins in the endogenous miRNA processing machinery. In all, nearly half of these miRNAs target one or more of the processing protein transcripts, and we confirmed significantly decreased expression of these mRNAs in RAS-active tumors. We therefore hypothesize that miRNA repression of processing proteins contributes to the observed down-regulation of some miRNAs in human tumors. This seems plausible as a similar phenomenon of Dicer regulation has been described . These observations require further work to examine whether the miRNA activity changes are contributory or causal in the RAS activation process, and to examine the link between miRNA processing machinery, RAS, and miRNA activity.
In addition to exploring miRNA activity, our study addresses the question of whether miRNA expression levels are reasonable surrogates for activity in sarcomas. Our data suggest that there is an imperfect relationship between activity and expression levels, and that it may be stronger for highly variant (in terms of expression level) miRNAs. It has been suggested that dramatic changes in miRNA levels may predictably result in activity changes, but activity can change even with small changes in expression level for various other reasons [46, 47]. For instance, functional alterations of proteins that have a role in the RNA-induced silencing complex (RISC), such as Argonaute, can cause activity changes without affecting miRNA levels . miRNA mutations can also cause altered miRNA activity while leaving the miRNA expression levels measured by microarray intact. Finally, it has been shown that certain transcripts may act as miRNA “sponges,” whereby miRNA regulatory effects may be modulated without changing their expression levels . Supporting these notions is a comparison of our findings and those of Subramanian et al. based on expression levels . We found only two synovial sarcoma-specific miRNAs, miR-126 and miR-129, that have both lower expression levels and decreased activity in both studies. While this question merits further study, these observations support the notion that expression levels and target mRNA levels capture different aspects of miRNA regulatory activity in sarcomas.
In conclusion, we present the first human specimen-based study using gene expression as a surrogate for miRNA activity patterns in sarcomas, while validating many of these miRNAs using a paraffin-embedded tissue cohort. Our analysis uncovers possible miRNA sequence alterations as a potential reason for differential activity, and we identify an association between RAS signaling and miRNA processing in which miRNAs may attenuate their own biogenesis. We show how relationships between miRNA activity and critical pathways can be assessed by high throughput genome-wide analysis. The logical next step would be a “Systems” level integration of miRNA, mRNA, and proteomic data, which would allow more comprehensive and definitive explorations of the role of miRNAs in mesenchymal tumors, and other malignancies.
Supported by the National Institutes of Health [U19 CA148065 to J.Q., K22 CA138716-01A2 to D.S.]; the Dana-Farber Cancer Institute Strategic Plan Fund [to M.C., J.Q.]; and a donation by Dr Richard and Virginia Clemmer to the sarcoma program at Beth Israel Deaconess Medical Center [to D.S., M.C.G.].
- Croce CM, Calin GA: miRNAs, cancer, and stem cell division. Cell. 2005, 13: 6-7.View ArticleGoogle Scholar
- Garzon R, Calin GA, Croce CM: MicroRNAs in Cancer. Annu Rev Med. 2009, 60: 167-179. 10.1146/annurev.med.59.053006.104707.View ArticlePubMedGoogle Scholar
- Fabbri M, Croce CM, Calin GA: MicroRNAs. Cancer J. 2008, 14: 1-6. 10.1097/PPO.0b013e318164145e.View ArticlePubMedGoogle Scholar
- Calin GA, Croce CM: MicroRNA signatures in human cancers. Nat Rev Cancer. 2006, 6: 857-866. 10.1038/nrc1997.View ArticlePubMedGoogle Scholar
- Zhang H, Li Y, Lai M: The microRNA network and tumor metastasis. Oncogene. 2010, 29: 937-948. 10.1038/onc.2009.406.View ArticlePubMedGoogle Scholar
- Thway K: Pathology of soft tissue sarcomas. Clin Oncol (R Coll Radiol). 2009, 21: 695-705. 10.1016/j.clon.2009.07.016.View ArticleGoogle Scholar
- Thway K, Fisher C: Histopathological diagnostic discrepancies in soft tissue tumours referred to a specialist centre. Sarcoma. 2009, 2009: 741975-Epub 2009 May 27PubMed CentralView ArticlePubMedGoogle Scholar
- Subramanian S, Lui WO, Lee CH, Espinosa I, Nielsen TO, Heinrich MC, Corless CL, Fire AZ, van de Rijn M: MicroRNA expression signature of human sarcomas. Oncogene. 2008, 27: 2015-2026. 10.1038/sj.onc.1210836.View ArticlePubMedGoogle Scholar
- Israel A, Sharan R, Ruppin E, Galun E: Increased microRNA activity in human cancers. PLoS One. 2009, 4: e6045-10.1371/journal.pone.0006045.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheng C, Fu X, Alves P, Gerstein M: mRNA expression profiles show differential regulatory effects of microRNAs between ER + and ER- breast cancer. Genome Biol. 2009, 10: R90-10.1186/gb-2009-10-9-r90.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheng C, Li LM: Inferring microRNA activities by combining gene expression with microRNA target prediction. PLoS One. 2008, 3: e1989-10.1371/journal.pone.0001989.PubMed CentralView ArticlePubMedGoogle Scholar
- Nakayama R, Nemoto T, Takahashi H, Ohta T, Kawai A, Seki K, Yoshida T, Toyama Y, Ichikawa H, Hasegawa T: Gene expression analysis of soft tissue sarcomas: characterization and reclassification of malignant fibrous histiocytoma. Mod Pathol. 2007, 20: 749-759. 10.1038/modpathol.3800794.View ArticlePubMedGoogle Scholar
- Detwiller KY, Fernando NT, Segal NH, Ryeom SW, D'Amore PA, Yoon SS: Analysis of hypoxia-related gene expression in sarcomas and effect of hypoxia on RNA interference of vascular endothelial cell growth factor A. Cancer Res. 2005, 65: 5881-5889. 10.1158/0008-5472.CAN-04-4078.View ArticlePubMedGoogle Scholar
- Henderson SR, Guiliano D, Presneau N, McLean S, Frow R, Vujovic S, Anderson J, Sebire N, Whelan J, Athanasou N, Flanagan AM, Boshoff C: A molecular map of mesenchymal tumors. Genome Biol. 2005, 6: R76-10.1186/gb-2005-6-9-r76.PubMed CentralView ArticlePubMedGoogle Scholar
- Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, Cooke MP, Walker JR, Hogenesch JB: A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci U S A. 2004, 101: 6062-6067. 10.1073/pnas.0400782101.PubMed CentralView ArticlePubMedGoogle Scholar
- Johnson WE, Li C, Rabinovic A: Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007, 8: 118-127. 10.1093/biostatistics/kxj037.View ArticlePubMedGoogle Scholar
- Berchuck A, Iversen ES, Lancaster JM, Pittman J, Luo J, Lee P, Murphy S, Dressman HK, Febbo PG, West M, Nevins JR, Marks JR: Patterns of gene expression that characterize long-term survival in advanced stage serous ovarian cancers. Clin Cancer Res. 2005, 11: 3686-3696. 10.1158/1078-0432.CCR-04-2398.View ArticlePubMedGoogle Scholar
- Hendrix ND, Wu R, Kuick R, Schwartz DR, Fearon ER, Cho KR: Fibroblast growth factor 9 has oncogenic activity and is a downstream target of Wnt signaling in ovarian endometrioid adenocarcinomas. Cancer Res. 2006, 66: 1354-1362. 10.1158/0008-5472.CAN-05-3694.View ArticlePubMedGoogle Scholar
- Zhang L, Volinia S, Bonome T, Calin GA, Greshock J, Yang N, Liu CG, Giannakakis A, Alexiou P, Hasegawa K, Johnstone CN, Megraw MS, Adams S, Lassus H, Huang J, Kaur S, Liang S, Sethupathy P, Leminen A, Simossis VA, Sandaltzopoulos R, Naomoto Y, Katsaros D, Gimotty PA, DeMichele A, Huang Q, Bützow R, Rustgi AK, Weber BL, Birrer MJ: Genomic and epigenetic alterations deregulate microRNA expression in human epithelial ovarian cancer. Proc Natl Acad Sci U S A. 2008, 105: 7004-7009. 10.1073/pnas.0801615105.PubMed CentralView ArticlePubMedGoogle Scholar
- O'Donnell RK, Kupferman M, Wei SJ, Singhal S, Weber R, O'Malley B, Cheng Y, Putt M, Feldman M, Ziober B, Muschel RJ: Gene expression signature predicts lymphatic metastasis in squamous cell carcinoma of the oral cavity. Oncogene. 2005, 24: 1244-1251. 10.1038/sj.onc.1208285.View ArticlePubMedGoogle Scholar
- Toruner GA, Ulger C, Alkan M, Galante AT, Rinaggio J, Wilk R, Tian B, Soteropoulos P, Hameed MR, Schwalb MN, Dermody JJ: Association between gene expression profile and tumor invasion in oral squamous cell carcinoma. Cancer Genet Cytogenet. 2004, 154: 27-35. 10.1016/j.cancergencyto.2004.01.026.View ArticlePubMedGoogle Scholar
- Pyeon D, Newton MA, Lambert PF, den Boon JA, Sengupta S, Marsit CJ, Woodworth CD, Connor JP, Haugen TH, Smith EM, Kelsey KT, Turek LP, Ahlquist P: Fundamental differences in cell cycle deregulation in human papillomavirus-positive and human papillomavirus-negative head/neck and cervical cancers. Cancer Res. 2007, 67: 4605-4619. 10.1158/0008-5472.CAN-06-3619.PubMed CentralView ArticlePubMedGoogle Scholar
- Abramovitz M, Ordanic-Kodani M, Wang Y, Li Z, Catzavelos C, Bouzyk M, Sledge GW: Moreno CS, Leyland-Jones B: Optimization of RNA extraction from FFPE tissues for expression profiling in the DASL assay. Biotechniques. 2008, 44: 417-423. 10.2144/000112703.PubMed CentralView ArticlePubMedGoogle Scholar
- April C, Klotzle B, Royce T, Wickham-Garcia E, Boyaniwsky T, Izzo J, Cox D, Jones W, Rubio R, Holton K, Matulonis U, Quackenbush J, Fan JB: Whole-genome gene expression profiling of formalin-fixed, paraffin-embedded tissue samples. PLoS One. 2009, 4: e8162-10.1371/journal.pone.0008162.PubMed CentralView ArticlePubMedGoogle Scholar
- Bibikova M, Yeakley JM, Wang-Rodriguez J, Fan JB: Quantitative expression profiling of RNA from formalin-fixed, paraffin-embedded tissues using randomly assembled bead arrays. Methods Mol Biol. 2008, 439: 159-177. 10.1007/978-1-59745-188-8_11.View ArticlePubMedGoogle Scholar
- Conway C, Mitra A, Jewell R, Randerson-Moor J, Lobo S, Nsengimana J, Edward S, Sanders DS, Cook M, Powell B, Boon A, Elliott F, de Kort F, Knowles MA, Bishop DT, Newton-Bishop J: Gene expression profiling of paraffin-embedded primary melanoma using the DASL assay identifies increased osteopontin expression as predictive of reduced relapse-free survival. Clin Cancer Res. 2009, 15: 6939-6946. 10.1158/1078-0432.CCR-09-1631.PubMed CentralView ArticlePubMedGoogle Scholar
- Waddell N, Cocciardi S, Johnson J, Healey S, Marsh A, Riley J, da Silva L, Vargas AC, Reid L, Simpson PT, Lakhani SR, Chenevix-Trench G, kConFab Investigators, kConFab Investigators: Gene expression profiling of formalin-fixed, paraffin-embedded familial breast tumours using the whole genome-DASL assay. J Pathol. 2010, 221: 452-461.PubMedGoogle Scholar
- Cunningham JM, Oberg AL, Borralho PM, Kren BT, French AJ, Wang L, Bot BM, Morlan BW, Silverstein KA, Staggs R, Zeng Y, Lamblin AF, Hilker CA, Fan JB, Steer CJ, Thibodeau SN: Evaluation of a new high-dimensional miRNA profiling platform. BMC Med Genomics. 2009, 2: 57-10.1186/1755-8794-2-57.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsao J, Yau P, Winegarden N: Profiling microRNA expression with the Illumina BeadChip platform. Methods Mol Biol. 2010, 632: 73-86. 10.1007/978-1-60761-663-4_5.View ArticlePubMedGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.PubMed CentralView ArticlePubMedGoogle Scholar
- Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39 (Database Issue): D152-D157.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang WC, Lin FM, Chang WC, Lin KY, Huang HD, Lin NS: miRExpress: Analyzing high-throughput sequencing data for profiling microRNA expression. BMC Bioinformatics. 2009, 10: 328-10.1186/1471-2105-10-328.PubMed CentralView ArticlePubMedGoogle Scholar
- Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4: R70-10.1186/gb-2003-4-10-r70.PubMed CentralView ArticlePubMedGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.PubMed CentralView ArticlePubMedGoogle Scholar
- Simon R, Lam A, Li MC, Ngan M, Menenzes D, Zhao Y: Analysis of Gene Expression Data Using BRB-Array Tools. Cancer Inform. 2007, 3: 11-17.PubMed CentralPubMedGoogle Scholar
- Bild AH, Yao G, Chang JT, Wang Q, Potti A, Chasse D, Joshi MB, Harpole D, Lancaster JM, Berchuck A, Olson JA, Marks JR, Dressman HK, West M, Nevins JR: Oncogenic pathway signatures in human cancers as a guide to targeted therapies. Nature. 2006, 439: 353-357. 10.1038/nature04296.View ArticlePubMedGoogle Scholar
- Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120: 15-20. 10.1016/j.cell.2004.12.035.View ArticlePubMedGoogle Scholar
- Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS: MicroRNA targets in Drosophila. Genome Biol. 2003, 5: R1-10.1186/gb-2003-5-1-r1.PubMed CentralView ArticlePubMedGoogle Scholar
- Konstantinopoulos PA, Fountzilas E, Goldsmith JD, Bhasin M, Pillay K, Francoeur N, Libermann TA, Gebhardt MC, Spentzos D: Analysis of multiple sarcoma expression datasets: implications for classification, oncogenic pathway activation and chemotherapy resistance. PLoS One. 2010, 5: e9747-10.1371/journal.pone.0009747.PubMed CentralView ArticlePubMedGoogle Scholar
- Liang Z, Zhou H, Zheng H, Wu J: Expression levels of microRNAs are not associated with their regulatory activities. Biol Direct. 2011, 6: 43-10.1186/1745-6150-6-43.PubMed CentralView ArticlePubMedGoogle Scholar
- Chou YT, Lin HH, Lien YC, Wang YH, Hong CF, Kao YR, Lin SC, Chang YC, Lin SY, Chen SJ, Chen HC, Yeh SD, Wu CW: EGFR Promotes Lung Tumorigenesis by Activating miR-7 through a Ras/ERK/Myc Pathway That Targets the Ets2 Transcriptional Repressor ERF. Cancer Res. 2010, 70: 8822-8831. 10.1158/0008-5472.CAN-10-0638.View ArticlePubMedGoogle Scholar
- Guo H, Ingolia NT, Weissman JS, Bartel DP: Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature. 2010, 466: 835-840. 10.1038/nature09267.PubMed CentralView ArticlePubMedGoogle Scholar
- Filipowicz W, Bhattacharyya SN, Sonenberg N: Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?. Nat Rev Genet. 2008, 9: 102-114.View ArticlePubMedGoogle Scholar
- Nagao K, Ito H, Yoshida H: Chromosomal translocation t(X;18) in human synovial sarcomas analyzed by fluorescence in situ hybridization using paraffin-embedded tissue. Am J Pathol. 1996, 148: 601-609.PubMed CentralPubMedGoogle Scholar
- Martello G, Rosato A, Ferrari F, Manfrin A, Cordenonsi M, Dupont S, Enzo E, Guzzardo V, Rondina M, Spruce T, Parenti AR, Daidone MG, Bicciato S, Piccolo S: A MicroRNA Targeting Dicer for Metastasis Control. Cell. 2010, 141: 1195-1207. 10.1016/j.cell.2010.05.017.View ArticlePubMedGoogle Scholar
- Zhang L, Huang J, Yang N, Greshock J, Megraw MS, Giannakakis A, Liang S, Naylor TL, Barchetti A, Ward MR, Yao G, Medina A, O'brien-Jenkins A, Katsaros D, Hatzigeorgiou A, Gimotty PA, Weber BL, Coukos G: microRNAs exhibit high frequency genomic alterations in human cancer. Proc Natl Acad Sci U S A. 2006, 103: 9136-9141. 10.1073/pnas.0508889103.PubMed CentralView ArticlePubMedGoogle Scholar
- Sumazin P, Yang X, Chiu HS, Chung WJ, Iyer A, Llobet-Navas D, Rajbhandari P, Bansal M, Guarnieri P, Silva J, Califano A: An extensive microRNA-mediated network of RNA-RNA interactions regulates established oncogenic pathways in glioblastoma. Cell. 2011, 147: 370-381. 10.1016/j.cell.2011.09.041.PubMed CentralView ArticlePubMedGoogle Scholar