Toutai Mituyama

MENε/β noncoding RNAs are essential for structural integrity of nuclear paraspeckles

Recent transcriptome analyses have shown that thousands of noncoding RNAs (ncRNAs) are transcribed from mammalian genomes. Although the number of functionally annotated ncRNAs is still limited, they are known to be frequently retained in the nucleus, where they coordinate regulatory networks of gene expression. Some subnuclear organelles or nuclear bodies include RNA species whose identity and structural roles are largely unknown. We identified 2 abundant overlapping ncRNAs, MENε and MENβ (MENε/β), which are transcribed from the corresponding site in the multiple endocrine neoplasia (MEN) I locus and which localize to nuclear paraspeckles. This finding raises the intriguing possibility that MENε/β are involved in paraspeckle organization, because paraspeckles are, reportedly, RNase-sensitive structures. Successful removal of MENε/β by a refined knockdown method resulted in paraspeckle disintegration. Furthermore, the reassembly of paraspeckles disassembled by transcriptional arrest appeared to be unsuccessful in the absence of MENε/β. RNA interference and immunoprecipitation further revealed that the paraspeckle proteins p54/nrb and PSF selectively associate with and stabilize the longer MENβ, thereby contributing to the organization of the paraspeckle structure. The paraspeckle protein PSP1 is not directly involved in either MENε/β stabilization or paraspeckle organization. We postulate a model for nuclear paraspeckle body organization where specific ncRNAs and RNA-binding proteins cooperate to maintain and, presumably, establish the structure.

A regulatory circuit for piwi by the large Maf gene traffic jam in Drosophila

PIWI-interacting RNAs (piRNAs) silence retrotransposons in Drosophila germ lines by associating with the PIWI proteins Argonaute 3 (AGO3), Aubergine (Aub) and Piwi. piRNAs in Drosophila are produced from intergenic repetitive genes and piRNA clusters by two systems: the primary processing pathway and the amplification loop. The amplification loop occurs in a Dicer-independent, PIWI-Slicer-dependent manner. However, primary piRNA processing remains elusive. Here we analysed piRNA processing in a Drosophila ovarian somatic cell line where Piwi, but not Aub or AGO3, is expressed; thus, only the primary piRNAs exist. In addition to flamenco, a Piwi-specific piRNA cluster, traffic jam (tj), a large Maf gene, was determined as a new piRNA cluster. piRNAs arising from tj correspond to the untranslated regions of tj messenger RNA and are sense-oriented. piRNA loading on to Piwi may occur in the cytoplasm. zucchini, a gene encoding a putative cytoplasmic nuclease, is required for tj-derived piRNA production. In tj and piwi mutant ovaries, somatic cells fail to intermingle with germ cells and Fasciclin III is overexpressed. Loss of tj abolishes Piwi expression in gonadal somatic cells. Thus, in gonadal somatic cells, tj gives rise simultaneously to two different molecules: the TJ protein, which activates Piwi expression, and piRNAs, which define the Piwi targets for silencing.

Characterization of endogenous human Argonautes and their miRNA partners in RNA silencing

Small RNAs triggering RNA silencing are loaded onto Argonautes and then sequence-specifically guide them to target transcripts. Epitope-tagged human Argonautes (hAgo1, hAgo2, hAgo3, and hAgo4) are associated with siRNAs and miRNAs, but only epitope-tagged hAgo2 has been shown to have Slicer activity. Contrarily, how endogenous hAgos behave with respect to small RNA association and target RNA destruction has remained unclear. Here, we produced monoclonal antibodies for individual hAgos. High-throughput pyrosequencing revealed that immunopurified endogenous hAgo2 and hAgo3 associated mostly with miRNAs. Endogenous hAgo3 did not show Slicer function but localized in P-bodies, suggesting that hAgo3 endogenously expressed is, like hAgo2, involved in the miRNA pathway but antagonizes the RNAi activity of hAgo2. Sequence variations of miRNAs were found at both 5′ and 3′ ends, suggesting that multiple mature miRNAs containing different “seed” sequences can arise from one miRNA precursor. The hAgo antibodies we raised are valuable tools for ascertaining the functional behavior of endogenous Argonautes and miRNAs in RNA silencing.

CentroidFold: a web server for RNA secondary structure prediction

The CentroidFold web server (http://www.ncrna.org/centroidfold/) is a web application for RNA secondary structure prediction powered by one of the most accurate prediction engine. The server accepts two kinds of sequence data: a single RNA sequence and a multiple alignment of RNA sequences. It responses with a prediction result shown as a popular base-pair notation and a graph representation. PDF version of the graph representation is also available. For a multiple alignment sequence, the server predicts a common secondary structure. Usage of the server is quite simple. You can paste a single RNA sequence (FASTA or plain sequence text) or a multiple alignment (CLUSTAL-W format) into the textarea then click on the ‘execute CentroidFold’ button. The server quickly responses with a prediction result. The major advantage of this server is that it employs our original CentroidFold software as its prediction engine which scores the best accuracy in our benchmark results. Our web server is freely available with no login requirement.

Prediction of RNA secondary structure using generalized centroid estimators

MOTIVATION Recent studies have shown that the methods for predicting secondary structures of RNAs on the basis of posterior decoding of the base-pairing probabilities has an advantage with respect to prediction accuracy over the conventionally utilized minimum free energy methods. However, there is room for improvement in the objective functions presented in previous studies, which are maximized in the posterior decoding with respect to the accuracy measures for secondary structures. RESULTS We propose novel estimators which improve the accuracy of secondary structure prediction of RNAs. The proposed estimators maximize an objective function which is the weighted sum of the expected number of the true positives and that of the true negatives of the base pairs. The proposed estimators are also improved versions of the ones used in previous works, namely CONTRAfold for secondary structure prediction from a single RNA sequence and McCaskill-MEA for common secondary structure prediction from multiple alignments of RNA sequences. We clarify the relations between the proposed estimators and the estimators presented in previous works, and theoretically show that the previous estimators include additional unnecessary terms in the evaluation measures with respect to the accuracy. Furthermore, computational experiments confirm the theoretical analysis by indicating improvement in the empirical accuracy. The proposed estimators represent extensions of the centroid estimators proposed in Ding et al. and Carvalho and Lawrence, and are applicable to a wide variety of problems in bioinformatics. AVAILABILITY Supporting information and the CentroidFold software are available online at: http://www.ncrna.org/software/centroidfold/.

Functional involvement of Tudor and dPRMT5 in the piRNA processing pathway in Drosophila germlines

In Drosophila, the PIWI proteins, Aubergine (Aub), AGO3, and Piwi are expressed in germlines and function in silencing transposons by associating with PIWI‐interacting RNAs (piRNAs). Recent studies show that PIWI proteins contain symmetric dimethyl‐arginines (sDMAs) and that dPRMT5/Capsuleen/DART5 is the modifying enzyme. Here, we show that Tudor (Tud), one of Tud domain‐containing proteins, associates with Aub and AGO3, specifically through their sDMA modifications and that these three proteins form heteromeric complexes. piRNA precursor‐like molecules are detected in these complexes. The expression levels of Aub and AGO3, along with their degree of sDMA modification, were not changed by tud mutations. However, the population of transposon‐derived piRNAs associated with Aub and AGO3 was altered by tud mutations, whereas the total amounts of small RNAs on Aub and AGO3 was increased. Loss of dprmt5 did not change the stability of Aub, but impaired its association with Tud and lowered piRNA association with Aub. Thus, in germline cells, piRNAs are quality‐controlled by dPRMT5 that modifies PIWI proteins, in tight association with Tud.

The Functional RNA Database 3.0: databases to support mining and annotation of functional RNAs

We developed a pair of databases that support two important tasks: annotation of anonymous RNA transcripts and discovery of novel non-coding RNAs. The database combo is called the Functional RNA Database and consists of two databases: a rewrite of the original version of the Functional RNA Database (fRNAdb) and the latest version of the UCSC GenomeBrowser for Functional RNA. The former is a sequence database equipped with a powerful search function and hosts a large collection of known/predicted non-coding RNA sequences acquired from existing databases as well as novel/predicted sequences reported by researchers of the Functional RNA Project. The latter is a UCSC Genome Browser mirror with large additional custom tracks specifically associated with non-coding elements. It also includes several functional enhancements such as a presentation of a common secondary structure prediction at any given genomic window ⩽500 bp. Our GenomeBrowser supports user authentication and user-specific tracks. The current version of the fRNAdb is a complete rewrite of the former version, hosting a larger number of sequences and with a much friendlier interface. The current version of UCSC GenomeBrowser for Functional RNA features a larger number of tracks and richer features than the former version. The databases are available at http://www.ncrna.org/.

CentroidAlign: fast and accurate aligner for structured RNAs by maximizing expected sum-of-pairs score.

MOTIVATION The importance of accurate and fast predictions of multiple alignments for RNA sequences has increased due to recent findings about functional non-coding RNAs. Recent studies suggest that maximizing the expected accuracy of predictions will be useful for many problems in bioinformatics. RESULTS We designed a novel estimator for multiple alignments of structured RNAs, based on maximizing the expected accuracy of predictions. First, we define the maximum expected accuracy (MEA) estimator for pairwise alignment of RNA sequences. This maximizes the expected sum-of-pairs score (SPS) of a predicted alignment under a probability distribution of alignments given by marginalizing the Sankoff model. Then, by approximating the MEA estimator, we obtain an estimator whose time complexity is O(L(3)+c(2)dL(2)) where L is the length of input sequences and both c and d are constants independent of L. The proposed estimator can handle uncertainty of secondary structures and alignments that are obstacles in Bioinformatics because it considers all the secondary structures and all the pairwise alignments as input sequences. Moreover, we integrate the probabilistic consistency transformation (PCT) on alignments into the proposed estimator. Computational experiments using six benchmark datasets indicate that the proposed method achieved a favorable SPS and was the fastest of many state-of-the-art tools for multiple alignments of structured RNAs. AVAILABILITY The software called CentroidAlign, which is an implementation of the algorithm in this article, is freely available on our website: http://www.ncrna.org/software/centroidalign/. CONTACT hamada-michiaki@aist.go.jp SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Bisulfighter: accurate detection of methylated cytosines and differentially methylated regions.

Analysis of bisulfite sequencing data usually requires two tasks: to call methylated cytosines (mCs) in a sample, and to detect differentially methylated regions (DMRs) between paired samples. Although numerous tools have been proposed for mC calling, methods for DMR detection have been largely limited. Here, we present Bisulfighter, a new software package for detecting mCs and DMRs from bisulfite sequencing data. Bisulfighter combines the LAST alignment tool for mC calling, and a novel framework for DMR detection based on hidden Markov models (HMMs). Unlike previous attempts that depend on empirical parameters, Bisulfighter can use the expectation-maximization algorithm for HMMs to adjust parameters for each data set. We conduct extensive experiments in which accuracy of mC calling and DMR detection is evaluated on simulated data with various mC contexts, read qualities, sequencing depths and DMR lengths, as well as on real data from a wide range of biological processes. We demonstrate that Bisulfighter consistently achieves better accuracy than other published tools, providing greater sensitivity for mCs with fewer false positives, more precise estimates of mC levels, more exact locations of DMRs and better agreement of DMRs with gene expression and DNase I hypersensitivity. The source code is available at http://epigenome.cbrc.jp/bisulfighter.

Predictions of RNA secondary structure by combining homologous sequence information

Motivation: Secondary structure prediction of RNA sequences is an important problem. There have been progresses in this area, but the accuracy of prediction from an RNA sequence is still limited. In many cases, however, homologous RNA sequences are available with the target RNA sequence whose secondary structure is to be predicted. Results: In this article, we propose a new method for secondary structure predictions of individual RNA sequences by taking the information of their homologous sequences into account without assuming the common secondary structure of the entire sequences. The proposed method is based on posterior decoding techniques, which consider all the suboptimal secondary structures of the target and homologous sequences and all the suboptimal alignments between the target sequence and each of the homologous sequences. In our computational experiments, the proposed method provides better predictions than those performed only on the basis of the formation of individual RNA sequences and those performed by using methods for predicting the common secondary structure of the homologous sequences. Remarkably, we found that the common secondary predictions sometimes give worse predictions for the secondary structure of a target sequence than the predictions from the individual target sequence, while the proposed method always gives good predictions for the secondary structure of target sequences in all tested cases. Availability: Supporting information and software are available online at: http://www.ncrna.org/software/centroidfold/ismb2009/. Contact: hamada-michiaki@aist.go.jp Supplementary information:Supplementary data are available at Bioinformatics online.