Search

Simple Summary: MicroRNAs (miRNAs) play a significant role as epigenetic regulators in cancer. IsomiRs are miRNA molecules that undergo small modifications during miRNA processing, which can affect their stability and their interaction with mRNA targets. While some isomiRs are linked to specific cancers, many, including those in the lung, remain understudied. To address this, small RNA sequencing data from lung adenocarcinoma (LUAD) and adult non-malignant lung (ANL) samples were analyzed to quantify isomiR expression. This analysis identified 16 A-to-I edited isomiRs, 213 5′ isomiRs, 128 3′ adenylated isomiRs, and 100 3′ uridylated isomiRs. A-to-I editing rates correlated with the expression of editing enzymes ADAR and ADARB1, both deregulated in LUAD. LUAD samples had lower A-to-I editing and 3′ adenylation rates compared to ANL. Machine learning models based on isomiR data effectively distinguished ANL from stage I/II LUAD samples, suggesting that isomiRs hold potential as cancer biomarkers. Background: MicroRNAs (miRNAs) modulate the expression of oncogenes and tumor suppressor genes, functioning as significant epigenetic regulators in cancer. IsomiRs are miRNA molecules that have undergone small modifications during miRNA processing. These modifications can alter an isomiR's binding stability with mRNA targets, and certain isomiRs have been implicated in the development of specific cancers. Still, the isomiRomes of many tissues, including the lung, have not been characterized; Methods: In this study, we analyzed small RNA sequencing data for three cohorts of lung adenocarcinoma (LUAD) and adult non-malignant lung (ANL) samples. Results: We quantified isomiR expression and found 16 A-to-I edited isomiRs expressed in multiple cohorts, as well as 213 5′ isomiRs, 128 3′ adenylated isomiRs, and 100 3′ uridylated isomiRs. Rates of A-to-I editing at editing hotspots correlated with mRNA expression of the editing enzymes ADAR and ADARB1, which were both observed to be deregulated in LUAD. LUAD samples displayed lower overall rates of A-to-I editing and 3′ adenylation than ANL samples. Support vector machines and random forest models were trained on one cohort to distinguish ANL and stage I/II LUAD samples using reads per million (RPM) and frequency data for different types of isomiRs. Models trained on A-to-I editing rates at editing hotspots displayed high accuracy when tested on the other two cohorts and compared favorably to classifiers trained on miRNA expression alone; Conclusions: We have identified isomiRs in the human lung and found that their expression differs between non-malignant and tumor tissues, suggesting they hold potential as cancer biomarkers. [ABSTRACT FROM AUTHOR]

The commonly used isomiR format (mirGFF3) combined with classification guidelines has been proposed to facilitate communication within the field. New technological and computational advances have validated certain isomiRs as biologically relevant. Various experimental studies have shown mediation of different regulation mechanisms through isomiRs as opposed to regulation mediated via archetypal miRNAs. In cancer, isomiRs have been identified as promising biomarkers in solid tissue for prognosis, and in circulatory fluids for diagnosis. Additionally, an isomiR has been identified as a promising therapeutic target. Similar studies are under way in the context of neurodegenerative, cardiovascular, and metabolic diseases. MicroRNAs (miRNAs) and isoforms of their archetype, called isomiRs, regulate gene expression via complementary base-pair binding to messenger RNAs (mRNAs). The partially evolutionarily conserved isomiR sequence variations are differentially expressed among tissues, populations, and genders, and between healthy and diseased states. Aiming towards the clinical use of isomiRs as diagnostic biomarkers and for therapeutic purposes, several challenges need to be addressed, including (i) clarification of isomiR definition, (ii) improved annotation in databases with new standardization (such as the mirGFF3 format), and (iii) improved methods of isomiR detection, functional verification, and in silico analysis. In this review we discuss the respective challenges, and highlight the opportunities for clinical use of isomiRs, especially in the light of increasing amounts of next-generation sequencing (NGS) data. [ABSTRACT FROM AUTHOR]

The analysis of small RNA sequencing data across a range of biofluids is a significant research area, given the diversity of RNA types that hold potential diagnostic, prognostic, and predictive value. The intricate task of segregating the complex mixture of small RNAs from both human and other species, including bacteria, fungi, and viruses, poses one of the most formidable challenges in the analysis of small RNA sequencing data, currently lacking satisfactory solutions. This study introduces sRNAflow, a user-friendly bioinformatic tool with a web interface designed for the analysis of small RNAs obtained from biological fluids. Tailored to the unique requirements of such samples, the proposed pipeline addresses various challenges, including filtering potential RNAs from reagents and environment, classifying small RNA types, managing small RNA annotation overlap, conducting differential expression assays, analysing isomiRs, and presenting an approach to identify the sources of small RNAs within samples. sRNAflow also encompasses an alternative alignment-free analysis of RNA-seq data, featuring clustering and initial RNA source identification using BLAST. This comprehensive approach facilitates meaningful comparisons of results between different analytical methods. [ABSTRACT FROM AUTHOR]

MiR-125b has therapeutic potential in the amelioration of myocardial ischemic injury. MicroRNA isomiRs, with either 5′ or 3′ addition or deletion of nucleotide(s), have been reported from next-generation sequencing data (NGS). However, due to technical challenges, validation and functional studies of isomiRs are few. In this study, we discovered using NGS, four 3′isomiRs of miR-125b, i.e., addition of A (adenosine), along with deletions of A, AG (guanosine) and AGU (uridine) from rat and sheep heart. These findings were validated using RT-qPCR. Comprehensive functional studies were carried out in the H9C2 hypoxia model. After miR-125b, isomiRs of Plus A, Trim A, AG and AGU mimic transfection, the H9C2 cells were subjected to hypoxic challenge. As assessed using cell viability, apoptosis, CCK-8 and LDH release, miR-125b and isomiRs were all protective against hypoxia. However, Plus A and Trim A were more effective than miR-125b, whilst Trim AG and Trim AGU had far weaker effects than miR-125b. Interestingly, both the gene regulation profile and apoptotic gene validation indicated a major overlap among miR-125b, Plus A and Trim A, whilst Trims AG and AGU revealed a different profile compared to miR-125b. Conclusions: miR-125b and its 3′ isomiRs are expressed stably in the heart. miR-125b and isomiRs with addition or deletion of A might function concurrently and concordantly under specific physiological and pathophysiological conditions. In-depth understanding of isomiRs' metabolism and function will contribute to better miRNA therapeutic drug design. [ABSTRACT FROM AUTHOR]

Simple Summary: This study investigated the cis-regulation of isomiR biogenesis in human lymphoblastoid cell lines. A total of 95 SNP–isomiR pairs demonstrated significant associations between SNPs and 5′-isomiRs, including base substitutions, trimmings, extensions, and additions. Notably, the study identified an association between rs6505162 and the 5′-extension of hsa-miR-423-3p, as well as the 5′-trimming of hsa-miR-423-5p. Additionally, the correlation of isomiR expression with breast cancer status in the TCGA dataset provided valuable insights into the genetic association with breast cancer tumorigenesis. The study also highlighted that canonical miRNAs may not be the most abundant isomiRs in human lymphoblastoid cell lines, emphasizing the role of isomiRs in biological processes. Furthermore, the presence of the allele-specific expression of miRNAs suggests the involvement of genetic variants in miRNA regulation. MicroRNAs play a critical role in regulating gene expression post-transcriptionally. Variations in mature microRNA sequences, known as isomiRs, arise from imprecise cleavage and nucleotide substitution or addition. These isomiRs can target different mRNAs or compete with their canonical counterparts, thereby expanding the scope of miRNA post-transcriptional regulation. Our study investigated the relationship between cis-acting single-nucleotide polymorphisms (SNPs) in precursor miRNA regions and isomiR composition, represented by the ratio of a specific 5′-isomiR subtype to all isomiRs identified for a particular mature miRNA. Significant associations between 95 SNP–isomiR pairs were identified. Of note, rs6505162 was significantly associated with both the 5′-extension of hsa-miR-423-3p and the 5′-trimming of hsa-miR-423-5p. Comparison of breast cancer and normal samples revealed that the expression of both isomiRs was significantly higher in tumors than in normal tissues. This study sheds light on the genetic regulation of isomiR maturation and advances our understanding of post-transcriptional regulation by microRNAs. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs), important regulators of gene expression, play critical roles in various biological processes and tumorigenesis. To reveal the potential relationships between multiple isomiRs and arm switching, we performed a comprehensive pan-cancer analysis to discuss their roles in tumorigenesis and cancer prognosis. Our results showed that many miR-#-5p and miR-#-3p pairs from the two arms of pre-miRNA may have abundant expression levels, and they are often involved in distinct functional regulatory networks by targeting different mRNAs, although they may also interact with common targets. The two arms may show diverse isomiR expression landscapes, and their expression ratio might vary, mainly depending on tissue type. Dominantly expressed isomiRs can be used to determine distinct cancer subtypes that are associated with clinical outcome, indicating that they may be potential prognostic biomarkers. Our findings indicate robust and flexible isomiR expression landscapes that will enrich the study of miRNAs/isomiRs and aid in revealing the potential roles of multiple isomiRs yielded by arm switching in tumorigenesis. [ABSTRACT FROM AUTHOR]

MiRNA isoforms (isomiRs) are single stranded small RNAs originating from the same pri-miRNA hairpin as a result of cleavage by Drosha and Dicer enzymes. Variations at the 5ʹ-end of a miRNA alter the seed region of the molecule, thus affecting the targetome of the miRNA. In this manuscript, we analysed the distribution of miRNA cleavage positions across 31 different cancers using miRNA sequencing data of TCGA project. As a result, we found that the processing positions are not tissue specific and that all miRNAs could be correctly classified as ones exhibiting homogeneous or heterogeneous cleavage at one of the four cleavage sites. In 42% of cases (42 out of 100 miRNAs), we observed imprecise 5ʹ-end Dicer cleavage, while this fraction was only 14% for Drosha (14 out of 99). To the contrary, almost all cleavage sites of 3ʹ-ends (either Drosha or Dicer) were heterogeneous. With the use of only four nucleotides surrounding a 5ʹ-end Dicer cleavage position we built a model which allowed us to distinguish between homogeneous and heterogeneous cleavage with the reliable quality (ROC AUC = 0.68). Finally, we showed the possible applications of the study by the analysis of two 5ʹ-end isoforms originating from the same exogeneous shRNA hairpin. It turned out that the less expressed shRNA variant was functionally active, which led to the increased off-targeting. Thus, the obtained results could be applied to the design of shRNAs whose processing will result in a single 5ʹ-variant. [ABSTRACT FROM AUTHOR]

A widely used procedure for selecting significant miRNA-mRNA or isomiR-mRNA pairs out of predicted interactions involves calculating the correlation between expression levels of miRNAs/isomiRs and mRNAs in a series of samples. In this manuscript, we aimed to assess the validity of this procedure by comparing isomiR-mRNA correlation profiles in sets of sequence-based predicted target mRNAs and non-target mRNAs (negative controls). Target prediction was carried out using RNA22 and TargetScan algorithms. Spearman’s correlation analysis was conducted using miRNA and mRNA sequencing data of The Cancer Genome Atlas Breast Invasive Carcinoma (TCGA-BRCA) project. Luminal A, luminal B, basal-like breast cancer subtypes, and adjacent normal tissue samples were analyzed separately. Using the sets of putative targets and non-targets, we introduced adjusted isomiR targeting activity (ITA)—the number of negatively correlated potential isomiR targets adjusted by the background (estimated using non-target mRNAs). We found that for most isomiRs a significant negative correlation between isomiR-mRNA expression levels appeared more often in a set of predicted targets compared to the non-targets. This trend was detected for both classical seed region binding types (8mer, 7mer-m8, 7mer-A1, 6mer) predicted by TargetScan and the non-classical ones (G:U wobbles and up to one mismatch or unpaired nucleotide within seed sequence) predicted by RNA22. Adjusted ITA distributions were similar for target sites located in 3′-UTRs and coding mRNA sequences, while 5′-UTRs had much lower scores. Finally, we observed strong cancer subtype-specific patterns of isomiR activity, highlighting the differences between breast cancer molecular subtypes and normal tissues. Surprisingly, our target prediction- and correlation-based estimates of isomiR activities were practically non-correlated with the average isomiR expression levels neither in cancerous nor in normal samples. [ABSTRACT FROM AUTHOR]

MiR-125b has therapeutic potential in the amelioration of myocardial ischemic injury. MicroRNA isomiRs, with either 5′ or 3′ addition or deletion of nucleotide(s), have been reported from next-generation sequencing data (NGS). However, due to technical challenges, validation and functional studies of isomiRs are few. In this study, we discovered using NGS, four 3′isomiRs of miR-125b, i.e., addition of A (adenosine), along with deletions of A, AG (guanosine) and AGU (uridine) from rat and sheep heart. These findings were validated using RT-qPCR. Comprehensive functional studies were carried out in the H9C2 hypoxia model. After miR-125b, isomiRs of Plus A, Trim A, AG and AGU mimic transfection, the H9C2 cells were subjected to hypoxic challenge. As assessed using cell viability, apoptosis, CCK-8 and LDH release, miR-125b and isomiRs were all protective against hypoxia. However, Plus A and Trim A were more effective than miR-125b, whilst Trim AG and Trim AGU had far weaker effects than miR-125b. Interestingly, both the gene regulation profile and apoptotic gene validation indicated a major overlap among miR-125b, Plus A and Trim A, whilst Trims AG and AGU revealed a different profile compared to miR-125b. Conclusions: miR-125b and its 3′ isomiRs are expressed stably in the heart. miR-125b and isomiRs with addition or deletion of A might function concurrently and concordantly under specific physiological and pathophysiological conditions. In-depth understanding of isomiRs’ metabolism and function will contribute to better miRNA therapeutic drug design.

A widely used procedure for selecting significant miRNA-mRNA or isomiR-mRNA pairs out of predicted interactions involves calculating the correlation between expression levels of miRNAs/isomiRs and mRNAs in a series of samples. In this manuscript, we aimed to assess the validity of this procedure by comparing isomiR-mRNA correlation profiles in sets of sequence-based predicted target mRNAs and non-target mRNAs (negative controls). Target prediction was carried out using RNA22 and TargetScan algorithms. Spearman’s correlation analysis was conducted using miRNA and mRNA sequencing data of The Cancer Genome Atlas Breast Invasive Carcinoma (TCGA-BRCA) project. Luminal A, luminal B, basal-like breast cancer subtypes, and adjacent normal tissue samples were analyzed separately. Using the sets of putative targets and non-targets, we introduced adjusted isomiR targeting activity (ITA)—the number of negatively correlated potential isomiR targets adjusted by the background (estimated using non-target mRNAs). We found that for most isomiRs a significant negative correlation between isomiR-mRNA expression levels appeared more often in a set of predicted targets compared to the non-targets. This trend was detected for both classical seed region binding types (8mer, 7mer-m8, 7mer-A1, 6mer) predicted by TargetScan and the non-classical ones (G:U wobbles and up to one mismatch or unpaired nucleotide within seed sequence) predicted by RNA22. Adjusted ITA distributions were similar for target sites located in 3′-UTRs and coding mRNA sequences, while 5′-UTRs had much lower scores. Finally, we observed strong cancer subtype-specific patterns of isomiR activity, highlighting the differences between breast cancer molecular subtypes and normal tissues. Surprisingly, our target prediction- and correlation-based estimates of isomiR activities were practically non-correlated with the average isomiR expression levels neither in cancerous nor in normal samples.

The identification of expression quantitative trait loci (eQTL) is an important component in efforts to understand how genetic variants influence disease risk. MicroRNAs (miRNAs) are short noncoding RNA molecules capable of regulating the expression of several genes simultaneously. Recently, several novel isomers of miRNAs (isomiRs) that differ slightly in length and sequence composition compared to their canonical miRNAs have been reported. Here we present isomiR-eQTL, a user-friendly database designed to help researchers find single nucleotide polymorphisms (SNPs) that can impact miRNA (miR-eQTL) and isomiR expression (isomiR-eQTL) in 30 cancer types. The isomiR-eQTL includes a total of 152,671 miR-eQTLs and 2,390,805 isomiR-eQTLs at a false discovery rate (FDR) of 0.05. It also includes 65,733 miR-eQTLs overlapping known cancer-associated loci identified through genome-wide association studies (GWAS). To the best of our knowledge, this is the first study investigating the impact of SNPs on isomiR expression at the genome-wide level. This database may pave the way for researchers toward finding a model for personalised medicine in which miRNAs, isomiRs, and genotypes are utilised. [ABSTRACT FROM AUTHOR]

Inaccurate cleavage of pri- and pre-miRNA hairpins by Drosha and Dicer results in the generation of miRNA isoforms known as isomiRs. isomiRs with 50-end variations (50-isomiRs) create a new dimension in miRNA research since they have different seed regions and distinct targetomes. We developed isomiRTar (https://isomirtar.hse.ru)--a comprehensive portal that allows one to analyze expression profiles and targeting activity of 50-isomiRs in cancer. Using the Cancer Genome Atlas sequencing data, we compiled the list of 1022 50-isomiRs expressed in 9282 tumor samples across 31 cancer types. Sequences of these isomiRs were used to predict target genes with miRDB and TargetScan. The putative interactions were then subjected to the co-expression analysis in each cancer type to identify isomiR-target pairs supported by significant negative correlations. Downstream analysis of the data deposited in isomiRTar revealed both cancer-specific and cancer-conserved 50-isomiR expression landscapes. Pairs of isomiRs differing in one nucleotide shift from 50-end had poorly overlapping targetomes with the median Jaccard index of 0.06. The analysis of colorectal cancer 50-isomiR-mediated regulatory networks revealed promising candidate tumor suppressor isomiRs: hsa- miR-203a-3p|+1, hsa-miR-192-5p|+1 and hsa-miR-148a-3p|0. In summary, we believe that isomiRTar will help researchers find novel mechanisms of isomiR-mediated gene silencing in different types of cancer. [ABSTRACT FROM AUTHOR]

Inaccurate cleavage of pri- and pre-miRNA hairpins by Drosha and Dicer results in the generation of miRNA isoforms known as isomiRs. isomiRs with 5′-end variations (5′-isomiRs) create a new dimension in miRNA research since they have different seed regions and distinct targetomes. We developed isomiRTar (https://isomirtar.hse.ru)—a comprehensive portal that allows one to analyze expression profiles and targeting activity of 5′-isomiRs in cancer. Using the Cancer Genome Atlas sequencing data, we compiled the list of 1022 5′-isomiRs expressed in 9282 tumor samples across 31 cancer types. Sequences of these isomiRs were used to predict target genes with miRDB and TargetScan. The putative interactions were then subjected to the co-expression analysis in each cancer type to identify isomiR-target pairs supported by significant negative correlations. Downstream analysis of the data deposited in isomiRTar revealed both cancer-specific and cancer-conserved 5′-isomiR expression landscapes. Pairs of isomiRs differing in one nucleotide shift from 5′-end had poorly overlapping targetomes with the median Jaccard index of 0.06. The analysis of colorectal cancer 5′-isomiR-mediated regulatory networks revealed promising candidate tumor suppressor isomiRs: hsa-miR-203a-3p—+1, hsa-miR-192-5p—+1 and hsa-miR-148a-3p—0. In summary, we believe that isomiRTar will help researchers find novel mechanisms of isomiR-mediated gene silencing in different types of cancer.

Circular RNAs (circRNAs), a class of new endogenous non-coding RNAs (ncRNAs), are closely related to the carcinogenic process and play a critical role in tumor metastasis. CircRNAs can lay the foundation for tumor metastasis via promoting tumor angiogenesis, make tumor cells gain the ability of migration and invasion by regulating epithelial-mesenchymal transition (EMT), interact with immune cells, cytokines, chemokines, and other non-cellular components in the tumor microenvironment, damage the normal immune function or escape the immunosuppressive network, and further promote cell survival and metastasis. Herein, based on the characteristics and biological functions of circRNA, we elaborated on the effect of circRNA via circRNA-associated competing endogenous RNA (ceRNA) network by acting as miRNA/isomiR sponges on tumor angiogenesis, cancer cell migration and invasion, and interaction with the tumor microenvironment (TME), then explored the potential interactions across different RNAs, and finally discussed the potential clinical value and application as a promising biomarker. These results provide a theoretical basis for the further application of metastasis-related circRNAs in cancer treatment. In summary, we briefly summarize the diverse roles of a circRNA-associated ceRNA network in cancer metastasis and the potential clinical application, especially the interaction of circRNA and miRNA/isomiR, which may complicate the RNA regulatory network and which will contribute to a novel insight into circRNA in the future. [ABSTRACT FROM AUTHOR]

ABSTRACT: In order to enrich the knowledge of chicken transcriptomic response to Salmonella enterica serovar Enteritidis infection, 2-day-old chicks were orally inoculated with this bacteria (1.0 × 108 cfu/mL), and then the cecum tissues of 3 days post-inoculation were utilized for RNA sequencing (6 replicates each for treatment group and control group). After analysis, we found a variety of inflammatory genes were triggered at the mRNA level upon infection. Notably, the expression profiles at the miRNA level and the isomiR level were heterogeneous. Certain isomiRs of chicken miR-146b-5p were significantly increased by more than 2 times compared to control (Padj < 0.05). Combining the bioinformatics prediction, transcriptome data and RT-qPCR results, we deduced that the isomiRs of chicken miR-146b-5p might act to sustain the RIG-I-like receptor signaling and type I interferon induction by repressing USP3 transcript. Our findings provide a new perspective on the regulatory function of miR-146b-5p and facilitate the study of isomiRs.

Abstract Background Next generation sequencing has allowed the discovery of miRNA isoforms, termed isomiRs. Some isomiRs are derived from imprecise processing of pre-miRNA precursors, leading to length variants. Additional variability is introduced by non-templated addition of bases at the ends or editing of internal bases, resulting in base differences relative to the template DNA sequence. We hypothesized that some component of the isomiR variation reported so far could be due to systematic technical noise and not real. Results We have developed the XICRA pipeline to analyze small RNA sequencing data at the isomiR level. We exploited its ability to use single or merged reads to compare isomiR results derived from paired-end (PE) reads with those from single reads (SR) to address whether detectable sequence differences relative to canonical miRNAs found in isomiRs are true biological variations or the result of errors in sequencing. We have detected non-negligible systematic differences between SR and PE data which primarily affect putative internally edited isomiRs, and at a much smaller frequency terminal length changing isomiRs. This is relevant for the identification of true isomiRs in small RNA sequencing datasets. Conclusions We conclude that potential artifacts derived from sequencing errors and/or data processing could result in an overestimation of abundance and diversity of miRNA isoforms. Efforts in annotating the isomiRnome should take this into account.

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

Abstract MicroRNA-derived structures play impressive roles in various biological processes. So dysregulation of miRNAs can lead to different human diseases. Recent studies have extended our comprehension of the control of miRNA function and features. Here, we overview some remarkable miRNA properties that have potential implications for the miRNA functions, including different variants of a miRNA called isomiRs, miRNA arm selection/arm switching, and the effect of these factors on miRNA target selection. Besides, we review some aspects of miRNA interactions such as the interaction between epigenetics and miRNA (different miRNAs and their related processing enzymes are epigenetically regulated by multiple DNA methylation enzymes. moreover, DNA methylation could be controlled by diverse mechanisms related to miRNAs), direct and indirect crosstalk between miRNA and lnc (Long Non-Coding) RNAs as a further approach to conduct intercellular regulation called “competing endogenous RNA” (ceRNA) that is involved in the pathogenesis of different diseases, and the interaction of miRNA activities and some Xeno-infectious (virus/bacteria/parasite) factors, which result in modulation of the pathogenesis of infections. This review provides some related studies to a better understanding of miRNA involvement mechanisms and overcoming the complexity of related diseases that may be applicable and useful to prognostic, diagnostic, therapeutic purposes and personalized medicine in the future.

Cholangiocarcinomas (CCAs) are tumors that arise from the cholangiocytes. Although some genes have been shown with important roles in pathological process, interactions or cross-talks among different RNAs are important to understand the detailed molecular mechanisms in cancer development, especially discussing cross-talks among isomiRs and other RNAs. Herein, to characterize crucial genes in CCA, the protein expression profile was performed to survey potential crucial mRNAs and related non-coding RNAs (ncRNAs) in mRNA-ncRNA network, mainly including miRNAs/isomiRs and lncRNAs. Deregulated mRNAs were firstly obtained if consistent expression patterns were found at protein and mRNA levels, and related miRNAs/isomiRs were screened according to regulatory relationships. Diverse isomiRs from a given miRNA locus also contributed to interactions between the small RNAs and target mRNAs, and miRNAs were further used to survey related lncRNAs to expand the interactions. Thus, several groups of RNAs were constructed as candidate competitive endogenous RNA (ceRNA) networks. Finally, we found that RAB11FIP1:miR-101-3p:MIR3142HG may be a potential ceRNA network, and the interactions among them may be more complex due to variety of isomiRs. Simultaneously, RAB11FIP1 and miR-194-5p were also detected other related lncRNAs (FBXL19-AS1, SNHG1 and PVT1) that may be crucial in coding-non-coding RNA regulatory network. Our results show that diverse isomiRs with sequence and expression heterogeneities contribute to ceRNA regulatory network that may have crucial roles in CCA, which will expand our understanding of interactions among diverse RNAs and their contributions in cancer development.

Abstract Background Mature sperm carry thousands of RNAs, including mRNAs, lncRNAs, tRNAs, rRNAs and sncRNAs, though their functional significance is still a matter of debate. Growing evidence suggests that sperm RNAs, especially sncRNAs, are selectively retained during spermiogenesis or specifically transferred during epididymis maturation, and are thus delivered to the oocyte at fertilization, providing resources for embryo development. However , a deep characterization of the sncRNA content of bull sperm and its expression profile across breeds is currently lacking. To fill this gap, we optimized a guanidinium–Trizol total RNA extraction protocol to prepare high-quality RNA from frozen bull sperm collected from 40 representative bulls from six breeds. Deep sequencing was performed (40 M single 50-bp reads per sample) to establish a comprehensive repertoire of cattle sperm sncRNA. Results Our study showed that it comprises mostly piRNAs (26%), rRNA fragments (25%), miRNAs (20%) and tRNA fragments (tsRNA, 14%). We identified 5p-halves as the predominant tsRNA subgroup in bull sperm, originating mostly from Gly and Glu isoacceptors. Our study also increased by ~ 50% the sperm repertoire of known miRNAs and identified 2022 predicted miRNAs. About 20% of sperm miRNAs were located within genomic clusters, expanding the list of known polycistronic pri-miRNA clusters and defining several networks of co-expressed miRNAs. Strikingly, our study highlighted the great diversity of isomiRs, resulting mainly from deletions and non-templated additions (A and U) at the 3p end. Substitutions within miRNA sequence accounted for 40% of isomiRs, with G>A, U>C and C>U substitutions being the most frequent variations. In addition, many sncRNAs were found to be differentially expressed across breeds. Conclusions Our study provides a comprehensive overview of cattle sperm sncRNA, and these findings will pave the way for future work on the role of sncRNAs in embryo development and their relevance as biomarkers of semen fertility.

Classical microRNA (miRNA) has been so far believed as a single sequence, but it indeed contains multiple miRNA isoforms (isomiR) with various sequences and expression patterns. It is not clear whether these diverse isomiRs have potential relationships and whether they contribute to miRNA:mRNA interactions. Here, we aimed to reveal the potential evolutionary and functional relationships of multiple isomiRs based on let-7 and miR-10 gene families that are prone to clustering together on chromosomes. Multiple isomiRs within gene families showed similar functions to their canonical miRNAs, indicating selection of the predominant sequence. IsomiRs containing novel seed regions showed increased/decreased biological function depending on whether they had more/less specific target mRNAs than their annotated seed. Few gene ontology(GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were shared among the target genes of the annotated seeds and the novel seeds. Various let-7 isomiRs with novel seed regions may cause opposing drug responses despite the fact that they are generated from the same miRNA locus and have highly similar sequences. IsomiRs, especially the dominant isomiRs with shifted seeds, may disturb the coding-non-coding RNA regulatory network. These findings provide insight into the multiple isomiRs and isomiR-mediated control of gene expression in the pathogenesis of cancer.

MicroRNA-seq data is produced by aligning small RNA sequencing reads of different miRNA transcript isoforms, called isomiRs, to known microRNAs. Aggregation to microRNA-level counts discards information and violates core assumptions of differential expression (DE) methods developed for mRNA-seq data. We establish miRglmm, a DE method for microRNA-seq data, that uses a generalized linear mixed model of isomiR-level counts, facilitating detection of miRNA with differential expression or differential isomiR usage. We demonstrate that miRglmm outperforms current DE methods in estimating DE for miRNA, whether or not there is significant isomiR variability, and simultaneously provides estimates of isomiR-level DE., Competing Interests: Competing interests The authors declare that they have no competing interests.

MiRNA isoforms (isomiRs) are single stranded small RNAs originating from the same pri-miRNA hairpin as a result of cleavage by Drosha and Dicer enzymes. Variations at the 5ʹ-end of a miRNA alter the seed region of the molecule, thus affecting the targetome of the miRNA. In this manuscript, we analysed the distribution of miRNA cleavage positions across 31 different cancers using miRNA sequencing data of TCGA project. As a result, we found that the processing positions are not tissue specific and that all miRNAs could be correctly classified as ones exhibiting homogeneous or heterogeneous cleavage at one of the four cleavage sites. In 42% of cases (42 out of 100 miRNAs), we observed imprecise 5ʹ-end Dicer cleavage, while this fraction was only 14% for Drosha (14 out of 99). To the contrary, almost all cleavage sites of 3ʹ-ends (either Drosha or Dicer) were heterogeneous. With the use of only four nucleotides surrounding a 5ʹ-end Dicer cleavage position we built a model which allowed us to distinguish between homogeneous and heterogeneous cleavage with the reliable quality (ROC AUC = 0.68). Finally, we showed the possible applications of the study by the analysis of two 5ʹ-end isoforms originating from the same exogeneous shRNA hairpin. It turned out that the less expressed shRNA variant was functionally active, which led to the increased off-targeting. Thus, the obtained results could be applied to the design of shRNAs whose processing will result in a single 5ʹ-variant. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) are approximately 20–22 nucleotides in length, which are well known to participate in the post-transcriptional modification. The mature miRNAs were observed to be varied on 5′ or 3′ that raise another term—the isoforms of mature miRNAs (isomiRs), which have been proven not the artifacts and discussed widely recently. In our research, we focused on studying the 5′ isomiRs in lung adenocarcinoma (LUAD) in The Cancer Genome Atlas (TCGA). We characterized 75 isomiRs significantly associated with better prognosis and 43 isomiRs with poor prognosis. The 75 protective isomiRs can successfully distinguish tumors from normal samples and are expressed differently between patients of early and late stages. We also found that most of the protective isomiRs tend to be with downstream shift and upregulated compared with those with upstream shift, implying that a possible selection occurs during cancer development. Among these protective isomiRs, we observed a highly positive and significant correlation, as well as in harmful isomiRs, suggesting cooperation within the group. However, between protective and harmful, there is no such a concordance but conversely more negative correlation, suggesting the possible antagonistic effect between protective and harmful isomiRs. We also identified that two isomiRs miR-181a-3p|-3 and miR-181a-3p|2, respectively, belong to the harmful and protective groups, suggesting a bidirectional regulation of their originated archetype—miR-181a-3p. Additionally, we found that the protective isomiRs of miR-21-5p, which is an oncomiR, may be evolved as the tumor suppressors through producing isomiRs to hinder metastasis. In summary, these results displayed the characteristics of the protective isomiRs and their potential for developing the treatment of lung cancer. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) are approximately 20–22 nucleotides in length, which are well known to participate in the post-transcriptional modification. The mature miRNAs were observed to be varied on 5′ or 3′ that raise another term—the isoforms of mature miRNAs (isomiRs), which have been proven not the artifacts and discussed widely recently. In our research, we focused on studying the 5′ isomiRs in lung adenocarcinoma (LUAD) in The Cancer Genome Atlas (TCGA). We characterized 75 isomiRs significantly associated with better prognosis and 43 isomiRs with poor prognosis. The 75 protective isomiRs can successfully distinguish tumors from normal samples and are expressed differently between patients of early and late stages. We also found that most of the protective isomiRs tend to be with downstream shift and upregulated compared with those with upstream shift, implying that a possible selection occurs during cancer development. Among these protective isomiRs, we observed a highly positive and significant correlation, as well as in harmful isomiRs, suggesting cooperation within the group. However, between protective and harmful, there is no such a concordance but conversely more negative correlation, suggesting the possible antagonistic effect between protective and harmful isomiRs. We also identified that two isomiRs miR-181a-3p|-3 and miR-181a-3p|2, respectively, belong to the harmful and protective groups, suggesting a bidirectional regulation of their originated archetype—miR-181a-3p. Additionally, we found that the protective isomiRs of miR-21-5p, which is an oncomiR, may be evolved as the tumor suppressors through producing isomiRs to hinder metastasis. In summary, these results displayed the characteristics of the protective isomiRs and their potential for developing the treatment of lung cancer.

The identification of expression quantitative trait loci (eQTL) is an important component in efforts to understand how genetic variants influence disease risk. MicroRNAs (miRNAs) are short noncoding RNA molecules capable of regulating the expression of several genes simultaneously. Recently, several novel isomers of miRNAs (isomiRs) that differ slightly in length and sequence composition compared to their canonical miRNAs have been reported. Here we present isomiR-eQTL, a user-friendly database designed to help researchers find single nucleotide polymorphisms (SNPs) that can impact miRNA (miR-eQTL) and isomiR expression (isomiR-eQTL) in 30 cancer types. The isomiR-eQTL includes a total of 152,671 miR-eQTLs and 2,390,805 isomiR-eQTLs at a false discovery rate (FDR) of 0.05. It also includes 65,733 miR-eQTLs overlapping known cancer-associated loci identified through genome-wide association studies (GWAS). To the best of our knowledge, this is the first study investigating the impact of SNPs on isomiR expression at the genome-wide level. This database may pave the way for researchers toward finding a model for personalised medicine in which miRNAs, isomiRs, and genotypes are utilised.

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. [ABSTRACT FROM AUTHOR]

Abstract Background MicroRNAs (miRNA) and other small RNAs are frequently dysregulated in cancer and are promising biomarkers for colon cancer. Here we profile human, virus and bacteria small RNAs in normal and tumor tissue from early stage colon cancer and correlate the expression with clinical parameters. Methods Small RNAs from colon cancer tissue and adjacent normal mucosa of 48 patients were sequenced using Illumina high-throughput sequencing. Clinical parameters were correlated with the small RNA expression data using linear models. We performed a meta-analysis by comparing publicly available small RNA sequencing datasets with our original sequencing data to confirm the main findings. Results We identified 331 differentially expressed miRNAs between tumor and normal samples. We found that the major changes in miRNA expression between left and right colon are due to miRNAs located within the Hox-developmental genes, including miR-10b, miR-196b and miR-615. Further, we identified new miRNAs associated with microsatellite instability (MSI), including miR-335, miR-26 and miR-625. We performed a meta-analysis on all publicly available miRNA-seq datasets and identified 117 common miRNAs that were differentially expressed between tumor and normal tissue. The miRNAs miR-135b and miR-31 were the most significant upregulated miRNA in tumor across all datasets. The miRNA miR-133a was the most strongly downregulated miRNA in our dataset and also showed consistent downregulation in the other datasets. The miRNAs associated with MSI and tumor location in our data showed similar changes in the other datasets. Finally, we show that small RNAs from Epstein-Barr virus and Fusobacterium nucleatum are differentially expressed between tumor and normal adjacent tissue. Conclusions Small RNA profiling in colon cancer tissue revealed novel RNAs associated with MSI and tumor location. We show that Fusobacterium nucleatum are detectable at the RNA-level in colon tissue, and that both Fusobacterium nucleatum and Epstein-Barr virus separate tumor and normal tissue.

Abstract Background The miRNA isoforms (isomiRs) have been suggested to regulate the same pathways as the canonical miRNA and play an important biological role in miRNA-mediated gene regulation. Recently, a study has demonstrated that the presence or absence of all isomiRs could efficiently discriminate amongst 32 TCGA cancer types. Besides, an effective reduction of distinguishing isomiR features for multiclass tumor discrimination must have a major impact on our understanding of the disease and treatment of cancer. Methods In this study, we have constructed a combination of the genetic algorithms (GA) with Random Forest (RF) algorithms to detect reliable sets of cancer-associated 5’isomiRs from TCGA isomiR expression data for multiclass tumor classification. Results We obtained 100 sets of the optimal predictive features, each of which comprised of 50–5’isomiRs that could effectively classify with an average sensitivity of 92% samples from 32 different tumor types. We calculated the frequency with which a 5’isomiR found in these sets as measuring its importance for tumor classification. Many highly frequent 5’isomiRs with different 5′ loci from canonical miRNAs were detected in these sets, supporting that the isomiRs play a significant role in the multiclass tumor classification. The further functional enrichment analysis showed that the target genes of the 10 most frequently appearing 5’isomiRs were involved in the activity of transcription activator and protein kinase and cell-cell adhesion. Conclusions The findings of the present study indicated that the 5’isomiRs might be employed for multiclass tumor classification and the suggested that GA/RF model could perform effective tumor classification by a series of largely independent optimal predictor 5′ isomiR sets.

Circular RNAs (circRNAs), a class of new endogenous non-coding RNAs (ncRNAs), are closely related to the carcinogenic process and play a critical role in tumor metastasis. CircRNAs can lay the foundation for tumor metastasis via promoting tumor angiogenesis, make tumor cells gain the ability of migration and invasion by regulating epithelial-mesenchymal transition (EMT), interact with immune cells, cytokines, chemokines, and other non-cellular components in the tumor microenvironment, damage the normal immune function or escape the immunosuppressive network, and further promote cell survival and metastasis. Herein, based on the characteristics and biological functions of circRNA, we elaborated on the effect of circRNA via circRNA-associated competing endogenous RNA (ceRNA) network by acting as miRNA/isomiR sponges on tumor angiogenesis, cancer cell migration and invasion, and interaction with the tumor microenvironment (TME), then explored the potential interactions across different RNAs, and finally discussed the potential clinical value and application as a promising biomarker. These results provide a theoretical basis for the further application of metastasis-related circRNAs in cancer treatment. In summary, we briefly summarize the diverse roles of a circRNA-associated ceRNA network in cancer metastasis and the potential clinical application, especially the interaction of circRNA and miRNA/isomiR, which may complicate the RNA regulatory network and which will contribute to a novel insight into circRNA in the future.

Abstract Background Gene expression-based profiling has been used to identify biomarkers for different breast cancer subtypes. However, this technique has many limitations. IsomiRs are isoforms of miRNAs that have critical roles in many biological processes and have been successfully used to distinguish various cancer types. Biomarker isomiRs for identifying different breast cancer subtypes has not been investigated. For the first time, we aim to show that isomiRs are better performing biomarkers and use them to explain molecular differences between breast cancer subtypes. Results In this study, a novel method is proposed to identify specific isomiRs that faithfully classify breast cancer subtypes. First, as a null hypothesis method we removed the lowly expressed isomiRs from small sequencing data generated from diverse breast cancers types. Second, we developed an improved mutual information-based feature selection method to calculate the weight of each isomiR expression. The weight of isomiR measures the importance of a given isomiR in classifying breast cancer subtypes. The improved mutual information enables to apply the dataset in which the feature is continuous data and label is discrete data; whereby, the traditional mutual information cannot be applied in this dataset. Finally, the support vector machine (SVM) classifier is applied to find isomiR biomarkers for subtyping. Conclusions Here we demonstrate that isomiRs can be used as biomarkers in the identification of different breast cancer subtypes, and in addition, they may provide new insights into the diverse molecular mechanisms of breast cancers. We have also shown that the classification of different subtypes of breast cancer based on isomiRs expression is more effective than using published gene expression profiling. The proposed method provides a better performance outcome than Fisher method and Hellinger method for discovering biomarkers to distinguish different breast cancer subtypes. This novel technique could be directly applied to identify biomarkers in other diseases.

Background: Next generation sequencing has allowed the discovery of miRNA isoforms, termed isomiRs. Some isomiRs are derived from imprecise processing of pre-miRNA precursors, leading to length variants. Additional variability is introduced by non-templated addition of bases at the ends or editing of internal bases, resulting in base differences relative to the template DNA sequence. We hypothesized that some component of the isomiR variation reported so far could be due to systematic technical noise and not real. Results: We have developed the XICRA pipeline to analyze small RNA sequencing data at the isomiR level. We exploited its ability to use single or merged reads to compare isomiR results derived from paired-end (PE) reads with those from single reads (SR) to address whether detectable sequence differences relative to canonical miRNAs found in isomiRs are true biological variations or the result of errors in sequencing. We have detected non-negligible systematic differences between SR and PE data which primarily affect putative internally edited isomiRs, and at a much smaller frequency terminal length changing isomiRs. This is relevant for the identification of true isomiRs in small RNA sequencing datasets. Conclusions: We conclude that potential artifacts derived from sequencing errors and/or data processing could result in an overestimation of abundance and diversity of miRNA isoforms. Efforts in annotating the isomiRnome should take this into account. [ABSTRACT FROM AUTHOR]

An increasing number of microRNAs (miRNAs) and miRNA‐related sequences produced during miRNA biogenesis, comprising the isomiRome, have been recently highlighted in different species as critical mediators of environmental stress responses. Conifers have some of the largest known genomes but an extensive characterization of the isomiRome from any conifer species has been lacking. We provide here a comprehensive overview of the Pinus pinaster isomiRome expressed in roots, stem and needles under well‐watered and drought conditions. From the 13,441 unique small RNA sequences identified, 2,980 were annotated as canonical miRNAs or miRNA* and the remaining were classified as isomiRNA or miRNA‐like sequences. A survey of their expression patterns highlighted roots as the most responsive organ under drought, where specific sequences of which a 24‐nt novel miRNA stood out, were strongly down‐regulated. Given the putative roles of the miRNA‐targeted transcripts validated specifically in root tissues, some of the miRNAs, conserved and novel, are shortlisted as potential regulators of drought response. These results provide a valuable resource for comparative studies between gymnosperms and angiosperms. Furthermore, it evidences high transferability of the isomiRome between pine species being a useful basis for further molecular regulation and physiological studies, and especially those focused on adaptation to drought conditions. The present work provides a first characterization of the rich and variable isomiRome encoded by a conifer genome. Results highlight a differential distribution of the isomiRome in different vegetative organs and under contrasting water availability, and potentially targeting transcripts controlling key physiological responses to drought. [ABSTRACT FROM AUTHOR]

MicroRNA-derived structures play impressive roles in various biological processes. So dysregulation of miRNAs can lead to different human diseases. Recent studies have extended our comprehension of the control of miRNA function and features. Here, we overview some remarkable miRNA properties that have potential implications for the miRNA functions, including different variants of a miRNA called isomiRs, miRNA arm selection/arm switching, and the effect of these factors on miRNA target selection. Besides, we review some aspects of miRNA interactions such as the interaction between epigenetics and miRNA (different miRNAs and their related processing enzymes are epigenetically regulated by multiple DNA methylation enzymes. moreover, DNA methylation could be controlled by diverse mechanisms related to miRNAs), direct and indirect crosstalk between miRNA and lnc (Long Non-Coding) RNAs as a further approach to conduct intercellular regulation called "competing endogenous RNA" (ceRNA) that is involved in the pathogenesis of different diseases, and the interaction of miRNA activities and some Xeno-infectious (virus/bacteria/parasite) factors, which result in modulation of the pathogenesis of infections. This review provides some related studies to a better understanding of miRNA involvement mechanisms and overcoming the complexity of related diseases that may be applicable and useful to prognostic, diagnostic, therapeutic purposes and personalized medicine in the future. [ABSTRACT FROM AUTHOR]

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/. [ABSTRACT FROM AUTHOR]

Cystic fibrosis (CF), a genetic disorder, is characterized by chronic lung disease. Small non‐coding RNAs are key regulators of gene expression and participate in various processes, which are dysregulated in CF; however, they remain poorly studied. Here, we determined the complete microRNAs (miRNAs) expression pattern in three CF ex vivo models. The miRNA profiles of air–liquid interface cultures of airway epithelia (bronchi, nasal cells, and nasal polyps) samples from patients with CF and non‐CF controls were obtained by deep sequencing. Compared with non‐CF controls, several miRNAs were deregulated in CF samples; for instance, miR‐181a‐5p and the miR‐449 family were upregulated. Moreover, mature miRNAs often showed variations (i.e. isomiRs) relative to their reference sequence, such as miR‐101, suggesting that miRNAs consist of heterogeneous repertoires of multiple isoforms with different effects on gene expression. Analysis of miR‐181a‐5p and miR‐101‐3p roles indicated that they regulate the expression of WISP1, a key component of cell proliferation/migration programs. We showed that miR‐101 and miR‐181a‐5p participated in aberrant recapitulation of wound healing programs by controlling WISP1 mRNA and protein level. Our miRNA expression data bring new insights into CF physiopathology and define new potential therapeutic targets in CF. © 2020 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

Regulatory changes occurring early in colorectal cancer development remain poorly investigated. Since the majority of cases develop from polyps in the adenoma-carcinoma transition, a search of early molecular features, such as aberrations in miRNA expression occurring prior to cancer development, would enable identification of potentially causal, rather than consequential, candidates in the progression of polyp to cancer. In the current study, by employing small RNA-seq profiling of colon biopsy samples, we described differentially expressed miRNAs and their isoforms in the adenoma-carcinoma transition. Analysis of healthy-adenoma-carcinoma sequence in an independent validation group enabled us to identify early deregulated miRNAs including hsa-miR-1246 and hsa-miR-215-5p, the expressions of which are, respectively, gradually increasing and decreasing. Loss-of-function experiments revealed that inhibition of hsa-miR-1246 lead to reduced cell viability, colony formation, and migration rate, thereby indicating an oncogenic effect of this miRNA in vitro. Subsequent western blot and luciferase reporter assay provided evidence of hsa-miR-1246 being involved in the regulation of target AXIN2 and CFTR genes’ expression. To conclude, the present study revealed possible involvement of hsa-miR-1246 in early colorectal cancer development and regulation of tumor suppressors AXIN2 and CFTR.

Neuroinflammation is a secondary injury mechanism that evolves in the brain for months after traumatic brain injury (TBI). We hypothesized that an altered small non-coding RNA (sncRNA) signature plays a key role in modulating post-TBI secondary injury and neuroinflammation. At 3threemonths post-TBI, messenger RNA sequencing (seq) and small RNAseq were performed on samples from the ipsilateral thalamus and perilesional cortex of selected rats with a chronic inflammatory endophenotype, and sham-operated controls. The small RNAseq identified dysregulation of 2 and 19 miRNAs in the thalamus and cortex, respectively. The two candidates from the thalamus and the top ten from the cortex were selected for validation. In the thalamus, miR-146a-5p and miR-155-5p levels were upregulated, and in the cortex, miR-375-3p and miR-211-5p levels were upregulated. Analysis of isomiRs of differentially expressed miRNAs identified 3′ nucleotide additions that were increased after TBI. Surprisingly, we found fragments originating from 16 and 13 tRNAs in the thalamus and cortex, respectively. We further analyzed two upregulated fragments, 3′tRF-IleAAT and 3′tRF-LysTTT. Increased expression of the full miR-146a profile, and 3′tRF-IleAAT and 3′tRF-LysTTT was associated with a worse behavioral outcome in animals with chronic neuroinflammation. Our results highlight the importance of understanding the regulatory roles of as-yet unknown sncRNAs for developing better strategies to treat TBI and neuroinflammation.

Abstract Background miRNAs play important roles in the regulation of gene expression. The rapidly developing field of microRNA sequencing (miRNA-seq; small RNA-seq) needs comprehensive, robust, user-friendly and standardized bioinformatics tools to analyze these large datasets. We present miRge 2.0, in which multiple enhancements were made towards these goals. Results miRge 2.0 has become more comprehensive with increased functionality including a novel miRNA detection method, A-to-I editing analysis, integrated standardized GFF3 isomiR reporting, and improved alignment to miRNAs. The novel miRNA detection method uniquely uses both miRNA hairpin sequence structure and composition of isomiRs resulting in higher specificity for potential miRNA identification. Using known miRNA data, our support vector machine (SVM) model predicted miRNAs with an average Matthews correlation coefficient (MCC) of 0.939 over 32 human cell datasets and outperformed miRDeep2 and miRAnalyzer regarding phylogenetic conservation. The A-to-I editing detection strongly correlated with a reference dataset with adjusted R2 = 0.96. miRge 2.0 is the most up-to-date aligner with custom libraries to both miRBase v22 and MirGeneDB v2.0 for 6 species: human, mouse, rat, fruit fly, nematode and zebrafish; and has a tool to create custom libraries. For user-friendliness, miRge 2.0 is incorporated into bcbio-nextgen and implementable through Bioconda. Conclusions miRge 2.0 is a redesigned, leading miRNA RNA-seq aligner with several improvements and novel utilities. miRge 2.0 is freely available at: https://github.com/mhalushka/miRge.

MicroRNAs (miRNAs) are small endogenous non-coding RNAs in eukaryotes which regulate the expression of numerous genes post-transcriptionally, thereby playing critical roles in cells and organismal development. The high-throughput sequencing technologies enable the effective detection and annotation of miRNAs. Several miRNA variants with heterogeneous ends, lengths, and sequences can be generated from a single miRNA locus. Discovery of these miRNA variants, also known as miRNA isoforms or isomiRs, has made our understanding of the cells' miRNome deeper than previously pictured. Despite their wide presence in multiple datasets, the different possible origins and true biological significance of isomiRs are yet to be uncovered. Several recent emerging studies suggest that isomiRs are biologically active and non-randomly formed. This review aims to provide a comprehensive insight into the origins and biological importance of isomiRs, highlighting the enormous complexity of miRNA regulatory networks which broadens our knowledge about the post-transcriptional gene regulation in plants. • miRNAs are known to regulate a wide range of key biological processes in plants • Several miRNA variants can be generated from a single miRNA locus • The origins and true biological significance of isomiRs are yet to be uncovered • studies suggest that isomiRs are biologically active and non-randomly formed • It broadens our knowledge about the post-transcriptional gene regulation in plants [ABSTRACT FROM AUTHOR]

Background: Mature sperm carry thousands of RNAs, including mRNAs, lncRNAs, tRNAs, rRNAs and sncRNAs, though their functional significance is still a matter of debate. Growing evidence suggests that sperm RNAs, especially sncRNAs, are selectively retained during spermiogenesis or specifically transferred during epididymis maturation, and are thus delivered to the oocyte at fertilization, providing resources for embryo development. However , a deep characterization of the sncRNA content of bull sperm and its expression profile across breeds is currently lacking. To fill this gap, we optimized a guanidinium–Trizol total RNA extraction protocol to prepare high-quality RNA from frozen bull sperm collected from 40 representative bulls from six breeds. Deep sequencing was performed (40 M single 50-bp reads per sample) to establish a comprehensive repertoire of cattle sperm sncRNA. Results: Our study showed that it comprises mostly piRNAs (26%), rRNA fragments (25%), miRNAs (20%) and tRNA fragments (tsRNA, 14%). We identified 5p-halves as the predominant tsRNA subgroup in bull sperm, originating mostly from Gly and Glu isoacceptors. Our study also increased by ~ 50% the sperm repertoire of known miRNAs and identified 2022 predicted miRNAs. About 20% of sperm miRNAs were located within genomic clusters, expanding the list of known polycistronic pri-miRNA clusters and defining several networks of co-expressed miRNAs. Strikingly, our study highlighted the great diversity of isomiRs, resulting mainly from deletions and non-templated additions (A and U) at the 3p end. Substitutions within miRNA sequence accounted for 40% of isomiRs, with G>A, U>C and C>U substitutions being the most frequent variations. In addition, many sncRNAs were found to be differentially expressed across breeds. Conclusions: Our study provides a comprehensive overview of cattle sperm sncRNA, and these findings will pave the way for future work on the role of sncRNAs in embryo development and their relevance as biomarkers of semen fertility. [ABSTRACT FROM AUTHOR]

Uveal melanoma (UVM) is the most common primary intraocular malignancy in adults. With over 50% of patients developing metastatic disease, there is an unmet need for improved diagnostic and therapeutic options. Efforts to understand the molecular biology of the disease have revealed several markers that correlate with patient prognosis, including the copy number of chromosome 3, genetic alterations in the BAP1, EIF1AX and SF3B1 genes, and other transcriptional features. Here, we expand upon previous reports by comprehensively characterizing the short RNA‐ome in 80 primary UVM tumor samples. In particular, we describe a previously unseen complex network involving numerous regulatory molecules that comprise microRNA (miRNAs), novel UVM‐specific miRNA loci, miRNA isoforms (isomiRs), and tRNA‐derived fragments (tRFs). Importantly, we show that the abundance profiles of isomiRs and tRFs associate with various molecular phenotypes, metastatic disease, and patient survival. Our findings suggest deep involvement of isomiRs and tRFs in the disease etiology of UVM. We posit that further study and characterization of these novel molecules will improve understanding of the mechanisms underlying UVM, and lead to the development of new diagnostic and therapeutic approaches. [ABSTRACT FROM AUTHOR]

After over a decade of nanosafety research, it is indisputable that the vast majority of nano-sized particles induce a plethora of adverse cellular responses – the severity of which is linked to the material's physicochemical properties. Differentiated THP-1 cells were previously exposed for 6 h and 24 h to silver, titanium dioxide, and zinc oxide nanoparticles at the maximum molar concentration at which no more than 15% cellular cytotoxicity was observed. All three nanoparticles differed in extent of induction of biological pathways corresponding to immune response signaling and metal ion homeostasis. In this study, we integrated gene and miRNA expression profiles from the same cells to propose miRNA biomarkers of adverse exposure to metal-based nanoparticles. We employed RNA sequencing together with a quantitative strategy that also enables analysis of the overlooked repertoire of length and sequence miRNA variants called isomiRs. Whilst only modest changes in expression were observed within the first 6 h of exposure, the miRNA/isomiR (miR) profiles of each nanoparticle were unique. Via canonical correlation and pathway enrichment analyses, we identified a co-regulated miR-mRNA cluster, predicted to be highly relevant for cellular response to metal ion homeostasis. These miRs were annotated to be canonical or variant isoforms of hsa-miR-142-5p, -342-3p, -5100, -6087, -6894-3p, and -7704. Hsa-miR-5100 was differentially expressed in response to each nanoparticle in both the 6 h and 24 h exposures. Taken together, this co-regulated miR-mRNA cluster could represent potential biomarkers of sub-toxic metal-based nanoparticle exposure. [ABSTRACT FROM AUTHOR]

This study aimed to discuss the potential roles of isomiRs of miR-27 family in metabolisms associated with disease via analyses of their evolution, expression, and function. miR-27b-3p was relatively highly expressed in liver cancer samples compared to miR-27a-3p and miR-27-5p loci. The diversity of isomiRs in miR-27-3p locus is similar to that of miRNAs among homologous genes. IsomiRs exhibited variable expression across different cancer tissue types, and some of them were abnormally expressed in ob/ob mice. Further experimental validation indicated that the protein expression of metabolism-related proteins, including PEPCK, G6Pase, FAS, and CPT1A, were significantly suppressed when canonical miR-27b was transfected into AML-12 cells. In contrast, the expression of these proteins was only slightly inhibited by isomiR-27b-1 or isomiR-27b-2 after transfection into AML-12 cells. These observations support that isomiRs exhibiting sequence divergence are functional regulatory molecules, and that they may contribute to biological processes via coordinated interactions in regulatory networks. [ABSTRACT FROM AUTHOR]

MicroRNAs as regulators of gene expression have been known for over a decade and have been correlated with various types of abiotic and biotic stresses. IsomiRs are modified forms of typical miRNAs altered by one or few nucleotides as a result of post transcriptional modifications of miRNAs at terminals or SNPs containing miRNA sequences producing 5′/3′ isomiRs or SNP_isomiRs. These isoforms function exactly like miRNAs, but sometimes alter the target gene preference particularly when they differ in sequences within seed region. Existence of isomiRs was considered earlier as errors in sequencing techniques due to unavailability of proper detection tools. Bread wheat is the most cultivated cereal crop globally with a strong hold on world's economy. Wheat production is frequently affected by rust diseases. Leaf rust, caused by Puccinia triticina, severely affects grain filling. Four small RNA libraries were prepared from leaves of two wheat Near Isogenic Lines, HD2329 (susceptible) and HD2329 + Lr24 (resistant) both under mock- and pathogen-inoculated conditions and sequenced using Illumina NGS. Several novel miRNAs were detected from these sequences. The present study focuses on identification of isomiRs derived from these miRNAs and their target genes. In total, 66 and 38 unique 5′ and 3′ isomiRs respectively were identified from 37 miRNAs. IsomiRs targeted genes with functions like amino acid metabolism, replication, transport, chelation of reactive oxygen species in wheat while genes with SunT and Structural Maintenance of Chromosome domain of Puccinia. This study will provide a basis for exploitation of isomiRs for genetic improvements in wheat. [ABSTRACT FROM AUTHOR]

Cutaneous melanoma (CM) is the most lethal tumor among skin cancers, and its incidence is constantly increasing. A deeper understanding of the molecular processes guiding melanoma pathogenesis could improve diagnosis, treatment and prognosis. MicroRNAs play a key role in melanoma biology. Recently, next generation sequencing (NGS) experiments, designed to assess small-RNA expression, revealed the existence of microRNA variants with different length and sequence. These microRNA isoforms are known as isomiRs and provide an additional layer to the complex non-coding RNA world. Here, we collected data from NGS experiments to provide a comprehensive characterization of miRNA and isomiR dysregulation in benign nevi (BN) and early-stage melanomas. We observed that melanoma and BN express different and specific isomiRs and have a different isomiR abundance distribution. Moreover, isomiRs from the same microRNA can have opposite expression trends between groups. Using The Cancer Genome Atlas (TCGA) dataset of skin cancers, we analyzed isomiR expression in primary melanoma and melanoma metastasis and tested their association with NF1, BRAF and NRAS mutations. IsomiRs differentially expressed were identified and catalogued with reference to the canonical form. The reported non-random dysregulation of specific isomiRs contributes to the understanding of the complex melanoma pathogenesis and serves as the basis for further functional studies.

MicroRNAs (miRNAs) are crucial regulators of gene expression in plants and a growing number of novel miRNA genes have been cloned in rice in recent years. However, there is no evidence that all miRNAs have been discovered, especially for those low expression ones which are difficult to be found by conventional methods. By taking advantage of the finding that DCL1-mediated cleavage signals for the processing of the miRNA precursors could be used as the clues for novel miRNAs' discovery, a genome-wide search for rice miRNA candidates was carried out. As a result, 51 previously validated miRNAs and 24 novel miRNA candidates were discovered. After target prediction and degradome sequencing data-based validation, coupled with reverse approach retest, 10 miRNA candidate–mRNA target pairs were further identified, providing a basis for in-depth functional analysis of these miRNA candidates. Besides, some isomiRs found in this study showed more likely to be the real miRNAs. We also found an exceptional example which did not obey the rule that 22-nt miRNAs have the ability to trigger the phased siRNAs production from the cleaved targets. [ABSTRACT FROM AUTHOR]

Background: MicroRNAs (miRNA) and other small RNAs are frequently dysregulated in cancer and are promising biomarkers for colon cancer. Here we profile human, virus and bacteria small RNAs in normal and tumor tissue from early stage colon cancer and correlate the expression with clinical parameters.Methods: Small RNAs from colon cancer tissue and adjacent normal mucosa of 48 patients were sequenced using Illumina high-throughput sequencing. Clinical parameters were correlated with the small RNA expression data using linear models. We performed a meta-analysis by comparing publicly available small RNA sequencing datasets with our original sequencing data to confirm the main findings.Results: We identified 331 differentially expressed miRNAs between tumor and normal samples. We found that the major changes in miRNA expression between left and right colon are due to miRNAs located within the Hox-developmental genes, including miR-10b, miR-196b and miR-615. Further, we identified new miRNAs associated with microsatellite instability (MSI), including miR-335, miR-26 and miR-625. We performed a meta-analysis on all publicly available miRNA-seq datasets and identified 117 common miRNAs that were differentially expressed between tumor and normal tissue. The miRNAs miR-135b and miR-31 were the most significant upregulated miRNA in tumor across all datasets. The miRNA miR-133a was the most strongly downregulated miRNA in our dataset and also showed consistent downregulation in the other datasets. The miRNAs associated with MSI and tumor location in our data showed similar changes in the other datasets. Finally, we show that small RNAs from Epstein-Barr virus and Fusobacterium nucleatum are differentially expressed between tumor and normal adjacent tissue.Conclusions: Small RNA profiling in colon cancer tissue revealed novel RNAs associated with MSI and tumor location. We show that Fusobacterium nucleatum are detectable at the RNA-level in colon tissue, and that both Fusobacterium nucleatum and Epstein-Barr virus separate tumor and normal tissue. [ABSTRACT FROM AUTHOR]

Background: The miRNA isoforms (isomiRs) have been suggested to regulate the same pathways as the canonical miRNA and play an important biological role in miRNA-mediated gene regulation. Recently, a study has demonstrated that the presence or absence of all isomiRs could efficiently discriminate amongst 32 TCGA cancer types. Besides, an effective reduction of distinguishing isomiR features for multiclass tumor discrimination must have a major impact on our understanding of the disease and treatment of cancer.Methods: In this study, we have constructed a combination of the genetic algorithms (GA) with Random Forest (RF) algorithms to detect reliable sets of cancer-associated 5'isomiRs from TCGA isomiR expression data for multiclass tumor classification.Results: We obtained 100 sets of the optimal predictive features, each of which comprised of 50-5'isomiRs that could effectively classify with an average sensitivity of 92% samples from 32 different tumor types. We calculated the frequency with which a 5'isomiR found in these sets as measuring its importance for tumor classification. Many highly frequent 5'isomiRs with different 5' loci from canonical miRNAs were detected in these sets, supporting that the isomiRs play a significant role in the multiclass tumor classification. The further functional enrichment analysis showed that the target genes of the 10 most frequently appearing 5'isomiRs were involved in the activity of transcription activator and protein kinase and cell-cell adhesion.Conclusions: The findings of the present study indicated that the 5'isomiRs might be employed for multiclass tumor classification and the suggested that GA/RF model could perform effective tumor classification by a series of largely independent optimal predictor 5' isomiR sets. [ABSTRACT FROM AUTHOR]