The RNA workbench 2.0: next generation RNA data analysis

Evolutionary Structure Conservation and Covariance Scores

Effective homology search for non-coding RNAs is frequently not possible via sequence similarity alone. Current methods leverage evolutionary information like structure conservation or covariance scores to identify homologs in organisms that are phylogenetically more distant. In this chapter, we introduce the theoretical background of evolutionary structure conservation and covariance score, and we show hands-on how current methods in the field are applied on example datasets.

IRF1 regulates self-renewal and stress responsiveness to support hematopoietic stem cell maintenance

Hematopoietic stem cells (HSCs) are tightly controlled to maintain a balance between blood cell production and self-renewal. While inflammation-related signaling is a critical regulator of HSC activity, the underlying mechanisms and the precise functions of specific factors under steady-state and stress conditions remain incompletely understood. We investigated the role of interferon regulatory factor 1 (IRF1), a transcription factor that is affected by multiple inflammatory stimuli, in HSC regulation. Our findings demonstrate that the loss of IRF1 from mouse HSCs significantly impairs self-renewal, increases stress-induced proliferation, and confers resistance to apoptosis. In addition, given the frequent abnormal expression of IRF1 in leukemia, we explored the potential of IRF1 expression level as a stratification marker for human acute myeloid leukemia. We show that IRF1-based stratification identifies distinct cancer-related signatures in patient subgroups. These findings establish IRF1 as a pivotal HSC controller and provide previously unknown insights into HSC regulation, with potential implications to IRF1 functions in the context of leukemia.

Interaction of Mesonivirus and Negevirus with arboviruses and the RNAi response in Culex tarsalis-derived cells

BACKGROUND

Mosquito-specific viruses (MSVs) comprise a variety of different virus families, some of which are known to interfere with infections of medically important arboviruses. Viruses belonging to the family Mesoniviridae or taxon Negevirus harbor several insect-specific viruses, including MSVs, which are known for their wide geographical distribution and extensive host ranges. Although these viruses are regularly identified in mosquitoes all over the world, their presence in mosquitoes in Germany had not yet been reported.

METHODS

A mix of three MSVs (Yichang virus [Mesoniviridae] and two negeviruses [Daeseongdong virus and Dezidougou virus]) in a sample that contained a pool of Coquillettidia richiardii mosquitoes collected in Germany was used to investigate the interaction of these viruses with different arboviruses in Culex-derived cells. In addition, small RNA sequencing and analysis of different mosquito-derived cells infected with this MSV mix were performed.

RESULTS

A strain of Yichang virus (Mesoniviridae) and two negeviruses (Daeseongdong virus and Dezidougou virus) were identified in the Cq. richiardii mosquitoes sampled in Germany, expanding current knowledge of their circulation in central Europe. Infection of mosquito-derived cells with these three viruses revealed that they are targeted by the small interfering RNA (siRNA) pathway. In Culex-derived cells, co-infection by these three viruses had varying effects on the representative arboviruses from different virus families (Togaviridae: Semliki forest virus [SFV]; Bunyavirales: Bunyamwera orthobunyavirus [BUNV]; or Flaviviridae: Usutu virus [USUV]). Specifically, persistent MSV co-infection inhibited BUNV infection, as well as USUV infection (but the latter only at specific time points). However, the impact on SFV infection was only noticeable at low multiplicity of infection (MOI 0.1) and at specific time points in combination with the infection status.

CONCLUSIONS

Taken together, these results are important findings that will lead to a better understanding of the complex interactions of MSVs, mosquitoes and arboviruses.

IRF1 regulates self-renewal and stress-responsiveness to support hematopoietic stem cell maintenance

Inflammatory mediators induce emergency myelopoiesis and cycling of adult hematopoietic stem cells (HSCs) through incompletely understood mechanisms. To suppress the unwanted effects of inflammation and preserve its beneficial outcomes, the mechanisms by which inflammation affects hematopoiesis need to be fully elucidated. Rather than focusing on specific inflammatory stimuli, we here investigated the role of transcription factor Interferon (IFN) regulatory factor 1 (IRF1), which receives input from several inflammatory signaling pathways. We identify IRF1 as a master HSC regulator. IRF1 loss impairs HSC self-renewal, increases stress-induced cell cycle activation, and confers apoptosis resistance. Transcriptomic analysis revealed an aged, inflammatory signature devoid of IFN signaling with reduced megakaryocytic/erythroid priming and antigen presentation in IRF1-deficient HSCs. Finally, we conducted IRF1-based AML patient stratification to identify groups with distinct proliferative, survival and differentiation features, overlapping with our murine HSC results. Our findings position IRF1 as a pivotal regulator of HSC preservation and stress-induced responses.

Antiviral RNAi Response in Culex quinquefasciatus-Derived HSU Cells

Culex spp. mosquitoes are important vectors of viruses, such as West Nile virus, Eastern equine encephalitis virus and Rift valley fever virus. However, their interactions with innate antiviral immunity, especially RNA interference (RNAi), are not well known. Most research on RNAi pathways in mosquitoes is focused on the tropical vector mosquito Aedes aegypti. Here, we investigated the production of arbovirus-specific small RNAs in Cx. quinquefasciatus-derived HSU cells. Furthermore, by silencing RNAi-related proteins, we investigated the antiviral role of these proteins for two different arboviruses: Semliki Forest virus (SFV) and Bunyamwera orthobunyavirus (BUNV). Our results showed an expansion of Ago2 and Piwi6 in Cx. quinquefasciatus compared to Ae. aegypti. While silencing Ago2a and Ago2b increased BUNV replication, only Ago2b showed antiviral activity against SFV. Our results suggest differences in the function of Cx. quinquefasciatus and Ae. aegypti RNAi proteins and highlight the virus-specific function of these proteins in Cx. quinquefasciatus.

High-Resolution Small RNAs Landscape Provides Insights into Alkane Adaptation in the Marine Alkane-Degrader Alcanivorax dieselolei B-5

Alkanes are widespread in the ocean, and Alcanivorax is one of the most ubiquitous alkane-degrading bacteria in the marine ecosystem. Small RNAs (sRNAs) are usually at the heart of regulatory pathways, but sRNA-mediated alkane metabolic adaptability still remains largely unknown due to the difficulties of identification. Here, differential RNA sequencing (dRNA-seq) modified with a size selection (~50-nt to 500-nt) strategy was used to generate high-resolution sRNAs profiling in the model species Alcanivorax dieselolei B-5 under alkane (n-hexadecane) and non-alkane (acetate) conditions. As a result, we identified 549 sRNA candidates at single-nucleotide resolution of 5'-ends, 63.4% of which are with transcription start sites (TSSs), and 36.6% of which are with processing sites (PSSs) at the 5'-ends. These sRNAs originate from almost any location in the genome, regardless of intragenic (65.8%), antisense (20.6%) and intergenic (6.2%) regions, and RNase E may function in the maturation of sRNAs. Most sRNAs locally distribute across the 15 reference genomes of Alcanivorax, and only 7.5% of sRNAs are broadly conserved in this genus. Expression responses to the alkane of several core conserved sRNAs, including 6S RNA, M1 RNA and tmRNA, indicate that they may participate in alkane metabolisms and result in more actively global transcription, RNA processing and stresses mitigation. Two novel CsrA-related sRNAs are identified, which may be involved in the translational activation of alkane metabolism-related genes by sequestering the global repressor CsrA. The relationships of sRNAs with the characterized genes of alkane sensing (ompS), chemotaxis (mcp, cheR, cheW2), transporting (ompT1, ompT2, ompT3) and hydroxylation (alkB1, alkB2, almA) were created based on the genome-wide predicted sRNA-mRNA interactions. Overall, the sRNA landscape lays the ground for uncovering cryptic regulations in critical marine bacterium, among which both the core and species-specific sRNAs are implicated in the alkane adaptive metabolisms.

Identification of Two Exosomal miRNAs in Circulating Blood of Cancer Patients by Using Integrative Transcriptome and Network Analysis

Exosomes carry molecules of great biological and clinical interest, such as miRNAs. The contents of exosomes vary between healthy controls and cancer patients. Therefore, miRNAs and other molecules transported in exosomes are considered a potential source of diagnostic and prognostic biomarkers in cancer. Many miRNAs have been detected in recent years. Consequently, a substantial amount of miRNA-related data comparing patients and healthy individuals is available, which contributes to a better understanding of the initiation, development, malignancy, and metastasis of cancer using non-invasive sampling procedures. However, a re-analysis of available ncRNA data is rare. This study used available data about miRNAs in exosomes comparing healthy individuals and cancer patients to identify possible global changes related to the presence of cancer. A robust transcriptomic analysis identified two common miRNAs (miR-495-3p and miR-543) deregulated in five cancer datasets. They had already been implicated in different cancers but not reported in exosomes circulating in blood. The study also examined their target genes and the implications of these genes for functional processes.

Micro-RNA-Regulated SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Expression and Cytokinin Accumulation Distinguish Early-Developing Male and Female Inflorescences in Oil Palm (Elaeis guineensis)

Sexual differentiation of inflorescences and flowers is important for reproduction and affects crop plant productivity. We report here on a molecular study of the process of sexual differentiation in the immature inflorescence of oil palm (Elaeis guineensis). This species is monoecious and exhibits gender diphasy, producing male and female inflorescences separately on the same plant in alternation. Three main approaches were used: small RNA-seq to characterise and study the expression of miRNA genes; RNA-seq to monitor mRNA accumulation patterns; hormone quantification to assess the role of cytokinins and auxins in inflorescence differentiation. Our study allowed the characterisation of 30 previously unreported palm MIRNA genes. In differential gene and miRNA expression studies, we identified a number of key developmental genes and miRNA-mRNA target modules previously described in relation to their developmental regulatory role in the cereal panicle, notably the miR156/529/535-SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE (SPL) gene regulatory module. Gene enrichment analysis highlighted the importance of hormone-related genes, and this observation was corroborated by the detection of much higher levels of cytokinins in the female inflorescence. Our data illustrate the importance of branching regulation within the developmental window studied, during which the female inflorescence, unlike its male counterpart, produces flower clusters on new successive axes by sympodial growth.

Bioinformatics and Machine Learning Approaches to Understand the Regulation of Mobile Genetic Elements

Simple Summary

Transposable elements (TEs) are DNA sequences that are, or were, able to move (transpose) within the genome of a single cell. They were first discovered by Barbara McClintock while working on maize, and they make up a large fraction of the genome. Transpositions can result in mutations and they can alter the genome size. Cells regulate the activity of TEs using a variety of mechanisms, such as chemical modifications of DNA and small RNAs. Machine learning (ML) is an interdisciplinary subject that studies computer algorithms that can improve through experience and by the use of data. ML has been successfully applied to a variety of problems in bioinformatics and has exhibited favorable precision and speed. Here, we provide a systematic and guided review on the ML and bioinformatic methods and tools that are used for the analysis of the regulation of TEs.

Abstract

Transposable elements (TEs, or mobile genetic elements, MGEs) are ubiquitous genetic elements that make up a substantial proportion of the genome of many species. The recent growing interest in understanding the evolution and function of TEs has revealed that TEs play a dual role in genome evolution, development, disease, and drug resistance. Cells regulate TE expression against uncontrolled activity that can lead to developmental defects and disease, using multiple strategies, such as DNA chemical modification, small RNA (sRNA) silencing, chromatin modification, as well as sequence-specific repressors. Advancements in bioinformatics and machine learning approaches are increasingly contributing to the analysis of the regulation mechanisms. A plethora of tools and machine learning approaches have been developed for prediction, annotation, and expression profiling of sRNAs, for methylation analysis of TEs, as well as for genome-wide methylation analysis through bisulfite sequencing data. In this review, we provide a guided overview of the bioinformatic and machine learning state of the art of fields closely associated with TE regulation and function.

ChiRA: an integrated framework for chimeric read analysis from RNA-RNA interactome and RNA structurome data

BACKGROUND

With the advances in next-generation sequencing technologies, it is possible to determine RNA-RNA interaction and RNA structure predictions on a genome-wide level. The reads from these experiments usually are chimeric, with each arm generated from one of the interaction partners. Owing to short read lengths, often these sequenced arms ambiguously map to multiple locations. Thus, inferring the origin of these can be quite complicated. Here we present ChiRA, a generic framework for sensitive annotation of these chimeric reads, which in turn can be used to predict the sequenced hybrids.

RESULTS

Grouping reference loci on the basis of aligned common reads and quantification improved the handling of the multi-mapped reads in contrast to common strategies such as the selection of the longest hit or a random choice among all hits. On benchmark data ChiRA improved the number of correct alignments to the reference up to 3-fold. It is shown that the genes that belong to the common read loci share the same protein families or similar pathways. In published data, ChiRA could detect 3 times more new interactions compared to existing approaches. In addition, ChiRAViz can be used to visualize and filter large chimeric datasets intuitively.

CONCLUSION

ChiRA tool suite provides a complete analysis and visualization framework along with ready-to-use Galaxy workflows and tutorials for RNA-RNA interactome and structurome datasets. Common read loci built by ChiRA can rescue multi-mapped reads on paralogous genes without requiring any information on gene relations. We showed that ChiRA is sensitive in detecting new RNA-RNA interactions from published RNA-RNA interactome datasets.

GLASSgo in Galaxy: high-throughput, reproducible and easy-to-integrate prediction of sRNA homologs

MOTIVATION

The correct prediction of bacterial sRNA homologs is a prerequisite for many downstream analyses based on comparative genomics, but it is frequently challenging due to the short length and distinct heterogeneity of such homologs. GLobal Automatic Small RNA Search go (GLASSgo) is an efficient tool for the prediction of sRNA homologs from a single input query. To make the algorithm available to a broader community, we offer a Docker container along with a free-access web service. For non-computer scientists, the web service provides a user-friendly interface. However, capabilities were lacking so far for batch processing, version control and direct interaction with compatible software applications as a workflow management system can provide.

RESULTS

Here, we present GLASSgo 1.5.2, an updated version that is fully incorporated into the workflow management system Galaxy. The improved version contains a new feature for extracting the upstream regions, allowing the search for conserved promoter elements. Additionally, it supports the use of accession numbers instead of the outdated GI numbers, which widens the applicability of the tool.

AVAILABILITY AND IMPLEMENTATION

GLASSgo is available at https://github.com/lotts/GLASSgo/ under the MIT license and is accompanied by instruction and application data. Furthermore, it can be installed into any Galaxy instance using the Galaxy ToolShed.

Dataset for de novo transcriptome assembly of the African bullfrog Pyxicephalus adspersus

In this article, we report the first de novo transcriptome assembly of the African bullfrog Pyxicephalus adspersus. In this data, 75,320,390 raw reads were acquired from African bullfrog mRNA using Illumina paired-end sequencing platform. De novo assembly resulted in a total of 136,958 unigenes. In the obtained unigenes, 30,039 open reading frames (ORFs) were detected. This dataset provides basic information for molecular level analysis of this species, which undergoes a state of dormancy under dry conditions at ordinary temperatures called estivation.