GraphProt: modeling binding preferences of RNA-binding proteins

We present GraphProt, a computational framework for learning sequence- and structure-binding preferences of RNA-binding proteins (RBPs) from high-throughput experimental data. We benchmark GraphProt, demonstrating that the modeled binding preferences conform to the literature, and showcase the biological relevance and two applications of GraphProt models. First, estimated binding affinities correlate with experimental measurements. Second, predicted Ago2 targets display higher levels of expression upon Ago2 knockdown, whereas control targets do not. Computational binding models, such as those provided by GraphProt, are essential for predicting RBP binding sites and affinities in all tissues. GraphProt is freely available at http://www.bioinf.uni-freiburg.de/Software/GraphProt.     

RNAstrand: reading direction of structured RNAs in multiple sequence alignments

Motivation

Genome-wide screens for structured ncRNA genes in mammals, urochordates, and nematodes have predicted thousands of putative ncRNA genes and other structured RNA motifs. A prerequisite for their functional annotation is to determine the reading direction with high precision.

Results

While folding energies of an RNA and its reverse complement are similar, the differences are sufficient at least in conjunction with substitution patterns to discriminate between structured RNAs and their complements. We present here a support vector machine that reliably classifies the reading direction of a structured RNA from a multiple sequence alignment and provides a considerable improvement in classification accuracy over previous approaches.

Software

RNAstrand is freely available as a stand-alone tool from http://www.bioinf.uni-leipzig.de/Software/RNAstrand and is also included in the latest release of RNAz, a part of the Vienna RNA Package.     

Enhancing protein-vitamin binding residues prediction by multiple heterogeneous subspace SVMs ensemble

Background

Vitamins are typical ligands that play critical roles in various metabolic processes. The accurate identification of the vitamin-binding residues solely based on a protein sequence is of significant importance for the functional annotation of proteins, especially in the post-genomic era, when large volumes of protein sequences are accumulating quickly without being functionally annotated.

Results

In this paper, a new predictor called TargetVita is designed and implemented for predicting protein-vitamin binding residues using protein sequences. In TargetVita, features derived from the position-specific scoring matrix (PSSM), predicted protein secondary structure, and vitamin binding propensity are combined to form the original feature space; then, several feature subspaces are selected by performing different feature selection methods. Finally, based on the selected feature subspaces, heterogeneous SVMs are trained and then ensembled for performing prediction.

Conclusions

The experimental results obtained with four separate vitamin-binding benchmark datasets demonstrate that the proposed TargetVita is superior to the state-of-the-art vitamin-specific predictor, and an average improvement of 10% in terms of the Matthews correlation coefficient (MCC) was achieved over independent validation tests. The TargetVita web server and the datasets used are freely available for academic use at http://csbio.njust.edu.cn/bioinf/TargetVita or http://www.csbio.sjtu.edu.cn/bioinf/TargetVita.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-297) contains supplementary material, which is available to authorized users.     

A low-latency,big database system and browser for storage,querying and visualization of 3D genomic data

Recent releases of genome three-dimensional (3D) structures have the potential to transform our understanding of genomes. Nonetheless, the storage technology and visualization tools need to evolve to offer to the scientific community fast and convenient access to these data. We introduce simultaneously a database system to store and query 3D genomic data (3DBG), and a 3D genome browser to visualize and explore 3D genome structures (3DGB). We benchmark 3DBG against state-of-the-art systems and demonstrate that it is faster than previous solutions, and importantly gracefully scales with the size of data. We also illustrate the usefulness of our 3D genome Web browser to explore human genome structures. The 3D genome browser is available at http://3dgb.cs.mcgill.ca/.     

RNAlien – Unsupervised RNA family model construction

Determining the function of a non-coding RNA requires costly and time-consuming wet-lab experiments. For this reason, computational methods which ascertain the homology of a sequence and thereby deduce functionality and family membership are often exploited. In this fashion, newly sequenced genomes can be annotated in a completely computational way. Covariance models are commonly used to assign novel RNA sequences to a known RNA family. However, to construct such models several examples of the family have to be already known. Moreover, model building is the work of experts who manually edit the necessary RNA alignment and consensus structure. Our method, RNAlien, starting from a single input sequence collects potential family member sequences by multiple iterations of homology search. RNA family models are fully automatically constructed for the found sequences. We have tested our method on a subset of the RfamRNA family database. RNAlien models are a starting point to construct models of comparable sensitivity and specificity to manually curated ones from the Rfam database. RNAlien Tool and web server are available at http://rna.tbi.univie.ac.at/rnalien/.     

Complete RNA inverse folding: computational design of functional hammerhead ribozymes

Nanotechnology and synthetic biology currently constitute one of the most innovative, interdisciplinary fields of research, poised to radically transform society in the 21st century. This paper concerns the synthetic design of ribonucleic acid molecules, using our recent algorithm, RNAiFold, which can determine all RNA sequences whose minimum free energy secondary structure is a user-specified target structure. Using RNAiFold, we design ten cis-cleaving hammerhead ribozymes, all of which are shown to be functional by a cleavage assay. We additionally use RNAiFold to design a functional cis-cleaving hammerhead as a modular unit of a synthetic larger RNA. Analysis of kinetics on this small set of hammerheads suggests that cleavage rate of computationally designed ribozymes may be correlated with positional entropy, ensemble defect, structural flexibility/rigidity and related measures. Artificial ribozymes have been designed in the past either manually or by SELEX (Systematic Evolution of Ligands by Exponential Enrichment); however, this appears to be the first purely computational design and experimental validation of novel functional ribozymes. RNAiFold is available at http://bioinformatics.bc.edu/clotelab/RNAiFold/.     

A New Supervised Over-Sampling Algorithm with Application to Protein-Nucleotide Binding Residue Prediction

Protein-nucleotide interactions are ubiquitous in a wide variety of biological processes. Accurately identifying interaction residues solely from protein sequences is useful for both protein function annotation and drug design, especially in the post-genomic era, as large volumes of protein data have not been functionally annotated. Protein-nucleotide binding residue prediction is a typical imbalanced learning problem, where binding residues are extremely fewer in number than non-binding residues. Alleviating the severity of class imbalance has been demonstrated to be a promising means of improving the prediction performance of a machine-learning-based predictor for class imbalance problems. However, little attention has been paid to the negative impact of class imbalance on protein-nucleotide binding residue prediction. In this study, we propose a new supervised over-sampling algorithm that synthesizes additional minority class samples to address class imbalance. The experimental results from protein-nucleotide interaction datasets demonstrate that the proposed supervised over-sampling algorithm can relieve the severity of class imbalance and help to improve prediction performance. Based on the proposed over-sampling algorithm, a predictor, called TargetSOS, is implemented for protein-nucleotide binding residue prediction. Cross-validation tests and independent validation tests demonstrate the effectiveness of TargetSOS. The web-server and datasets used in this study are freely available at http://www.csbio.sjtu.edu.cn/bioinf/TargetSOS/.     

Translational Control of the Oogenic Program by Components of OMA Ribonucleoprotein Particles in Caenorhabditis elegans

The oocytes of most sexually reproducing animals arrest in meiotic prophase I. Oocyte growth, which occurs during this period of arrest, enables oocytes to acquire the cytoplasmic components needed to produce healthy progeny and to gain competence to complete meiosis. In the nematode Caenorhabditis elegans, the major sperm protein hormone promotes meiotic resumption (also called meiotic maturation) and the cytoplasmic flows that drive oocyte growth. Prior work established that two related TIS11 zinc-finger RNA-binding proteins, OMA-1 and OMA-2, are redundantly required for normal oocyte growth and meiotic maturation. We affinity purified OMA-1 and identified associated mRNAs and proteins using genome-wide expression data and mass spectrometry, respectively. As a class, mRNAs enriched in OMA-1 ribonucleoprotein particles (OMA RNPs) have reproductive functions. Several of these mRNAs were tested and found to be targets of OMA-1/2-mediated translational repression, dependent on sequences in their 3′-untranslated regions (3′-UTRs). Consistent with a major role for OMA-1 and OMA-2 in regulating translation, OMA-1-associated proteins include translational repressors and activators, and some of these proteins bind directly to OMA-1 in yeast two-hybrid assays, including OMA-2. We show that the highly conserved TRIM-NHL protein LIN-41 is an OMA-1-associated protein, which also represses the translation of several OMA-1/2 target mRNAs. In the accompanying article in this issue, we show that LIN-41 prevents meiotic maturation and promotes oocyte growth in opposition to OMA-1/2. Taken together, these data support a model in which the conserved regulators of mRNA translation LIN-41 and OMA-1/2 coordinately control oocyte growth and the proper spatial and temporal execution of the meiotic maturation decision.     

Integrity and Function of the Saccharomyces cerevisiae Spindle Pole Body Depends on Connections Between the Membrane Proteins Ndc1, Rtn1, and Yop1

The nuclear envelope in Saccharomyces cerevisiae harbors two essential macromolecular protein assemblies: the nuclear pore complexes (NPCs) that enable nucleocytoplasmic transport, and the spindle pole bodies (SPBs) that mediate chromosome segregation. Previously, based on metazoan and budding yeast studies, we reported that reticulons and Yop1/DP1 play a role in the early steps of de novo NPC assembly. Here, we examined if Rtn1 and Yop1 are required for SPB function in S. cerevisiae. Electron microscopy of rtn1Δ yop1Δ cells revealed lobular abnormalities in SPB structure. Using an assay that monitors lateral expansion of the SPB central layer, we found that rtn1Δ yop1Δ SPBs had decreased connections to the NE compared to wild type, suggesting that SPBs are less stable in the NE. Furthermore, large budded rtn1Δ yop1Δ cells exhibited a high incidence of short mitotic spindles, which were frequently misoriented with respect to the mother–daughter axis. This correlated with cytoplasmic microtubule defects. We found that overexpression of the SPB insertion factors NDC1, MPS2, or BBP1 rescued the SPB defects observed in rtn1Δ yop1Δ cells. However, only overexpression of NDC1, which is also required for NPC biogenesis, rescued both the SPB and NPC associated defects. Rtn1 and Yop1 also physically interacted with Ndc1 and other NPC membrane proteins. We propose that NPC and SPB biogenesis are altered in cells lacking Rtn1 and Yop1 due to competition between these complexes for Ndc1, an essential common component of both NPCs and SPBs.     

Online database for mosquito (Diptera,Culicidae) occurrence records in French Guiana

A database providing information on mosquito specimens (Arthropoda: Diptera: Culicidae) collected in French Guiana is presented. Field collections were initiated in 2013 under the auspices of the CEnter for the study of Biodiversity in Amazonia (CEBA: http://www.labexceba.fr/en/). This study is part of an ongoing process aiming to understand the distribution of mosquitoes, including vector species, across French Guiana. Occurrences are recorded after each collecting trip in a database managed by the laboratory Evolution et Diversité Biologique (EDB), Toulouse, France. The dataset is updated monthly and is available online. Voucher specimens and their associated DNA are stored at the laboratory Ecologie des Forêts de Guyane (Ecofog), Kourou, French Guiana. The latest version of the dataset is accessible through EDB’s Integrated Publication Toolkit at http://130.120.204.55:8080/ipt/resource.do?r=mosquitoes_of_french_guiana or through the Global Biodiversity Information Facility data portal at http://www.gbif.org/dataset/5a8aa2ad-261c-4f61-a98e-26dd752fe1c5 It can also be viewed through the Guyanensis platform at http://guyanensis.ups-tlse.fr     

leeHom: adaptor trimming and merging for Illumina sequencing reads

The sequencing of libraries containing molecules shorter than the read length, such as in ancient or forensic applications, may result in the production of reads that include the adaptor, and in paired reads that overlap one another. Challenges for the processing of such reads are the accurate identification of the adaptor sequence and accurate reconstruction of the original sequence most likely to have given rise to the observed read(s). We introduce an algorithm that removes the adaptors and reconstructs the original DNA sequences using a Bayesian maximum a posteriori probability approach. Our algorithm is faster, and provides a more accurate reconstruction of the original sequence for both simulated and ancient DNA data sets, than other approaches. leeHom is released under the GPLv3 and is freely available from: https://bioinf.eva.mpg.de/leehom/     

The Abundant Histone Chaperones Spt6 and FACT Collaborate to Assemble,Inspect, and Maintain Chromatin Structure in Saccharomyces cerevisiae

Saccharomyces cerevisiae Spt6 protein is a conserved chromatin factor with several distinct functional domains, including a natively unstructured 30-residue N-terminal region that binds competitively with Spn1 or nucleosomes. To uncover physiological roles of these interactions, we isolated histone mutations that suppress defects caused by weakening Spt6:Spn1 binding with the spt6-F249K mutation. The strongest suppressor was H2A-N39K, which perturbs the point of contact between the two H2A-H2B dimers in an assembled nucleosome. Substantial suppression also was observed when the H2A-H2B interface with H3-H4 was altered, and many members of this class of mutations also suppressed a defect in another essential histone chaperone, FACT. Spt6 is best known as an H3-H4 chaperone, but we found that it binds with similar affinity to H2A-H2B or H3-H4. Like FACT, Spt6 is therefore capable of binding each of the individual components of a nucleosome, but unlike FACT, Spt6 did not produce endonuclease-sensitive reorganized nucleosomes and did not displace H2A-H2B dimers from nucleosomes. Spt6 and FACT therefore have distinct activities, but defects can be suppressed by overlapping histone mutations. We also found that Spt6 and FACT together are nearly as abundant as nucleosomes, with ∼24,000 Spt6 molecules, ∼42,000 FACT molecules, and ∼75,000 nucleosomes per cell. Histone mutations that destabilize interfaces within nucleosomes therefore reveal multiple spatial regions that have both common and distinct roles in the functions of these two essential and abundant histone chaperones. We discuss these observations in terms of different potential roles for chaperones in both promoting the assembly of nucleosomes and monitoring their quality.     

ChiloKey,an interactive identification tool for the geophilomorph centipedes of Europe (Chilopoda,Geophilomorpha)

ChiloKey is a matrix-based, interactive key to all 179 species of Geophilomorpha (Chilopoda) recorded from Europe, including species of uncertain identity and those whose morphology is known partially only. The key is intended to assist in identification of subadult and adult specimens, by means of microscopy and simple dissection techniques whenever necessary. The key is freely available through the web at: http://www.biologia.unipd.it/chilokey/ and at http://www.interactive-keys.eu/chilokey/.