Duplicated RNA genes in teleost fish genomes

Phylogenetic and Chemical Probing Information as Soft Constraints in RNA Secondary Structure Prediction

Extrinsic, experimental information can be incorporated into thermodynamics-based RNA folding algorithms in the form of pseudo-energies. Evolutionary conservation of RNA secondary structure elements is detectable in alignments of phylogenetically related sequences and provides evidence for the presence of certain base pairs that can also be converted into pseudo-energy contributions. We show that the centroid base pairs computed from a consensus folding model such as RNAalifold result in a substantial improvement of the prediction accuracy for single sequences. Evidence for specific base pairs turns out to be more informative than a position-wise profile for the conservation of the pairing status. A comparison with chemical probing data, furthermore, strongly suggests that phylogenetic base pairing data are more informative than position-specific data on (un)pairedness as obtained from chemical probing experiments. In this context we demonstrate, in addition, that the conversion of signal from probing data into pseudo-energies is possible using thermodynamic structure predictions as a reference instead of known RNA structures.

Under control: The innate immunity of fish from the inhibitors' perspective

The innate immune response involves a concerted network of induced gene products, preformed immune effectors, biochemical signalling cascades and specialised cells. However, the multifaceted activation of these defensive measures can derail or overshoot and, if left unchecked, overwhelm the host. A plenty of regulatory devices therefore mediate the fragile equilibrium between pathogen defence and pathophysiological manifestations. Over the past decade in particular, an almost complete set of teleostean sequences orthologous to mammalian immunoregulatory factors has been identified in various fish species, which prove the remarkable conservation of innate immune-control concepts among vertebrates. This review will present the current knowledge on more than 50 teleostean regulatory factors (plus additional fish-specific paralogs) that are of paramount importance for controlling the clotting cascade, the complement system, pattern-recognition pathways and cytokine-signalling networks. A special focus lies on those immunoregulatory features that have emerged as potential biomarker genes in transcriptome-wide research studies. Moreover, we report on the latest progress in elucidating control elements that act directly with immune-gene-encoding nucleic acids, such as transcription factors, hormone receptors and micro- and long noncoding RNAs. Investigations into the function of teleostean inhibitory factors are still mainly based on gene-expression profiling or overexpression studies. However, in support of structural and in-vitro analyses, evidence from in-vivo trials is also available and revealed many biochemical details on piscine immune regulation. The presence of multiple gene copies in fish adds a degree of complexity, as it is so far hardly understood if they might play distinct roles during inflammation. The present review addresses this and other open questions that should be tackled by fish immunologists in future.

Metabolic consequences of microRNA-122 inhibition in rainbow trout, Oncorhynchus mykiss

Background

MicroRNAs (miRNAs) are small regulatory molecules which post-transcriptionally regulate mRNA stability and translation. Several microRNAs have received attention due to their role as key metabolic regulators. In spite of the high evolutionary conservation of several miRNAs, the role of miRNAs in lower taxa of vertebrates has not been studied with regard to metabolism. The liver-specific and highly abundant miRNA-122 is one of the most widely studied miRNA in mammals, where it has been implicated in the control of hepatic lipid metabolism. Following our identification of acute postprandial, nutritional and endocrine regulation of hepatic miRNA-122 isomiRNA expression in rainbow trout, we used complementary in silico and in vivo approaches to study the role of miRNA-122 in rainbow trout metabolism. We hypothesized that the role of miRNA-122 in regulating lipid metabolism in rainbow trout is conserved to that in mammals and that modulation of miRNA-122 function would result in altered lipid homeostasis and secondarily altered glucose homeostasis, since lipogenesis has been suggested to act as glucose sink in trout.

Results

Our results show that miRNA-122 was functionally inhibited in vivo in the liver. Postprandial glucose concentrations increased significantly in rainbow trout injected with a miRNA-122 inhibitor, and this effect correlated with decreases in hepatic FAS protein abundance, indicative of altered lipogenic potential. Additionally, miRNA-122 inhibition resulted in a 20% decrease in plasma cholesterol concentration, an effect associated with increased expression of genes involved in cholesterol degradation and excretion.

Conclusions

Overall evidence suggests that miRNA-122 may have evolved in early vertebrates to support liver-specific metabolic functions. Nevertheless, our data also indicate that metabolic consequences of miRNA-122 inhibition may differ quantitatively between vertebrate species and that distinct direct molecular targets of miRNA-122 may mediate metabolic effects between vertebrate species, indicating that miRNA-122 - mRNA target relationships may have undergone species-specific evolutionary changes.

CRISPRmap: an automated classification of repeat conservation in prokaryotic adaptive immune systems

Central to Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-Cas systems are repeated RNA sequences that serve as Cas-protein-binding templates. Classification is based on the architectural composition of associated Cas proteins, considering repeat evolution is essential to complete the picture. We compiled the largest data set of CRISPRs to date, performed comprehensive, independent clustering analyses and identified a novel set of 40 conserved sequence families and 33 potential structure motifs for Cas-endoribonucleases with some distinct conservation patterns. Evolutionary relationships are presented as a hierarchical map of sequence and structure similarities for both a quick and detailed insight into the diversity of CRISPR-Cas systems. In a comparison with Cas-subtypes, I-C, I-E, I-F and type II were strongly coupled and the remaining type I and type III subtypes were loosely coupled to repeat and Cas1 evolution, respectively. Subtypes with a strong link to CRISPR evolution were almost exclusive to bacteria; nevertheless, we identified rare examples of potential horizontal transfer of I-C and I-E systems into archaeal organisms. Our easy-to-use web server provides an automated assignment of newly sequenced CRISPRs to our classification system and enables more informed choices on future hypotheses in CRISPR-Cas research: http://rna.informatik.uni-freiburg.de/CRISPRmap.

GraphClust: alignment-free structural clustering of local RNA secondary structures

Motivation: Clustering according to sequence–structure similarity has now become a generally accepted scheme for ncRNA annotation. Its application to complete genomic sequences as well as whole transcriptomes is therefore desirable but hindered by extremely high computational costs.

Results: We present a novel linear-time, alignment-free method for comparing and clustering RNAs according to sequence and structure. The approach scales to datasets of hundreds of thousands of sequences. The quality of the retrieved clusters has been benchmarked against known ncRNA datasets and is comparable to state-of-the-art sequence–structure methods although achieving speedups of several orders of magnitude. A selection of applications aiming at the detection of novel structural ncRNAs are presented. Exemplarily, we predicted local structural elements specific to lincRNAs likely functionally associating involved transcripts to vital processes of the human nervous system. In total, we predicted 349 local structural RNA elements.

Availability: The GraphClust pipeline is available on request.

Contact:backofen@informatik.uni-freiburg.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Molecular evolution of the non-coding eosinophil granule ontogeny transcript

Eukaryotic genomes are pervasively transcribed. A large fraction of the transcriptional output consists of long, mRNA-like, non-protein-coding transcripts (mlncRNAs). The evolutionary history of mlncRNAs is still largely uncharted territory. In this contribution, we explore in detail the evolutionary traces of the eosinophil granule ontogeny transcript (EGOT), an experimentally confirmed representative of an abundant class of totally intronic non-coding transcripts (TINs). EGOT is located antisense to an intron of the ITPR1 gene. We computationally identify putative EGOT orthologs in the genomes of 32 different amniotes, including orthologs from primates, rodents, ungulates, carnivores, afrotherians, and xenarthrans, as well as putative candidates from basal amniotes, such as opossum or platypus. We investigate the EGOT gene phylogeny, analyze patterns of sequence conservation, and the evolutionary conservation of the EGOT gene structure. We show that EGO-B, the spliced isoform, may be present throughout the placental mammals, but most likely dates back even further. We demonstrate here for the first time that the whole EGOT locus is highly structured, containing several evolutionary conserved, and thermodynamic stable secondary structures. Our analyses allow us to postulate novel functional roles of a hitherto poorly understood region at the intron of EGO-B which is highly conserved at the sequence level. The region contains a novel ITPR1 exon and also conserved RNA secondary structures together with a conserved TATA-like element, which putatively acts as a promoter of an independent regulatory element.

Hox clusters of the bichir (Actinopterygii, Polypterus senegalus) highlight unique patterns of sequence evolution in gnathostome phylogeny

Teleost fishes have extra Hox gene clusters owing to shared or lineage-specific genome duplication events in rayfinned fish (actinopterygian) phylogeny. Hence, extrapolating between genome function of teleosts and human or even between different fish species is difficult. We have sequenced and analyzed Hox gene clusters of the Senegal bichir (Polypterus senegalus), an extant representative of the most basal actinopterygian lineage. Bichir possesses four Hox gene clusters (A, B, C, D); phylogenetic analysis supports their orthology to the four Hox gene clusters of the gnathostome ancestor. We have generated a comprehensive database of conserved Hox noncoding sequences that include cartilaginous, lobe-finned, and ray-finned fishes (bichir and teleosts). Our analysis identified putative and known Hox cis-regulatory sequences with differing depths of conservation in Gnathostoma. We found that although bichir possesses four Hox gene clusters, its pattern of conservation of noncoding sequences is mosaic between outgroups, such as human, coelacanth, and shark, with four Hox gene clusters and teleosts, such as zebrafish and pufferfish, with seven or eight Hox gene clusters. Notably, bichir Hox gene clusters have been invaded by DNA transposons and this trend is further exemplified in teleosts, suggesting an as yet unrecognized mechanism of genome evolution that may explain Hox cluster plasticity in actinopterygians. Taken together, our results suggest that actinopterygian Hox gene clusters experienced a reduction in selective constraints that surprisingly predates the teleost-specific genome duplication.

Bioinformatics for RNomics

Rapid improvements in high-throughput experimental technologies make it nowadays possible to study the expression, as well as changes in expression, of whole transcriptomes under different environmental conditions in a detailed view. We describe current approaches to identify genome-wide functional RNA transcripts (experimentally as well as computationally), and focus on computational methods that may be utilized to disclose their function. While genome databases offer a wealth of information about known and putative functions for protein-coding genes, functional information for novel non-coding RNA genes is almost nonexistent. This is mainly explained by the lack of established software tools to efficiently reveal the function and evolutionary origin of non-coding RNA genes. Here, we describe in detail computational approaches one may follow to annotate and classify an RNA transcript.

Non-coding RNA annotation of the genome of Trichoplax adhaerens

A detailed annotation of non-protein coding RNAs is typically missing in initial releases of newly sequenced genomes. Here we report on a comprehensive ncRNA annotation of the genome of Trichoplax adhaerens, the presumably most basal metazoan whose genome has been published to-date. Since blast identified only a small fraction of the best-conserved ncRNAs--in particular rRNAs, tRNAs and some snRNAs--we developed a semi-global dynamic programming tool, GotohScan, to increase the sensitivity of the homology search. It successfully identified the full complement of major and minor spliceosomal snRNAs, the genes for RNase P and MRP RNAs, the SRP RNA, as well as several small nucleolar RNAs. We did not find any microRNA candidates homologous to known eumetazoan sequences. Interestingly, most ncRNAs, including the pol-III transcripts, appear as single-copy genes or with very small copy numbers in the Trichoplax genome.