Divergent target recognition by coexpressed 5′-isomiRs of miR-142-3p and selective viral mimicry

Sequence heterogeneity at the ends of mature microRNAs (miRNAs) is well documented, but its effects on miRNA function are largely unexplored. Here we studied the impact of miRNA 5′-heterogeneity, which affects the seed region critical for target recognition. Using the example of miR-142-3p, an emerging regulator of the hematopoietic lineage in vertebrates, we show that naturally coexpressed 5′-variants (5′-isomiRs) can recognize largely distinct sets of binding sites. Despite this, both miR-142-3p isomiRs regulate exclusive and shared targets involved in actin dynamics. Thus, 5′-heterogeneity can substantially broaden and enhance regulation of one pathway. Other 5′-isomiRs, in contrast, recognize largely overlapping sets of binding sites. This is exemplified by two herpesviral 5′-isomiRs that selectively mimic one of the miR-142-3p 5′-isomiRs. We hypothesize that other cellular and viral 5′-isomiRs can similarly be grouped into those with divergent or convergent target repertoires, based on 5′-sequence features. Taken together, our results provide a detailed characterization of target recognition by miR-142-3p and its 5′-isomiR-specific viral mimic. We furthermore demonstrate that miRNA 5′-end variation leads to differential targeting and can thus broaden the target range of miRNAs.     

A Novel miR-451a isomiR,Associated with Amelanotypic Phenotype,Acts as a Tumor Suppressor in Melanoma by Retarding Cell Migration and Invasion

miRNAs are key regulatory small non-coding RNAs involved in critical steps of melanoma tumorigenesis; however, the relationship between sequence specific variations at the 5′ or 3′ termini (isomiR) of a miRNA and cancer phenotype remains unclear. Deep-sequencing and qRT-PCR showed reduced expression of miR-144/451a cluster and most abundant isomiR (miR451a.1) in dysplastic nevi, in-situ and invasive melanomas compared to common nevi and normal skin (n = 101). miRNA in situ hybridization reproducibly confirmed lost miR-451a.1 in melanoma compared to nevus cells or adjacent keratinocytes. Significantly higher expression of miR-451a.1 was associated with amelanotic phenotype in melanomas (n = 47). In contrast, miR-451a was associated with melanotic phenotype, absent pagetoid scatter of intraepidermal melanocytes, superficial spreading histological subtype and tumor inflammation. Sequencing miRNAs from cultured melanocytes with cytoplasmic melanin gradient (light, medium to dark) showed absent miR-451a while revealing other melanin-associated miRNAs, e.g. miR-30b, miR-100 and miR-590 in darkly and let-7a, let-7i and let-7f in lightly to moderately pigmented cultured melanocytes. Ectopic expression of miR-144/451a in melanoma cell lines resulted in markedly higher levels of mature miR-451a.1 than miR451a or miR-144; and significantly retarded cell migration and inhibited invasion in a glucose-sensitive manner. Surprisingly, these effects were not mediated by calcium binding protein 39 (CAB39), a proven miR451a gene target. miR-144/miR-451a cluster is a novel miRNA locus with tumor suppressive activity in melanoma.     

miRNAome analysis of the mammalian neuronal nicotinic acetylcholine receptor gene family

Nicotine binds to and activates a family of ligand-gated ion channels, neuronal nicotinic acetylcholine receptors (nAChRs). Chronic nicotine exposure alters the expression of various nAChR subtypes, which likely contributes to nicotine dependence; however, the underlying mechanisms regulating these changes remain unclear. A growing body of evidence indicates that microRNAs (miRNAs) may be involved in nAChR regulation. Using bioinformatics, miRNA library screening, site-directed mutagenesis, and gene expression analysis, we have identified a limited number of miRNAs that functionally interact with the 3′-untranslated regions (3′ UTRs) of mammalian neuronal nAChR subunit genes. In silico analyses revealed specific, evolutionarily conserved sites within the 3′ UTRs through which the miRNAs regulate gene expression. Mutating these sites disrupted miRNA regulation confirming the in silico predictions. In addition, the miRNAs that target nAChR 3′ UTRs are expressed in mouse brain and are regulated by chronic nicotine exposure. Furthermore, we show that expression of one of these miRNAs, miR-542-3p, is modulated by nicotine within the mesocorticolimbic reward pathway. Importantly, overexpression of miR-542-3p led to a decrease in the protein levels of its target, the nAChR β2 subunit. Bioinformatic analysis suggests that a number of the miRNAs play a general role in regulating cholinergic signaling. Our results provide evidence for a novel mode of nicotine-mediated regulation of the mammalian nAChR gene family.     

An Exploration of Evolution,Maturation, Expression and Function Relationships in Mir-23～27～24 Cluster

The study aims to explore the potential relationships of evolution, maturation, expression and function between homologous/clustered miRNAs. mir-23∼27∼24 gene cluster, including the two gene clusters (mir-23a and mir-23b) and the three miRNA gene families (mir-23, mir-27 and mir-24), was typically selected as an example. These related miRNAs show similar evolutionary patterns and various expression patterns. Most of them show consistent isomiR expression pattern, and the “switching” phenomenon can be found between different abundant isomiR species. These findings suggest that these sequence or location related miRNAs show the similar miRNA processing and maturation processes, and the robust selection of the most dominant isomiR exists in specific tissues. Functional analysis show that these miRNAs show similar distributions of enriched gene categories, suggesting the close functional prelateships via direct or indirect coordinate regulation in biological processes. The study reveals the close evolutionary, expression and functional relationships between related homologous/clustered miRNAs, which will further enrich miRNA studies and understand direct or indirect interactions between miRNAs.     

miR-1207-5p and miR-1266 suppress gastric cancer growth and invasion by targeting telomerase reverse transcriptase

hTERT is the catalytic subunit of the telomerase complex. Elevated expression of hTERT is associated with the expansion and metastasis of gastric tumor. In this study, we aimed to identify novel tumor suppressor miRNAs that restrain hTERT expression. We began our screen for hTERT-targeting miRNAs with a miRNA microarray. miRNA candidates were further filtered by bioinformatic analysis, general expression pattern in different cell lines, gain-of-function effects on hTERT protein and the potential of these effects to suppress hTERT 3′ untranslated region (3′UTR) luciferase activity. The clinical relevance of two miRNAs (miR-1207-5p and miR-1266) was evaluated by real-time RT-PCR. The effects of these miRNAs on cell growth, cell cycle and invasion of gastric cancer cells were measured with CCK-8, flow cytometry and transwell assays. Finally, the ability of these miRNAs to suppress the transplanted tumors was also investigated. Fourteen miRNAs were identified using a combination of bioinformatics and miRNA microarray analysis. Of these fourteen miRNAs, nine were expressed at significantly lower levels in hTERT-positive cell lines compared with hTERT-negative cell lines and five could downregulate hTERT protein expression. Only miR-1207-5p and miR-1266 interacted with the 3′ UTR of hTERT and the expression levels of these two miRNAs were significantly decreased in gastric cancer tissues. These two miRNAs also inhibited gastric tumor growth in vitro and in vivo. Altogether, miR-1207-5p and miR-1266 were determined to be hTERT suppressors in gastric cancer, and the delivery of these two miRNAs represents a novel therapeutic strategy for gastric cancer treatment.     

3′-Biotin-tagged microRNA-27 does not associate with Argonaute proteins in cells

Synthetic 3′-biotin-tagged microRNAs (miRNAs) have often been used to select interacting messenger RNA (mRNA) and noncoding RNA (ncRNA) targets. Here, we examined the extent of association of 3′-end biotinylated miR-27 with Argonaute (Ago) proteins in transfected human cells using a coimmunoprecipitation assay followed by Northern blot analysis. We report that biotinylated miR-27 does not efficiently associate with Ago compared to unmodified miR-27. These results suggest that 3′-end biotin-modified miRNAs are questionable monitors of miRNA function in cells.     

ADAR2 induces reproducible changes in sequence and abundance of mature microRNAs in the mouse brain

Adenosine deaminases that act on RNA (ADARs) deaminate adenosines to inosines in double-stranded RNAs including miRNA precursors. A to I editing is widespread and required for normal life. By comparing deep sequencing data of brain miRNAs from wild-type and ADAR2 deficient mouse strains, we detect editing sites and altered miRNA processing at high sensitivity. We detect 48 novel editing events in miRNAs. Some editing events reach frequencies of up to 80%. About half of all editing events depend on ADAR2 while some miRNAs are preferentially edited by ADAR1. Sixty-four percent of all editing events are located within the seed region of mature miRNAs. For the highly edited miR-3099, we experimentally prove retargeting of the edited miRNA to novel 3′ UTRs. We show further that an abundant editing event in miR-497 promotes processing by Drosha of the corresponding pri-miRNA. We also detect reproducible changes in the abundance of specific miRNAs in ADAR2-deficient mice that occur independent of adjacent A to I editing events. This indicates that ADAR2 binding but not editing of miRNA precursors may influence their processing. Correlating with changes in miRNA abundance we find misregulation of putative targets of these miRNAs in the presence or absence of ADAR2.     

Selected isomiR expression profiles via arm switching?

A mature miRNA may be generated from 5p or 3p arm of a hairpin precursor. The selection may be flexible via “arm switching”. However, accumulating evidences suggest that both arms of many pre-miRNAs can yield mature functional miRNAs. Herein, we attempted to compare the isomiR expression profiles between the two arms through analyzing in-house and published small RNA deep sequencing datasets. Although many miR-#-5p and miR-#-3p have been reported as functional miRNAs, fewer miRNA pairs (11 and 6 pairs are collected in tumor and normal cells, respectively) are simultaneously identified as abundant miRNA species. According to isomiR types and dominant isomiR species, miR-#-5p and miR-#-3p show various isomiR expression profiles as well as diverse enrichment levels. IsomiR profiles of non-dominant arm are not well-conserved in 5′ ends as well as isomiR profiles of dominant arm. If both the miR-#-5p and miR-#-3p are abundantly expressed, their isomiR expression profiles are always stable across different samples. Similar to diverse enrichment levels of miR-#-5p and miR-#-3p, the isomiR expression patterns may also be influenced by the phenomenon of “arm switching”. The diverged isomiR expression profiles further enrich the complexity of multiple isomiRs, and complicate the coding-non-coding RNA regulatory network.     

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.     

A meta-analysis revealed insights into the sources,conservation and impact of microRNA 5′-isoforms in four model species

MicroRNA (miRNA) 5′-isoforms, or 5′-isomiRs, are small-RNA species that originate from the same genomic loci as the major miRNAs with their 5′ ends shifted from the 5′ ends of the miRNAs by a few nucleotides. Although 5′-isomiRs have been reported, their origins, properties and potential functions remain to be examined. We systematically studied 5′-isomiRs in human, mouse, fruitfly and worm by analysing a large collection of small non-coding RNA and mRNA profiling data. The results revealed a broad existence of 5′-isomiRs in the four species, many of which were conserved and could arise from genomic loci of canonical and non-canonical miRNAs. The well-conserved 5′-isomiRs have several features, including a preference of the 3p over the 5p arms of hairpins of conserved mammalian miRNAs, altered 5′-isomiRs across species and across tissues, and association with structural variations of miRNA hairpins. Importantly, 5′-isomiRs and their major miRNAs may have different mRNA targets and thus potentially play distinct roles of gene regulation, as shown by an integrative analysis combining miRNA and mRNA profiling data from psoriatic and normal human skin and from murine miRNA knockout assays. Indeed, 18 5′-isomiRs had aberrant expression in psoriatic human skin, suggesting their potential function in psoriasis pathogenesis. The results of the current study deepened our understanding of the diversity and conservation of miRNAs, their plasticity in gene regulation and potential broad function in complex diseases.     

Human miR-3145 inhibits influenza A viruses replication by targeting and silencing viral PB1 gene

MicroRNAs (miRNAs) play an important role in the regulation of gene expression and are involved in many cellular processes including inhibition of viral replication in infected cells. In this study, three subtypes of influenza A viruses (pH1N1, H5N1 and H3N2) were analyzed to identify candidate human miRNAs targeting and silencing viral genes expression. Candidate human miRNAs were predicted by miRBase and RNAhybrid based on minimum free energy (MFE) and hybridization patterns between human miRNAs and viral target genes. In silico analysis presented 76 miRNAs targeting influenza A viruses, including 70 miRNAs that targeted specific subtypes (21 for pH1N1, 27 for H5N1 and 22 for H3N2) and 6 miRNAs (miR-216b, miR-3145, miR-3682, miR-4513, miR-4753 and miR-5693) that targeted multiple subtypes of influenza A viruses. Interestingly, miR-3145 is the only candidate miRNA targeting all three subtypes of influenza A viruses. The miR-3145 targets to PB1 encoding polymerase basic protein 1, which is the main component of the viral polymerase complex. The silencing effect of miR-3145 was validated by 3′-UTR reporter assay and inhibition of influenza viral replication in A549 cells. In 3′-UTR reporter assay, results revealed that miR-3145 triggered significant reduction of the luciferase activity. Moreover, expression of viral PB1 genes was also inhibited considerably (P value < 0.05) in viral infected cells expressing mimic miR-3145. In conclusion, this study demonstrated that human miR-3145 triggered silencing of viral PB1 genes and lead to inhibition of multiple subtypes of influenza viral replication. Therefore, hsa-miR-3145 might be useful for alternative treatment of influenza A viruses in the future.     

Genetic Analysis of Loop Sequences in the Let-7 Gene Family Reveal a Relationship between Loop Evolution and Multiple IsomiRs

While mature miRNAs have been widely studied, the terminal loop sequences are rarely examined despite regulating both primary and mature miRNA functions. Herein, we attempted to understand the evolutionary pattern of loop sequences by analyzing loops in the let-7 gene family. Compared to the stable miRNA length distributions seen in most metazoans, higher metazoan species exhibit a longer length distribution. Examination of these loop sequence length distributions, in addition to phylogenetic tree construction, implicated loop sequences as the main evolutionary drivers in miRNA genes. Moreover, loops from relevant clustered miRNA gene families showed varying length distributions and higher levels of nucleotide divergence, even between homologous pre-miRNA loops. Furthermore, we found that specific nucleotides were dominantly distributed in the 5′ and 3′ terminal loop ends, which may contribute to the relatively precise cleavage that leads to a stable isomiR expression profile. Overall, this study provides further insight into miRNA processing and maturation and further enriches our understanding of miRNA biogenesis.     

Unraveling the role of microRNA/isomiR network in multiple primary melanoma pathogenesis

Malignant cutaneous melanoma (CM) is a potentially lethal form of skin cancer whose worldwide incidence has been constantly increasing over the past decades. During their lifetime, about 8% of CM patients will develop multiple primary melanomas (MPMs), usually at a young age and within 3 years from the first tumor/diagnosis. With the aim of improving our knowledge on MPM biology and pathogenesis, we explored the miRNome of 24 single and multiple primary melanomas, including multiple tumors from the same patient, using a small RNA-sequencing approach. From a supervised analysis, 22 miRNAs were differentially expressed in MPM compared to single CM, including key miRNAs involved in epithelial–mesenchymal transition. The first and second melanoma from the same patient presented a different miRNA profile. Ten miRNAs, including miR-25-3p, 149-5p, 92b-3p, 211-5p, 125a-5p, 125b-5p, 205-5p, 200b-3p, 21-5p, and 146a-5p, were further validated in 47 single and multiple melanoma samples. Pathway enrichment analysis of miRNA target genes revealed a more differentiated and less invasive status of MPMs compared to CMs. Bioinformatic analyses at the miRNA isoform (isomiR) level detected a panel of highly expressed isomiRs belonging to miRNA families implicated in human tumorigenesis, including miR-200, miR-30, and miR-10 family. Moreover, we identified hsa-miR-125a-5p|0|−2 isoform as tenfold over-represented in melanoma than the canonical form and differentially expressed in MPMs arising in the same patient. Target prediction analysis revealed that the miRNA shortening could change the pattern of target gene regulation, specifically in genes implicated in cell adhesion and neuronal differentiation. Overall, we provided a putative and comprehensive characterization of the miRNA/isomiR regulatory network of MPMs, highlighting mechanisms of tumor development and molecular features differentiating this subtype from single melanomas.Subject terms: Small RNAs, Melanoma