Abstract:Background RNA-binding proteins (RBPs) play diverse roles in eukaryotic RNA processing. Despite their pervasive functions in coding and noncoding RNA biogenesis and regulation, elucidating the sequence specificities that define protein-RNA interactions remains a major challenge. Recently, CLIP-seq (Cross-linking immunoprecipitation followed by high-throughput sequencing) has been successfully implemented to study the transcriptome-wide binding patterns of SRSF1, PTBP1, NOVA and fox2 proteins. These studies either adopted traditional methods like Multiple EM for Motif Elicitation (MEME) to discover the sequence consensus of RBP's binding sites or used Z-score statistics to search for the overrepresented nucleotides of a certain size. We argue that most of these methods are not well-suited for RNA motif identification, as they are unable to incorporate the RNA structural context of protein-RNA interactions, which may affect to binding specificity. Here, we describe a novel model-based approach-- RNAMotifModeler to identify the consensus of protein-RNA binding regions by integrating sequence features and RNA secondary structures. Results As an example, we implemented RNAMotifModeler on SRSF1 (SF2/ASF) CLIP-seq data. The sequence-structural consensus we identified is a purine-rich octamer 'AGAAGAAG' in a highly single-stranded RNA context. The unpaired probabilities, the probabilities of not forming pairs, are significantly higher than negative controls and the flanking sequence surrounding the binding site, indicating that SRSF1 proteins tend to bind on single-stranded RNA. Further statistical evaluations revealed that the second and fifth bases of SRSF1octamer motif have much stronger sequence specificities, but weaker single-strandedness, while the third, fourth, sixth and seventh bases are far more likely to be single-stranded, but have more degenerate sequence specificities. Therefore, we hypothesize that nucleotide specificity and secondary structure play complementary roles during binding site recognition by SRSF1. Conclusion In this study, we presented a computational model to predict the sequence consensus and optimal RNA secondary structure for protein-RNA binding regions. The successful implementation on SRSF1 CLIP-seq data demonstrates great potential to improve our understanding on the binding specificity of RNA binding proteins.

Identification of SRSF3 target mRNAs using inducible TRIBE

Predicting Sequence and Structural Specificities of RNA Binding Regions Recognized by Splicing Factor SRSF1

Serine and Arginine-Rich Splicing Factor 3 Promotes the Activation of Quiescent Mouse Neural Stem Cells.

Genome-wide analysis reveals SR protein cooperation and competition in regulated splicing.

A Genome Landscape of SRSF3-regulated Splicing Events and Gene Expression in Human Osteosarcoma U2OS Cells.

Revisiting and Corrections to the Annotated SRSF3 (srp20) Gene Structure and RefSeq Sequences from the Human and Mouse Genomes

Oncogenic SRSF3 in Health and Diseases.

Systematic characterization of short intronic splicing-regulatory elements

Single-cell profiling identifies LIN28A mRNA targets in the mouse pluripotent-to-2C-like transition and somatic cell reprogramming

Systematic characterization of short intronic splicing-regulatory elements in SMN2 pre-mRNA

Rbm24 Regulates Alternative Splicing Switch in Embryonic Stem Cell Cardiac Lineage Differentiation.

Emerging Roles of SRSF3 As a Therapeutic Target for Cancer

Defining The Regulatory Network Of The Tissue-Specific Splicing Factors Fox-1 And Fox-2

Factors affecting the expression and stability of full-length and truncated SRSF3 proteins in human cancer cells

Integrated Analysis of a Compendium of RNA-Seq Datasets for Splicing Factors.

Systematic Identification of Cis-Regulatory Sequences Active in Mouse and Human Embryonic Stem Cells

Computational Identification of Tissue-Specific Splicing Regulatory Elements in Human Genes from RNA-Seq Data

SRSF3 and SRSF7 modulate 3′UTR length through suppression or activation of proximal polyadenylation sites and regulation of CFIm levels

Towards Understandings of Serine/arginine-Rich Splicing Factors

A Trna-Derived Small RNA Regulates Ribosomal Protein S28 Protein Levels after Translation Initiation in Humans and Mice.

Transcriptome analysis of alternative splicing events regulated by SRSF10 reveals position-dependent splicing modulation