Concordant and discordant regulation of target genes by miR-31 and its isoforms

Circulating isomiRs May Be Superior Biomarkers Compared to Their Corresponding miRNAs: A Pilot Biomarker Study of Using isomiR-Ome to Detect Coronary Calcium-Based Cardiovascular Risk in Patients with NAFLD

Circulating miRNAs are increasingly being considered as biomarkers in various medical contexts, but the value of analyzing isomiRs (isoforms of canonical miRNA sequences) has not frequently been assessed. Here we hypothesize that an in-depth analysis of the full circulating miRNA landscape could identify specific isomiRs that are stronger biomarkers, compared to their corresponding miRNA, for identifying increased CV risk in patients with non-alcoholic fatty liver disease (NAFLD)-a clinical unmet need. Plasma miRNAs were sequenced with next-generation sequencing (NGS). Liver fat content was measured with magnetic-resonance spectrometry (MRS); CV risk was determined, beyond using traditional biomarkers, by a CT-based measurement of coronary artery calcium (CAC) score and the calculation of a CAC score-based CV-risk percentile (CAC-CV%). This pilot study included n = 13 patients, age > 45 years, with an MRS-measured liver fat content of ≥5% (wt/wt), and free of overt CVD. NGS identified 1103 miRNAs and 404,022 different isomiRs, of which 280 (25%) and 1418 (0.35%), respectively, passed an abundance threshold. Eighteen (sixteen/two) circulating miRNAs correlated positively/negatively, respectively, with CAC-CV%, nine of which also significantly discriminated between high/low CV risk through ROC-AUC analysis. IsomiR-ome analyses uncovered 67 isomiRs highly correlated (R ≥ 0.55) with CAC-CV%. Specific isomiRs of miRNAs 101-3p, 144-3p, 421, and 484 exhibited stronger associations with CAC-CV% compared to their corresponding miRNA. Additionally, while miRNAs 140-3p, 223-3p, 30e-5p, and 342-3p did not correlate with CAC-CV%, specific isomiRs with altered seed sequences exhibited a strong correlation with coronary atherosclerosis burden. Their predicted isomiRs-specific targets were uniquely enriched (compared to their canonical miRNA sequence) in CV Disease (CVD)-related pathways. Two of the isomiRs exhibited discriminative ROC-AUC, and another two showed a correlation with reverse cholesterol transport from cholesterol-loaded macrophages to ApoB-depleted plasma. In summary, we propose a pipeline for exploring circulating isomiR-ome as an approach to uncover novel and strong CVD biomarkers.

Tumor-derived miRNAs as tumor microenvironment regulators for synergistic therapeutic options

MicroRNAs (miRNAs) are a class of non-coding RNAs that perform post-transcriptional gene regulation. This review focuses on the role of tumor cell-derived miRNAs in the regulation of the tumor microenvironment (TME) via receptor cell recoding, including angiogenesis, expression of immunosuppressive molecules, formation of radiation resistance, and chemoresistance. Furthermore, we discuss the potential of these molecules as adjuvant therapies in combination with chemotherapy, radiotherapy, or immunotherapy, as well as their advantages as efficacy predictors for personalized therapy. MiRNA-based therapeutic agents for tumors are currently in clinical trials, and while challenges remain, additional research on tumor-derived miRNAs is warranted, which may provide significant clinical benefits to cancer patients.

Dynamic and static circulating cancer microRNA biomarkers - a validation study

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

IsomiR-eQTL: A Cancer-Specific Expression Quantitative Trait Loci Database of miRNAs and Their Isoforms

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.

Obesity and Roux-en-Y gastric bypass drive changes in miR-31 and miR-215 expression in the human rectal mucosa

Background/Objectives

Obesity increases colorectal cancer (CRC) risk. However, the effects of weight loss on CRC risk are unclear. Epigenetic mechanisms involving microRNAs that lead to dysregulated gene expression may mediate the effects of obesity and weight loss on CRC risk. We examined the effects of obesity and weight loss following Roux-en-Y gastric bypass (RYGB) on microRNA expression in the human rectal mucosa.

Subjects/Methods

We collected rectal mucosal biopsies from obese patients (n = 22) listed for RYGB and age- and sex-matched healthy non-obese Controls (n = 20), at baseline and six months post-surgery. We quantified microRNA expression in rectal mucosal biopsies using Next Generation Sequencing and bioinformatics analysis to investigate the likely functional consequences of these epigenetic changes.

Results

Compared with non-obese individuals, obese individuals showed differential expression of 112 microRNAs (p < 0.05). At six-months post-RYGB, when mean body mass had fallen by 27 kg, 60 microRNAs were differentially expressed, compared with baseline (p < 0.05). The expression of 36 microRNAs differed significantly between both i) obese and non-obese individuals and ii) obese individuals pre- and post-RYGB. Quantitative polymerase chain reaction (qPCR) demonstrated that expression of miR-31 and miR-215 was significantly (p < 0.05) higher, 143-fold and 15-fold respectively, in obese than in non-obese individuals. Weight loss, following RYGB, reduced expression of miR-31 and miR-215 to levels comparable with Controls. These differentially expressed microRNAs are implicated in pathways linked with inflammation, obesity and cancer.

Conclusion

Our findings show, for the first time, that obesity is associated with dysregulated microRNA expression in the human rectal mucosa. Further, surgically-induced weight loss may normalise microRNA expression in this tissue.

Identification and Elucidation of the Protective isomiRs in Lung Cancer Patient Prognosis

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.

Molecular Mechanisms of Chemoresistance Induced by Cisplatin in NSCLC Cancer Therapy

Cancer cells utilise several mechanisms to increase their survival and progression as well as their resistance to anticancer therapy: deregulation of growth regulatory pathways by acquiring grow factor independence, immune system suppression, reducing the expression of antigens activating T lymphocyte cells (mimicry), induction of anti-apoptotic signals to counter the action of drugs, activation of several DNA repair mechanisms and driving the active efflux of drugs from the cell cytoplasm, and epigenetic regulation by microRNAs (miRNAs). Because it is commonly diagnosed late, lung cancer remains a major malignancy with a low five-year survival rate; when diagnosed, the cancer is often highly advanced, and the cancer cells may have acquired drug resistance. This review summarises the main mechanisms involved in cisplatin resistance and interactions between cisplatin-resistant cancer cells and the tumour microenvironment. It also analyses changes in the gene expression profile of cisplatin sensitive vs. cisplatin-resistant non-small cell lung cancer (NSCLC) cellular model using the GSE108214 Gene Expression Omnibus database. It describes a protein-protein interaction network that indicates highly dysregulated TP53, MDM2, and CDKN1A genes as they encode the top networking proteins that may be involved in cisplatin tolerance, these all being upregulated in cisplatin-resistant cells. Furthermore, it illustrates the multifactorial nature of cisplatin resistance by examining the diversity of dysregulated pathways present in cisplatin-resistant NSCLC cells based on KEGG pathway analysis.

Regulation of breast cancer metastasis signaling by miRNAs

Despite the decline in death rate from breast cancer and recent advances in targeted therapies and combinations for the treatment of metastatic disease, metastatic breast cancer remains the second leading cause of cancer-associated death in U.S. women. The invasion-metastasis cascade involves a number of steps and multitudes of proteins and signaling molecules. The pathways include invasion, intravasation, circulation, extravasation, infiltration into a distant site to form a metastatic niche, and micrometastasis formation in a new environment. Each of these processes is regulated by changes in gene expression. Noncoding RNAs including microRNAs (miRNAs) are involved in breast cancer tumorigenesis, progression, and metastasis by post-transcriptional regulation of target gene expression. miRNAs can stimulate oncogenesis (oncomiRs), inhibit tumor growth (tumor suppressors or miRsupps), and regulate gene targets in metastasis (metastamiRs). The goal of this review is to summarize some of the key miRNAs that regulate genes and pathways involved in metastatic breast cancer with an emphasis on estrogen receptor α (ERα+) breast cancer. We reviewed the identity, regulation, human breast tumor expression, and reported prognostic significance of miRNAs that have been documented to directly target key genes in pathways, including epithelial-to-mesenchymal transition (EMT) contributing to the metastatic cascade. We critically evaluated the evidence for metastamiRs and their targets and miRNA regulation of metastasis suppressor genes in breast cancer progression and metastasis. It is clear that our understanding of miRNA regulation of targets in metastasis is incomplete.

Key MicroRNA's and Their Targetome in Adrenocortical Cancer

Adrenocortical Carcinoma (ACC) is a rare but aggressive malignancy with poor prognosis and limited response to available systemic therapies. Although complete surgical resection gives the best chance for long-term survival, ACC has a two-year recurrence rate of 50%, which poses a therapeutic challenge. High throughput analyses focused on characterizing the molecular signature of ACC have revealed specific micro-RNAs (miRNAs) that are associated with aggressive tumor phenotypes. MiRNAs are small non-coding RNA molecules that regulate gene expression by inhibiting mRNA translation or degrading mRNA transcripts and have been generally implicated in carcinogenesis. This review summarizes the current insights into dysregulated miRNAs in ACC tumorigenesis, their known functions, and specific targetomes. In addition, we explore the possibility of particular miRNAs to be exploited as clinical biomarkers in ACC and as potential therapeutics.

Circulating microRNA/isomiRs as novel biomarkers of esophageal squamous cell carcinoma

BACKGROUND

MicroRNA (miR)s are promising diagnostic biomarkers of cancer. Recent next generation sequencer (NGS) studies have found that isoforms of micro RNA (isomiR) circulate in the bloodstream similarly to mature micro RNA (miR). We hypothesized that combination of circulating miR and isomiRs detected by NGS are potentially powerful cancer biomarker. The present study aimed to investigate their application in esophageal cancer.

METHODS

Serum samples from patients with esophageal squamous cell carcinoma (ESCC) and age and sex matched healthy control (HC) individuals were investigated for the expression of miR/isomiRs using NGS. Candidate miR/isomiRs which met the criteria in the 1st group (ESCC = 18 and HC = 12) were validated in the 2nd group (ESCC = 30 and HC = 30). A diagnostic panel was generated using miR/isomiRs that were consistently confirmed in the 1st and 2nd groups. Accuracy of the panel was tested then in the 3rd group (ESCC = 18 and HC = 18). Their use was also investigated in 22 paired samples obtained pre- and post-treatment, and in patients with esophageal adenocarcinoma (EAD) and high-grade dysplasia (HGD).

RESULTS

Twenty-four miR/isomiRs met the criteria for diagnostic biomarker in the 1st and 2nd group. A multiple regression model selected one mature miR (miR-30a-5p) and two isomiRs (isoform of miR-574-3p and miR-205-5p). The index calculated from the diagnostic panel was significantly higher in ESCC patients than in the HCs (13.3±8.9 vs. 3.1±1.3, p<0.001). The area under the receiver operating characteristics (ROC) curves of the panel index was 0.95. Sensitivity and specificity were 93.8%, and 81% in the 1st and 2nd groups, and 88.9% and 72.3% in the 3rd group, respectively. The panel index was significantly lower in patients with EAD (6.2±4.5) and HGD (4.2±1.7) than in those with ESCC and was significantly decreased at post-treatment compared with pre-treatment (6.2±5.6 vs 11.6±11.5, p = 0.03).

CONCLUSION

Our diagnostic panel had high accuracy in the diagnosis of ESCC. MiR/isomiRs detected by NGS could serve as novel biomarkers of ESCC.

The role and mechanisms of action of microRNAs in cancer drug resistance

MicroRNAs (miRNAs) are small non-coding RNAs with a length of about 19–25 nt, which can regulate various target genes and are thus involved in the regulation of a variety of biological and pathological processes, including the formation and development of cancer. Drug resistance in cancer chemotherapy is one of the main obstacles to curing this malignant disease. Statistical data indicate that over 90% of the mortality of patients with cancer is related to drug resistance. Drug resistance of cancer chemotherapy can be caused by many mechanisms, such as decreased antitumor drug uptake, modified drug targets, altered cell cycle checkpoints, or increased DNA damage repair, among others. In recent years, many studies have shown that miRNAs are involved in the drug resistance of tumor cells by targeting drug-resistance-related genes or influencing genes related to cell proliferation, cell cycle, and apoptosis. A single miRNA often targets a number of genes, and its regulatory effect is tissue-specific. In this review, we emphasize the miRNAs that are involved in the regulation of drug resistance among different cancers and probe the mechanisms of the deregulated expression of miRNAs. The molecular targets of miRNAs and their underlying signaling pathways are also explored comprehensively. A holistic understanding of the functions of miRNAs in drug resistance will help us develop better strategies to regulate them efficiently and will finally pave the way toward better translation of miRNAs into clinics, developing them into a promising approach in cancer therapy.

Baptista P, Mullin M, Dalby MJ, Berry CC. Nanoparticle-antagomiR based targeting of miR-31 to induce osterix and osteocalcin expression in mesenchymal stem cells

Mesenchymal stem cells are multipotent adult stem cells capable of generating bone, cartilage and fat, and are thus currently being exploited for regenerative medicine. When considering osteogenesis, developments have been made with regards to chemical induction (e.g. differentiation media) and physical induction (e.g. material stiffness, nanotopography), targeting established early transcription factors or regulators such as runx2 or bone morphogenic proteins and promoting increased numbers of cells committing to osteo-specific differentiation. Recent research highlighted the involvement of microRNAs in lineage commitment and terminal differentiation. Herein, gold nanoparticles that confer stability to short single stranded RNAs were used to deliver MiR-31 antagomiRs to both pre-osteoblastic cells and primary human MSCs in vitro. Results showed that blocking miR-31 led to an increase in osterix protein in both cell types at day 7, with an increase in osteocalcin at day 21, suggesting MSC osteogenesis. In addition, it was noted that antagomiR sequence direction was important, with the 5 prime reading direction proving more effective than the 3 prime. This study highlights the potential that miRNA antagomiR-tagged nanoparticles offer as novel therapeutics in regenerative medicine.

MicroRNAs as regulators of cisplatin-resistance in non-small cell lung carcinomas

With more than 80% of all diagnosed lung cancer cases, non-small cell lung cancer (NSCLC) remains the leading cause of cancer death worldwide. Exact diagnosis is mostly very late and advanced-stage NSCLCs are inoperable at admission. Tailored therapies with tyrosine kinase inhibitors are only available for a minority of patients. Thus, chemotherapy is often the treatment of choice. As first-line chemotherapy for NSCLCs, platinum-based substances (e.g. cisplatin, CDDP) are mainly used. Unfortunately, the positive effects of CDDP are frequently diminished due to development of drug resistance and negative influence of microenvironmental factors like hypoxia. MicroRNAs (miRNAs) are small, non-coding molecules involved in the regulation of gene expression and modification of biological processes like cell proliferation, apoptosis and cell response to chemotherapeutics. Expression of miRNAs is often deregulated in lung cancer compared to corresponding non-malignant tissue. In this review we summarize the present knowledge about the effects of miRNAs on CDDP-resistance in NSCLCs. Further, we focus on miRNAs deregulated by hypoxia, which is an important factor in the development of CDDP-resistance in NSCLCs. This review will contribute to the general understanding of miRNA-regulated biological processes in NSCLC, with special focus on the role of miRNA in CDDP-resistance.

Naturally existing isoforms of miR-222 have distinct functions

Deep-sequencing reveals extensive variation in the sequence of endogenously expressed microRNAs (termed ‘isomiRs’) in human cell lines and tissues, especially in relation to the 3′ end. From the immunoprecipitation of the microRNA-binding protein Argonaute and the sequencing of associated small RNAs, we observe extensive 3′-isomiR variation, including for miR-222 where the majority of endogenously expressed miR-222 is extended by 1–5 nt compared to the canonical sequence. We demonstrate this 3′ heterogeneity has dramatic implications for the phenotype of miR-222 transfected cells, with longer isoforms promoting apoptosis in a size (but not 3′ sequence)-dependent manner. The transfection of longer miR-222 isomiRs did not induce an interferon response, but did downregulate the expression of many components of the pro-survival PI3K-AKT pathway including PIK3R3, a regulatory subunit whose knockdown phenocopied the expression of longer 222 isoforms in terms of apoptosis and the inhibition of other PI3K-AKT genes. As this work demonstrates the capacity for 3′ isomiRs to mediate differential functions, we contend more attention needs to be given to 3′ variance given the prevalence of this class of isomiR.

IsomiR processing during differentiation of myelogenous leukemic cell line K562 by phorbol ester PMA

Chronic myelocytic leukemia cell line K562 undergoes differentiation by phorbol esters to megakaryocytes and we have used this system to understand miRNA processing leading to isomiR generation. PMA treatment significantly altered the production of miRNA in K562 cells. Expression of 24.4% of miRNAs were found to be stimulated whereas expression of 10% miRNAs were inhibited by PMA treatment. Our results suggest that miRNA precursors are processed into isomiRs in a deterministic manner. The relative levels of different isomiRs of a miRNA remained mainly unchanged even after PMA treatment irrespective of overall changes in expression (either up-regulation or down-regulation). However, not all miRNAs behave in the same way, about 7% showed a variation of isomiR profiles after PMA treatment. Most of the later class of miRNAs were found to be oncogenic miRNAs. Further, it was also found that number of isomiRs was independent of abundance of a miRNA. Functional importance of different isomiRs was demonstrated using three different isomiRs of miR-22. Our results showed that different isomiRs could inhibit expression of targets genes with different efficiencies. Our study suggests that the heterogeneity of a miRNA population generated during processing is in general regulated and that variation in the generation of an isomiR can be a functionally important regulatory feature.

Identification of aminosulfonylarylisoxazole as microRNA-31 regulators

The discovery of small-molecule regulators of microRNAs remains challenging, but a few have been reported. Herein, we describe small-molecule inhibitors of miR-31, a tumor-associated microRNA (miRNA), identified by high-throughput screening using a cell-based reporter assay. Aminosulfonylarylisoxazole compounds exhibited higher specificity for miR-31 than for six other miRNAs, i.e., miR-15a, miR-16, miR-21, miR-92a-1, miR-146a, and miR-155, and increased the expression of miR-31 target genes. The down-regulation of mature miR-31 was observed, while its precursor form increased following treatment with the compounds. Thus, the compounds may target the processing of pre-miR-31 into mature miR-31 and thereby inhibit the production of mature miR-31.

The highly expressed 5'isomiR of hsa-miR-140-3p contributes to the tumor-suppressive effects of miR-140 by reducing breast cancer proliferation and migration

Background

miRNAs are small noncoding RNA molecules that play an important role in post-transcriptional regulation of gene expression. Length and/or sequence variants of the same miRNA are termed isomiRs. While most isomiRs are functionally redundant compared to their canonical counterparts, the so-called 5’isomiRs exhibit a shifted 5’ end and therefore a shifted seed sequence resulting in a different target spectrum. However, not much is known about the functional relevance of these isoforms.

Results

Analysis of miRNA-seq data from breast cancer cell lines identified six pairs of highly expressed miRNAs and associated 5’isomiRs. Among them, hsa-miR-140-3p was of particular interest because its 5’isomiR showed higher expression compared to the canonical miRNA annotated in miRbase. This miRNA has previously been shown to control stemness of breast cancer cells. miRNAseq data of breast cancer patients (TCGA dataset) showed that both the canonical hsa-miR-140-3p and its 5’isomiR-140-3p were highly expressed in patients’ tumors compared to normal breast tissue. In the current work, we present the functional characterization of 5’isomiR-140-3p and the cellular phenotypes associated with its overexpression in MCF10A, MDA-MB-468 and MDA-MB-231 cell lines in comparison to the canonical hsa-miR-140-3p. Contrary to the effect of the canonical hsa-miR-140-3p, overexpression of the 5’isomiR-140-3p led to a decrease in cell viability. The latter observation was supported by cell cycle analysis, where the 5’isomiR-140-3p but not the hsa-miR-140-3p caused cell cycle arrest in G₀/G₁-phase. Additionally, 5’ismoiR-140-3p overexpression was found to cause a decrease in cell migration in the three cell lines. We identified three novel direct target genes of the 5’isomiR-140-3p; COL4A1, ITGA6 and MARCKSL1. Finally, we have shown that knocking down these genes partially phenocopied the effects of the 5’isomiR-140-4p overexpression, where COL4A1 and ITGA6 knockdown led to reduced cell viability and cell cycle arrest, while MARCKSL1 knockdown resulted in a decrease in the migratory potential of cells.

Conclusions

In summary, this work presents evidence that there is functional synergy between the canonical hsa-miR-140-3p and the newly identified 5’isomiR-140-3p in suppressing growth and progression of breast cancer by simultaneously targeting genes related to differentiation, proliferation, and migration.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2869-x) contains supplementary material, which is available to authorized users.

Dual regulation by microRNA-200b-3p and microRNA-200b-5p in the inhibition of epithelial-to-mesenchymal transition in triple-negative breast cancer

Epithelial to mesenchymal transition (EMT) involves loss of an epithelial phenotype and activation of a mesenchymal one. Enhanced expression of genes associated with a mesenchymal transition includes ZEB1/2, TWIST, and FOXC1. miRNAs are known regulators of gene expression and altered miRNA expression is known to enhance EMT in breast cancer. Here we demonstrate that the tumor suppressive miRNA family, miR-200, is not expressed in triple negative breast cancer (TNBC) cell lines and that miR-200b-3p over-expression represses EMT, which is evident through decreased migration and increased CDH1 expression. Despite the loss of migratory capacity following re-expression of miR-200b-3p, no subsequent loss of the conventional miR-200 family targets and EMT markers ZEB1/2 was observed. Next generation RNA-sequencing analysis showed that enhanced expression of pri-miR-200b lead to ectopic expression of both miR-200b-3p and miR-200b-5p with multiple isomiRs expressed for each of these miRNAs. Furthermore, miR-200b-5p was expressed in the receptor positive, epithelial breast cancer cell lines but not in the TNBC (mesenchymal) cell lines. In addition, a compensatory mechanism for miR-200b-3p/200b-5p targeting, where both miRNAs target the RHOGDI pathway leading to non-canonical repression of EMT, was demonstrated. Collectively, these data are the first to demonstrate dual targeting by miR-200b-3p and miR-200b-5p and a previously undescribed role for microRNA processing and strand expression in EMT and TNBC, the most aggressive breast cancer subtype.

Analysis of Biobanked Serum from a Mycobacterium avium subsp paratuberculosis Bovine Infection Model Confirms the Remarkable Stability of Circulating miRNA Profiles and Defines a Bovine Serum miRNA Repertoire

Johne's Disease (JD) is a chronic enteritis of ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP). Current disease control strategies are hampered by the lack of sensitive and specific diagnostic modalities. Therefore, novel diagnostic and prognostic tools are needed, and circulating microRNAs (miRNAs) may hold potential in this area. The aims of this study were twofold: (i) to address the stability of miRNA in bovine sera from biobanked samples, and (ii) to assess the potential of miRNAs as biomarkers for JD disease progression. To address these aims we used bovine sera from an experimental MAP infection model that had been stored at -20°C for over a decade, allowing us to also assess the stability of miRNA profiles in biobanked serum samples through comparison with fresh sera. Approximately 100-200 intact miRNAs were identified in each sample with 83 of these being consistently detected across all 57 samples. The miRNA profile of the biobanked sera stored at -20°C for over 10 years was highly similar to the profile of <1 year-old sera stored at -80°C, with an overlap of 73 shared miRNAs. IsomiR analysis also indicated a distinct bovine serum-specific isomiR profile as compared to previously reported bovine macrophage miRNA profiles. To explore the prognostic potential of miRNA profiles cattle defined as seropositive for anti-MAP antibodies (n = 5) were compared against seronegative cattle (n = 7). No significant differential expressed miRNAs were detected at either the early (6 months) or late (43, 46 and 49 months) intervals (FDR≤0.05, fold-change≥1.5) across seropositive or seronegative animals. However, comparing pre-infection sera to the early and late time-points identified increased miR-29a and miR-92b abundance (2-fold) that may be due to blood-cell population changes over time (P<0.001). In conclusion our study has demonstrated that bovine circulating miRNAs retain their integrity under long-term sub-optimal storage temperatures opening the way for increased miRNA analyses from biobanked samples for a range of infectious and non-infectious diseases.

Dengue virus infection alters post-transcriptional modification of microRNAs in the mosquito vector Aedes aegypti

Recent discoveries regarding the importance of isomiRs have increased our understanding of the regulatory complexities of the miRNAome. Observed changes in the miRNA profiles in mosquitoes infected with flaviviruses have implicated small RNAs in the interactions between viruses and their vectors. Here we analysed the isomiR profiles of both uninfected and infected Aedes aegypti mosquitoes with the major human pathogen dengue virus (DENV). We found that several specific isomiRs were significantly altered in their abundance patterns in response to DENV infection potentially affecting their target repertoire. Notable among these were isomiR variants which displayed arm-switching. We also demonstrate that modifications to the 3p end of miRNAs are vastly more prevalent than those at the 5p ends. We also observed that in only 45% of Ae. aegypti miRNAs the most abundant read matches the exact sequence reported in miRBase. Further, we found positive correlations between the number of mature miRNA reads, pre-miRNA length, GC content and secondary structure minimum free energy with the number of isomiRs. The findings presented here provide some evidence that isomiR production is not a random phenomenon and may be important in DENV replication in its vector.

Small RNA Sequencing Uncovers New miRNAs and moRNAs Differentially Expressed in Normal and Primary Myelofibrosis CD34+ Cells

Myeloproliferative neoplasms (MPN) are chronic myeloid cancers thought to arise at the level of CD34+ hematopoietic stem/progenitor cells. They include essential thrombocythemia (ET), polycythemia vera (PV) and primary myelofibrosis (PMF). All can progress to acute leukemia, but PMF carries the worst prognosis. Increasing evidences indicate that deregulation of microRNAs (miRNAs) might plays an important role in hematologic malignancies, including MPN. To attain deeper knowledge of short RNAs (sRNAs) expression pattern in CD34+ cells and of their possible role in mediating post-transcriptional regulation in PMF, we sequenced with Illumina HiSeq2000 technology CD34+ cells from healthy subjects and PMF patients. We detected the expression of 784 known miRNAs, with a prevalence of miRNA up-regulation in PMF samples, and discovered 34 new miRNAs and 99 new miRNA-offset RNAs (moRNAs), in CD34+ cells. Thirty-seven small RNAs were differentially expressed in PMF patients compared with healthy subjects, according to microRNA sequencing data. Five miRNAs (miR-10b-5p, miR-19b-3p, miR-29a-3p, miR-379-5p, and miR-543) were deregulated also in PMF granulocytes. Moreover, 3’-moR-128-2 resulted consistently downregulated in PMF according to RNA-seq and qRT-PCR data both in CD34+ cells and granulocytes. Target predictions of these validated small RNAs de-regulated in PMF and functional enrichment analyses highlighted many interesting pathways involved in tumor development and progression, such as signaling by FGFR and DAP12 and Oncogene Induced Senescence. As a whole, data obtained in this study deepened the knowledge of miRNAs and moRNAs altered expression in PMF CD34+ cells and allowed to identify and validate a specific small RNA profile that distinguishes PMF granulocytes from those of normal subjects. We thus provided new information regarding the possible role of miRNAs and, specifically, of new moRNAs in this disease.

Genetic variants in microRNA genes: impact on microRNA expression, function, and disease

MicroRNAs (miRNAs) are important regulators of gene expression and like any other gene, their coding sequences are subject to genetic variation. Variants in miRNA genes can have profound effects on miRNA functionality at all levels, including miRNA transcription, maturation, and target specificity, and as such they can also contribute to disease. The impact of variants in miRNA genes is the focus of the present review. To put these effects into context, we first discuss the requirements of miRNA transcripts for maturation. In the last part an overview of available databases and tools and experimental approaches to investigate miRNA variants related to human disease is presented.

MicroRNA and Breast Cancer: Understanding Pathogenesis, Improving Management

The advent of the microRNAs in the early 1990s has proven to be a tremendously significant development within the purview of gene regulation. They participate in the regulation of a broad assembly of processes vital to proper cell function and the perturbation of these pathways following alteration of miRNA expression is strongly believed to contribute to the pathogenesis of cancer. This review provides a comprehensive overview of the miRNAs that have to date been well-characterized in the context of human breast neoplasia. Detailed discussion will center around their role in tumor initiation and progression, control of epithelial-mesenchymal transition (EMT), cancer stem cell formation, use as biomarkers in tissues and circulation, as well as their role in cancer treatment. In addition, attention will be given to topics which remain underexplored, such as miRNA control of cancer cell metabolism and the genomic/epigenetic origins underlying the preliminary disruption of miRNA expression in disease. This review will also address and attempt to resolve instances where discordant, inter-study findings have been reported (examples of which are replete in the literature) while also identifying bottlenecks hampering progress in miRNA research and other challenges that confront this fledgling but promising field of biomedical research.

Expression of dicer and its related miRNAs in the progression of prostate cancer

Dicer is aberrantly expressed in several types of malignancies. Cleaved by Dicer, the small noncoding microRNAs (miRNAs) are considered potential tools for the diagnosis and prognosis of cancer. This study investigated the expression of miRNAs thought to target Dicer. Expression of 1,205 human miRNAs and miRNA*s were examined in four patients with prostate cancer (PCa) by miRNA array in which the threshold was set as two-fold. Seventy-three miRNAs and miRNA*s were significantly down-regulated while 10 were up-regulated in PCa tissues compared with matched histologically normal glands. Of these, miR-29b-1, miR-200a, miR-370, and miR-31, which were the most down/up-regulated and closely potentially target to the Dicer 3' UTR, were investigated further. Tissues of primary tumors and matched normal prostate glands from 185 patients with PCa were collected for further investigation. Dicer mRNA levels were negatively correlated with miR-29b-1 (ρs = -0.177, p = 0.017), miR-200a (ρs = -0.489, p < 0.0001) and miR-31 (ρs = -0.314, p < 0.0001) expression. Compared with adjacent normal glands, PCa tissues showed significantly lower miR-200a and miR-31 expression levels. Furthermore, in metastatic PCa, the expression levels of miR-200a, miR-370, and miR-31 were dramatically higher than in localized PCa. Additionally, elevated expression levels of miR-200a and miR-31 appeared to be associated with castration-resistant PCa. These findings suggest possibilities that miR-200a and miR-31 target Dicer and are involved in the carcinogenesis, migration, and behavior of castration-resistant PCa, indicating that they could be potential biomarkers for monitoring PCa progression.

Comparisons of isomiR patterns and classification performance using the rank-based MANOVA and 10-fold cross-validation

Next generation sequencing technology has identified a series of miRNA variants (named "isomiRs"), which might be associated with cancer progression. We provide a new strategy to reanalyze the miR-seq datasets through a view of the isomiR spectrum. Firstly, differentially expressed (DE) isomiRs were detected with the DESeq algorithm based on negative binomial distribution. Secondly, the rank-based MANOVA was adopted to compare the isomiR patterns between normal and tumor tissues. Moreover, a comprehensive survey on classification performance of three features was conducted, including the logistic regression, k-nearest neighbors and Random Forest. Finally, functional enrichment analysis was performed with the putative targets of specific isomiRs to elucidate their biological functions. Furthermore, the methods were applied to the downloaded miR-seq datasets of breast invasive carcinoma from TCGA. We found that the expression levels of multiple isomiRs derived from the same miRNA locus showed significant inconsistency between normal and tumor samples. In most cases, logistic regression with multiple DE isomiRs was superior to the others, with highest AUC and lowest AIC. Similarly, DE isomiRs performed best in the average accuracy of standard classifiers. Integrated targets were significantly enriched in some cancer-related pathways, including MAPK signaling pathway, and focal adhesion. Collectively, we could recommend the rank-based MANOVA for comparing different isomiR patterns, and further investigation on isomiRs needs to be considered in miRNA sequencing research.

Insight into small RNA abundance and expression in high- and low-temperature stress response using deep sequencing in Arabidopsis

Small RNA profiling and assessing its dependence on changing environmental factors have expanded our understanding of the transcriptional and post-transcriptional regulation of plant stress responses. Insufficient data have been documented earlier to depict the profiling of small RNA classes in temperature-associated stress which has a wide implication for climate change biology. In the present study, we report a comparative assessment of the genome-wide profiling of small RNAs in Arabidopsis thaliana using two conditional responses, induced by high- and low-temperature. Genome-wide profiling of small RNAs revealed an abundance of 21 nt small RNAs at low temperature, while high temperature showed an abundance of 21 nt and 24 nt small RNAs. The two temperature treatments altered the expression of a specific subset of mature miRNAs and displayed differential expression of a number of miRNA isoforms (isomiRs). Comparative analysis demonstrated that a large number of protein-coding genes can give rise to differentially expressed small RNAs following temperature shifts. Low temperature caused accumulation of small RNAs, corresponding to the sense strand of a number of cold-responsive genes. In contrast, high temperature stimulated the production of small RNAs of both polarities from genes encoding functionally diverse proteins.

MicroRNA-5p and -3p co-expression and cross-targeting in colon cancer cells

Background

Two mature miRNA species may be generated from the 5’ and 3’ arms of a pre-miRNA precursor. In most cases, only one species remains while the complementary species is degraded. However, co-existence of miRNA-5p and -3p species is increasingly being reported. In this work, we aimed to systematically investigate co-expression of miRNA-5p/3p in colon cancer cells in a genome-wide analysis, and to examine cross-targeting of the dysregulated miRNAs and 5p/3p species.

Results

Four colon cancer cell lines were examined relative to two normal colon tissues. Of the 1,190 miRNAs analyzed, 92 and 36 were found to be up- or down-regulated, respectively, in cancer cells. Nineteen co-expressed miRNA-5p/3p pairs were further identified suggesting frequent 5p/3p co-accumulation in colon cancer cells. Of these, 14 pairs were co-up-regulated and 3 pairs were co-down-regulated indicating concerted 5p/3p dysregulation. Nine dysregulated miRNA pairs fell into three miRNA gene families, namely let-7, mir-8/200 and mir-17, which showed frequent cross-targeting in the metastasis process. Focusing on the let-7d-5p/3p pair, the respectively targeted IGF1R and KRAS were shown to be in a reverse relationship with expression of the respective miRNA, which was confirmed in transient transfection assays using let-7d mimic or inhibitor. Targeting of KRAS by let-7d was previous reported; targeting of IGF1R by let-7d-5p was confirmed in luciferase assays in this study. The findings of let-7d-5p/3p and multiple other miRNAs targeting IGF1R, KRAS and other metastasis-related factors suggest that 5p/3p miRNAs contribute to cross-targeting of multiple cancer-associated factors and processes possibly to evade functional abolishment when any one of the crucial factors are inactivated.

Conclusions

miRNA-5p/3p species are frequently co-expressed and are coordinately regulated in colon cancer cells. In cancer cells, multiple cross-targeting by the miRNAs, including the co-existing 5p/3p species, frequently occurs in an apparent safe-proof scheme of miRNA regulation of important tumorigenesis processes. Further systematic analysis of co-existing miRNA-5p/3p pairs in clinical tissues is important in elucidating 5p/3p contributions to cancer pathogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-014-0095-x) contains supplementary material, which is available to authorized users.

IsomiR expression profiles in human lymphoblastoid cell lines exhibit population and gender dependencies

For many years it was believed that each mature microRNA (miRNA) existed as a single entity with fixed endpoints and a 'static' and unchangeable primary sequence. However, recent evidence suggests that mature miRNAs are more 'dynamic' and that each miRNA precursor arm gives rise to multiple isoforms, the isomiRs. Here we report on our identification of numerous and abundant isomiRs in the lymphoblastoid cell lines (LCLs) of 452 men and women from five different population groups. Unexpectedly, we find that these isomiRs exhibit an expression profile that is population-dependent and gender-dependent. This is important as it indicates that the LCLs of each gender/population combination have their own unique collection of mature miRNA transcripts. Moreover, each identified isomiR has its own characteristic abundance that remains consistent across biological replicates indicating that these are not degradation products. The primary sequences of identified isomiRs differ from the known miRBase miRNA either at their 5´-endpoint (leads to a different 'seed' sequence and suggests a different targetome), their 3´-endpoint, or both simultaneously. Our analysis of Argonaute PAR-CLIP data from LCLs supports the association of many of these newly identified isomiRs with the Argonaute silencing complex and thus their functional roles through participation in the RNA interference pathway.

Role of microRNAs and DNA Methyltransferases in Transmitting Induced Genomic Instability between Cell Generations

There is limited understanding of how radiation or chemicals induce genomic instability, and how the instability is epigenetically transmitted to the progeny of exposed cells or organisms. Here, we measured the expression of microRNAs (miRNAs) and DNA methyltransferases (DNMTs) in murine embryonal fibroblasts exposed to ionizing radiation or 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), which were previously shown to induce genomic instability in this cell line. Cadmium was used as a reference agent that does not induce genomic instability in our experimental model. Measurements at 8 and 15 days after exposure did not identify any such persistent changes that could be considered as signals transmitting genomic instability to the progeny of exposed cells. However, measurements at 2 days after exposure revealed findings that may reflect initial stages of genomic instability. Changes that were common to TCDD and two doses of radiation (but not to cadmium) included five candidate signature miRNAs and general up-regulation of miRNA expression. Expression of DNMT3a, DNMT3b, and DNMT2 was suppressed by cadmium but not by TCDD or radiation, consistently with the hypothesis that sufficient expression of DNMTs is necessary in the initial phase of induced genomic instability.

K14-EGFP-miR-31 transgenic mice have high susceptibility to chemical-induced squamous cell tumorigenesis that is associating with Ku80 repression

Squamous cell carcinoma (SCC) occurring in the head and neck region and the esophagus causes tremendous cancer mortality around the world. miR-31 is among the most eminently upregulated MicroRNAs in SCC, when it occurs in the head and neck region and the esophagus. We established miR-31 transgenic mouse lines, in which miR-31 is under the control of the K14 promoter. 4-nitroquinoline 1-oxide (4NQO) is a mutagen that causes double strand breaks. The transgenic mice exhibited a higher potential for tumor induction than wild-type (Wt) mice of the tongue and esophagus after 4NQO treatment. After 4NQO treatment or irradiation, p-γH2AX expression in squamous epithelium of transgenic mice was increased more than in Wt mice. Exogenous expression of miR-31 was also found to be associated with the higher p-γH2AX expression induced by 4NQO in human oral SCC (OSCC) cell lines. The repair genes PARP1 and Ku80 were validated as new targets of miR-31 in human OSCC cell lines, and were found to be downregulated in the squamous epithelium of the tongue in transgenic mice. However, only the downregulation of Ku80 was essential for maintaining the high level of p-γH2AX induced by 4NQO in OSCC cells. Inverse expression profiles for miR-31 and Ku80 were noted in human OSCC tissue. Our study identifies the high sensitivity of K14-EGFP-miR-31 transgenic mice to chemical carcinogen-induced squamous cell tumorigenesis and shows that this seems to be associated with the downregulation of Ku80 and an impairment of repair activity in squamous cells, which are mediated by miR-31.

A challenge for miRNA: multiple isomiRs in miRNAomics

Accumulating evidence suggests that a single microRNA (miRNA) locus can generate a series of sequences during miRNA maturation process. These multiple sequences, called miRNA variants, or isomiRs, have different lengths and different 5' and 3' ends. Some of these isomiRs are detected as varied nucleotides and 3' additional non-template nucleotides. As physiological miRNA isoforms, they have drawn attention for possible regulatory biological roles. The present work mainly reviews miRNA/isomiR biogenesis, isomiR expression patterns, and functional and evolutionary implications, especially between isomiRs from homologous and clustered miRNA loci. The phenomenon of multiple isomiRs and their biological roles indicates that analysis performed at the miRNA and isomiR levels should be included in miRNA studies. This may enrich and complicate miRNA biogenesis and coding-non-coding RNA regulatory networks.

The panorama of miRNA-mediated mechanisms in mammalian cells

MicroRNAs comprise a large family of short, non-coding RNAs that are present in most eukaryotic organisms and are typically involved in downregulating the expression of protein-coding genes. The detailed mechanisms of miRNA functioning in animals and plants have been under investigation for more than decade. In mammalian cells, miRNA guides the effector complex miRISC to bind with partially complementary sequences, usually within the 3′UTR of mRNAs, and inhibit protein synthesis with or without transcript degradation. In addition to these main mechanisms, several other modes of miRNA-mediated gene expression regulation have been described, but their scale and importance remain a matter of debate. In this review, we briefly summarize the pathway of miRNA precursor processing during miRNA biogenesis and continue with the description of the miRISC assembly process. Then, we present the miRNA-mediated mechanisms of gene expression regulation in detail, and we gather information concerning the proteins involved in these processes. In addition, we briefly refer to the current applications of miRNA mechanisms in therapeutic strategies. Finally, we highlight some of the remaining controversies surrounding the regulation of mammalian gene expression by miRNAs.

Regulation of breast cancer and bone metastasis by microRNAs

Breast cancer progression including bone metastasis is a complex process involving numerous changes in gene expression and function. MicroRNAs (miRNAs) are small endogenous noncoding RNAs that regulate gene expression by targeting protein-coding mRNAs posttranscriptionally, often affecting a number of gene targets simultaneously. Alteration in expression of miRNAs is common in human breast cancer, possessing with either oncogenic or tumor suppressive activity. The expression and the functional role of several miRNAs (miR-206, miR-31, miR-27a/b, miR-21, miR-92a, miR-205, miR-125a/b, miR-10b, miR-155, miR-146a/b, miR-335, miR-204, miR-211, miR-7, miR-22, miR-126, and miR-17) in breast cancer has been identified. In this review we summarize the experimentally validated targets of up- and downregulated miRNAs and their regulation in breast cancer and bone metastasis for diagnostic and therapeutic purposes.

The diverse role of miR-31 in regulating cancer associated phenotypes

In the past 10 years research on miRNAs has demonstrated their central role in regulating gene expression both in normal and diseased tissue. The expression of miRNAs is widely altered in cancer, leading to abnormal expression of the genes regulated by these miRNAs, and subsequently alterations in entire molecular networks and pathways. One especially interesting cancer-related miRNA is miR-31 which is frequently altered in a large variety of cancers. The functional role of miR-31 is extremely complex and miR-31 can hold both tumor suppressive and oncogenic roles in different tumor types. The phenotype caused by aberrant miR-31 expression seems to be strongly dependent on the endogenous expression levels. For example, in breast cancer loss of miR-31 expression is associated with high risk of metastases, whereas in colorectal cancer high miR-31 expression correlates with advanced disease stage. This review summarizes the complex expression patterns of miR-31 in human cancers, describes the variable phenotypes caused by altered miR-31 expression, and highlights the current knowledge on the genes targeted by miR-31.

Regulation of huntingtin gene expression by miRNA-137, -214, -148a, and their respective isomiRs

With the advent of deep sequencing technology, a variety of miRNA length and sequence variants, termed isomiRNAs (isomiRs), have been discovered. However, the functional roles of these commonly detected isomiRs remain unknown. In this paper, we demonstrated that miRNAs regulate the expression of the HTT gene, whose mutation leads to Huntington’s disease (HD), a hereditary degenerative disorder. Specifically, we validated the interactions of canonical miRNAs, miR-137, miR-214, and miR-148a, with the HTT 3′UTR using a luciferase assay. Moreover, we applied synthetic miRNA mimics to examine whether a slight shifting of miRNA seed regions might alter the regulation of the HTT transcript. We also examined miR-137, miR-214, and miR-148a isomiRs and showed the activity of these isoforms on reporter constructs bearing appropriate sequences from the HTT 3′UTR. Hence, we demonstrated that certain 5′-end variants of miRNAs might be functional for the regulation of the same targets as canonical miRNAs.

Distinctive profile of IsomiR expression and novel microRNAs in rat heart left ventricle

MicroRNAs (miRNAs) are single-stranded non-coding RNAs that negatively regulate target gene expression through mRNA cleavage or translational repression. There is mounting evidence that they play critical roles in heart disease. The expression of known miRNAs in the heart has been studied at length by microarray and quantitative PCR but it is becoming evident that microRNA isoforms (isomiRs) are potentially physiologically important. It is well known that left ventricular (patho)physiology is influenced by transmural heterogeneity of cardiomyocyte phenotype, and this likely reflects underlying heterogeneity of gene expression. Given the significant role of miRNAs in regulating gene expression, knowledge of how the miRNA profile varies across the ventricular wall will be crucial to better understand the mechanisms governing transmural physiological heterogeneity. To determinine miRNA/isomiR expression profiles in the rat heart we investigated tissue from different locations across the left ventricular wall using deep sequencing. We detected significant quantities of 145 known rat miRNAs and 68 potential novel orthologs of known miRNAs, in mature, mature* and isomiR formation. Many isomiRs were detected at a higher frequency than their canonical sequence in miRBase and have different predicted targets. The most common miR-133a isomiR was more effective at targeting a construct containing a sequence from the gelsolin gene than was canonical miR-133a, as determined by dual-fluorescence assay. We identified a novel rat miR-1 homolog from a second miR-1 gene; and a novel rat miRNA similar to miR-676. We also cloned and sequenced the rat miR-486 gene which is not in miRBase (v18). Signalling pathways predicted to be targeted by the most highly detected miRNAs include Ubiquitin-mediated Proteolysis, Mitogen-Activated Protein Kinase, Regulation of Actin Cytoskeleton, Wnt signalling, Calcium Signalling, Gap junctions and Arrhythmogenic Right Ventricular Cardiomyopathy. Most miRNAs are not expressed in a gradient across the ventricular wall, with exceptions including miR-10b, miR-21, miR-99b and miR-486.