Ticks are ectoparasites that modulate host responses to sustain prolonged blood feeding, and in Ixodes ricinus, salivary microRNAs (miRNAs) represent promising candidates for manipulating host gene expression. Using phylogenetic footprinting combined with cooperative targeting analysis, we predicted deeply conserved miRNA-mRNA interactions that appear to contribute to tick lifecycle completion and tick's ability to parasitize diverse vertebrate hosts. Our analysis suggests that conserved tick miRNAs may mimic host miRNAs by exploiting shared target sites on host transcripts, potentially influencing regulatory circuits active during tick feeding and cross-kingdom communication. We identified 12 core tick salivary miRNAs-highly expressed and enriched in conserved target sites-predicted to cooperatively target 22 human genes, several of which are associated with processes relevant to parasitism establishment. These predicted cooperative interactions involve hub genes such as PDGFRA and NRG1, which are components of MAPK and PI3K-AKT pathways that regulate immune defense, tissue repair, and sensory responses. Our results suggest that I. ricinus miRNAs could potentially influence host homeostasis through cooperative targeting mechanisms, and they introduce a broadly applicable framework that leverages evolutionary conservation to detect meaningful cross-species miRNA-mRNA interactions.

MicroRNAs are by far the most extensively studied and best-characterized class of small RNAs. This depth of knowledge may partly explain the persistent eagerness in the field to classify new sequences as microRNAs-since labeling a molecule as such often implies immediate functional insight and biological relevance. However, this enthusiasm has led to the publication and deposition of thousands of spurious microRNA entries in public repositories, increasing the risk of misinterpretation-particularly when these sequences are linked to human disease. To address this, we present a concise, four-part checklist for evaluating putative microRNAs, based on structural, genomic, functional, and expression-related criteria. Applying it to "Mir-690" and nine other examples, which are not listed in the curated microRNA database MirGeneDB, but highly cited, we find no support for their classification as microRNAs. They fail to meet core microRNA annotation criteria, lack structural features, conservation of targets, and reproducible or microRNA-typical expression evidence. Our findings underscore the need for heightened scrutiny and standardized validation in microRNA research. The checklist we provide offers a practical tool for authors, editors, reviewers, and readers to critically assess microRNA claims-particularly when such sequences are proposed as biomarkers or therapeutic targets in translational studies.

Background: Platelet-derived Extracellular Vesicles, or “Platelet Dust” (PD), are reported as the most-abundant extracellular vesicles in plasma. However, the PD molecular content, especially the small RNA profile, is still poorly characterized. This study aims to characterize PD and other extracellular vesicles (EVs) in patients with non-small-cell lung cancer (NSCLC), focusing on their small RNA signatures and diagnostic potential. Methods: The EVs were isolated directly from the plasma of healthy donors and patients with NSCLC using the surface markers CD9, CD63, CD81 (overall EVs), and CD61 (PD). Small RNA sequencing was then performed to comprehensively profile the miRNAs. Results: Our analysis revealed distinct small RNA profiles in the EVs and the PD from the patients with NSCLC. The EVs (CD9-, CD63-, and CD81-positive) showed the enrichment of four miRNAs and the depletion of ten miRNAs, while the PD (CD61-positive) exhibited a more complex profile, with nineteen miRNAs enriched and nine miRNAs depleted in the patients with NSCLC compared to those of the healthy controls. Conclusions: This exploratory study enhances our understanding of miRNA composition within different plasma vesicle populations, shedding light on the biology of plasma vesicles and their contents. Furthermore, utilizing an extracellular vesicle isolation method with potential clinical applicability offers the prospect of improved cancer characterization and detection by selecting the most informative subpopulation of plasma vesicles.

BackgroundTicks, hematophagous Acari, pose a significant threat by transmitting various pathogens to their vertebrate hosts during feeding. Despite advances in tick genomics, high-quality genomes were lacking until recently, particularly in the genus Ixodes, which includes the main vectors of Lyme disease.ResultsHere, we present the genome sequences of four tick species, derived from a single female individual, with a particular focus on the European species Ixodes ricinus, achieving a chromosome-level assembly. Additionally, draft assemblies were generated for the three other Ixodes species, I. persulcatus, I. pacificus, and I. hexagonus. The quality of the four genomes and extensive annotation of several important gene families have allowed us to study the evolution of gene repertoires at the level of the genus Ixodes and of the tick group. We have determined gene families that have undergone major amplifications during the evolution of ticks, while an expression atlas obtained for I. ricinus reveals striking patterns of specialization both between and within gene families. Notably, several gene family amplifications are associated with a proliferation of single-exon genes-most strikingly for fatty acid elongases and sulfotransferases.ConclusionsThe integration of our data with existing genomes establishes a solid framework for the study of gene evolution, improving our understanding of tick biology. In addition, our work lays the foundations for applied research and innovative control targeting these organisms.

Plasmodium is an obligate intracellular parasite that has numerous interactions with different hosts during its elaborate life cycle. This is also the case for the other parasites belonging to the same phylum Apicomplexa. In this study, we bioinformatically identified the components of the multi-synthetase complexes (MSCs) of several Apicomplexa parasites and modelled their assembly using AlphaFold2. It appears that none of these MSCs resemble the two MSCs that we have identified and characterized in Plasmodium. Indeed, tRip, the central protein involved in the association of the two Plasmodium MSCs is different from its homologues, suggesting also that the tRip-dependent import of exogenous tRNAs is not conserved in other apicomplexan parasites. Based on this observation, we searched for obvious differences that could explain the singularity of Plasmodium protein synthesis by comparing tRNA genes and amino acid usage in the different genomes. We noted a contradiction between the large number of asparagine residues used in Plasmodium proteomes and the single gene encoding the tRNA that inserts them into proteins. This observation remains true for all the Plasmodia strains studied, even those that do not contain long asparagine homorepeats.

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.

Liquid biopsy-derived RNA sequencing (lbRNA-seq) exhibits significant promise for clinic-oriented cancer diagnostics due to its non-invasiveness and ease of repeatability. Despite substantial advancements, obstacles like technical artefacts and process standardisation impede seamless clinical integration. Alongside addressing technical aspects such as normalising fluctuating low-input material and establishing a standardised clinical workflow, the lack of result validation using independent datasets remains a critical factor contributing to the often low reproducibility of liquid biopsy-detected biomarkers.Considering the outlined drawbacks, our objective was to establish a workflow/methodology characterised by: 1. Harness the rich diversity of biological features accessible through lbRNA-seq data, encompassing a holistic range of molecular and functional attributes. These components are seamlessly integrated via a Machine Learning-based Ensemble Classification framework, enabling a unified and comprehensive analysis of the intricate information encoded within the data. 2. Implementing and rigorously benchmarking intra-sample normalisation methods to heighten their relevance within clinical settings. 3. Thoroughly assessing its efficacy across independent test sets to ascertain its robustness and potential utility.Using ten datasets from several studies comprising three different sources of biological material, we first show that while the best-performing normalisation methods depend strongly on the dataset and coupled Machine Learning method, the rather simple Counts Per Million method is generally very robust, showing comparable performance to cross-sample methods. Subsequently, we demonstrate that the innovative biofeature types introduced in this study, such as the Fraction of Canonical Transcript, harbour complementary information. Consequently, their inclusion consistently enhances prediction power compared to models relying solely on gene expression-based biofeatures. Finally, we demonstrate that the workflow is robust on completely independent datasets, generally from different labs and/or different protocols. Taken together, the workflow presented here outperforms generally employed methods in prediction accuracy and may hold potential for clinical diagnostics application due to its specific design.

We present a major update of MirGeneDB (3.0), the manually curated animal microRNA gene database. Beyond moving to a new server and the creation of a computational mirror, we have expanded the database with the addition of 33 invertebrate species, including representatives of 5 previously unsampled phyla, and 6 mammal species. MirGeneDB now contains entries for 21 822 microRNA genes (5160 of these from the new species) belonging to 1743 microRNA families. The inclusion of these new species allowed us to refine both the evolutionary node of appearance of a number of microRNA genes/families, as well as MirGeneDB's phylogenetically informed nomenclature system. Updated covariance models of all microRNA families, along with all smallRNA read data are now downloadable. These enhanced annotations will allow researchers to analyze microRNA properties such as secondary structure and features of their biogenesis within a robust phylogenetic context and without the database plagued with numerous false positives and false negatives. In light of these improvements, MirGeneDB 3.0 will assume the responsibility for naming conserved novel metazoan microRNAs. MirGeneDB is part of RNAcentral and Elixir Norway and is publicly and freely available at mirgenedb.org.

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. Gomez-Marti ' n et al. (2023)1 and Van Eijndhoven et al. (2021).2

Michael Hackenberg: Books in Artisan Homes of Sixteenth-Century Germany. In: Journal of Library History 21, 1986, S. 72-91.

RNA-sequencing has become one of the most used high-throughput approaches to gain knowledge about the expression of all different RNA subpopulations. However, technical artifacts, either introduced during library preparation and/or data analysis, can influence the detected RNA expression levels. A critical step, especially in large and low input datasets or studies, is data normalization, which aims at eliminating the variability in data that is not related to biology. Many normalization methods have been developed, each of them relying on different assumptions, making the selection of the appropriate normalization strategy key to preserve biological information. To address this, we developed NormSeq, a free web-server tool to systematically assess the performance of normalization methods in a given dataset. A key feature of NormSeq is the implementation of information gain to guide the selection of the best normalization method, which is crucial to eliminate or at least reduce non-biological variability. Altogether, NormSeq provides an easy-to-use platform to explore different aspects of gene expression data with a special focus on data normalization to help researchers, even without bioinformatics expertise, to obtain reliable biological inference from their data. NormSeq is freely available at: https://arn.ugr.es/normSeq.

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 drought-upregulated 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.

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.

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.

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.

Malignant melanoma (MM) can spread to other organs and is resistant in part due to the presence of cancer stem cell subpopulations (CSCs). While a controversial high dose of interferon-alpha (IFN-α) has been used to treat non-metastatic high-risk melanoma, it comes with undesirable side effects. In this study, we evaluated the effect of low and high doses of IFN-α on CSCs by analyzing ALDH activity, side population and specific surface markers in established and patient-derived primary cell lines. We also assessed the clonogenicity, migration and tumor initiation capacities of IFN-α treated CSCs. Additionally, we investigated genomic modulations related to stemness properties using microRNA sequencing and microarrays. The effect of IFN-α on CSCs-derived exosomes was also analyzed using NanoSight and liquid chromatography (LC-HRMS)-based metabolomic analysis, among others. Our results showed that even low doses of IFN-α reduced CSC formation and stemness properties, and led to a significant decrease in the ability to form tumors in mice xenotransplants. IFN-α also modulated the expression of genes and microRNAs involved in several cancer processes and metabolomics of released exosomes. Our work suggests the utility of low doses of interferon, combined with the analysis of metabolic biomarkers, as a potential clinical approach against the aggressiveness of CSCs in melanoma.

Abstract 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.

Viruses use miRNAs to enable efficient replication, control host defense mechanisms, and regulate latent infection. Herpes simplex virus 1 (HSV-1) expresses multiple miRNAs, whose functions are largely unknown. The evolutionary conservation of many HSV-1 miRNAs in the closely related HSV-2 suggests their functional importance. miRNAs, similar to other transcripts, can undergo various posttranscriptional modifications that may affect their biogenesis, stability and targeting. To investigate whether editing occurs in HSV-1 miRNAs, we sequenced samples from latently infected human ganglia. We show that one of the six HSV-1 miRNAs (miR-H2 to -H8) that define HSV-1 latency, miR-H2, exhibits A-to-I hyperediting within the miRNA seed sequence. We observed the same specific miR-H2 hyperediting phenomenon in miRNAs isolated from the ganglia of latently infected mice and, to a lesser extent, during productive infection in cultured cells. Curiously, we found no evidence of editing of the encoded HSV-2 homolog in latently infected mice or in cultured cells. The efficient loading of the edited miRNAs onto the RISC complex, indicates their ability to function as miRNAs. To investigate the potential of the edited miRNA to alter mRNA targeting, we predicted the host and viral targets for the modified miRNAs. Nucleotide substitution in the seed region significantly increased the number of potential host and viral targets. Most notably, ICP4, an essential viral protein, was predicted to be an additional target. Using transfection assays, we demonstrated that edited miRNAs have the potential to regulate ICP4 in addition to the previously identified target ICP0. Our study identifies a specific hyperedited HSV-1 mRNA, miR-H2, and highlights how the virus can use a single miRNA to target multiple transcripts during persistent, latent infection.Importance Herpes simplex virus 1 is an important human pathogen and 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 (miR-H2-e) can target ICP4, an essential viral protein. Interestingly, we found no evidence of hyperdating of its homolog, miR-H2, which is expressed by the closely related virus HSV-2. The discovery of posttranslational 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.

Despite the diversity of liquid biopsy transcriptomic repertoire, numerous studies often exploit only a single RNA type signature for diagnostic biomarker potential. This frequently results in insufficient sensitivity and specificity necessary to reach diagnostic utility. Combinatorial biomarker approaches may offer a more reliable diagnosis. Here, we investigated the synergistic contributions of circRNA and mRNA signatures derived from blood platelets as biomarkers for lung cancer detection. We developed a comprehensive bioinformatics pipeline permitting an analysis of platelet-circRNA and mRNA derived from non-cancer individuals and lung cancer patients. An optimal selected signature is then used to generate the predictive classification model using machine learning algorithm. Using an individual signature of 21 circRNA and 28 mRNA, the predictive models reached an area under the curve (AUC) of 0.88 and 0.81, respectively. Importantly, combinatorial analysis including both types of RNAs resulted in an 8-target signature (6 mRNA and 2 circRNA), enhancing the differentiation of lung cancer from controls (AUC of 0.92). Additionally, we identified five biomarkers potentially specific for early-stage detection of lung cancer. Our proof-of-concept study presents the first multi-analyte-based approach for the analysis of platelets-derived biomarkers, providing a potential combinatorial diagnostic signature for lung cancer detection.

Trichomes are specialised epidermal cells developed in the aerial surface of almost every terrestrial plant. These structures form physical barriers, which combined with their capability of synthesis of complex molecules, prevent plagues from spreading and confer trichomes a key role in the defence against herbivores. In this work, the tomato gene HAIRPLUS (HAP) that controls glandular trichome density in tomato plants was characterised. HAP belongs to a group of proteins involved in histone tail modifications although some also bind methylated DNA. HAP loss of function promotes epigenomic modifications in the tomato genome reflected in numerous differentially methylated cytosines and causes transcriptomic changes in hap mutant plants. Taken together, these findings demonstrate that HAP links epigenome remodelling with multicellular glandular trichome development and reveal that HAP is a valuable genomic tool for pest resistance in tomato breeding.