ID	Duplex structure	Position	Score	MFE
1	miRNA 3' acguccgucuuccccaCCACUCc 5' || ||| Target 5' tcttttgactttttgcGGGGAGg 3'	135 - 157	88.00	-10.35
2	miRNA 3' acGUCCGUCUUCCCCACCacucc 5' |||| | | ||| Target 5' gtCAGGGCTATGCTATGGgcctt 3'	31 - 53	69.00	-5.50
3	miRNA 3' acguccgucuuccccaccACUcc 5' ||| Target 5' ctttcttcccaggataacTGAag 3'	112 - 134	60.00	-7.67

Mutant ID	Mutant Position	Mutant Source
COSN31486338	11	COSMIC
COSN30454765	36	COSMIC
COSN8646064	38	COSMIC
COSN15900125	101	COSMIC
COSN30513898	102	COSMIC
COSN20609487	135	COSMIC
COSN19749875	151	COSMIC
COSN31485761	189	COSMIC
COSN30113416	198	COSMIC
COSN31556743	210	COSMIC
COSN30126216	216	COSMIC
COSN30111185	230	COSMIC
COSN30124990	256	COSMIC
COSN30132489	268	COSMIC
COSN31530164	298	COSMIC

Mutant ID	Mutant Position	Mutant Source
rs1455611035	2	dbSNP
rs775967539	4	dbSNP
rs761412401	7	dbSNP
rs1415606586	10	dbSNP
rs766922629	11	dbSNP
rs1348667041	12	dbSNP
rs1465142681	13	dbSNP
rs760176269	15	dbSNP
rs1226521050	17	dbSNP
rs1048617819	25	dbSNP
rs1325833858	34	dbSNP
rs767967823	35	dbSNP
rs1217747245	37	dbSNP
rs1417511007	41	dbSNP
rs1161846335	42	dbSNP
rs1254000453	46	dbSNP
rs753144068	48	dbSNP
rs12562680	58	dbSNP
rs1474012981	62	dbSNP
rs1387311591	64	dbSNP
rs1181776250	71	dbSNP
rs888591388	76	dbSNP
rs1235857379	77	dbSNP
rs1200061206	79	dbSNP
rs60372173	92	dbSNP
rs60515741	94	dbSNP
rs879359847	98	dbSNP
rs1288776732	99	dbSNP
rs1445517104	109	dbSNP
rs1358250214	110	dbSNP
rs548995950	114	dbSNP
rs1037250594	116	dbSNP
rs1445985670	116	dbSNP
rs12073773	120	dbSNP
rs1333454941	123	dbSNP
rs1378571283	130	dbSNP
rs1349506698	134	dbSNP
rs1321906418	136	dbSNP
rs1407408701	138	dbSNP
rs1013958839	143	dbSNP
rs41299262	150	dbSNP
rs1042881	151	dbSNP
rs41299264	152	dbSNP
rs999881047	154	dbSNP
rs1290370004	155	dbSNP
rs1239255573	158	dbSNP
rs1339180096	160	dbSNP
rs1032822646	161	dbSNP
rs1239041211	165	dbSNP
rs1306388620	181	dbSNP
rs1352251782	184	dbSNP
rs955839684	191	dbSNP
rs1224201262	192	dbSNP
rs866607150	202	dbSNP
rs1440698704	208	dbSNP
rs1205298762	209	dbSNP
rs375069063	212	dbSNP
rs1461907993	216	dbSNP
rs1018902637	229	dbSNP
rs966347992	246	dbSNP
rs1423885198	248	dbSNP
rs10919261	251	dbSNP
rs12128362	263	dbSNP
rs1159277260	277	dbSNP
rs1247013697	280	dbSNP
rs1215613449	281	dbSNP
rs1425074755	282	dbSNP
rs1413306972	292	dbSNP
rs976808057	297	dbSNP
rs1217271077	298	dbSNP
rs1311634926	304	dbSNP

- Cell, 2010

RNA transcripts are subject to posttranscriptional gene regulation involving hundreds of RNA-binding proteins (RBPs) and microRNA-containing ribonucleoprotein complexes (miRNPs) expressed in a cell-type dependent fashion. We developed a cell-based crosslinking approach to determine at high resolution and transcriptome-wide the binding sites of cellular RBPs and miRNPs. The crosslinked sites are revealed by thymidine to cytidine transitions in the cDNAs prepared from immunopurified RNPs of 4-thiouridine-treated cells. We determined the binding sites and regulatory consequences for several intensely studied RBPs and miRNPs, including PUM2, QKI, IGF2BP1-3, AGO/EIF2C1-4 and TNRC6A-C. Our study revealed that these factors bind thousands of sites containing defined sequence motifs and have distinct preferences for exonic versus intronic or coding versus untranslated transcript regions. The precise mapping of binding sites across the transcriptome will be critical to the interpretation of the rapidly emerging data on genetic variation between individuals and how these variations contribute to complex genetic diseases.

- Nature methods, 2011

Cross-linking and immunoprecipitation (CLIP) is increasingly used to map transcriptome-wide binding sites of RNA-binding proteins. We developed a method for CLIP data analysis, and applied it to compare CLIP with photoactivatable ribonucleoside-enhanced CLIP (PAR-CLIP) and to uncover how differences in cross-linking and ribonuclease digestion affect the identified sites. We found only small differences in accuracies of these methods in identifying binding sites of HuR, which binds low-complexity sequences, and Argonaute 2, which has a complex binding specificity. We found that cross-link-induced mutations led to single-nucleotide resolution for both PAR-CLIP and CLIP. Our results confirm the expectation from original CLIP publications that RNA-binding proteins do not protect their binding sites sufficiently under the denaturing conditions used during the CLIP procedure, and we show that extensive digestion with sequence-specific RNases strongly biases the recovered binding sites. This bias can be substantially reduced by milder nuclease digestion conditions.

- Genes & development, 2011

MicroRNAs are important regulators in many cellular processes, including stem cell self-renewal. Recent studies demonstrated their function as pluripotency factors with the capacity for somatic cell reprogramming. However, their role in human embryonic stem (ES) cells (hESCs) remains poorly understood, partially due to the lack of genome-wide strategies to identify their targets. Here, we performed comprehensive microRNA profiling in hESCs and in purified neural and mesenchymal derivatives. Using a combination of AGO cross-linking and microRNA perturbation experiments, together with computational prediction, we identified the targets of the miR-302/367 cluster, the most abundant microRNAs in hESCs. Functional studies identified novel roles of miR-302/367 in maintaining pluripotency and regulating hESC differentiation. We show that in addition to its role in TGF-beta signaling, miR-302/367 promotes bone morphogenetic protein (BMP) signaling by targeting BMP inhibitors TOB2, DAZAP2, and SLAIN1. This study broadens our understanding of microRNA function in hESCs and is a valuable resource for future studies in this area.

- PLoS pathogens, 2012

KSHV is the etiological agent of Kaposi's sarcoma (KS), primary effusion lymphoma (PEL), and a subset of multicentricCastleman's disease (MCD). The fact that KSHV-encoded miRNAs are readily detectable in all KSHV-associated tumors suggests a potential role in viral pathogenesis and tumorigenesis. MiRNA-mediated regulation of gene expression is a complex network with each miRNA having many potential targets, and to date only few KSHV miRNA targets have been experimentally determined. A detailed understanding of KSHV miRNA functions requires high-through putribonomics to globally analyze putative miRNA targets in a cell type-specific manner. We performed Ago HITS-CLIP to identify viral and cellular miRNAs and their cognate targets in two latently KSHV-infected PEL cell lines. Ago HITS-CLIP recovered 1170 and 950 cellular KSHV miRNA targets from BCBL-1 and BC-3, respectively. Importantly, enriched clusters contained KSHV miRNA seed matches in the 3'UTRs of numerous well characterized targets, among them THBS1, BACH1, and C/EBPbeta. KSHV miRNA targets were strongly enriched for genes involved in multiple pathways central for KSHV biology, such as apoptosis, cell cycle regulation, lymphocyte proliferation, and immune evasion, thus further supporting a role in KSHV pathogenesis and potentially tumorigenesis. A limited number of viral transcripts were also enriched by HITS-CLIP including vIL-6 expressed only in a subset of PEL cells during latency. Interestingly, Ago HITS-CLIP revealed extremely high levels of Ago-associated KSHV miRNAs especially in BC-3 cells where more than 70% of all miRNAs are of viral origin. This suggests that in addition to seed match-specific targeting of cellular genes, KSHV miRNAs may also function by hijacking RISCs, thereby contributing to a global de-repression of cellular gene expression due to the loss of regulation by human miRNAs. In summary, we provide an extensive list of cellular and viral miRNA targets representing an important resource to decipher KSHV miRNA function.

- Cell, 2013

The induction of pluripotency or trans-differentiation of one cell type to another can be accomplished with cell-lineage-specific transcription factors. Here, we report that repression of a single RNA binding polypyrimidine-tract-binding (PTB) protein, which occurs during normal brain development via the action of miR-124, is sufficient to induce trans-differentiation of fibroblasts into functional neurons. Besides its traditional role in regulated splicing, we show that PTB has a previously undocumented function in the regulation of microRNA functions, suppressing or enhancing microRNA targeting by competitive binding on target mRNA or altering local RNA secondary structure. A key event during neuronal induction is the relief of PTB-mediated blockage of microRNA action on multiple components of the REST complex, thereby derepressing a large array of neuronal genes, including miR-124 and multiple neuronal-specific transcription factors, in nonneuronal cells. This converts a negative feedback loop to a positive one to elicit cellular reprogramming to the neuronal lineage.

- Nature, 2013

Circular RNAs (circRNAs) in animals are an enigmatic class of RNA with unknown function. To explore circRNAs systematically, we sequenced and computationally analysed human, mouse and nematode RNA. We detected thousands of well-expressed, stable circRNAs, often showing tissue/developmental-stage-specific expression. Sequence analysis indicated important regulatory functions for circRNAs. We found that a human circRNA, antisense to the cerebellar degeneration-related protein 1 transcript (CDR1as), is densely bound by microRNA (miRNA) effector complexes and harbours 63 conserved binding sites for the ancient miRNA miR-7. Further analyses indicated that CDR1as functions to bind miR-7 in neuronal tissues. Human CDR1as expression in zebrafish impaired midbrain development, similar to knocking down miR-7, suggesting that CDR1as is a miRNA antagonist with a miRNA-binding capacity ten times higher than any other known transcript. Together, our data provide evidence that circRNAs form a large class of post-transcriptional regulators. Numerous circRNAs form by head-to-tail splicing of exons, suggesting previously unrecognized regulatory potential of coding sequences.

- mBio, 2013

UNLABELLED: The question of how HIV-1 interfaces with cellular microRNA (miRNA) biogenesis and effector mechanisms has been highly controversial. Here, we first used deep sequencing of small RNAs present in two different infected cell lines (TZM-bl and C8166) and two types of primary human cells (CD4(+) peripheral blood mononuclear cells [PBMCs] and macrophages) to unequivocally demonstrate that HIV-1 does not encode any viral miRNAs. Perhaps surprisingly, we also observed that infection of T cells by HIV-1 has only a modest effect on the expression of cellular miRNAs at early times after infection. Comprehensive analysis of miRNA binding to the HIV-1 genome using the photoactivatable ribonucleoside-induced cross-linking and immunoprecipitation (PAR-CLIP) technique revealed several binding sites for cellular miRNAs, a subset of which were shown to be capable of mediating miRNA-mediated repression of gene expression. However, the main finding from this analysis is that HIV-1 transcripts are largely refractory to miRNA binding, most probably due to extensive viral RNA secondary structure. Together, these data demonstrate that HIV-1 neither encodes viral miRNAs nor strongly influences cellular miRNA expression, at least early after infection, and imply that HIV-1 transcripts have evolved to avoid inhibition by preexisting cellular miRNAs by adopting extensive RNA secondary structures that occlude most potential miRNA binding sites. IMPORTANCE: MicroRNAs (miRNAs) are a ubiquitous class of small regulatory RNAs that serve as posttranscriptional regulators of gene expression. Previous work has suggested that HIV-1 might subvert the function of the cellular miRNA machinery by expressing viral miRNAs or by dramatically altering the level of cellular miRNA expression. Using very sensitive approaches, we now demonstrate that neither of these ideas is in fact correct. Moreover, HIV-1 transcripts appear to largely avoid regulation by cellular miRNAs by adopting an extensive RNA secondary structure that occludes the ability of cellular miRNAs to interact with viral mRNAs. Together, these data suggest that HIV-1, rather than seeking to control miRNA function in infected cells, has instead evolved a mechanism to become largely invisible to cellular miRNA effector mechanisms.

- Genome biology, 2013

BACKGROUND: In recent years, a variety of small RNAs derived from other RNAs with well-known functions such as tRNAs and snoRNAs, have been identified. The functional relevance of these RNAs is largely unknown. To gain insight into the complexity of snoRNA processing and the functional relevance of snoRNA-derived small RNAs, we sequence long and short RNAs, small RNAs that co-precipitate with the Argonaute 2 protein and RNA fragments obtained in photoreactive nucleotide-enhanced crosslinking and immunoprecipitation (PAR-CLIP) of core snoRNA-associated proteins. RESULTS: Analysis of these data sets reveals that many loci in the human genome reproducibly give rise to C/D box-like snoRNAs, whose expression and evolutionary conservation are typically less pronounced relative to the snoRNAs that are currently cataloged. We further find that virtually all C/D box snoRNAs are specifically processed inside the regions of terminal complementarity, retaining in the mature form only 4-5 nucleotides upstream of the C box and 2-5 nucleotides downstream of the D box. Sequencing of the total and Argonaute 2-associated populations of small RNAs reveals that despite their cellular abundance, C/D box-derived small RNAs are not efficiently incorporated into the Ago2 protein. CONCLUSIONS: We conclude that the human genome encodes a large number of snoRNAs that are processed along the canonical pathway and expressed at relatively low levels. Generation of snoRNA-derived processing products with alternative, particularly miRNA-like, functions appears to be uncommon.

- Genome biology, 2014

BACKGROUND: Various microRNAs (miRNAs) are up- or downregulated in tumors. However, the repression of cognate miRNA targets responsible for the phenotypic effects of this dysregulation in patients remains largely unexplored. To define miRNA targets and associated pathways, together with their relationship to outcome in breast cancer, we integrated patient-paired miRNA-mRNA expression data with a set of validated miRNA targets and pathway inference. RESULTS: To generate a biochemically-validated set of miRNA-binding sites, we performed argonaute-2 photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (AGO2-PAR-CLIP) in MCF7 cells. We then defined putative miRNA-target interactions using a computational model, which ranked and selected additional TargetScan-predicted interactions based on features of our AGO2-PAR-CLIP binding-site data. We subselected modeled interactions according to the abundance of their constituent miRNA and mRNA transcripts in tumors, and we took advantage of the variability of miRNA expression within molecular subtypes to detect miRNA repression. Interestingly, our data suggest that miRNA families control subtype-specific pathways; for example, miR-17, miR-19a, miR-25, and miR-200b show high miRNA regulatory activity in the triple-negative, basal-like subtype, whereas miR-22 and miR-24 do so in the HER2 subtype. An independent dataset validated our findings for miR-17 and miR-25, and showed a correlation between the expression levels of miR-182 targets and overall patient survival. Pathway analysis associated miR-17, miR-19a, and miR-200b with leukocyte transendothelial migration. CONCLUSIONS: We combined PAR-CLIP data with patient expression data to predict regulatory miRNAs, revealing potential therapeutic targets and prognostic markers in breast cancer.

- Breast cancer research and treatment, 2014

miRNAs regulate the expression of genes in both normal physiology and disease. While miRNAs have been demonstrated to play a pivotal role in aspects of cancer biology, these reports have generally focused on the regulation of single genes. Such single-gene approaches have significant limitations, relying on miRNA expression levels and heuristic predictions of mRNA-binding sites. This results in only circumstantial evidence of miRNA-target interaction and typically leads to large numbers of false positive predictions. Here, we used a genome-wide approach (high-throughput sequencing of RNA isolated by crosslinking immunoprecipitation, HITS-CLIP) to define direct miRNA-mRNA interactions in three breast cancer subtypes (estrogen receptor positive, Her2 amplified, and triple negative). Focusing on steroid receptor signaling, we identified two novel regulators of the ER pathway (miR-9-5p and miR-193a/b-3p), which together target multiple genes involved in ER signaling. Moreover, this approach enabled the definition of miR-9-5p as a global regulator of steroid receptor signaling in breast cancer. We show that miRNA targets and networks defined by HITS-CLIP under physiologic conditions are predictive of patient outcomes and provide global insight into miRNA regulation in breast cancer.

- Neoplasia (New York, N.Y.), 2016

MicroRNA (miRNA) deregulation in prostate cancer (PCa) contributes to PCa initiation and metastatic progression. To comprehensively define the cancer-associated changes in miRNA targeting and function in commonly studied models of PCa, we performed photoactivatable ribonucleoside-enhanced cross-linking immunoprecipitation of the Argonaute protein in a panel of PCa cell lines modeling different stages of PCa progression. Using this comprehensive catalogue of miRNA targets, we analyzed miRNA targeting on known drivers of PCa and examined tissue-specific and stage-specific pathway targeting by miRNAs. We found that androgen receptor is the most frequently targeted PCa oncogene and that miR-148a targets the largest number of known PCa drivers. Globally, tissue-specific and stage-specific changes in miRNA targeting are driven by homeostatic response to active oncogenic pathways. Our findings indicate that, even in advanced PCa, the miRNA pool adapts to regulate continuing alterations in the cancer genome to balance oncogenic molecular changes. These findings are important because they are the first to globally characterize miRNA changes in PCa and demonstrate how the miRNA target spectrum responds to staged tumorigenesis.