BACKGROUND
microrna regulate mrna and protein levels by cleavage and/or translation/transcriptional repression in a tissue specific manner  <cit> .
by modulating key cellular processes such as metabolism, division, differentiation, development and apoptosis, they can simultaneously regulate both oncogenes and tumor suppressor genes  <cit> .
aberrant microrna profiles have been noted in many cancers  <cit> , including renal cell carcinoma  <cit> .
almost half the known micrornas are in cancer-associated chromosomal fragile sites, susceptible to point mutation, amplification, deletion, or translocation  <cit> .
recent evidence demonstrates that microrna play an important role in the patho-physiology of many cancers  <cit>  and they are believed to be involved in pathogenesis in ccrcc  <cit> .
microrna are also being studied in various tumors to understand their significance for drug resistance  <cit> , diagnosis and prognosis  <cit>  and for their therapeutic potential  <cit> .
their secondary structure preserves them better in ffpe samples than mrna, making them easier to extract in intact form, resulting in higher identification accuracy in the analysis of archived clinical material  <cit> .
their tissue specificity and tight regulation makes them more reliable identifiers of tissue of origin in highly differentiated tumors  <cit> .
single microrna can regulate multiple mrna and are therefore both better identifiers of mechanism and possibly better drug targets  <cit> .
however, while it is clear that microrna play an important role in the biology of many cancers, their complex biology and tissue specificity makes it difficult to understand the precise role they play in the disease process and the genes affected by their dysregulation  <cit> .
mature micrornas are produced in a multi-stage process.
after transcription, they are processed by rna pol ii or pol iii to create capped and polyadenylated primary transcripts , which are further processed in the nucleus by the enzyme drosha/pasha  or by dgcr <dig>  to produce ~60-nucleotide pre-microrna stem-loop molecules.
these are then exported to the cytoplasm by exportin and ran-gtp where they are further processed by dicer to ~ <dig> nt double-stranded rna duplexes, which form complexes with risc  leading to unwinding of the duplexes to form single-stranded micrornas.
micrornas bound to risc can down-regulate protein levels using at least two alternative pathways: 1) if the microrna has imperfect complementarity with a matching sequence in the 3'utr of its target mrna, the microrna-risc complex can combine with the complement mrna sequence and cause translational repression.
2) on the other hand, if the microrna and its mrna target have perfect or near perfect complementarity, the microrna-risc complex binding to its target mrnas can result in the cleavage and degradation of the mrna by argonaute <dig>   <cit> .
although many studies have identified signatures of microrna dysregulation, the identification of tissue specific targets of aberrantly regulated microrna is difficult.
putative identification using seed sequence complementarity and free energy predictions of rna-rna duplexes  <cit>  are available in databases such as targetscan: http://www.targetscan.org.
however, the false positive rate for such matches is unacceptably high, with different algorithms identifying different mrna targets for the same microrna  <cit> .
the tissue specificity of microrna regulation is known only in some specific cases  and a general methodology for target identification, tissue specificity of action and specific biological role of microrna in the initiation and progression of most cancers remains an open problem.
we describe a novel method to identify "direct mrna targets" of microrna in any cancer based on measuring an anti-correlation signal between differentially expressed microrna and mrna in patient matched tumor and normal samples.
in this paper, the words "direct mrna targets" is used in a very strict and limited sense.
a direct target is one which: a) has an exact seed sequence match in its 3'utr to the corresponding microrna, b) the seed sequence match is conserved across mouse, human, rat and dog genomes, c) the expression levels of both the microrna and the mrna can distinguish tumor from normal with high statistical significance and d) the mrna and microrna levels are strongly and significantly anti-correlated in tumor and/or normal.
these requirements could be relaxed to find additional targets or eliminated altogether to find indirect regulations .
the method proceeds as follows: a) identify significantly differentially expressed microrna and mrna between the two classes ; b) for each microrna which is differentially expressed, identify all its putative target mrna by restricting to those differentially expressed mrna with a matching seed sequence in their 3'utr, with the further requirement that it be conserved in human, mouse, rat and dog genomes; c) compute the pearson correlation between microrna and mrna expression levels for samples in each class  and d) retain only those microrna/mrna pairs whose expression levels are highly anti-correlated.
these constraints remove spurious matches, reducing relatively speculative "putative" seed match based mrna targets in databases to a highly robust subset of direct functional targets.
note that our method can be extended  by removing constraint b) and looking for a correlation  signature in c).
this allows the identification of indirect regulation.
for example, if a microrna up-regulated in cancer down-regulates a gene which is a transcriptional repressor of an oncogene, then the expression level of the microrna will be correlated with the level of the oncogene without a seed sequence match.
since the direct gene target  of the microrna should already be identified using our method, such an analysis would extend the regulation network beyond first order interactions.
note that, although the method as described above does not identify regulation by translation inhibition , if protein levels were also measured, the method could easily be extended to identify such regulation.
we demonstrate the use of our method on expression data from clear cell renal cell carcinoma  and matched normal kidney samples.
renal cell carcinoma  represents ~3% of all malignancies in the us, with  <dig>  new cases and  <dig>  deaths each year http://www.nci.nih.gov/cancertopics/types/kidney.
the most common histological class is ccrcc, accounting for ~75% of kidney cancers.
ccrcc is known to be characterized by the loss of the vhl gene, which under normal oxygen pressure, binds to the α subunits of hypoxia-inducible factors , inducing their poly-ubiquitinylation and subsequent degradation in the proteasome.
in hypoxic conditions, or if hif regulation is lost because of vhl inactivation, hif accumulates to high levels and promotes the transcription of genes such as vegf, pdgf-β, tgf-α, epo etc which trigger angiogenesis, cell growth, migration and proliferation  <cit> .
the spectrum of hif target genes expressed in individual tumors and the factors which influence them are the object of active ongoing research.
ccrcc tumors have a wide range of natural histories and varied responses to vegf-targeted therapy  <cit> .
early stage, fuhrman grade  <dig>  tumors tend to have significantly better disease free survival after resection than higher stage and grade   <cit> .
although vhl mutation is associated with all grades of ccrcc, the other molecular factors associated with ccrcc initiation and progression are largely unknown.
the molecular basis of the diversity in histologic grade, clinical behavior, and response to vegf-targeted is also unclear, and makes ccrcc a ripe target for studies investigating the molecular and genetic nature of these heterogeneities.
in rcc, various studies have identified panels of microrna and mrna that are differentially expressed between normal renal tissue and tumor or between histological subtypes of tumor  <cit> .
the present study extends these previous studies by linking the microrna to some of their mrna targets, thus elucidating a hitherto unknown part of the biology of ccrcc disease.
some of the identified microrna/mrna anti-correlations were validated on a new cohort of ccrcc/normal samples.
sema6a was confirmed as a direct target of mir- <dig> by over-expressing mir- <dig> in a ccrcc cell line and showing strong down-regulation of the sema6a transcript.
RESULTS
the underlying hypothesis in our method is that the expression levels of microrna and their direct mrna targets should be strongly anti-correlated when averaged over matched samples in either tumor or normal tissue.
the stepwise procedure is as follows:
step 1: identify significantly up/down regulated micrornas in ccrcc samples vs normal samples.
step 2: identify significantly up/down regulated mrnas in ccrcc samples vs normal samples.
step 3: using targetscan, retain only the mrna in step  <dig> which have a conserved seed sequence in their 3' utr for at least one of the microrna from step  <dig>
step 4: find anti-correlated pairs of up-regulated microrna and down-regulated mrna in ccrcc samples using a strict cutoff  in pearson correlation coefficient in ccrcc samples.
similarly, find anti-correlated pairs of down-regulated microrna and up-regulated mrnas in normal kidney samples using a strict cutoff.
in step  <dig>  putative target mrna were identified using targetscan version  <dig>  http://www.targetscan.org, which identifies possible regulatory targets of mammalian micrornas as those with conserved sequences of matching seed regions for each microrna.
the term "conserved" means that the sequence is conserved in human, mouse, rat and dog.
in step  <dig>  we ran  <dig> permutations in brb-arraytools' http://linus.nci.nih.gov/brb-arraytools.html and multivariate/univariate analysis at p <  <dig> , fdr <  <dig>  to assess significance of discovered pairs and to find the appropriate cutoff p <dig> for significance of the measured pearson correlation.
in the primary dataset, because of the high accuracy of qrt-pcr, we were able to set a strict cutoff p <dig> = - <dig> ; thus in step  <dig> only microrna/mrna pairs with p < p <dig>  were considered to be significant.
step 1: microrna significantly differentially expressed in ccrcc versus normal kidney tissue
<dig> microrna were identified as differentially expressed  in ccrcc versus normal kidney,  <dig> down-regulated and  <dig> up-regulated.
the microrna down-regulated in ccrcc were mir- <dig>  mir-10b, mir-125b, mir-26a+, mir-133b, mir-135a, mir-135b, mir- <dig>  mir- <dig>  mir- <dig>  mir- <dig>  mir-199a, mir-200a, mir-200b, mir-200c, mir- <dig>  mir- <dig>  mir- <dig>  mir-30a-3p, mir-30a-5p, mir- <dig>  mir- <dig>  mir- <dig>  mir- <dig>  mir- <dig>  mir- <dig> and the microrna up-regulated in ccrcc were mir-142-3p, mir- <dig>  mir- <dig>  mir- <dig>  mir- <dig>  mir- <dig>  mir-34a, mir-34b, mir- <dig> , in agreement with recent studies  <cit> .
the chromosomal location of these microrna  shows that several are proximal to fragile regions and regions commonly deleted in ccrcc  <cit> .
figure 1a shows a heat-map of the microrna expression levels for samples in the primary dataset showing that they can robustly separate ccrcc from normal kidney.
cancer-related regions listed are either deleted regions , amplified regions  or breakpoint regions.
step 2: mrna differentially expressed in ccrcc versus normal kidney tissue
using sam http://www-stat.stanford.edu/~tibs/sam/ at fdr =  <dig>  and  <dig> permutation tests,  <dig> mrna probes were found up-regulated and  <dig> were found down-regulated in ccrcc.
the genes represented by these probes included many genes known to be dysregulated in ccrcc  <cit> .
we saw significant upregulation of oncogenes vegfa  and egfr , integrin-mediated cell adhesion pathway genes , cell adhesion related genes  and cytokines .
similarly we noted down-regulation of tumor suppressor genes  and of members of the metallothionein family .
steps  <dig> and 4: identifying direct mrna targets of dysregulated microrna
pearson correlation analysis with p <dig> = - <dig>  was applied to each of the  <dig> differentially expressed microrna and its putative mrna targets .
this procedure identified  <dig> mrna targets for the  <dig> up-regulated microrna.
this list included several important tumor suppressor genes, such as pten, erbb <dig> and sfrp <dig>  known to be mutated or down-regulated in many tumors, including ccrcc  <cit> .
the  <dig> down-regulated microrna had  <dig> direct up-regulated mrna targets, including oncogenes vegfa, lox, loxl <dig> and fas, well known to be involved in kidney cancer  <cit> .
the nine most significantly down-regulated and twelve most significantly up-regulated mrna are listed in table  <dig> and their heatmap is shown in figure 1b.
in figure  <dig> we plot microrna and mrna levels for mir-200c and its target vegfa.
note that the levels of mir-200c and its target vegfa are not only anti-correlated overall, but are also anti-correlated separately in both ccrcc and normal tissue.
additional files  <dig> and  <dig> contain the full list of microrna/mrna regulations identified by our analysis.
the full list of identified microrna/mrna pairs is given in additional files  <dig> and  <dig>  + marks microrna/mrna anti-correlations that were tested/validated on an independent set of matched ccrcc/normal kidney samples.
validation on a cohort of ccrcc/normal kidney samples
a. validation of differentially expressed microrna and mrna
seventeen validation samples were collected from a new cohort of patients and analyzed by microrna and mrna profiling on agilent arrays .
figure 1c shows a heat map of twenty-seven microrna  of the thirty five differentially expressed microrna identified previously.
a weighted voting classifier on binarized microrna expression data in the validation set had 100% accuracy in leave-one-out  cross validation in distinguishing ccrcc from normal kidney.
figure 1d shows a heatmap of mrna expression levels for the genes in figure 1b.
these were also found to be 100% accurate at discriminating ccrcc from normal kidney using weighted voting and loo cross validation analysis.
mrna levels of three genes  which were down-regulated in ccrcc and one gene  up-regulated in ccrcc were also measured by quantitative rt-pcr in twelve of the test samples.
the results, shown in figure 3a, demonstrate that mrna levels of erbb <dig>  sfrp <dig>  slc12a <dig> and vegfa were quantitatively and significantly down/up-regulated as expected.
b. validation of anti-correlation signature between some identified microrna/mrna pairs
in figure 3b-e, we plot the qrt-pcr expression levels of erbb <dig>  sfrp <dig>  slc12a <dig> and vegfa versus agilent chip measured levels of their regulatory microrna  for the twelve samples of figure 3a.
the overall strong anti-correlation signature between microrna and mrna levels is clearly visible in these plots.
figure  <dig> summarizes our validation analysis of a number of anti-correlation measurements between several identified microrna/mrna pairs in the agilent chip data.
the measured correlations between predicted microrna/mrna pairs are shown in the figure.
because of the higher level of noise in the agilent chip data compared to qrt-pcr, we cannot apply the strict criterion  used in the primary dataset.
instead, the significance of the correlation  was computed using permutation tests as follows: a large number of datasets were obtained by permuting the sample labels in the microrna or mrna measurements.
for each microrna/mrna pair, these permuted datasets were used to compute the null distribution for p and the significance of the measured value of p was estimated in this null distribution.
as the p-values in figure  <dig> indicate, we validate a strong anti-correlation signature between mrna levels of , vegf, sema6a, , sfrp <dig>  erbb <dig>  slc12a <dig> and tcf <dig>  and their identified regulators: mir- <dig>  mir-200c, mir- <dig>  mir-142-3p, mir- <dig>  mir-34a, mir- <dig> and mir- <dig> respectively.
c. in vitro validation of sema6a as a target of mir- <dig> in an rcc cell line
finally, to establish that this method can accurately predict functional and direct microrna/mrna regulation, we performed an in vitro analysis of one microrna , and its identified direct target sema6a.
the rcc cell line crl- <dig> was transfected  with either pre-mir- <dig> or a control pre-mir, and levels of sema6a were measured on the case/control cell lines by semi-quantitative rt-pcr.
the results  showed that introduction of pre-mir- <dig> produced a significant reduction in the level of sema6a mrna, validating sema6a as a functional and direct target of mir- <dig>
discussion
we have developed and validated a simple method to identify direct functional mrna targets of microrna in ccrcc using patient matched tumor/normal samples.
the method looks for the subset of anti-correlated microrna/mrna pairs from a larger set of microrna and mrna able to distinguish tumor from normal tissue, with the additional requirement of a highly conserved seed sequences for the corresponding microrna in the 3'utr of the corresponding gene.
our method can not only identify microrna/mrna pairs which discriminate normal from cancer tissue, but also dysregulated regulation mechanisms between them that may initiate and/or drive the disease process.
we used an rt-pcr panel for our initial discovery experiments and microrna and microarray chips both for validation of our previously discovered microrna/mrna pairings and to expand the available pool of micrornas to all currently known micrornas for a comprehensive examination of microrna/mrna pairs.
in comparing the two data modalities, we noticed that in the rt-pcr data, the sample to sample variation in the normal samples  was comparable to the change in level between normal and tumor  which caused a statistically significant up/down regulation of the target mrna.
to find reliable matches, it was therefore crucial to minimize stochastic error.
we observed that the rt-pcr data had lower stochastic variance than the microarray data .
this noise effect is also reflected in the measured p <dig> values in the validation analysis on the microarray data , which are smaller in magnitude  than the strict cutoff |p0| >  <dig>  used in the discovery phase on the rt-pcr data.
since the present study is mainly to demonstrate "proof of concept", we limited the rt-pcr data for discovery and the microarray data for validation.
in a more extensive study, with sufficiently large sample sizes and more accurate data from high throughput technologies  it may be feasible and cost effective to use a high throughput platform for discovery and rt-pcr for validation.
such an analysis might also identify a much bigger set of microrna/mrna relationships.
although our method is very robust, as described, it cannot find all mrna targets.
its principal merit is the ability to reduce the large number of relatively speculative matches from seed sequences alone to a smaller set of functional, tissue specific targets.
while this smaller set is perhaps incomplete, it is more reliable.
additionally, the use of perfect seed sequence matches can easily be relaxed.
a more serious issue is that the method as described cannot find targets affected by translation inhibition.
for such targets, the changes in microrna levels would leave mrna levels unaffected, but reduce protein levels.
such targets might be identified if protein levels in tumor/normal samples were also measured along with the microrna and mrna levels.
the method could then be adapted to identify targets affected by translation inhibition by using an appropriate joint correlation/anti-correlation analysis of microrna/mrna/protein levels.
our method, as described here, would also miss regulation that proceeds via intermediate states.
in these cases, there would be no seed sequence match  between a microrna and the mrna whose level was affected by it.
we could attempt to identify such secondary effects by eliminating the seed sequence match requirement but keeping the requirement of a high correlation/anti-correlation signal between microrna/mrna levels.
whether such a procedure would work would depend on the degree to which the intermediate state gene/protein is affected and measurable.
we briefly discuss some of these issues below.
however, it should be noted that these studies would require significantly larger sample sizes and are beyond the scope of the present paper.
the next round of tcga http://cancergenome.nih.gov/ may make these types of analyses more feasible.
although here, we focused on microrna/mrna pairs with exact seed sequence match and an anti-correlation signal in both tumor and normal samples, as noted above, with additional samples it should be possible to use this method and simple extensions to identify more subtle types of dysregulation.
for example, if a mutation in a microrna in the tumor samples causes loss of its function , then although the microrna/mrna levels would be anti-correlated in normal tissue, they would not be anti-correlated in tumor samples.
similarly a de-novo gain of microrna function would be signaled by an anti-correlated signal in tumor samples which was absent in the normal samples.
in our data, we did observe several correlations and anti-correlations between microrna/mrna pairs in tumor or normal samples without a corresponding seed sequence match.
as noted above, these most likely represent regulations which proceed via intermediate states and hence cannot be identified by seed sequence matches.
for example, if increased levels of microrna x down-regulates expression of a protein which is a transcriptional repressor of gene y, there will be a strong correlation between x and y levels but no seed sequence match.
conversely, if microrna x regulates a protein which is a transcription factor for gene y then the levels of x and y will be anti-correlated without a seed sequence match.
the measurement of such correlations would extend the network of microrna control beyond first level regulators but would require significant increases in the number of samples  for statistical significance.
since the method finds functional relationships, it should be useful for identifying pharmaceutically relevant mechanisms which suggest drug targets for therapy.
we describe below some of the regulations we identified which have pharmaceutical relevance.
microrna which shut down multiple tumor suppressor genes in ccrcc
we found several dysregulated microrna which targeted multiple tumor suppressor genes.
for example, the oncogenic mir- <dig> was significantly up-regulated in ccrcc and anti-correlated with the tumor suppressor gene pten, suggesting that its gain of function shuts down pten in ccrcc.
mutated or down-regulated in many advanced cancers  <cit> , pten loss activates the pi3k-akt  <cit>  signaling pathway and its downstream target mtor, with important implications in rcc development and therapeutic selection  <cit> .
another identified target of mir- <dig> was ptpn <dig> , which is a fas-associated protein tyrosine phosphatase and putative tumor suppressor gene that can inhibit pi3k/akt signaling, suppress the influence of insulin-like growth factor-i on cell survival and induce apoptosis  <cit> .
kcnj <dig>  member of the potassium channel subfamily of membrane proteins, was also identified as a target of mir- <dig>  such membrane proteins have been suggested in anti-cancer therapies because of their important role in cell growth  <cit>  and are known to be down-regulated in rcc  <cit> .
mir-34a is known to be over-expressed in various tumors and associated with cell proliferation  <cit> .
in our data, it was up-regulated in ccrcc and predicted to target sfrp <dig>  a known regulator of the wnt signaling pathway and a tumor suppressor gene whose loss has been observed in a majority of rcc patients  <cit> .
another tumor suppressor regulated by mir-34a was camta <dig>  <cit> , a reduction in whose levels correlates with poor outcome in neuroblastoma  <cit> .
finally, we find mir- <dig>  associated with a human chromosome fragile site on chr-xq <dig>  <cit> , is up-regulated in ccrcc and predicted to target erbb <dig>  a member of the egfr family, a potential tumor suppressor known to be strongly down-regulated in ccrcc  <cit> .
hypoxia induced micrornas
a hypoxic tumor microenvironment can directly activate the expression of several microrna  <cit> .
for example, mir- <dig>  mir- <dig> and mir- <dig> reduce pro-apoptotic signaling in response to a hypoxic environment and are consistently over-expressed in a variety of human tumors .
ccrcc is a unique setting in which to study these micrornas, given that vhl loss constitutively stabilizes one or more hif factors, thereby creating a pseudo-hypoxic scenario in ccrcc tumor cells  <cit> .
we found that hypoxia related microrna had the most significant fold changes, with mir- <dig>  mir- <dig> and mir- <dig> being amongst the top., suggesting a major role for them in renal carcinogenesis.
we found  <dig> genes down-regulated in tumor, inversely correlated with mir- <dig> and enriched  in the cell adhesion  pathway at p =  <dig> .
loss of cam degrades the intra and extracellular matrix, leading to abnormal cell growth patterns.
among these  <dig> genes, two were identified as direct targets of mir-21: slc12a <dig> and tcf <dig>  both of which have been reported down-regulated in ccrcc  <cit> .
identification of microrna family interactions
we found evidence that families of micrornas may be coordinately participating in microrna/mrna interactions.
one example in ccrcc is the mir-204/ <dig> family, which was significantly down-regulated in ccrcc samples.
we identified thirty five significant mrna targets  for this family .
among these, eight were on chromosome 3q, a common amplicon region in many epithelial tumors  <cit> .
these genes include c3orf <dig>  ccdc <dig>  dtx3l, pld <dig>  trim <dig>  two oncogenes ect <dig> and rap2b, and a hypoxia associated protein serp <dig>  <cit> .
the gain of this chromosome arm was previously associated with papillary rccs  <cit> , and our observation in clear cell rcc implies a possible regulatory relation between mir-204/ <dig> and the genes in this region as an alternate mechanism of up-regulation of this group of genes.
another example is the mir- <dig> family which includes two microrna clusters, one on chromosome 1p <dig>   and another on chromosome 12p <dig> .
five mir- <dig> family members contain very similar seed sequences - aauacu for mir-200b/200c/ <dig> and aacacu for mir-200a*/ <dig>  <cit> .
recently, several other groups have reported a role for the mir- <dig> family in the epithelial-mesenchymal-transition  and in cancer cell migration, the latter by directly targeting the transcription factors zeb <dig> and zeb <dig>  which regulate e-cadherin, a mediator of cell-cell adhesion  <cit> .
another study  <cit>  identified a regulatory loop between these micrornas and zeb transcription factors as well as the emt inducer tgfβ.
in epithelial cells, mir- <dig> family micrornas and e-cadherin maintain higher level expression by repressing zeb <dig>  zeb <dig> and tgfβ; on the other hand, in mesenchymal cells and tumors, the up-regulation of zeb factors is triggered by tgfβ and suppresses the transcription of mir-141/200c by binding to their putative common promoter region.
in our primary dataset, zeb <dig> and zeb <dig> were both up-regulated in six out of our eight ccrcc samples and, their expression levels were highly anti-correlated with the mir- <dig> family in both tumor and normal samples.
as confirmation of these results, down-regulation of mir- <dig> and mir-200c and their function on zeb <dig> in ccrcc has recently been reported  <cit> .
we also noted that in our data, the anti-correlation between vegfa and the mir- <dig> family was strongest in normal kidney tissue, suggesting that loss of this regulation may be an important factor providing a permissive environment for hif transcriptional signaling.
our hypothesis  from these various observations is that in normal kidney, the expression level of hif2α and its downstream targets  are regulated by mir- <dig>  200a*, 200b and 200c and the loss of this microrna regulation, in concert with vhl loss, is responsible for activation of the hif pathway.
one intriguing association which we have identified  is highly significant for therapy in ccrcc.
this is because the soluble extracellular domain of sema6a has been engineered to effectively inhibit vegf-mediated tumor formation  <cit> .
hence, our results imply that mir- <dig> may have a role in gene therapy.
a model which summarizes our observations and integrates it with mechanisms for ccrcc dysregulation from the literature is shown in figure  <dig>  this model integrates our measured microrna regulatory mechanisms with known transcriptional activity resulting from vhl loss and activation of the hypoxia response pathway.
the pathways shown in figure  <dig> are likely highly interconnecting, and this model, and the individual functional interactions it suggests, need to be validated  by direct experimental targeting of microrna levels and measurements of each of the predicted target genes in a larger cohort  in its next phase).
CONCLUSIONS
the main innovation in this paper is the use of an anti-correlation analysis of microrna/mrna levels in two cell types from the same patient  to identify functional mrna targets of altered microrna.
the method can obviously extend to any tissue type and might be useful in other contexts: for instance, by using microdissection to harvest cells from different cellular compartments from the same breast cancer patient, it might reveal the microrna/mrna regulation program which causes progression of normal epithelium to hyperplasia to dcis to invasive disease.
the method might also be useful in non-cancer studies, such as in identifying the role of microrna and their target genes in the transition from stem cells to differentiated cells or in embryogenesis.
in summary, in this paper, we have demonstrated a simple method to identify tissue specific mrna targets of microrna, which is expandable to many study types.

 microrna regulate mrna and protein levels by cleavage and/or translation/transcriptional repression in a tissue specific manner. by modulating key cellular processes such as metabolism, division, differentiation, development and apoptosis, they can simultaneously regulate both oncogenes and tumor suppressor genes. aberrant microrna profiles have been noted in many cancers, including renal cell carcinoma. almost half the known micrornas are in cancer-associated chromosomal fragile sites, susceptible to point mutation, amplification, deletion, or translocation. recent evidence demonstrates that microrna play an important role in the patho-physiology of many cancers and they are believed to be involved in pathogenesis in ccrcc. microrna are also being studied in various tumors to understand their significance for drug resistance, diagnosis and prognosis and for their therapeutic potential. their secondary structure preserves them better in ffpe samples than mrna, making them easier to extract in intact form, resulting in higher identification accuracy in the analysis of archived clinical material. their tissue specificity and tight regulation makes them more reliable identifiers of tissue of origin in highly differentiated tumors. single microrna can regulate multiple mrna and are therefore both better identifiers of mechanism and possibly better drug targets. however, while it is clear that microrna play an important role in the biology of many cancers, their complex biology and tissue specificity makes it difficult to understand the precise role they play in the disease process and the genes affected by their dysregulation. mature micrornas are produced in a multi-stage process. after transcription, they are processed by rna pol ii or pol iii to create capped and polyadenylated primary transcripts, which are further processed in the nucleus by the enzyme drosha/pasha or by dgcr to produce ~60-nucleotide pre-microrna stem-loop molecules. these are then exported to the cytoplasm by exportin and ran-gtp where they are further processed by dicer to ~ nt double-stranded rna duplexes, which form complexes with risc leading to unwinding of the duplexes to form single-stranded micrornas. micrornas bound to risc can down-regulate protein levels using at least two alternative pathways: 1) if the microrna has imperfect complementarity with a matching sequence in the 3'utr of its target mrna, the microrna-risc complex can combine with the complement mrna sequence and cause translational repression. 2) on the other hand, if the microrna and its mrna target have perfect or near perfect complementarity, the microrna-risc complex binding to its target mrnas can result in the cleavage and degradation of the mrna by argonaute. although many studies have identified signatures of microrna dysregulation, the identification of tissue specific targets of aberrantly regulated microrna is difficult. putative identification using seed sequence complementarity and free energy predictions of rna-rna duplexes are available in databases such as targetscan: http://www.targetscan.org. however, the false positive rate for such matches is unacceptably high, with different algorithms identifying different mrna targets for the same microrna. the tissue specificity of microrna regulation is known only in some specific cases and a general methodology for target identification, tissue specificity of action and specific biological role of microrna in the initiation and progression of most cancers remains an open problem. we describe a novel method to identify "direct mrna targets" of microrna in any cancer based on measuring an anti-correlation signal between differentially expressed microrna and mrna in patient matched tumor and normal samples. in this paper, the words "direct mrna targets" is used in a very strict and limited sense. a direct target is one which: a) has an exact seed sequence match in its 3'utr to the corresponding microrna, b) the seed sequence match is conserved across mouse, human, rat and dog genomes, c) the expression levels of both the microrna and the mrna can distinguish tumor from normal with high statistical significance and d) the mrna and microrna levels are strongly and significantly anti-correlated in tumor and/or normal. these requirements could be relaxed to find additional targets or eliminated altogether to find indirect regulations. the method proceeds as follows: a) identify significantly differentially expressed microrna and mrna between the two classes; b) for each microrna which is differentially expressed, identify all its putative target mrna by restricting to those differentially expressed mrna with a matching seed sequence in their 3'utr, with the further requirement that it be conserved in human, mouse, rat and dog genomes; c) compute the pearson correlation between microrna and mrna expression levels for samples in each class and d) retain only those microrna/mrna pairs whose expression levels are highly anti-correlated. these constraints remove spurious matches, reducing relatively speculative "putative" seed match based mrna targets in databases to a highly robust subset of direct functional targets. note that our method can be extended by removing constraint b) and looking for a correlation signature in c). this allows the identification of indirect regulation. for example, if a microrna up-regulated in cancer down-regulates a gene which is a transcriptional repressor of an oncogene, then the expression level of the microrna will be correlated with the level of the oncogene without a seed sequence match. since the direct gene target of the microrna should already be identified using our method, such an analysis would extend the regulation network beyond first order interactions. note that, although the method as described above does not identify regulation by translation inhibition, if protein levels were also measured, the method could easily be extended to identify such regulation. we demonstrate the use of our method on expression data from clear cell renal cell carcinoma and matched normal kidney samples. renal cell carcinoma represents ~3% of all malignancies in the us, with new cases and deaths each year http://www.nci.nih.gov/cancertopics/types/kidney. the most common histological class is ccrcc, accounting for ~75% of kidney cancers. ccrcc is known to be characterized by the loss of the vhl gene, which under normal oxygen pressure, binds to the α subunits of hypoxia-inducible factors, inducing their poly-ubiquitinylation and subsequent degradation in the proteasome. in hypoxic conditions, or if hif regulation is lost because of vhl inactivation, hif accumulates to high levels and promotes the transcription of genes such as vegf, pdgf-β, tgf-α, epo etc which trigger angiogenesis, cell growth, migration and proliferation. the spectrum of hif target genes expressed in individual tumors and the factors which influence them are the object of active ongoing research. ccrcc tumors have a wide range of natural histories and varied responses to vegf-targeted therapy. early stage, fuhrman grade tumors tend to have significantly better disease free survival after resection than higher stage and grade. although vhl mutation is associated with all grades of ccrcc, the other molecular factors associated with ccrcc initiation and progression are largely unknown. the molecular basis of the diversity in histologic grade, clinical behavior, and response to vegf-targeted is also unclear, and makes ccrcc a ripe target for studies investigating the molecular and genetic nature of these heterogeneities. in rcc, various studies have identified panels of microrna and mrna that are differentially expressed between normal renal tissue and tumor or between histological subtypes of tumor. the present study extends these previous studies by linking the microrna to some of their mrna targets, thus elucidating a hitherto unknown part of the biology of ccrcc disease. some of the identified microrna/mrna anti-correlations were validated on a new cohort of ccrcc/normal samples. sema6a was confirmed as a direct target of mir- by over-expressing mir- in a ccrcc cell line and showing strong down-regulation of the sema6a transcript.  the underlying hypothesis in our method is that the expression levels of microrna and their direct mrna targets should be strongly anti-correlated when averaged over matched samples in either tumor or normal tissue. the stepwise procedure is as follows: step 1: identify significantly up/down regulated micrornas in ccrcc samples vs normal samples. step 2: identify significantly up/down regulated mrnas in ccrcc samples vs normal samples. step 3: using targetscan, retain only the mrna in step which have a conserved seed sequence in their 3' utr for at least one of the microrna from step step 4: find anti-correlated pairs of up-regulated microrna and down-regulated mrna in ccrcc samples using a strict cutoff in pearson correlation coefficient in ccrcc samples. similarly, find anti-correlated pairs of down-regulated microrna and up-regulated mrnas in normal kidney samples using a strict cutoff. in step putative target mrna were identified using targetscan version http://www.targetscan.org, which identifies possible regulatory targets of mammalian micrornas as those with conserved sequences of matching seed regions for each microrna. the term "conserved" means that the sequence is conserved in human, mouse, rat and dog. in step we ran permutations in brb-arraytools' http://linus.nci.nih.gov/brb-arraytools.html and multivariate/univariate analysis at p <, fdr < to assess significance of discovered pairs and to find the appropriate cutoff p for significance of the measured pearson correlation. in the primary dataset, because of the high accuracy of qrt-pcr, we were able to set a strict cutoff p = -; thus in step only microrna/mrna pairs with p < p were considered to be significant. step 1: microrna significantly differentially expressed in ccrcc versus normal kidney tissue microrna were identified as differentially expressed in ccrcc versus normal kidney, down-regulated and up-regulated. the microrna down-regulated in ccrcc were mir- mir-10b, mir-125b, mir-26a+, mir-133b, mir-135a, mir-135b, mir- mir- mir- mir- mir-199a, mir-200a, mir-200b, mir-200c, mir- mir- mir- mir-30a-3p, mir-30a-5p, mir- mir- mir- mir- mir- mir- and the microrna up-regulated in ccrcc were mir-142-3p, mir- mir- mir- mir- mir- mir-34a, mir-34b, mir-, in agreement with recent studies. the chromosomal location of these microrna shows that several are proximal to fragile regions and regions commonly deleted in ccrcc. figure 1a shows a heat-map of the microrna expression levels for samples in the primary dataset showing that they can robustly separate ccrcc from normal kidney. cancer-related regions listed are either deleted regions, amplified regions or breakpoint regions. step 2: mrna differentially expressed in ccrcc versus normal kidney tissue using sam http://www-stat.stanford.edu/~tibs/sam/ at fdr = and permutation tests, mrna probes were found up-regulated and were found down-regulated in ccrcc. the genes represented by these probes included many genes known to be dysregulated in ccrcc. we saw significant upregulation of oncogenes vegfa and egfr, integrin-mediated cell adhesion pathway genes, cell adhesion related genes and cytokines. similarly we noted down-regulation of tumor suppressor genes and of members of the metallothionein family. steps and 4: identifying direct mrna targets of dysregulated microrna pearson correlation analysis with p = - was applied to each of the differentially expressed microrna and its putative mrna targets. this procedure identified mrna targets for the up-regulated microrna. this list included several important tumor suppressor genes, such as pten, erbb and sfrp known to be mutated or down-regulated in many tumors, including ccrcc. the down-regulated microrna had direct up-regulated mrna targets, including oncogenes vegfa, lox, loxl and fas, well known to be involved in kidney cancer. the nine most significantly down-regulated and twelve most significantly up-regulated mrna are listed in table and their heatmap is shown in figure 1b. in figure we plot microrna and mrna levels for mir-200c and its target vegfa. note that the levels of mir-200c and its target vegfa are not only anti-correlated overall, but are also anti-correlated separately in both ccrcc and normal tissue. additional files and contain the full list of microrna/mrna regulations identified by our analysis. the full list of identified microrna/mrna pairs is given in additional files and + marks microrna/mrna anti-correlations that were tested/validated on an independent set of matched ccrcc/normal kidney samples. validation on a cohort of ccrcc/normal kidney samples a. validation of differentially expressed microrna and mrna seventeen validation samples were collected from a new cohort of patients and analyzed by microrna and mrna profiling on agilent arrays. figure 1c shows a heat map of twenty-seven microrna of the thirty five differentially expressed microrna identified previously. a weighted voting classifier on binarized microrna expression data in the validation set had 100% accuracy in leave-one-out cross validation in distinguishing ccrcc from normal kidney. figure 1d shows a heatmap of mrna expression levels for the genes in figure 1b. these were also found to be 100% accurate at discriminating ccrcc from normal kidney using weighted voting and loo cross validation analysis. mrna levels of three genes which were down-regulated in ccrcc and one gene up-regulated in ccrcc were also measured by quantitative rt-pcr in twelve of the test samples. the results, shown in figure 3a, demonstrate that mrna levels of erbb sfrp slc12a and vegfa were quantitatively and significantly down/up-regulated as expected. b. validation of anti-correlation signature between some identified microrna/mrna pairs in figure 3b-e, we plot the qrt-pcr expression levels of erbb sfrp slc12a and vegfa versus agilent chip measured levels of their regulatory microrna for the twelve samples of figure 3a. the overall strong anti-correlation signature between microrna and mrna levels is clearly visible in these plots. figure summarizes our validation analysis of a number of anti-correlation measurements between several identified microrna/mrna pairs in the agilent chip data. the measured correlations between predicted microrna/mrna pairs are shown in the figure. because of the higher level of noise in the agilent chip data compared to qrt-pcr, we cannot apply the strict criterion used in the primary dataset. instead, the significance of the correlation was computed using permutation tests as follows: a large number of datasets were obtained by permuting the sample labels in the microrna or mrna measurements. for each microrna/mrna pair, these permuted datasets were used to compute the null distribution for p and the significance of the measured value of p was estimated in this null distribution. as the p-values in figure indicate, we validate a strong anti-correlation signature between mrna levels of, vegf, sema6a,, sfrp erbb slc12a and tcf and their identified regulators: mir- mir-200c, mir- mir-142-3p, mir- mir-34a, mir- and mir- respectively. c. in vitro validation of sema6a as a target of mir- in an rcc cell line finally, to establish that this method can accurately predict functional and direct microrna/mrna regulation, we performed an in vitro analysis of one microrna, and its identified direct target sema6a. the rcc cell line crl- was transfected with either pre-mir- or a control pre-mir, and levels of sema6a were measured on the case/control cell lines by semi-quantitative rt-pcr. the results showed that introduction of pre-mir- produced a significant reduction in the level of sema6a mrna, validating sema6a as a functional and direct target of mir- discussion we have developed and validated a simple method to identify direct functional mrna targets of microrna in ccrcc using patient matched tumor/normal samples. the method looks for the subset of anti-correlated microrna/mrna pairs from a larger set of microrna and mrna able to distinguish tumor from normal tissue, with the additional requirement of a highly conserved seed sequences for the corresponding microrna in the 3'utr of the corresponding gene. our method can not only identify microrna/mrna pairs which discriminate normal from cancer tissue, but also dysregulated regulation mechanisms between them that may initiate and/or drive the disease process. we used an rt-pcr panel for our initial discovery experiments and microrna and microarray chips both for validation of our previously discovered microrna/mrna pairings and to expand the available pool of micrornas to all currently known micrornas for a comprehensive examination of microrna/mrna pairs. in comparing the two data modalities, we noticed that in the rt-pcr data, the sample to sample variation in the normal samples was comparable to the change in level between normal and tumor which caused a statistically significant up/down regulation of the target mrna. to find reliable matches, it was therefore crucial to minimize stochastic error. we observed that the rt-pcr data had lower stochastic variance than the microarray data. this noise effect is also reflected in the measured p values in the validation analysis on the microarray data, which are smaller in magnitude than the strict cutoff |p0| > used in the discovery phase on the rt-pcr data. since the present study is mainly to demonstrate "proof of concept", we limited the rt-pcr data for discovery and the microarray data for validation. in a more extensive study, with sufficiently large sample sizes and more accurate data from high throughput technologies it may be feasible and cost effective to use a high throughput platform for discovery and rt-pcr for validation. such an analysis might also identify a much bigger set of microrna/mrna relationships. although our method is very robust, as described, it cannot find all mrna targets. its principal merit is the ability to reduce the large number of relatively speculative matches from seed sequences alone to a smaller set of functional, tissue specific targets. while this smaller set is perhaps incomplete, it is more reliable. additionally, the use of perfect seed sequence matches can easily be relaxed. a more serious issue is that the method as described cannot find targets affected by translation inhibition. for such targets, the changes in microrna levels would leave mrna levels unaffected, but reduce protein levels. such targets might be identified if protein levels in tumor/normal samples were also measured along with the microrna and mrna levels. the method could then be adapted to identify targets affected by translation inhibition by using an appropriate joint correlation/anti-correlation analysis of microrna/mrna/protein levels. our method, as described here, would also miss regulation that proceeds via intermediate states. in these cases, there would be no seed sequence match between a microrna and the mrna whose level was affected by it. we could attempt to identify such secondary effects by eliminating the seed sequence match requirement but keeping the requirement of a high correlation/anti-correlation signal between microrna/mrna levels. whether such a procedure would work would depend on the degree to which the intermediate state gene/protein is affected and measurable. we briefly discuss some of these issues below. however, it should be noted that these studies would require significantly larger sample sizes and are beyond the scope of the present paper. the next round of tcga http://cancergenome.nih.gov/ may make these types of analyses more feasible. although here, we focused on microrna/mrna pairs with exact seed sequence match and an anti-correlation signal in both tumor and normal samples, as noted above, with additional samples it should be possible to use this method and simple extensions to identify more subtle types of dysregulation. for example, if a mutation in a microrna in the tumor samples causes loss of its function, then although the microrna/mrna levels would be anti-correlated in normal tissue, they would not be anti-correlated in tumor samples. similarly a de-novo gain of microrna function would be signaled by an anti-correlated signal in tumor samples which was absent in the normal samples. in our data, we did observe several correlations and anti-correlations between microrna/mrna pairs in tumor or normal samples without a corresponding seed sequence match. as noted above, these most likely represent regulations which proceed via intermediate states and hence cannot be identified by seed sequence matches. for example, if increased levels of microrna x down-regulates expression of a protein which is a transcriptional repressor of gene y, there will be a strong correlation between x and y levels but no seed sequence match. conversely, if microrna x regulates a protein which is a transcription factor for gene y then the levels of x and y will be anti-correlated without a seed sequence match. the measurement of such correlations would extend the network of microrna control beyond first level regulators but would require significant increases in the number of samples for statistical significance. since the method finds functional relationships, it should be useful for identifying pharmaceutically relevant mechanisms which suggest drug targets for therapy. we describe below some of the regulations we identified which have pharmaceutical relevance. microrna which shut down multiple tumor suppressor genes in ccrcc we found several dysregulated microrna which targeted multiple tumor suppressor genes. for example, the oncogenic mir- was significantly up-regulated in ccrcc and anti-correlated with the tumor suppressor gene pten, suggesting that its gain of function shuts down pten in ccrcc. mutated or down-regulated in many advanced cancers, pten loss activates the pi3k-akt signaling pathway and its downstream target mtor, with important implications in rcc development and therapeutic selection. another identified target of mir- was ptpn, which is a fas-associated protein tyrosine phosphatase and putative tumor suppressor gene that can inhibit pi3k/akt signaling, suppress the influence of insulin-like growth factor-i on cell survival and induce apoptosis. kcnj member of the potassium channel subfamily of membrane proteins, was also identified as a target of mir- such membrane proteins have been suggested in anti-cancer therapies because of their important role in cell growth and are known to be down-regulated in rcc. mir-34a is known to be over-expressed in various tumors and associated with cell proliferation. in our data, it was up-regulated in ccrcc and predicted to target sfrp a known regulator of the wnt signaling pathway and a tumor suppressor gene whose loss has been observed in a majority of rcc patients. another tumor suppressor regulated by mir-34a was camta, a reduction in whose levels correlates with poor outcome in neuroblastoma. finally, we find mir- associated with a human chromosome fragile site on chr-xq, is up-regulated in ccrcc and predicted to target erbb a member of the egfr family, a potential tumor suppressor known to be strongly down-regulated in ccrcc. hypoxia induced micrornas a hypoxic tumor microenvironment can directly activate the expression of several microrna. for example, mir- mir- and mir- reduce pro-apoptotic signaling in response to a hypoxic environment and are consistently over-expressed in a variety of human tumors. ccrcc is a unique setting in which to study these micrornas, given that vhl loss constitutively stabilizes one or more hif factors, thereby creating a pseudo-hypoxic scenario in ccrcc tumor cells. we found that hypoxia related microrna had the most significant fold changes, with mir- mir- and mir- being amongst the top., suggesting a major role for them in renal carcinogenesis. we found genes down-regulated in tumor, inversely correlated with mir- and enriched in the cell adhesion pathway at p =. loss of cam degrades the intra and extracellular matrix, leading to abnormal cell growth patterns. among these genes, two were identified as direct targets of mir-21: slc12a and tcf both of which have been reported down-regulated in ccrcc. identification of microrna family interactions we found evidence that families of micrornas may be coordinately participating in microrna/mrna interactions. one example in ccrcc is the mir-204/ family, which was significantly down-regulated in ccrcc samples. we identified thirty five significant mrna targets for this family. among these, eight were on chromosome 3q, a common amplicon region in many epithelial tumors. these genes include c3orf ccdc dtx3l, pld trim two oncogenes ect and rap2b, and a hypoxia associated protein serp. the gain of this chromosome arm was previously associated with papillary rccs, and our observation in clear cell rcc implies a possible regulatory relation between mir-204/ and the genes in this region as an alternate mechanism of up-regulation of this group of genes. another example is the mir- family which includes two microrna clusters, one on chromosome 1p and another on chromosome 12p. five mir- family members contain very similar seed sequences - aauacu for mir-200b/200c/ and aacacu for mir-200a*/. recently, several other groups have reported a role for the mir- family in the epithelial-mesenchymal-transition and in cancer cell migration, the latter by directly targeting the transcription factors zeb and zeb which regulate e-cadherin, a mediator of cell-cell adhesion. another study identified a regulatory loop between these micrornas and zeb transcription factors as well as the emt inducer tgfβ. in epithelial cells, mir- family micrornas and e-cadherin maintain higher level expression by repressing zeb zeb and tgfβ; on the other hand, in mesenchymal cells and tumors, the up-regulation of zeb factors is triggered by tgfβ and suppresses the transcription of mir-141/200c by binding to their putative common promoter region. in our primary dataset, zeb and zeb were both up-regulated in six out of our eight ccrcc samples and, their expression levels were highly anti-correlated with the mir- family in both tumor and normal samples. as confirmation of these results, down-regulation of mir- and mir-200c and their function on zeb in ccrcc has recently been reported. we also noted that in our data, the anti-correlation between vegfa and the mir- family was strongest in normal kidney tissue, suggesting that loss of this regulation may be an important factor providing a permissive environment for hif transcriptional signaling. our hypothesis from these various observations is that in normal kidney, the expression level of hif2α and its downstream targets are regulated by mir- 200a*, 200b and 200c and the loss of this microrna regulation, in concert with vhl loss, is responsible for activation of the hif pathway. one intriguing association which we have identified is highly significant for therapy in ccrcc. this is because the soluble extracellular domain of sema6a has been engineered to effectively inhibit vegf-mediated tumor formation. hence, our results imply that mir- may have a role in gene therapy. a model which summarizes our observations and integrates it with mechanisms for ccrcc dysregulation from the literature is shown in figure this model integrates our measured microrna regulatory mechanisms with known transcriptional activity resulting from vhl loss and activation of the hypoxia response pathway. the pathways shown in figure are likely highly interconnecting, and this model, and the individual functional interactions it suggests, need to be validated by direct experimental targeting of microrna levels and measurements of each of the predicted target genes in a larger cohort in its next phase).  the main innovation in this paper is the use of an anti-correlation analysis of microrna/mrna levels in two cell types from the same patient to identify functional mrna targets of altered microrna. the method can obviously extend to any tissue type and might be useful in other contexts: for instance, by using microdissection to harvest cells from different cellular compartments from the same breast cancer patient, it might reveal the microrna/mrna regulation program which causes progression of normal epithelium to hyperplasia to dcis to invasive disease. the method might also be useful in non-cancer studies, such as in identifying the role of microrna and their target genes in the transition from stem cells to differentiated cells or in embryogenesis. in summary, in this paper, we have demonstrated a simple method to identify tissue specific mrna targets of microrna, which is expandable to many study types. 
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. these are further pruned by requiring a strong anti-correlation signature in tumor and/or normal samples. putative mrna targets of microrna identified from seed sequence matches are available in many databases. we apply and validate the method using clear cell renal cell carcinoma 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. 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. the constraints we impose remove spurious matches and identify a subset of robust, tissue specific, functional mrna targets of dysregulated microrna. however, such matches have a high false positive rate and cannot identify tissue specificity of regulation. microrna regulate mrna levels in a tissue specific way, either by inducing degradation of the transcript or by inhibiting translation or transcription. dysregulated microrna and mrna are first identified by comparing their expression levels in tumor vs normal samples. we also found strong anti-correlation between vegfa and the mir- family of microrna: mir-200a*, 200b, 200c and mir- several identified microrna/mrna pairs were validated on an independent set of matched ccrcc/normal samples.
 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.  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 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- mir-200c and mir- causes gain of function of oncogenes, vegfa and sema6a respectively and increased levels of mir-142-3p, mir- mir-34a, mir- mir- cause loss of function of tumor suppressors lrrc ptpn sfrp erbb and respectively. we also found strong anti-correlation between vegfa and the mir- family of microrna: mir-200a*, 200b, 200c and mir- several identified microrna/mrna pairs were validated on an independent set of matched ccrcc/normal samples. the regulation of sema6a by mir- was verified by a transfection assay.  we describe a simple and reliable method to identify direct gene targets of microrna in any cancer. the constraints we impose remove spurious matches and identify a subset of robust, tissue specific, functional mrna targets of dysregulated microrna. 