mirPRo-a novel standalone program for differential expression and variation analysis of miRNAs

From follicle to blastocyst: microRNA-34c from follicular fluid-derived extracellular vesicles modulates blastocyst quality

BACKGROUND

Within the follicular fluid, extracellular vesicles (EVs) guide oocyte growth through their cargo microRNAs (miRNAs). Here, we investigated the role of EVs and their cargo miRNAs by linking the miRNAs found in EVs, derived from the fluid of an individual follicle, to the ability of its oocyte to become a blastocyst (competent) or not (non-competent).

METHODS

Bovine antral follicles were dissected, categorized as small (2-4 mm) or large (5-8 mm) and the corresponding oocytes were subjected to individual maturation, fertilization and embryo culture to the blastocyst stage. Follicular fluid was pooled in 4 groups (4 replicates) based on follicle size and competence of the corresponding oocyte to produce a blastocyst. Follicular fluid-derived EVs were isolated, characterized, and subjected to miRNA-sequencing (Illumina Miseq) to assess differential expression (DE) in the 4 groups. Functional validation of the effect of miR-34c on embryo development was performed by supplementation of mimics and inhibitors during in vitro maturation (IVM).

RESULTS

We identified 16 DE miRNAs linked to oocyte competence when follicular size was not considered. Within the large and small follicles, 46 DE miRNAs were driving blastocyst formation in each group. Comparison of EVs from competent small and large follicles revealed 90 DE miRNAs. Cell regulation, cell differentiation, cell cycle, and metabolic process regulation were the most enriched pathways targeted by the DE miRNAs from competent oocytes. We identified bta-miR-34c as the most abundant in follicular fluid containing competent oocytes. Supplementation of miR-34c mimic and inhibitor during IVM did not affect embryo development. However, blastocyst quality, as evidenced by higher cell numbers, was significantly improved following oocyte IVM in the presence of miR-34c mimics, while miR-34c inhibitors resulted in the opposite effect.

CONCLUSION

This study demonstrates the regulatory effect of miRNAs from follicular fluid-derived EVs on oocyte competence acquisition, providing a further basis for understanding the significance of miRNAs in oocyte maturation and embryonic development. Up-regulation of miR-34c in EVs from follicular fluid containing competent oocytes and the positive impact of miR-34c mimics added during IVM on the resulting blastocysts indicate its pivotal role in oocyte competence.

Expression profile and gap-junctional transfer of microRNAs in the bovine cumulus-oocyte complex

MicroRNAs (miRNA) are important regulators of oocyte maturation, playing a key role in modulating gene expression both in a temporal- and spatial-specific manner. These small non-coding RNAs are involved in important processes during oocyte maturation, acting as messengers between the oocyte and its surrounding cumulus cells. Despite its significance, the bidirectional communication mechanism is still unknown. To test miRNA communication between oocyte and surrounding cumulus cells through the gap junctions the gap junctions were either blocked with carbenoxolone or not. MiRNA sequencing of oocytes at 1, 6, and 22 h of in vitro maturation was then performed. Among the differentially expressed miRNAs, bta-miR-21-5p, a regulator of cumulus cell viability and oocyte maturation, was the only previously known miRNA. Furthermore, by labeling a bta-miR-21-5p mimic with FAM, crossing of this miRNA through the gap junctions within the cumulus-oocyte complex could be visualized and internalization in the oocyte was confirmed by RT-qPCR. In conclusion, this study provides, for the first time, evidence that miRNA communication within the bovine cumulus-oocyte complex is enabled through the gap junctional network.

A Unified Computational Framework for a Robust, Reliable, and Reproducible Identification of Novel miRNAs From the RNA Sequencing Data

In eukaryotic cells, miRNAs regulate a plethora of cellular functionalities ranging from cellular metabolisms, and development to the regulation of biological networks and pathways, both under homeostatic and pathological states like cancer.Despite their immense importance as key regulators of cellular processes, accurate and reliable estimation of miRNAs using Next Generation Sequencing is challenging, largely due to the limited availability of robust computational tools/methods/pipelines. Here, we introduce miRPipe, an end-to-end computational framework for the identification, characterization, and expression estimation of small RNAs, including the known and novel miRNAs and previously annotated pi-RNAs from small-RNA sequencing profiles. Our workflow detects unique novel miRNAs by incorporating the sequence information of seed and non-seed regions, concomitant with clustering analysis. This approach allows reliable and reproducible detection of unique novel miRNAs and functionally same miRNAs (paralogues). We validated the performance of miRPipe with the available state-of-the-art pipelines using both synthetic datasets generated using the newly developed miRSim tool and three cancer datasets (Chronic Lymphocytic Leukemia, Lung cancer, and breast cancer). In the experiment over the synthetic dataset, miRPipe is observed to outperform the existing state-of-the-art pipelines (accuracy: 95.23% and F1-score: 94.17%). Analysis on all the three cancer datasets shows that miRPipe is able to extract more number of known dysregulated miRNAs or piRNAs from the datasets as compared to the existing pipelines.

Identification and profiling of stable microRNAs in hemolymph of young and old Locusta migratoria fifth instars

Since the discovery of the first microRNA (miRNA) in the nematode Caenorhabditis elegans, numerous novel miRNAs have been identified which can regulate presumably every biological process in a wide range of metazoan species. In accordance, several insect miRNAs have been identified and functionally characterized. While regulatory RNA pathways are traditionally described at an intracellular level, studies reporting on the presence and potential role of extracellular (small) sRNAs have been emerging in the last decade, mainly in mammalian systems. Interestingly, evidence in several species indicates the functional transfer of extracellular RNAs between donor and recipient cells, illustrating RNA-based intercellular communication. In insects, however, reports on extracellular small RNAs are emerging but the number of detailed studies is still very limited. Here, we demonstrate the presence of stable sRNAs in the hemolymph of the migratory locust, Locusta migratoria. Moreover, the levels of several extracellular miRNAs (ex-miRNAs) present in locust hemolymph differed significantly between young and old fifth nymphal instars. In addition, we performed a 'proof of principle' experiment which suggested that extracellularly delivered miRNA molecules are capable of affecting the locusts' development.

Hatching is modulated by microRNA-378a-3p derived from extracellular vesicles secreted by blastocysts

SignificanceHatching from the zona pellucida is a prerequisite for embryo implantation and is less likely to occur in vitro for reasons unknown. Extracellular vesicles (EVs) are secreted by the embryo into the culture medium. Yet the role that embryonic EVs and their cargo microRNAs (miRNAs) play in blastocyst hatching has not been elucidated, partially due to the difficulties of isolating them from low amounts of culture medium. Here, we optimized EV-miRNA isolation from medium conditioned by individually cultured bovine embryos and subsequently showed that miR-378a-3p, which was up-regulated in EVs secreted by blastocysts, plays a crucial role in promoting blastocyst hatching. This demonstrates the regulatory effect of miR-378-3p on hatching, which is an established embryo quality parameter linked with implantation.

The DGCR8 E518K mutation found in Wilms tumors leads to a partial miRNA processing defect that alters gene expression patterns and biological processes

Wilms tumor (WT) is the most common renal tumor in childhood. We and others have previously identified oncogenic driver mutations affecting the microprocessor genes DROSHA and DGCR8 that lead to altered miRNA expression patterns. In the case of DGCR8, a single recurrent hotspot mutation (E518K) was found in the RNA binding domain. To functionally assess this mutation in vitro, we generated mouse Dgcr8-KO embryonic stem cell (mESC) lines with an inducible expression of wild-type or mutant DGCR8, mirroring the hemizygous mutant expression seen in WT. RNA-seq analysis revealed significant differences of miRNA expression profiles in DGCR8-E518K compared to DGCR8-wild-type mESCs. The E518K mutation only led to a partial rescue of the reported miRNA processing defect in Dgcr8-KO, with selectively reduced expression of numerous canonical miRNAs. Nevertheless, DGCR8-E518K retained significant activity given its ability to still process many miRNAs. Subsequent to altered miRNA levels, the expression of mRNA targets was likewise changed. Functional assays showed that DGCR8-E518K cells still have a partial proliferation and differentiation defect but were able to rescue critical biological processes in embryoid body development. The stem cell program could be shut down and all three germ layers were formed. These findings suggest that the E518K mutation leads to a partial reduction of microprocessor activity and altered specificity with selective impairment only in certain developmental contexts, apparently including nephrogenesis.

Distinct and Coordinated Regulation of Small Non-coding RNAs by E2f1 and p53 During Drosophila Development and in Response to DNA Damage

Small non-coding RNAs (ncRNAs), including microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs), play a pivotal role in biological processes. A comprehensive quantitative reference of small ncRNAs expression during development and in DNA damage response (DDR) would significantly advance our understanding of their roles. In this study, we systemically analyzed the expression profile of miRNAs and piRNAs in wild-type flies, e2f1 mutant, p53 mutant and e2f1 p53 double mutant during development and after X-ray irradiation. By using small RNA sequencing and bioinformatic analysis, we found that both miRNAs and piRNAs were expressed in a dynamic mode and formed 4 distinct clusters during development. Notably, the expression pattern of miRNAs and piRNAs was changed in e2f1 mutant at multiple developmental stages, while retained in p53 mutant, indicating a critical role of E2f1 played in mediating small ncRNAs expression. Moreover, we identified differentially expressed (DE) small ncRNAs in e2f1 mutant and p53 mutant after X-ray irradiation. Furthermore, we mapped the binding motif of E2f1 and p53 around the small ncRNAs. Our data suggested that E2f1 and p53 work differently yet coordinately to regulate small ncRNAs expression, and E2f1 may play a major role to regulate miRNAs during development and after X-ray irradiation. Collectively, our results provide comprehensive characterization of small ncRNAs, as well as the regulatory roles of E2f1 and p53 in small ncRNAs expression, during development and in DNA damage response, which reveal new insights into the small ncRNAs biology.

miRMaster 2.0: multi-species non-coding RNA sequencing analyses at scale

Analyzing all features of small non-coding RNA sequencing data can be demanding and challenging. To facilitate this process, we developed miRMaster. After the analysis of over 125 000 human samples and 1.5 trillion human small RNA reads over 4 years, we present miRMaster 2 with a wide range of updates and new features. We extended our reference data sets so that miRMaster 2 now supports the analysis of eight species (e.g. human, mouse, chicken, dog, cow) and 10 non-coding RNA classes (e.g. microRNAs, piRNAs, tRNAs, rRNAs, circRNAs). We also incorporated new downstream analysis modules such as batch effect analysis or sample embeddings using UMAP, and updated annotation data bases included by default (miRBase, Ensembl, GtRNAdb). To accommodate the increasing popularity of single cell small-RNA sequencing data, we incorporated a module for unique molecular identifier (UMI) processing. Further, the output tables and graphics have been improved based on user feedback and new output formats that emerged in the community are now supported (e.g. miRGFF3). Finally, we integrated differential expression analysis with the miRNA enrichment analysis tool miEAA. miRMaster is freely available at https://www.ccb.uni-saarland.de/mirmaster2.

Identification of Novel miRNAs and Their Target Genes in the Response to Abscisic Acid in Arabidopsis

miRNAs are involved in various biological processes, including adaptive responses to abiotic stress. To understand the role of miRNAs in the response to ABA, ABA-responsive miRNAs were identified by small RNA sequencing in wild-type Arabidopsis, as well as in abi1td, mkkk17, and mkkk18 mutants. We identified 10 novel miRNAs in WT after ABA treatment, while in abi1td, mkkk17, and mkkk18 mutants, three, seven, and nine known miRNAs, respectively, were differentially expressed after ABA treatment. One novel miRNA (miRn-8) was differentially expressed in the mkkk17 mutant. Potential target genes of the miRNA panel were identified using psRNATarget. Sequencing results were validated by quantitative RT-PCR of several known and novel miRNAs in all genotypes. Of the predicted targets of novel miRNAs, seven target genes of six novel miRNAs were further validated by 5' RLM-RACE. Gene ontology analyses showed the potential target genes of ABA-responsive known and novel miRNAs to be involved in diverse cellular processes in plants, including development and stomatal movement. These outcomes suggest that a number of the identified miRNAs have crucial roles in plant responses to environmental stress, as well as in plant development, and might have common regulatory roles in the core ABA signaling pathway.

Different Expressions of Pericardial Fluid MicroRNAs in Patients With Arrhythmogenic Right Ventricular Cardiomyopathy and Ischemic Heart Disease Undergoing Ventricular Tachycardia Ablation

Introduction: Pericardial fluid is enriched with biologically active molecules of cardiovascular origin including microRNAs. Investigation of the disease-specific extracellular microRNAs could shed light on the molecular processes underlying disease development. Arrhythmogenic right ventricular cardiomyopathy (ARVC) is an inherited heart disease characterized by life-threatening arrhythmias and progressive heart failure development. The current data about the association between microRNAs and ARVC development are limited.

Methods and Results: We performed small RNA sequence analysis of microRNAs of pericardial fluid samples obtained during transcutaneous epicardial access for ventricular tachycardia (VT) ablation of six patients with definite ARVC and three post-infarction VT patients. Disease-associated microRNAs of pericardial fluid were identified. Five microRNAs (hsa-miR-1-3p, hsa-miR-21-5p, hsa-miR-122-5p, hsa-miR-206, and hsa-miR-3679-5p) were found to be differentially expressed between patients with ARVC and patients with post-infarction VT. Enrichment analysis of differentially expressed microRNAs revealed their close linkage to cardiac diseases.

Conclusion: Our data extend the knowledge of pericardial fluid microRNA composition and highlight five pericardial fluid microRNAs potentially linked to ARVC pathogenesis. Further studies are required to confirm the use of pericardial fluid RNA sequencing in differential diagnosis of ARVC.

Pi-starvation induced transcriptional changes in barley revealed by a comprehensive RNA-Seq and degradome analyses

BACKGROUND

Small RNAs (sRNAs) are 20-30 nt regulatory elements which are responsible for plant development regulation and participate in many plant stress responses. Insufficient inorganic phosphate (Pi) concentration triggers plant responses to balance the internal Pi level.

RESULTS

In this study, we describe Pi-starvation-responsive small RNAs and transcriptome changes in barley (Hordeum vulgare L.) using Next-Generation Sequencing (NGS) RNA-Seq data derived from three different types of NGS libraries: (i) small RNAs, (ii) degraded RNAs, and (iii) functional mRNAs. We find that differentially and significantly expressed miRNAs (DEMs, Bonferroni adjusted p-value < 0.05) are represented by 15 molecules in shoot and 13 in root; mainly various miR399 and miR827 isomiRs. The remaining small RNAs (i.e., those without perfect match to reference sequences deposited in miRBase) are considered as differentially expressed other sRNAs (DESs, p-value Bonferroni correction < 0.05). In roots, a more abundant and diverse set of other sRNAs (DESs, 1796 unique sequences, 0.13% from the average of the unique small RNA expressed under low-Pi) contributes more to the compensation of low-Pi stress than that in shoots (DESs, 199 unique sequences, 0.01%). More than 80% of differentially expressed other sRNAs are up-regulated in both organs. Additionally, in barley shoots, up-regulation of small RNAs is accompanied by strong induction of two nucleases (S1/P1 endonuclease and 3'-5' exonuclease). This suggests that most small RNAs may be generated upon nucleolytic cleavage to increase the internal Pi pool. Transcriptomic profiling of Pi-starved barley shoots identifies 98 differentially expressed genes (DEGs). A majority of the DEGs possess characteristic Pi-responsive cis-regulatory elements (P1BS and/or PHO element), located mostly in the proximal promoter regions. GO analysis shows that the discovered DEGs primarily alter plant defense, plant stress response, nutrient mobilization, or pathways involved in the gathering and recycling of phosphorus from organic pools.

CONCLUSIONS

Our results provide comprehensive data to demonstrate complex responses at the RNA level in barley to maintain Pi homeostasis and indicate that barley adapts to Pi-starvation through elicitation of RNA degradation. Novel P-responsive genes were selected as putative candidates to overcome low-Pi stress in barley plants.

Identification of microRNAs and their gene targets in cytoplasmic male sterile and fertile maintainer lines of pigeonpea

Comparative analysis of genome-wide miRNAs and their gene targets between cytoplasmic male sterile (CMS) and fertile lines of pigeonpea suggests a possible role of miRNA-regulated pathways in reproductive development. Exploitation of hybrid vigor using CMS technology has delivered nearly 50% yield gain in pigeonpea. Among various sterility-inducing cytoplasms (A₁-A₉) reported so far in pigeonpea, A₂ and A₄ are the two major sources that facilitate hybrid seed production. Recent evidence suggests involvement of micro RNA in vast array of biological processes including plant reproductive development. In pigeonpea, information about the miRNAs is insufficient. In view of this, we sequenced six small RNA libraries of CMS line UPAS 120A and isogenic fertile line UPAS 120B using Illumina technology. Results revealed 316 miRNAs including 248 known and 68 novel types. A total of 637 gene targets were predicted for known miRNAs, while 324 genes were associated with novel miRNAs. Degradome analysis revealed 77 gene targets of predicted miRNAs, which included a variety of transcription factors playing key roles in plant reproduction such as F-box family proteins, apetala 2, auxin response factors, ethylene-responsive factors, homeodomain-leucine zipper proteins etc. Differential expression of both known and novel miRNAs implied roles for both conserved as well as species-specific players. We also obtained several miRNA families such as miR156, miR159, miR167 that are known to influence crucial aspects of plant fertility. Gene ontology and pathway level analyses of the target genes showed their possible implications for crucial events during male reproductive development such as tapetal degeneration, pollen wall formation, retrograde signaling etc. To the best of our knowledge, present study is first to combine deep sequencing of small RNA and degradome for elucidating the role of miRNAs in flower and male reproductive development in pigeonpea.

MicroRNAs Regulating Autophagy in Neurodegeneration

Social and economic impacts of neurodegenerative diseases (NDs) become more prominent in our constantly aging population. Currently, due to the lack of knowledge about the aetiology of most NDs, only symptomatic treatment is available for patients. Hence, researchers and clinicians are in need of solid studies on pathological mechanisms of NDs. Autophagy promotes degradation of pathogenic proteins in NDs, while microRNAs post-transcriptionally regulate multiple signalling networks including autophagy. This chapter will critically discuss current research advancements in the area of microRNAs regulating autophagy in NDs. Moreover, we will introduce basic strategies and techniques used in microRNA research. Delineation of the mechanisms contributing to NDs will result in development of better approaches for their early diagnosis and effective treatment.

isomiRs-Hidden Soldiers in the miRNA Regulatory Army, and How to Find Them?

Numerous studies on microRNAs (miRNA) in cancer and other diseases have been accompanied by diverse computational approaches and experimental methods to predict and validate miRNA biological and clinical significance as easily accessible disease biomarkers. In recent years, the application of the next-generation deep sequencing for the analysis and discovery of novel RNA biomarkers has clearly shown an expanding repertoire of diverse sequence variants of mature miRNAs, or isomiRs, resulting from alternative post-transcriptional processing events, and affected by (patho)physiological changes, population origin, individual's gender, and age. Here, we provide an in-depth overview of currently available bioinformatics approaches for the detection and visualization of both mature miRNA and cognate isomiR sequences. An attempt has been made to present in a systematic way the advantages and downsides of in silico approaches in terms of their sensitivity and accuracy performance, as well as used methods, workflows, and processing steps, and end output dataset overlapping issues. The focus is given to the challenges and pitfalls of isomiR expression analysis. Specifically, we address the availability of tools enabling research without extensive bioinformatics background to explore this fascinating corner of the small RNAome universe that may facilitate the discovery of new and more reliable disease biomarkers.

Encyclopedia of tools for the analysis of miRNA isoforms

RNA sequencing data sets rapidly increase in quantity. For microRNAs (miRNAs), frequently dozens to hundreds of billion reads are generated per study. The quantification of annotated miRNAs and the prediction of new miRNAs are leading computational tasks. Now, the increased depth of coverage allows to gain deeper insights into the variability of miRNAs. The analysis of isoforms of miRNAs (isomiRs) is a trending topic, and a range of computational tools for the analysis of isomiRs has been developed. We provide an overview on 27 available computational solutions for the analysis of isomiRs. These include both stand-alone programs (17 tools) and web-based solutions (10 tools) and span a publication time range from 2010 to 2020. Seven of the tools were published in 2019 and 2020, confirming the rising importance of the topic. While most of the analyzed tools work for a broad range of organisms or are completely independent of a reference organism, several tools have been tailored for the analysis of human miRNA data or for plants. While 14 of the tools are general analysis tools of miRNAs, and isomiR analysis is one of their features, the remaining 13 tools have specifically been developed for isomiR analysis. A direct comparison on 20 deep sequencing data sets for selected tools provides insights into the heterogeneity of results. With our work, we provide users a comprehensive overview on the landscape of isomiR analysis tools and in that support the selection of the most appropriate tool for their respective research task.

Ross River Virus Provokes Differentially Expressed MicroRNA and RNA Interference Responses in Aedes aegypti Mosquitoes

Alphaviruses are globally distributed and predominately transmitted by mosquitoes. Aedes species are common vectors for the clinically important alphaviruses-Chikungunya, Sindbis, and Ross River (RRV) viruses-with Aedes aegypti also being a vector for the flaviviruses dengue, Yellow Fever, and Zika viruses. Ae. aegypti was putatively implicated in the large 1979-1980 South Pacific Islands outbreak of RRV-the leading cause of arboviral disease in Australia today. The RNA interference (RNAi) defense response in mosquitoes involves a number of small RNAs, with their kinetics induced by alphaviruses being poorly understood, particularly at the tissue level. We compared the small RNA profiles between RRV-infected and noninfected Ae. aegypti midgut and fat body tissues at 2, 6, and 12 days post-inoculation (dpi). RRV induced an incremental RNAi response, yielding short interfering and P-element-induced-wimpy-testis (PIWI)-interacting RNAs. Fourteen host microRNAs were differentially expressed due to RRV with the majority in the fat body at 2 dpi. The largely congruent pattern of microRNA regulation with previous reports for alphaviruses and divergence from those for flaviviruses suggests a degree of conservation, whereas patterns of microRNA expression unique to this study provide novel insights into the tissuespecific hostvirus attributes of Ae. aegypti responses to this previously unexplored oldworld alphavirus.

Maternal Recognition of Pregnancy in the Horse: Are MicroRNAs the Secret Messengers?

The signal for maternal recognition of pregnancy (MRP) has still not been identified in the horse. High-throughput molecular biology at the embryo-maternal interface has substantially contributed to the knowledge on pathways affected during MRP, but an integrated study in which proteomics, transcriptomics and miRNA expression can be linked directly is currently lacking. The aim of this study was to provide such analysis. Endometrial biopsies, uterine fluid, embryonic tissues, and yolk sac fluid were collected 13 days after ovulation during pregnant and control cycles from the same mares. Micro-RNA-Sequencing was performed on all collected samples, mRNA-Sequencing on the same tissue samples and mass spectrometry was conducted previously on the same fluid samples. Differential expression of miRNA, mRNA and proteins showed high conformity with literature and confirmed involvement in pregnancy establishment, embryo quality, steroid synthesis and prostaglandin regulation, but the link between differential miRNAs and their targets was limited and did not indicate the identity of an unequivocal signal for MRP in the horse. Differential expression at the embryo-maternal interface was prominent, highlighting a potential role of miRNAs in embryo-maternal communication during early pregnancy in the horse. These data provide a strong basis for future targeted studies.

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.

Deciphering miRNAs' Action through miRNA Editing

MicroRNAs (miRNAs) are small non-coding RNAs with the capability of modulating gene expression at the post-transcriptional level either by inhibiting messenger RNA (mRNA) translation or by promoting mRNA degradation. The outcome of a myriad of physiological processes and pathologies, including cancer, cardiovascular and metabolic diseases, relies highly on miRNAs. However, deciphering the precise roles of specific miRNAs in these pathophysiological contexts is challenging due to the high levels of complexity of their actions. Indeed, regulation of mRNA expression by miRNAs is frequently cell/organ specific; highly dependent on the stress and metabolic status of the organism; and often poorly correlated with miRNA expression levels. Such biological features of miRNAs suggest that various regulatory mechanisms control not only their expression, but also their activity and/or bioavailability. Several mechanisms have been described to modulate miRNA action, including genetic polymorphisms, methylation of miRNA promoters, asymmetric miRNA strand selection, interactions with RNA-binding proteins (RBPs) or other coding/non-coding RNAs. Moreover, nucleotide modifications (A-to-I or C-to-U) within the miRNA sequences at different stages of their maturation are also critical for their functionality. This regulatory mechanism called "RNA editing" involves specific enzymes of the adenosine/cytidine deaminase family, which trigger single nucleotide changes in primary miRNAs. These nucleotide modifications greatly influence a miRNA's stability, maturation and activity by changing its specificity towards target mRNAs. Understanding how editing events impact miRNA's ability to regulate stress responses in cells and organs, or the development of specific pathologies, e.g., metabolic diseases or cancer, should not only deepen our knowledge of molecular mechanisms underlying complex diseases, but can also facilitate the design of new therapeutic approaches based on miRNA targeting. Herein, we will discuss the current knowledge on miRNA editing and how this mechanism regulates miRNA biogenesis and activity.

Computational Resources for Prediction and Analysis of Functional miRNA and Their Targetome

microRNAs are evolutionarily conserved, endogenously produced, noncoding RNAs (ncRNAs) of approximately 19-24 nucleotides (nts) in length known to exhibit gene silencing of complementary target sequence. Their deregulated expression is reported in various disease conditions and thus has therapeutic implications. In the last decade, various computational resources are published in this field. In this chapter, we have reviewed bioinformatics resources, i.e., miRNA-centered databases, algorithms, and tools to predict miRNA targets. First section has enlisted more than 75 databases, which mainly covers information regarding miRNA registries, targets, disease associations, differential expression, interactions with other noncoding RNAs, and all-in-one resources. In the algorithms section, we have compiled about 140 algorithms from eight subcategories, viz. for the prediction of precursor (pre-) and mature miRNAs. These algorithms are developed on various sequence, structure, and thermodynamic based features incorporated into different machine learning techniques (MLTs). In addition, computational identification of miRNAs from high-throughput next generation sequencing (NGS) data and their variants, viz. isomiRs, differential expression, miR-SNPs, and functional annotation, are discussed. Prediction and analysis of miRNAs and their associated targets are also evaluated under miR-targets section providing knowledge regarding novel miRNA targets and complex host-pathogen interactions. In conclusion, we have provided comprehensive review of in silico resources published in miRNA research to help scientific community be updated and choose the appropriate tool according to their needs.

Altered microRNA expression profiles in large offspring syndrome and Beckwith-Wiedemann syndrome

The use of assisted reproductive technologies (ART) can induce a congenital overgrowth condition in humans and ruminants, namely Beckwith-Wiedemann syndrome (BWS) and large offspring syndrome (LOS), respectively. Shared phenotypes and epigenotypes have been found between BWS and LOS. We have observed global misregulation of transcripts in bovine foetuses with LOS. microRNAs (miRNAs) are important post-transcriptional gene expression regulators. We hypothesize that there is miRNA misregulation in LOS and that this misregulation is shared with BWS. In this study, small RNA sequencing was conducted to investigate miRNA expression profiles in bovine and human samples. We detected 407 abundant known miRNAs and predicted 196 putative miRNAs from the bovine sequencing results and identified 505 abundant miRNAs in human tongue. Differentially expressed miRNAs (DE-miRNAs) were identified between control and LOS groups in all tissues analysed as well as between BWS and control human samples. DE-miRNAs were detected from several miRNA clusters including DLK1-DIO3 genomic imprinted cluster in LOS and BWS. DNA hypermethylation was associated with downregulation of miRNAs in the DLK1-DIO3. mRNA targets of the DE-miRNAs were predicted and signalling pathways associated with control of organ size (including the Hippo signalling pathway), cell proliferation, apoptosis, cell survival, cell cycle, and cell adhesion were found to be enriched with these genes. Yes associated protein 1 (YAP1) is the core effector of the Hippo signalling pathway, and increased level of active (non-phosphorylated) YAP1 protein was detected in skeletal muscle of LOS foetuses. Overall, our data provide evidence of miRNA misregulation in LOS and BWS.

Bovine Embryo-Secreted microRNA-30c Is a Potential Non-invasive Biomarker for Hampered Preimplantation Developmental Competence

Recently, secreted microRNAs (miRNAs) have received a lot of attention since they may act as autocrine factors. However, how secreted miRNAs influence embryonic development is still poorly understood. We identified 294 miRNAs, 114 known, and 180 novel, in the conditioned medium of individually cultured bovine embryos. Of these miRNAs, miR-30c and miR-10b were much more abundant in conditioned medium of slow cleaving embryos compared to intermediate cleaving ones. MiR-10b, miR-novel-44, and miR-novel-45 were higher expressed in the conditioned medium of degenerate embryos compared to blastocysts, while the reverse was observed for miR-novel-113 and miR-novel-139. Supplementation of miR-30c mimics into the culture medium confirmed the uptake of miR-30c mimics by embryos and resulted in increased cell apoptosis, as also shown after delivery of miR-30c mimics in Madin-Darby bovine kidney cells (MDBKs). We also demonstrated that miR-30c directly targets Cyclin-dependent kinase 12 (CDK12) through its 3′ untranslated region (3′-UTR) and inhibits its expression. Overexpression and downregulation of CDK12 revealed the opposite results of the delivery of miRNA-30c mimics and inhibitor. The significant down-regulation of several tested DNA damage response (DDR) genes, after increasing miR-30c or reducing CDK12 expression, suggests a possible role for miR-30c in regulating embryo development through DDR pathways.

MicroRNA Sequencing Identifies a Serum MicroRNA Panel, Which Combined With Aspartate Aminotransferase to Platelet Ratio Index Can Detect and Monitor Liver Disease in Pediatric Cystic Fibrosis

Cystic fibrosis (CF)-associated liver disease (CFLD) is a hepatobiliary complication of CF. Current diagnostic modalities rely on nonspecific assessments, whereas liver biopsy is the gold standard to assess severity of fibrosis. MicroRNAs (miRNAs) regulate liver disease pathogenesis and are proposed as diagnostic biomarkers. We investigated the combined use of serum miRNAs and aspartate aminotransferase (AST) to platelet ratio (APRI) to diagnose and assess CFLD severity. This was a cross-sectional cohort study of the circulatory miRNA signature of 124 children grouped by clinical, biochemical, and imaging assessments as follows: CFLD (n = 44), CF patients with no evidence of liver disease (CFnoLD; n = 40), and healthy controls (n = 40). Serum miRNAs were analyzed using miRNA sequencing (miRNA-Seq). Selected differentially expressed serum miRNA candidates were further validated by qRT-PCR and statistical analysis performed to evaluate utility to predict CFLD and fibrosis severity validated by liver biopsy, alone or in combination with APRI. Serum miR-122-5p, miR-365a-3p, and miR-34a-5p levels were elevated in CFLD compared to CFnoLD, whereas miR-142-3p and let-7g-5p were down-regulated in CFLD compared to CFnoLD. Logistic regression analysis combining miR-365a-3p, miR-142-3p, and let-7g-5p with APRI showed 21 times greater odds of accurately predicting liver disease in CF with an area under the receiver operating characteristics curve (AUROC) = 0.91 (sensitivity = 83%, specificity = 92%; P < 0.0001). Expression levels of serum miR-18a-5p were correlated with increasing hepatic fibrosis (HF) stage in CFLD (r_s  = 0.56; P < 0.0001), showing good diagnostic accuracy for distinguishing severe (F3-F4) from mild/moderate fibrosis (F0-F2). A unit increase of miR-18a-5p showed a 7-fold increased odds of having severe fibrosis with an AUROC = 0.82 (sensitivity = 93%, specificity = 73%; P = 0.004), indicating its potential to predict fibrosis severity. Conclusion: We identified a distinct circulatory miRNA profile in pediatric CFLD with potential to accurately discriminate liver disease and fibrosis severity in children with CF.

Genome-wide transcriptome analysis reveals small RNA profiles involved in early stages of stolon-to-tuber transitions in potato under photoperiodic conditions

BACKGROUND

Small RNAs (sRNAs), especially miRNAs, act as crucial regulators of plant growth and development. Two other sRNA groups, trans-acting short-interfering RNAs (tasiRNAs) or phased siRNAs (phasiRNAs), are also emerging as potential regulators of plant development. Stolon-to-tuber transition in potato is an important developmental phase governed by many environmental, biochemical and hormonal cues. Among different environmental factors, photoperiod has a major influence on tuberization. Several mobile signals, mRNAs, proteins and transcription factors have been widely studied for their role in tuber formation in potato, however, no information is yet available that describes the molecular signals governing the early stages of stolon transitions or cell-fate changes at the stolon tip before it matures to potato. Stolon could be an interesting model for studying below ground organ development and we hypothesize that small RNAs might be involved in regulation of stolon-to-tuber transition process in potato. Also, there is no literature that describes the phased siRNAs in potato development.

RESULTS

We performed sRNA profiling of early stolon stages (4, 7 and 10 d) under long-day (LD; 16 h light, 8 h dark) and short-day (SD; 8 h light, 16 h dark) photoperiodic conditions. Altogether, 7 (out of 324) conserved and 12 (out of 311) novel miRNAs showed differential expression in early stolon stages under SD vs LD photoperiodic conditions. Key target genes (StGRAS, StTCP2/4 and StPTB6) exhibited differential expression in early stolon stages under SD vs LD photoperiodic conditions, indicative of their potential role in tuberization. Out of 830 TAS-like loci identified, 24 were cleaved by miRNAs to generate 190 phased siRNAs. Some of them targeted crucial tuberization genes such as StPTB1, POTH1 and StCDPKs. Two conserved TAS loci, referred as StTAS3 and StTAS5, which share close conservation with members of the Solanaceae family, were identified in our analysis. One TAS-like locus (StTm2) was validated for phased siRNA generation and one of its siRNA was predicted to cleave an important tuber marker gene StGA2ox1.

CONCLUSION

Our study suggests that sRNAs and their selective target genes could be associated with the regulation of early stages of stolon-to-tuber transitions in a photoperiod-dependent manner in potato.

Extracellular vesicle-associated miRNAs in ovarian cancer – design of an integrated NGS-based workflow for the identification of blood-based biomarkers for platinum-resistance

Background

Extracellular vesicle (EV)-associated microRNAs (miRNAs) have been suggested as promising biomarkers for blood-based cancer diagnosis. However, one of the major limitations for the use of EVs with diagnostic purpose is the lack of standardized EV-profiling techniques. In this regard, the objective of our study was to design an integrated next-generation sequencing (NGS)-based workflow for analyzing the signature of EV-associated miRNA in the plasma of platinum-resistant ovarian cancer patients.

Methods

For EV-extraction, different enrichment methods were compared (ExoQuick vs. exoRNeasy). NGS was performed with the Illumina platform.

Results

We established an integrated NGS-based workflow, including EV-enrichment with the ExoQuick system, which resulted in an optimal RNA-yield and consistent small RNA libraries. We applied this workflow in a pilot cohort of clinically documented platinum-sensitive (n=15) vs. platinum-resistant (n=15) ovarian cancer patients, resulting in a panel of mature EV-associated miRNAs (including ovarian cancer associated miR-181a, miR-1908, miR-21, miR-486 and miR-223), which were differentially abundant in the plasma of platinum-resistant patients.

Conclusions

This is the first study, analyzing the profile of EV-associated miRNAs in platinum-resistant ovarian cancer patients. We provide rationale to further validate these miRNA candidates in an independent set of patients, in order to characterize their biomarker potential as predictors for platinum-resistance.

Small RNA-Seq reveals novel miRNAs shaping the transcriptomic identity of rat brain structures

Small RNA-Seq of the rat central nervous system reveals known and novel miRNAs specifically regulated in brain structures and correlated with the expression of their predicted target genes, suggesting a critical role in the transcriptomic identity of brain structures.

In the central nervous system (CNS), miRNAs are involved in key functions, such as neurogenesis and synaptic plasticity. Moreover, they are essential to define specific transcriptomes in tissues and cells. However, few studies were performed to determine the miRNome of the different structures of the rat CNS, although a major model in neuroscience. Here, we determined by small RNA-Seq, the miRNome of the olfactory bulb, the hippocampus, the cortex, the striatum, and the spinal cord and showed the expression of 365 known miRNAs and 90 novel miRNAs. Differential expression analysis showed that several miRNAs were specifically enriched/depleted in these CNS structures. Transcriptome analysis by mRNA-Seq and correlation based on miRNA target predictions suggest that the specifically enriched/depleted miRNAs have a strong impact on the transcriptomic identity of the CNS structures. Altogether, these results suggest the critical role played by these enriched/depleted miRNAs, in particular the novel miRNAs, in the functional identities of CNS structures.

miRmapper: A Tool for Interpretation of miRNA⁻mRNA Interaction Networks

It is estimated that 30% of all genes in the mammalian cells are regulated by microRNA (miRNAs). The most relevant miRNAs in a cellular context are not necessarily those with the greatest change in expression levels between healthy and diseased tissue. Differentially expressed (DE) miRNAs that modulate a large number of messenger RNA (mRNA) transcripts ultimately have a greater influence in determining phenotypic outcomes and are more important in a global biological context than miRNAs that modulate just a few mRNA transcripts. Here, we describe the development of a tool, “miRmapper”, which identifies the most dominant miRNAs in a miRNA–mRNA network and recognizes similarities between miRNAs based on commonly regulated mRNAs. Using a list of miRNA–target gene interactions and a list of DE transcripts, miRmapper provides several outputs: (1) an adjacency matrix that is used to calculate miRNA similarity utilizing the Jaccard distance; (2) a dendrogram and (3) an identity heatmap displaying miRNA clusters based on their effect on mRNA expression; (4) a miRNA impact table and (5) a barplot that provides a visual illustration of this impact. We tested this tool using nonmetastatic and metastatic bladder cancer cell lines and demonstrated that the most relevant miRNAs in a cellular context are not necessarily those with the greatest fold change. Additionally, by exploiting the Jaccard distance, we unraveled novel cooperative interactions between miRNAs from independent families in regulating common target mRNAs; i.e., five of the top 10 miRNAs act in synergy.

Identifying and characterizing functional 3' nucleotide addition in the miRNA pathway

Over the past decade, modifications to microRNAs (miRNAs) via 3' end nucleotide addition have gone from a deep-sequencing curiosity to experimentally confirmed drivers of a range of regulatory activities. Here we overview the methods that have been deployed by researchers seeking to untangle these diverse functional roles and include characterizing not only the nucleotidyl transferases catalyzing the additions but also the nucleotides being added, and the timing of their addition during the miRNA pathway. These methods and their further development are key to clarifying the diverse and sometimes contradictory functional findings presently attributed to these nucleotide additions.

Large-scale screening of circulating microRNAs in individuals with HIV-1 mono-infections reveals specific liver damage signatures

Human immunodeficiency virus type 1 (HIV-1)-induced inflammation and/or long-term antiretroviral drug toxicity may contribute to the evolution of liver disease. We investigated circulating plasma microRNAs (miRNAs) as potential biomarkers of liver injury in patients mono-infected with HIV-1. We performed large-scale deep sequencing analyses of small RNA level on plasma samples from patients with HIV-1 mono-infection that had elevated or normal levels of alanine aminotransferase (ALT) or focal nodular hyperplasia (FNH). Hepatitis C virus (HCV) mono-infected patients were also studied. Compared to healthy donors, patients with HIV-1 or HCV mono-infections showed significantly altered (fold change >2, adjusted p < 0.05) level of 25 and 70 miRNAs, respectively. Of the 25 altered miRNAs found in patients with HIV-1, 19 were also found in patients mono-infected with HCV. Moreover, 13 of the 14 most up-regulated miRNAs (range: 9.3-3.4-fold increase) in patients with HCV mono-infections were also up-regulated in patients with HIV-1 mono-infections. Importantly, most of these miRNAs significantly and positively correlated with ALT and aspartate aminotransferase (AST) levels, and liver fibrosis stage (p < 0.05). MiR-122-3p and miR-193b-5p were highly up-regulated HIV-1 mono-infected patients with elevated ALT or FNH, but not in HIV-1 patients with normal levels of ALT. These results reveal that HIV-1 infections impacted liver-related miRNA levels in the absence of an HCV co-infection, which highlights the potential of miRNAs as biomarkers for the progression of liver injury in HIV-1 infected patients.

Bioinformatic Analysis of MicroRNA Sequencing Data

The vital role of microRNAs (miRNAs) involved in gene expression regulation has been confirmed in many biological processes. With the growing power and reducing cost of next-generation sequencing, more and more researchers turn to apply this high-throughput method to solve their biological problems. For miRNAs with known sequences, their expression profiles can be generated from the sequencing data. It also allows us to identify some novel miRNAs and explore the sequence variations under different conditions. Currently, there are a handful of tools available to analyze the miRNA sequencing data with separated or combined features, such as reads preprocessing, mapping and differential expression analysis. However, to our knowledge, a hands-on guideline for miRNA sequencing data analysis covering all steps is not available. Here we will utilize a set of published tools to perform the miRNA analysis with detailed explanation. Particularly, the miRNA target prediction and annotation may provide useful information for further experimental verification.

Association of microRNAs with Argonaute proteins in the malaria mosquito Anopheles gambiae after blood ingestion

Drastic changes in gene expression occur after adult female mosquitoes take a blood meal and use the nutrients for egg maturation. A growing body of evidence indicates that microRNAs (miRNAs) contribute to this tightly controlled tissue- and stage-specific gene expression. To investigate the role of miRNAs, we monitored miRNA expression in the mosquito Anopheles gambiae during the 72-h period immediately after blood feeding. We also measured the association of miRNAs with Argonaute 1 (Ago1) and Argonaute 2 (Ago2) to assess the functional status of individual miRNA species. Overall, 173 mature miRNAs were precipitated with Ago1 and Ago2, including 12 new miRNAs, the orthologs of which are found thus far only in other Anopheles species. Ago1 is the predominant carrier of miRNAs in Anopheles gambiae. The abundance and Ago loading of most of the mature miRNAs were relatively stable after blood ingestion. However, miRNAs of the miR-309/286/2944 cluster were considerably upregulated after blood feeding. Injection of the specific antagomir for miR-309 resulted in smaller developing oocytes and ultimately fewer eggs. In addition, the Ago association of some miRNAs was not proportional to their cellular abundance, suggesting that integration of miRNAs into the Ago complexes is regulated by additional mechanisms.