A max-margin model for efficient simultaneous alignment and folding of RNA sequences

Motivation: The need for accurate and efficient tools for computational RNA structure analysis has become increasingly apparent over the last several years: RNA folding algorithms underlie numerous applications in bioinformatics, ranging from microarray probe selection to de novo non-coding RNA gene prediction

Characterization of eight novel proteins with male germ cell-specific expression in mouse

<p>Abstract</p> <p>Background</p> <p>Spermatogenesis and fertilization are highly unique processes. Discovery and characterization of germ cell-specific genes are important for the understanding of these reproductive processes. We investigated eight proteins encoded by novel spermatogenic cell-specific genes previously identified from the mouse round spermatid UniGene library.</p> <p>Methods</p> <p>Polyclonal antibodies were generated against the novel proteins and western blot analysis was performed with various protein samples. Germ cell specificity was investigated using testes from germ cell-less mutant mice. Developmental expression pattern was examined in testicular germ cells, testicular sperm and mature sperm. Subcellular localization was assessed by cell surface biotin labeling and trypsinization. Protein localization and properties in sperm were investigated by separation of head and tail fractions, and extractabilities by a non-ionic detergent and urea.</p> <p>Results</p> <p>The authenticity of the eight novel proteins and their specificity to spermatogenic cells were confirmed. In examining the developmental expression patterns, we found the presence of four proteins only in testicular germ cells, a single protein in testicular germ cells and testicular sperm, and three proteins in the testicular stages and mature sperm from the epididymis. Further analysis of the three proteins present in sperm disclosed that one is located at the surface of the acrosomal region and the other two are associated with cytoskeletal structures in the sperm flagellum. We name the genes for these sperm proteins Shsp1 (Sperm head surface protein 1), Sfap1 (Sperm flagellum associated protein 1) and Sfap2 (Sperm flagellum associated protein 2).</p> <p>Conclusion</p> <p>We analyzed eight novel germ cell-specific proteins, providing new and inclusive information about their developmental and cellular characteristics. Our findings will facilitate future investigation into the biological roles of these novel proteins in spermatogenesis and sperm functions.</p

Integrative characterization of germ cell-specific genes from mouse spermatocyte UniGene library

<p>Abstract</p> <p>Background</p> <p>The primary regulator of spermatogenesis, a highly ordered and tightly regulated developmental process, is an intrinsic genetic program involving male germ cell-specific genes.</p> <p>Results</p> <p>We analyzed the mouse spermatocyte UniGene library containing 2155 gene-oriented transcript clusters. We predict that 11% of these genes are testis-specific and systematically identified 24 authentic genes specifically and abundantly expressed in the testis via in <it>silico </it>and <it>in vitro </it>approaches. Northern blot analysis disclosed various transcript characteristics, such as expression level, size and the presence of isoform. Expression analysis revealed developmentally regulated and stage-specific expression patterns in all of the genes. We further analyzed the genes at the protein and cellular levels. Transfection assays performed using GC-2 cells provided information on the cellular characteristics of the gene products. In addition, antibodies were generated against proteins encoded by some of the genes to facilitate their identification and characterization in spermatogenic cells and sperm. Our data suggest that a number of the gene products are implicated in transcriptional regulation, nuclear integrity, sperm structure and motility, and fertilization. In particular, we found for the first time that Mm.333010, predicted to contain a trypsin-like serine protease domain, is a sperm acrosomal protein.</p> <p>Conclusion</p> <p>We identify 24 authentic genes with spermatogenic cell-specific expression, and provide comprehensive information about the genes. Our findings establish a new basis for future investigation into molecular mechanisms underlying male reproduction.</p

PARTS: Probabilistic Alignment for RNA joinT Secondary structure prediction

A novel method is presented for joint prediction of alignment and common secondary structures of two RNA sequences. The joint consideration of common secondary structures and alignment is accomplished by structural alignment over a search space defined by the newly introduced motif called matched helical regions. The matched helical region formulation generalizes previously employed constraints for structural alignment and thereby better accommodates the structural variability within RNA families. A probabilistic model based on pseudo free energies obtained from precomputed base pairing and alignment probabilities is utilized for scoring structural alignments. Maximum a posteriori (MAP) common secondary structures, sequence alignment and joint posterior probabilities of base pairing are obtained from the model via a dynamic programming algorithm called PARTS. The advantage of the more general structural alignment of PARTS is seen in secondary structure predictions for the RNase P family. For this family, the PARTS MAP predictions of secondary structures and alignment perform significantly better than prior methods that utilize a more restrictive structural alignment model. For the tRNA and 5S rRNA families, the richer structural alignment model of PARTS does not offer a benefit and the method therefore performs comparably with existing alternatives. For all RNA families studied, the posterior probability estimates obtained from PARTS offer an improvement over posterior probability estimates from a single sequence prediction. When considering the base pairings predicted over a threshold value of confidence, the combination of sensitivity and positive predictive value is superior for PARTS than for the single sequence prediction. PARTS source code is available for download under the GNU public license at http://rna.urmc.rochester.edu

Towards Reliable Automatic Protein Structure Alignment

A variety of methods have been proposed for structure similarity calculation, which are called structure alignment or superposition. One major shortcoming in current structure alignment algorithms is in their inherent design, which is based on local structure similarity. In this work, we propose a method to incorporate global information in obtaining optimal alignments and superpositions. Our method, when applied to optimizing the TM-score and the GDT score, produces significantly better results than current state-of-the-art protein structure alignment tools. Specifically, if the highest TM-score found by TMalign is lower than (0.6) and the highest TM-score found by one of the tested methods is higher than (0.5), there is a probability of (42%) that TMalign failed to find TM-scores higher than (0.5), while the same probability is reduced to (2%) if our method is used. This could significantly improve the accuracy of fold detection if the cutoff TM-score of (0.5) is used. In addition, existing structure alignment algorithms focus on structure similarity alone and simply ignore other important similarities, such as sequence similarity. Our approach has the capacity to incorporate multiple similarities into the scoring function. Results show that sequence similarity aids in finding high quality protein structure alignments that are more consistent with eye-examined alignments in HOMSTRAD. Even when structure similarity itself fails to find alignments with any consistency with eye-examined alignments, our method remains capable of finding alignments highly similar to, or even identical to, eye-examined alignments.Comment: Peer-reviewed and presented as part of the 13th Workshop on Algorithms in Bioinformatics (WABI2013

RNAG: a new Gibbs sampler for predicting RNA secondary structure for unaligned sequences

Motivation: RNA secondary structure plays an important role in the function of many RNAs, and structural features are often key to their interaction with other cellular components. Thus, there has been considerable interest in the prediction of secondary structures for RNA families. In this article, we present a new global structural alignment algorithm, RNAG, to predict consensus secondary structures for unaligned sequences. It uses a blocked Gibbs sampling algorithm, which has a theoretical advantage in convergence time. This algorithm iteratively samples from the conditional probability distributions P(Structure | Alignment) and P(Alignment | Structure). Not surprisingly, there is considerable uncertainly in the high-dimensional space of this difficult problem, which has so far received limited attention in this field. We show how the samples drawn from this algorithm can be used to more fully characterize the posterior space and to assess the uncertainty of predictions

Rfam: Wikipedia, clans and the “decimal” release

The Rfam database aims to catalogue non-coding RNAs through the use of sequence alignments and statistical profile models known as covariance models. In this contribution, we discuss the pros and cons of using the online encyclopedia, Wikipedia, as a source of community‐derived annotation. We discuss the addition of groupings of related RNA families into clans and new developments to the website. Rfam is available on the Web at http://rfam.sanger.ac.uk

Mechanisms of Epithelial-Mesenchymal Transition of Peritoneal Mesothelial Cells During Peritoneal Dialysis

A growing body of evidence indicates that epithelial-mesenchymal transition (EMT) of human peritoneal mesothelial cells (HPMC) may play an important role in the development and progression of peritoneal fibrosis during long-term peritoneal dialysis (PD) leading to failure of peritoneal membrane function. Here, we review our own observations and those of others on the mechanisms of EMT of HPMC and suggest potential therapeutic strategies to prevent EMT and peritoneal fibrosis during long-term PD. We found that high glucose and H2O2 as well as transforming growth factor-β1 (TGF-β1) induced EMT in HPMC and that high glucose-induced EMT was blocked not only by inhibition of TGF-β1 but also by antioxidants or inhibitors of mitogen-activated protein kinases (MAPK). Since MAPKs are downstream target molecules of reactive oxygen species (ROS), these data suggest that high glucose-induced generation of ROS and subsequent MAPK activation mediate high glucose-induced EMT in HPMC. We and others also observed that bone morphogenetic protein-7 (BMP-7) prevented EMT in HPMC. Glucose degradation products (GDP) were shown to play a role in inducing EMT. Involvement of a mammalian target of rapamycin (mTOR) in TGF-β1-induced EMT has also been proposed in cultured HPMC. A better understanding of the precise mechanisms involved in EMT of HPMC may provide new therapeutic strategies for inhibiting peritoneal fibrosis in long-term PD patients

Evolutionary distances in the twilight zone -- a rational kernel approach

Phylogenetic tree reconstruction is traditionally based on multiple sequence alignments (MSAs) and heavily depends on the validity of this information bottleneck. With increasing sequence divergence, the quality of MSAs decays quickly. Alignment-free methods, on the other hand, are based on abstract string comparisons and avoid potential alignment problems. However, in general they are not biologically motivated and ignore our knowledge about the evolution of sequences. Thus, it is still a major open question how to define an evolutionary distance metric between divergent sequences that makes use of indel information and known substitution models without the need for a multiple alignment. Here we propose a new evolutionary distance metric to close this gap. It uses finite-state transducers to create a biologically motivated similarity score which models substitutions and indels, and does not depend on a multiple sequence alignment. The sequence similarity score is defined in analogy to pairwise alignments and additionally has the positive semi-definite property. We describe its derivation and show in simulation studies and real-world examples that it is more accurate in reconstructing phylogenies than competing methods. The result is a new and accurate way of determining evolutionary distances in and beyond the twilight zone of sequence alignments that is suitable for large datasets.Comment: to appear in PLoS ON

First report of generalized face processing difficulties in möbius sequence.

Reverse simulation models of facial expression recognition suggest that we recognize the emotions of others by running implicit motor programmes responsible for the production of that expression. Previous work has tested this theory by examining facial expression recognition in participants with Möbius sequence, a condition characterized by congenital bilateral facial paralysis. However, a mixed pattern of findings has emerged, and it has not yet been tested whether these individuals can imagine facial expressions, a process also hypothesized to be underpinned by proprioceptive feedback from the face. We investigated this issue by examining expression recognition and imagery in six participants with Möbius sequence, and also carried out tests assessing facial identity and object recognition, as well as basic visual processing. While five of the six participants presented with expression recognition impairments, only one was impaired at the imagery of facial expressions. Further, five participants presented with other difficulties in the recognition of facial identity or objects, or in lower-level visual processing. We discuss the implications of our findings for the reverse simulation model, and suggest that facial identity recognition impairments may be more severe in the condition than has previously been noted