Functional and evolutionary significance of human microRNA seed region mutations

MicroRNA-125a-3p, -4530, and -92a as a Potential Circulating MicroRNA Panel for Noninvasive Pancreatic Cancer Diagnosis

MicroRNA (miRNA) expression dysregulations in pancreatic ductal adenocarcinoma (PDAC) have been studied widely for their diagnostic and prognostic utility. By the use of bioinformatics-based methods, in our previous study, we identified some potential miRNA panels for diagnosis of pancreatic cancer patients from noncancerous controls (the screening stage). In this report, we used 142 plasma samples from people with and without pancreatic cancer (PC) to conduct RT-qPCR differential expression analysis to assess the strength of the first previously proposed diagnostic panel (consisting of miR-125a-3p, miR-4530, and miR-92a-2-5p). As the result, we identified significant upregulation for all the three considered miRNAs in the serum of PC patients. After that, a three-miRNA panel in serum was developed. The area under the receiver operating characteristic curves (AUC) for the panel were 0.850, 0.910, and 0.86, respectively, indicating that it had a higher diagnostic value than individual miRNAs. Therefore, we detected a promising three-miRNA panel in the plasma for noninvasive PC diagnosis (miR-125a-3p, miR-4530, and miR-92a-2-5p).

Signatures of genetic variation in human microRNAs point to processes of positive selection and population-specific disease risks

The occurrence of natural variation in human microRNAs has been the focus of numerous studies during the last 20 years. Most of them have been focused on the role of specific mutations in disease, while a minor proportion seek to analyse microRNA diversity in the genomes of human populations. We analyse the latest human microRNA annotations in the light of the most updated catalogue of genetic variation provided by the 1000 Genomes Project. By means of the in silico analysis of microRNA genetic variation we show that the level of evolutionary constraint of these sequences is governed by the interplay of different factors, like their evolutionary age or genomic location. The role of mutations in the shaping of microRNA-driven regulatory interactions is emphasized with the acknowledgement that, while the whole microRNA sequence is highly conserved, the seed region shows a pattern of higher genetic diversity that appears to be caused by the dramatic frequency shifts of a fraction of human microRNAs. We highlight the participation of these microRNAs in population-specific processes by identifying that not only the seed, but also the loop, are particularly differentiated regions among human populations. The quantitative computational comparison of signatures of population differentiation showed that candidate microRNAs with the largest differences are enriched in variants implicated in gene expression levels (eQTLs), selective sweeps and pathological processes. We explore the implication of these evolutionary-driven microRNAs and their SNPs in human diseases, such as different types of cancer, and discuss their role in population-specific disease risk.

A novel microRNA, novel-m009C, regulates methamphetamine rewarding effects

Methamphetamine (METH) is a widely abused psychostimulant, whose hyper-rewarding property is believed to underlie its addictive effect, but the molecular mechanism regulating this effect remains unclear. We previously reported that decreased expression of a novel microRNA (miRNA), novel-m009C, is implicated in the regulation of METH hyperlocomotion. Here, we found that novel-m009C may be homologous to hsa-miR-604. Its expression is consistently downregulated in the nucleus accumbens (NAc) of mice when exposed to METH and cocaine, whereas significant alterations in novel-m009C expression were not observed in the NAc of mice subjected to other rewarding and psychiatric stimuli, such as sucrose, morphine and MK-801. We further found the substantial reduction in novel-m009C expression may be regulated by both dopamine receptor D1 (D1R) and D2 (D2R). Increasing novel-m009C levels in the NAc attenuated METH-induced conditioned place preference (CPP) and hyperlocomotion, whereas inhibiting novel-m009C expression in the NAc enhanced these effects but did not change the preference of mice for a natural reward, i.e., sucrose. These effects may involve targeting of genes important for the synaptic transmission, such as Grin1 (NMDAR subunit 1). Our findings demonstrate an important role for NAc novel-m009C in regulating METH reward, reveal a novel molecular regulator of the actions of METH on brain reward circuitries and provide a new strategy for treating METH addiction based on the modulation of small non-coding RNAs.

Circ_0009035 regulates the progression of cervical cancer by targeting miR-1305/CREBRF axis

Circular RNAs (circRNAs) have a crucial role in the occurrence of many diseases, such as tumors. Yet the roles of circ_0009035 (circRACGAP1) in cervical cancer are not fully characterized. The expression levels of circRACGAP1, miR-1305 and cAMP-responsive element-binding protein 3 regulatory factor (CREBRF) were detected by using real-time quantitative PCR or western blot. Cell counting kit-8 assay, 5-ethynyl-2'-deoxyuridine, colony formation assay, transwell assay and tube formation assay were used to detect cell proliferation, migration and invasion and angiogenesis, respectively. Flow cytometry assay was used to analyze the cell apoptosis. Dual-luciferase reporter assay and RNA immunoprecipitation assay were performed to analyze the targeting about miR-1305 and circ_0009035 or CREBRF. Xenograft model was built to study the role of circ_0009035 in vivo. Immunohistochemistry was used to detect the expression of Ki67, epithelial cadherin and vimentin. First, we found that circ_0009035 expression was significantly upregulated in tumor cells and tissues; second, knockdown of circ_0009035 could inhibit cell proliferation, migration and invasion and promote cell apoptosis. Subsequently, circ_0009035 was found to be able to target miR-1305, and the expression of miR-1305 in tumor tissues and cells was significantly lower. MiR-1305 inhibitor could restore cell-related progression of cervical cancer inhibited by si-circ_0009035. Finally, miR-1305 could target CREBRF, and circ_0009035 could regulate CREBRF expression by targeting miR-1305, thereby affecting cervical cancer tumorigenesis. In summary, our study confirmed that circ_0009035 could influence the development of cervical cancer through the targeted regulation of miR-1305/CREBRF.

A Mutation in Endogenous saRNA miR-23a Influences Granulosa Cells Response to Oxidative Stress

Phenotypes are the result of the interaction between the gene and the environment, so the response of individuals with different genotypes to an environment is variable. Here, we reported that a mutation in miR-23a influences granulosa cells (GCs) response to oxidative stress, a common mechanism of environmental factors affecting female reproduction. We showed that nuclear miR-23a is a pro-apoptotic miRNA in porcine GCs through the activation of the transcription and function of NORHA, a long non-coding RNA (lncRNA) induces GC apoptosis and responses to oxidative stress. Mechanistically, miR-23a acts as an endogenous small activating RNA (saRNA) to alter histone modifications of the NORHA promoter through the direct binding to its core promoter. A C > T mutation was identified at −398 nt of the miR-23a core promoter, which created a novel binding site for the transcription factor SMAD4 and recruited the transcription repressor SMAD4 to inhibit miR-23a transcription and function in GCs. Notably, g.−398C > T mutation in the miR-23a promoter reduced GCs response to oxidative stress. In addition, g.−398C > T mutation was significantly associated with sow fertility traits. In short, our findings preliminarily revealed the genetic basis of individual differences in the response to oxidative stress from the perspective of a single mutation and identified miR-23a as a candidate gene for the environmental adaptation to oxidative stress.

Genetic Variants of MIR27A, MIR196A2 May Impact the Risk for the Onset of Coronary Artery Disease in the Pakistani Population

Genetic variants in microRNA genes have a detrimental effect on miRNA-mediated regulation of gene expression and may contribute to coronary artery disease (CAD). CAD is the primary cause of mortality worldwide. Several environmental, genetic, and epigenetic factors are responsible for CAD susceptibility. The contribution of protein-coding genes is extensively studied. However, the role of microRNA genes in CAD is at infancy. The study is aimed to investigate the impact of rs895819, rs11614913, and rs2168518 variants in MIR27A, MIR196A2, and MIR4513, respectively, in CAD using allele-specific PCR. Results: For variant rs11614913, significant distribution of the genotypes among the cases and controls was determined by co-dominant [χ² = 54.4; p value ≤ 0.0001], dominant (C/C vs. C/T + T/T) [OR = 0.257 (0.133-0.496); p value ≤ 0.0001], recessive (T/T vs. C/T + C/C) [OR = 1.56 (0.677-0.632); p value = 0.398], and additive models [OR = 0.421 (0.262-0.675); p value = 0.0004]. Similarly, a significant association of rs895819 was determined by co-dominant [χ² = 9.669; p value ≤ 0.008], dominant (A/A vs. A/G + G/G) [OR = 0.285 (0.1242-0.6575); p value ≤ 0.0034], recessive (G/G vs. A/G + A/A) [OR = 0.900 (0.3202-3.519); p value = 1.000], and additive models [OR = 0.604 (0.3640-1.002); p value = 0.05] while no significant association of rs2168518 with CAD was found. Conclusion: The variants rs895819 and rs11614913 are the susceptibility factors for CAD.

In silico Analysis of Single Nucleotide Polymorphisms Associated with MicroRNA Regulating 5-fluorouracil Resistance in Colorectal Cancer

Background:

Due to the broad influence and reversible nature of microRNA (miRNA) on the expression and regulation of target genes, researchers suggest that miRNAs and single nucleotide polymorphisms (SNPs) in miRNA genes interfere with 5-fluorouracil (5-FU) drug resistance in colorectal cancer chemotherapy.

Methods:

Computational assessment and cataloging of miRNA gene polymorphisms that target mRNA transcripts directly or indirectly through regulation of 5-FU chemoresistance in CRC were screened out by applying various universally accessible datasets such as miRNA SNP3.0 software.

Results:

1255 SNPs in 85 miRNAs affecting 5-FU resistance (retrieved from literature) were detected. Computational analysis showed that 167 from 1255 SNPs alter microRNA expression levels leading to inadequate response to 5-FU resistance in CRC. Among these 167 SNPs, 39 were located in the seed region of 25/85 miRNA and were more critical than other SNPs. Has-miR-320a-5p with 4 SNP in seed region was miRNA with the most number of SNPs. On the other hand, it has been identified that proteoglycan in cancer, adherents junction, ECM-receptor interaction, Hippo signaling pathway, TGF-beta signaling cascade, biosynthesis of fatty acid, and fatty acid metabolism were the most important pathways targeted by these 85 predicted miRNAs.

Conclusion:

Our data suggest 39 SNPs in the seed region of 25 miRNAs as catalog in miRNA genes that control the 5-FU resistance in CRC. These data also identify the most important pathways regulated by miRNA.

MicroRNA-210 regulates the metabolic and inflammatory status of primary human astrocytes

Background

Astrocytes are the most numerous glial cell type with important roles in maintaining homeostasis and responding to diseases in the brain. Astrocyte function is subject to modulation by microRNAs (miRs), which are short nucleotide strands that regulate protein expression in a post-transcriptional manner. Understanding the miR expression profile of astrocytes in disease settings provides insight into the cellular stresses present in the microenvironment and may uncover pathways of therapeutic interest.

Methods

Laser-capture microdissection was used to isolate human astrocytes surrounding stroke lesions and those from neurological control tissue. Astrocytic miR expression profiles were examined using quantitative reverse transcription polymerase chain reaction (RT-qPCR). Primary human fetal astrocytes were cultured under in vitro stress conditions and transfection of a miR mimic was used to better understand how altered levels of miR-210 affect astrocyte function. The astrocytic response to stress was studied using qPCR, enzyme-linked immunosorbent assays (ELISAs), measurement of released lactate, and Seahorse.

Results

Here, we measured miR expression levels in astrocytes around human ischemic stroke lesions and observed differential expression of miR-210 in chronic stroke astrocytes compared to astrocytes from neurological control tissue. We also identified increased expression of miR-210 in mouse white matter tissue around middle cerebral artery occlusion (MCAO) brain lesions. We aimed to understand the role of miR-210 in primary human fetal astrocytes by developing an in vitro assay of hypoxic, metabolic, and inflammatory stresses. A combination of hypoxic and inflammatory stresses was observed to upregulate miR-210 expression. Transfection with miR-210-mimic (210M) increased glycolysis, enhanced lactate export, and promoted an anti-inflammatory transcriptional and translational signature in astrocytes. Additionally, 210M transfection resulted in decreased expression of complement 3 (C3) and semaphorin 5b (Sema5b).

Conclusions

We conclude that miR-210 expression in human astrocytes is modulated in response to ischemic stroke disease and under in vitro stress conditions, supporting a role for miR-210 in the astrocytic response to disease conditions. Further, the anti-inflammatory and pro-glycolytic impact of miR-210 on astrocytes makes it a potential candidate for further research as a neuroprotective agent.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12974-021-02373-y.

In silico identification of variations in microRNAs with a potential impact on dairy traits using whole ruminant genome SNP datasets

MicroRNAs are small noncoding RNAs that have important roles in the lactation process and milk biosynthesis. Some polymorphisms have been studied in various livestock species from the perspective of pathology or production traits. To target variants that could be the causal variants of dairy traits, genetic variants of microRNAs expressed in the mammary gland or present in milk and localized in dairy quantitative trait loci (QTLs) were investigated in bovine, caprine, and ovine species. In this study, a total of 59,124 (out of 28 millions), 13,427 (out of 87 millions), and 4761 (out of 38 millions) genetic variants in microRNAs expressed in the mammary gland or present in milk were identified in bovine, caprine, and ovine species, respectively. A total of 4679 of these detected bovine genetic variants are located in dairy QTLs. In caprine species, 127 genetic variants are localized in dairy QTLs. In ovine species, no genetic variant was identified in dairy QTLs. This study leads to the detection of microRNA genetic variants of interest in the context of dairy production, taking advantage of whole genome data to identify microRNA genetic variants expressed in the mammary gland and localized in dairy QTLs.

Osteoclast signaling-targeting miR-146a-3p and miR-155-5p are downregulated in Paget's disease of bone

MicroRNA (miRs) are small, non-coding RNA that post-transcriptionally regulate DNA expression. We hypothesized that specific miR profiles may be a feature of overactive osteoclasts in Paget's disease of bone (PDB), a disorder characterized by an increased and disorganized bone remodeling that typically begins with excessive bone resorption. We compared the expression profile of 13 miRs in human osteoclasts differentiated in vitro from peripheral blood mononuclear cells (PBMCs) of patients with PDB (n = 10) or age- and sex- matched healthy subjects (n = 10). We selected 13 miRs for testing, on the basis of their previously reported roles either in human osteoclast differentiation, in bone diseases, or in osteoclast important signaling pathways. From those expression results, 3 miRNAs were further selected for in-vitro studies aiming at modulating miR expression in human cord blood monocyte derived osteoclasts: 2 miRs (miR-146a-3p and miR-155-5p) whose expression was significantly reduced in pagetic osteoclasts, as well as miRNA-133a-3p, stable in PDB relative to controls, but with known regulatory importance within osteoclasts. We demonstrated a positive (miR-133a-3p) or negative (miR-155-5p, miR-146a-3p) impact of those miRs on the formation of osteoclasts and/or their bone resorption capacity in this human model. Signaling pathways were significantly affected, including p38 MAP-kinase (miR-133a-3p), RANKL-induced TRAF6/NFκB signaling (miR-146a-3p), and MITF expression (miR-155-5p). Osteoclast miRNA profiles might have an important value to yield significant new insights into the osteoclast phenotype in PDB and in other bone diseases with hyperactive osteoclasts.

The ability of miRNAs to induce mesenchymal-to-epithelial transition (MET) in cancer cells is highly dependent upon genetic background

Understanding of the molecular basis of host cell-miRNA interactions is prerequisite to the successful application of miRNAs as potential therapeutic agents. We studied the morphological and molecular consequences of over expression of three sequence divergent miRNAs previously implicated in the mesenchymal-to-epithelial transition process (MET) in three distinct mesenchymal-like cancer cell lines. The ability of miRNAs to induce morphological changes characteristic of MET positively correlated with induced changes in the expression of genes previously implicated in the process. Variability in the responses of different mesenchymal-like cells to over expression of the same miRNAs was attributable to inherent differences in trans-regulatory profiles pre-disposing these cells to miRNA-induced MET. Collectively our results indicate that miRNA-mediated regulation of MET is a highly integrated process that is significantly modulated by the molecular background of individual cells.

MicroRNAs in Pancreatic Cancer: biomarkers, prognostic, and therapeutic modulators

A severe lack of early diagnosis coupled with resistance to most available therapeutic options renders pancreatic cancer as a major clinical concern. The limited efficacy of current treatments necessitates the development of novel therapeutic strategies that are based on an understanding of the molecular mechanisms involved in pancreatic cancer progression. MicroRNAs (miRNAs) are non-coding small RNAs that regulate the expression of multiple proteins in the post-translation process and thus have promise as biomarkers, prognostic agents, and as advanced pancreatic therapies. Profiling of deregulated miRNAs in pancreatic cancer can correlate to diagnosis, indicate optimal treatment and predict response to therapy. Furthermore, understanding the main effector genes in pancreatic cancer along with downstream pathways can identify possible miRNAs as therapeutic candidates. Additionally, obstacles to the translation of miRNAs into the clinic are also considered. Distinct miRNA expression profiles can correlate to stages of malignant pancreatic disease, and hold potential as biomarkers, prognostic markers and clinical targets. However, a limited understanding and validation of the specific role of such miRNAs stunts clinical application. Target prediction using algorithms provides a wide range of possible targets, but these miRNAs still require validation through pre-clinical studies to determine the knock-on genetic effects.

Inhibition of miR-378a-3p by Inflammation Enhances IL-33 Levels: A Novel Mechanism of Alarmin Modulation in Ulcerative Colitis

Ulcerative colitis (UC) is an inflammatory bowel disease (IBD) characterized by mucosa damage associated with an uncontrolled inflammatory response. This immunological impairment leads to altered inflammatory mediators such as IL-33, which is shown to increase in the mucosa of active UC (aUC) patients. MicroRNAs present a distorted feature in inflamed colonic mucosa and are potential IL-33 regulating candidates in UC. Therefore, we studied the microRNA and mRNA profiles in inflamed colonic samples of UC patients, evaluating the effect of a microRNA (selected by in silico analysis and its expression in UC patients), on IL-33 under inflammatory conditions. We found that inflamed mucosa (n = 8) showed increased expression of 40 microRNAs and 2,120 mRNAs, while 49 microRNAs and 1,734 mRNAs were decreased, as determined by microarrays. In particular, IL-33 mRNA showed a 3.8-fold increase and eight members of a microRNA family (miR-378), which targets IL-33 mRNA in the 3'UTR, were decreased (-3.9 to -3.0 times). We selected three members of the miR-378 family (miR-378a-3p, miR-422a, and miR-378c) according to background information and interaction energy analysis, for further correlation analyses with IL-33 expression through qPCR and ELISA, respectively. We determined that aUC (n = 24) showed high IL-33 levels, and decreased expression of miR-378a-3p and miR-422a compared to inactive UC (n = 10) and controls (n = 6). Moreover, both microRNAs were inversely correlated with IL-33 expression, while miR-378c does not show a significant difference. To evaluate the effect of TNFα on the studied microRNAs, aUC patients with anti-TNF therapy were compared to aUC receiving other treatments. The levels of miR-378a-3p and miR-378c were higher in aUC patients with anti-TNF. Based on these findings, we selected miR-378a-3p to exploring the molecular mechanism involved by in vitro assays, showing that over-expression of miR-378a-3p decreased the levels of an IL-33 target sequence β-gal-reporter gene in HEK293 cells. Stable miR-378a-3p over-expression/inhibition inversely modulated IL-33 content and altered viability of HT-29 cells. Additionally, in an inflammatory context, TNFα decreased miR-378a-3p levels in HT-29 cells enhancing IL-33 expression. Together, our results propose a regulatory mechanism of IL-33 expression exerted by miR-378a-3p in an inflammatory environment, contributing to the understanding of UC pathogenesis.

Functional Analysis of Genetic Variants and Somatic Mutations Impacting MicroRNA-Target Recognition: Bioinformatics Resources

MicroRNAs are small noncoding RNA molecules with great importance in regulating a large number of diverse biological processes in health and disease. MicroRNAs can bind to both coding and noncoding RNAs and regulate their stability and expression. Genetic variants and somatic mutations may alter microRNA sequences and their target sites and therefore impact microRNA-target recognition. Aberrant microRNA-target interactions have been associated with many diseases. In recent years, computational resources have been developed for retrieving, annotating, and analyzing the impact of mutations on microRNA-target recognition. In this chapter, we provide an overview on the computational analysis of mutations impacting microRNA target recognition, followed by a detailed tutorial on how to use three major Web-based bioinformatics resources: PolymiRTS ( http://compbio.uthsc.edu/miRSNP ), a database of genetic variants impacting microRNA target recognition; SomamiR ( http://compbio.uthsc.edu/SomamiR ), a database of somatic mutations affecting the interactions between microRNAs and their targets in mRNAs and noncoding RNAs; and miR2GO ( http://compbio.uthsc.edu/miR2GO ), a computational tool for knowledge-based functional analysis of genetic variants and somatic mutations in microRNA seed regions.

Framework for microRNA variant annotation and prioritization using human population and disease datasets

MicroRNA (miRNA) expression is frequently deregulated in human disease, in contrast, disease-associated miRNA mutations are understudied. We developed Annotative Database of miRNA Elements, ADmiRE, which combines multiple existing and new biological annotations to aid prioritization of causal miRNA variation. We annotated 10,206 mature (3,257 within seed region) miRNA variants from multiple large sequencing datasets including gnomAD (15,496 genomes; 123,136 exomes). The pattern of miRNA variation closely resembles protein-coding exonic regions, with no difference between intragenic and intergenic miRNAs (P = 0.56), and high confidence miRNAs demonstrate higher sequence constraint (P < 0.001). Conservation analysis across 100 vertebrates identified 765 highly conserved miRNAs that also have limited genetic variation in gnomAD. We applied ADmiRE to the TCGA PanCancerAtlas WES dataset containing over 10,000 individuals across 33 adult cancers and annotated 1,267 germline (rare in gnomAD) and 1,492 somatic miRNA variants. Several miRNA families with deregulated gene expression in cancer have low levels of both somatic and germline variants, e.g., let-7 and miR-10. In addition to known somatic miR-142 mutations in hematologic cancers, we describe novel somatic miR-21 mutations in esophageal cancers impacting downstream miRNA targets. Through the development of ADmiRE, we present a framework for annotation and prioritization of miRNA variation in disease datasets.

Expanding the horizons of microRNA bioinformatics

MicroRNA regulation of key biological and developmental pathways is a rapidly expanding area of research, accompanied by vast amounts of experimental data. This data, however, is not widely available in bioinformatic resources, making it difficult for researchers to find and analyze microRNA-related experimental data and define further research projects. We are addressing this problem by providing two new bioinformatics data sets that contain experimentally verified functional information for mammalian microRNAs involved in cardiovascular-relevant, and other, processes. To date, our resource provides over 4400 Gene Ontology annotations associated with over 500 microRNAs from human, mouse, and rat and over 2400 experimentally validated microRNA:target interactions. We illustrate how this resource can be used to create microRNA-focused interaction networks with a biological context using the known biological role of microRNAs and the mRNAs they regulate, enabling discovery of associations between gene products, biological pathways and, ultimately, diseases. This data will be crucial in advancing the field of microRNA bioinformatics and will establish consistent data sets for reproducible functional analysis of microRNAs across all biological research areas.

Clustering Pattern and Functional Effect of SNPs in Human miRNA Seed Regions

miRNAs are a class of noncoding RNAs important in posttranscriptional repressors and involved in the regulation of almost every biological process by base paring with target genes through sequence in their seed regions. Genetic variations in the seed regions have vital effects on gene expression, phenotypic variation, and disease susceptibility in humans. The distribution pattern of genetic variation in miRNA seed regions might be related to miRNA function and is worth paying more attention to. We here employed computational analyses to explore the clustering pattern and functional effect of SNPs in human miRNA seed regions. A total of 1879 SNPs were mapped to 1226 human miRNA seed regions. We found that miRNAs with SNPs in their seed region are significantly enriched in miRNA clusters. We also found that SNPs in clustered miRNA seed regions have a lower functional effect than have SNPs in nonclustered miRNA seed regions. Additionally, we found that clustered miRNAs with SNPs in seed regions are involved in more pathways. Overall, our results demonstrate that SNPs in clustered miRNA seed regions can take part in more intricate and complex gene-regulating networks with lower functional cost by functional complementarity. Moreover, our results also broaden current knowledge on the genetic variation in human miRNA seed regions.

microRNA single polynucleotide polymorphism influences on microRNA biogenesis and mRNA target specificity

microRNAs (miRNAs) are nowadays recognized as an essential component of gene regulatory networks. Furthermore, deregulation of miRNAs expression often contributes to human pathologies. Recently, a substantial number of single nucleotide polymorphism (SNPs) and rare mutations within pri-, pre- and mature miRNA sequences have been reported. These miRNA SNPs have often been associated with human disease. However, due to the complexity of miRNA biogenesis and the genome-wide functional effects of miRNAs, the determination of biological consequences of these miRNA SNPs remains challenging. Despite an increasing number of reports linking miRNA SNPs with human pathologies, few reports have analyzed the mechanism by which miRNA-SNPs contribute to disease pathogenesis. In this review, we discuss how single polynucleotide polymorphisms in miRNAs genes may influence miRNAs expression and function and thus potentially alter disease pathogenesis.

The Enigma of miRNA Regulation in Cancer

MicroRNAs (miRNAs or miRs) are small 19-22 nucleotide long, noncoding, single-stranded, and multifunctional RNAs that regulate a diverse assortment of gene and protein functions that impact on a vast network of pathways. Lin-4, a noncoding transcript discovered in 1993 and named miRNA, initiated the exploration of research into these intriguing molecules identified in almost all organisms. miRNAs interfere with translation or posttranscriptional regulation of their target gene and regulate multiple biological actions exerted by these target genes. In cancer, they function as both oncogenes and tumor suppressor genes displaying differential activity in various cellular contexts. Although the role of miRNAs on target gene functions has been extensively investigated, less is currently known about the upstream regulatory molecules that regulate miRNAs. This chapter focuses on the factors and processes involved in miRNA regulation.

MicroRNA Genetic Variation: From Population Analysis to Functional Implications of Three Allele Variants Associated with Cancer

Nucleotide variants in microRNA regions have been associated with disease; nevertheless, few studies still have addressed the allele-dependent effect of these changes. We studied microRNA genetic variation in human populations and found that while low-frequency variants accumulate indistinctly in microRNA regions, the mature and seed regions tend to be depleted of high-frequency variants, probably as a result of purifying selection. Comparison of pairwise population fixation indexes among regions showed that the seed had higher population fixation indexes than the other regions, suggesting the existence of local adaptation in the seed region. We further performed functional studies of three microRNA variants associated with cancer (rs2910164:C > G in MIR146A, rs11614913:C > T in MIR196A2, and rs3746444:A > G in both MIR499A and MIR499B). We found differences in the expression between alleles and in the regulation of several genes involved in cancer, such as TP53, KIT, CDH1, CLH, and TERT, which may result in changes in regulatory networks related to tumorigenesis. Furthermore, luciferase-based assays showed that MIR499A could be regulating the cadherin CDH1 and the cell adhesion molecule CLH1 in an allele-dependent fashion. A better understanding of the effect of microRNA variants associated with disease could be key in our way to a more personalized medicine.

Knowledge-based analysis of functional impacts of mutations in microRNA seed regions

MicroRNAs are a class of important post-transcriptional regulators. Genetic and somatic mutations in miRNAs, especially those in the seed regions, have profound and broad impacts on gene expression and physiological and pathological processes. Over 500 SNPs were mapped to the miRNA seeds, which are located at position 2-8 of the mature miRNA sequences. We found that the central positions of the miRNA seeds contain fewer genetic variants and therefore are more evolutionary conserved than the peripheral positions in the seeds. We developed a knowledgebased method to analyse the functional impacts of mutations in miRNA seed regions. We computed the gene ontology-based similarity score GOSS and the GOSS percentile score for all 517 SNPs in miRNA seeds. In addition to the annotation of SNPs for their functional effects, in the present article we also present a detailed analysis pipeline for finding the key functional changes for seed SNPs. We performed a detailed gene ontology graph-based analysis of enriched functional categories for miRNA target gene sets. In the analysis of a SNP in the seed region of hsa-miR-96 we found that two key biological processes for progressive hearing loss 'Neurotrophin TRK receptor signaling pathway' and 'Epidermal growth factor receptor signaling pathway' were significantly and differentially enriched by the two sets of allele-specific target genes of miRNA hsa-miR-96.

LncRNA-MIF, a c-Myc-activated long non-coding RNA, suppresses glycolysis by promoting Fbxw7-mediated c-Myc degradation

The c-Myc proto-oncogene is activated in more than half of all human cancers. However, the precise regulation of c-Myc protein stability is unknown. Here, we show that the lncRNA-MIF (c-Myc inhibitory factor), a c-Myc-induced long non-coding RNA, is a competing endogenous RNA for miR-586 and attenuates the inhibitory effect of miR-586 on Fbxw7, an E3 ligase for c-Myc, leading to increased Fbxw7 expression and subsequent c-Myc degradation. Our data reveal the existence of a feedback loop between c-Myc and lncRNA-MIF, through which c-Myc protein stability is finely controlled. Additionally, we show that the lncRNA-MIF inhibits aerobic glycolysis and tumorigenesis by suppressing c-Myc and miR-586.

microRNA editing in seed region aligns with cellular changes in hypoxic conditions

RNA editing is a finely tuned, dynamic mechanism for post-transcriptional gene regulation that has been thoroughly investigated in the last decade. Nevertheless, RNA editing in non-coding RNA, such as microRNA (miRNA), have caused great debate and have called for deeper investigation. Until recently, in fact, inadequate methodologies and experimental contexts have been unable to provide detailed insights for further elucidation of RNA editing affecting miRNAs, especially in cancer.In this work, we leverage on recent innovative bioinformatics approaches applied to a more informative experimental context in order to analyze the variations in miRNA seed region editing activity during a time course of a hypoxia-exposed breast cancer cell line. By investigating its behavior in a dynamic context, we found that miRNA editing events in the seed region are not depended on miRNA expression, unprecedentedly providing insights on the targetome shifts derived from these modifications. This reveals that miRNA editing acts under the influence of environmentally induced stimuli.Our results show a miRNA editing activity trend aligning with cellular pathways closely associated to hypoxia, such as the VEGF and PI3K/Akt pathways, providing important novel insights on this poorly elucidated phenomenon.

SomamiR 2.0: a database of cancer somatic mutations altering microRNA-ceRNA interactions

SomamiR 2.0 (http://compbio.uthsc.edu/SomamiR) is a database of cancer somatic mutations in microRNAs (miRNA) and their target sites that potentially alter the interactions between miRNAs and competing endogenous RNAs (ceRNA) including mRNAs, circular RNAs (circRNA) and long noncoding RNAs (lncRNA). Here, we describe the recent major updates to the SomamiR database. We expanded the scope of the database by including somatic mutations that impact the interactions between miRNAs and two classes of non-coding RNAs, circRNAs and lncRNAs. Recently, a large number of miRNA target sites have been discovered by newly emerged high-throughput technologies for mapping the miRNA interactome. We have mapped 388 247 somatic mutations to the experimentally identified miRNA target sites. The updated database also includes a list of somatic mutations in the miRNA seed regions, which contain the most important guiding information for miRNA target recognition. A recently developed webserver, miR2GO, was integrated with the database to provide a seamless pipeline for assessing functional impacts of somatic mutations in miRNA seed regions. Data and functions from multiple sources including biological pathways and genome-wide association studies were updated and integrated with SomamiR 2.0 to make it a better platform for functional analysis of somatic mutations altering miRNA–ceRNA interactions.

A Uniform System for the Annotation of Vertebrate microRNA Genes and the Evolution of the Human microRNAome

Although microRNAs (miRNAs) are among the most intensively studied molecules of the past 20 years, determining what is and what is not a miRNA has not been straightforward. Here, we present a uniform system for the annotation and nomenclature of miRNA genes. We show that less than a third of the 1,881 human miRBase entries, and only approximately 16% of the 7,095 metazoan miRBase entries, are robustly supported as miRNA genes. Furthermore, we show that the human repertoire of miRNAs has been shaped by periods of intense miRNA innovation and that mature gene products show a very different tempo and mode of sequence evolution than star products. We establish a new open access database--MirGeneDB ( http://mirgenedb.org )--to catalog this set of miRNAs, which complements the efforts of miRBase but differs from it by annotating the mature versus star products and by imposing an evolutionary hierarchy upon this curated and consistently named repertoire.

Genetic polymorphism of 3' untranslated region of zeta-chain associated protein kinase 70 kDa in southern Taiwanese patients with rheumatoid arthritis

T cell activation participates in the pathogenesis of rheumatoid arthritis (RA), and the signaling molecule zeta-chain-associated protein kinase 70 kDa (ZAP-70) plays a crucial role in this process. Different mutations in the coding sequence of ZAP-70 are involved in a variety of immunological phenotypes, and recent evidence indicates that genetic variations within the 3' untranslated regions (UTR) of microRNA binding sites may affect the hybridization with target mRNAs, leading to phenotype changes with disease status. In this study, we evaluated the possible effect of ZAP-70 polymorphism as a genetic risk factor in RA by examining the single-nucleotide polymorphism in 100 patients and 100 ethnicity- and sex-matched healthy individuals from southern Taiwan. In both groups, the genotype distribution of rs2278699 in the 3' UTR was in the Hardy-Weinberg equilibrium. In RA, there were higher frequencies of the G allele (15.5 versus 8.0 %, odds ratio 2.1, P = 0.020) and significant differences in the trend of various genotypes (P = 0.024). The results suggest that genetic polymorphism in the 3' UTR of ZAP-70 is associated with RA susceptibility in southern Taiwanese.

MiR-200 modulates coelomocytes antibacterial activities and LPS priming via targeting Tollip in Apostichopus japonicus

In order to explore the potential roles of microRNAs (miRNAs) in regulating Toll-like receptor (TLR) pathways, we identified Toll interacting protein as a putative target of miR-200 in Apostichopus japonicus coelomocytes by RNA-seq screening (denoted as AjTollip). The positive expression profiles of miR-200 and AjTollip were detected in both LPS exposure primary coelomocytes and Vibrio splendidus challenge sea cucumber. Co-infection miR-200 mimics significantly elevated the expression of AjTollip and its down-stream molecules. In contrast, miR-200 inhibitor significantly repressed the expression of these TLR-pathway members. More importantly, miR-200 displayed not only to enhance coelomocytes antibacterial activities, but to suppress LPS priming in vitro. Overall, all these results will enhance our understanding on miR-200 regulatory roles in anti-bacterial process in sea cucumber.