Base-pair ambiguity and the kinetics of RNA folding.

BACKGROUND: A pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between "bound" molecules, which are active as part of ribonucleoprotein (RNP) complexes, and "unbound," with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an "ambiguity index"-a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative ("gold standard") and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules. RESULTS: We have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule-an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure. CONCLUSIONS: A statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.

Mapping the RNA-Seq trash bin: unusual transcripts in prokaryotic transcriptome sequencing data.

Prokaryotic transcripts constitute almost always uninterrupted intervals when mapped back to the genome. Split reads, i.e., RNA-seq reads consisting of parts that only map to discontiguous loci, are thus disregarded in most analysis pipelines. There are, however, some well-known exceptions, in particular, tRNA splicing and circularized small RNAs in Archaea as well as self-splicing introns. Here, we reanalyze a series of published RNA-seq data sets, screening them specifically for non-contiguously mapping reads. We recover most of the known cases together with several novel archaeal ncRNAs associated with circularized products. In Eubacteria, only a handful of interesting candidates were obtained beyond a few previously described group I and group II introns. Most of the atypically mapping reads do not appear to correspond to well-defined, specifically processed products. Whether this diffuse background is, at least in part, an incidental by-product of prokaryotic RNA processing or whether it consists entirely of technical artifacts of reverse transcription or amplification remains unknown.

Predicting Human miRNA-like Sequences within Human Papillomavirus Genomes.

BACKGROUND: This study presents a prediction of putative miRNA within several Human Papillomavirus (HPV) types by using bioinformatics tools and a strategy based on sequence and structure alignment. Currently, little is known about HPV miRNAs. METHODS: Computational methods have been widely applied in the identification of novel miRNAs when analyzing genome sequences. Here, ten whole-genome sequences from HPV-6, -11, -16, -18, -31, -33, -35, -45, -52, and -58 were analyzed. Software based on local contiguous structure-sequence features and support vector machine (SVM), as well as additional bioinformatics tools, were utilized for identification and classification of real and pseudo microRNA precursors. RESULTS: An initial analysis predicted 200 putative pre-miRNAs for all the ten HPV genome variants. To derive a smaller set of pre-miRNAs candidates, stringent validation criteria was conducted by applying <‒10 ΔG value (Gibbs Free Energy). Thus, only pre-miRNAs with total scores above the cut-off points of 90% were considered as putative pre-miRNAs. As a result of this strategy, 19 pre-miRNAs were selected (hpv-pre-miRNAs). These novel pre-miRNAs were located in different clusters within HPV genomes and some of them were positioned at splice regions. Additionally, the 19 identified pre-miRNAs sequences varied between HPV genotypes. Interestingly, the newly identified miRNAs, 297, 27b, 500, 501-5, and 509-3-5p, were closely implicated in carcinogenesis participating in cellular longevity, cell cycle, metastasis, apoptosis evasion, tissue invasion and cellular growth pathways. CONCLUSIONS: The novel putative miRNAs candidates could be promising biomarkers of HPV infection and furthermore, could be targeted for potential therapeutic interventions in HPV-induced malignancies.

Mobile Introns Shape the Genetic Diversity of Their Host Genes.

Self-splicing introns populate several highly conserved protein-coding genes in fungal and plant mitochondria. In fungi, many of these introns have retained their ability to spread to intron-free target sites, often assisted by intron-encoded endonucleases that initiate the homing process. Here, leveraging population genomic data from Saccharomyces cerevisiae, Schizosaccharomyces pombe, and Lachancea kluyveri, we expose nonrandom patterns of genetic diversity in exons that border self-splicing introns. In particular, we show that, in all three species, the density of single nucleotide polymorphisms increases as one approaches a mobile intron. Through multiple lines of evidence, we rule out relaxed purifying selection as the cause of uneven nucleotide diversity. Instead, our findings implicate intron mobility as a direct driver of host gene diversity. We discuss two mechanistic scenarios that are consistent with the data: either endonuclease activity and subsequent error-prone repair have left a mutational footprint on the insertion environment of mobile introns or nonrandom patterns of genetic diversity are caused by exonic coconversion, which occurs when introns spread to empty target sites via homologous recombination. Importantly, however, we show that exonic coconversion can only explain diversity gradients near intron-exon boundaries if the conversion template comes from outside the population. In other words, there must be pervasive and ongoing horizontal gene transfer of self-splicing introns into extant fungal populations.