In Planta Determination of the mRNA-Binding Proteome of Arabidopsis Etiolated Seedlings.

RNA binding proteins (RBPs) control the fate and expression of a transcriptome. Despite this fundamental importance, our understanding of plant RBPs is rudimentary, being mainly derived via bioinformatic extrapolation from other kingdoms. Here, we adapted the mRNA-protein interactome capture method to investigate the RNA binding proteome in planta. From Arabidopsis thaliana etiolated seedlings, we captured more than 700 proteins, including 300 with high confidence that we have defined as the At-RBP set. Approximately 75% of these At-RBPs are bioinformatically linked with RNA biology, containing a diversity of canonical RNA binding domains (RBDs). As no prior experimental RNA binding evidence exists for the majority of these proteins, their capture now authenticates them as RBPs. Moreover, we identified protein families harboring emerging and potentially novel RBDs, including WHIRLY, LIM, ALBA, DUF1296, and YTH domain-containing proteins, the latter being homologous to animal RNA methylation readers. Other At-RBP set proteins include major signaling proteins, cytoskeleton-associated proteins, membrane transporters, and enzymes, suggesting the scope and function of RNA-protein interactions within a plant cell is much broader than previously appreciated. Therefore, our foundation data set has provided an unbiased insight into the RNA binding proteome of plants, on which future investigations into plant RBPs can be based.

Optimization of signal-to-noise ratio for efficient microarray probe design.

MOTIVATION: Target-specific hybridization depends on oligo-probe characteristics that improve hybridization specificity and minimize genome-wide cross-hybridization. Interplay between specific hybridization and genome-wide cross-hybridization has been insufficiently studied, despite its crucial role in efficient probe design and in data analysis. RESULTS: In this study, we defined hybridization specificity as a ratio between oligo target-specific hybridization and oligo genome-wide cross-hybridization. A microarray database, derived from the Genomic Comparison Hybridization (GCH) experiment and performed using the Affymetrix platform, contains two different types of probes. The first type of oligo-probes does not have a specific target on the genome and their hybridization signals are derived from genome-wide cross-hybridization alone. The second type includes oligonucleotides that have a specific target on the genomic DNA and their signals are derived from specific and cross-hybridization components combined together in a total signal. A comparative analysis of hybridization specificity of oligo-probes, as well as their nucleotide sequences and thermodynamic features was performed on the database. The comparison has revealed that hybridization specificity was negatively affected by low stability of the fully-paired oligo-target duplex, stable probe self-folding, G-rich content, including GGG motifs, low sequence complexity and nucleotide composition symmetry. CONCLUSION: Filtering out the probes with defined 'negative' characteristics significantly increases specific hybridization and dramatically decreasing genome-wide cross-hybridization. Selected oligo-probes have two times higher hybridization specificity on average, compared to the probes that were filtered from the analysis by applying suggested cutoff thresholds to the described parameters. A new approach for efficient oligo-probe design is described in our study. CONTACT: shabalin@ncbi.nlm.nih.gov or olga.matveeva@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Hybridization-based reconstruction of small non-coding RNA transcripts from deep sequencing data.

Recent advances in RNA sequencing technology (RNA-Seq) enables comprehensive profiling of RNAs by producing millions of short sequence reads from size-fractionated RNA libraries. Although conventional tools for detecting and distinguishing non-coding RNAs (ncRNAs) from reference-genome data can be applied to sequence data, ncRNA detection can be improved by harnessing the full information content provided by this new technology. Here we present NorahDesk, the first unbiased and universally applicable method for small ncRNAs detection from RNA-Seq data. NorahDesk utilizes the coverage-distribution of small RNA sequence data as well as thermodynamic assessments of secondary structure to reliably predict and annotate ncRNA classes. Using publicly available mouse sequence data from brain, skeletal muscle, testis and ovary, we evaluated our method with an emphasis on the performance for microRNAs (miRNAs) and piwi-interacting small RNA (piRNA). We compared our method with Dario and mirDeep2 and found that NorahDesk produces longer transcripts with higher read coverage. This feature makes it the first method particularly suitable for the prediction of both known and novel piRNAs.

Quantitative prediction of miRNA-mRNA interaction based on equilibrium concentrations.

MicroRNAs (miRNAs) suppress gene expression by forming a duplex with a target messenger RNA (mRNA), blocking translation or initiating cleavage. Computational approaches have proven valuable for predicting which mRNAs can be targeted by a given miRNA, but currently available prediction methods do not address the extent of duplex formation under physiological conditions. Some miRNAs can at low concentrations bind to target mRNAs, whereas others are unlikely to bind within a physiologically relevant concentration range. Here we present a novel approach in which we find potential target sites on mRNA that minimize the calculated free energy of duplex formation, compute the free energy change involved in unfolding these sites, and use these energies to estimate the extent of duplex formation at specified initial concentrations of both species. We compare our predictions to experimentally confirmed miRNA-mRNA interactions (and non-interactions) in Drosophila melanogaster and in human. Although our method does not predict whether the targeted mRNA is degraded and/or its translation to protein inhibited, our quantitative estimates generally track experimentally supported results, indicating that this approach can be used to predict whether an interaction occurs at specified concentrations. Our approach offers a more-quantitative understanding of post-translational regulation in different cell types, tissues, and developmental conditions.

Insights into the biogenesis and potential functions of exonic circular RNA.

Circular RNAs (circRNAs) exhibit unique properties due to their covalently closed nature. Models of circRNAs synthesis and function are emerging but much remains undefined about this surprisingly prevalent class of RNA. Here, we identified exonic circRNAs from human and mouse RNA-sequencing datasets, documenting multiple new examples. Addressing function, we found that many circRNAs co-sediment with ribosomes, indicative of their translation potential. By contrast, circRNAs with potential to act as microRNA sponges were scarce, with some support for a collective sponge function by groups of circRNAs. Addressing circRNA biogenesis, we delineated several features commonly associated with circRNA occurrence. CircRNA-producing genes tend to be longer and to contain more exons than average. Back-splice acceptor exons are strongly enriched at ordinal position 2 within genes, and circRNAs typically have a short exon span with two exons being the most prevalent. The flanking introns either side of circRNA loci are exceptionally long. Of note also, single-exon circRNAs derive from unusually long exons while multi-exon circRNAs are mostly generated from exons of regular length. These findings independently validate and extend similar observations made in a number of prior studies. Furthermore, we analysed high-resolution RNA polymerase II occupancy data from two separate human cell lines to reveal distinctive transcription dynamics at circRNA-producing genes. Specifically, RNA polymerase II traverses the introns of these genes at above average speed concomitant with an accentuated slow-down at exons. Collectively, these features indicate how a perturbed balance between transcription and linear splicing creates important preconditions for circRNA production. We speculate that these preconditions need to be in place so that looping interactions between flanking introns can promote back-splicing to raise circRNA production to appreciable levels.

Small non-coding RNA expression from anterior cingulate cortex in schizophrenia shows sex specific regulation.

MicroRNAs (miRNAs) are known to regulate the expression of genes that are important for brain development and function, but the roles of other classes of small non-coding RNAs (sncRNAs) are less well understood. Additionally, although miRNA expression studies have been conducted in post-mortem brain samples from schizophrenia (SCZ) patients, other classes of sncRNAs are yet to be investigated in SCZ. We profiled the expression of miRNAs, piwi-interacting RNAs (piRNAs), small nucleolar RNAs (snoRNAs) and small nuclear RNAs (snRNAs) in SCZ by applying small RNA sequencing (RNA-Seq) to sncRNA isolated from post-mortem anterior cingulate cortex (ACC) of SCZ-affected individuals (n=22) and matched controls (n=22). We identified about one-third of annotated miRNAs, one-quarter of snoRNAs and a small proportion of piRNAs and snRNAs. No sncRNAs were significantly differentially expressed between SCZ and controls, but there was evidence for an interaction between disease status and sex on the expression level of a number of miRNAs and snoRNAs. Many of these transcripts exhibited differential expression between male and female cases, and/or between female cases and controls, suggesting sex based dysregulation in ACC of SCZ. These findings require replication in an independent sample, but our study provides further insights into the potential involvement of sncRNAs in brain function and SCZ.

The Cardiomyocyte RNA-Binding Proteome: Links to Intermediary Metabolism and Heart Disease.

RNA functions through the dynamic formation of complexes with RNA-binding proteins (RBPs) in all clades of life. We determined the RBP repertoire of beating cardiomyocytic HL-1 cells by jointly employing two in vivo proteomic methods, mRNA interactome capture and RBDmap. Together, these yielded 1,148 RBPs, 391 of which are shared with all other available mammalian RBP repertoires, while 393 are thus far unique to cardiomyocytes. RBDmap further identified 568 regions of RNA contact within 368 RBPs. The cardiomyocyte mRNA interactome composition reflects their unique biology. Proteins with roles in cardiovascular physiology or disease, mitochondrial function, and intermediary metabolism are all highly represented. Notably, we identified 73 metabolic enzymes as RBPs. RNA-enzyme contacts frequently involve Rossmann fold domains with examples in evidence of both, mutual exclusivity of, or compatibility between RNA binding and enzymatic function. Our findings raise the prospect of previously hidden RNA-mediated regulatory interactions among cardiomyocyte gene expression, physiology, and metabolism.

Long Noncoding RNA-Directed Epigenetic Regulation of Gene Expression Is Associated With Anxiety-like Behavior in Mice.

BACKGROUND: RNA-directed regulation of epigenetic processes has recently emerged as an important feature of mammalian differentiation and development. Perturbation of this regulatory system in the brain may contribute to the development of neuropsychiatric disorders. METHODS: RNA sequencing was used to identify changes in the experience-dependent expression of long noncoding RNAs (lncRNAs) within the medial prefrontal cortex of adult mice. Transcripts were validated by real-time quantitative polymerase chain reaction and a candidate lncRNA, Gomafu, was selected for further investigation. The functional role of this schizophrenia-related lncRNA was explored in vivo by antisense oligonucleotide-mediated gene knockdown in the medial prefrontal cortex, followed by behavioral training and assessment of fear-related anxiety. Long noncoding RNA-directed epigenetic regulation of gene expression was investigated by chromatin and RNA immunoprecipitation assays. RESULTS: RNA sequencing analysis revealed changes in the expression of a significant number of genes related to neural plasticity and stress, as well as the dynamic regulation of lncRNAs. In particular, we detected a significant downregulation of Gomafu lncRNA. Our results revealed that Gomafu plays a role in mediating anxiety-like behavior and suggest that this may occur through an interaction with a key member of the polycomb repressive complex 1, BMI1, which regulates the expression of the schizophrenia-related gene beta crystallin (Crybb1). We also demonstrated a novel role for Crybb1 in mediating fear-induced anxiety-like behavior. CONCLUSIONS: Experience-dependent expression of lncRNAs plays an important role in the epigenetic regulation of adaptive behavior, and the perturbation of Gomafu may be related to anxiety and the development of neuropsychiatric disorders.

Transcriptome-wide prediction of miRNA targets in human and mouse using FASTH.

Transcriptional regulation by microRNAs (miRNAs) involves complementary base-pairing at target sites on mRNAs, yielding complex secondary structures. Here we introduce an efficient computational approach and software (FASTH) for genome-scale prediction of miRNA target sites based on minimizing the free energy of duplex structure. We apply our approach to identify miRNA target sites in the human and mouse transcriptomes. Our results show that short sequence motifs in the 5' end of miRNAs frequently match mRNAs perfectly, not only at validated target sites but additionally at many other, energetically favourable sites. High-quality matching regions are abundant and occur at similar frequencies in all mRNA regions, not only the 3'UTR. About one-third of potential miRNA target sites are reassigned to different mRNA regions, or gained or lost altogether, among different transcript isoforms from the same gene. Many potential miRNA target sites predicted in human are not found in mouse, and vice-versa, but among those that do occur in orthologous human and mouse mRNAs most are situated in corresponding mRNA regions, i.e. these sites are themselves orthologous. Using a luciferase assay in HEK293 cells, we validate four of six predicted miRNA-mRNA interactions, with the mRNA level reduced by an average of 73%. We demonstrate that a thermodynamically based computational approach to prediction of miRNA binding sites on mRNAs can be scaled to analyse complete mammalian transcriptome datasets. These results confirm and extend the scope of miRNA-mediated species- and transcript-specific regulation in different cell types, tissues and developmental conditions.