LocARNA-P: accurate boundary prediction and improved detection of structural RNAs

Genome-wide identification of Wig-1 mRNA targets by RIP-Seq analysis

RNA-binding proteins (RBPs) play important roles in the regulation of gene expression through a variety of post-transcriptional mechanisms. The p53-induced RBP Wig-1 (Zmat3) binds RNA through its zinc finger domains and enhances stability of p53 and N-Myc mRNAs and decreases stability of FAS mRNA. To identify novel Wig-1-bound RNAs, we performed RNA-immunoprecipitation followed by high-throughput sequencing (RIP-Seq) in HCT116 and Saos-2 cells. We identified 286 Wig-1-bound mRNAs common between the two cell lines. Sequence analysis revealed that AU-rich elements (AREs) are highly enriched in the 3′UTR of these Wig-1-bound mRNAs. Network enrichment analysis showed that Wig-1 preferentially binds mRNAs involved in cell cycle regulation. Moreover, we identified a 2D Wig-1 binding motif in HIF1A mRNA. Our findings confirm that Wig-1 is an ARE-BP that regulates cell cycle-related processes and provide a novel view of how Wig-1 may bind mRNA through a putative structural motif. We also significantly extend the repertoire of Wig-1 target mRNAs. Since Wig-1 is a transcriptional target of the tumor suppressor p53, these results have implications for our understanding of p53-dependent stress responses and tumor suppression.

Phylogenetic distribution of plant snoRNA families

Background

Small nucleolar RNAs (snoRNAs) are one of the most ancient families amongst non-protein-coding RNAs. They are ubiquitous in Archaea and Eukarya but absent in bacteria. Their main function is to target chemical modifications of ribosomal RNAs. They fall into two classes, box C/D snoRNAs and box H/ACA snoRNAs, which are clearly distinguished by conserved sequence motifs and the type of chemical modification that they govern. Similarly to microRNAs, snoRNAs appear in distinct families of homologs that affect homologous targets. In animals, snoRNAs and their evolution have been studied in much detail. In plants, however, their evolution has attracted comparably little attention.

Results

In order to chart the phylogenetic distribution of individual snoRNA families in plants, we applied a sophisticated approach for identifying homologs of known plant snoRNAs across the plant kingdom. In response to the relatively fast evolution of snoRNAs, information on conserved sequence boxes, target sequences, and secondary structure is combined to identify additional snoRNAs. We identified 296 families of snoRNAs in 24 species and traced their evolution throughout the plant kingdom. Many of the plant snoRNA families comprise paralogs. We also found that targets are well-conserved for most snoRNA families.

Conclusions

The sequence conservation of snoRNAs is sufficient to establish homologies between phyla. The degree of this conservation tapers off, however, between land plants and algae. Plant snoRNAs are frequently organized in highly conserved spatial clusters. As a resource for further investigations we provide carefully curated and annotated alignments for each snoRNA family under investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3301-2) contains supplementary material, which is available to authorized users.

Complete mitochondrial genome and evolutionary analysis of Turritopsis dohrnii, the "immortal" jellyfish with a reversible life-cycle

Turritopsis dohrnii (Cnidaria, Hydrozoa, Hydroidolina, Anthoathecata) is the only known metazoan that is capable of reversing its life cycle via morph rejuvenation from the adult medusa stage to the juvenile polyp stage. Here, we present a complete mitochondrial (mt) genome sequence of T. dohrnii, which harbors genes for 13 proteins, two transfer RNAs, and two ribosomal RNAs. The T. dohrnii mt genome is characterized by typical features of species in the Hydroidolina subclass, such as a high A+T content (71.5%), reversed transcriptional orientation for the large rRNA subunit gene, and paucity of CGN codons. An incomplete complementary duplicate of the cox1 gene was found at the 5' end of the T. dohrnii mt chromosome, as were variable repeat regions flanking the chromosome. We identified species-specific variations (nad5, nad6, cob, and cox1 genes) and putative selective constraints (atp8, nad1, nad2, and nad5 genes) in the mt genes of T. dohrnii, and predicted alterations in tertiary structures of respiratory chain proteins (NADH4, NADH5, and COX1 proteins) of T. dohrnii. Based on comparative analyses of available hydrozoan mt genomes, we also determined the taxonomic relationships of T. dohrnii, recovering Filifera IV as a paraphyletic taxon, and assessed intraspecific diversity of various Hydrozoa species.

The Conserved Dcw Gene Cluster of R. sphaeroides Is Preceded by an Uncommonly Extended 5' Leader Featuring the sRNA UpsM

Cell division and cell wall synthesis mechanisms are similarly conserved among bacteria. Consequently some bacterial species have comparable sets of genes organized in the dcw (division andcellwall) gene cluster. Dcw genes, their regulation and their relative order within the cluster are outstandingly conserved among rod shaped and gram negative bacteria to ensure an efficient coordination of growth and division. A well studied representative is the dcw gene cluster of E. coli. The first promoter of the gene cluster (mraZ1p) gives rise to polycistronic transcripts containing a 38 nt long 5’ UTR followed by the first gene mraZ. Despite reported conservation we present evidence for a much longer 5’ UTR in the gram negative and rod shaped bacterium Rhodobacter sphaeroides and in the family of Rhodobacteraceae. This extended 268 nt long 5’ UTR comprises a Rho independent terminator, which in case of termination gives rise to a non-coding RNA (UpsM). This sRNA is conditionally cleaved by RNase E under stress conditions in an Hfq- and very likely target mRNA-dependent manner, implying its function in trans. These results raise the question for the regulatory function of this extended 5’ UTR. It might represent the rarely described case of a trans acting sRNA derived from a riboswitch with exclusive presence in the family of Rhodobacteraceae.

Plasticity of archaeal C/D box sRNA biogenesis

Archaeal and eukaryotic organisms contain sets of C/D box s(no)RNAs with guide sequences that determine ribose 2'-O-methylation sites of target RNAs. The composition of these C/D box sRNA sets is highly variable between organisms and results in varying RNA modification patterns which are important for ribosomal RNA folding and stability. Little is known about the genomic organization of C/D box sRNA genes in archaea. Here, we aimed to obtain first insights into the biogenesis of these archaeal C/D box sRNAs and analyzed the genetic context of more than 300 archaeal sRNA genes. We found that the majority of these genes do not possess independent promoters but are rather located at positions that allow for co-transcription with neighboring genes and their start or stop codons were frequently incorporated into the conserved boxC and D motifs. The biogenesis of plasmid-encoded C/D box sRNA variants was analyzed in vivo in Sulfolobus acidocaldarius. It was found that C/D box sRNA maturation occurs independent of their genetic context and relies solely on the presence of intact RNA kink-turn structures. The observed plasticity of C/D box sRNA biogenesis is suggested to enable their accelerated evolution and, consequently, allow for adjustments of the RNA modification landscape.

Mycoplasma non-coding RNA: identification of small RNAs and targets

Background

Bacterial non-coding RNAs act by base-pairing as regulatory elements in crucial biological processes. We performed the identification of trans-encoded small RNAs (sRNA) from the genomes of Mycoplama hyopneumoniae, Mycoplasma flocculare and Mycoplasma hyorhinis, which are Mycoplasma species that have been identified in the porcine respiratory system.

Results

A total of 47, 15 and 11 putative sRNAs were predicted in M. hyopneumoniae, M. flocculare and M. hyorhinis, respectively. A comparative genomic analysis revealed the presence of species or lineage specific sRNA candidates. Furthermore, the expression profile of some M. hyopneumoniae sRNAs was determined by a reverse transcription amplification approach, in three different culture conditions. All tested sRNAs were transcribed in at least one condition. A detailed investigation revealed a differential expression profile for two M. hyopneumoniae sRNAs in response to oxidative and heat shock stress conditions, suggesting that their expression is influenced by environmental signals. Moreover, we analyzed sRNA-mRNA hybrids and accessed putative target genes for the novel sRNA candidates. The majority of the sRNAs showed interaction with multiple target genes, some of which could be linked to pathogenesis and cell homeostasis activity.

Conclusion

This study contributes to our knowledge of Mycoplasma sRNAs and their response to environmental changes. Furthermore, the mRNA target prediction provides a perspective for the characterization and comprehension of the function of the sRNA regulatory mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3061-z) contains supplementary material, which is available to authorized users.

Elaborate uORF/IRES features control expression and localization of human glycyl-tRNA synthetase

The canonical activity of glycyl-tRNA synthetase (GARS) is to charge glycine onto its cognate tRNAs. However, outside translation, GARS also participates in many other functions. A single gene encodes both the cytosolic and mitochondrial forms of GARS but 2 mRNA isoforms were identified. Using immunolocalization assays, in vitro translation assays and bicistronic constructs we provide experimental evidence that one of these mRNAs tightly controls expression and localization of human GARS. An intricate regulatory domain was found in its 5'-UTR which displays a functional Internal Ribosome Entry Site and an upstream Open Reading Frame. Together, these elements hinder the synthesis of the mitochondrial GARS and target the translation of the cytosolic enzyme to ER-bound ribosomes. This finding reveals a complex picture of GARS translation and localization in mammals. In this context, we discuss how human GARS expression could influence its moonlighting activities and its involvement in diseases.

Coronavirus cis-Acting RNA Elements

Coronaviruses have exceptionally large RNA genomes of approximately 30 kilobases. Genome replication and transcription is mediated by a multisubunit protein complex comprised of more than a dozen virus-encoded proteins. The protein complex is thought to bind specific cis-acting RNA elements primarily located in the 5'- and 3'-terminal genome regions and upstream of the open reading frames located in the 3'-proximal one-third of the genome. Here, we review our current understanding of coronavirus cis-acting RNA elements, focusing on elements required for genome replication and packaging. Recent bioinformatic, biochemical, and genetic studies suggest a previously unknown level of conservation of cis-acting RNA structures among different coronavirus genera and, in some cases, even beyond genus boundaries. Also, there is increasing evidence to suggest that individual cis-acting elements may be part of higher-order RNA structures involving long-range and dynamic RNA-RNA interactions between RNA structural elements separated by thousands of nucleotides in the viral genome. We discuss the structural and functional features of these cis-acting RNA elements and their specific functions in coronavirus RNA synthesis.

High diversity of picornaviruses in rats from different continents revealed by deep sequencing

Outbreaks of zoonotic diseases in humans and livestock are not uncommon, and an important component in containment of such emerging viral diseases is rapid and reliable diagnostics. Such methods are often PCR-based and hence require the availability of sequence data from the pathogen. Rattus norvegicus (R. norvegicus) is a known reservoir for important zoonotic pathogens. Transmission may be direct via contact with the animal, for example, through exposure to its faecal matter, or indirectly mediated by arthropod vectors. Here we investigated the viral content in rat faecal matter (n=29) collected from two continents by analyzing 2.2 billion next-generation sequencing reads derived from both DNA and RNA. Among other virus families, we found sequences from members of the Picornaviridae to be abundant in the microbiome of all the samples. Here we describe the diversity of the picornavirus-like contigs including near-full-length genomes closely related to the Boone cardiovirus and Theiler's encephalomyelitis virus. From this study, we conclude that picornaviruses within R. norvegicus are more diverse than previously recognized. The virome of R. norvegicus should be investigated further to assess the full potential for zoonotic virus transmission.

Fly Models of Human Diseases: Drosophila as a Model for Understanding Human Mitochondrial Mutations and Disease

Mitochondrial diseases are a prevalent, heterogeneous class of diseases caused by defects in oxidative phosphorylation, whose severity depends upon particular genetic mutations. These diseases can be difficult to diagnose, and current therapeutics have limited efficacy, primarily treating only symptoms. Because mitochondria play a pivotal role in numerous cellular functions, especially ATP production, their diminished activity has dramatic physiological consequences. While this in and of itself makes treating mitochondrial disease complex, these organelles contain their own DNA, mtDNA, whose products are required for ATP production, in addition to the hundreds of nucleus-encoded proteins. Drosophila offers a tractable whole-animal model to understand the mechanisms underlying loss of mitochondrial function, the subsequent cellular and tissue damage that results, and how these organelles are inherited. Human and Drosophila mtDNAs encode the same set of products, and the homologous nucleus-encoded genes required for mitochondrial function are conserved. In addition, Drosophila contain sufficiently complex organ systems to effectively recapitulate many basic symptoms of mitochondrial diseases, yet are relatively easy and fast to genetically manipulate. There are several Drosophila models for specific mitochondrial diseases, which have been recently reviewed (Foriel, Willems, Smeitink, Schenck, & Beyrath, 2015). In this review, we highlight the conservation between human and Drosophila mtDNA, the present and future techniques for creating mtDNA mutations for further study, and how Drosophila has contributed to our current understanding of mitochondrial inheritance.

PR2S2Clust: Patched RNA-seq read segments' structure-oriented clustering

RNA-seq, the next generation sequencing platform, enables researchers to explore deep into the transcriptome of organisms, such as identifying functional non-coding RNAs (ncRNAs), and quantify their expressions on tissues. The functions of ncRNAs are mostly related to their secondary structures. Thus by exploring the clustering in terms of structural profiles of the corresponding read-segments would be essential and this fuels in our motivation behind this research. In this manuscript we proposed PR2S2Clust, Patched RNA-seq Read Segments' Structure-oriented Clustering, which is an analysis platform to extract features to prepare the secondary structure profiles of the RNA-seq read segments. It provides a strategy to employ the profiles to annotate the segments into ncRNA classes using several clustering strategies. The system considers seven pairwise structural distance metrics by considering short-read mappings onto each structure, which we term as the "patched structure" while clustering the segments. In this regard, we show applications of both classical and ensemble clusterings of the partitional and hierarchical variations. Extensive real-world experiments over three publicly available RNA-seq datasets and a comparative analysis over four competitive systems confirm the effectiveness and superiority of the proposed system. The source codes and dataset of PR2S2Clust are available at the http://biomecis.uta.edu/~ashis/res/PR2S2Clust-suppl/ .

Evolutionary clues in lncRNAs

The diversity of long non-coding RNAs (lncRNAs) in the human transcriptome is in stark contrast to the sparse exploration of their functions concomitant with their conservation and evolution. The pervasive transcription of the largely non-coding human genome makes the evolutionary age and conservation patterns of lncRNAs to a topic of interest. Yet it is a fairly unexplored field and not that easy to determine as for protein-coding genes. Although there are a few experimentally studied cases, which are conserved at the sequence level, most lncRNAs exhibit weak or untraceable primary sequence conservation. Recent studies shed light on the interspecies conservation of secondary structures among lncRNA homologs by using diverse computational methods. This highlights the importance of structure on functionality of lncRNAs as opposed to the poor impact of primary sequence changes. Further clues in the evolution of lncRNAs are given by selective constraints on non-coding gene structures (e.g., promoters or splice sites) as well as the conservation of prevalent spatio-temporal expression patterns. However, a rapid evolutionary turnover is observable throughout the heterogeneous group of lncRNAs. This still gives rise to questions about its functional meaning. WIREs RNA 2017, 8:e1376. doi: 10.1002/wrna.1376 For further resources related to this article, please visit the WIREs website.

IF-combined smRNA FISH reveals interaction of MCPIP1 protein with IER3 mRNA

MCPIP1 and IER3 are recently described proteins essential for maintenance of immune homeostasis. IER3 is involved in the regulation of apoptosis and differentiation and has been shown lately to protect activated T cells and macrophages from apoptosis. MCPIP1 is an RNase critical for controlling inflammation-related mRNAs. MCPIP1 interacts with and degrades a set of stem-loop-containing mRNAs (including IL-6). Our results demonstrate the involvement of MCPIP1 in the regulation of IER3 mRNA levels. A dual luciferase assay revealed that over-expression of MCPIP1 resulted in a decrease of luciferase activity in the samples co-transfected with constructs containing luciferase CDS attached to IER3 3′UTR. We identified a stem-loop structure similar to that described to be important for destabilization of the IL-6 mRNA by MCPIP1. Examination of IER3 3′UTR sequence, structure and evolutionary conservation revealed that the identified stem-loop is buried within a bigger element. Deletion of this fragment abolished the regulation of IER3 3′UTR-containing transcript by MCPIP1. Finally, using immunofluorescence-combined single-molecule RNA FISH we have shown that the MCPIP1 protein co-localizes with IER3 mRNA. By this method we also proved that the presence of the wild-type NYN/PIN-like domain of MCPIP1 correlated with the decreased level of IER3 mRNA. RNA immunoprecipitation further confirmed the interaction of MCPIP1 with IER3 transcripts in vivo.

Summary: We identify IER3 mRNA as a newly discovered MCPIP1 target using recently described IF-based procedures, and also identify a conserved element involved in MCPIP1-dependent IER3 transcript destabilization.

Argonaute-associated short introns are a novel class of gene regulators

MicroRNAs (miRNAs) are short (∼22 nucleotides) regulators of gene expression acting by direct base pairing to 3'-UTR target sites in messenger RNAs. Mature miRNAs are produced by two sequential endonucleolytic cleavages facilitated by Drosha in the nucleus and Dicer in the cytoplasm. A subclass of miRNAs, termed mirtrons, derives from short introns and enters the miRNA biogenesis pathway as Dicer substrates. Here we uncover a third biogenesis strategy that, similar to mirtron biogenesis, initiates from short introns but bypasses Dicer cleavage. These short introns (80-100 nucleotides), coined agotrons, are associated with and stabilized by Argonaute (Ago) proteins in the cytoplasm. Some agotrons are completely conserved in mammalian species, suggesting that they are functionally important. Furthermore, we demonstrate that the agotrons are capable of repressing mRNAs with seed-matching target sequences in the 3'-UTR. These data provide evidence for a novel RNA regulator of gene expression, which bypasses the canonical miRNA biogenesis machinery.

Predicting RNA secondary structures from sequence and probing data

RNA secondary structures have proven essential for understanding the regulatory functions performed by RNA such as microRNAs, bacterial small RNAs, or riboswitches. This success is in part due to the availability of efficient computational methods for predicting RNA secondary structures. Recent advances focus on dealing with the inherent uncertainty of prediction by considering the ensemble of possible structures rather than the single most stable one. Moreover, the advent of high-throughput structural probing has spurred the development of computational methods that incorporate such experimental data as auxiliary information.

Novel modes of RNA editing in mitochondria

Gene structure and expression in diplonemid mitochondria are unparalleled. Genes are fragmented in pieces (modules) that are separately transcribed, followed by the joining of module transcripts to contiguous RNAs. Some instances of unique uridine insertion RNA editing at module boundaries were noted, but the extent and potential occurrence of other editing types remained unknown. Comparative analysis of deep transcriptome and genome data from Diplonema papillatum mitochondria reveals ∼220 post-transcriptional insertions of uridines, but no insertions of other nucleotides nor deletions. In addition, we detect in total 114 substitutions of cytosine by uridine and adenosine by inosine, amassed into unusually compact clusters. Inosines in transcripts were confirmed experimentally. This is the first report of adenosine-to-inosine editing of mRNAs and ribosomal RNAs in mitochondria. In mRNAs, editing causes mostly amino-acid additions and non-synonymous substitutions; in ribosomal RNAs, it permits formation of canonical secondary structures. Two extensively edited transcripts were compared across four diplonemids. The pattern of uridine-insertion editing is strictly conserved, whereas substitution editing has diverged dramatically, but still rendering diplonemid proteins more similar to other eukaryotic orthologs. We posit that RNA editing not only compensates but also sustains, or even accelerates, ultra-rapid evolution of genome structure and sequence in diplonemid mitochondria.

Contrasting evolutionary patterns of 28S and ITS rRNA genes reveal high intragenomic variation in Cephalenchus (Nematoda): Implications for species delimitation

Concerted evolution is often assumed to be the evolutionary force driving multi-family genes, including those from ribosomal DNA (rDNA) repeat, to complete homogenization within a species, although cases of non-concerted evolution have been also documented. In this study, sequence variation of 28S and ITS ribosomal RNA (rRNA) genes in the genus Cephalenchus is assessed at three different levels, intragenomic, intraspecific, and interspecific. The findings suggest that not all Cephalenchus species undergo concerted evolution. High levels of intraspecific polymorphism, mostly due to intragenomic variation, are found in Cephalenchus sp1 (BRA-01). Secondary structure analyses of both rRNA genes and across different species show a similar substitution pattern, including mostly compensatory (CBC) and semi-compensatory (SBC) base changes, thus suggesting the functionality of these rRNA copies despite the variation found in some species. This view is also supported by low sequence variation in the 5.8S gene in relation to the flanking ITS-1 and ITS-2 as well as by the existence of conserved motifs in the former gene. It is suggested that potential cross-fertilization in some Cephalenchus species, based on inspection of female reproductive system, might contribute to both intragenomic and intraspecific polymorphism of their rRNA genes. These results reinforce the potential implications of intragenomic and intraspecific genetic diversity on species delimitation, especially in biodiversity studies based solely on metagenetic approaches. Knowledge of sequence variation will be crucial for accurate species diversity estimation using molecular methods.

An assessment of bacterial small RNA target prediction programs

Most bacterial regulatory RNAs exert their function through base-pairing with target RNAs. Computational prediction of targets is a busy research field that offers biologists a variety of web sites and software. However, it is difficult for a non-expert to evaluate how reliable those programs are. Here, we provide a simple benchmark for bacterial sRNA target prediction based on trusted E. coli sRNA/target pairs. We use this benchmark to assess the most recent RNA target predictors as well as earlier programs for RNA-RNA hybrid prediction. Moreover, we consider how the definition of mRNA boundaries can impact overall predictions. Recent algorithms that exploit both conservation of targets and accessibility information offer improved accuracy over previous software. However, even with the best predictors, the number of true biological targets with low scores and non-targets with high scores remains puzzling.

Identification of Conserved and Potentially Regulatory Small RNAs in Heterocystous Cyanobacteria

Small RNAs (sRNAs) are a growing class of non-protein-coding transcripts that participate in the regulation of virtually every aspect of bacterial physiology. Heterocystous cyanobacteria are a group of photosynthetic organisms that exhibit multicellular behavior and developmental alternatives involving specific transcriptomes exclusive of a given physiological condition or even a cell type. In the context of our ongoing effort to understand developmental decisions in these organisms we have undertaken an approach to the global identification of sRNAs. Using differential RNA-Seq we have previously identified transcriptional start sites for the model heterocystous cyanobacterium Nostoc sp. PCC 7120. Here we combine this dataset with a prediction of Rho-independent transcriptional terminators and an analysis of phylogenetic conservation of potential sRNAs among 89 available cyanobacterial genomes. In contrast to predictive genome-wide approaches, the use of an experimental dataset comprising all active transcriptional start sites (differential RNA-Seq) facilitates the identification of bona fide sRNAs. The output of our approach is a dataset of predicted potential sRNAs in Nostoc sp. PCC 7120, with different degrees of phylogenetic conservation across the 89 cyanobacterial genomes analyzed. Previously described sRNAs appear among the predicted sRNAs, demonstrating the performance of the algorithm. In addition, new predicted sRNAs are now identified that can be involved in regulation of different aspects of cyanobacterial physiology, including adaptation to nitrogen stress, the condition that triggers differentiation of heterocysts (specialized nitrogen-fixing cells). Transcription of several predicted sRNAs that appear exclusively in the genomes of heterocystous cyanobacteria is experimentally verified by Northern blot. Cell-specific transcription of one of these sRNAs, NsiR8 (nitrogen stress-induced RNA 8), in developing heterocysts is also demonstrated.

Loss of Conserved Noncoding RNAs in Genomes of Bacterial Endosymbionts

The genomes of intracellular symbiotic or pathogenic bacteria, such as of Buchnera, Mycoplasma, and Rickettsia, are typically smaller compared with their free-living counterparts. Here we showed that noncoding RNA (ncRNA) families, which are conserved in free-living bacteria, frequently could not be detected by computational methods in the small genomes. Statistical tests demonstrated that their absence is not an artifact of low GC content or small deletions in these small genomes, and thus it was indicative of an independent loss of ncRNAs in different endosymbiotic lineages. By analyzing the synteny (conservation of gene order) between the reduced and nonreduced genomes, we revealed instances of protein-coding genes that were preserved in the reduced genomes but lost cis-regulatory elements. We found that the loss of cis-regulatory ncRNA sequences, which regulate the expression of cognate protein-coding genes, is characterized by the reduction of secondary structure formation propensity, GC content, and length of the corresponding genomic regions.

Multiple sequence alignment modeling: methods and applications

This review provides an overview on the development of Multiple sequence alignment (MSA) methods and their main applications. It is focused on progress made over the past decade. The three first sections review recent algorithmic developments for protein, RNA/DNA and genomic alignments. The fourth section deals with benchmarks and explores the relationship between empirical and simulated data, along with the impact on method developments. The last part of the review gives an overview on available MSA local reliability estimators and their dependence on various algorithmic properties of available methods.

Survey of Natural Language Processing Techniques in Bioinformatics

Informatics methods, such as text mining and natural language processing, are always involved in bioinformatics research. In this study, we discuss text mining and natural language processing methods in bioinformatics from two perspectives. First, we aim to search for knowledge on biology, retrieve references using text mining methods, and reconstruct databases. For example, protein-protein interactions and gene-disease relationship can be mined from PubMed. Then, we analyze the applications of text mining and natural language processing techniques in bioinformatics, including predicting protein structure and function, detecting noncoding RNA. Finally, numerous methods and applications, as well as their contributions to bioinformatics, are discussed for future use by text mining and natural language processing researchers.

Optimizing RNA structures by sequence extensions using RNAcop

A key aspect of RNA secondary structure prediction is the identification of novel functional elements. This is a challenging task because these elements typically are embedded in longer transcripts where the borders between the element and flanking regions have to be defined. The flanking sequences impact the folding of the functional elements both at the level of computational analyses and when the element is extracted as a transcript for experimental analysis. Here, we analyze how different flanking region lengths impact folding into a constrained structure by computing probabilities of folding for different sizes of flanking regions. Our method, RNAcop (RNA context optimization by probability), is tested on known and de novo predicted structures. In vitro experiments support the computational analysis and suggest that for a number of structures, choosing proper lengths of flanking regions is critical. RNAcop is available as web server and stand-alone software via http://rth.dk/resources/rnacop.

RC3H1 post-transcriptionally regulates A20 mRNA and modulates the activity of the IKK/NF-κB pathway

The RNA-binding protein RC3H1 (also known as ROQUIN) promotes TNFα mRNA decay via a 3′UTR constitutive decay element (CDE). Here we applied PAR-CLIP to human RC3H1 to identify ∼3,800 mRNA targets with >16,000 binding sites. A large number of sites are distinct from the consensus CDE and revealed a structure-sequence motif with U-rich sequences embedded in hairpins. RC3H1 binds preferentially short-lived and DNA damage-induced mRNAs, indicating a role of this RNA-binding protein in the post-transcriptional regulation of the DNA damage response. Intriguingly, RC3H1 affects expression of the NF-κB pathway regulators such as IκBα and A20. RC3H1 uses ROQ and Zn-finger domains to contact a binding site in the A20 3′UTR, demonstrating a not yet recognized mode of RC3H1 binding. Knockdown of RC3H1 resulted in increased A20 protein expression, thereby interfering with IκB kinase and NF-κB activities, demonstrating that RC3H1 can modulate the activity of the IKK/NF-κB pathway.

The RNA-binding protein RC3H1/ROQUIN1 promotes the degradation of mRNA by binding to a consensus CDE present in the 3′UTR. Here the authors expand the set of consensus sequences through which RCH31 binds and regulates mRNA encoding members of the DNA damage response and IKK/NF-κB pathway.

Conserved RNA secondary structures and long-range interactions in hepatitis C viruses

Hepatitis C virus (HCV) is a hepatotropic virus with a plus-strand RNA genome of ∼9.600 nt. Due to error-prone replication by its RNA-dependent RNA polymerase (RdRp) residing in nonstructural protein 5B (NS5B), HCV isolates are grouped into seven genotypes with several subtypes. By using whole-genome sequences of 106 HCV isolates and secondary structure alignments of the plus-strand genome and its minus-strand replication intermediate, we established refined secondary structures of the 5' untranslated region (UTR), the cis-acting replication element (CRE) in NS5B, and the 3' UTR. We propose an alternative structure in the 5' UTR, conserved secondary structures of 5B stem-loop (SL)1 and 5BSL2, and four possible structures of the X-tail at the very 3' end of the HCV genome. We predict several previously unknown long-range interactions, most importantly a possible circularization interaction between distinct elements in the 5' and 3' UTR, reminiscent of the cyclization elements of the related flaviviruses. Based on analogy to these viruses, we propose that the 5'-3' UTR base-pairing in the HCV genome might play an important role in viral RNA replication. These results may have important implications for our understanding of the nature of the cis-acting RNA elements in the HCV genome and their possible role in regulating the mutually exclusive processes of viral RNA translation and replication.

ITS2 Secondary Structure Improves Discrimination between Medicinal "Mu Tong" Species when Using DNA Barcoding

DNA barcoding is a promising species identification method, but it has proved difficult to find a standardized DNA marker in plant. Although the ITS/ITS2 RNA transcript has been proposed as the core barcode for seed plants, it has been criticized for being too conserved in some species to provide enough information or too variable in some species to align it within the different taxa ranks. We selected 30 individuals, representing 16 species and four families, to explore whether ITS2 can successfully resolve species in terms of secondary structure. Secondary structure was predicted using Mfold software and sequence-structure was aligned by MARNA. RNAstat software transformed the secondary structures into 28 symbol code data for maximum parsimony (MP) analysis. The results showed that the ITS2 structures in our samples had a common four-helix folding type with some shared motifs. This conserved structure facilitated the alignment of ambiguous sequences from divergent families. The structure alignment yielded a MP tree, in which most topological relationships were congruent with the tree constructed using nucleotide sequence data. When the data was combined, we obtained a well-resolved and highly supported phylogeny, in which individuals of a same species were clustered together into a monophyletic group. As a result, the different species that are often referred to as the herb “Mu tong” were successfully identified using short fragments of 250 bp ITS2 sequences, together with their secondary structure. Thus our analysis strengthens the potential of ITS2 as a promising DNA barcode because it incorporates valuable secondary structure information that will help improve discrimination between species.

Identifying novel sequence variants of RNA 3D motifs

Predicting RNA 3D structure from sequence is a major challenge in biophysics. An important sub-goal is accurately identifying recurrent 3D motifs from RNA internal and hairpin loop sequences extracted from secondary structure (2D) diagrams. We have developed and validated new probabilistic models for 3D motif sequences based on hybrid Stochastic Context-Free Grammars and Markov Random Fields (SCFG/MRF). The SCFG/MRF models are constructed using atomic-resolution RNA 3D structures. To parameterize each model, we use all instances of each motif found in the RNA 3D Motif Atlas and annotations of pairwise nucleotide interactions generated by the FR3D software. Isostericity relations between non-Watson–Crick basepairs are used in scoring sequence variants. SCFG techniques model nested pairs and insertions, while MRF ideas handle crossing interactions and base triples. We use test sets of randomly-generated sequences to set acceptance and rejection thresholds for each motif group and thus control the false positive rate. Validation was carried out by comparing results for four motif groups to RMDetect. The software developed for sequence scoring (JAR3D) is structured to automatically incorporate new motifs as they accumulate in the RNA 3D Motif Atlas when new structures are solved and is available free for download.

In silico discovery and modeling of non-coding RNA structure in viruses

This review covers several computational methods for discovering structured non-coding RNAs in viruses and modeling their putative secondary structures. Here we will use examples from two target viruses to highlight these approaches: influenza A virus-a relatively small, segmented RNA virus; and Epstein-Barr virus-a relatively large DNA virus with a complex transcriptome. Each system has unique challenges to overcome and unique characteristics to exploit. From these particular cases, generically useful approaches can be derived for the study of additional viral targets.

Multiple 3D RNA Structure Superposition Using Neighbor Joining

Recent advances in RNA research and the steady growth of available RNA structures call for bioinformatics methods for handling and analyzing RNA structural data. Recently, we introduced SETTER-a fast and accurate method for RNA pairwise structure alignment. In this paper, we describe MultiSETTER, SETTER extension for multiple RNA structure alignment. MultiSETTER combines SETTER's decomposition of RNA structures into non-overlapping structural subunits with the multiple sequence alignment algorithm ClustalW adapted for the structure alignment. The accuracy of MultiSETTER was assessed by the automatic classification of RNA structures and its comparison to SCOR annotations. In addition, MultiSETTER classification was also compared to multiple sequence alignment-based and secondary structure alignment-based classifications provided by LocARNA and RNADistance tools, respectively. MultiSETTER precompiled Windows libraries, as well as the C++ source code, are freely available from http://siret.cz/multisetter.

SPARSE: quadratic time simultaneous alignment and folding of RNAs without sequence-based heuristics

Motivation: RNA-Seq experiments have revealed a multitude of novel ncRNAs. The gold standard for their analysis based on simultaneous alignment and folding suffers from extreme time complexity of O(n6). Subsequently, numerous faster ‘Sankoff-style’ approaches have been suggested. Commonly, the performance of such methods relies on sequence-based heuristics that restrict the search space to optimal or near-optimal sequence alignments; however, the accuracy of sequence-based methods breaks down for RNAs with sequence identities below 60%. Alignment approaches like LocARNA that do not require sequence-based heuristics, have been limited to high complexity (≥ quartic time).

Results: Breaking this barrier, we introduce the novel Sankoff-style algorithm ‘sparsified prediction and alignment of RNAs based on their structure ensembles (SPARSE)’, which runs in quadratic time without sequence-based heuristics. To achieve this low complexity, on par with sequence alignment algorithms, SPARSE features strong sparsification based on structural properties of the RNA ensembles. Following PMcomp, SPARSE gains further speed-up from lightweight energy computation. Although all existing lightweight Sankoff-style methods restrict Sankoff’s original model by disallowing loop deletions and insertions, SPARSE transfers the Sankoff algorithm to the lightweight energy model completely for the first time. Compared with LocARNA, SPARSE achieves similar alignment and better folding quality in significantly less time (speedup: 3.7). At similar run-time, it aligns low sequence identity instances substantially more accurate than RAF, which uses sequence-based heuristics.

Availability and implementation: SPARSE is freely available at http://www.bioinf.uni-freiburg.de/Software/SPARSE.

Contact:backofen@informatik.uni-freiburg.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Small regulatory RNA-induced growth rate heterogeneity of Bacillus subtilis

Isogenic bacterial populations can consist of cells displaying heterogeneous physiological traits. Small regulatory RNAs (sRNAs) could affect this heterogeneity since they act by fine-tuning mRNA or protein levels to coordinate the appropriate cellular behavior. Here we show that the sRNA RnaC/S1022 from the Gram-positive bacterium Bacillus subtilis can suppress exponential growth by modulation of the transcriptional regulator AbrB. Specifically, the post-transcriptional abrB-RnaC/S1022 interaction allows B. subtilis to increase the cell-to-cell variation in AbrB protein levels, despite strong negative autoregulation of the abrB promoter. This behavior is consistent with existing mathematical models of sRNA action, thus suggesting that induction of protein expression noise could be a new general aspect of sRNA regulation. Importantly, we show that the sRNA-induced diversity in AbrB levels generates heterogeneity in growth rates during the exponential growth phase. Based on these findings, we hypothesize that the resulting subpopulations of fast- and slow-growing B. subtilis cells reflect a bet-hedging strategy for enhanced survival of unfavorable conditions.

Bacterial cells that share the same genetic information can display very different phenotypes, even if they grow under identical conditions. Despite the relevance of this population heterogeneity for processes like drug resistance and development, the molecular players that induce heterogenic phenotypes are often not known. Here we report that in the Gram-positive model bacterium Bacillus subtilis a small regulatory RNA (sRNA) can induce heterogeneity in growth rates by increasing cell-to-cell variation in the levels of the transcriptional regulator AbrB, which is important for rapid growth. Remarkably, the observed variation in AbrB levels is induced post-transcriptionally because of AbrB’s negative autoregulation, and is not observed at the abrB promoter level. We show that our observations are consistent with mathematical models of sRNA action, thus suggesting that induction of protein expression noise could be a new general aspect of sRNA regulation. Since a low growth rate can be beneficial for cellular survival, we propose that the observed subpopulations of fast- and slow-growing B. subtilis cells reflect a bet-hedging strategy for enhanced survival of unfavorable conditions.

Effective alignment of RNA pseudoknot structures using partition function posterior log-odds scores

Background

RNA pseudoknots play important roles in many biological processes. Previous methods for comparative pseudoknot analysis mainly focus on simultaneous folding and alignment of RNA sequences. Little work has been done to align two known RNA secondary structures with pseudoknots taking into account both sequence and structure information of the two RNAs.

Results

In this article we present a novel method for aligning two known RNA secondary structures with pseudoknots. We adopt the partition function methodology to calculate the posterior log-odds scores of the alignments between bases or base pairs of the two RNAs with a dynamic programming algorithm. The posterior log-odds scores are then used to calculate the expected accuracy of an alignment between the RNAs. The goal is to find an optimal alignment with the maximum expected accuracy. We present a heuristic to achieve this goal. The performance of our method is investigated and compared with existing tools for RNA structure alignment. An extension of the method to multiple alignment of pseudoknot structures is also discussed.

Conclusions

The method described here has been implemented in a tool named RKalign, which is freely accessible on the Internet. As more and more pseudoknots are revealed, collected and stored in public databases, we anticipate a tool like RKalign will play a significant role in data comparison, annotation, analysis, and retrieval in these databases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0464-9) contains supplementary material, which is available to authorized users.

The incidence and pattern of copollinator diversification in dioecious and monoecious figs

Differences in breeding system are associated with correlated ecological and morphological changes in plants. In Ficus, dioecy and monoecy are strongly associated with different suites of traits (tree height, population density, fruiting frequency, pollinator dispersal ecology). Although approximately 30% of fig species are pollinated by multiple species of fig-pollinating wasps, it has been suggested that copollinators are rare in dioecious figs. Here, we test whether there is a connection between the fig breeding system and copollinator incidence and diversification by conducting a meta-analysis of molecular data from pollinators of 119 fig species that includes new data from 15 Asian fig species. We find that the incidence of copollinators is not significantly different between monoecious and dioecious Ficus. Surprisingly, while all copollinators in dioecious figs are sister taxa, only 32.1% in monoecious figs are sister taxa. We present hypotheses to explain those patterns and discuss their consequences on the evolution of this mutualism.

An internal ribosome entry site in the 5' untranslated region of epidermal growth factor receptor allows hypoxic expression

The expression of epidermal growth factor receptor (EGFR/ERBB1/HER1) is implicated in the progress of numerous cancers, a feature that has been exploited in the development of EGFR antibodies and EGFR tyrosine kinase inhibitors as anti-cancer drugs. However, EGFR also has important normal cellular functions, leading to serious side effects when EGFR is inhibited. One damaging characteristic of many oncogenes is the ability to be expressed in the hypoxic conditions associated with the tumour interior. It has previously been demonstrated that expression of EGFR is maintained in hypoxic conditions via an unknown mechanism of translational control, despite global translation rates generally being attenuated under hypoxic conditions. In this report, we demonstrate that the human EGFR 5′ untranslated region (UTR) sequence can initiate the expression of a downstream open reading frame via an internal ribosome entry site (IRES). We show that this effect is not due to either cryptic promoter activity or splicing events. We have investigated the requirement of the EGFR IRES for eukaryotic initiation factor 4A (eIF4A), which is an RNA helicase responsible for processing RNA secondary structure as part of translation initiation. Treatment with hippuristanol (a potent inhibitor of eIF4A) caused a decrease in EGFR 5′ UTR-driven reporter activity and also a reduction in EGFR protein level. Importantly, we show that expression of a reporter gene under the control of the EGFR IRES is maintained under hypoxic conditions despite a fall in global translation rates.

CRISPRstrand: predicting repeat orientations to determine the crRNA-encoding strand at CRISPR loci

MOTIVATION

The discovery of CRISPR-Cas systems almost 20 years ago rapidly changed our perception of the bacterial and archaeal immune systems. CRISPR loci consist of several repetitive DNA sequences called repeats, inter-spaced by stretches of variable length sequences called spacers. This CRISPR array is transcribed and processed into multiple mature RNA species (crRNAs). A single crRNA is integrated into an interference complex, together with CRISPR-associated (Cas) proteins, to bind and degrade invading nucleic acids. Although existing bioinformatics tools can recognize CRISPR loci by their characteristic repeat-spacer architecture, they generally output CRISPR arrays of ambiguous orientation and thus do not determine the strand from which crRNAs are processed. Knowledge of the correct orientation is crucial for many tasks, including the classification of CRISPR conservation, the detection of leader regions, the identification of target sites (protospacers) on invading genetic elements and the characterization of protospacer-adjacent motifs.

RESULTS

We present a fast and accurate tool to determine the crRNA-encoding strand at CRISPR loci by predicting the correct orientation of repeats based on an advanced machine learning approach. Both the repeat sequence and mutation information were encoded and processed by an efficient graph kernel to learn higher-order correlations. The model was trained and tested on curated data comprising >4500 CRISPRs and yielded a remarkable performance of 0.95 AUC ROC (area under the curve of the receiver operator characteristic). In addition, we show that accurate orientation information greatly improved detection of conserved repeat sequence families and structure motifs. We integrated CRISPRstrand predictions into our CRISPRmap web server of CRISPR conservation and updated the latter to version 2.0.

AVAILABILITY

CRISPRmap and CRISPRstrand are available at http://rna.informatik.uni-freiburg.de/CRISPRmap.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

ExpaRNA-P: simultaneous exact pattern matching and folding of RNAs

Background

Identifying sequence-structure motifs common to two RNAs can speed up the comparison of structural RNAs substantially. The core algorithm of the existent approach ExpaRNA solves this problem for a priori known input structures. However, such structures are rarely known; moreover, predicting them computationally is no rescue, since single sequence structure prediction is highly unreliable.

Results

The novel algorithm ExpaRNA-P computes exactly matching sequence-structure motifs in entire Boltzmann-distributed structure ensembles of two RNAs; thereby we match and fold RNAs simultaneously, analogous to the well-known “simultaneous alignment and folding” of RNAs. While this implies much higher flexibility compared to ExpaRNA, ExpaRNA-P has the same very low complexity (quadratic in time and space), which is enabled by its novel structure ensemble-based sparsification. Furthermore, we devise a generalized chaining algorithm to compute compatible subsets of ExpaRNA-P’s sequence-structure motifs. Resulting in the very fast RNA alignment approach ExpLoc-P, we utilize the best chain as anchor constraints for the sequence-structure alignment tool LocARNA. ExpLoc-P is benchmarked in several variants and versus state-of-the-art approaches. In particular, we formally introduce and evaluate strict and relaxed variants of the problem; the latter makes the approach sensitive to compensatory mutations. Across a benchmark set of typical non-coding RNAs, ExpLoc-P has similar accuracy to LocARNA but is four times faster (in both variants), while it achieves a speed-up over 30-fold for the longest benchmark sequences (≈400nt). Finally, different ExpLoc-P variants enable tailoring of the method to specific application scenarios. ExpaRNA-P and ExpLoc-P are distributed as part of the LocARNA package. The source code is freely available at http://www.bioinf.uni-freiburg.de/Software/ExpaRNA-P.

Conclusions

ExpaRNA-P’s novel ensemble-based sparsification reduces its complexity to quadratic time and space. Thereby, ExpaRNA-P significantly speeds up sequence-structure alignment while maintaining the alignment quality. Different ExpaRNA-P variants support a wide range of applications.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0404-0) contains supplementary material, which is available to authorized users.

ITS1, 5.8S and ITS2 secondary structure modelling for intra-specific differentiation among species of the Colletotrichum gloeosporioides sensu lato species complex

The Colletotrichum gloeosporioides species complex is among the most destructive fungal plant pathogens in the world, however, identification of member species which are of quarantine importance is impacted by a number of factors that negatively affect species identification. Structural information of the rRNA marker may be considered to be a conserved marker which can be used as supplementary information for possible species identification. The difficulty in using ITS rDNA sequences for identification lies in the low level of sequence variation at the intra-specific level and the generation of artificially-induced sequence variation due to errors in polymerization of the ITS array during DNA replication. Type and query ITS sequences were subjected to sequence analyses prior to generation of predicted consensus secondary structures, including the pattern of nucleotide polymorphisms and number of indel haplotypes, GC content, and detection of artificially-induced sequence variation. Data pertaining to structure stability, the presence of conserved motifs in secondary structures and mapping of all sequences onto the consensus C. gloeosporioides sensu stricto secondary structure for ITS1, 5.8S and ITS2 markers was then carried out. Motifs that are evolutionarily conserved among eukaryotes were found for all ITS1, 5.8S and ITS2 sequences. The sequences exhibited conserved features typical of functional rRNAs. Generally, polymorphisms occurred within less conserved regions and were seen as bulges, internal and terminal loops or non-canonical G–U base-pairs within regions of the double stranded helices. Importantly, there were also taxonomic motifs and base changes that were unique to specific taxa and which may be used to support intra-specific identification of members of the C. gloeosporioides sensu lato species complex.

RNA structure analysis of alphacoronavirus terminal genome regions

•

Review of current knowledge of cis-acting RNA elements essential to coronavirus replication.

•

Identification of RNA structural elements in alphacoronavirus terminal genome regions.

•

Discussion of intra- and intergeneric conservation of genomic cis-acting RNA elements in alpha- and betacoronaviruses.

Coronavirus genome replication is mediated by a multi-subunit protein complex that is comprised of more than a dozen virally encoded and several cellular proteins. Interactions of the viral replicase complex with cis-acting RNA elements located in the 5′ and 3′-terminal genome regions ensure the specific replication of viral RNA. Over the past years, boundaries and structures of cis-acting RNA elements required for coronavirus genome replication have been extensively characterized in betacoronaviruses and, to a lesser extent, other coronavirus genera. Here, we review our current understanding of coronavirus cis-acting elements located in the terminal genome regions and use a combination of bioinformatic and RNA structure probing studies to identify and characterize putative cis-acting RNA elements in alphacoronaviruses. The study suggests significant RNA structure conservation among members of the genus Alphacoronavirus but also across genus boundaries. Overall, the conservation pattern identified for 5′ and 3′-terminal RNA structural elements in the genomes of alpha- and betacoronaviruses is in agreement with the widely used replicase polyprotein-based classification of the Coronavirinae, suggesting co-evolution of the coronavirus replication machinery with cognate cis-acting RNA elements.

Structure of the complete bacterial SRP Alu domain

The Alu domain of the signal recognition particle (SRP) arrests protein biosynthesis by competition with elongation factor binding on the ribosome. The mammalian Alu domain is a protein-RNA complex, while prokaryotic Alu domains are protein-free with significant extensions of the RNA. Here we report the crystal structure of the complete Alu domain of Bacillus subtilis SRP RNA at 2.5 Å resolution. The bacterial Alu RNA reveals a compact fold, which is stabilized by prokaryote-specific extensions and interactions. In this 'closed' conformation, the 5' and 3' regions are clamped together by the additional helix 1, the connecting 3-way junction and a novel minor groove interaction, which we term the 'minor-saddle motif' (MSM). The 5' region includes an extended loop-loop pseudoknot made of five consecutive Watson-Crick base pairs. Homology modeling with the human Alu domain in context of the ribosome shows that an additional lobe in the pseudoknot approaches the large subunit, while the absence of protein results in the detachment from the small subunit. Our findings provide the structural basis for purely RNA-driven elongation arrest in prokaryotes, and give insights into the structural adaption of SRP RNA during evolution.

Analysis of new aphid lethal paralysis virus (ALPV) isolates suggests evolution of two ALPV species

Aphid lethal paralysis virus (ALPV; family Dicistroviridae) was first isolated from the bird cherry-oat aphid, Rhopalosiphum padi. ALPV-like virus sequences have been reported from many insects and insect predators. We identified a new isolate of ALPV (ALPV-AP) from the pea aphid, Acyrthosiphon pisum, and a new isolate (ALPV-DvV) from western corn rootworm, Diabrotica virgifera virgifera. ALPV-AP has an ssRNA genome of 9940 nt. Based on phylogenetic analysis, ALPV-AP was closely related to ALPV-AM, an ALPV isolate from honeybees, Apis mellifera, in Spain and Brookings, SD, USA. The distinct evolutionary branches suggested the existence of two lineages of the ALPV virus. One consisted of ALPV-AP and ALPV-AM, whilst all other isolates of ALPV grouped into the other lineage. The similarity of ALPV-AP and ALPV-AM was up to 88 % at the RNA level, compared with 78-79 % between ALPV-AP and other ALPV isolates. The sequence identity of proteins between ALPV-AP and ALPV-AM was 98-99 % for both ORF1 and ORF2, whilst only 85-87 % for ORF1 and 91-92 % for ORF2 between ALPV-AP and other ALPV isolates. Sequencing of RACE (rapid amplification of cDNA ends) products and cDNA clones of the virus genome revealed sequence variation in the 5' UTRs and in ORF1, indicating that ALPV may be under strong selection pressure, which could have important biological implications for ALPV host range and infectivity. Our results indicated that ALPV-like viruses infect insects in the order Coleoptera, in addition to the orders Hemiptera and Hymenoptera, and we propose that ALPV isolates be classified as two separate viral species.

Human immunodeficiency virus Tat associates with a specific set of cellular RNAs

BACKGROUND

Human Immunodeficiency Virus 1 (HIV-1) exhibits a wide range of interactions with the host cell but whether viral proteins interact with cellular RNA is not clear. A candidate interacting factor is the trans-activator of transcription (Tat) protein. Tat is required for expression of virus genes but activates transcription through an unusual mechanism; binding to an RNA stem-loop, the transactivation response element (TAR), with the host elongation factor P-TEFb. HIV-1 Tat has also been shown to alter the expression of host genes during infection, contributing to viral pathogenesis but, whether Tat also interacts with cellular RNAs is unknown.

RESULTS

Using RNA immunoprecipitation coupled with microarray analysis, we have discovered that HIV-1 Tat is associated with a specific set of human mRNAs in T cells. mRNAs bound by Tat share a stem-loop structural element and encode proteins with common biological roles. In contrast, we do not find evidence that Tat associates with microRNAs or the RNA-induced silencing complex (RISC). The interaction of Tat with cellular RNA requires an intact RNA binding domain and Tat RNA binding is linked to an increase in RNA abundance in cell lines and during infection of primary CD4+ T cells by HIV.

CONCLUSIONS

We conclude that Tat interacts with a specific set of human mRNAs in T cells, many of which show changes in abundance in response to Tat and HIV infection. This work uncovers a previously unrecognised interaction between HIV and its host that may contribute to viral alteration of the host cellular environment.

ITS2, 18S, 16S or any other RNA - simply aligning sequences and their individual secondary structures simultaneously by an automatic approach

Secondary structures of RNA sequences are increasingly being used as additional information in reconstructing phylogenies and/or in distinguishing species by compensatory base change (CBC) analyses. However, in most cases just one secondary structure is used in manually correcting an automatically generated multiple sequence alignment and/or just one secondary structure is used in guiding a sequence alignment still completely generated by hand. With the advent of databases and tools offering individual RNA secondary structures, here we re-introduce a twelve letter code already implemented in 4SALE - a tool for synchronous sequence and secondary structure alignment and editing - that enables one to align RNA sequences and their individual secondary structures synchronously and fully automatic, while dramatically increasing the phylogenetic information content. We further introduce a scaled down non-GUI version of 4SALE particularly designed for big data analysis, and available at: http://4sale.bioapps.biozentrum.uni-wuerzburg.de.

Identification of non-coding RNAs with a new composite feature in the Hybrid Random Forest Ensemble algorithm

To identify non-coding RNA (ncRNA) signals within genomic regions, a classification tool was developed based on a hybrid random forest (RF) with a logistic regression model to efficiently discriminate short ncRNA sequences as well as long complex ncRNA sequences. This RF-based classifier was trained on a well-balanced dataset with a discriminative set of features and achieved an accuracy, sensitivity and specificity of 92.11%, 90.7% and 93.5%, respectively. The selected feature set includes a new proposed feature, SCORE. This feature is generated based on a logistic regression function that combines five significant features-structure, sequence, modularity, structural robustness and coding potential-to enable improved characterization of long ncRNA (lncRNA) elements. The use of SCORE improved the performance of the RF-based classifier in the identification of Rfam lncRNA families. A genome-wide ncRNA classification framework was applied to a wide variety of organisms, with an emphasis on those of economic, social, public health, environmental and agricultural significance, such as various bacteria genomes, the Arthrospira (Spirulina) genome, and rice and human genomic regions. Our framework was able to identify known ncRNAs with sensitivities of greater than 90% and 77.7% for prokaryotic and eukaryotic sequences, respectively. Our classifier is available at http://ncrna-pred.com/HLRF.htm.

Bioinformatics of prokaryotic RNAs

The genome of most prokaryotes gives rise to surprisingly complex transcriptomes, comprising not only protein-coding mRNAs, often organized as operons, but also harbors dozens or even hundreds of highly structured small regulatory RNAs and unexpectedly large levels of anti-sense transcripts. Comprehensive surveys of prokaryotic transcriptomes and the need to characterize also their non-coding components is heavily dependent on computational methods and workflows, many of which have been developed or at least adapted specifically for the use with bacterial and archaeal data. This review provides an overview on the state-of-the-art of RNA bioinformatics focusing on applications to prokaryotes.

Vicinal: a method for the determination of ncRNA ends using chimeric reads from RNA-seq experiments

Non-coding (nc)RNAs are important structural and regulatory molecules. Accurate determination of the primary sequence and secondary structure of ncRNAs is important for understanding their functions. During cDNA synthesis, RNA 3' end stem-loops can self-prime reverse transcription, creating RNA-cDNA chimeras. We found that chimeric RNA-cDNA fragments can also be detected at 5' end stem-loops, although at much lower frequency. Using the Gubler-Hoffman method, both types of chimeric fragments can be converted to cDNA during library construction, and they are readily detectable in high-throughput RNA sequencing (RNA-seq) experiments. Here, we show that these chimeric reads contain valuable information about the boundaries of ncRNAs. We developed a bioinformatic method, called Vicinal, to precisely map the ends of numerous fruitfly, mouse and human ncRNAs. Using this method, we analyzed chimeric reads from over 100 RNA-seq datasets, the results of which we make available for users to find RNAs of interest. In summary, we show that Vicinal is a useful tool for determination of the precise boundaries of uncharacterized ncRNAs, facilitating further structure/function studies.

Small non-coding RNAs in plant-pathogenic Xanthomonas spp

The genus Xanthomonas comprises a large group of plant-pathogenic bacteria. The infection and bacterial multiplication in the plant tissue depends on the type III secretion system and other virulence determinants. Recent studies revealed that bacterial virulence is also controlled at the post-transcriptional level by small non-coding RNAs (sRNAs). In this review, we highlight our current knowledge about sRNAs and RNA-binding proteins in Xanthomonas species.

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.