miR-isomiRExp: a web-server for the analysis of expression of miRNA at the miRNA/isomiR levels

LncRNA-miRNA-mRNA Network Analysis Reveals the Potential Biomarkers in Crohn’s Disease Rats Treated with Herb-Partitioned Moxibustion

Background

Long noncoding RNA (lncRNA) is receiving growing attention in Crohn’s disease (CD). However, the mechanism by which herb-partitioned moxibustion (HPM) regulates the expression and functions of lncRNAs in CD rats is still unclear. The aim of our study is to identify lncRNA-miRNA-mRNA network potential biological functions in CD.

Methods

RNA sequencing and microRNA (miRNA) sequencing were carried out to analyze lncRNA, miRNA and mRNA expression profiles among the CD rats, normal control rats, and CD rats after HPM treatment and constructed the potential related lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) networks. Then, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, protein–protein interaction (PPI) analysis and quantitative real-time polymerase chain reaction (qRT-PCR) were performed to explore potentially important genes in ceRNA networks.

Results

A total of 189 lncRNAs, 32 miRNAs and 463 mRNAs were determined as differentially expressed (DE) genes in CD rats compared to normal control rats, and 161 lncRNAs, 12 miRNAs and 130 mRNAs were identified as remarkably DE genes in CD rats after HPM treatment compared to CD rats. GO analysis indicated that the target genes were most enriched in cAMP and in KEGG pathway analysis the main pathways included adipocytokine, PPAR, AMPK, FoxO and PI3K-Akt signaling pathway. Finally, qRT-PCR results confirmed that lncRNA LOC102550026 sponged miRNA-34c-5p to regulate the intestinal immune inflammatory response by targeting Pck1.

Conclusion

By constructing a ceRNA network with lncRNA-miRNA-mRNA, PCR verification, and KEGG analysis, we revealed that LOC102550026/miRNA-34c-5p/Pck1 axis and adipocytokine, PPAR, AMPK, FoxO, and PI3K-Akt signaling pathways might regulate the intestinal immune-inflammatory response, and HPM may regulate the lncRNA LOC102550026/miR-34c-5p/Pck1 axis and adipocytokine, PPAR, AMPK, FoxO, and PI3K-Akt signaling pathways, thus improving intestinal inflammation in CD. These findings may be novel potential targets in CD.

MicroRNA Databases and Tools

In this era of big data, sets of methodologies and strategies are designed to extract knowledge from huge volumes of data. However, the cost of where and how to get this information accurately and quickly is extremely important, given the diversity of genomes and the different ways of representing that information. Among the huge set of information and relationships that the genome carries, there are sequences called miRNAs (microRNAs). These sequences were described in the 1990s and are mainly involved in mechanisms of regulation and gene expression. Having this in mind, this chapter focuses on exploring the available literature and providing useful and practical guidance on the miRNA database and tools topic. For that, we organized and present this text in two ways: (a) the update reviews and articles, which best summarize and discuss the theme; and (b) our update investigation on miRNA literature and portals about databases and tools. Finally, we present the main challenge and a possible solution to improve resources and tools.

Systems analysis of miRNA biomarkers to inform drug safety

microRNAs (miRNAs or miRs) are short non-coding RNA molecules which have been shown to be dysregulated and released into the extracellular milieu as a result of many drug and non-drug-induced pathologies in different organ systems. Consequently, circulating miRs have been proposed as useful biomarkers of many disease states, including drug-induced tissue injury. miRs have shown potential to support or even replace the existing traditional biomarkers of drug-induced toxicity in terms of sensitivity and specificity, and there is some evidence for their improved diagnostic and prognostic value. However, several pre-analytical and analytical challenges, mainly associated with assay standardization, require solutions before circulating miRs can be successfully translated into the clinic. This review will consider the value and potential for the use of circulating miRs in drug-safety assessment and describe a systems approach to the analysis of the miRNAome in the discovery setting, as well as highlighting standardization issues that at this stage prevent their clinical use as biomarkers. Highlighting these challenges will hopefully drive future research into finding appropriate solutions, and eventually circulating miRs may be translated to the clinic where their undoubted biomarker potential can be used to benefit patients in rapid, easy to use, point-of-care test systems.

Comparative Mucous miRomics in Cynoglossus semilaevis Related to Vibrio harveyi Caused Infection

Epidermal mucus is an important barrier and regulating mediator in fish. MicroRNAs (miRNAs) are proved to be involved in various biological processes, also as promising biomarkers for disease diagnosis. Vibrio harveyi has long been a noticeable bacterial pathogen in Cynoglossus semilaevis aquaculture. To find the evidence whether there are indicating miRNAs in mucus and whether the miRNAs are related to infections caused by V. harveyi, miRNA profiles of mucus from V. harveyi infected fish and healthy controls were screened by small RNA sequencing and verified by quantitative real-time PCR. This is the first report about miRNA profiling of flatfish mucus, aiming at illustrating the pathogenesis of V. harveyi caused infection and developing disease-related biomarkers. The results revealed significant differences in expression levels of some miRNAs between infected fish and healthy ones. Three hundred differentially expressed miRNAs were obtained after filtering through FC > 2 or FC < 0.5 and most of the differential miRNAs were downregulated. After verification through qRT-PCR, four unique miRNAs, dre-miR-451, dre-miR-184, dre-miR-205-5p > ssa-miR-205b-5p, and dre-miR-181a-5p > ssa-miR-181a-5p, were identified as V. harveyi infection-related signatures, consistent with sequencing trend. The expression levels of these four miRNAs in the infected fish were all significantly lower than controls. These miRNAs in mucus could be used to differentiate diseased and healthy fish in a non-invasive way with practical value for large-scale disease screening. They also provided new insights into the mechanism underlying the bacterial infections in fish.

Role of SNPs in the Biogenesis of Mature miRNAs

Single nucleotide polymorphisms (SNPs) play a significant role in microRNA (miRNA) generation, processing, and function and contribute to multiple phenotypes and diseases. Therefore, whole-genome analysis of how SNPs affect miRNA maturation mechanisms is important for precision medicine. The present study established an SNP-associated pre-miRNA (SNP-pre-miRNA) database, named miRSNPBase, and constructed SNP-pre-miRNA sequences. We also identified phenotypes and disease biomarker-associated isoform miRNA (isomiR) based on miRFind, which was developed in our previous study. We identified functional SNPs and isomiRs. We analyzed the biological characteristics of functional SNPs and isomiRs and studied their distribution in different ethnic groups using whole-genome analysis. Notably, we used individuals from Great Britain (GBR) as examples and identified isomiRs and isomiR-associated SNPs (iso-SNPs). We performed sequence alignments of isomiRs and miRNA sequencing data to verify the identified isomiRs and further revealed GBR ethnographic epigenetic dominant biomarkers. The SNP-pre-miRNA database consisted of 886 pre-miRNAs and 2640 SNPs. We analyzed the effects of SNP type, SNP location, and SNP-mediated free energy change during mature miRNA biogenesis and found that these factors were closely associated to mature miRNA biogenesis. Remarkably, 158 isomiRs were verified in the miRNA sequencing data for the 18 GBR samples. Our results indicated that SNPs affected the mature miRNA processing mechanism and contributed to the production of isomiRs. This mechanism may have important significance for epigenetic changes and diseases.

Methodological considerations for measuring biofluid-based microRNA biomarkers

MicroRNAs (miRNAs) are small non-coding RNA that regulate the expression of messenger RNA and are implicated in almost all cellular processes. Importantly, miRNAs can be released extracellularly and are stable in these matrices where they may serve as indicators of organ or cell-specific toxicity, disease, and biological status. There has thus been great enthusiasm for developing miRNAs as biomarkers of adverse outcomes for scientific, regulatory, and clinical purposes. Despite advances in measurement capabilities for miRNAs, miRNAs are still not routinely employed as noninvasive biomarkers. This is in part due to the lack of standard approaches for sample preparation and miRNA measurement and uncertainty in their biological interpretation. Members of the microRNA Biomarkers Workgroup within the Health and Environmental Sciences Institute's (HESI) Committee on Emerging Systems Toxicology for the Assessment of Risk (eSTAR) are a consortium of private- and public-sector scientists dedicated to developing miRNAs as applied biomarkers. Here, we explore major impediments to routine acceptance and use of miRNA biomarkers and case examples of successes and deficiencies in development. Finally, we provide insight on miRNA measurement, collection, and analysis tools to provide solid footing for addressing knowledge gaps toward routine biomarker use.

Extracellular vesicles piwi-interacting RNAs from skin mucus for identification of infected Cynoglossus semilaevis with Vibrio harveyi

Extracellular vesicles play a regulatory role in intracellular and intercellular transmission through a variety of biological information molecules, including mRNA, small RNAs and proteins. piRNAs are one kind of regulatory small RNAs in the vesicles at the post transcriptional level. Hereby, we isolated the extracellular vesicles from skin mucus and screened the piRNA profiles of these vesicles, aiming at developing biomarkers related to bacterial infections in Cynoglossus semilaevis. The different profilings of piRNAs in mucous extracellular vesicles of C. semilaevis were compared through small RNA sequencing, between fish infected with Vibrio harveyi and healthy ones. The number of clean reads on the alignment of exosome sick (ES) group was 105, 345 and that of exosome control (EC) group was 455, 144. GO and KEGG pathway enrichment analysis showed that most of the target genes were involved in cellular process, response to stimulus, biological regulation, immune system process and signal transduction, signal molecular and interaction, transport and catabolism. The 45 final candidate piRNAs related to immunity or infectious diseases included 20 piRNAs with high expression in the ES group and 25 piRNAs with a low expression in the ES group. After verification by qRT-PCR, there was significant difference of five piRNAs expression level between infected fish and healthy fish, in line with the sequencing. The expression level of piR-mmu-16401212, piR-mmu-26829319 and piR-gga-244092 in infected fish were significantly lower than that of control group, while piR-gga-71717 and piR-gga-99034 were higher, which implying that these piRNAs in mucous extracellular vesicles can be used to identify diseased fish from normal ones. This work supplied a novel class of biomarker for infection diagnosis in fish, and it will be benefit for screening disease resistant breeding of C. semilaevis.

IsomiR_Window: a system for analyzing small-RNA-seq data in an integrative and user-friendly manner

BACKGROUND

IsomiRs are miRNA variants that vary in length and/or sequence when compared to their canonical forms. These variants display differences in length and/or sequence, including additions or deletions of one or more nucleotides (nts) at the 5' and/or 3' end, internal editings or untemplated 3' end additions. Most available tools for small RNA-seq data analysis do not allow the identification of isomiRs and often require advanced knowledge of bioinformatics. To overcome this, we have developed IsomiR Window, a platform that supports the systematic identification, quantification and functional exploration of isomiR expression in small RNA-seq datasets, accessible to users with no computational skills.

METHODS

IsomiR Window enables the discovery of isomiRs and identification of all annotated non-coding RNAs in RNA-seq datasets from animals and plants. It comprises two main components: the IsomiR Window pipeline for data processing; and the IsomiR Window Browser interface. It integrates over ten third-party softwares for the analysis of small-RNA-seq data and holds a new algorithm that allows the detection of all possible types of isomiRs. These include 3' and 5'end isomiRs, 3' end tailings, isomiRs with single nucleotide polymorphisms (SNPs) or potential RNA editings, as well as all possible fuzzy combinations. IsomiR Window includes all required databases for analysis and annotation, and is freely distributed as a Linux virtual machine, including all required software.

RESULTS

IsomiR Window processes several datasets in an automated manner, without restrictions of input file size. It generates high quality interactive figures and tables which can be exported into different formats. The performance of isomiR detection and quantification was assessed using simulated small-RNA-seq data. For correctly mapped reads, it identified different types of isomiRs with high confidence and 100% accuracy. The analysis of a small RNA-seq data from Basal Cell Carcinomas (BCCs) using isomiR Window confirmed that miR-183-5p is up-regulated in Nodular BCCs, but revealed that this effect was predominantly due to a novel 5'end variant. This variant displays a different seed region motif and 1756 isoform-exclusive mRNA targets that are significantly associated with disease pathways, underscoring the biological relevance of isomiR-focused analysis. IsomiR Window is available at https://isomir.fc.ul.pt/ .

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.

Detecting Cynoglossus semilaevis infected with Vibrio harveyi using micro RNAs from mucous exosomes

Exosomes are important mediators of vesicle transportation and contain microRNAs (miRNAs) that mediate transcriptional gene knockout and silencing in biological processes. Moreover, exosomic miRNAs are promising biomarkers for disease diagnosis and physiological status indication in many species, including fish. The impact of the Vibrio harveyi pathogen on Cynoglossus semilaevis aquaculture is becoming more and more serious as the industry expands. To overcome this challenge, miRNAs in mucous exosomes were screened by small RNA sequencing and verified by quantitative real-time PCR to develop biomarkers. This is the first capture of exosomes from flatfish mucus coupled with miRNA profiling. The results revealed significant differences in expression levels of some miRNAs between infected and healthy fish. Three unique miRNAs were identified for V. harveyi infection diagnosis; expression levels of dre-miR-205-5p and dre-miR-205-5p in infected fish were significantly lower than controls, while dre-miR-100-5p expression was higher. These miRNAs in mucous exosomes could be used to differentiate diseased and healthy fish in an early screening method with practical value for breeding disease-resistant C. semilaevis.

An Approach to Identify Individual Functional Single Nucleotide Polymorphisms and Isoform MicroRNAs

MicroRNAs (miRNAs) and single nucleotide polymorphisms (SNPs) play important roles in disease risk and development, especially cancer. Importantly, when SNPs are located in pre-miRNAs, they affect their splicing mechanism and change the function of miRNAs. To improve disease risk assessment, we propose an approach and developed a software tool, IsomiR_Find, to identify disease/phenotype-related SNPs and isomiRs in individuals. Our approach is based on the individual's samples, with SNP information extracted from the 1000 Genomes Project. SNPs were mapped to pre-miRNAs based on whole-genome coordinates and then SNP-pre-miRNA sequences were constructed. Moreover, we developed matpred2, a software tool to identify the four splicing sites of mature miRNAs. Using matpred2, we identified isomiRs and then verified them by searching within individual miRNA sequencing data. Our approach yielded biomarkers for biological experiments, mined functions of miRNAs and SNPs, improved disease risk assessment, and provided a way to achieve individualized precision medicine.

miRNA analysis with Prost! reveals evolutionary conservation of organ-enriched expression and post-transcriptional modifications in three-spined stickleback and zebrafish

MicroRNAs (miRNAs) can have organ-specific expression and functions; they can originate from dedicated miRNA genes, from non-canonical miRNA genes, or from mirror-miRNA genes and can also experience post-transcriptional variation. It remains unclear, however, which mechanisms of miRNA production or modification are organ-specific and the extent of their evolutionary conservation. To address these issues, we developed the software Prost! (PRocessing Of Short Transcripts), which, among other features, helps quantify mature miRNAs, accounts for post-transcriptional processing, such as nucleotide editing, and identifies mirror-miRNAs. Here, we applied Prost! to annotate and analyze miRNAs in three-spined stickleback (Gasterosteus aculeatus), a model fish for evolutionary biology reported to have a miRNome larger than most teleost fish. Zebrafish (Danio rerio), a distantly related teleost with a well-known miRNome, served as comparator. Our results provided evidence for the existence of 286 miRNA genes and 382 unique mature miRNAs (excluding mir430 gene duplicates and the vaultRNA-derived mir733), which doesn't represent a miRNAome larger than other teleost miRNomes. In addition, small RNA sequencing data from brain, heart, testis, and ovary in both stickleback and zebrafish identified suites of mature miRNAs that display organ-specific enrichment, many of which are evolutionarily-conserved in the brain and heart in both species. These data also supported the hypothesis that evolutionarily-conserved, organ-specific mechanisms may regulate post-transcriptional variations in miRNA sequence. In both stickleback and zebrafish, miR2188-5p was edited frequently with similar nucleotide changes in the seed sequence with organ specific editing rates, highest in the brain. In summary, Prost! is a new tool to identify and understand small RNAs, to help clarify a species' miRNA biology as shown here for an important model for the evolution of developmental mechanisms, and to provide insight into organ-enriched expression and the evolutionary conservation of miRNA post-transcriptional modifications.

Genome-wide miRNA profiling reinforces the importance of miR-9 in human papillomavirus associated oral and oropharyngeal head and neck cancer

Head and neck cancer is the sixth most common malignancy worldwide, predominantly developing from squamous cell epithelia (HNSCC). The main HNSCC risk factors are tobacco, excessive alcohol use, and the presence of human papillomavirus (HPV). HPV positive (+) cancers are etiologically different from other HNSCC and often show better prognosis. The current knowledge regarding HNSCC miRNA profiles is still incomplete especially in the context of HPV+ cancer. Thus, we analyzed 61 freshly collected primary oral (OSCC) and oropharyngeal (OPSCC) SCC samples. HPV DNA and RNA was found in 21% cases. The Illumina whole-genome small-RNA profiling by next-generation sequencing was done on 22 samples and revealed 7 specific miRNAs to HPV+ OSCC, 77 to HPV+ OPSCC, and additional 3 shared with both; 51 miRNAs were specific to HPV− OPSCC, 62 to HPV− OSCC, and 31 shared with both. The results for 9 miRNAs (miR-9, -21, -29a, -100, -106b, -143 and -145) were assessed by reverse transcription-quantitative polymerase chain reaction on the whole study population. The data was additionally confirmed by reanalyzing publicly available miRNA sequencing Cancer Genome Atlas consortium (TCGA) HNSCC data. Cell signaling pathway analysis revealed differences between HPV+ and HPV− HNSCC. Our findings compared with literature data revealed extensive heterogeneity of miRNA deregulation with only several miRNAs consistently affected, and miR-9 being the most likely HPV related miRNA.

3' RNA Uridylation in Epitranscriptomics, Gene Regulation, and Disease

Emerging evidence implicates a wide range of post-transcriptional RNA modifications that play crucial roles in fundamental biological processes including regulating gene expression. Collectively, they are known as epitranscriptomics. Recent studies implicate 3' RNA uridylation, the non-templated addition of uridine(s) to the terminal end of RNA, as a key player in epitranscriptomics. In this review, we describe the functional roles and significance of 3' terminal RNA uridylation that has diverse functions in regulating both mRNAs and non-coding RNAs. In mammals, three Terminal Uridylyl Transferases (TUTases) are primarily responsible for 3' RNA uridylation. These enzymes are also referred to as polyU polymerases. TUTase 1 (TUT1) is implicated in U6 snRNA maturation via uridylation. The TUTases TUT4 and/or TUT7 are the predominant mediators of all other cellular uridylation. Terminal uridylation promotes turnover for many polyadenylated mRNAs, replication-dependent histone mRNAs that lack polyA-tails, and aberrant structured noncoding RNAs. In addition, uridylation regulates biogenesis of a subset of microRNAs and generates isomiRs, sequent variant microRNAs that have altered function in specific cases. For example, the RNA binding protein and proto-oncogene LIN28A and TUT4 work together to polyuridylate pre-let-7, thereby blocking biogenesis and function of the tumor suppressor let-7 microRNA family. In contrast, monouridylation of Group II pre-miRNAs creates an optimal 3' overhang that promotes recognition and subsequent cleavage by the Dicer-TRBP complex that then yields the mature microRNA. Also, uridylation may play a role in non-canonical microRNA biogenesis. The overall significance of 3' RNA uridylation is discussed with an emphasis on mammalian development, gene regulation, and disease, including cancer and Perlman syndrome. We also introduce recent changes to the HUGO-approved gene names for multiple terminal nucleotidyl transferases that affects in part TUTase nomenclature (TUT1/TENT1, TENT2/PAPD4/GLD2, TUT4/ZCCHC11/TENT3A, TUT7/ZCCHC6/TENT3B, TENT4A/PAPD7, TENT4B/PAPD5, TENT5A/FAM46A, TENT5B/FAM46B, TENT5C/FAM46C, TENT5D/FAM46D, MTPAP/TENT6/PAPD1).

Evidence for plant-derived xenomiRs based on a large-scale analysis of public small RNA sequencing data from human samples

In recent years, an increasing number of studies have reported the presence of plant miRNAs in human samples, which resulted in a hypothesis asserting the existence of plant-derived exogenous microRNA (xenomiR). However, this hypothesis is not widely accepted in the scientific community due to possible sample contamination and the small sample size with lack of rigorous statistical analysis. This study provides a systematic statistical test that can validate (or invalidate) the plant-derived xenomiR hypothesis by analyzing 388 small RNA sequencing data from human samples in 11 types of body fluids/tissues. A total of 166 types of plant miRNAs were found in at least one human sample, of which 14 plant miRNAs represented more than 80% of the total plant miRNAs abundance in human samples. Plant miRNA profiles were characterized to be tissue-specific in different human samples. Meanwhile, the plant miRNAs identified from microbiome have an insignificant abundance compared to those from humans, while plant miRNA profiles in human samples were significantly different from those in plants, suggesting that sample contamination is an unlikely reason for all the plant miRNAs detected in human samples. This study also provides a set of testable synthetic miRNAs with isotopes that can be detected in situ after being fed to animals.

Estradiol and the Development of the Cerebral Cortex: An Unexpected Role?

The cerebral cortex undergoes rapid folding in an "inside-outside" manner during embryonic development resulting in the establishment of six discrete cortical layers. This unique cytoarchitecture occurs via the coordinated processes of neurogenesis and cell migration. In addition, these processes are fine-tuned by a number of extracellular cues, which exert their effects by regulating intracellular signaling pathways. Interestingly, multiple brain regions have been shown to develop in a sexually dimorphic manner. In many cases, estrogens have been demonstrated to play an integral role in mediating these sexual dimorphisms in both males and females. Indeed, 17β-estradiol, the main biologically active estrogen, plays a critical organizational role during early brain development and has been shown to be pivotal in the sexually dimorphic development and regulation of the neural circuitry underlying sex-typical and socio-aggressive behaviors in males and females. However, whether and how estrogens, and 17β-estradiol in particular, regulate the development of the cerebral cortex is less well understood. In this review, we outline the evidence that estrogens are not only present but are engaged and regulate molecular machinery required for the fine-tuning of processes central to the cortex. We discuss how estrogens are thought to regulate the function of key molecular players and signaling pathways involved in corticogenesis, and where possible, highlight if these processes are sexually dimorphic. Collectively, we hope this review highlights the need to consider how estrogens may influence the development of brain regions directly involved in the sex-typical and socio-aggressive behaviors as well as development of sexually dimorphic regions such as the cerebral cortex.

Which statistical significance test best detects oncomiRNAs in cancer tissues? An exploratory analysis

MicroRNAs(miRNAs) often exert their oncogenic and tumor suppressor functions by suppressing protein-coding genes expressions in cancers and thus have a strong association with cancers' generation, development and metastasis. Through comprehensively understanding differentially expressed miRNAs (oncomiRNA) in tumor tissues, we can elucidate the underlying molecular mechanisms in tumorigenesis and develop novel strategies for cancer diagnosis and treatment. The differential expression of miRNAs can now be analyzed through numerous statistical significance tests based on different principles, which are also available in various R packages. However, the results can be notably different. In this study, we compared miRNAs obtained from 6 common significance tests/R packages (t-test, Limma, DESeq, edgeR, LRT and MARS) with the miRNAs archived in two databases; HMDD 2.0 database, which collects experimentally validated differentially expressed miRNAs, and Infer microRNA-disease association database, which contains the potential disease-associated miRNAs by network forecasting. Finally, we sought the MARS method in DEGseq package more effectively searched out differentially expressed miRNAs than other common methods.

MiR-93-5p Promotes Cell Proliferation through Down-Regulating PPARGC1A in Hepatocellular Carcinoma Cells by Bioinformatics Analysis and Experimental Verification

Peroxisome proliferator-activated receptor gamma coactivator-1 alpha (PPARGC1A, formerly known as PGC-1a) is a transcriptional coactivator and metabolic regulator. Previous studies are mainly focused on the association between PPARGC1A and hepatoma. However, the regulatory mechanism remains unknown. A microRNA associated with cancer (oncomiR), miR-93-5p, has recently been found to play an essential role in tumorigenesis and progression of various carcinomas, including liver cancer. Therefore, this paper aims to explore the regulatory mechanism underlying these two proteins in hepatoma cells. Firstly, an integrative analysis was performed with miRNA–mRNA modules on microarray and The Cancer Genome Atlas (TCGA) data and obtained the core regulatory network and miR-93-5p/PPARGC1A pair. Then, a series of experiments were conducted in hepatoma cells with the results including miR-93-5p upregulated and promoted cell proliferation. Thirdly, the inverse correlation between miR-93-5p and PPARGC1A expression was validated. Finally, we inferred that miR-93-5p plays an essential role in inhibiting PPARGC1A expression by directly targeting the 3′-untranslated region (UTR) of its mRNA. In conclusion, these results suggested that miR-93-5p overexpression contributes to hepatoma development by inhibiting PPARGC1A. It is anticipated to be a promising therapeutic strategy for patients with liver cancer in the future.

Plant IsomiR Atlas: Large Scale Detection, Profiling, and Target Repertoire of IsomiRs in Plants

microRNAs (miRNAs) play an important role as key regulators controlling the post-transcriptional events in plants across development, abiotic and biotic stress, tissue polarity and also in defining the evolutionary basis of the origin of the post-transcriptional machinery. Identifying patterns of regulated and co-regulated small RNAs, in particular miRNAs and their sequence variants with the availability of next generation sequencing approaches has widely demonstrated the role of miRNAs and their temporal regulation in maintaining plant development and their response to stress conditions. Although the role of canonical miRNAs has been widely explored and functional diversity is revealed, those works for isomiRs are still limited and urgent to be carried out across plants. This relative lack of information with respect to isomiRs might be attributed to the non-availability of large-scale detection of isomiRs across wide plant species. In the present research, we addressed this by developing Plant isomiR Atlas, which provides large-scale detection of isomiRs across 23 plant species utilizing 677 smallRNAs datasets and reveals a total of 98,374 templated and non-templated isomiRs from 6,167 precursors. Plant isomiR Atlas provides several visualization features such as species specific isomiRs, isomiRs and canonical miRNAs overlap, terminal modification classifications, target identification using psRNATarget and TargetFinder and also canonical miRNAs:target interactions. Plant isomiR Atlas will play a key role in understanding the regulatory nature of miRNAome and will accelerate to understand the functional role of isomiRs. Plant isomiR Atlas is available at www.mcr.org.in/isomir. One Sentence Summary  Plant isomiR Atlas will play a key role in understanding the regulatory nature of miRNAome and will accelerate the understanding and diversity of functional targets of plants isomiRs.

Web-based NGS data analysis using miRMaster: a large-scale meta-analysis of human miRNAs

The analysis of small RNA NGS data together with the discovery of new small RNAs is among the foremost challenges in life science. For the analysis of raw high-throughput sequencing data we implemented the fast, accurate and comprehensive web-based tool miRMaster. Our toolbox provides a wide range of modules for quantification of miRNAs and other non-coding RNAs, discovering new miRNAs, isomiRs, mutations, exogenous RNAs and motifs. Use-cases comprising hundreds of samples are processed in less than 5 h with an accuracy of 99.4%. An integrative analysis of small RNAs from 1836 data sets (20 billion reads) indicated that context-specific miRNAs (e.g. miRNAs present only in one or few different tissues / cell types) still remain to be discovered while broadly expressed miRNAs appear to be largely known. In total, our analysis of known and novel miRNAs indicated nearly 22 000 candidates of precursors with one or two mature forms. Based on these, we designed a custom microarray comprising 11 872 potential mature miRNAs to assess the quality of our prediction. MiRMaster is a convenient-to-use tool for the comprehensive and fast analysis of miRNA NGS data. In addition, our predicted miRNA candidates provided as custom array will allow researchers to perform in depth validation of candidates interesting to them.

Small RNAs in Circulating Exosomes of Cancer Patients: A Minireview

Extracellular vesicles (EVs) secreted from many cell types play important roles in intercellular communication, both as paracrine and endocrine factors, as they can circulate in biological fluids, including plasma. Amid EVs, exosomes are actively secreted vesicles that contain proteins, lipids, soluble factors, and nucleic acids, including microRNAs (miRNAs) and other classes of small RNAs (sRNA). miRNAs are prominent post-transcriptional regulators of gene expression and epigenetic silencers of transcription. We concisely review the roles of miRNAs in cell-fate determination and development and their regulatory activity on almost all the processes and pathways controlling tumor formation and progression. Next, we consider the evidence linking exosomes to tumor progression, particularly to the setting-up of permissive pre-metastatic niches. The study of exosomes in patients with different survival and therapy response can inform on the possible correlations between exosomal cargo and disease features. Moreover, the exploration of circulating exosomes as possible sources of non-invasive biomarkers could give new implements for anti-cancer therapy and metastasis prevention. Since the characterization of sRNAs in exosomes of cancer patients sparks opportunities to better understand their roles in cancer, we briefly present current experimental and computational protocols for sRNAs analysis in circulating exosomes by RNA-seq.

Evaluation of high-throughput isomiR identification tools: illuminating the early isomiRome of Tribolium castaneum

Background

MicroRNAs carry out post-transcriptional gene regulation in animals by binding to the 3' untranslated regions of mRNAs, causing their degradation or translational repression. MicroRNAs influence many biological functions, and dysregulation can therefore disrupt development or even cause death. High-throughput sequencing and the mining of animal small RNA data has shown that microRNA genes can yield differentially expressed isoforms, known as isomiRs. Such isoforms are particularly relevant during early development, and the extension or truncation of the 5' end can change the profile of mRNA targets compared to the original mature sequence. We used the publicly available small RNA dataset of the model beetle Tribolium castaneum to create the first comparative isomiRome of early developmental stages in this species. Standard microRNA analysis software does not specifically account for isomiRs. We therefore carried out the first comparative evaluation of the specialized tools isomiRID, isomiR-SEA and miraligner, which can be downloaded for local use and can handle next generation sequencing data.

Results

We compared the performance of isomiRID, isomiR-SEA and miraligner using simulated Illumina HiSeq2000 and MiSeq data to test the impact of technical errors. We also created artificial microRNA isoforms to determine the effect of biological variants on the performance of each algorithm. We found that isomiRID achieved the best true positive rate among the three algorithms, but only accounted for one mutation at a time. In contrast, miraligner reported all variations simultaneously but with 78% sensitivity, yielding isomiRs with 3' or 5' deletions. Finally, isomiR-SEA achieved a sensitivity of 25–33% when the seed region was mutated or partly deleted, but was the only tool that could accommodate more than one mismatch. Using the best tool, we performed a complete isomiRome analysis of the early developmental stages of T. castaneum.

Conclusions

Our findings will help researchers to select the most suitable isomiR analysis tools for their experiments. We confirmed the dynamic expression of 3′ non-template isomiRs and expanded the isomiRome by all known isomiR modifications during the early development of T. castaneum.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1772-z) contains supplementary material, which is available to authorized users.

Identification of exosomes and its signature miRNAs of male and female Cynoglossus semilaevis

Exosomes are small membrane particles which are widely found in various cell lines and physiological fluids in mammalian. MicroRNAs (miRNAs) enclosed in exosomes have been identified as proper signatures for many diseases and response to therapies. However, the composition of exosomes and enclosed miRNAs in fishes has not been investigated. Cynoglossus semilaevis is an important commercial flatfish with ambiguous distinction between males and females before sex maturation, which leads to screening difficulty in reproduction and cultivation. An effective detection method was required for sex differentiation of C. semilaevis. In this work, we successfully identified exosomes in C. semilaevis serum. The analysis of nucleotide composition showed that miRNA dominated in exosomes. Thereafter the miRNA profiles in exosomes from males and females were sequenced and compared to identify the signature miRNAs corresponding to sex differentiation. The functions of signature miRNAs were analyzed by target matching and annotation. Furthermore, 7 miRNAs with high expression in males were selected from signature miRNAs as the markers for sex identification with their expression profiles verified by real time quantitative PCR. Exosomes were first found in fish serum in this work. Investigation of marker miRNAs supplies an effective index for the filtration of male and female C. semilaevis in cultivation.

miRNA-based heavy metal homeostasis and plant growth

Plants have been naturally gifted with mechanisms to adjust under very high or low nutrient concentrations. Heavy metal toxicity is considered as a major growth and yield-limiting factor for plants. This stress includes essential as well as non-essential metals. MicroRNAs (miRNAs) are known for mediating post-transcriptional regulation by cleaving transcripts or translational inhibition. It is commonly agreed that an extensive understanding of plant miRNAs will significantly help in the induction of tolerance against environmental stresses. With the introduction of the latest technology like next generation sequencing (NGS), a growing figure of miRNAs has been productively recognized in several plants for their diverse roles. These miRNAs are well-known modulators of plant responses to heavy metal (HM) stress. Data regarding metal-responsive miRNAs point out the vital role of plant miRNAs in supplementing metal detoxification by means of transcription factors (TF) or gene regulation. Acting as systemic signals, miRNAs also synchronize different physiological processes for plant responses to metal toxicities. In contrast to practicing techniques, using miRNA is a greatly helpful, pragmatic, and feasible approach. The earlier findings point towards miRNAs as a prospective target to engineer heavy metal tolerance in plants. Therefore, there is a need to augment our knowledge about the orchestrated functions of miRNAs during HM stress. We reviewed the deterministic significance of plant miRNAs in heavy metal tolerance and their role in mediating plant responses to HM toxicities. This review also summarized the topical developments by identification and validation of different metal stress-responsive miRNAs.

Extracting microRNA-gene relations from biomedical literature using distant supervision

Many biomedical relation extraction approaches are based on supervised machine learning, requiring an annotated corpus. Distant supervision aims at training a classifier by combining a knowledge base with a corpus, reducing the amount of manual effort necessary. This is particularly useful for biomedicine because many databases and ontologies have been made available for many biological processes, while the availability of annotated corpora is still limited. We studied the extraction of microRNA-gene relations from text. MicroRNA regulation is an important biological process due to its close association with human diseases. The proposed method, IBRel, is based on distantly supervised multi-instance learning. We evaluated IBRel on three datasets, and the results were compared with a co-occurrence approach as well as a supervised machine learning algorithm. While supervised learning outperformed on two of those datasets, IBRel obtained an F-score 28.3 percentage points higher on the dataset for which there was no training set developed specifically. To demonstrate the applicability of IBRel, we used it to extract 27 miRNA-gene relations from recently published papers about cystic fibrosis. Our results demonstrate that our method can be successfully used to extract relations from literature about a biological process without an annotated corpus. The source code and data used in this study are available at https://github.com/AndreLamurias/IBRel.

IsomiR expression patterns in canonical and Dicer‑independent microRNAs

Multiple microRNA (miRNA) variants, known as isomiRs, are extensively distributed in miRNA loci and predominantly derive from the alternative cleavage of Drosha/Dicer and 3′addition events. The present study aimed to investigate the expression patterns of multiple isomiRs in typical miRNA and Dicer-independent miRNA loci by conducting evolutionary and expression analysis using public datasets. Although different miRNA maturation processes exist, multiple isomiRs can be detected by similar expression distributions. However, isomiR expression in Dicer-independent miRNA loci tends to be at a moderate level, particularly for random distribution in the ends that are split by Dicer in the typical miRNA loci. Compared with the mature miRNA locus (dominant miRNA locus), the non-dominant miRNA locus indicates an expression distribution similar to that of the Dicer-independent miRNA locus. These results increase the understanding of multiple isomiRs in the progression of diseases.

Computational identification of putative lincRNAs in mouse embryonic stem cell

As the regulatory factors, lncRNAs play critical roles in embryonic stem cells. And lincRNAs are most widely studied lncRNAs, however, there might still might exist a large member of uncovered lncRNAs. In this study, we constructed the de novo assembly of transcriptome to detect 6,701 putative long intergenic non-coding transcripts (lincRNAs) expressed in mouse embryonic stem cells (ESCs), which might be incomplete with the lack coverage of 5' ends assessed by CAGE peaks. Comparing the TSS proximal regions between the known lincRNAs and their closet protein coding transcripts, our results revealed that the lincRNA TSS proximal regions are associated with the characteristic genomic and epigenetic features. Subsequently, 1,293 lincRNAs were corrected at their 5' ends using the putative lincRNA TSS regions predicted by the TSS proximal region prediction model based on genomic and epigenetic features. Finally, 43 putative lincRNAs were annotated by Gene Ontology terms. In conclusion, this work provides a novel catalog of mouse ESCs-expressed lincRNAs with the relatively complete transcript length, which might be useful for the investigation of transcriptional and post-transcriptional regulation of lincRNA in mouse ESCs and even mammalian development.

lncScore: alignment-free identification of long noncoding RNA from assembled novel transcripts

RNA-Seq based transcriptome assembly has been widely used to identify novel lncRNAs. However, the best-performing transcript reconstruction methods merely identified 21% of full-length protein-coding transcripts from H. sapiens. Those partial-length protein-coding transcripts are more likely to be classified as lncRNAs due to their incomplete CDS, leading to higher false positive rate for lncRNA identification. Furthermore, potential sequencing or assembly error that gain or abolish stop codons also complicates ORF-based prediction of lncRNAs. Therefore, it remains a challenge to identify lncRNAs from the assembled transcripts, particularly the partial-length ones. Here, we present a novel alignment-free tool, lncScore, which uses a logistic regression model with 11 carefully selected features. Compared to other state-of-the-art alignment-free tools (e.g. CPAT, CNCI, and PLEK), lncScore outperforms them on accurately distinguishing lncRNAs from mRNAs, especially partial-length mRNAs in the human and mouse datasets. In addition, lncScore also performed well on transcripts from five other species (Zebrafish, Fly, C. elegans, Rat, and Sheep). To speed up the prediction, multithreading is implemented within lncScore, and it only took 2 minute to classify 64,756 transcripts and 54 seconds to train a new model with 21,000 transcripts with 12 threads, which is much faster than other tools. lncScore is available at https://github.com/WGLab/lncScore.

BP Neural Network Could Help Improve Pre-miRNA Identification in Various Species

MicroRNAs (miRNAs) are a set of short (21–24 nt) noncoding RNAs that play significant regulatory roles in cells. In the past few years, research on miRNA-related problems has become a hot field of bioinformatics because of miRNAs' essential biological function. miRNA-related bioinformatics analysis is beneficial in several aspects, including the functions of miRNAs and other genes, the regulatory network between miRNAs and their target mRNAs, and even biological evolution. Distinguishing miRNA precursors from other hairpin-like sequences is important and is an essential procedure in detecting novel microRNAs. In this study, we employed backpropagation (BP) neural network together with 98-dimensional novel features for microRNA precursor identification. Results show that the precision and recall of our method are 95.53% and 96.67%, respectively. Results further demonstrate that the total prediction accuracy of our method is nearly 13.17% greater than the state-of-the-art microRNA precursor prediction software tools.

Dynamic mRNA and miRNA expression analysis in response to intermuscular bone development of blunt snout bream (Megalobrama amblycephala)

Intermuscular bone (IB), which occurs only in the myosepta of lower teleosts, is attracting more attention because they are difficult to remove and make the fish unpleasant to eat. By gaining a better understanding of the genetic regulation of IB development, an integrated analysis of miRNAs and mRNAs expression profiling was performed on Megalobrama amblycephala. Four key development stages were selected for transcriptome and small RNA sequencing. A number of significantly differentially expressed miRNAs/genes associated with bone formation and differentiation were identified and the functional characteristics of these miRNAs/genes were revealed by GO function and KEGG pathway analysis. These were involved in TGF-β, ERK and osteoclast differentiation pathways known in the literature to affect bone formation and differentiation. MiRNA-mRNA interaction pairs were detected from comparison of expression between different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating bone development and differentiation. A negative regulation effect of two miRNAs was verified through dual luciferase reporter assay. As a unique public resource for gene expression and regulation during the IB development, this study is expected to provide forwards ideas and resources for further biological researches to understand the IBs’ development.