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

Towards understanding sex differences in autism spectrum disorders

Autism Spectrum Disorder (ASD) is a neurodevelopmental condition characterized by social deficits, repetitive behaviours and lack of empathy. Its significant genetic heritability and potential comorbidities often lead to diagnostic and therapeutic challenges. This review addresses the biological basis of ASD, focusing on the sex differences in gene expression and hormonal influences. ASD is more commonly diagnosed in males at a ratio of 4:1, indicating a potential oversight in female-specific ASD research and a risk of underdiagnosis in females. We consider how ASD manifests differently across sexes by exploring differential gene expression in female and male brains and consider how variations in steroid hormones influence ASD characteristics. Synaptic function, including excitation/inhibition ratio imbalance, is influenced by gene mutations and this is explored as a key factor in the cognitive and behavioural manifestations of ASD. We also discuss the role of micro RNAs (miRNAs) and highlight a novel mutation in miRNA-873, which affects a suite of key synaptic genes, neurexin, neuroligin, SHANK and post-synaptic density proteins, implicated in the pathology of ASD. Our review suggests that genetic predisposition, sex differences in brain gene expression, and hormonal factors significantly contribute to the presentation, identification and severity of ASD, necessitating sex-specific considerations in diagnosis and treatments. These findings advocate for personalized interventions to improve the outcomes for individuals with ASD.

Biosensor capability of the endometrium is mediated in part, by altered miRNA cargo from conceptus-derived extracellular vesicles

We tested the hypothesis that the biosensor capability of the endometrium is mediated in part, by the effect of different cargo contained in the extracellular vesicles secreted by the conceptus during the peri-implantation period of pregnancy. We transferred Bos taurus taurus embryos of different origin, in vivo (high developmental potential (IV)), in vitro (intermediate developmental potential (IVF)), or cloned (low developmental potential (NT)), into Bos taurus indicus recipients. Extracellular vesicles (EVs) recovered from Day 16 conceptus-conditioned medium were characterized and their microRNA (miRNA) cargo sequenced alongside RNA sequencing of their respective endometria. There were substantial differences in the endometrial response to in vivo versus in vitro and in vivo versus cloned conceptuses (1153 and 334DEGs respectively) with limited differences between in vitro Vs cloned conceptuses (36 DEGs). The miRNA cargo contained in conceptus-derived EVs was similar between all three groups (426 miRNA in common). Only 8 miRNAs were different between in vivo and cloned conceptuses, while only 6 miRNAs were different between in vivo and in vitro-derived conceptuses. Treatment of endometrial epithelial cells with mimic or inhibitors for miR-128 and miR-1298 changed the proteomic content of target cells (96 and 85, respectively) of which mRNAs are altered in the endometrium in vivo (PLXDC2, COPG1, HSPA12A, MCM5, TBL1XR1, and TTF). In conclusion, we have determined that the biosensor capability of the endometrium is mediated in part, by its response to different EVs miRNA cargo produced by the conceptus during the peri-implantation period of pregnancy.

Analysis of Amblyomma americanum microRNAs in response to Ehrlichia chaff eensis 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 chaf feensis, 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.

The Evolution and Characterization of the RNA Interference Pathways in Lophotrochozoa

In animals, three main RNA interference mechanisms have been described so far, which respectively maturate three types of small noncoding RNAs (sncRNAs): miRNAs, piRNAs, and endo-siRNAs. The diversification of these mechanisms is deeply linked with the evolution of the Argonaute gene superfamily since each type of sncRNA is typically loaded by a specific Argonaute homolog. Moreover, other protein families play pivotal roles in the maturation of sncRNAs, like the DICER ribonuclease family, whose DICER1 and DICER2 paralogs maturate respectively miRNAs and endo-siRNAs. Within Metazoa, the distribution of these families has been only studied in major groups, and there are very few data for clades like Lophotrochozoa. Thus, we here inferred the evolutionary history of the animal Argonaute and DICER families including 43 lophotrochozoan species. Phylogenetic analyses along with newly sequenced sncRNA libraries suggested that in all Trochozoa, the proteins related to the endo-siRNA pathway have been lost, a part of them in some phyla (i.e. Nemertea, Bryozoa, Entoprocta), while all of them in all the others. On the contrary, early diverging phyla, Platyhelminthes and Syndermata, showed a complete endo-siRNA pathway. On the other hand, miRNAs were revealed the most conserved and ubiquitous mechanism of the metazoan RNA interference machinery, confirming their pivotal role in animal cell regulation.

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.

Elucidating the ecophysiology of soybean pod-sucking stinkbug Riptortus pedestris (Hemiptera: Alydidae) based on de novo genome assembly and transcriptome analysis

Food security is important for the ever-growing global population. Soybean, Glycine max (L.) Merr., is cultivated worldwide providing a key source of food, protein and oil. Hence, it is imperative to maintain or to increase its yield under different conditions including challenges caused by abiotic and biotic stresses. In recent years, the soybean pod-sucking stinkbug Riptortus pedestris has emerged as an important agricultural insect pest in East, South and Southeast Asia. Here, we present a genomics resource for R. pedestris including its genome assembly, messenger RNA (mRNA) and microRNA (miRNA) transcriptomes at different developmental stages and from different organs. As insect hormone biosynthesis genes (genes involved in metamorphosis) and their regulators such as miRNAs are potential targets for pest control, we analyzed the sesquiterpenoid (juvenile) and ecdysteroid (molting) hormone biosynthesis pathway genes including their miRNAs and relevant neuropeptides. Temporal gene expression changes of these insect hormone biosynthesis pathways were observed at different developmental stages. Similarly, a diet-specific response in gene expression was also observed in both head and salivary glands. Furthermore, we observed that microRNAs (bantam, miR-14, miR-316, and miR-263) of R. pedestris fed with different types of soybeans were differentially expressed in the salivary glands indicating a diet-specific response. Interestingly, the opposite arms of miR-281 (-5p and -3p), a miRNA involved in regulating development, were predicted to target Hmgs genes of R. pedestris and soybean, respectively. These observations among others highlight stinkbug's responses as a function of its interaction with soybean. In brief, the results of this study not only present salient findings that could be of potential use in pest management and mitigation but also provide an invaluable resource for R. pedestris as an insect model to facilitate studies on plant-pest interactions.

Convergent and divergent evolution of microRNA-mediated regulation in metazoans

The evolution of microRNAs (miRNAs) has been studied extensively to understand their roles in gene regulation and evolutionary processes. This review focuses on how miRNA-mediated regulation has evolved in bilaterian animals, highlighting both convergent and divergent evolution. Since animals and plants display significant differences in miRNA biogenesis and target recognition, the 'independent origin' hypothesis proposes that miRNA pathways in these groups independently evolved from the RNA interference (RNAi) pathway, leading to modern miRNA repertoires through convergent evolution. However, recent evidence raises the alternative possibility that the miRNA pathway might have already existed in the last common ancestor of eukaryotes, and that the differences in miRNA pathway and miRNA repertoires among animal and plant lineages arise from lineage-specific innovations and losses of miRNA pathways, miRNA acquisition, and loss of miRNAs after eukaryotic divergence. The repertoire of miRNAs has considerably expanded during bilaterian evolution, primarily through de novo creation and duplication processes, generating new miRNAs. Although ancient functionally established miRNAs are rarely lost, many newly emerged miRNAs are transient and lineage specific, following a birth-death evolutionary pattern aligning with the 'out-of-the-testis' and 'transcriptional control' hypotheses. Our focus then shifts to the convergent molecular evolution of miRNAs. We summarize how miRNA clustering and seed mimicry contribute to this phenomenon, and we review how miRNAs from different sources converge to degrade maternal messenger RNAs (mRNAs) during animal development. Additionally, we describe how miRNAs evolve across species due to changes in sequence, seed shifting, arm switching, and spatiotemporal expression patterns, which can result in variations in target sites among orthologous miRNAs across distant strains or species. We also provide a summary of the current understanding regarding how the target sites of orthologous miRNAs can vary across strains or distantly related species. Although many paralogous miRNAs retain their seed or mature sequences after duplication, alterations can occur in the seed or mature sequences or expression patterns of paralogous miRNAs, leading to functional diversification. We discuss our current understanding of the functional divergence between duplicated miRNAs, and illustrate how the functional diversification of duplicated miRNAs impacts target site evolution. By investigating these topics, we aim to enhance our current understanding of the functions and evolutionary dynamics of miRNAs. Additionally, we shed light on the existing challenges in miRNA evolutionary studies, particularly the complexity of deciphering the role of miRNA-mediated regulatory network evolution in shaping gene expression divergence and phenotypic differences among species.

BioNexusSentinel: a visual tool for bioregulatory network and cytohistological RNA-seq genetic expression profiling within the context of multicellular simulation research using ChatGPT-augmented software engineering

Summary

Motivated by the need to parameterize ongoing multicellular simulation research, this paper documents the culmination of a ChatGPT augmented software engineering cycle resulting in an integrated visual platform for efficient cytohistological RNA-seq and bioregulatory network exploration. As contrasted to other systems and synthetic biology tools, BioNexusSentinel was developed de novo to uniquely combine these features. Reactome served as the primary source of remotely accessible biological models, accessible using BioNexusSentinel's novel search engine and REST API requests. The innovative, feature-rich gene expression profiler component was developed to enhance the exploratory experience for the researcher, culminating in the cytohistological RNA-seq explorer based on Human Protein Atlas data. A novel cytohistological classifier would be integrated via pre-processed analysis of the RNA-seq data via R statistical language, providing for useful analytical functionality and good performance for the end-user. Implications of the work span prospects for model orthogonality evaluations, gap identification in network modelling, prototyped automatic kinetics parameterization, and downstream simulation and cellular biological state analysis. This unique computational biology software engineering collaboration with generative natural language processing artificial intelligence was shown to enhance worker productivity, with evident benefits in terms of accelerating coding and machine-human intelligence transfer.

Availability and implementation

BioNexusSentinel project releases, with corresponding data and installation instructions, are available at https://github.com/RichardMatzko/BioNexusSentinel.

sRNA-Effector: A tool to expedite discovery of small RNA regulators

microRNAs (miRNAs) are small regulatory RNAs that repress target mRNA transcripts through base pairing. Although the mechanisms of miRNA production and function are clearly established, new insights into miRNA regulation or miRNA-mediated gene silencing are still emerging. In order to facilitate the discovery of miRNA regulators or effectors, we have developed sRNA-Effector, a machine learning algorithm trained on enhanced crosslinking and immunoprecipitation sequencing and RNA sequencing data following knockdown of specific genes. sRNA-Effector can accurately identify known miRNA biogenesis and effector proteins and identifies 9 putative regulators of miRNA function, including serine/threonine kinase STK33, splicing factor SFPQ, and proto-oncogene BMI1. We validated the role of STK33, SFPQ, and BMI1 in miRNA regulation, showing that sRNA-Effector is useful for identifying new players in small RNA biology. sRNA-Effector will be a web tool available for all researchers to identify potential miRNA regulators in any cell line of interest.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

Two-motif model illuminates DICER cleavage preferences

In humans, DICER is a key regulator of gene expression through its production of miRNAs and siRNAs by processing miRNA precursors (pre-miRNAs), short-hairpin RNAs (shRNAs), and long double-stranded RNAs (dsRNAs). To advance our understanding of this process, we employed high-throughput dicing assays using various shRNA variants and both wild-type and mutant DICER. Our analysis revealed that DICER predominantly cleaves shRNAs at two positions, specifically at 21 (DC21) and 22 (DC22) nucleotides from their 5'-end. Our investigation identified two different motifs, mWCU and YCR, that determine whether DICER cleaves at DC21 or DC22, depending on their locations in shRNAs/pre-miRNAs. These motifs can work together or independently to determine the cleavage sites of DICER. Furthermore, our findings indicate that dsRNA-binding domain (dsRBD) of DICER enhances its cleavage, and mWCU strengthens the interaction between dsRBD and RNA, leading to an even greater enhancement of the cleavage. Conversely, YCR functions independently of dsRBD. Our study proposes a two-motif model that sheds light on the intricate regulatory mechanisms involved in gene expression by elucidating how DICER recognizes its substrates, providing valuable insights into this critical biological process.

The miRNA-target interactions: An underestimated intricacy

MicroRNAs (miRNAs) play indispensable roles in posttranscriptional gene regulation. Their cellular regulatory impact is determined not solely by their sheer number, which likely amounts to >2000 individual miRNAs in human, than by the regulatory effectiveness of single miRNAs. Although, one begins to develop an understanding of the complex mechanisms underlying miRNA-target interactions (MTIs), the overall knowledge of MTI functionality is still rather patchy. In this critical review, we summarize key features of mammalian MTIs. We especially highlight latest insights on (i) the dynamic make-up of miRNA binding sites including non-canonical binding sites, (ii) the cooperativity between miRNA binding sites, (iii) the adaptivity of MTIs through sequence modifications, (iv) the bearing of intra-cellular miRNA localization changes and (v) the role of cell type and cell status specific miRNA interaction partners. The MTI biology is discussed against the background of state-of-the-art approaches with particular emphasis on experimental strategies for evaluating miRNA functionality.

Discovery of novel microRNA mimic repressors of ribosome biogenesis

While microRNAs and other non-coding RNAs are the next frontier of novel regulators of mammalian ribosome biogenesis (RB), a systematic exploration of microRNA-mediated RB regulation has not yet been undertaken. We carried out a high-content screen in MCF10A cells for changes in nucleolar number using a library of 2603 mature human microRNA mimics. Following a secondary screen for nucleolar rRNA biogenesis inhibition, we identified 72 novel microRNA negative regulators of RB after stringent hit calling. Hits included 27 well-conserved microRNAs present in MirGeneDB, and were enriched for mRNA targets encoding proteins with nucleolar localization or functions in cell cycle regulation. Rigorous selection and validation of a subset of 15 microRNA hits unexpectedly revealed that most of them caused dysregulated pre-rRNA processing, elucidating a novel role for microRNAs in RB regulation. Almost all hits impaired global protein synthesis and upregulated CDKN1A (p21) levels, while causing diverse effects on RNA Polymerase 1 (RNAP1) transcription and TP53 protein levels. We provide evidence that the MIR-28 siblings, hsa-miR-28-5p and hsa-miR-708-5p, potently target the ribosomal protein mRNA RPS28 via tandem primate-specific 3' UTR binding sites, causing a severe pre-18S pre-rRNA processing defect. Our work illuminates novel microRNA attenuators of RB, forging a promising new path for microRNA mimic chemotherapeutics.

Rapid adaptation of cellular metabolic rate to the MicroRNA complements of mammals and its relevance to the evolution of endothermy

The metabolic efficiency of mammalian cells depends on the attenuation of intrinsic translation noise by microRNAs. We devised a metric of cellular metabolic rate (cMR), rMR/Mexp optimally fit to the number of microRNA families (mirFam), that is robust to variation in mass and sensitive to body temperature (Tb), consistent with the heat dissipation limit theory of Speakman and Król (2010). Using mirFam as predictor, an Ornstein-Uhlenbeck process of stabilizing selection, with an adaptive shift at the divergence of Boreoeutheria, accounted for 95% of the variation in cMR across mammals. Branchwise rates of evolution of cMR, mirFam and Tb concurrently increased 6- to 7-fold at the divergence of Boreoeutheria, independent of mass. Cellular MR variation across placental mammals was also predicted by the sum of model conserved microRNA-target interactions, revealing an unexpected degree of integration of the microRNA-target apparatus into the energy economy of the mammalian cell.

Recent advances in understanding microRNA function and regulation in C. elegans

MicroRNAs (miRNAs) were first discovered in C. elegans as essential post-transcriptional regulators of gene expression. Since their initial discovery, miRNAs have been implicated in numerous areas of physiology and disease in all animals examined. In recent years, the C. elegans model continues to contribute important advances to all areas of miRNA research. Technological advances in tissue-specific miRNA profiling and genome editing have driven breakthroughs in understanding biological functions of miRNAs, mechanism of miRNA action, and regulation of miRNAs. In this review, we highlight these new C. elegans findings from the past five to seven years.

NCAD v1.0: a database for non-coding variant annotation and interpretation

The application of whole genome sequencing is expanding in clinical diagnostics across various genetic disorders, and the significance of non-coding variants in penetrant diseases is increasingly being demonstrated. Therefore, it is urgent to improve the diagnostic yield by exploring the pathogenic mechanisms of variants in non-coding regions. However, the interpretation of non-coding variants remains a significant challenge, due to the complex functional regulatory mechanisms of non-coding regions and the current limitations of available databases and tools. Hence, we develop the non-coding variant annotation database (NCAD, http://www.ncawdb.net/), encompassing comprehensive insights into 665,679,194 variants, regulatory elements, and element interaction details. Integrating data from 96 sources, spanning both GRCh37 and GRCh38 versions, NCAD v1.0 provides vital information to support the genetic diagnosis of non-coding variants, including allele frequencies of 12 diverse populations, with a particular focus on the population frequency information for 230,235,698 variants in 20,964 Chinese individuals. Moreover, it offers prediction scores for variant functionality, five categories of regulatory elements, and four types of non-coding RNAs. With its rich data and comprehensive coverage, NCAD serves as a valuable platform, empowering researchers and clinicians with profound insights into non-coding regulatory mechanisms while facilitating the interpretation of non-coding variants.

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.

Differential expression of miRNAs associated with pectoral myopathies in young broilers: insights from a comparative transcriptome analysis

INTRODUCTION

White Striping (WS) and Wooden Breast (WB) pectoral myopathies are relevant disorders for contemporary broiler production worldwide. Several studies aimed to elucidate the genetic components associated with the occurrence of these myopathies. However, epigenetic factors that trigger or differentiate these two conditions are still unclear. The aim of this study was to identify miRNAs differentially expressed (DE) between normal and WS and WB-affected broilers, and to verify the possible role of these miRNAs in metabolic pathways related to the manifestation of these pectoral myopathies in 28-day-old broilers.

RESULTS

Five miRNAs were DE in the WS vs control (gga-miR-375, gga-miR-200b-3p, gga-miR-429-3p, gga-miR-1769-5p, gga-miR-200a-3p), 82 between WB vs control and 62 between WB vs WS. Several known miRNAs were associated with WB, such as gga-miR-155, gga-miR-146b, gga-miR-222, gga-miR-146-5p, gga-miR- 29, gga-miR-21-5p, gga-miR-133a-3p and gga-miR-133b. Most of them had not previously been associated with the development of this myopathy in broilers. We also have predicted 17 new miRNAs expressed in the broilers pectoral muscle. DE miRNA target gene ontology analysis enriched 6 common pathways for WS and WB compared to control: autophagy, insulin signaling, FoxO signaling, endocytosis, and metabolic pathways. The WS vs control contrast had two unique pathways, ERBB signaling and the mTOR signaling, while WB vs control had 14 unique pathways, with ubiquitin-mediated proteolysis and endoplasmic reticulum protein processing being the most significant.

CONCLUSIONS

We found miRNAs DE between normal broilers and those affected with breast myopathies at 28 days of age. Our results also provide novel evidence of the miRNAs role on the regulation of WS and in the differentiation of both WS and WB myopathies. Overall, our study provides insights into miRNA-mediated and pathways involved in the occurrence of WS and WB helping to better understand these chicken growth disorders in an early age. These findings can help developing new approaches to reduce these complex issues in poultry production possibly by adjustments in nutrition and management conditions. Moreover, the miRNAs and target genes associated with the initial stages of WS and WB development could be potential biomarkers to be used in selection to reduce the occurrence of these myopathies in broiler production.

Hepatitis B Virus and microRNAs: A Bioinformatics Approach

In recent decades, microRNAs (miRNAs) have emerged as key regulators of gene expression, and the identification of viral miRNAs (v-miRNAs) within some viruses, including hepatitis B virus (HBV), has attracted significant attention. HBV infections often progress to chronic states (CHB) and may induce fibrosis/cirrhosis and hepatocellular carcinoma (HCC). The presence of HBV can dysregulate host miRNA expression, influencing several biological pathways, such as apoptosis, innate and immune response, viral replication, and pathogenesis. Consequently, miRNAs are considered a promising biomarker for diagnostic, prognostic, and treatment response. The dynamics of miRNAs during HBV infection are multifaceted, influenced by host variability and miRNA interactions. Given the ability of miRNAs to target multiple messenger RNA (mRNA), understanding the viral-host (human) interplay is complex but essential to develop novel clinical applications. Therefore, bioinformatics can help to analyze, identify, and interpret a vast amount of miRNA data. This review explores the bioinformatics tools available for viral and host miRNA research. Moreover, we introduce a brief overview focusing on the role of miRNAs during HBV infection. In this way, this review aims to help the selection of the most appropriate bioinformatics tools based on requirements and research goals.

Generation of de novo miRNAs from template switching during DNA replication

The mechanisms generating novel genes and genetic information are poorly known, even for microRNA (miRNA) genes with an extremely constrained design. All miRNA primary transcripts need to fold into a stem-loop structure to yield short gene products ([Formula: see text]22 nt) that bind and repress their mRNA targets. While a substantial number of miRNA genes are ancient and highly conserved, short secondary structures coding for entirely novel miRNA genes have been shown to emerge in a lineage-specific manner. Template switching is a DNA-replication-related mutation mechanism that can introduce complex changes and generate perfect base pairing for entire hairpin structures in a single event. Here, we show that the template-switching mutations (TSMs) have participated in the emergence of over 6,000 suitable hairpin structures in the primate lineage to yield at least 18 new human miRNA genes, that is 26% of the miRNAs inferred to have arisen since the origin of primates. While the mechanism appears random, the TSM-generated miRNAs are enriched in introns where they can be expressed with their host genes. The high frequency of TSM events provides raw material for evolution. Being orders of magnitude faster than other mechanisms proposed for de novo creation of genes, TSM-generated miRNAs enable near-instant rewiring of genetic information and rapid adaptation to changing environments.

Beyond ribosome biogenesis: noncoding nucleolar RNAs in physiology and tumor biology

The nucleolus, the largest subcompartment of the nucleus, stands out from the nucleoplasm due to its exceptionally high local RNA and low DNA concentrations. Within this central hub of nuclear RNA metabolism, ribosome biogenesis is the most prominent ribonucleoprotein (RNP) biogenesis process, critically determining the structure and function of the nucleolus. However, recent studies have shed light on other roles of the nucleolus, exploring the interplay with various noncoding RNAs that are not directly involved in ribosome synthesis. This review focuses on this intriguing topic and summarizes the techniques to study and the latest findings on nucleolar long noncoding RNAs (lncRNAs) as well as microRNAs (miRNAs) in the context of nucleolus biology beyond ribosome biogenesis. We particularly focus on the multifaceted roles of the nucleolus and noncoding RNAs in physiology and tumor biology.

Quadruplexes and aging: G4-binding proteins regulate the presence of miRNA in small extracellular vesicles (sEVs)

The interaction between proteins and nucleic acids is a core element of life. Many proteins bind nucleic acids via a sequence-specific manner, but there are also many types of proteins that recognize various structural motifs. Researchers have recently found that proteins that can recognize DNA and RNA G-quadruplexes (G4s) are very important for basic cellular processes, particularly in eukaryotes. Some of these proteins are located outside the nucleus and interact with RNA, potentially affecting miRNA functions in intercellular communication, which is facilitated by small extracellular vesicles (sEVs). Imbalances in the production of sEVs are associated with various pathologies and senescence in humans. The distribution of miRNA into sEVs is regulated by two RNA-binding proteins, Alyref and FUS. Both proteins possess G-rich recognition motifs that are compatible with the formation of RNA parallel G4 structures. This lends credence to the new hypothesis that G4-formation in RNAs and their interaction with G4-binding proteins can affect the fate of miRNAs and control their distribution in sEVs that are associated with senescence and aging.

Discovery of novel microRNA mimic repressors of ribosome biogenesis

While microRNAs and other non-coding RNAs are the next frontier of novel regulators of mammalian ribosome biogenesis (RB), a systematic exploration of microRNA-mediated RB regulation has not yet been undertaken. We carried out a high-content screen in MCF10A cells for changes in nucleolar number using a library of 2,603 mature human microRNA mimics. Following a secondary screen for nucleolar rRNA biogenesis inhibition, we identified 72 novel microRNA negative regulators of RB after stringent hit calling. Hits included 27 well-conserved microRNAs present in MirGeneDB, and were enriched for mRNA targets encoding proteins with nucleolar localization or functions in cell cycle regulation. Rigorous selection and validation of a subset of 15 microRNA hits unexpectedly revealed that most of them caused dysregulated pre-rRNA processing, elucidating a novel role for microRNAs in RB regulation. Almost all hits impaired global protein synthesis and upregulated CDKN1A ( p21 ) levels, while causing diverse effects on RNA Polymerase 1 (RNAP1) transcription and TP53 protein levels. We discovered that the MIR-28 siblings, hsa-miR-28-5p and hsa-miR-708-5p, directly and potently target the ribosomal protein mRNA RPS28 via tandem primate-specific 3' UTR binding sites, causing a severe pre-18S pre-rRNA processing defect. Our work illuminates novel microRNA attenuators of RB, forging a promising new path for microRNA mimic chemotherapeutics.

Detection of transcriptome-wide microRNA-target interactions in single cells with agoTRIBE

MicroRNAs (miRNAs) exert their gene regulatory effects on numerous biological processes based on their selection of target transcripts. Current experimental methods available to identify miRNA targets are laborious and require millions of cells. Here we have overcome these limitations by fusing the miRNA effector protein Argonaute2 to the RNA editing domain of ADAR2, allowing the detection of miRNA targets transcriptome-wide in single cells. miRNAs guide the fusion protein to their natural target transcripts, causing them to undergo A>I editing, which can be detected by sensitive single-cell RNA sequencing. We show that agoTRIBE identifies functional miRNA targets, which are supported by evolutionary sequence conservation. In one application of the method we study microRNA interactions in single cells and identify substantial differential targeting across the cell cycle. AgoTRIBE also provides transcriptome-wide measurements of RNA abundance and allows the deconvolution of miRNA targeting in complex tissues at the single-cell level.

Virus-Derived Small RNAs and microRNAs in Health and Disease

MicroRNAs (miRNAs) are short noncoding RNAs that can regulate all steps of gene expression (induction, transcription, and translation). Several virus families, primarily double-stranded DNA viruses, encode small RNAs (sRNAs), including miRNAs. These virus-derived miRNAs (v-miRNAs) help the virus evade the host's innate and adaptive immune system and maintain an environment of chronic latent infection. In this review, the functions of the sRNA-mediated virus-host interactions are highlighted, delineating their implication in chronic stress, inflammation, immunopathology, and disease. We provide insights into the latest viral RNA-based research-in silico approaches for functional characterization of v-miRNAs and other RNA types. The latest research can assist toward the identification of therapeutic targets to combat viral infections.

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.

High-quality genome assemblies provide clues on the evolutionary advantage of blue peafowl over green peafowl

An intriguing example of differential adaptability is the case of two Asian peafowl species, Pavo cristatus (blue peafowl) and Pavo muticus (green peafowl), where the former has a "Least Concern" conservation status and the latter is an "Endangered" species. To understand the genetic basis of this differential adaptability of the two peafowl species, a comparative analysis of these species is much needed to gain the genomic and evolutionary insights. Thus, we constructed a high-quality genome assembly of blue peafowl with an N50 value of 84.81 Mb (pseudochromosome-level assembly), and a high-confidence coding gene set to perform the genomic and evolutionary analyses of blue and green peafowls with 49 other avian species. The analyses revealed adaptive evolution of genes related to neuronal development, immunity, and skeletal muscle development in these peafowl species. Major genes related to axon guidance such as NEO1 and UNC5, semaphorin (SEMA), and ephrin receptor showed adaptive evolution in peafowl species. However, blue peafowl showed the presence of 42% more coding genes compared to the green peafowl along with a higher number of species-specific gene clusters, segmental duplicated genes and expanded gene families, and comparatively higher evolution in neuronal and developmental pathways. Blue peafowl also showed longer branch length compared to green peafowl in the species phylogenetic tree. These genomic insights obtained from the high-quality genome assembly of P. cristatus constructed in this study provide new clues on the superior adaptability of the blue peafowl over green peafowl despite having a recent species divergence time.

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.

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

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

Sites of transcription initiation drive mRNA isoform selection

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

MicroRNA profiling of cerebrospinal fluid from dogs with steroid responsive meningitis-arteritis and meningoencephalitis of unknown origin

Introduction

Non-infectious inflammatory diseases of the central nervous system in dogs, such as steroid responsive meningitis-arteritis (SRMA) and meningoencephalitis of unknown origin (MUO), represent a common clinical challenge that needs extensive and multimodal work-up to reach a presumptive diagnosis. Both diseases are presumably caused by dysregulations of the immune system, but further research is needed in order to understand the molecular mechanisms behind each disease and to optimize treatment.

Methods

By next-generation sequencing and subsequent quantitative real-time PCR (qPCR) verification, we designed a prospective case-control pilot study to analyze the small RNA profiles of cerebrospinal fluid from dogs suffering from MUO (N = 5), dogs suffering from SRMA (N = 8), and healthy dogs (N = 5) presented for elective euthanasia used as the Control group.

Results

Our results showed an overall enrichment in Y-RNA fragments across all samples, followed by microRNAs (miRNAs) and ribosomal RNAs as the major findings. Additional traces of short RNA reads mapped to long non-coding RNAs and protein-coding genes were also found. From the detected canine miRNAs, miR-21, miR-486, miR-148a, miR-99a, miR-191 and miR-92a were among the most abundant. Dogs with SRMA showed higher differences in miRNA abundance than dogs with MUO when compared to healthy dogs, and miR-142-3p was consistently detected as differentially upregulated in both diseases, although at a low concentration. Moreover, miR-405-5p and miR-503-5p showed different profiles between SRMA and MUO dogs. Subsequent qPCR analyses confirmed miR-142-5p, miR-191-5p and miR-92a-3p as significantly upregulated miRNAs in dogs with SRMA and/or MUO.

Discussion

Cerebrospinal fluid is a challenging biological material to use for profiling miRNAs due to the low content of circulating RNAs. Despite this, we could confirm several miRNAs being differentially abundant when comparing healthy dogs and dogs with MUO and SRMA, respectively. The results of this study indicate a potential role of miRNAs in the underlying molecular mechanisms of these diseases and establish the basis for further studies.

Direct observation of the evolution of cell-type-specific microRNA expression signatures supports the hematopoietic origin model of endothelial cells

The evolution of specialized cell-types is a long-standing interest of biologists, but given the deep time-scales very difficult to reconstruct or observe. microRNAs have been linked to the evolution of cellular complexity and may inform on specialization. The endothelium is a vertebrate-specific specialization of the circulatory system that enabled a critical new level of vasoregulation. The evolutionary origin of these endothelial cells is unclear. We hypothesized that Mir-126, an endothelial cell-specific microRNA may be informative. We here reconstruct the evolutionary history of Mir-126. Mir-126 likely appeared in the last common ancestor of vertebrates and tunicates, which was a species without an endothelium, within an intron of the evolutionary much older EGF Like Domain Multiple (Egfl) locus. Mir-126 has a complex evolutionary history due to duplications and losses of both the host gene and the microRNA. Taking advantage of the strong evolutionary conservation of the microRNA among Olfactores, and using RNA in situ hybridization, we localized Mir-126 in the tunicate Ciona robusta. We found exclusive expression of the mature Mir-126 in granular amebocytes, supporting a long-proposed scenario that endothelial cells arose from hemoblasts, a type of proto-endothelial amoebocyte found throughout invertebrates. This observed change of expression of Mir-126 from proto-endothelial amoebocytes in the tunicate to endothelial cells in vertebrates is the first direct observation of the evolution of a cell-type in relation to microRNA expression indicating that microRNAs can be a prerequisite of cell-type evolution.

DIANA-microT 2023: including predicted targets of virally encoded miRNAs

DIANA-microT-CDS is a state-of-the-art miRNA target prediction algorithm catering the scientific community since 2009. It is one of the first algorithms to predict miRNA binding sites in both the 3' Untranslated Region (3'-UTR) and the coding sequence (CDS) of transcripts, with increased performance. Its current version, DIANA-microT 2023 (www.microrna.gr/microt_webserver/), brings forward a significantly updated set of interactions. DIANA-microT-CDS has been executed utilizing annotation information from Ensembl v102, miRBase 22.1 and, for the first time, MirGeneDB 2.1, yielding more than 83 million interactions in human, mouse, rat, chicken, fly and worm species. Additionally, this version delivers predicted interactions of miRNAs encoded from 20 viruses against host transcripts from human, mouse and chicken species. Numerous resources have been interconnected into DIANA-microT, including DIANA-TarBase, plasmiR, HMDD, UCSC, dbSNP, ClinVar, as well as miRNA/gene abundance values for 369 distinct cell-lines/tissues. The server interface has been redesigned allowing users to use smart filtering options, identify abundance patterns of interest, pinpoint known SNPs residing on binding sites and obtain miRNA-disease information. The contents of DIANA-microT webserver are freely accessible and can also be locally downloaded without any login requirements.

Small RNAs in Cnidaria: A review

As fundamental components of RNA silencing, small RNA (sRNA) molecules ranging from 20 to 32 nucleotides in length have been found as potent regulators of gene expression and genome stability in many biological processes of eukaryotes. Three major small RNAs are active in animals, including the microRNA (miRNA), short interfering RNA (siRNA), and PIWI-interacting RNA (piRNA). Cnidarians, the sister group to bilaterians, are at a critical phylogenetic node to better model eukaryotic small RNA pathway evolution. To date, most of our understanding of sRNA regulation and its potential contribution to evolution has been limited to a few triploblastic bilaterian and plant models. The diploblastic nonbilaterians, including the cnidarians, are understudied in this regard. Therefore, this review will present the current-known small RNA information in cnidarians to enhance our understanding of the development of the small RNA pathways in early branch animals.

Shared and organ-specific gene-expression programs during the development of the cochlea and the superior olivary complex

The peripheral and central auditory subsystems together form a complex sensory network that allows an organism to hear. The genetic programs of the two subsystems must therefore be tightly coordinated during development. Yet, their interactions and common expression pathways have never been systematically explored. MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene expression and are essential for normal development of the auditory system. We performed mRNA and small-RNA sequencing of organs from both auditory subsystems at three critical developmental timepoints (E16, P0, P16) to obtain a comprehensive and unbiased insight of their expression profiles. Our analysis reveals common and organ-specific expression patterns for differentially regulated mRNAs and miRNAs, which could be clustered with a particular selection of functions such as inner ear development, Wnt signalling, K+ transport, and axon guidance, based on gene ontology. Bioinformatics detected enrichment of predicted targets of specific miRNAs in the clusters and predicted regulatory interactions by monitoring opposite trends of expression of miRNAs and their targets. This approach identified six miRNAs as strong regulatory candidates for both subsystems. Among them was miR-96, an established critical factor for proper development in both subsystems, demonstrating the strength of our approach. We suggest that other miRNAs identified by this analysis are also common effectors of proper hearing acquirement. This first combined comprehensive analysis of the developmental program of the peripheral and central auditory systems provides important data and bioinformatics insights into the shared genetic program of the two sensory subsystems and their regulation by miRNAs.

Dynamic and static circulating cancer microRNA biomarkers - a validation study

For cancers and other pathologies, early diagnosis remains the most promising path to survival. Profiling of longitudinal cohorts facilitates insights into trajectories of biomarkers. We measured microRNA expression in 240 serum samples from patients with colon, lung, and breast cancer and from cancer-free controls. Each patient provided at least two serum samples, one prior to diagnosis and one following diagnosis. The median time interval between the samples was 11.6 years. Using computational models, we evaluated the circulating profiles of 21 microRNAs. The analysis yielded two sets of biomarkers, static ones that show an absolute difference between certain cancer types and controls and dynamic ones where the level over time provided higher diagnostic information content. In the first group, miR-99a-5p stands out for all three cancer types. In the second group, miR-155-5p allows to predict lung cancers and colon cancers. Classification in samples from cancer and non-cancer patients using gradient boosted trees reached an average accuracy of 79.9%. The results suggest that individual change over time or an absolute value at one time point may predict a disease with high specificity and sensitivity.

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.

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.

Systematic characterization of seed overlap microRNA cotargeting associated with lupus pathogenesis

BACKGROUND

Combinatorial gene regulation by multiple microRNAs (miRNAs) is widespread and closely spaced target sites often act cooperatively to achieve stronger repression ("neighborhood" miRNA cotargeting). While miRNA cotarget sites are suggested to be more conserved and implicated in developmental control, the pathological significance of miRNA cotargeting remains elusive.

RESULTS

Here, we report the pathogenic impacts of combinatorial miRNA regulation on inflammation in systemic lupus erythematosus (SLE). In the SLE mouse model, we identified the downregulation of two miRNAs, miR-128 and miR-148a, by TLR7 stimulation in plasmacytoid dendritic cells. Functional analyses using human cell lines demonstrated that miR-128 and miR-148a additively target KLF4 via extensively overlapping target sites ("seed overlap" miRNA cotargeting) and suppress the inflammatory responses. At the transcriptome level, "seed overlap" miRNA cotargeting increases susceptibility to downregulation by two miRNAs, consistent with additive but not cooperative recruitment of two miRNAs. Systematic characterization further revealed that extensive "seed overlap" is a prevalent feature among broadly conserved miRNAs. Highly conserved target sites of broadly conserved miRNAs are largely divided into two classes-those conserved among eutherian mammals and from human to Coelacanth, and the latter, including KLF4-cotargeting sites, has a stronger association with both "seed overlap" and "neighborhood" miRNA cotargeting. Furthermore, a deeply conserved miRNA target class has a higher probability of haplo-insufficient genes.

CONCLUSIONS

Our study collectively suggests the complexity of distinct modes of miRNA cotargeting and the importance of their perturbations in human diseases.

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.

miRSCAPE - inferring miRNA expression from scRNA-seq data

Our understanding of miRNA activity at cellular resolution is thwarted by the inability of standard scRNA-seq protocols to capture miRNAs. We introduce a novel tool, miRSCAPE, to infer miRNA expression in a sample from its RNA-seq profile. We establish miRSCAPE’s accuracy in 10 tumor and normal cohorts demonstrating its superiority over alternatives. miRSCAPE accurately infers cell type-specific miRNA activities (predicted versus observed fold-difference correlation ∼0.81) in two independent scRNA-seq datasets. We apply miRSCAPE to infer miRNA activities in scRNA clusters in pancreatic and lung adenocarcinomas, as well as in 56 cell types in the human cell landscape (HCL). In pancreatic and breast cancer scRNA-seq data, miRSCAPE recapitulates miRNAs associated with stemness and epithelial-mesenchymal transition (EMT) cell states, respectively. Overall, miRSCAPE recapitulates and refines miRNA biology at cellular resolution. miRSCAPE is freely available and is easily applicable to scRNA-seq data to infer miRNA activities at cellular resolution.

•

Novel machine learning-based tool to infer miRNA expression at single cell level

•

Predicts miRNA activity with high accuracy in various contexts

•

Recaps miRNAs associated with specific cellular states and suggests novel candidates

•

Provides a general framework to predict other types of molecular data at a single cell

Bioinformatic Analysis of Ixodes ricinus Long Non-Coding RNAs Predicts Their Binding Ability of Host miRNAs

Ixodes ricinus ticks are distributed across Europe and are a vector of tick-borne diseases. Although I. ricinus transcriptome studies have focused exclusively on protein coding genes, the last decade witnessed a strong increase in long non-coding RNA (lncRNA) research and characterization. Here, we report for the first time an exhaustive analysis of these non-coding molecules in I. ricinus based on 131 RNA-seq datasets from three different BioProjects. Using this data, we obtained a consensus set of lncRNAs and showed that lncRNA expression is stable among different studies. While the length distribution of lncRNAs from the individual data sets is biased toward short length values, implying the existence of technical artefacts, the consensus lncRNAs show a more homogeneous distribution emphasizing the importance to incorporate data from different sources to generate a solid reference set of lncRNAs. KEGG enrichment analysis of host miRNAs putatively targeting lncRNAs upregulated upon feeding showed that these miRNAs are involved in several relevant functions for the tick-host interaction. The possibility that at least some tick lncRNAs act as host miRNA sponges was further explored by identifying lncRNAs with many target regions for a given host miRNA or sets of host miRNAs that consistently target lncRNAs together. Overall, our findings suggest that lncRNAs that may act as sponges have diverse biological roles related to the tick–host interaction in different tissues.

A curated human cellular microRNAome based on 196 primary cell types

Background

An incomplete picture of the expression distribution of microRNAs (miRNAs) across human cell types has long hindered our understanding of this important regulatory class of RNA. With the continued increase in available public small RNA sequencing datasets, there is an opportunity to more fully understand the general distribution of miRNAs at the cell level.

Results

From the NCBI Sequence Read Archive, we obtained 6,054 human primary cell datasets and processed 4,184 of them through the miRge3.0 small RNA sequencing alignment software. This dataset was curated down, through shared miRNA expression patterns, to 2,077 samples from 196 unique cell types derived from 175 separate studies. Of 2,731 putative miRNAs listed in miRBase (v22.1), 2,452 (89.8%) were detected. Among reasonably expressed miRNAs, 108 were designated as cell specific/near specific, 59 as infrequent, 52 as frequent, 54 as near ubiquitous, and 50 as ubiquitous. The complexity of cellular microRNA expression estimates recapitulates tissue expression patterns and informs on the miRNA composition of plasma.

Conclusions

This study represents the most complete reference, to date, of miRNA expression patterns by primary cell type. The data are available through the human cellular microRNAome track at the UCSC Genome Browser (https://genome.ucsc.edu/cgi-bin/hgHubConnect) and an R/Bioconductor package (https://bioconductor.org/packages/microRNAome/).

Urinary microRNAome in healthy cats and cats with pyelonephritis or other urological conditions

MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene expression at the post-transcriptional level. miRNAs have been found in urine and have shown diagnostic potential in human nephropathies. Here, we aimed to characterize, for the first time, the feline urinary miRNAome and explore the use of urinary miRNA profiles as non-invasive biomarkers for feline pyelonephritis (PN). Thirty-eight cats were included in a prospective case-control study and classified in five groups: healthy Control cats (n = 11), cats with PN (n = 10), cats with subclinical bacteriuria or cystitis (SB/C, n = 5), cats with ureteral obstruction (n = 7) and cats with chronic kidney disease (n = 5). By small RNA sequencing we identified 212 miRNAs in cat urine, including annotated (n = 137) and putative novel (n = 75) miRNAs. The 15 most highly abundant urinary miRNAs accounted for nearly 71% of all detected miRNAs, most of which were previously identified in feline kidney. Ninety-nine differentially abundant (DA) miRNAs were identified when comparing Control cats to cats with urological conditions and 102 DA miRNAs when comparing PN to other urological conditions. Tissue clustering analysis revealed that the majority of urine samples clustered close to kidney, which confirm the likely cellular origin of the secreted urinary miRNAs. Relevant DA miRNAs were verified by quantitative real-time PCR (qPCR). Eighteen miRNAs discriminated Control cats from cats with a urological condition. Of those, seven miRNAs were DA by both RNAseq and qPCR methods between Control and PN cats (miR-125b-5p, miR-27a-3p, miR-21-5p, miR-27b-3p, miR-125a-5p, miR-17-5p and miR-23a-3p) or DA between Control and SB/C cats (miR-125b-5p). Six additional miRNAs (miR-30b-5p, miR-30c, miR-30e-5p, miR-27a-3p, miR-27b-39 and miR-222) relevant for discriminating PN from other urological conditions were identified by qPCR alone (n = 4) or by both methods (n = 2) (P<0.05). This panel of 13 miRNAs has potential as non-invasive urinary biomarkers for diagnostic of PN and other urological conditions in cats.

A deep learning method for miRNA/isomiR target detection

Accurate identification of microRNA (miRNA) targets at base-pair resolution has been an open problem for over a decade. The recent discovery of miRNA isoforms (isomiRs) adds more complexity to this problem. Despite the existence of many methods, none considers isomiRs, and their performance is still suboptimal. We hypothesize that by taking the isomiR–mRNA interactions into account and applying a deep learning model to study miRNA–mRNA interaction features, we may improve the accuracy of miRNA target predictions. We developed a deep learning tool called DMISO to capture the intricate features of miRNA/isomiR–mRNA interactions. Based on tenfold cross-validation, DMISO showed high precision (95%) and recall (90%). Evaluated on three independent datasets, DMISO had superior performance to five tools, including three popular conventional tools and two recently developed deep learning-based tools. By applying two popular feature interpretation strategies, we demonstrated the importance of the miRNA regions other than their seeds and the potential contribution of the RNA-binding motifs within miRNAs/isomiRs and mRNAs to the miRNA/isomiR–mRNA interactions.

The limits of human microRNA annotation have been met

Over the last few years, the number of microRNAs in the human genome has become a controversially debated issue. Several publications reported thousands of putative novel microRNAs not included in the curated microRNA gene database MirGeneDB and the repository miRBase. Recently, by using sequencing of ∼300 human tissues and cell lines, the human RNA atlas, an expanded inventory of human RNA annotations, was published, reporting thousands of putative microRNAs. We, the developers of established microRNA prediction tools and hosts of MirGeneDB, raise concerns about the frequently applied prediction and functional validation strategies, briefly discussing the drawbacks of false positive detections. By means of quantifying well-established biogenesis-derived features, we show that the reported novel microRNAs essentially represent false-positives and argue that the human microRNA complement, at about 550 microRNA genes, is already near complete. Output of available tools must be curated as false predictions will misguide scientists looking for biomarkers or therapeutic targets.