This pipeline infers putative direct gene targets of miRs based on miRseq and mRNAseq expression profiles across multiple samples. This pipeline has the following steps:
The CLR approach is applied to infer putative miR:gene regulatory connections.
Filtering miR:gene pairs based on Pearson correlation (<= -0.3 ).
Filtering miR:gene pairs based on predicted interactions in three sequence prediction databases (Miranda, Pictar, Targetscan)
The CLR algorithm was applied on 815 miRs and 18514 mRNAs across 342 samples. After 2 filtering steps, the number of 125 miR:genes pairs were detected.
|hsa-miR-96-5p||15||CPEB4, EPAS1, GPR124, ITPR2, LAT2, LPPR4, MITF, MRAS, PDGFRA, PTGER3, RFTN1, SMAD7, TNS3, VCL,CACNA1C|
|hsa-miR-141-3p||13||CXCL12, DIXDC1, DLC1, FOXP1, MYH10, NRP1, PRKACB, TMEM110, TNS3, WIPF1, ZCCHC24, ZEB2,APBB2|
|hsa-miR-200c-3p||12||CYP1B1, ETS1, LHFP, PALM2-AKAP2, PRDM1, RECK, TLN1, WIPF1, ZCCHC24, ZEB2, ZFPM2,CDH11|
|hsa-miR-200b-3p||11||ETS1, PALM2-AKAP2, RECK, RHOA, SEC23A, SYDE1, TLN1, ZCCHC24, ZEB1, ZEB2,CFL2|
|hsa-miR-205-5p||6||LHFPL2, MGRN1, SORBS1, TRAK2, VASN,AXIN2|
|hsa-miR-16-5p||6||HSPG2, ISLR, LAMB2, N4BP1, PAPPA,COL12A1|
|hsa-miR-32-5p||6||ITGA5, ITGAV, LUZP1, PMEPA1, TRIO,COL1A2|
|hsa-miR-25-3p||5||CXXC5, EDEM1, LHFPL2, MCL1,CPEB4|
|hsa-miR-182-5p||5||CBFA2T3, CD47, MBNL2, RFTN1,BHLHE22|
|hsa-miR-9-5p||4||DUSP6, ERG, RHOJ,CPEB2|
This section should list the files that were used as input.
miRseq (at precursor level) of RPM value (reads per million reads aligned to miRBase precursor) with log2 transformed = /xchip/cga/gdac-prod/tcga-gdac/jobResults/miRseq_mature_preprocess/LUSC-TP/22315347/LUSC-TP.miRseq_mature_RPM_log2.txt
mRNAseq of RSEM/RPKM value with log2 transformed = /xchip/cga/gdac-prod/tcga-gdac/jobResults/mRNAseq_preprocessor/LUSC-TP/22315346/LUSC-TP.uncv2.mRNAseq_RSEM_normalized_log2.txt
miR:gene predicted interactions file = /xchip/cga_home/hailei/FH/CLR/human_interactions.predicted.v2.mirbase21.txt
Miranda = microrna.org Aug 2010 release, Microcosm version 5
Pictar = version 1
Targetscan = release 5.2
The CLR (Context Likelihood of Relatedness) algorithm builds upon the relevance network strategies, by applying a background correction step. After computing the mutual information between regulators and their potential target genes, CLR calculates the statistical likelihood of each mutual information value within its network context. The algorithm compares the mutual information between a miR/gene pair to the background distribution of mutual information scores for all possible miR/gene pairs that include either the miR or its target. After this background correction, the most probable interactions are those whose mutual information scores stand significantly above the background distribution of mutual information scores .
Pairwise Pearson correlations coefficients between all miR:gene pairs are first computed. All genes that have correlation values less than the user-defined threshold (-0.3) with a particular miR and have been predicted as targets of that miR in three sequence based prediction databases: Miranda Pictar, TargetScan  are identified as putative direct targets of that miR. We infer a direct target miR:gene network which comprises all such putative direct associations.
threshold = -0.3
In addition to the links below, the full results of the analysis summarized in this report can also be downloaded programmatically using firehose_get, or interactively from either the Broad GDAC website or TCGA Data Coordination Center Portal.