SEAweb: the small RNA Expression Atlas web application

Pathogen detection in RNA-seq data with Pathonoia

BACKGROUND

Bacterial and viral infections may cause or exacerbate various human diseases and to detect microbes in tissue, one method of choice is RNA sequencing. The detection of specific microbes using RNA sequencing offers good sensitivity and specificity, but untargeted approaches suffer from high false positive rates and a lack of sensitivity for lowly abundant organisms.

RESULTS

We introduce Pathonoia, an algorithm that detects viruses and bacteria in RNA sequencing data with high precision and recall. Pathonoia first applies an established k-mer based method for species identification and then aggregates this evidence over all reads in a sample. In addition, we provide an easy-to-use analysis framework that highlights potential microbe-host interactions by correlating the microbial to the host gene expression. Pathonoia outperforms state-of-the-art methods in microbial detection specificity, both on in silico and real datasets.

CONCLUSION

Two case studies in human liver and brain show how Pathonoia can support novel hypotheses on microbial infection exacerbating disease. The Python package for Pathonoia sample analysis and a guided analysis Jupyter notebook for bulk RNAseq datasets are available on GitHub.

isomiRdb: microRNA expression at isoform resolution

A significant fraction of mature miRNA transcripts carries sequence and/or length variations, termed isomiRs. IsomiRs are differentially abundant in cell types, tissues, body fluids or patients' samples. Not surprisingly, multiple studies describe a physiological and pathophysiological role. Despite their importance, systematically collected and annotated isomiR information available in databases remains limited. We thus developed isomiRdb, a comprehensive resource that compiles miRNA expression data at isomiR resolution from various sources. We processed 42 499 human miRNA-seq datasets (5.9 × 1011 sequencing reads) and consistently analyzed them using miRMaster and sRNAbench. Our database provides online access to the 90 483 most abundant isomiRs (>1 RPM in at least 1% of the samples) from 52 tissues and 188 cell types. Additionally, the full set of over 3 million detected isomiRs is available for download. Our resource can be queried at the sample, miRNA or isomiR level so users can quickly answer common questions about the presence/absence of a particular miRNA/isomiR in tissues of interest. Further, the database facilitates to identify whether a potentially interesting new isoform has been detected before and its frequency. In addition to expression tables, isomiRdb can generate multiple interactive visualisations including violin plots and heatmaps. isomiRdb is free to use and publicly available at: https://www.ccb.uni-saarland.de/isomirdb.

A Novel Tissue Atlas and Online Tool for the Interrogation of Small RNA Expression in Human Tissues and Biofluids

One promising goal for utilizing the molecular information circulating in biofluids is the discovery of clinically useful biomarkers. Extracellular RNAs (exRNAs) are one of the most diverse classes of molecular cargo, easily assayed by sequencing and with expressions that rapidly change in response to subject status. Despite diverse exRNA cargo, most evaluations from biofluids have focused on small RNA sequencing and analysis, specifically on microRNAs (miRNAs). Another goal of characterizing circulating molecular information, is to correlate expression to injuries associated with specific tissues of origin. Biomarker candidates are often described as being specific, enriched in a particular tissue or associated with a disease process. Likewise, miRNA data is often reported to be specific, enriched for a tissue, without rigorous testing to support the claim. Here we provide a tissue atlas of small RNAs from 30 different tissues and three different blood cell types. We analyzed the tissues for enrichment of small RNA sequences and assessed their expression in biofluids: plasma, cerebrospinal fluid, urine, and saliva. We employed published data sets representing physiological (resting vs. acute exercise) and pathologic states (early- vs. late-stage liver fibrosis, and differential subtypes of stroke) to determine differential tissue-enriched small RNAs. We also developed an online tool that provides information about exRNA sequences found in different biofluids and tissues. The data can be used to better understand the various types of small RNA sequences in different tissues as well as their potential release into biofluids, which should help in the validation or design of biomarker studies.

piRNAs and PIWI Proteins as Diagnostic and Prognostic Markers of Genitourinary Cancers

piRNAs (PIWI-interacting RNAs) are small non-coding RNAs capable of regulation of transposon and gene expression. piRNAs utilise multiple mechanisms to affect gene expression, which makes them potentially more powerful regulators than microRNAs. The mechanisms by which piRNAs regulate transposon and gene expression include DNA methylation, histone modifications, and mRNA degradation. Genitourinary cancers (GC) are a large group of neoplasms that differ by their incidence, clinical course, biology, and prognosis for patients. Regardless of the GC type, metastatic disease remains a key therapeutic challenge, largely affecting patients’ survival rates. Recent studies indicate that piRNAs could serve as potentially useful biomarkers allowing for early cancer detection and therapeutic interventions at the stage of non-advanced tumour, improving patient’s outcomes. Furthermore, studies in prostate cancer show that piRNAs contribute to cancer progression by affecting key oncogenic pathways such as PI3K/AKT. Here, we discuss recent findings on biogenesis, mechanisms of action and the role of piRNAs and the associated PIWI proteins in GC. We also present tools that may be useful for studies on the functioning of piRNAs in cancers.

Bioinformatic Tools for the Analysis and Prediction of ncRNA Interactions

Noncoding RNAs (ncRNAs) play prominent roles in the regulation of gene expression via their interactions with other biological molecules such as proteins and nucleic acids. Although much of our knowledge about how these ncRNAs operate in different biological processes has been obtained from experimental findings, computational biology can also clearly substantially boost this knowledge by suggesting possible novel interactions of these ncRNAs with other molecules. Computational predictions are thus used as an alternative source of new insights through a process of mutual enrichment because the information obtained through experiments continuously feeds through into computational methods. The results of these predictions in turn shed light on possible interactions that are subsequently validated experimentally. This review describes the latest advances in databases, bioinformatic tools, and new in silico strategies that allow the establishment or prediction of biological interactions of ncRNAs, particularly miRNAs and lncRNAs. The ncRNA species described in this work have a special emphasis on those found in humans, but information on ncRNA of other species is also included.

miRNATissueAtlas2: an update to the human miRNA tissue atlas

Small non-coding RNAs (sncRNAs) are pervasive regulators of physiological and pathological processes. We previously developed the human miRNA Tissue Atlas, detailing the expression of miRNAs across organs in the human body. Here, we present an updated resource containing sequencing data of 188 tissue samples comprising 21 organ types retrieved from six humans. Sampling the organs from the same bodies minimizes intra-individual variability and facilitates the making of a precise high-resolution body map of the non-coding transcriptome. The data allow shedding light on the organ- and organ system-specificity of piwi-interacting RNAs (piRNAs), transfer RNAs (tRNAs), microRNAs (miRNAs) and other non-coding RNAs. As use case of our resource, we describe the identification of highly specific ncRNAs in different organs. The update also contains 58 samples from six tissues of the Tabula Muris collection, allowing to check if the tissue specificity is evolutionary conserved between Homo sapiens and Mus musculus. The updated resource of 87 252 non-coding RNAs from nine non-coding RNA classes for all organs and organ systems is available online without any restrictions (https://www.ccb.uni-saarland.de/tissueatlas2).

DIANA-miTED: a microRNA tissue expression database

microRNAs (miRNAs) are short (∼23nt) single-stranded non-coding RNAs that act as potent post-transcriptional gene expression regulators. Information about miRNA expression and distribution across cell types and tissues is crucial to the understanding of their function and for their translational use as biomarkers or therapeutic targets. DIANA-miTED is the most comprehensive and systematic collection of miRNA expression values derived from the analysis of 15 183 raw human small RNA-Seq (sRNA-Seq) datasets from the Sequence Read Archive (SRA) and The Cancer Genome Atlas (TCGA). Metadata quality maximizes the utility of expression atlases, therefore we manually curated SRA and TCGA-derived information to deliver a comprehensive and standardized set, incorporating in total 199 tissues, 82 anatomical sublocations, 267 cell lines and 261 diseases. miTED offers rich instant visualizations of the expression and sample distributions of requested data across variables, as well as study-wide diagrams and graphs enabling efficient content exploration. Queries also generate links towards state-of-the-art miRNA functional resources, deeming miTED an ideal starting point for expression retrieval, exploration, comparison, and downstream analysis, without requiring bioinformatics support or expertise. DIANA-miTED is freely available at http://www.microrna.gr/mited.

Single-cell biology to decode the immune cellular composition of kidney inflammation

Single-cell biology is transforming the ability of researchers to understand cellular signaling and identity across medical and biological disciplines. Especially for immune-mediated diseases, a single-cell look at immune cell subtypes, signaling, and activity might yield fundamental insights into the disease etiology, mechanisms, and potential therapeutic interventions. In this review, we highlight recent advances in the field of single-cell RNA profiling and their application to understand renal function in health and disease. With a focus on the immune system, in particular on T cells, we propose some key directions of understanding renal inflammation using single-cell approaches. We detail the benefits and shortcomings of the various technological approaches outlined and give advice on potential pitfalls and challenges in experimental setup and computational analysis. Finally, we conclude with a brief outlook into a promising future for single-cell technologies to elucidate kidney function.

Roles of host small RNAs in the evolution and host tropism of coronaviruses

Human coronaviruses (CoVs) can cause respiratory infection epidemics that sometimes expand into globally relevant pandemics. All human CoVs have sister strains isolated from animal hosts and seem to have an animal origin, yet the process of host jumping is largely unknown. RNA interference (RNAi) is an ancient mechanism in many eukaryotes to defend against viral infections through the hybridization of host endogenous small RNAs (miRNAs) with target sites in invading RNAs. Here, we developed a method to identify potential RNAi-sensitive sites in the viral genome and discovered that human-adapted coronavirus strains had deleted some of their sites targeted by miRNAs in human lungs when compared to their close zoonic relatives. We further confirmed using a phylogenetic analysis that the loss of RNAi-sensitive target sites could be a major driver of the host-jumping process, and adaptive mutations that lead to the loss-of-target might be as simple as point mutation. Up-to-date genomic data of severe acute respiratory syndrome coronavirus 2 and Middle-East respiratory syndromes-CoV strains demonstrate that the stress from host miRNA milieus sustained even after their epidemics in humans. Thus, this study illustrates a new mechanism about coronavirus to explain its host-jumping process and provides a novel avenue for pathogenesis research, epidemiological modeling, and development of drugs and vaccines against coronavirus, taking into consideration these findings.

Two MicroRNAs Are Sufficient for Embryonic Patterning in C. elegans

MicroRNAs (miRNAs) are a class of post-transcriptional repressors with diverse roles in animal development and physiology [1]. The Microprocessor complex, composed of Drosha and Pasha/DGCR8, is necessary for the biogenesis of all canonical miRNAs and essential for the early stages of animal embryogenesis [2, 3, 4, 5, 6, 7, 8]. However, the cause for this requirement is largely unknown. Animals often express hundreds of miRNAs, and it remains unclear whether the Microprocessor is required to produce one or few essential miRNAs or many individually non-essential miRNAs. Additionally, both Drosha and Pasha/DGCR8 bind and cleave a variety of non-miRNA substrates [9, 10, 11, 12, 13, 14, 15], and it is unknown whether these activities account for the Microprocessor’s essential requirement. To distinguish between these possibilities, we developed a system in C. elegans to stringently deplete embryos of Microprocessor activity. Using a combination of auxin-inducible degradation (AID) and RNA interference (RNAi), we achieved Drosha and Pasha/DGCR8 depletion starting in the maternal germline, resulting in Microprocessor and miRNA-depleted embryos, which fail to undergo morphogenesis or form organs. Using a Microprocessor-bypass strategy, we show that this early embryonic arrest is rescued by the addition of just two miRNAs, one miR-35 and one miR-51 family member, resulting in morphologically normal larvae. Thus, just two out of ∼150 canonical miRNAs are sufficient for morphogenesis and organogenesis, and the processing of these miRNAs accounts for the essential requirement for Drosha and Pasha/DGCR8 during the early stages of C. elegans embryonic development.

Video Abstract

•

Depletion of Drosha and Pasha results in embryos that fail to undergo morphogenesis

•

The mirtron pathway enables expression of miRNAs in the absence of Drosha and Pasha

•

Two miRNAs are sufficient to rescue embryogenesis in the absence of Drosha and Pasha

•

miR-35 and miR-51 play an unexplored, likely conserved role in animal development

Explainable Deep Learning for Augmentation of Small RNA Expression Profiles

The lack of well-structured metadata annotations complicates the reusability and interpretation of the growing amount of publicly available RNA expression data. The machine learning-based prediction of metadata (data augmentation) can considerably improve the quality of expression data annotation. In this study, we systematically benchmark deep learning (DL) and random forest (RF)-based metadata augmentation of tissue, age, and sex using small RNA (sRNA) expression profiles. We use 4243 annotated sRNA-Seq samples from the sRNA expression atlas database to train and test the augmentation performance. In general, the DL machine learner outperforms the RF method in almost all tested cases. The average cross-validated prediction accuracy of the DL algorithm for tissues is 96.5%, for sex is 77%, and for age is 77.2%. The average tissue prediction accuracy for a completely new data set is 83.1% (DL) and 80.8% (RF). To understand which sRNAs influence DL predictions, we employ backpropagation-based feature importance scores using the DeepLIFT method, which enable us to obtain information on biological relevance of sRNAs.