Glomerular-tubular crosstalk via cold shock Y-box binding protein-1 in the kidney.

Glomerular-tubular crosstalk within the kidney has been proposed, but the paracrine signals enabling this remain largely unknown. The cold-shock protein Y-box binding protein 1 (YBX1) is known to regulate inflammation and kidney diseases but its role in podocytes remains undetermined. Therefore, we analyzed mice with podocyte specific Ybx1 deletion (Ybx1ΔPod). Albuminuria was increased in unchallenged Ybx1ΔPod mice, which surprisingly was associated with reduced glomerular, but enhanced tubular damage. Tubular toll-like receptor 4 (TLR4) expression, node-like receptor protein 3 (NLRP3) inflammasome activation and kidney inflammatory cell infiltrates were all increased in Ybx1ΔPod mice. In vitro, extracellular YBX1 inhibited NLRP3 inflammasome activation in tubular cells. Co-immunoprecipitation, immunohistochemical analyses, microscale cell-free thermophoresis assays, and blunting of the YBX1-mediated TLR4-inhibition by a unique YBX1-derived decapeptide suggests a direct interaction of YBX1 and TLR4. Since YBX1 can be secreted upon post-translational acetylation, we hypothesized that YBX1 secreted from podocytes can inhibit TLR4 signaling in tubular cells. Indeed, mice expressing a non-secreted YBX1 variant specifically in podocytes (Ybx1PodK2A mice) phenocopied Ybx1ΔPod mice, demonstrating a tubular-protective effect of YBX1 secreted from podocytes. Lipopolysaccharide-induced tubular injury was aggravated in Ybx1ΔPod and Ybx1PodK2A mice, indicating a pathophysiological relevance of this glomerular-tubular crosstalk. Thus, our data show that YBX1 is physiologically secreted from podocytes, thereby negatively modulating sterile inflammation in the tubular compartment, apparently by binding to and inhibiting tubular TLR4 signaling. Hence, we have uncovered an YBX1-dependent molecular mechanism of glomerular-tubular crosstalk.

RNA interaction format: a general data format for RNA interactions.

SUMMARY: RNA molecules play crucial roles in various biological processes. They mediate their function mainly by interacting with other RNAs or proteins. At present, information about these interactions is distributed over different resources, often providing the data in simple tab-delimited formats that differ between the databases. There is no standardized data format that can capture the nature of all these different interactions in detail. AVAILABILITY AND IMPLEMENTATION: Here, we propose the RNA interaction format (RIF) for the detailed representation of RNA-RNA and RNA-Protein interactions and provide reference implementations in C/C++, Python, and JavaScript. RIF is released under licence GNU General Public License version 3 (GNU GPLv3) and is available on https://github.com/RNABioInfo/rna-interaction-format.

Application of RtcB ligase to monitor self-cleaving ribozyme activity by RNA-seq.

Self-cleaving ribozymes are catalytic RNAs and can be found in all domains of life. They catalyze a site-specific cleavage that results in a 5' fragment with a 2',3' cyclic phosphate (2',3' cP) and a 3' fragment with a 5' hydroxyl (5' OH) end. Recently, several strategies to enrich self-cleaving ribozymes by targeted biochemical methods have been introduced by us and others. Here, we develop an alternative strategy in which 5' OH RNAs are specifically ligated by RtcB ligase, which first guanylates the 3' phosphate of the adapter and then ligates it directly to RNAs with 5' OH ends. Our results demonstrate that adapter ligation to highly structured ribozyme fragments is much more efficient using the thermostable RtcB ligase from Pyrococcus horikoshii than the broadly applied Escherichia coli enzyme. Moreover, we investigated DNA, RNA and modified RNA adapters for their suitability in RtcB ligation reactions. We used the optimized RtcB-mediated ligation to produce RNA-seq libraries and captured a spiked 3' twister ribozyme fragment from E. coli total RNA. This RNA-seq-based method is applicable to detect ribozyme fragments as well as other cellular RNAs with 5' OH termini from total RNA.

Alzheimer-related genes show accelerated evolution.

Alzheimer's disease (AD) is a neurodegenerative disorder of unknown cause with complex genetic and environmental traits. While AD is extremely prevalent in human elderly, it hardly occurs in non-primate mammals and even non-human-primates develop only an incomplete form of the disease. This specificity of AD to human clearly implies a phylogenetic aspect. Still, the evolutionary dimension of AD pathomechanism remains difficult to prove and has not been established so far. To analyze the evolutionary age and dynamics of AD-associated-genes, we established the AD-associated genome-wide RNA-profile comprising both protein-coding and non-protein-coding transcripts. We than applied a systematic analysis on the conservation of splice-sites as a measure of gene-structure based on multiple alignments across vertebrates of homologs of AD-associated-genes. Here, we show that nearly all AD-associated-genes are evolutionarily old and did not originate later in evolution than not-AD-associated-genes. However, the gene-structures of loci, that exhibit AD-associated changes in their expression, evolve faster than the genome at large. While protein-coding-loci exhibit an enhanced rate of small changes in gene structure, non-coding loci show even much larger changes. The accelerated evolution of AD-associated-genes indicates a more rapid functional adaptation of these genes. In particular AD-associated non-coding-genes play an important, as yet largely unexplored, role in AD. This phylogenetic trait indicates that recent adaptive evolution of human brain is causally involved in basic principles of neurodegeneration. It highlights the necessity for a paradigmatic change of our disease-concepts and to reconsider the appropriateness of current animal-models to develop disease-modifying strategies that can be translated to human.

In vitro iCLIP-based modeling uncovers how the splicing factor U2AF2 relies on regulation by cofactors.

Alternative splicing generates distinct mRNA isoforms and is crucial for proteome diversity in eukaryotes. The RNA-binding protein (RBP) U2AF2 is central to splicing decisions, as it recognizes 3' splice sites and recruits the spliceosome. We establish "in vitro iCLIP" experiments, in which recombinant RBPs are incubated with long transcripts, to study how U2AF2 recognizes RNA sequences and how this is modulated by trans-acting RBPs. We measure U2AF2 affinities at hundreds of binding sites and compare in vitro and in vivo binding landscapes by mathematical modeling. We find that trans-acting RBPs extensively regulate U2AF2 binding in vivo, including enhanced recruitment to 3' splice sites and clearance of introns. Using machine learning, we identify and experimentally validate novel trans-acting RBPs (including FUBP1, CELF6, and PCBP1) that modulate U2AF2 binding and affect splicing outcomes. Our study offers a blueprint for the high-throughput characterization of in vitro mRNP assembly and in vivo splicing regulation.

cyPhyRNA-seq: a genome-scale RNA-seq method to detect active self-cleaving ribozymes by capturing RNAs with 2',3' cyclic phosphates and 5' hydroxyl ends.

Self-cleaving ribozymes are catalytically active RNAs that cleave themselves into a 5'-fragment with a 2',3'-cyclic phosphate and a 3'-fragment with a 5'-hydroxyl. They are widely applied for the construction of synthetic RNA devices and RNA-based therapeutics. However, the targeted discovery of self-cleaving ribozymes remains a major challenge. We developed a transcriptome-wide method, called cyPhyRNA-seq, to screen for ribozyme cleavage fragments in total RNA extract. This approach employs the specific ligation-based capture of ribozyme 5'-fragments using a variant of the Arabidopsis thaliana tRNA ligase we engineered. To capture ribozyme 3'-fragments, they are enriched from total RNA by enzymatic treatments. We optimized and enhanced the individual steps of cyPhyRNA-seq in vitro and in spike-in experiments. Then, we applied cyPhyRNA-seq to total RNA isolated from the bacterium Desulfovibrio vulgaris and detected self-cleavage of the three predicted type II hammerhead ribozymes, whose activity had not been examined to date. cyPhyRNA-seq can be used for the global analysis of active self-cleaving ribozymes with the advantage to capture both ribozyme cleavage fragments from total RNA. Especially in organisms harbouring many self-cleaving RNAs, cyPhyRNA-seq facilitates the investigation of cleavage activity. Moreover, this method has the potential to be used to discover novel self-cleaving ribozymes in different organisms. [Figure: see text].

MCPIP1 ribonuclease can bind and cleave AURKA mRNA in MYCN-amplified neuroblastoma cells.

The role of the inflammation-silencing ribonuclease, MCPIP1 (monocyte chemoattractant protein-induced protein 1), in neoplasia continuous to emerge. The ribonuclease can cleave not only inflammation-related transcripts but also some microRNAs (miRNAs) and viral RNAs. The suppressive effect of the protein has been hitherto suggested in breast cancer, clear cell renal cell carcinoma, osteosarcoma, and neuroblastoma. Our previous results have demonstrated a reduced levels of several oncogenes, as well as inhibited growth of neuroblastoma cells upon MCPIP1 overexpression. Here, we investigate the mechanisms underlying the suppression of MYCN proto-oncogene, bHLH transcription factor (MYCN)-amplified neuroblastoma cells overexpressing the MCPIP1 protein. We showed that the levels of several transcripts involved in cell cycle progression decreased in BE(2)-C and KELLY cells overexpressing MCPIP1 in a ribonucleolytic activity-dependent manner. However, RNA immunoprecipitation indicated that only AURKA mRNA (encoding for Aurora A kinase) interacts with the ribonuclease. Furthermore, the application of a luciferase assay suggested MCPIP1-dependent destabilization of the transcript. Further analyses demonstrated that the entire conserved region of AURKA seems to be indispensable for the interaction with the MCPIP1 protein. Additionally, we examined the effect of the ribonuclease overexpression on the miRNA expression profile in MYCN-amplified neuroblastoma cells. However, no significant alterations were observed. Our data indicate a key role of the binding and cleavage of the AURKA transcript in an MCPIP1-dependent suppressive effect on neuroblastoma cells.

The RNA workbench 2.0: next generation RNA data analysis.

RNA has become one of the major research topics in molecular biology. As a central player in key processes regulating gene expression, RNA is in the focus of many efforts to decipher the pathways that govern the transition of genetic information to a fully functional cell. As more and more researchers join this endeavour, there is a rapidly growing demand for comprehensive collections of tools that cover the diverse layers of RNA-related research. However, increasing amounts of data, from diverse types of experiments, addressing different aspects of biological questions need to be consolidated and integrated into a single framework. Only then is it possible to connect findings from e.g. RNA-Seq experiments and methods for e.g. target predictions. To address these needs, we present the RNA Workbench 2.0 , an updated online resource for RNA related analysis. With the RNA Workbench we created a comprehensive set of analysis tools and workflows that enables researchers to analyze their data without the need for sophisticated command-line skills. This update takes the established framework to the next level, providing not only a containerized infrastructure for analysis, but also a ready-to-use platform for hands-on training, analysis, data exploration, and visualization. The new framework is available at https://rna.usegalaxy.eu , and login is free and open to all users. The containerized version can be found at https://github.com/bgruening/galaxy-rna-workbench.

Adipose tissue derived bacteria are associated with inflammation in obesity and type 2 diabetes.

OBJECTIVE: Bacterial translocation to various organs including human adipose tissue (AT) due to increased intestinal permeability remains poorly understood. We hypothesised that: (1) bacterial presence is highly tissue specific and (2) related in composition and quantity to immune inflammatory and metabolic burden. DESIGN: We quantified and sequenced the bacterial 16S rRNA gene in blood and AT samples (omental, mesenteric and subcutaneous) of 75 subjects with obesity with or without type 2 diabetes (T2D) and used catalysed reporter deposition (CARD) - fluorescence in situ hybridisation (FISH) to detect bacteria in AT. RESULTS: Under stringent experimental and bioinformatic control for contaminants, bacterial DNA was detected in blood and omental, subcutaneous and mesenteric AT samples in the range of 0.1 to 5 pg/µg DNA isolate. Moreover, CARD-FISH allowed the detection of living, AT-borne bacteria. Proteobacteria and Firmicutes were the predominant phyla, and bacterial quantity was associated with immune cell infiltration, inflammatory and metabolic parameters in a tissue-specific manner. Bacterial composition differed between subjects with and without T2D and was associated with related clinical measures, including systemic and tissues-specific inflammatory markers. Finally, treatment of adipocytes with bacterial DNA in vitro stimulated the expression of TNFA and IL6. CONCLUSIONS: Our study provides contaminant aware evidence for the presence of bacteria and bacterial DNA in several ATs in obesity and T2D and suggests an important role of bacteria in initiating and sustaining local AT subclinical inflammation and therefore impacting metabolic sequelae of obesity.

Automatic curation of large comparative animal MicroRNA datasets.

MOTIVATION: MicroRNAs form an important class of RNA regulators that has been studied extensively. The miRBase and Rfam database provide rich, frequently updated information on both pre-miRNAs and their mature forms. These data sources, however, rely on individual data submission and thus are neither complete nor consistent in their coverage across different miRNA families. Quantitative studies of miRNA evolution therefore are difficult or impossible on this basis. RESULTS: We present here a workflow and a corresponding implementation, MIRfix, that automatically curates miRNA datasets by improving alignments of their precursors, the consistency of the annotation of mature miR and miR* sequence, and the phylogenetic coverage. MIRfix produces alignments that are comparable across families and sets the stage for improved homology search as well as quantitative analyses. AVAILABILITY AND IMPLEMENTATION: MIRfix can be downloaded from https://github.com/Bierinformatik/MIRfix. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

LOTTE-seq (Long hairpin oligonucleotide based tRNA high-throughput sequencing): specific selection of tRNAs with 3'-CCA end for high-throughput sequencing.

Transfer RNAs belong to the most abundant type of ribonucleic acid in the cell, and detailed investigations revealed correlations between alterations in the tRNA pool composition and certain diseases like breast cancer. However, currently available methods do not sample the entire tRNA pool or lack specificity for tRNAs. A specific disadvantage of such methods is that only full-length tRNAs are analysed, while tRNA fragments or incomplete cDNAs due to RT stops at modified nucleosides are lost. Another drawback in certain approaches is that the tRNA fraction has to be isolated and separated from high molecular weight RNA, resulting in considerable labour costs and loss of material. Based on a hairpin-shaped adapter oligonucleotide selective for tRNA transcripts, we developed a highly specific protocol for efficient and comprehensive high-throughput analysis of tRNAs that combines the benefits of existing methods and eliminates their disadvantages. Due to a 3'-TGG overhang, the adapter is specifically ligated to the tRNA 3'-CCA end. Reverse transcription prior to the ligation of a second adapter allows to include prematurely terminated cDNA products, increasing the number of tRNA reads. This strategy renders this approach a powerful and universal tool to analyse the tRNA pool of cells and organisms under different conditions in health and disease.

Unusual Occurrence of Two Bona-Fide CCA-Adding Enzymes in Dictyostelium discoideum.

Dictyostelium discoideum, the model organism for the evolutionary supergroup of Amoebozoa, is a social amoeba that, upon starvation, undergoes transition from a unicellular to a multicellular organism. In its genome, we identified two genes encoding for tRNA nucleotidyltransferases. Such pairs of tRNA nucleotidyltransferases usually represent collaborating partial activities catalyzing CC- and A-addition to the tRNA 3'-end, respectively. In D. discoideum, however, both enzymes exhibit identical activities, representing bona-fide CCA-adding enzymes. Detailed characterization of the corresponding activities revealed that both enzymes seem to be essential and are regulated inversely during different developmental stages of D. discoideum. Intriguingly, this is the first description of two functionally equivalent CCA-adding enzymes using the same set of tRNAs and showing a similar distribution within the cell. This situation seems to be a common feature in Dictyostelia, as other members of this phylum carry similar pairs of tRNA nucleotidyltransferase genes in their genome.

Accurate mapping of tRNA reads.

Motivation: Many repetitive DNA elements are transcribed at appreciable expression levels. Mapping the corresponding RNA sequencing reads back to a reference genome is notoriously difficult and error-prone task, however. This is in particular true if chemical modifications introduce systematic mismatches, while at the same time the genomic loci are only approximately identical, as in the case of tRNAs. Results: We therefore developed a dedicated mapping strategy to handle RNA-seq reads that map to tRNAs relying on a modified target genome in which known tRNA loci are masked and instead intronless tRNA precursor sequences are appended as artificial 'chromosomes'. In a first pass, reads that overlap the boundaries of mature tRNAs are extracted. In the second pass, the remaining reads are mapped to a tRNA-masked target that is augmented by representative mature tRNA sequences. Using both simulated and real life data we show that our best-practice workflow removes most of the mapping artefacts introduced by simpler mapping schemes and makes it possible to reliably identify many of chemical tRNA modifications in generic small RNA-seq data. Using simulated data the FDR is only 2%. We find compelling evidence for tissue specific differences of tRNA modification patterns. Availability and implementation: The workflow is available both as a bash script and as a Galaxy workflow from https://github.com/AnneHoffmann/tRNA-read-mapping. Contact: fabian@tbi.univie.ac.at. Supplementary information: Supplementary data are available at Bioinformatics online.

The RNA workbench: best practices for RNA and high-throughput sequencing bioinformatics in Galaxy.

RNA-based regulation has become a major research topic in molecular biology. The analysis of epigenetic and expression data is therefore incomplete if RNA-based regulation is not taken into account. Thus, it is increasingly important but not yet standard to combine RNA-centric data and analysis tools with other types of experimental data such as RNA-seq or ChIP-seq. Here, we present the RNA workbench, a comprehensive set of analysis tools and consolidated workflows that enable the researcher to combine these two worlds. Based on the Galaxy framework the workbench guarantees simple access, easy extension, flexible adaption to personal and security needs, and sophisticated analyses that are independent of command-line knowledge. Currently, it includes more than 50 bioinformatics tools that are dedicated to different research areas of RNA biology including RNA structure analysis, RNA alignment, RNA annotation, RNA-protein interaction, ribosome profiling, RNA-seq analysis and RNA target prediction. The workbench is developed and maintained by experts in RNA bioinformatics and the Galaxy framework. Together with the growing community evolving around this workbench, we are committed to keep the workbench up-to-date for future standards and needs, providing researchers with a reliable and robust framework for RNA data analysis. AVAILABILITY: The RNA workbench is available at https://github.com/bgruening/galaxy-rna-workbench.

AREsite2: an enhanced database for the comprehensive investigation of AU/GU/U-rich elements.

AREsite2 represents an update for AREsite, an on-line resource for the investigation of AU-rich elements (ARE) in human and mouse mRNA 3'UTR sequences. The new updated and enhanced version allows detailed investigation of AU, GU and U-rich elements (ARE, GRE, URE) in the transcriptome of Homo sapiens, Mus musculus, Danio rerio, Caenorhabditis elegans and Drosophila melanogaster. It contains information on genomic location, genic context, RNA secondary structure context and conservation of annotated motifs. Improvements include annotation of motifs not only in 3'UTRs but in the whole gene body including introns, additional genomes, and locally stable secondary structures from genome wide scans. Furthermore, we include data from CLIP-Seq experiments in order to highlight motifs with validated protein interaction. Additionally, we provide a REST interface for experienced users to interact with the database in a semi-automated manner. The database is publicly available at: http://rna.tbi.univie.ac.at/AREsite.

Tristetraprolin binding site atlas in the macrophage transcriptome reveals a switch for inflammation resolution.

Precise regulation of mRNA decay is fundamental for robust yet not exaggerated inflammatory responses to pathogens. However, a global model integrating regulation and functional consequences of inflammation-associated mRNA decay remains to be established. Using time-resolved high-resolution RNA binding analysis of the mRNA-destabilizing protein tristetraprolin (TTP), an inflammation-limiting factor, we qualitatively and quantitatively characterize TTP binding positions in the transcriptome of immunostimulated macrophages. We identify pervasive destabilizing and non-destabilizing TTP binding, including a robust intronic binding, showing that TTP binding is not sufficient for mRNA destabilization. A low degree of flanking RNA structuredness distinguishes occupied from silent binding motifs. By functionally relating TTP binding sites to mRNA stability and levels, we identify a TTP-controlled switch for the transition from inflammatory into the resolution phase of the macrophage immune response. Mapping of binding positions of the mRNA-stabilizing protein HuR reveals little target and functional overlap with TTP, implying a limited co-regulation of inflammatory mRNA decay by these proteins. Our study establishes a functionally annotated and navigable transcriptome-wide atlas (http://ttp-atlas.univie.ac.at) of cis-acting elements controlling mRNA decay in inflammation.

AREsite: a database for the comprehensive investigation of AU-rich elements.

AREsite is an online resource for the detailed investigation of AU-rich elements (ARE) in vertebrate mRNA 3'-untranslated regions (UTRs). AREs are one of the most prominent cis-acting regulatory elements found in 3'-UTRs of mRNAs. Various ARE-binding proteins that possess RNA stabilizing or destabilizing functions are recruited by sequence-specific motifs. Recent findings suggest an essential role of the structural mRNA context in which these sequence motifs are embedded. AREsite is the first database that allows to quantify the structuredness of ARE motif sites in terms of opening energies and accessibility probabilities. Moreover, we also provide a detailed phylogenetic analysis of ARE motifs and incorporate information about experimentally validated targets of the ARE-binding proteins TTP, HuR and Auf1. The database is publicly available at: http://rna.tbi.univie.ac.at/AREsite.

Toward a Systematic Assessment of Sex Differences in Cystic Fibrosis.

(1) Background: Cystic fibrosis (CF) is a disease with well-documented clinical differences between female and male patients. However, this gender gap is very poorly studied at the molecular level. (2) Methods: Expression differences in whole blood transcriptomics between female and male CF patients are analyzed in order to determine the pathways related to sex-biased genes and assess their potential influence on sex-specific effects in CF patients. (3) Results: We identify sex-biased genes in female and male CF patients and provide explanations for some sex-specific differences at the molecular level. (4) Conclusion: Genes in key pathways associated with CF are differentially expressed between sexes, and thus may account for the gender gap in morbidity and mortality in CF.

Tailored machine learning models for functional RNA detection in genome-wide screens.

The in silico prediction of non-coding and protein-coding genetic loci has received considerable attention in comparative genomics aiming in particular at the identification of properties of nucleotide sequences that are informative of their biological role in the cell. We present here a software framework for the alignment-based training, evaluation and application of machine learning models with user-defined parameters. Instead of focusing on the one-size-fits-all approach of pervasive in silico annotation pipelines, we offer a framework for the structured generation and evaluation of models based on arbitrary features and input data, focusing on stable and explainable results. Furthermore, we showcase the usage of our software package in a full-genome screen of Drosophila melanogaster and evaluate our results against the well-known but much less flexible program RNAz.