Systemic CLIP-seq analysis and game theory approach to model microRNA mode of binding

microRNAs (miRNAs) associate with Ago proteins to post-transcriptionally silence gene expression by targeting mRNAs. To characterize the modes of miRNA-binding, we developed a novel computational framework, called optiCLIP, which considers the reproducibility of the identified peaks among replicates based on the peak overlap. We identified 98999 binding sites for mouse and human miRNAs, from eleven Ago2 CLIP-seq datasets. Clustering the binding preferences, we found heterogeneity of the mode of binding for different miRNAs. Finally, we set up a quantitative model, named miRgame, based on an adaptation of the game theory. We have developed a new algorithm to translate the miRgame into a score that corresponds to a miRNA degree of occupancy for each Ago2 peak. The degree of occupancy summarizes the number of miRNA-binding sites and miRNAs targeting each binding site, and binding energy of each miRNA::RNA heteroduplex in each peak. Ago peaks were stratified accordingly to the degree of occupancy. Target repression correlates with higher score of degree of occupancy and number of miRNA-binding sites within each Ago peak. We validated the biological performance of our new method on miR-155-5p. In conclusion, our data demonstrate that miRNA-binding sites within each Ago2 CLIP-seq peak synergistically interplay to enhance target repression.

Automated summary

The study established a new computational framework, optiCLIP, to characterize the modes of miRNA-binding and their role in gene silencing. Clustering the binding preferences revealed heterogeneity in the modes of miRNA binding. A quantitative model, miRgame, was developed to demonstrate the relationships between miRNA binding sites and miRNAs.