Identifying mRNA targets of microRNA dysregulated in cancer: with application to clear cell Renal Cell Carcinoma

BACKGROUND. MicroRNA regulate mRNA levels in a tissue specific way, either by inducing degradation of the transcript or by inhibiting translation or transcription. Putative mRNA targets of microRNA identified from seed sequence matches are available in many databases. However, such matches have a high false positive rate and cannot identify tissue specificity of regulation. RESULTS. We describe a simple method to identify direct mRNA targets of microRNA dysregulated in cancers from expression level measurements in patient matched tumor/normal samples. The word "direct" is used here in a strict sense to: a) represent mRNA which have an exact seed sequence match to the microRNA in their 3'UTR, b) the seed sequence match is strictly conserved across mouse, human, rat and dog genomes, c) the mRNA and microRNA expression levels can distinguish tumor from normal with high significance and d) the microRNA/mRNA expression levels are strongly and significantly anti-correlated in tumor and/or normal samples. We apply and validate the method using clear cell Renal Cell Carcinoma (ccRCC) and matched normal kidney samples, limiting our analysis to mRNA targets which undergo degradation of the mRNA transcript because of a perfect seed sequence match. Dysregulated microRNA and mRNA are first identified by comparing their expression levels in tumor vs normal samples. Putative dysregulated microRNA/mRNA pairs are identified from these using seed sequence matches, requiring that the seed sequence be conserved in human/dog/rat/mouse genomes. These are further pruned by requiring a strong anti-correlation signature in tumor and/or normal samples. The method revealed many new regulations in ccRCC. For instance, loss of miR-149, miR-200c and mir-141 causes gain of function of oncogenes (KCNMA1, LOX), VEGFA and SEMA6A respectively and increased levels of miR-142-3p, miR-185, mir-34a, miR-224, miR-21 cause loss of function of tumor suppressors LRRC2, PTPN13, SFRP1, ERBB4, and (SLC12A1, TCF21) respectively. We also found strong anti-correlation between VEGFA and the miR-200 family of microRNA: miR-200a*, 200b, 200c and miR-141. Several identified microRNA/mRNA pairs were validated on an independent set of matched ccRCC/normal samples. The regulation of SEMA6A by miR-141 was verified by a transfection assay. CONCLUSIONS. We describe a simple and reliable method to identify direct gene targets of microRNA in any cancer. The constraints we impose (strong dysregulation signature for microRNA and mRNA levels between tumor/normal samples, evolutionary conservation of seed sequence and strong anti-correlation of expression levels) remove spurious matches and identify a subset of robust, tissue specific, functional mRNA targets of dysregulated microRNA.Cancer Institute of New Jersy; New Jersey Commission for Cacner Research; Lineberger Comprehensive Cancer Center Tissue Procurement and Genomics Core Facility; Crawford Fun

Previously Unidentified Changes in Renal Cell Carcinoma Gene Expression Identified by Parametric Analysis of Microarray Data

BACKGROUND. Renal cell carcinoma is a common malignancy that often presents as a metastatic-disease for which there are no effective treatments. To gain insights into the mechanism of renal cell carcinogenesis, a number of genome-wide expression profiling studies have been performed. Surprisingly, there is very poor agreement among these studies as to which genes are differentially regulated. To better understand this lack of agreement we profiled renal cell tumor gene expression using genome-wide microarrays (45,000 probe sets) and compare our analysis to previous microarray studies. METHODS. We hybridized total RNA isolated from renal cell tumors and adjacent normal tissue to Affymetrix U133A and U133B arrays. We removed samples with technical defects and removed probesets that failed to exhibit sequence-specific hybridization in any of the samples. We detected differential gene expression in the resulting dataset with parametric methods and identified keywords that are overrepresented in the differentially expressed genes with the Fisher-exact test. RESULTS. We identify 1,234 genes that are more than three-fold changed in renal tumors by t-test, 800 of which have not been previously reported to be altered in renal cell tumors. Of the only 37 genes that have been identified as being differentially expressed in three or more of five previous microarray studies of renal tumor gene expression, our analysis finds 33 of these genes (89%). A key to the sensitivity and power of our analysis is filtering out defective samples and genes that are not reliably detected. CONCLUSIONS. The widespread use of sample-wise voting schemes for detecting differential expression that do not control for false positives likely account for the poor overlap among previous studies. Among the many genes we identified using parametric methods that were not previously reported as being differentially expressed in renal cell tumors are several oncogenes and tumor suppressor genes that likely play important roles in renal cell carcinogenesis. This highlights the need for rigorous statistical approaches in microarray studies.National Institutes of Healt

Up-regulation of microRNA-21 correlates with lower kidney cancer survival.

BackgroundMicroRNA-21 is up-regulated in a variety of cancers like, breast, colorectal, lung, head and neck etc. However, the regulation of miR-21 in renal cell carcinoma (RCC) has not yet been studied systematically.Methods and resultsWe measured miR-21 levels in 54 pairs of kidney cancers and their normal matched tissues by real-time PCR. The expression level of miR-21 was correlated with 5 year survival and the pathological stage. Functional studies were done after inhibiting miR-21 in RCC cell lines. We studied in vitro and in vivo effects of the chemo preventive agent genistein on miR-21 expression. In 48 cases (90%), miR-21 was increased. All patients with low miR-21 expression survived 5 years, while with high miR-21 expression, only 50% survived. Higher expression of miR-21 is associated with an increase in the stage of renal cancer. Functional studies after inhibiting miRNA-21 in RCC cell lines show cell cycle arrest, induction of apoptosis and reduced invasive and migratory capabilities. Western blot analysis showed an increase in the expression of p21 and p38 MAP kinase genes and a reduction in cyclin E2. Genistein inhibited the expression of miR-21 in A-498 cells and in the tumors formed after injecting genistein treated A-498 cells in nude mice besides inhibiting tumor formation.ConclusionsThe current study shows a clear correlation between miR-21 expression and clinical characteristics of renal cancer. Thus we believe that miR-21 can be used as a tumor marker and its inhibition may prove to be useful in controlling cancers with up-regulated miR-21

Glutathione Metabolism in Renal Cell Carcinoma Progression and Implications for Therapies

A significantly increased level of the reactive oxygen species (ROS) scavenger glutathione (GSH) has been identified as a hallmark of renal cell carcinoma (RCC). The proposed mechanism for increased GSH levels is to counteract damaging ROS to sustain the viability and growth of the malignancy. Here, we review the current knowledge about the three main RCC subtypes, namely clear cell RCC (ccRCC), papillary RCC (pRCC), and chromophobe RCC (chRCC), at the genetic, transcript, protein, and metabolite level and highlight their mutual influence on GSH metabolism. A further discussion addresses the question of how the manipulation of GSH levels can be exploited as a potential treatment strategy for RCC

Global isoform-specific transcript alterations and deregulated networks in clear cell renal cell carcinoma.

Extensive genome-wide analyses of deregulated gene expression have now been performed for many types of cancer. However, most studies have focused on deregulation at the gene-level, which may overlook the alterations of specific transcripts for a given gene. Clear cell renal cell carcinoma (ccRCC) is one of the best-characterized and most pervasive renal cancers, and ccRCCs are well-documented to have aberrant RNA processing. In the present study, we examine the extent of aberrant isoform-specific RNA expression by reporting a comprehensive transcript-level analysis, using the new kallisto-sleuth-RATs pipeline, investigating coding and non-coding differential transcript expression in ccRCC. We analyzed 50 ccRCC tumors and their matched normal samples from The Cancer Genome Altas datasets. We identified 7,339 differentially expressed transcripts and 94 genes exhibiting differential transcript isoform usage in ccRCC. Additionally, transcript-level coexpression network analyses identified vasculature development and the tricarboxylic acid cycle as the most significantly deregulated networks correlating with ccRCC progression. These analyses uncovered several uncharacterized transcripts, including lncRNAs FGD5-AS1 and AL035661.1, as potential regulators of the tricarboxylic acid cycle associated with ccRCC progression. As ccRCC still presents treatment challenges, our results provide a new resource of potential therapeutics targets and highlight the importance of exploring alternative methodologies in transcriptome-wide studies

SWIM: A computational tool to unveiling crucial nodes in complex biological networks

SWItchMiner (SWIM) is a wizard-like software implementation of a procedure, previously described, able to extract information contained in complex networks. Specifically, SWIM allows unearthing the existence of a new class of hubs, called "fight-club hubs", characterized by a marked negative correlation with their first nearest neighbors. Among them, a special subset of genes, called "switch genes", appears to be characterized by an unusual pattern of intra- and inter-module connections that confers them a crucial topological role, interestingly mirrored by the evidence of their clinic-biological relevance. Here, we applied SWIM to a large panel of cancer datasets from The Cancer Genome Atlas, in order to highlight switch genes that could be critically associated with the drastic changes in the physiological state of cells or tissues induced by the cancer development. We discovered that switch genes are found in all cancers we studied and they encompass protein coding genes and non-coding RNAs, recovering many known key cancer players but also many new potential biomarkers not yet characterized in cancer context. Furthermore, SWIM is amenable to detect switch genes in different organisms and cell conditions, with the potential to uncover important players in biologically relevant scenarios, including but not limited to human cancer

CRISPR-Mediated VHL Knockout Generates an Improved Model for Metastatic Renal Cell Carcinoma.

Metastatic renal cell carcinoma (mRCC) is nearly incurable and accounts for most of the mortality associated with RCC. Von Hippel Lindau (VHL) is a tumour suppressor that is lost in the majority of clear cell RCC (ccRCC) cases. Its role in regulating hypoxia-inducible factors-1α (HIF-1α) and -2α (HIF-2α) is well-studied. Recent work has demonstrated that VHL knock down induces an epithelial-mesenchymal transition (EMT) phenotype. In this study we showed that a CRISPR/Cas9-mediated knock out of VHL in the RENCA model leads to morphologic and molecular changes indicative of EMT, which in turn drives increased metastasis to the lungs. RENCA cells deficient in HIF-1α failed to undergo EMT changes upon VHL knockout. RNA-seq revealed several HIF-1α-regulated genes that are upregulated in our VHL knockout cells and whose overexpression signifies an aggressive form of ccRCC in the cancer genome atlas (TCGA) database. Independent validation in a new clinical dataset confirms the upregulation of these genes in ccRCC samples compared to adjacent normal tissue. Our findings indicate that loss of VHL could be driving tumour cell dissemination through stabilization of HIF-1α in RCC. A better understanding of the mechanisms involved in this phenomenon can guide the search for more effective treatments to combat mRCC

Expression of LDL receptor-related proteins (LRPs) in common solid malignancies correlates with patient survival

LDL receptor-related proteins (LRPs) are transmembrane receptors involved in endocytosis, cell-signaling, and trafficking of other cellular proteins. Considerable work has focused on LRPs in the fields of vascular biology and neurobiology. How these receptors affect cancer progression in humans remains largely unknown. Herein, we mined provisional data-bases in The Cancer Genome Atlas (TCGA) to compare expression of thirteen LRPs in ten common solid malignancies in patients. Our first goal was to determine the abundance of LRP mRNAs in each type of cancer. Our second goal was to determine whether expression of LRPs is associated with improved or worsened patient survival. In total, data from 4,629 patients were mined. In nine of ten cancers studied, the most abundantly expressed LRP was LRP1; however, a correlation between LRP1 mRNA expression and patient survival was observed only in bladder urothelial carcinoma. In this malignancy, high levels of LRP1 mRNA were associated with worsened patient survival. High levels of LDL receptor (LDLR) mRNA were associated with decreased patient survival in pancreatic adenocarcinoma. High levels of LRP10 mRNA were associated with decreased patient survival in hepatocellular carcinoma, lung adenocarcinoma, and pancreatic adenocarcinoma. LRP2 was the only LRP for which high levels of mRNA expression correlated with improved patient survival. This correlation was observed in renal clear cell carcinoma. Insights into LRP gene expression in human cancers and their effects on patient survival should guide future research