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

In search of epigenetic hallmarks of different tissues: an integrative omics study of horse liver, lung, and heart

DNA methylation and microRNA (miRNA) expression are epigenetic mechanisms essential for regulating tissue-specific gene expression and metabolic processes. However, high-resolution transcriptome, methylome, or miRNAome data is only available for a few model organisms and selected tissues. Up to date, only a few studies have reported on gene expression, DNA methylation, or miRNA expression in adult equine tissues at the genome-wide level. In the present study, we used RNA-Seq, miRNA-seq, and reduced representation bisulfite sequencing (RRBS) data from the heart, lung, and liver tissues of healthy cold-blooded horses to identify differentially expressed genes (DEGs), differentially expressed miRNA (DE miRNA) and differentially methylated sites (DMSs) between three types of horse tissues. Additionally, based on integrative omics analysis, we described the observed interactions of epigenetic mechanisms with tissue-specific gene expression alterations. The obtained data allowed identification from 4067 to 6143 DMSs, 9733 to 11,263 mRNAs, and 155 to 185 microRNAs, differentially expressed between various tissues. We pointed out specific genes whose expression level displayed a negative correlation with the level of CpG methylation and miRNA expression and revealed biological processes that they enrich. Furthermore, we confirmed and validated the accuracy of the Next-Generation Sequencing (NGS) results with bisulfite sequencing PCR (BSP) and quantitative PCR (qPCR). This comprehensive analysis forms a strong foundation for exploring the epigenetic mechanisms involved in tissue differentiation, especially the growth and development of the equine heart, lungs, and liver.

miRglmm: a generalized linear mixed model of isomiR-level counts improves estimation of miRNA-level differential expression and uncovers variable differential expression between isomiRs

MicroRNA-seq data is produced by aligning small RNA sequencing reads of different miRNA transcript isoforms, called isomiRs, to known microRNAs. Aggregation to microRNA-level counts discards information and violates core assumptions of differential expression (DE) methods developed for mRNA-seq data. We establish miRglmm, a DE method for microRNA-seq data, that uses a generalized linear mixed model of isomiR-level counts, facilitating detection of miRNA with differential expression or differential isomiR usage. We demonstrate that miRglmm outperforms current DE methods in estimating DE for miRNA, whether or not there is significant isomiR variability, and simultaneously provides estimates of isomiR-level DE.

Analysis of Amblyomma americanum microRNAs in response to Ehrlichia chaffeensis infection and their potential role in vectorial capacity

Background

MicroRNAs (miRNAs) represent a subset of small noncoding RNAs and carry tremendous potential for regulating gene expression at the post-transcriptional level. They play pivotal roles in distinct cellular mechanisms including inhibition of bacterial, parasitic, and viral infections via immune response pathways. Intriguingly, pathogens have developed strategies to manipulate the host's miRNA profile, fostering environments conducive to successful infection. Therefore, changes in an arthropod host's miRNA profile in response to pathogen invasion could be critical in understanding host-pathogen dynamics. Additionally, this area of study could provide insights into discovering new targets for disease control and prevention. The main objective of the present study is to investigate the functional role of differentially expressed miRNAs upon Ehrlichia chaffeensis, a tick-borne pathogen, infection in tick vector, Amblyomma americanum.

Methods

Small RNA libraries from uninfected and E. chaffeensis-infected Am. americanum midgut and salivary gland tissues were prepared using the Illumina Truseq kit. Small RNA sequencing data was analyzed using miRDeep2 and sRNAtoolbox to identify novel and known miRNAs. The differentially expressed miRNAs were validated using a quantitative PCR assay. Furthermore, a miRNA inhibitor approach was used to determine the functional role of selected miRNA candidates.

Results

The sequencing of small RNA libraries generated >147 million raw reads in all four libraries and identified a total of >250 miRNAs across the four libraries. We identified 23 and 14 differentially expressed miRNAs in salivary glands, and midgut tissues infected with E. chaffeensis, respectively. Three differentially expressed miRNAs (miR-87, miR-750, and miR-275) were further characterized to determine their roles in pathogen infection. Inhibition of target miRNAs significantly decreased the E. chaffeensis load in tick tissues, which warrants more in-depth mechanistic studies.

Conclusions

The current study identified known and novel miRNAs and suggests that interfering with these miRNAs may impact the vectorial capacity of ticks to harbor Ehrlichia. This study identified several new miRNAs for future analysis of their functions in tick biology and tick-pathogen interaction studies.

SingmiR: a single-cell miRNA alignment and analysis tool

Single-cell RNA sequencing (RNA-seq) has revolutionized our understanding of cell biology, developmental and pathophysiological molecular processes, paving the way toward novel diagnostic and therapeutic approaches. However, most of the gene regulatory processes on the single-cell level are still unknown, including post-transcriptional control conferred by microRNAs (miRNAs). Like the established single-cell gene expression analysis, advanced computational expertise is required to comprehensively process newly emerging single-cell miRNA-seq datasets. A web server providing a workflow tailored for single-cell miRNA-seq data with a self-explanatory interface is currently not available. Here, we present SingmiR, enabling the rapid (pre-)processing and quantification of human miRNAs from noncoding single-cell samples. It performs read trimming for different library preparation protocols, generates automated quality control reports and provides feature-normalized count files. Numerous standard and advanced analyses such as dimension reduction, clustered feature heatmaps, sample correlation heatmaps and differential expression statistics are implemented. We aim to speed up the prototyping pipeline for biologists developing single-cell miRNA-seq protocols on small to medium-sized datasets. SingmiR is freely available to all users without the need for a login at https://www.ccb.uni-saarland.de/singmir.

The intricacies of isomiRs: from classification to clinical relevance

MicroRNAs (miRNAs) and isoforms of their archetype, called isomiRs, regulate gene expression via complementary base-pair binding to messenger RNAs (mRNAs). The partially evolutionarily conserved isomiR sequence variations are differentially expressed among tissues, populations, and genders, and between healthy and diseased states. Aiming towards the clinical use of isomiRs as diagnostic biomarkers and for therapeutic purposes, several challenges need to be addressed, including (i) clarification of isomiR definition, (ii) improved annotation in databases with new standardization (such as the mirGFF3 format), and (iii) improved methods of isomiR detection, functional verification, and in silico analysis. In this review we discuss the respective challenges, and highlight the opportunities for clinical use of isomiRs, especially in the light of increasing amounts of next-generation sequencing (NGS) data.

Early life interventions metformin and trodusquemine metabolically reprogram the developing mouse liver through transcriptomic alterations

Recent studies have demonstrated the remarkable potential of early life intervention strategies at influencing the course of postnatal development, thereby offering exciting possibilities for enhancing longevity and improving overall health. Metformin (MF), an FDA-approved medication for type II diabetes mellitus, has recently gained attention for its promising anti-aging properties, acting as a calorie restriction mimetic, and delaying precocious puberty. Additionally, trodusquemine (MSI-1436), an investigational drug, has been shown to combat obesity and metabolic disorders by inhibiting the enzyme protein tyrosine phosphatase 1b (Ptp1b), consequently reducing hepatic lipogenesis and counteracting insulin and leptin resistance. In this study, we aimed to further explore the effects of these compounds on young, developing mice to uncover biomolecular signatures that are central to liver metabolic processes. We found that MSI-1436 more potently alters mRNA and miRNA expression in the liver compared with MF, with bioinformatic analysis suggesting that cohorts of differentially expressed miRNAs inhibit the action of phosphoinositide 3-kinase (Pi3k), protein kinase B (Akt), and mammalian target of rapamycin (Mtor) to regulate the downstream processes of de novo lipogenesis, fatty acid oxidation, very-low-density lipoprotein transport, and cholesterol biosynthesis and efflux. In summary, our study demonstrates that administering these compounds during the postnatal window metabolically reprograms the liver through induction of potent epigenetic changes in the transcriptome, potentially forestalling the onset of age-related diseases and enhancing longevity. Future studies are necessary to determine the impacts on lifespan and overall quality of life.

Full-length tRNAs lacking a functional CCA tail are selectively sorted into the lumen of extracellular vesicles

Small extracellular vesicles (sEVs) are heterogenous lipid membrane particles typically less than 200 nm in size and secreted by most cell types either constitutively or upon activation signals. sEVs isolated from biofluids contain RNAs, including small non-coding RNAs (ncRNAs), that can be either encapsulated within the EV lumen or bound to the EV surface. EV-associated microRNAs (miRNAs) are, despite a relatively low abundance, extensively investigated for their selective incorporation and their role in cell-cell communication. In contrast, the sorting of highly-structured ncRNA species is understudied, mainly due to technical limitations of traditional small RNA sequencing protocols. Here, we adapted ALL-tRNAseq to profile the relative abundance of highly structured and potentially methylated small ncRNA species, including transfer RNAs (tRNAs), small nucleolar RNAs (snoRNAs), and Y RNAs in bulk EV preparations. We determined that full-length tRNAs, typically 75 to 90 nucleotides in length, were the dominant small ncRNA species (>60% of all reads in the 18-120 nucleotides size-range) in all cell culture-derived EVs, as well as in human plasma-derived EV samples, vastly outnumbering 21 nucleotides-long miRNAs. Nearly all EV-associated tRNAs were protected from external RNAse treatment, indicating a location within the EV lumen. Strikingly, the vast majority of luminal-sorted, full-length, nucleobase modification-containing EV-tRNA sequences, harbored a dysfunctional 3' CCA tail, 1 to 3 nucleotides truncated, rendering them incompetent for amino acid loading. In contrast, in non-EV associated extracellular particle fractions (NVEPs), tRNAs appeared almost exclusively fragmented or 'nicked' into tRNA-derived small RNAs (tsRNAs) with lengths between 18 to 35 nucleotides. We propose that in mammalian cells, tRNAs that lack a functional 3' CCA tail are selectively sorted into EVs and shuttled out of the producing cell, offering a new perspective into the physiological role of secreted EVs and luminal cargo-selection.

Analysis of Amblyomma americanum microRNAs in response to Ehrlichia chaffeensis infection and their potential role in vectorial capacity

Background

MicroRNAs (miRNAs) represent a subset of small noncoding RNAs and carry tremendous potential for regulating gene expression at the post-transcriptional level. They play pivotal roles in distinct cellular mechanisms including inhibition of bacterial, parasitic, and viral infections via immune response pathways. Intriguingly, pathogens have developed strategies to manipulate the host's miRNA profile, fostering environments conducive to successful infection. Therefore, changes in an arthropod host's miRNA profile in response to pathogen invasion could be critical in understanding host-pathogen dynamics. Additionally, this area of study could provide insights into discovering new targets for disease control and prevention. The main objective of the present study is to investigate the functional role of differentially expressed miRNAs upon Ehrlichia chaffeensis, a tick-borne pathogen, infection in tick vector, Amblyomma americanum.

Methods

Small RNA libraries from uninfected and E. chaffeensis-infected Am. americanum midgut and salivary gland tissues were prepared using the Illumina Truseq kit. Small RNA sequencing data was analyzed using miRDeep2 and sRNAtoolbox to identify novel and known miRNAs. The differentially expressed miRNAs were validated using a quantitative PCR assay. Furthermore, a miRNA inhibitor approach was used to determine the functional role of selected miRNA candidates.

Results

The sequencing of small RNA libraries generated >147 million raw reads in all four libraries and identified a total of >250 miRNAs across the four libraries. We identified 23 and 14 differentially expressed miRNAs in salivary glands, and midgut tissues infected with E. chaffeensis, respectively. Three differentially expressed miRNAs (miR-87, miR-750, and miR-275) were further characterized to determine their roles in pathogen infection. Inhibition of target miRNAs significantly decreased the E. chaffeensis load in tick tissues, which warrants more in-depth mechanistic studies.

Conclusions

The current study identified known and novel miRNAs and suggests that interfering with these miRNAs may impact the vectorial capacity of ticks to harbor Ehrlichia. This study identified several new miRNAs for future analysis of their functions in tick biology and tick-pathogen interaction studies.

Spitting in the wind?-The challenges of RNA sequencing for biomarker discovery from saliva

Forensic trace contextualization, i.e., assessing information beyond who deposited a biological stain, has become an issue of great and steadily growing importance in forensic genetic casework and research. The human transcriptome encodes a wide variety of information and thus has received increasing interest for the identification of biomarkers for different aspects of forensic trace contextualization over the past years. Massively parallel sequencing of reverse-transcribed RNA ("RNA sequencing") has emerged as the gold standard technology to characterize the transcriptome in its entirety and identify RNA markers showing significant expression differences not only between different forensically relevant body fluids but also within a single body fluid between forensically relevant conditions of interest. Here, we analyze the quality and composition of four RNA sequencing datasets (whole transcriptome as well as miRNA sequencing) from two different research projects (the RNAgE project and the TrACES project), aiming at identifying contextualizing forensic biomarker from the forensically relevant body fluid saliva. We describe and characterize challenges of RNA sequencing of saliva samples arising from the presence of oral bacteria, the heterogeneity of sample composition, and the confounding factor of degradation. Based on these observations, we formulate recommendations that might help to improve RNA biomarker discovery from the challenging but forensically relevant body fluid saliva.

Characterisation of the in-vivo miRNA landscape in Drosophila ribonuclease mutants reveals Pacman-mediated regulation of the highly conserved let-7 cluster during apoptotic processes

The control of gene expression is a fundamental process essential for correct development and to maintain homeostasis. Many post-transcriptional mechanisms exist to maintain the correct levels of each RNA transcript within the cell. Controlled and targeted cytoplasmic RNA degradation is one such mechanism with the 5'-3' exoribonuclease Pacman (XRN1) and the 3'-5' exoribonuclease Dis3L2 playing crucial roles. Loss of function mutations in either Pacman or Dis3L2 have been demonstrated to result in distinct phenotypes, and both have been implicated in human disease. One mechanism by which gene expression is controlled is through the function of miRNAs which have been shown to be crucial for the control of almost all cellular processes. Although the biogenesis and mechanisms of action of miRNAs have been comprehensively studied, the mechanisms regulating their own turnover are not well understood. Here we characterise the miRNA landscape in a natural developing tissue, the Drosophila melanogaster wing imaginal disc, and assess the importance of Pacman and Dis3L2 on the abundance of miRNAs. We reveal a complex landscape of miRNA expression and show that whilst a null mutation in dis3L2 has a minimal effect on the miRNA expression profile, loss of Pacman has a profound effect with a third of all detected miRNAs demonstrating Pacman sensitivity. We also reveal a role for Pacman in regulating the highly conserved let-7 cluster (containing miR-100, let-7 and miR-125) and present a genetic model outlining a positive feedback loop regulated by Pacman which enhances our understanding of the apoptotic phenotype observed in Pacman mutants.

Computational tools supporting known miRNA identification

The study of small RNAs is a field that is expanding quickly. Other functional short RNA molecules other than microRNAs, and gene expression regulators, have been found in animals and plants. MicroRNAs play a significant role in host-microbe interactions, and parasite microRNAs may affect the host's innate immunity. Furthermore, short RNAs are intriguing non-invasive biomarker possibilities because they can be found in physiological fluids. These trends suggest that for many researchers, quick and simple techniques for expression profiling and subsequent downstream analysis of miRNA-seq data are crucial. We selected sRNAtoolbox to make integrated sRNA research easier. Each tool can be used separately or to explore and analyze sRNAbench results in further depth. A special focus was placed on the tools' usability. We review available miRNA research tools to have an overview of the evaluation of the tools. Mainly we evaluate the tool sRNAtoolbox.

sRNAflow: A Tool for the Analysis of Small RNA-Seq Data

The analysis of small RNA sequencing data across a range of biofluids is a significant research area, given the diversity of RNA types that hold potential diagnostic, prognostic, and predictive value. The intricate task of segregating the complex mixture of small RNAs from both human and other species, including bacteria, fungi, and viruses, poses one of the most formidable challenges in the analysis of small RNA sequencing data, currently lacking satisfactory solutions. This study introduces sRNAflow, a user-friendly bioinformatic tool with a web interface designed for the analysis of small RNAs obtained from biological fluids. Tailored to the unique requirements of such samples, the proposed pipeline addresses various challenges, including filtering potential RNAs from reagents and environment, classifying small RNA types, managing small RNA annotation overlap, conducting differential expression assays, analysing isomiRs, and presenting an approach to identify the sources of small RNAs within samples. sRNAflow also encompasses an alternative alignment-free analysis of RNA-seq data, featuring clustering and initial RNA source identification using BLAST. This comprehensive approach facilitates meaningful comparisons of results between different analytical methods.

IsoSeek for unbiased and UMI-informed sequencing of miRNAs from low input samples at single-nucleotide resolution

Besides canonical microRNAs (miRNAs), sequence-based variations called isomiRs have biological relevance and diagnostic potential; however, accurate calling of these post-transcriptional modifications is challenging, especially for low input samples. Here, we present IsoSeek, a sequencing protocol that reduces ligation and PCR amplification bias and improves the accuracy of miRNA detection in low input samples. We describe steps for using randomized adapters combined with unique molecular identifiers (UMI), library quantification, and sequencing, followed by detailed procedures for data processing and analysis. For complete details on the use and execution of this protocol, please refer to C. Gómez-Martín et al. (2023)1 and Van Eijndhoven et al. (2021).2.

Benchmarking of bioinformatics tools for NGS-based microRNA profiling with RT-qPCR method

MicroRNAs are vital gene expression regulators, extensively studied worldwide. The large-scale characterization of miRNAomes is possible using next-generation sequencing (NGS). This technology offers great opportunities, but these cannot be fully exploited without proper and comprehensive bioinformatics analysis. This may be achieved by the use of reliable dedicated software; however, different programs may generate divergent results, leading to additional discrepancies. Thus, the aim of this study was to compare three bioinformatic algorithms dedicated to NGS-based microRNA profiling and validate them using an alternative method, namely RT-qPCR. The comparison analysis revealed differences in the number and sets of identified miRNAs. The qPCR confirmed the expression of the investigated microRNAs. The correlation analysis of NGS and qPCR measurements showed strong and significant coefficients for a subset of the tested miRNAs, including those detected by all three algorithms. Single miRNA variants (isomiRs) showed different levels of correlation with the qPCR data. The obtained results revealed the good performance of all tested programs, despite the observed differences. Moreover, they implied that some specific miRNAs may be differentially estimated using NGS technology and the qPCR method, regardless of the used bioinformatics software. These discrepancies may stem from many factors, including the composition of the isomiR profile, their abundance, length, and investigated species. In conclusion, in this study, we shed light on the bioinformatics aspects of miRNAome profiling, elucidating its complexity and pinpointing potential features influencing validation. Thus, qPCR validation results should be open to interpretation when not fully concordant with NGS results until further, additional analyses are conducted.

Unraveling the Potential of miRNAs from CSCs as an Emerging Clinical Tool for Breast Cancer Diagnosis and Prognosis

Breast cancer (BC) is the most diagnosed cancer in women and the second most common cancer globally. Significant advances in BC research have led to improved early detection and effective therapies. One of the key challenges in BC is the presence of BC stem cells (BCSCs). This small subpopulation within the tumor possesses unique characteristics, including tumor-initiating capabilities, contributes to treatment resistance, and plays a role in cancer recurrence and metastasis. In recent years, microRNAs (miRNAs) have emerged as potential regulators of BCSCs, which can modulate gene expression and influence cellular processes like BCSCs' self-renewal, differentiation, and tumor-promoting pathways. Understanding the miRNA signatures of BCSCs holds great promise for improving BC diagnosis and prognosis. By targeting BCSCs and their associated miRNAs, researchers aim to develop more effective and personalized treatment strategies that may offer better outcomes for BC patients, minimizing tumor recurrence and metastasis. In conclusion, the investigation of miRNAs as regulators of BCSCs opens new directions for advancing BC research through the use of bioinformatics and the development of innovative therapeutic approaches. This review summarizes the most recent and innovative studies and clinical trials on the role of BCSCs miRNAs as potential tools for early diagnosis, prognosis, and resistance.

HSV-1 miRNAs are post-transcriptionally edited in latently infected human ganglia

IMPORTANCE

Herpes simplex virus 1 is an important human pathogen that has been intensively studied for many decades. Nevertheless, the molecular mechanisms regulating its establishment, maintenance, and reactivation from latency are poorly understood. Here, we show that HSV-1-encoded miR-H2 is post-transcriptionally edited in latently infected human tissues. Hyperediting of viral miRNAs increases the targeting potential of these miRNAs and may play an important role in regulating latency. We show that the edited miR-H2 can target ICP4, an essential viral protein. Interestingly, we found no evidence of hyperediting of its homolog, miR-H2, which is expressed by the closely related virus HSV-2. The discovery of post-translational modifications of viral miRNA in the latency phase suggests that these processes may also be important for other non-coding viral RNA in the latency phase, including the intron LAT, which in turn may be crucial for understanding the biology of this virus.

Starvation induces changes in abundance and small RNA cargo of extracellular vesicles released from Plasmodium falciparum infected red blood cells

The lethal malaria parasite Plasmodium falciparum needs to constantly respond and adapt to changes within the human host in order to survive and transmit. One such change is composed of nutritional limitation, which is augmented with increased parasite loads and intimately linked to severe disease development. Extracellular vesicles released from infected red blood cells have been proposed as important mediators of disease pathogenesis and intercellular communication but whether important for the parasite response to nutritional availability is unknown. Therefore, we investigated the abundance and small RNA cargo of extracellular vesicles released upon short-term nutritional starvation of P. falciparum in vitro cultures. We show that primarily ring-stage parasite cultures respond to glucose and amino acid deprivation with an increased release of extracellular vesicles. Small RNA sequencing of these extracellular vesicles further revealed human miRNAs and parasitic tRNA fragments as the main constituent biotypes. Short-term starvations led to alterations in the transcriptomic profile, most notably in terms of the over-represented biotypes. These data suggest a potential role for extracellular vesicles released from P. falciparum infected red blood cells in the response to nutritional perturbations, their potential as prognostic biomarkers and point towards an evolutionary conserved role among protozoan parasites.

Insights into the microRNA landscape of Rhodnius prolixus, a vector of Chagas disease

The growing interest in microRNAs (miRNAs) over recent years has led to their characterization in numerous organisms. However, there is currently a lack of data available on miRNAs from triatomine bugs (Reduviidae: Triatominae), which are the vectors of the protozoan parasite Trypanosoma cruzi, the causative agent of Chagas disease. A comprehensive understanding of the molecular biology of vectors provides new insights into insect-host interactions and insect control approaches, which are key methods to prevent disease incidence in endemic areas. In this work, we describe the miRNome profiles from gut, hemolymph, and salivary gland tissues of the Rhodnius prolixus triatomine. Small RNA sequencing data revealed abundant expression of miRNAs, along with tRNA- and rRNA-derived fragments. Fifty-two mature miRNAs, previously reported in Ecdysozoa, were identified, including 39 ubiquitously expressed in the three tissues. Additionally, 112, 73, and 78 novel miRNAs were predicted in the gut, hemolymph, and salivary glands, respectively. In silico prediction showed that the top eight most highly expressed miRNAs from salivary glands potentially target human blood-expressed genes, suggesting that R. prolixus may modulate the host's gene expression at the bite site. This study provides the first characterization of miRNAs in a Triatominae species, shedding light on the role of these crucial regulatory molecules.

Accurate microRNA annotation of animal genomes using trained covariance models of curated microRNA complements in MirMachine

The annotation of microRNAs depends on the availability of transcriptomics data and expert knowledge. This has led to a gap between the availability of novel genomes and high-quality microRNA complements. Using >16,000 microRNAs from the manually curated microRNA gene database MirGeneDB, we generated trained covariance models for all conserved microRNA families. These models are available in our tool MirMachine, which annotates conserved microRNAs within genomes. We successfully applied MirMachine to a range of animal species, including those with large genomes and genome duplications and extinct species, where small RNA sequencing is hard to achieve. We further describe a microRNA score of expected microRNAs that can be used to assess the completeness of genome assemblies. MirMachine closes a long-persisting gap in the microRNA field by facilitating automated genome annotation pipelines and deeper studies into the evolution of genome regulation, even in extinct organisms.

Integrative Transcriptomic Profiling of the Wilms Tumor

Our study aimed to identify relevant transcriptomic biomarkers for the Wilms tumor, the most common pediatric kidney cancer, independent of the histological type and stage. Using next-generation sequencing, we analyzed the miRNA profiles of 74 kidney samples, which were divided into two independent groups: fresh frozen tissue and formalin-fixed paraffin-embedded tissue samples. Subsequent mRNA expression profiling and pathway analysis were performed to establish the interplay and potential involvement of miRNAs and mRNA in the Wilms tumor. Comparative analysis, irrespective of post-dissection tissue processing, revealed 41 differentially expressed miRNAs, with 27 miRNAs having decreased expression and 14 miRNAs having increased expression in the Wilms tumor tissue compared to healthy kidney tissue. Among global mRNA transcriptomic profile differences, cross-sectional analysis suggested a limited list of genes potentially regulated by differentially expressed miRNAs in the Wilms tumor. This study identified the comprehensive miRNA and mRNA profile of the Wilms tumor using next-generation sequencing and bioinformatics approach, providing better insights into the pathogenesis of the Wilms tumor. The identified Wilms tumor miRNAs have potential as biomarkers for the diagnosis and treatment of the Wilms tumor, regardless of histological subtype and disease stage.

Small RNA sequencing of field Culex mosquitoes identifies patterns of viral infection and the mosquito immune response

Mosquito-borne disease remains a significant burden on global health. In the United States, the major threat posed by mosquitoes is transmission of arboviruses, including West Nile virus by mosquitoes of the Culex genus. Virus metagenomic analysis of mosquito small RNA using deep sequencing and advanced bioinformatic tools enables the rapid detection of viruses and other infecting organisms, both pathogenic and non-pathogenic to humans, without any precedent knowledge. In this study, we sequenced small RNA samples from over 60 pools of Culex mosquitoes from two major areas of Southern California from 2017 to 2019 to elucidate the virome and immune responses of Culex. Our results demonstrated that small RNAs not only allowed the detection of viruses but also revealed distinct patterns of viral infection based on location, Culex species, and time. We also identified miRNAs that are most likely involved in Culex immune responses to viruses and Wolbachia bacteria, and show the utility of using small RNA to detect antiviral immune pathways including piRNAs against some pathogens. Collectively, these findings show that deep sequencing of small RNA can be used for virus discovery and surveillance. One could also conceive that such work could be accomplished in various locations across the world and over time to better understand patterns of mosquito infection and immune response to many vector-borne diseases in field samples.

Reassessment of miRNA variant (isomiRs) composition by small RNA sequencing

IsomiRs, sequence variants of mature microRNAs, are usually detected and quantified using high-throughput sequencing. Many examples of their biological relevance have been reported, but sequencing artifacts identified as artificial variants might bias biological inference and therefore need to be ideally avoided. We conducted a comprehensive evaluation of 10 different small RNA sequencing protocols, exploring both a theoretically isomiR-free pool of synthetic miRNAs and HEK293T cells. We calculated that, with the exception of two protocols, less than 5% of miRNA reads can be attributed to library preparation artifacts. Randomized-end adapter protocols showed superior accuracy, with 40% of true biological isomiRs. Nevertheless, we demonstrate concordance across protocols for selected miRNAs in non-templated uridyl additions. Notably, NTA-U calling and isomiR target prediction can be inaccurate when using protocols with poor single-nucleotide resolution. Our results highlight the relevance of protocol choice for biological isomiRs detection and annotation, which has key potential implications for biomedical applications.

ALL-tRNAseq enables robust tRNA profiling in tissue samples

Transfer RNAs (tRNAs) are small adaptor RNAs essential for mRNA translation. Alterations in the cellular tRNA population can directly affect mRNA decoding rates and translational efficiency during cancer development and progression. To evaluate changes in the composition of the tRNA pool, multiple sequencing approaches have been developed to overcome reverse transcription blocks caused by the stable structures of these molecules and their numerous base modifications. However, it remains unclear whether current sequencing protocols faithfully capture tRNAs existing in cells or tissues. This is specifically challenging for clinical tissue samples that often present variable RNA qualities. For this reason, we developed ALL-tRNAseq, which combines the highly processive MarathonRT and RNA demethylation for the robust assessment of tRNA expression, together with a randomized adapter ligation strategy prior to reverse transcription to assess tRNA fragmentation levels in both cell lines and tissues. Incorporation of tRNA fragments not only informed on sample integrity but also significantly improved tRNA profiling of tissue samples. Our data showed that our profiling strategy effectively improves classification of oncogenic signatures in glioblastoma and diffuse large B-cell lymphoma tissues, particularly for samples presenting higher levels of RNA fragmentation, further highlighting the utility of ALL-tRNAseq for translational research.

sRNAtoolbox: Dockerized Analysis of Small RNA Sequencing Data in Model and Non-model Species

The current versions of the microRNA databases MiRgeneDB, miRBase, and PmiREN contain annotations for a total of 358 different species. Public repositories, however, host small RNA sequencing data for over 800 species. This discrepancy implies that microRNA research is also very active in species that neither have an available high-quality genome assembly nor annotations for microRNAs or other types of noncoding genes. These cases are particularly challenging to analyze because reference sequences need to be collected from different sources and processed and formatted appropriately so that the dedicated small RNA analysis tools can make use of them. In this protocol we describe how small RNA sequencing data can be easily analyzed by means of a dockerized version of the well-established sRNAtoolbox/sRNAbench small RNA tools. We outline the analysis of two publicly available datasets to demonstrate basic aspects like the preparation of the local database, expression profiling, or differential expression analysis as well as more advanced features such as quantification of exogenous RNA content and data analysis in non-model species.

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.