Validation and genomic interrogation of the MET variant rs11762213 as a predictor of adverse outcomes in clear cell renal cell carcinoma.

BACKGROUND: The exonic single-nucleotide variant rs11762213 located in the MET oncogene has recently been identified as a prognostic marker in clear cell renal cell carcinoma (ccRCC). This finding was validated with The Cancer Genome Atlas (TCGA) cohort, and the biologic implications were explored. METHODS: The genotype status for rs11762213 was available for 272 patients. Paired tumor-normal data, genomic data, and clinical information were acquired from ccRCC TCGA data sets. Cancer-specific survival (CSS) was analyzed with the competing risk method, and Cox proportional hazards regression was used for the analysis of the time to recurrence (TTR). Multivariate competing risk models were fitted to adjust for the validated Mayo Clinic Stage, Size, Grade, and Necrosis (SSIGN) score. RESULTS: The variant allele of rs11762213 was detected in 10.3% of the cohort. After adjustments for the SSIGN score, the risk allele remained a significant predictor for adverse CSS (hazard ratio [HR], 3.88; 95% confidence interval [CI], 1.99-7.56; P < .0001) and for TTR (OR, 2.97; 95% CI, 1.43-6.2; P = .003). The mapping of rs11762213 to regulatory regions within the genome suggested that it might affect a DNA enhancer region. RNA and protein sequencing data for MET did not reveal differences in steady-state expression with stratification by risk allele. CONCLUSIONS: The exonic MET variant rs11762213 is an independent predictor of adverse CSS and TTR in ccRCC and should be integrated into clinical practice for prognostic stratification. Genomic analysis suggests that the single-nucleotide polymorphism may affect an enhancer region located in the coding region of MET. Further biological mechanistic interrogation is currently underway.

TCEB1-mutated renal cell carcinoma: a distinct genomic and morphological subtype.

Integrated sequencing analysis identified a group of tumors among clear cell renal cell carcinomas characterized by hotspot mutations in TCEB1 (a gene that contributes to the VHL complex to ubiquitinate hypoxia-inducible factor). We analyzed 11 tumors from two distinct cohorts with TCEB1 mutations along with an expanded cohort to assess whether these should be considered an entity distinct from clear cell renal cell carcinoma and clear cell papillary renal cell carcinoma. All tumors were characterized by hotspot mutations in TCEB1 Y79C/S/F/N or A100P. Morphological and immunohistochemical characteristics of the tumors were assessed by two experienced genitourinary pathologists. Clinical and pathological variables, copy number alterations, mutations, and expression signatures were compared with a cohort of TCEB1 wild-type tumors. All TCEB1-mutated tumors were VHL and PBRM1 wild type and contained distinct copy number profiles including loss of heterozygosity of chromosome 8, the location of TCEB1 (8q21.11). All tumors lacked the clear cell renal cell carcinoma signature 3p loss and contained distinct gene expression signatures. None of the clear cell papillary tumors harbored TCEB1 mutations. Pathologically, all TCEB1-mutated tumors shared characteristic features including thick fibromuscular bands transecting the tumor, pure clear cell cytology frequently with cells showing voluminous cytoplasm, and clear cell renal cell carcinoma-like acinar areas associated with infolding tubular and focally papillary architecture. The presence of voluminous cytoplasm, absence of luminal polarization of tumor nuclei, and lack of extensive cup-like distribution of carbonic anhydrase-IX expression distinguish it from clear cell papillary carcinoma. None of the patients developed metastases at last follow-up (median 48 months). In sum, TCEB1-mutated renal cell carcinoma is a distinct entity with recurrent hotspot mutations, specific copy number alterations, pathway activation, and characteristic morphological features. Further clinical follow-up is needed to determine whether these tumors are more indolent compared with the conventional clear cell renal cell carcinoma.

Signatures of RNA binding proteins globally coupled to effective microRNA target sites.

MicroRNAs (miRNAs) and small interfering RNAs (siRNAs), bound to Argonaute proteins (RISC), destabilize mRNAs through base-pairing with the mRNA. However, the gene expression changes after perturbations of these small RNAs are only partially explained by predicted miRNA/siRNA targeting. Targeting may be modulated by other mRNA sequence elements such as binding sites for the hundreds of RNA binding proteins (RNA-BPs) expressed in any cell, and this aspect has not been systematically explored. Across a panel of published experiments, we systematically investigated to what extent sequence motifs in 3' untranslated regions (UTRs) correlate with expression changes following transfection of small RNAs. The most significantly overrepresented motifs in down-regulated mRNAs are two novel U-rich motifs (URMs), UUUUAAA and UUUGUUU, recently discovered as binding sites for the ELAVL4 (also known as HuD) RNA-BP. Surprisingly, the most significantly overrepresented motif in up-regulated mRNAs is the heptanucleotide AU-rich element (ARE), UAUUUAU, which is known to affect mRNA stability via at least 20 different RNA-BPs. We show that destabilization mediated by the transfected miRNA is generally attenuated by ARE motifs and augmented by URM motifs. These ARE and URM signatures were confirmed in different types of published experiments covering eight different cell lines. Finally, we show that both ARE and URM motifs couple to presumed endogenous miRNA binding sites in mRNAs bound by Argonaute proteins. This is the first systematic investigation of 3' UTR motifs that globally couple to regulation by miRNAs and may potentially antagonize or cooperate with miRNA/siRNA regulation. Our results suggest that binding sites of miRNAs and RNA-BPs should be considered in combination when interpreting and predicting miRNA regulation in vivo.

The molecular diversity of Luminal A breast tumors.

Breast cancer is a collection of diseases with distinct molecular traits, prognosis, and therapeutic options. Luminal A breast cancer is the most heterogeneous, both molecularly and clinically. Using genomic data from over 1,000 Luminal A tumors from multiple studies, we analyzed the copy number and mutational landscape of this tumor subtype. This integrated analysis revealed four major subtypes defined by distinct copy-number and mutation profiles. We identified an atypical Luminal A subtype characterized by high genomic instability, TP53 mutations, and increased Aurora kinase signaling; these genomic alterations lead to a worse clinical prognosis. Aberrations of chromosomes 1, 8, and 16, together with PIK3CA, GATA3, AKT1, and MAP3K1 mutations drive the other subtypes. Finally, an unbiased pathway analysis revealed multiple rare, but mutually exclusive, alterations linked to loss of activity of co-repressor complexes N-Cor and SMRT. These rare alterations were the most prevalent in Luminal A tumors and may predict resistance to endocrine therapy. Our work provides for a further molecular stratification of Luminal A breast tumors, with potential direct clinical implications.

MicroRNA transfection and AGO-bound CLIP-seq data sets reveal distinct determinants of miRNA action.

Microarray expression analyses following miRNA transfection/inhibition and, more recently, Argonaute cross-linked immunoprecipitation (CLIP)-seq assays have been used to detect miRNA target sites. CLIP and expression approaches measure differing stages of miRNA functioning-initial binding of the miRNP complex and subsequent message repression. We use nonparametric predictive models to characterize a large number of known target and flanking features, utilizing miRNA transfection, HITS-CLIP, and PAR-CLIP data. In particular, we utilize the precise spatial information provided by CLIP-seq to analyze the predictive effect of target flanking features. We observe distinct target determinants between expression-based and CLIP-based data. Target flanking features such as flanking region conservation are an important AGO-binding determinant-we hypothesize that CLIP experiments have a preference for strongly bound miRNP-target interactions involving adjacent RNA-binding proteins that increase the strength of cross-linking. In contrast, seed-related features are major determinants in expression-based studies, but less so for CLIP-seq studies, and increased miRNA concentrations typical of transfection studies contribute to this difference. While there is a good overlap between miRNA targets detected by miRNA transfection and CLIP-seq, the detection of CLIP-seq targets is largely independent of the level of subsequent mRNA degradation. Also, models built using CLIP-seq data show strong predictive power between independent CLIP-seq data sets, but are not strongly predictive for expression change. Similarly, models built from expression data are not strongly predictive for CLIP-seq data sets, supporting the finding that the determinants of miRNA binding and mRNA degradation differ. Predictive models and results are available at http://servers.binf.ku.dk/antar/.

microRNA-146a inhibits G protein-coupled receptor-mediated activation of NF-κB by targeting CARD10 and COPS8 in gastric cancer.

BACKGROUND: Gastric cancer is the second most common cause of cancer-related death in the world. Inflammatory signals originating from gastric cancer cells are important for recruiting inflammatory cells and regulation of metastasis of gastric cancer. Several microRNAs (miRNA) have been shown to be involved in development and progression of gastric cancer. miRNA-146a (miR-146a) is a modulator of inflammatory signals, but little is known about its importance in gastric cancer. We therefore wanted to identify targets of miR-146a in gastric cancer and examine its biological roles. RESULTS: The expression of miR-146a was evaluated by quantitative PCR (qPCR) and found up-regulated in the gastrin knockout mice, a mouse model of gastric cancer, and in 73% of investigated human gastric adenocarcinomas. Expression of miR-146a by gastric cancer cells was confirmed by in situ hybridization. Global analysis of changes in mRNA levels after miR-146a transfection identified two transcripts, caspase recruitment domain-containing protein 10 (CARD10) and COP9 signalosome complex subunit 8 (COPS8), as new miR-146a targets. qPCR, Western blotting and luciferase assays confirmed these transcripts as direct miR-146a targets. CARD10 and COPS8 were shown to be part of the G protein-coupled receptor (GPCR) pathway of nuclear factor-kappaB (NF-kappaB) activation. Lysophosphatidic acid (LPA) induces NF-kappaB activation via this pathway and over-expression of miR-146a inhibited LPA-induced NF-kappaB activation, reduced LPA-induced expression of tumor-promoting cytokines and growth factors and inhibited monocyte attraction. CONCLUSIONS: miR-146a expression is up-regulated in a majority of gastric cancers where it targets CARD10 and COPS8, inhibiting GPCR-mediated activation of NF-kappaB, thus reducing expression of NF-kappaB-regulated tumor-promoting cytokines and growth factors. By targeting components of several NF-kappaB-activating pathways, miR-146a is a key component in the regulation of NF-kappaB activity.

MicroRNA-143 down-regulates Hexokinase 2 in colon cancer cells.

BACKGROUND: MicroRNAs (miRNAs) are well recognized as gene regulators and have been implicated in the regulation of development as well as human diseases. miR-143 is located at a fragile site on chromosome 5 frequently deleted in cancer, and has been reported to be down-regulated in several cancers including colon cancer. METHODS: To gain insight into the role of miR-143 in colon cancer, we used a microarray-based approach in combination with seed site enrichment analysis to identify miR-143 targets. RESULTS: As expected, transcripts down-regulated upon miR-143 overexpression had a significant enrichment of miR-143 seed sites in their 3'UTRs. Here we report the identification of Hexokinase 2 (HK2) as a direct target of miR-143. We show that re-introduction of miR-143 in the colon cancer cell line DLD-1 results in a decreased lactate secretion. CONCLUSION: We have identified and validated HK2 as a miR-143 target. Furthermore, our results indicate that miR-143 mediated down-regulation of HK2 affects glucose metabolism in colon cancer cells. We hypothesize that loss of miR-143-mediated repression of HK2 can promote glucose metabolism in cancer cells, contributing to the shift towards aerobic glycolysis observed in many tumors.

miR-449 inhibits cell proliferation and is down-regulated in gastric cancer.

BACKGROUND: Gastric cancer is the fourth most common cancer in the world and the second most prevalent cause of cancer related death. The development of gastric cancer is mainly associated with H. Pylori infection leading to a focus in pathology studies on bacterial and environmental factors, and to a lesser extent on the mechanistic development of the tumour. MicroRNAs are small non-coding RNA molecules involved in post-transcriptional gene regulation. They are found to regulate genes involved in diverse biological functions and alterations in microRNA expression have been linked to the pathogenesis of many malignancies. The current study is focused on identifying microRNAs involved in gastric carcinogenesis and to explore their mechanistic relevance by characterizing their targets. RESULTS: Invitrogen NCode miRNA microarrays identified miR-449 to be decreased in 1-year-old Gastrin KO mice and in H. Pylori infected gastric tissues compared to tissues from wild type animals. Growth rate of gastric cell lines over-expressing miR-449 was inhibited by 60% compared to controls. FACS cell cycle analysis of miR-449 over-expressing cells showed a significant increase in the sub-G1 fraction indicative of apoptosis. ß-Gal assays indicated a senescent phenotype of gastric cell lines over-expressing miR-449. Affymetrix 133v2 arrays identified GMNN, MET, CCNE2, SIRT1 and CDK6 as miR-449 targets. Luciferase assays were used to confirm GMNN, MET, CCNE2 and SIRT1 as direct targets. We also show that miR-449 over-expression activated p53 and its downstream target p21 as well as the apoptosis markers cleaved CASP3 and PARP. Importantly, qPCR analyses showed a loss of miR-449 expression in human clinical gastric tumours compared to normal tissues. CONCLUSIONS: In this study, we document a diminished expression of miR-449 in Gastrin KO mice and further confirmed its loss in human gastric tumours. We investigated the function of miR-449 by identifying its direct targets. Furthermore we show that miR-449 induces senescence and apoptosis by activating the p53 pathway.

miRMaid: a unified programming interface for microRNA data resources.

BACKGROUND: MicroRNAs (miRNAs) are endogenous small RNAs that play a key role in post-transcriptional regulation of gene expression in animals and plants. The number of known miRNAs has increased rapidly over the years. The current release (version 14.0) of miRBase, the central online repository for miRNA annotation, comprises over 10.000 miRNA precursors from 115 different species. Furthermore, a large number of decentralized online resources are now available, each contributing with important miRNA annotation and information. RESULTS: We have developed a software framework, designated here as miRMaid, with the goal of integrating miRNA data resources in a uniform web service interface that can be accessed and queried by researchers and, most importantly, by computers. miRMaid is built around data from miRBase and is designed to follow the official miRBase data releases. It exposes miRBase data as inter-connected web services. Third-party miRNA data resources can be modularly integrated as miRMaid plugins or they can loosely couple with miRMaid as individual entities in the World Wide Web. miRMaid is available as a public web service but is also easily installed as a local application. The software framework is freely available under the LGPL open source license for academic and commercial use. CONCLUSION: miRMaid is an intuitive and modular software platform designed to unify miRBase and independent miRNA data resources. It enables miRNA researchers to computationally address complex questions involving the multitude of miRNA data resources. Furthermore, miRMaid constitutes a basic framework for further programming in which microRNA-interested bioinformaticians can readily develop their own tools and data sources.

Intragenomic matching reveals a huge potential for miRNA-mediated regulation in plants.

microRNAs (miRNAs) are important post-transcriptional regulators, but the extent of this regulation is uncertain, both with regard to the number of miRNA genes and their targets. Using an algorithm based on intragenomic matching of potential miRNAs and their targets coupled with support vector machine classification of miRNA precursors, we explore the potential for regulation by miRNAs in three plant genomes: Arabidopsis thaliana, Populus trichocarpa, and Oryza sativa. We find that the intragenomic matching in conjunction with a supervised learning approach contains enough information to allow reliable computational prediction of miRNA candidates without requiring conservation across species. Using this method, we identify approximately 1,200, approximately 2,500, and approximately 2,100 miRNA candidate genes capable of extensive base-pairing to potential target mRNAs in A. thaliana, P. trichocarpa, and O. sativa, respectively. This is more than five times the number of currently annotated miRNAs in the plants. Many of these candidates are derived from repeat regions, yet they seem to contain the features necessary for correct processing by the miRNA machinery. Conservation analysis indicates that only a few of the candidates are conserved between the species. We conclude that there is a large potential for miRNA-mediated regulatory interactions encoded in the genomes of the investigated plants. We hypothesize that some of these interactions may be realized under special environmental conditions, while others can readily be recruited when organisms diverge and adapt to new niches.

Long-term results of a phase II trial of high-dose radiotherapy (60 Gy) and UFT/l-leucovorin in patients with non-resectable locally advanced rectal cancer (LARC).

BACKGROUND: Preoperative radiochemotherapy is a cornerstone in patients with non- resectable locally advanced rectal cancer (LARC). To improve outcome (number of R0 resections and survival) high-dose radiotherapy (RT) was combined with oral UFT/l-leucovorin to allow tumour regression before radical intended surgery. METHODS: Pelvic RT was delivered with megavoltage photons using a 5 field technique. RT was CT-based, given 5 days a week through one posterior field and two lateral fields (48.6 Gy/27 fractions) to encompass the primary tumour and the regional lymph nodes. In addition, the tumour bed received a concurrent boost (5.4 Gy/27 fractions) and a final boost (6 Gy/3 fractions); thus GTV received 60 Gy/30 fractions. Concurrent with RT patients received a daily dose of oral UFT 300 mg/m(2) plus l-leucovorin 22.5 mg 5/7 days (divided in three doses). RESULTS: From September 2000 to November 2004, 52 patients (median age 60 years (32-83); median PS 0 (0-2)) with LARC (36 primary, 16 recurrent) were included in this phase II study. All but three patients received the planned 60 Gy, median duration of RT was 42 days (25-49). Toxicity was very modest; only four patients had a dose reduction of UFT. No hematological toxicity of clinical significance was seen. Non-hematological toxicity grade 1 (GI-toxicity, fatigue and/or dysuria) was frequently observed but only four patients experienced grade 3 toxicity (diarrhoea and/or nausea/vomiting). Forty patients (77%) were operated (30 R0, 5 R1, 5 R2) median 55 days (27-112) after completion of RT. Seven (13%) patients had a pathological complete response (pCR). Thirty-one patients (60%) died after median 25.4 months (1.6-45.2 months). Twenty-one patients (40%) are still alive June 2007. CONCLUSIONS: Preoperative high-dose RT and concomitant UFT produces major regression in most patients with non-resectable LARC and thus a good chance of cure.

The lncRNA MIR31HG regulates p16(INK4A) expression to modulate senescence.

Oncogene-induced senescence (OIS) can occur in response to oncogenic insults and is considered an important tumour suppressor mechanism. Here we identify the lncRNA MIR31HG as upregulated in OIS and find that knockdown of MIR31HG promotes a strong p16(INK4A)-dependent senescence phenotype. Under normal conditions, MIR31HG is found in both nucleus and cytoplasm, but following B-RAF expression MIR31HG is located mainly in the cytoplasm. We show that MIR31HG interacts with both INK4A and MIR31HG genomic regions and with Polycomb group (PcG) proteins, and that MIR31HG is required for PcG-mediated repression of the INK4A locus. We further identify a functional enhancer, located between MIR31HG and INK4A, which becomes activated during OIS and interacts with the MIR31HG promoter. Data from melanoma patients show a negative correlation between MIR31HG and p16(INK4A) expression levels, suggesting a role for this transcript in cancer. Hence, our data provide a new lncRNA-mediated regulatory mechanism for the tumour suppressor p16(INK4A).

Molecular subtypes of uterine leiomyosarcoma and correlation with clinical outcome.

The molecular etiology of uterine leiomyosarcoma (ULMS) is poorly understood, which accounts for the wide disparity in outcomes among women with this disease. We examined and compared the molecular profiles of ULMS and normal myometrium (NL) to identify clinically relevant molecular subtypes. Discovery cases included 29 NL and 23 ULMS specimens. RNA was hybridized to Affymetrix U133A 2.0 transcription microarrays. Differentially expressed genes and pathways were identified using standard methods. Fourteen NL and 44 ULMS independent archival samples were used for external validation. Molecular subgroups were correlated with clinical outcome. Pathway analyses of differentially expressed genes between ULMS and NL samples identified overrepresentation of cell cycle regulation, DNA repair, and genomic integrity. External validation confirmed differential expression in 31 genes (P < 4.4 × 10(-4), Bonferroni corrected), with 84% of the overexpressed genes, including CDC7, CDC20, GTSE1, CCNA2, CCNB1, and CCNB2, participating in cell cycle regulation. Unsupervised clustering of ULMS identified two clades that were reproducibly associated with progression-free (median, 4.0 vs 26.0 months; P = .02; HR, 0.33) and overall (median, 18.2 vs 77.2 months; P = .04; HR, 0.33) survival. Cell cycle genes play a key role in ULMS sarcomagenesis, providing opportunities for therapeutic targeting. Reproducible molecular subtypes associated with clinical outcome may permit individualized adjuvant treatment after clinical trial validation.

Genome-wide analysis of noncoding regulatory mutations in cancer.

Cancer primarily develops because of somatic alterations in the genome. Advances in sequencing have enabled large-scale sequencing studies across many tumor types, emphasizing the discovery of alterations in protein-coding genes. However, the protein-coding exome comprises less than 2% of the human genome. Here we analyze the complete genome sequences of 863 human tumors from The Cancer Genome Atlas and other sources to systematically identify noncoding regions that are recurrently mutated in cancer. We use new frequency- and sequence-based approaches to comprehensively scan the genome for noncoding mutations with potential regulatory impact. These methods identify recurrent mutations in regulatory elements upstream of PLEKHS1, WDR74 and SDHD, as well as previously identified mutations in the TERT promoter. SDHD promoter mutations are frequent in melanoma and are associated with reduced gene expression and poor prognosis. The non-protein-coding cancer genome remains widely unexplored, and our findings represent a step toward targeting the entire genome for clinical purposes.

Down-Regulation of miR-129-5p and the let-7 Family in Neuroendocrine Tumors and Metastases Leads to Up-Regulation of Their Targets Egr1, G3bp1, Hmga2 and Bach1.

Expression of miRNAs in Neuroendocrine Neoplasms (NEN) is poorly characterized. We therefore wanted to examine the miRNA expression in Neuroendocrine Tumors (NETs), and identify their targets and importance in NET carcinogenesis. miRNA expression in six NEN primary tumors, six NEN metastases and four normal intestinal tissues was characterized using miRNA arrays, and validated by in-situ hybridization and qPCR. Among the down-regulated miRNAs miR-129-5p and the let-7f/let-7 family, were selected for further characterization. Transfection of miR-129-5p inhibited growth of a pulmonary and an intestinal carcinoid cell line. Analysis of mRNA expression changes identified EGR1 and G3BP1 as miR-129-5p targets. They were validated by luciferase assay and western blotting, and found robustly expressed in NETs by immunohistochemistry. Knockdown of EGR1 and G3BP1 mimicked the growth inhibition induced by miR-129-5p. let-7 overexpression inhibited growth of carcinoid cell lines, and let-7 inhibition increased protein content of the transcription factor BACH1 and its targets MMP1 and HMGA2, all known to promote bone metastases. Immunohistochemistry analysis revealed that let-7 targets are highly expressed in NETs and metastases. We found down-regulation of miR-129-5p and the let-7 family, and identified new neuroendocrine specific targets for these miRNAs, which contributes to the growth and metastatic potential of these tumors.

cWords - systematic microRNA regulatory motif discovery from mRNA expression data.

BACKGROUND: Post-transcriptional regulation of gene expression by small RNAs and RNA binding proteins is of fundamental importance in development of complex organisms, and dysregulation of regulatory RNAs can influence onset, progression and potentially be target for treatment of many diseases. Post-transcriptional regulation by small RNAs is mediated through partial complementary binding to messenger RNAs leaving nucleotide signatures or motifs throughout the entire transcriptome. Computational methods for discovery and analysis of sequence motifs in high-throughput mRNA expression profiling experiments are becoming increasingly important tools for the identification of post-transcriptional regulatory motifs and the inference of the regulators and their targets. RESULTS: cWords is a method designed for regulatory motif discovery in differential case-control mRNA expression datasets. We have improved the algorithms and statistical methods of cWords, resulting in at least a factor 100 speed gain over the previous implementation. On a benchmark dataset of 19 microRNA (miRNA) perturbation experiments cWords showed equal or better performance than two comparable methods, miReduce and Sylamer. We have developed rigorous motif clustering and visualization that accompany the cWords analysis for more intuitive and effective data interpretation. To demonstrate the versatility of cWords we show that it can also be used for identification of potential siRNA off-target binding. Moreover, cWords analysis of an experiment profiling mRNAs bound by Argonaute ribonucleoprotein particles discovered endogenous miRNA binding motifs. CONCLUSIONS: cWords is an unbiased, flexible and easy-to-use tool designed for regulatory motif discovery in differential case-control mRNA expression datasets. cWords is based on rigorous statistical methods that demonstrate comparable or better performance than other existing methods. Rich visualization of results promotes intuitive and efficient interpretation of data. cWords is available as a stand-alone Open Source program at Github https://github.com/simras/cWords and as a web-service at: http://servers.binf.ku.dk/cwords/.

Analysis of microRNA-target interactions across diverse cancer types.

Little is known about the extent to which individual microRNAs (miRNAs) regulate common processes of tumor biology across diverse cancer types. Using molecular profiles of >3,000 tumors from 11 human cancer types in The Cancer Genome Atlas, we systematically analyzed expression of miRNAs and mRNAs across cancer types to infer recurrent cancer-associated miRNA-target relationships. As we expected, the inferred relationships were consistent with sequence-based predictions and published data from miRNA perturbation experiments. Notably, miRNAs with recurrent target relationships were frequently regulated by genetic and epigenetic alterations across the studied cancer types. We also identify new examples of miRNAs that coordinately regulate cancer pathways, including the miR-29 family, which recurrently regulates active DNA demethylation pathway members TET1 and TDG. The online resource http://cancerminer.org allows exploration and prioritization of miRNA-target interactions that potentially regulate tumorigenesis.

Comprehensive analysis of long non-coding RNAs in ovarian cancer reveals global patterns and targeted DNA amplification.

Long non-coding RNAs (lncRNAs) are emerging as potent regulators of cell physiology, and recent studies highlight their role in tumor development. However, while established protein-coding oncogenes and tumor suppressors often display striking patterns of focal DNA copy-number alteration in tumors, similar evidence is largely lacking for lncRNAs. Here, we report on a genomic analysis of GENCODE lncRNAs in high-grade serous ovarian adenocarcinoma, based on The Cancer Genome Atlas (TCGA) molecular profiles. Using genomic copy-number data and deep coverage transcriptome sequencing, we derived dual copy-number and expression data for 10,419 lncRNAs across 407 primary tumors. We describe global correlations between lncRNA copy-number and expression, and associate established expression subtypes with distinct lncRNA signatures. By examining regions of focal copy-number change that lack protein-coding targets, we identified an intergenic lncRNA on chromosome 1, OVAL, that shows narrow focal genomic amplification in a subset of tumors. While weakly expressed in most tumors, focal amplification coincided with strong OVAL transcriptional activation. Screening of 16 other cancer types revealed similar patterns in serous endometrial carcinomas. This shows that intergenic lncRNAs can be specifically targeted by somatic copy-number amplification, suggestive of functional involvement in tumor initiation or progression. Our analysis provides testable hypotheses and paves the way for further study of lncRNAs based on TCGA and other large-scale cancer genomics datasets.

Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal.

The cBioPortal for Cancer Genomics (http://cbioportal.org) provides a Web resource for exploring, visualizing, and analyzing multidimensional cancer genomics data. The portal reduces molecular profiling data from cancer tissues and cell lines into readily understandable genetic, epigenetic, gene expression, and proteomic events. The query interface combined with customized data storage enables researchers to interactively explore genetic alterations across samples, genes, and pathways and, when available in the underlying data, to link these to clinical outcomes. The portal provides graphical summaries of gene-level data from multiple platforms, network visualization and analysis, survival analysis, patient-centric queries, and software programmatic access. The intuitive Web interface of the portal makes complex cancer genomics profiles accessible to researchers and clinicians without requiring bioinformatics expertise, thus facilitating biological discoveries. Here, we provide a practical guide to the analysis and visualization features of the cBioPortal for Cancer Genomics.

An epidemiologic and genomic investigation into the obesity paradox in renal cell carcinoma.

BACKGROUND: Obesity increases risk for clear-cell renal cell carcinoma (ccRCC), yet obese patients appear to experience longer survival than nonobese patients. We examined body mass index (BMI) in relation to stage, grade, and cancer-specific mortality (CSM) while considering detection bias, nutritional status, and molecular tumor features. METHODS: Data were available from 2119 ccRCC patients who underwent renal mass surgery at Memorial Sloan-Kettering Cancer Center between 1995 and 2012. Logistic regression models produced associations between BMI and advanced disease. Multivariable competing risks regression models estimated associations between BMI and CSM. Somatic mutation, copy number, methylation, and expression data were examined by BMI among a subset of 126 patients who participated in the Cancer Genome Atlas Project for ccRCC using the Kruskal-Wallis or Fisher exact tests. All statistical tests were two-sided. RESULTS: Obese and overweight patients were less likely to present with advanced-stage disease compared with normal-weight patients (odds ratio [OR] = 0.61, 95% confidence interval [CI] = 0.48 to 0.79 vs OR = 0.65, 95% CI = 0.51 to 0.83, respectively). Higher BMI was associated with reduced CSM in univariable analyses (P < .005). It remained statistically significant after adjustment for comorbidities and albumin level, but it became non-statistically significant after adjusting for stage and grade (P > .10). Genome-wide interrogation by BMI suggested differences in gene expression of metabolic and fatty acid genes, including fatty acid synthase (FASN), consistent with the obesity paradox. CONCLUSIONS: Our findings suggest that although BMI is not an independent prognostic factor for CSM after controlling for stage and grade, tumors developing in an obesogenic environment may be more indolent.