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 usingmiRMaster 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.Saarland Universit

MirGeneDB 2.1: toward a complete sampling of all major animal phyla

B.F. is supported by the Tromso forskningsstiftelse (TFS) [20 SG BF `MIRevolution']; Strategic Research Area (SFO) program of the Swedish Research Council (to V.R.) through Stockholm University (to B.F., W.K., E.M.-S. and M.R.F.); M.R.F. is additionally supported by ERC [758397 `miRCell']; South-Eastern Norway Regional Health Authority support is acknowledged [2018014 to E.H.]; P.J. Chabot is supported by the Junior Scholars Program (Dartmouth College); V.O.'s research funding was awarded to Dr Mary J. O'Connell (Associate Professor) from the School of Life Sciences University of Nottingham; M.H. is supported by the Spanish Government [AGL2017-88702C2-2-R]; University of Granada [A-BIO-481-UGR18, FEDER 18]; K.J.P. has been supported by the National Science Foundation; NASA Ames; Dartmouth College.We describe an update of MirGeneDB, the manually curated microRNA gene database. Adhering to uniform and consistent criteria for microRNA annotation and nomenclature, we substantially expanded MirGeneDB with 30 additional species representing previously missing metazoan phyla such as sponges, jellyfish, rotifers and flatworms. MirGeneDB 2.1 now consists of 75 species spanning over ∼800 million years of animal evolution, and contains a total number of 16 670 microRNAs from 1549 families. Over 6000 microRNAs were added in this update using ∼550 datasets with ∼7.5 billion sequencing reads. By adding new phylogenetically important species, especially those relevant for the study of whole genome duplication events, and through updating evolutionary nodes of origin for many families and genes, we were able to substantially refine our nomenclature system. All changes are traceable in the specifically developed MirGeneDB version tracker. The performance of read-pages is improved and microRNA expression matrices for all tissues and species are now also downloadable. Altogether, this update represents a significant step toward a complete sampling of all major metazoan phyla, and a widely needed foundation for comparative microRNA genomics and transcriptomics studies. MirGeneDB 2.1 is part of RNAcentral and Elixir Norway, publicly and freely available at http://www.mirgenedb.org/.Tromso forskningsstiftelse (TFS) 20_SG_BFStrategic Research Area (SFO) program of the Swedish Research Council through Stockholm UniversityEuropean Research Council (ERC) European Commission 758397South-Eastern Norway Regional Health Authority 2018014Junior Scholars Program (Dartmouth College)School of Life Sciences University of NottinghamSpanish GovernmentEuropean Commission AGL2017-88702-C2-2-RUniversity of Granada A-BIO-481-UGR18 FEDER 18National Science Foundation (NSF)National Aeronautics & Space Administration (NASA)Dartmouth Colleg

Transcriptomic analysis of the tick midgut and salivary gland responses upon repeated blood-feeding on a vertebrate host

Ticks are blood-feeding arthropods that use the components of their salivary glands to counter the host’s hemostatic, inflammatory, and immune responses. The tick midgut also plays a crucial role in hematophagy. It is responsible for managing blood meals (storage and digestion) and protecting against host immunity and pathogen infections. Previous transcriptomic studies revealed the complexity of tick sialomes (salivary gland transcriptomes) and mialomes (midgut transcriptomes) which encode for protease inhibitors, lipocalins (histamine-binding proteins), disintegrins, enzymes, and several other tickspecific proteins. Several studies have demonstrated that mammalian hosts acquire tick resistance against repeated tick bites. Consequently, there is an urgent need to uncover how tick sialomes and mialomes respond to resistant hosts, as they may serve to develop novel tick control strategies and applications. Here, we mimicked natural repeated tick bites in a laboratory setting and analyzed gene expression dynamics in the salivary glands and midguts of adult female ticks. Rabbits were subjected to a primary (feeding on a naive host) and a secondary infestation of the same host (we re-exposed the hosts but to other ticks). We used single salivary glands and midguts dissected from individual siblings adult pathogen-free female Ixodes ricinus to reduce genetic variability between individual ticks. The comprehensive analysis of 88 obtained RNA-seq data sets allows us to provide high-quality annotated sialomes and mialomes from individual ticks. Comparisons between fed/ unfed, timepoints, and exposures yielded as many as 3000 putative differentially expressed genes (DEG). Interestingly, when classifying the exposure DEGs by means of a clustering approach we observed that the majority of these genes show increased expression at early feeding timepoints in the mid-gut of re-exposed ticks. The existence of clearly defined groups of genes with highly similar responses to re-exposure suggests the existence of molecular swiches. In silico functional analysis shows that these early feeding reexposure response genes form a dense interaction network at protein level being related to virtually all aspects of gene expression regulation and glycosylation. The processed data is available through an easy-to-use database-associated webpage (https://arn.ugr.es/IxoriDB/) that can serve as a valuable resource for tick research.Grant Agency of the Czech Republic 19-382 07247SERD Funds 384 CZ.02.1.01/0.0/0.0/16_019/0000759Programa Operativo FEDER de Andalucia A-BIO-481-UGR18European Union within ESIF in frame of Operational Programme Research, Development and Education CZ.02.2.69/0.0/0.0/20_079/001780

miEAA 2023: updates, new functional microRNA sets and improved enrichment visualizations

MicroRNAs (miRNAs) are small non-coding RNAs that play a critical role in regulating diverse biological processes. Extracting functional insights from a list of miRNAs is challenging, as each miRNA can potentially interact with hundreds of genes. To address this challenge, we developed miEAA, a flexible and comprehensive miRNA enrichment analysis tool based on direct and indirect miRNA annotation. The latest release of miEAA includes a data warehouse of 19 miRNA repositories, covering 10 different organisms and 139 399 functional categories. We have added information on the cellular context of miRNAs, isomiRs, and high-confidence miRNAs to improve the accuracy of the results. We have also improved the representation of aggregated results, including interactive Upset plots to aid users in understanding the interaction among enriched terms or categories. Finally, we demonstrate the functionality of miEAA in the context of ageing and highlight the importance of carefully considering the miRNA input list. MiEAA is free to use and publicly available at https://www.ccb.uni-saarland.de/mieaa/

Genome-Wide Analysis of microRNA Expression Profile in Roots and Leaves of Three Wheat Cultivars under Water and Drought Conditions

Wheat is one of the most important food sources on Earth. MicroRNAs (miRNAs) play important roles in wheat productivity. To identify wheat miRNAs as well as their expression profiles under drought condition, we constructed and sequenced small RNA (sRNA) libraries from the leaves and roots of three wheat cultivars (Kukri, RAC875 and Excalibur) under water and drought conditions. A total of 636 known miRNAs and 294 novel miRNAs were identified, of which 34 miRNAs were tissue- or cultivar-specific. Among these, 314 were significantly regulated under drought conditions. miRNAs that were drought-regulated in all cultivars displayed notably higher expression than those that responded in a cultivar-specific manner. Cultivar-specific drought response miRNAs were mainly detected in roots and showed significantly different drought regulations between cultivars. By using wheat degradome library, 6619 target genes were identified. Many target genes were strongly enriched for protein domains, such as MEKHLA, that play roles in drought response. Targeting analysis showed that drought-downregulated miRNAs targeted more genes than droughtupregulated miRNAs. Furthermore, such genes had more important functions. Additionally, the genes targeted by drought-downregulated miRNAs had multiple interactions with each other, while the genes targeted by drought-upregulated miRNAs had no interactions. Our data provide valuable information on wheat miRNA expression profiles and potential functions in different tissues, cultivars and drought conditions

liqDB: a small-RNAseq knowledge discovery database for liquid biopsy studies

MiRNAs are important regulators of gene expression and are frequently deregulated under pathologic conditions. They are highly stable in bodily fluids which makes them feasible candidates to become minimally invasive biomarkers. In fact, several studies already proposed circulating miRNA-based biomarkers for different types of neoplastic, cardiovascular and degenerative diseases. However, many of these studies rely on small RNA sequencing experiments that are based on different RNA extraction and processing protocols, rendering results incomparable. We generated liqDB, a database for liquid biopsy small RNA sequencing profiles that provides users with meaningful information to guide their small RNA liquid biopsy research and to overcome technical and conceptual problems. By means of a user-friendly web interface, miRNA expression profiles from 1607 manually annotated samples can be queried and explored at different levels. Result pages include downloadable expression matrices, differential expression analysis, most stably expressed miRNAs, cluster analysis and relevant visualizations by means of boxplots and heatmaps. We anticipate that liqDB will be a useful tool in liquid biopsy research as it provides a consistently annotated large compilation of experiments together with tools for reproducible analysis, comparison and hypothesis generation.European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie [grant agreement ELBA, 765492 to M.H. and D.K.L.]; Spanish Government [AGL2017-88702-C2-2-R to M.H.]; Instituto de Salud Carlos III, FEDER funds [PIE16/00045 to J.A.M.]; Chair ‘Doctors Galera-Requena in cancer stem cell research’ (to J.A.M.); Ministry of Education of Spain [FPU13/05662 to R.L. and IFI16/00041 to E.A.]. Funding for open access charge: Instituto de Salud Carlos III (FEDER funds) [PIE16/00045]

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

sRNAbench and sRNAtoolbox 2022 update: accurate miRNA and sncRNA profiling for model and non-model organisms

European Union [765492 to M.H.]; Spanish Government [AGL2017-88702-C2-2-R to M.H.]; Chair 'Doctors Galera-Requena in cancer stem cell research' (to J.A.M.); Tromsoforskningsstiftelse (TFS) [20 SG BF 'MIRevolution' to B.F.]; Stichting Cancer Center Amsterdam [CCA2021-9-77 to C.G]; TKI-Health Holland ['AQrate' project to C.G. and M.P.]. This publication is part of a project that has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 765492.The NCBI Sequence Read Archive currently hosts microRNA sequencing data for over 800 different species, evidencing the existence of a broad taxonomic distribution in the field of small RNA research. Simultaneously, the number of samples per miRNA-seq study continues to increase resulting in a vast amount of data that requires accurate, fast and user-friendly analysis methods. Since the previous release of sRNAtoolbox in 2019, 55 000 sRNAbench jobs have been submitted which has motivated many improvements in its usability and the scope of the underlying annotation database. With this update, users can upload an unlimited number of samples or import them from Google Drive, Dropbox or URLs. Micro- and small RNA profiling can now be carried out using high-confidence Metazoan and plant specific databases, MirGeneDB and PmiREN respectively, together with genome assemblies and libraries from 441 Ensembl species. The new results page includes straightforward sample annotation to allow downstream differential expression analysis with sRNAde. Unassigned reads can also be explored by means of a new tool that performsmapping to microbial references, which can reveal contamination events or biologically meaningful findings as we describe in the example. sRNAtoolbox is available at: https://arn.ugr.es/srnatoolbox/.European Commission 765492Spanish GovernmentEuropean Commission AGL2017-88702-C2-2-RChair 'Doctors Galera-Requena in cancer stem cell research'Stichting Cancer Center Amsterdam CCA2021-9-77Tromsoforskningsstiftelse (TFS) ['MIRevolution'] 20 SG BFTKI-Health Holland ['AQrate' project

Functional Enrichment Analysis of Regulatory Elements

This work has been partially supported by FEDER/Junta de Andalucia-Consejeria de Economia y Conocimiento/(grant CV20-36723), grant PID2020-119032RB-I00, MCIN/AEI/10.13039/501100011033 and FEDER/Junta de Andalucia-Consejeria de Transformacion Economica, Industria, Conocimiento y Universidades (Grant P20_00335).Statistical methods for enrichment analysis are important tools to extract biological information from omics experiments. Although these methods have been widely used for the analysis of gene and protein lists, the development of high-throughput technologies for regulatory elements demands dedicated statistical and bioinformatics tools. Here, we present a set of enrichment analysis methods for regulatory elements, including CpG sites, miRNAs, and transcription factors. Statistical significance is determined via a power weighting function for target genes and tested by theWallenius noncentral hypergeometric distribution model to avoid selection bias. These new methodologies have been applied to the analysis of a set of miRNAs associated with arrhythmia, showing the potential of this tool to extract biological information from a list of regulatory elements. These new methods are available in GeneCodis 4, a web tool able to perform singular and modular enrichment analysis that allows the integration of heterogeneous information.FEDER/Junta de Andalucia-Consejeria de Economia y Conocimiento CV20-36723MCIN/AEI PID2020-119032RB-I00FEDER/Junta de Andalucia-Consejeria de Transformacion Economica, Industria, Conocimiento y Universidades P20_0033

sRNAbench and sRNAtoolbox 2022 update: accurate miRNA and sncRNA profiling for model and non-model organisms

The NCBI Sequence Read Archive currently hosts microRNA sequencing data for over 800 different species, evidencing the existence of a broad taxonomic distribution in the field of small RNA research. Simultaneously, the number of samples per miRNA-seq study continues to increase resulting in a vast amount of data that requires accurate, fast and user-friendly analysis methods. Since the previous release of sRNAtoolbox in 2019, 55 000 sRNAbench jobs have been submitted which has motivated many improvements in its usability and the scope of the underlying annotation database. With this update, users can upload an unlimited number of samples or import them from Google Drive, Dropbox or URLs. Micro- and small RNA profiling can now be carried out using high-confidence Metazoan and plant specific databases, MirGeneDB and PmiREN respectively, together with genome assemblies and libraries from 441 Ensembl species. The new results page includes straightforward sample annotation to allow downstream differential expression analysis with sRNAde. Unassigned reads can also be explored by means of a new tool that performs mapping to microbial references, which can reveal contamination events or biologically meaningful findings as we describe in the example. sRNAtoolbox is available at: https://arn.ugr.es/srnatoolbox/</a