miRTrace reveals the organismal origins of microRNA sequencing data

Analysis of miRNAs in milk of four livestock species

BACKGROUND

Milk is essential for mammalian nutrition because it provides vital nutrients for growth and development. Milk composition, which is influenced by genetic and environmental factors, supports lactation, a complex process crucial for milk production and quality. Recent research has focused on noncoding RNAs, particularly microRNAs (miRNAs), which are present in body fluids and regulate gene expression post-transcriptionally. This study comprehensively characterizes miRNAs in milk of four livestock species, namely Bubalus bubalis, Capra hircus, Equus asinus, and Ovis aries and identifies potential target genes.

RESULTS

High-throughput sequencing of milk RNA resulted in distinct read counts across species: B. bubalis (8,790,441 reads), C. hircus (12,976,275 reads), E. asinus (9,385,067 reads), and O. aries (7,295,297 reads). E. asinus had the highest RNA mapping rate (94.6%) and O. aries the lowest (84.8%). A substantially greater proportion of miRNAs over other small RNAs was observed for the donkey milk sample (7.74%) compared to buffalo (0.87%), goat (1.57%), and sheep (1.12%). Shared miRNAs, which included miR-200a, miR-200b, miR-200c, and miR-23a among others, showed varying expression levels across species, confirmed by qPCR analysis. Functional annotation of predicted miRNA target genes highlighted diverse roles, with an enrichment in functions linked to metabolism and immunity. Pathway analysis identified immune response pathways as significant, with several miRNAs targeting specific genes across species, suggesting their regulatory function in milk.

CONCLUSIONS

Both conserved and species-specific miRNAs were detected in milk of the investigated species. The identified target genes of these miRNAs have important roles in neonatal development, adaptation, growth, and immune response. Furthermore, they influence milk and meat production traits in livestock.

Insights into early animal evolution from the genome of the xenacoelomorph worm Xenoturbella bocki

The evolutionary origins of Bilateria remain enigmatic. One of the more enduring proposals highlights similarities between a cnidarian-like planula larva and simple acoel-like flatworms. This idea is based in part on the view of the Xenacoelomorpha as an outgroup to all other bilaterians which are themselves designated the Nephrozoa (protostomes and deuterostomes). Genome data can provide important comparative data and help understand the evolution and biology of enigmatic species better. Here, we assemble and analyze the genome of the simple, marine xenacoelomorph Xenoturbella bocki, a key species for our understanding of early bilaterian evolution. Our highly contiguous genome assembly of X. bocki has a size of ~111 Mbp in 18 chromosome-like scaffolds, with repeat content and intron, exon, and intergenic space comparable to other bilaterian invertebrates. We find X. bocki to have a similar number of genes to other bilaterians and to have retained ancestral metazoan synteny. Key bilaterian signaling pathways are also largely complete and most bilaterian miRNAs are present. Overall, we conclude that X. bocki has a complex genome typical of bilaterians, which does not reflect the apparent simplicity of its body plan that has been so important to proposals that the Xenacoelomorpha are the simple sister group of the rest of the Bilateria.

Space radiation damage rescued by inhibition of key spaceflight associated miRNAs

Our previous research revealed a key microRNA signature that is associated with spaceflight that can be used as a biomarker and to develop countermeasure treatments to mitigate the damage caused by space radiation. Here, we expand on this work to determine the biological factors rescued by the countermeasure treatment. We performed RNA-sequencing and transcriptomic analysis on 3D microvessel cell cultures exposed to simulated deep space radiation (0.5 Gy of Galactic Cosmic Radiation) with and without the antagonists to three microRNAs: miR-16-5p, miR-125b-5p, and let-7a-5p (i.e., antagomirs). Significant reduction of inflammation and DNA double strand breaks (DSBs) activity and rescue of mitochondria functions are observed after antagomir treatment. Using data from astronaut participants in the NASA Twin Study, Inspiration4, and JAXA missions, we reveal the genes and pathways implicated in the action of these antagomirs are altered in humans. Our findings indicate a countermeasure strategy that can potentially be utilized by astronauts in spaceflight missions to mitigate space radiation damage.

Spatiotemporal expression and control of haemoglobin in space

It is now widely recognised that the environment in space activates a diverse set of genes involved in regulating fundamental cellular pathways. This includes the activation of genes associated with blood homoeostasis and erythropoiesis, with a particular emphasis on those involved in globin chain production. Haemoglobin biology provides an intriguing model for studying space omics, as it has been extensively explored at multiple -omic levels, spanning DNA, RNA, and protein analyses, in both experimental and clinical contexts. In this study, we examined the developmental expression of haemoglobin over time and space using a unique suite of multi-omic datasets available on NASA GeneLab, from the NASA Twins Study, the JAXA CFE study, and the Inspiration4 mission. Our findings reveal significant variations in globin gene expression corresponding to the distinct spatiotemporal characteristics of the collected samples. This study sheds light on the dynamic nature of globin gene regulation in response to the space environment and provides valuable insights into the broader implications of space omics research.

Grandfathers-to-Grandsons Transgenerational Transmission of Exercise Positive Effects on Cognitive Performance

Physical exercise is a robust lifestyle intervention known for its enhancement of cognitive abilities. Nevertheless, the extent to which these benefits can be transmitted across generations (intergenerational inheritance to F1, and transgenerational to F2 and beyond) remains a topic of limited comprehension. We have already shown that cognitive improvements resulting from physical exercise can be inherited from parents to their offspring, proving intergenerational effects. So, we set out to explore whether these enhancements might extend transgenerationally, impacting the F2 generation. In this study, we initially examined the behavioral traits of second generation (F2) male mice, whose grandfathers (F0) had an exercise intervention. Our findings revealed that F2 mice with physically active grandpaternal F0 progenitors displayed significantly improved memory recall, encompassing both spatial and non-spatial information when compared to their counterparts from sedentary F0 progenitors, and proving for the first time the transgenerational inheritance of physical exercise induced cognitive enhancement. Surprisingly, while F2 memory improved (as was the case with F1), adult hippocampal neurogenesis remained unchanged between experimental and control groups (unlike in F1). Additionally, our analysis of small RNA sequences in the hippocampus identified 35 differentially expressed miRNAs linked to important brain function categories. Notably, two of these miRNAs, miRNA-144 and miRNA-298, displayed a robust negative correlation with cognitive performance. These findings highlight the enduring transgenerational transmission of cognitive benefits associated with exercise, even after two generations, suggesting that moderate exercise training can have lasting positive effects, possibly orchestrated by a specific set of miRNAs that exert their influence across multiple generations.

A novel multi-stage enrichment workflow and comprehensive characterization for HEK293F-derived extracellular vesicles

Extracellular vesicles (EVs) are emerging as a promising drug delivery vehicle as they are biocompatible and capable of targeted delivery. However, clinical translation of EVs remains challenging due to the lack of standardized and scalable manufacturing protocols to consistently isolate small EVs (sEVs) with both high yield and high purity. The heterogenous nature of sEVs leading to unknown composition of biocargos causes further pushback due to safety concerns. In order to address these issues, we developed a robust quality-controlled multi-stage process to produce and isolate sEVs from human embryonic kidney HEK293F cells. We then compared different 2-step and 3-step workflows for eliminating protein impurities and cell-free nucleic acids to meet acceptable limits of regulatory authorities. Our results showed that sEV production was maximized when HEK293F cells were grown at high-density stationary phase in semi-continuous culture. The novel 3-step workflow combining tangential flow filtration, sucrose-cushion ultracentrifugation and bind-elute size-exclusion chromatography outperformed other methods in sEV purity while still preserved high yield and particle integrity. The purified HEK293F-derived sEVs were thoroughly characterized for identity including sub-population analysis, content profiling including proteomics and miRNA sequencing, and demonstrated excellent preclinical safety profile in both in-vitro and in-vivo testing. Our rigorous enrichment workflow and comprehensive characterization will help advance the development of EVs, particularly HEK293F-derived sEVs, to be safe and reliable drug carriers for therapeutic applications.

Evolution of microRNAs in Amoebozoa and implications for the origin of multicellularity

MicroRNAs (miRNAs) are important and ubiquitous regulators of gene expression in both plants and animals. They are thought to have evolved convergently in these lineages and hypothesized to have played a role in the evolution of multicellularity. In line with this hypothesis, miRNAs have so far only been described in few unicellular eukaryotes. Here, we investigate the presence and evolution of miRNAs in Amoebozoa, focusing on species belonging to Acanthamoeba, Physarum and dictyostelid taxonomic groups, representing a range of unicellular and multicellular lifestyles. miRNAs that adhere to both the stringent plant and animal miRNA criteria were identified in all examined amoebae, expanding the total number of protists harbouring miRNAs from 7 to 15. We found conserved miRNAs between closely related species, but the majority of species feature only unique miRNAs. This shows rapid gain and/or loss of miRNAs in Amoebozoa, further illustrated by a detailed comparison between two evolutionary closely related dictyostelids. Additionally, loss of miRNAs in the Dictyostelium discoideum drnB mutant did not seem to affect multicellular development and, hence, demonstrates that the presence of miRNAs does not appear to be a strict requirement for the transition from uni- to multicellular life.

A time course analysis through diapause reveals dynamic temporal patterns of microRNAs associated with endocrine regulation in the butterfly Pieris napi

Organisms inhabiting highly seasonal environments must cope with a wide range of environmentally induced challenges. Many seasonal challenges require extensive physiological modification to survive. In winter, to survive extreme cold and limited resources, insects commonly enter diapause, which is an endogenously derived dormant state associated with minimized cellular processes and low energetic expenditure. Due to the high degree of complexity involved in diapause, substantial cellular regulation is required, of which our understanding primarily derives from the transcriptome via messenger RNA expression dynamics. Here we aim to advance our understanding of diapause by investigating microRNA (miRNA) expression in diapausing and direct developing pupae of the butterfly Pieris napi. We identified coordinated patterns of miRNA expression throughout diapause in both head and abdomen tissues of pupae, and via miRNA target identification, found several expression patterns to be enriched for relevant diapause-related physiological processes. We also identified two candidate miRNAs, miR-14-5p and miR-2a-3p, that are likely involved in diapause progression through their activity in the ecdysone pathway, a critical regulator of diapause termination. miR-14-5p targets phantom, a gene in the ecdysone synthesis pathway, and is upregulated early in diapause. miR-2a-3p has been found to be expressed in response to ecdysone, and is upregulated during diapause termination. Together, the expression patterns of these two miRNAs match our current understanding of the timing of hormonal regulation of diapause in P. napi and provide interesting candidates to further explore the mechanistic role of microRNAs in diapause regulation.

Tissue-specific profiling of age-dependent miRNAomic changes in Caenorhabditis elegans

Ageing exhibits common and distinct features in various tissues, making it critical to decipher the tissue-specific ageing mechanisms. MiRNAs are essential regulators in ageing and are recently highlighted as a class of intercellular messengers. However, little is known about the tissue-specific transcriptomic changes of miRNAs during ageing. C. elegans is a well-established model organism in ageing research. Here, we profile the age-dependent miRNAomic changes in five isolated worm tissues. Besides the diverse ageing-regulated miRNA expression across tissues, we discover numerous miRNAs in the tissues without their transcription. We further profile miRNAs in the extracellular vesicles and find that worm miRNAs undergo inter-tissue trafficking via these vesicles in an age-dependent manner. Using these datasets, we uncover the interaction between body wall muscle-derived mir-1 and DAF-16/FOXO in the intestine, suggesting mir-1 as a messenger in inter-tissue signalling. Taken together, we systematically investigate worm miRNAs in the somatic tissues and extracellular vesicles during ageing, providing a valuable resource to study tissue-autonomous and nonautonomous functions of miRNAs in ageing.

Identification and characterization of microRNAs in Biomphalaria tenagophila and comparative analysis of their expression in Schistosoma mansoni-resistant and -susceptible snail populations

BACKGROUND

Schistosomiasis is a neglected tropical disease caused by Schistosoma and affects over 240 million people worldwide. One of the most prominent causative agents is Schistosoma mansoni, which develops inside the intermediate host. Biomphalaria tenagophila is the second most important vector of schistosomiasis in Brazil and the Taim population is completely resistant to infection by S. mansoni.

OBJECTIVE

This study aims to identify and characterize B. tenagophila microRNAs (miRNAs) and evaluate their differential expression in S. mansoni-susceptible and -resistant populations of B. tenagophila.

METHODS

Two populations of B. tenagophila snails, susceptible and resistant to S. mansoni infection, were used to investigate the small RNA response of these snails after being infected with the parasite. Small RNA sequencing and quantitative real-time PCR were employed to identify and validate differentially expressed miRNAs. Bioinformatics analysis were performed to identify miRNA precursors and mature and evaluate their differential expression.

FINDINGS

The study predicted 173 mature miRNAs and 123 precursors. Among them were six Lophotrochozoa-specific miRNAs, three mollusk-specific miRNAs, and six pre-miRNAs in a cluster. The small RNA sequencing and RT-PCR of B. tenagophila samples allowed assessing the expression patterns of miRNAs.

MAIN CONCLUSIONS

The results obtained may support future studies in Biomphalaria spp., generating a global impact on disease control.

Age-Dependent Alterations in Semen Parameters and Human Sperm MicroRNA Profile

The trend to delay parenthood is increasing, impacting fertility and reproductive outcomes. Advanced paternal age (APA), defined as men's age above 40 years at conception, has been linked with testicular impairment, abnormal semen parameters, and poor reproductive and birth outcomes. Recently, the significance of sperm microRNA for fertilization and embryonic development has emerged. This work aimed to investigate the effects of men's age on semen parameters and sperm microRNA profiles. The ejaculates of 333 Portuguese men were collected between 2018 and 2022, analyzed according to WHO guidelines, and a density gradient sperm selection was performed. For microRNA expression analysis, 16 normozoospermic human sperm samples were selected and divided into four age groups: ≤30, 31-35, 36-40, and >40 years. microRNA target genes were retrieved from the miRDB and TargetScan databases and Gene Ontology analysis was performed using the DAVID tool. No significant correlation was found between male age and conventional semen parameters, except for volume. Fifteen differentially expressed microRNAs (DEMs) between groups were identified. Enrichment analysis suggested the involvement of DEMs in the sperm of men with advanced age in critical biological processes like embryonic development, morphogenesis, and male gonad development. Targets of DEMs were involved in signaling pathways previously associated with the ageing process, including cellular senescence, autophagy, insulin, and mTOR pathways. These results suggest that although conventional semen parameters were not affected by men's age, alterations in microRNA regulation may occur and be responsible for poor fertility and reproductive outcomes associated with APA.

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.

Historical RNA expression profiles from the extinct Tasmanian tiger

Paleogenomics continues to yield valuable insights into the evolution, population dynamics, and ecology of our ancestors and other extinct species. However, DNA sequencing cannot reveal tissue-specific gene expression, cellular identity, or gene regulation, which are only attainable at the transcriptional level. Pioneering studies have shown that useful RNA can be extracted from ancient specimens preserved in permafrost and historical skins from extant canids, but no attempts have been made so far on extinct species. We extract, sequence, and analyze historical RNA from muscle and skin tissue of a ∼130-year-old Tasmanian tiger (Thylacinus cynocephalus) preserved in desiccation at room temperature in a museum collection. The transcriptional profiles closely resemble those of extant species, revealing specific anatomical features such as slow muscle fibers or blood infiltration. Metatranscriptomic analysis, RNA damage, tissue-specific RNA profiles, and expression hotspots genome-wide further confirm the thylacine origin of the sequences. RNA sequences are used to improve protein-coding and noncoding annotations, evidencing missing exonic loci and the location of ribosomal RNA genes while increasing the number of annotated thylacine microRNAs from 62 to 325. We discover a thylacine-specific microRNA isoform that could not have been confirmed without RNA evidence. Finally, we detect traces of RNA viruses, suggesting the possibility of profiling viral evolution. Our results represent the first successful attempt to obtain transcriptional profiles from an extinct animal species, providing thought-to-be-lost information on gene expression dynamics. These findings hold promising implications for the study of RNA molecules across the vast collections of natural history museums and from well-preserved permafrost remains.

ncOrtho: efficient and reliable identification of miRNA orthologs

MicroRNAs (miRNAs) are post-transcriptional regulators that finetune gene expression via translational repression or degradation of their target mRNAs. Despite their functional relevance, frameworks for the scalable and accurate detection of miRNA orthologs are missing. Consequently, there is still no comprehensive picture of how miRNAs and their associated regulatory networks have evolved. Here we present ncOrtho, a synteny informed pipeline for the targeted search of miRNA orthologs in unannotated genome sequences. ncOrtho matches miRNA annotations from multi-tissue transcriptomes in precision, while scaling to the analysis of hundreds of custom-selected species. The presence-absence pattern of orthologs to 266 human miRNA families across 402 vertebrate species reveals four bursts of miRNA acquisition, of which the most recent event occurred in the last common ancestor of higher primates. miRNA families are rarely modified or lost, but notable exceptions for both events exist. miRNA co-ortholog numbers faithfully indicate lineage-specific whole genome duplications, and miRNAs are powerful markers for phylogenomic analyses. Their exceptionally low genetic diversity makes them suitable to resolve clades where the phylogenetic signal is blurred by incomplete lineage sorting of ancestral alleles. In summary, ncOrtho allows to routinely consider miRNAs in evolutionary analyses that were thus far reserved to protein-coding genes.

Role of miRNAs in regulation of SA-mediated upregulation of genes involved in folate and methionine metabolism in foxtail millet

The effect of exogenous salicylic acid (SA) on folate metabolism and the related gene regulatory mechanisms is still unclear. In this study, the panicle of foxtail millet treated with different SA concentrations showed that 6 mM SA doubled the 5-methyltetrahydrofolate content compared to that of the control. An untargeted metabolomic analysis revealed that 275 metabolites were enriched in amino acid metabolic pathways. Significantly, the relative content of methionine (Met) after 6 mM SA treatment was 3.14 times higher than the control. Transcriptome analysis revealed that differentially expressed genes were mainly enriched in the folate and amino acid biosynthesis pathways (including Met, Cys, Pro, Ser et al.). The miRNA-mRNA interactions related to the folate and Met metabolic pathways were analyzed and several likely structural gene targets for miRNAs were identified, miRNA-seq analysis revealed that 33 and 51 miRNAs targeted 11 and 15 genes related to the folate and Met pathways, respectively. Eight key genes in the folate metabolism pathway were likely to be up-regulated by 14 new miRNAs and 20 new miRNAs up-regulated the 9 key genes in the Met metabolism pathway. The 6 miRNA-mRNA interactions related to the folate and Met metabolism pathways were verified by qRT-PCR, and consistent with the prediction. The results showed that DHFR1 gene expression level related to folate synthesis was directly up-regulated by Nov-m0139-3p with 3.8 times, but DHFR2 was down-regulated by Nov-m0731-5p with 0.62 times. The expression level of CYSC1 and APIP related to Met synthesis were up-regulated by Nov-m0461-5p and Nov-m0664-3p with 4.27 and 1.32 times, respectively. Our results suggested that exogenous SA could induce the folate and Met accumulated in the panicle of foxtail millet. The higher expression level of DHFR1, FTHFD, CYSC1 and APIP in the folate and Met metabolism pathway and their regulators, including Nov-m0139-3p, Nov-m0717-5p, Nov-m0461-5p and Nov-m0664-3p, could be responsible for these metabolites accumulation. This study lays the theoretical foundation for elucidating the post-transcription regulatory mechanisms of folate and Met metabolism.

MicroRNAs are deeply linked to the emergence of the complex octopus brain

Soft-bodied cephalopods such as octopuses are exceptionally intelligent invertebrates with a highly complex nervous system that evolved independently from vertebrates. Because of elevated RNA editing in their nervous tissues, we hypothesized that RNA regulation may play a major role in the cognitive success of this group. We thus profiled messenger RNAs and small RNAs in three cephalopod species including 18 tissues of the Octopus vulgaris. We show that the major RNA innovation of soft-bodied cephalopods is an expansion of the microRNA (miRNA) gene repertoire. These evolutionarily novel miRNAs were primarily expressed in adult neuronal tissues and during the development and had conserved and thus likely functional target sites. The only comparable miRNA expansions happened, notably, in vertebrates. Thus, we propose that miRNAs are intimately linked to the evolution of complex animal brains.

Inhibition of Respiratory Syncytial Virus Infection by Small Non-Coding RNA Fragments

Respiratory syncytial virus (RSV) causes acute lower respiratory tract infection in infants, immunocompromised individuals and the elderly. As the only current specific treatment options for RSV are monoclonal antibodies, there is a need for efficacious antiviral treatments against RSV to be developed. We have previously shown that a group of synthetic non-coding single-stranded DNA oligonucleotides with lengths of 25–40 nucleotides can inhibit RSV infection in vitro and in vivo. Based on this, herein, we investigate whether naturally occurring single-stranded small non-coding RNA (sncRNA) fragments present in the airways have antiviral effects against RSV infection. From publicly available sequencing data, we selected sncRNA fragments such as YRNAs, tRNAs and rRNAs present in human bronchoalveolar lavage fluid (BALF) from healthy individuals. We utilized a GFP-expressing RSV to show that pre-treatment with the selected sncRNA fragments inhibited RSV infection in A549 cells in vitro. Furthermore, by using a flow cytometry-based binding assay, we demonstrate that these naturally occurring sncRNAs fragments inhibit viral infection most likely by binding to the RSV entry receptor nucleolin and thereby preventing the virus from binding to host cells, either directly or via steric hindrance. This finding highlights a new function of sncRNAs and displays the possibility of using naturally occurring sncRNAs as treatments against RSV.

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

Paired miRNA and RNA sequencing provides a first insight into molecular defense mechanisms of Scapharca broughtonii during ostreid herpesvirus-1 infection

Ostreid herpesvirus 1 (OsHV-1) infection caused mortalities with relevant economic losses in bivalve aquaculture industry worldwide. Initially described as an oyster pathogen, OsHV-1 can infect other bivalve species, like the blood clam Scapharca broughtonii. However, at present, little is known about the molecular interactions during OsHV-1 infection in the blood clam. We produced paired miRNA and total RNA-seq data to investigate the blood clam transcriptional changes from 0 to 72 h after experimental infection with OsHV-1. High-throughput miRNA sequencing of 24 libraries revealed 580 conserved and 270 new blood clam miRNAs, whereas no genuine miRNA was identified for OsHV-1. Total 88-203 differently expressed miRNAs were identified per time point, mostly up-regulated and mainly targeting metabolic pathways. Most of the blood clam mRNAs, in contrast, were down-regulated up to 60 h post-injection, with the trend analysis revealing the activation of immune genes only when comparing the early and latest stage of infection. Taken together, paired short and long RNA data suggested a miRNA-mediated down-regulation of host metabolic and energetic processes as a possible antiviral strategy during early infection stages, whereas antiviral pathways appeared upregulated only at late infection.

The Small RNA Universe of Capitella teleta

RNAi is an evolutionarily fluid mechanism with dramatically different activities across animal phyla. One major group where there has been little investigation is annelid worms. Here, the small RNAs of the polychaete developmental model Capitella teleta are profiled across development. As is seen with nearly all animals, nearly 200 microRNAs were found with 58 high-confidence novel species. Greater miRNA diversity was associated with later stages consistent with differentiation of tissues. Outside miRNA, a distinct composition of other small RNA pathways was found. Unlike many invertebrates, an endogenous siRNA pathway was not observed, indicating pathway loss relative to basal planarians. No processively generated siRNA-class RNAs could be found arising from dsRNA precursors. This has a significant impact on RNAi technology development for this group of animals. Unlike the apparent absence of siRNAs, a significant population of piRNAs was observed. For many piRNAs, phasing and ping-pong biogenesis pathways were identified. Interestingly, piRNAs were found to be highly expressed during early development, suggesting a potential role in regulation in metamorphosis. Critically, the configuration of RNAi factors in C. teleta is found in other annelids and mollusks, suggesting that similar biology is likely to be present in the wider clade. This study is the first in providing comprehensive analysis of small RNAs in annelids.

Dissecting miRNA–Gene Networks to Map Clinical Utility Roads of Pharmacogenomics-Guided Therapeutic Decisions in Cardiovascular Precision Medicine

MicroRNAs (miRNAs) create systems networks and gene-expression circuits through molecular signaling and cell interactions that contribute to health imbalance and the emergence of cardiovascular disorders (CVDs). Because the clinical phenotypes of CVD patients present a diversity in their pathophysiology and heterogeneity at the molecular level, it is essential to establish genomic signatures to delineate multifactorial correlations, and to unveil the variability seen in therapeutic intervention outcomes. The clinically validated miRNA biomarkers, along with the relevant SNPs identified, have to be suitably implemented in the clinical setting in order to enhance patient stratification capacity, to contribute to a better understanding of the underlying pathophysiological mechanisms, to guide the selection of innovative therapeutic schemes, and to identify innovative drugs and delivery systems. In this article, the miRNA–gene networks and the genomic signatures resulting from the SNPs will be analyzed as a method of highlighting specific gene-signaling circuits as sources of molecular knowledge which is relevant to CVDs. In concordance with this concept, and as a case study, the design of the clinical trial GESS (NCT03150680) is referenced. The latter is presented in a manner to provide a direction for the improvement of the implementation of pharmacogenomics and precision cardiovascular medicine trials.

A comprehensive framework for analysis of microRNA sequencing data in metastatic colorectal cancer

Although microRNAs (miRNAs) contribute to all hallmarks of cancer, miRNA dysregulation in metastasis remains poorly understood. The aim of this work was to reliably identify miRNAs associated with metastatic progression of colorectal cancer (CRC) using novel and previously published next-generation sequencing (NGS) datasets generated from 268 samples of primary (pCRC) and metastatic CRC (mCRC; liver, lung and peritoneal metastases) and tumor adjacent tissues. Differential expression analysis was performed using a meticulous bioinformatics pipeline, including only bona fide miRNAs, and utilizing miRNA-tailored quality control and processing. Five miRNAs were identified as up-regulated at multiple metastatic sites Mir-210_3p, Mir-191_5p, Mir-8-P1b_3p [mir-141–3p], Mir-1307_5p and Mir-155_5p. Several have previously been implicated in metastasis through involvement in epithelial-to-mesenchymal transition and hypoxia, while other identified miRNAs represent novel findings. The use of a publicly available pipeline facilitates reproducibility and allows new datasets to be added as they become available. The set of miRNAs identified here provides a reliable starting-point for further research into the role of miRNAs in metastatic progression.

Expression Profiles of Estrogen-Regulated MicroRNAs in Cancer Cells

MicroRNAs play critical roles through their impact on posttranscriptional gene regulation. In cancer, they can act as oncogenes or tumor suppressors and can also function as biomarkers. Here, we describe a method for robust characterization of estrogen-regulated microRNA profiles. The activity of estrogen is mediated by two nuclear receptors, estrogen receptor alpha and estrogen receptor beta, and a transmembrane G-protein coupled estrogen receptor 1. This chapter details how to prepare cells for optimal estrogen response, directions for estrogen treatment, RNA extraction, different microRNA profiling approaches, and subsequent confirmations.

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

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

MapToCleave: High-throughput profiling of microRNA biogenesis in living cells

Previous large-scale studies have uncovered many features that determine the processing of microRNA (miRNA) precursors; however, they have been conducted in vitro. Here, we introduce MapToCleave, a method to simultaneously profile processing of thousands of distinct RNA structures in living cells. We find that miRNA precursors with a stable lower basal stem are more efficiently processed and also have higher expression in vivo in tissues from 20 animal species. We systematically compare the importance of known and novel sequence and structural features and test biogenesis of miRNA precursors from 10 animal and plant species in human cells. Lastly, we provide evidence that the GHG motif better predicts processing when defined as a structure rather than sequence motif, consistent with recent cryogenic electron microscopy (cryo-EM) studies. In summary, we apply a screening assay in living cells to reveal the importance of lower basal stem stability for miRNA processing and in vivo expression.

Large-scale identification of extracellular plant miRNAs in mammals implicates their dietary intake

Extracellular microRNAs (miRNAs) have been proposed to function in cross-kingdom gene regulation. Among these, plant-derived miRNAs of dietary origin have been reported to survive the harsh conditions of the human digestive system, enter the circulatory system, and regulate gene expression and metabolic function. However, definitive evidence supporting the presence of plant-derived miRNAs of dietary origin in mammals has been difficult to obtain due to limited sample sizes. We have developed a bioinformatics pipeline (ePmiRNA_finder) that provides strident miRNA classification and applied it to analyze 421 small RNA sequencing data sets from 10 types of human body fluids and tissues and comparative samples from carnivores and herbivores. A total of 35 miRNAs were identified that map to plants typically found in the human diet and these miRNAs were found in at least one human blood sample and their abundance was significantly different when compared to samples from human microbiome or cow. The plant-derived miRNA profiles were body fluid/tissue-specific and highly abundant in the brain and the breast milk samples, indicating selective absorption and/or the ability to be transported across tissue/organ barriers. Our data provide conclusive evidence for the presence of plant-derived miRNAs as a consequence of dietary intake and their cross-kingdom regulatory function within human circulating system.

Digging into bivalve miRNAomes: between conservation and innovation

Bivalves are a diverse mollusc group of economic and ecological importance. An evident resilience to pollution, parasites and extreme environments makes some bivalve species important models for studying adaptation and immunity. Despite substantial progress in sequencing projects of bivalves, information on non-coding genes and gene-regulatory aspects is still lacking. Here, we review the current repertoire of bivalve microRNAs (miRNAs), important regulators of gene expression in Metazoa. We exploited available short non-coding RNA (sncRNA) data for Pinctada martensii, Crassostrea gigas, Corbicula fluminea, Tegillarca granosa and Ruditapes philippinarum, and we produced new sncRNA data for two additional bivalves, the Mediterranean mussel Mytilus galloprovincialis and the blood clam Scapharca broughtonii. We found substantial heterogeneity and incorrect annotations of miRNAs; hence, we reannotated conserved miRNA families using recently established criteria for bona fide microRNA annotation. We found 106 miRNA families missing in the previously published bivalve datasets and 89 and 87 miRNA complements were identified in the two additional species. The overall results provide a homogeneous and evolutionarily consistent picture of miRNAs in bivalves and enable future comparative studies. The identification of two bivalve-specific miRNA families sheds further light on the complexity of transcription and its regulation in bivalve molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Forensic botany: time to embrace natural history collections, large scale environmental data and environmental DNA

Forensic botany is a diverse discipline that spans many aspects of plant sciences, particularly taxonomy, field botany, anatomy, and ecology. Internationally, there is a significant opportunity to expand the application of forensic botany in criminal investigations, especially war crimes, genocide, homicide, sexual violence, serious physical assault, illegal trade in endangered species and wildlife crime. In civil proceedings, forensic botany may, for example, be called upon in trade disputes such as accidental contamination of commodities. Despite the potential, there are barriers to the wider application of forensic botany in criminal cases; there is a widespread need to improve the efficiency of botanical trace evidence identification. This could partly be addressed by embracing innovations in image recognition and by accessing the huge quantity of specimens and images housed in natural history collections worldwide. Additionally, the recent advances in DNA sequencing technologies and the expansion of environmental DNA (eDNA) and forensic ecogenomics, offers opportunities to more rapidly provide species-level identifications. The impact of taphonomic processes upon vegetation, and vice versa, remains poorly understood; improved understanding of these interactions and their ecological impacts may be invaluable in improving clandestine burial search protocols.

Differential expression of microRNAs and tRNA fragments mediate the adaptation of the liver fluke Fasciola gigantica to its intermediate snail and definitive mammalian hosts

The tropical liver fluke Fasciola gigantica affects livestock and humans in many Asian countries, large parts of Africa, and parts of Europe. Despite the public health and economic impacts of F. gigantica, understanding of F. gigantica biology and how the complex lifecycle of this liver fluke is transcriptionally regulated remain unknown. Here, we tested the hypothesis that the regulatory small non-coding RNAs (sncRNAs), microRNAs (miRNAs) and tRNA-derived fragments (tRFs) play roles in the adaptation of F. gigantica to its intermediate and definitive hosts. We sequenced sncRNAs of eight lifecycle stages of F. gigantica. In total, 56 miRNAs from 33 conserved families and four Fasciola-specific miRNAs were identified. Expression analysis of miRNAs suggested clear stage-related patterns. By leveraging the existing transcriptomic data, we predicted a miRNA-based regulation of metabolism, transport, growth and developmental processes. Also, by comparing miRNA complement of F. gigantica with that of Fasciola hepatica, we detected a high level of conservation and identified differences in some miRNAs, which can be used to distinguish the two species. Moreover, we found that tRFs at each lifecycle stage were predominantly derived by tRNA^-Lys and tRNA^-Gly at 5' half sites, but relatively high expression was related to the buffalo-infecting stages. Taken together, we provided a comprehensive overview of the dynamic transcriptional changes of small RNAs that occur during the developmental stages of F. gigantica. This global analysis of F. gigantica lifecycle stages revealed new roles of miRNAs and tRFs in parasite development and will facilitate future research into understanding of fasciolosis pathobiology.

Ancient microRNA profiles of a 14,300-year-old canid samples confirm taxonomic origin and give glimpses into tissue-specific gene regulation from the Pleistocene

DNA sequencing is the current key technology for historic or ancient biological samples and has led to many exciting discoveries in the field of paleogenomics. However, functional insights into tissue identity, cellular composition or gene regulation cannot be gained from DNA. Recent analyses have shown that, under favorable conditions, RNA can also be sequenced from ancient samples, enabling studies at the transcriptomic and regulatory level. Analyzing ancient RNA data from a Pleistocene canid, we find hundreds of intact microRNAs that are taxonomically informative, show tissue-specificity and have functionally predictive characteristics. With an extraordinary age of 14,300 years, these microRNA sequences are by far the oldest ever reported. The authenticity of the sequences is further supported by a) the presence of canid / Caniformia-specific sequences that never evolved outside of this clade, b) tissue-specific expression patterns (cartilage, liver and muscle) that resemble those of modern dogs and c) RNA damage patterns that are clearly distinct from those of fresh samples. By performing computational microRNA-target enrichment analyses on the ancient sequences, we predict microRNA functions consistent with their tissue pattern of expression. For instance, we find a liver-specific microRNA that regulates carbohydrate metabolism and starvation responses in canids. In summary, we show that straightforward paleotranscriptomic microRNA analyses can give functional glimpses into tissue identity, cellular composition and gene regulatory activity of ancient samples and biological processes that took place in the Pleistocene, thus holding great promise for deeper insights into gene regulation in extinct animals based on ancient RNA sequencing. .

Circulating miRNA Spaceflight Signature Reveals Targets for Countermeasure Development

We have identified and validated a spaceflight-associated microRNA (miRNA) signature that is shared by rodents and humans in response to simulated, short-duration and long-duration spaceflight. Previous studies have identified miRNAs that regulate rodent responses to spaceflight in low-Earth orbit, and we have confirmed the expression of these proposed spaceflight-associated miRNAs in rodents reacting to simulated spaceflight conditions. Moreover, astronaut samples from the NASA Twins Study confirmed these expression signatures in miRNA sequencing, single-cell RNA sequencing (scRNA-seq), and single-cell assay for transposase accessible chromatin (scATAC-seq) data. Additionally, a subset of these miRNAs (miR-125, miR-16, and let-7a) was found to regulate vascular damage caused by simulated deep space radiation. To demonstrate the physiological relevance of key spaceflight-associated miRNAs, we utilized antagomirs to inhibit their expression and successfully rescue simulated deep-space-radiation-mediated damage in human 3D vascular constructs.

mirnaQC: a webserver for comparative quality control of miRNA-seq data

Although miRNA-seq is extensively used in many different fields, its quality control is frequently restricted to a PhredScore-based filter. Other important quality related aspects like microRNA yield, the fraction of putative degradation products (such as rRNA fragments) or the percentage of adapter-dimers are hard to assess using absolute thresholds. Here we present mirnaQC, a webserver that relies on 34 quality parameters to assist in miRNA-seq quality control. To improve their interpretability, quality attributes are ranked using a reference distribution obtained from over 36 000 publicly available miRNA-seq datasets. Accepted input formats include FASTQ and SRA accessions. The results page contains several sections that deal with putative technical artefacts related to library preparation, sequencing, contamination or yield. Different visualisations, including PCA and heatmaps, are available to help users identify underlying issues. Finally, we show the usefulness of this approach by analysing two publicly available datasets and discussing the different quality issues that can be detected using mirnaQC.

sRNAbench and sRNAtoolbox 2019: intuitive fast small RNA profiling and differential expression

Since the original publication of sRNAtoolbox in 2015, small RNA research experienced notable advances in different directions. New protocols for small RNA sequencing have become available to address important issues such as adapter ligation bias, PCR amplification artefacts or to include internal controls such as spike-in sequences. New microRNA reference databases were developed with different foci, either prioritizing accuracy (low number of false positives) or completeness (low number of false negatives). Additionally, other small RNA molecules as well as microRNA sequence and length variants (isomiRs) have continued to gain importance. Finally, the number of microRNA sequencing studies deposited in GEO nearly triplicated from 2014 (280) to 2018 (764). These developments imply that fast and easy-to-use tools for expression profiling and subsequent downstream analysis of miRNA-seq data are essential to many researchers. Key features in this sRNAtoolbox release include addition of all major RNA library preparation protocols to sRNAbench and improvements in sRNAde, a tool that summarizes several aspects of small RNA sequencing studies including the detection of consensus differential expression. A special emphasis was put on the user-friendliness of the tools, for instance sRNAbench now supports parallel launching of several jobs to improve reproducibility and user time efficiency.

MirGeneDB 2.0: the metazoan microRNA complement

Small non-coding RNAs have gained substantial attention due to their roles in animal development and human disorders. Among them, microRNAs are special because individual gene sequences are conserved across the animal kingdom. In addition, unique and mechanistically well understood features can clearly distinguish bona fide miRNAs from the myriad other small RNAs generated by cells. However, making this distinction is not a common practice and, thus, not surprisingly, the heterogeneous quality of available miRNA complements has become a major concern in microRNA research. We addressed this by extensively expanding our curated microRNA gene database - MirGeneDB - to 45 organisms, encompassing a wide phylogenetic swath of animal evolution. By consistently annotating and naming 10,899 microRNA genes in these organisms, we show that previous microRNA annotations contained not only many false positives, but surprisingly lacked >2000 bona fide microRNAs. Indeed, curated microRNA complements of closely related organisms are very similar and can be used to reconstruct ancestral miRNA repertoires. MirGeneDB represents a robust platform for microRNA-based research, providing deeper and more significant insights into the biology and evolution of miRNAs as well as biomedical and biomarker research. MirGeneDB is publicly and freely available at http://mirgenedb.org/.

Disentangling sRNA-Seq data to study RNA communication between species

Many organisms exchange small RNAs (sRNAs) during their interactions, that can target or bolster defense strategies in host-pathogen systems. Current sRNA-Seq technology can determine the sRNAs present in any symbiotic system, but there are very few bioinformatic tools available to interpret the results. We show that one of the biggest challenges comes from sequences that map equally well to the genomes of both interacting organisms. This arises due to the small size of the sRNAs compared to large genomes, and because a large portion of sequenced sRNAs come from genomic regions that encode highly conserved miRNAs, rRNAs or tRNAs. Here, we present strategies to disentangle sRNA-Seq data from samples of communicating organisms, developed using diverse plant and animal species that are known to receive or exchange RNA with their symbionts. We show that sequence assembly, both de novo and genome-guided, can be used for these sRNA-Seq data, greatly reducing the ambiguity of mapping reads. Even confidently mapped sequences can be misleading, so we further demonstrate the use of differential expression strategies to determine true parasite-derived sRNAs within host cells. We validate our methods on new experiments designed to probe the nature of the extracellular vesicle sRNAs from the parasitic nematode Heligmosomoides bakeri that get into mouse intestinal epithelial cells.

Small noncoding RNA discovery and profiling with sRNAtools based on high-throughput sequencing

Small noncoding RNAs (sRNA/sncRNAs) are generated from different genomic loci and play important roles in biological processes, such as cell proliferation and the regulation of gene expression. Next-generation sequencing (NGS) has provided an unprecedented opportunity to discover and quantify diverse kinds of sncRNA, such as tRFs (tRNA-derived small RNA fragments), phasiRNAs (phased, secondary, small-interfering RNAs), Piwi-interacting RNA (piRNAs) and plant-specific 24-nt short interfering RNAs (siRNAs). However, currently available web-based tools do not provide approaches to comprehensively analyze all of these diverse sncRNAs. This study presents a novel integrated platform, sRNAtools (https://bioinformatics.caf.ac.cn/sRNAtools), that can be used in conjunction with high-throughput sequencing to identify and functionally annotate sncRNAs, including profiling microRNAss, piRNAs, tRNAs, small nuclear RNAs, small nucleolar RNAs and rRNAs and discovering isomiRs, tRFs, phasiRNAs and plant-specific 24-nt siRNAs for up to 21 model organisms. Different modules, including single case, batch case, group case and target case, are developed to provide users with flexible ways of studying sncRNA. In addition, sRNAtools supports different ways of uploading small RNA sequencing data in a very interactive queue system, while local versions based on the program package/Docker/virtureBox are also available. We believe that sRNAtools will greatly benefit the scientific community as an integrated tool for studying sncRNAs.

Evolutionary Implications of the microRNA- and piRNA Complement of Lepidodermella squamata (Gastrotricha)

Gastrotrichs-'hairy bellies'-are microscopic free-living animals inhabiting marine and freshwater habitats. Based on morphological and early molecular analyses, gastrotrichs were placed close to nematodes, but recent phylogenomic analyses have suggested their close relationship to flatworms (Platyhelminthes) within Spiralia. Small non-coding RNA data on e.g., microRNAs (miRNAs) and PIWI-interacting RNAs (piRNA) may help to resolve this long-standing question. MiRNAs are short post-transcriptional gene regulators that together with piRNAs play key roles in development. In a 'multi-omics' approach we here used small-RNA sequencing, available transcriptome and genomic data to unravel the miRNA- and piRNA complements along with the RNAi (RNA interference) protein machinery of Lepidodermella squamata (Gastrotricha, Chaetonotida). We identified 52 miRNA genes representing 35 highly conserved miRNA families specific to Eumetazoa, Bilateria, Protostomia, and Spiralia, respectively, with overall high similarities to platyhelminth miRNA complements. In addition, we found four large piRNA clusters that also resemble flatworm piRNAs but not those earlier described for nematodes. Congruently, transcriptomic annotation revealed that the Lepidodermella protein machinery is highly similar to flatworms, too. Taken together, miRNA, piRNA, and protein data support a close relationship of gastrotrichs and flatworms.