Unravelling the multiple roles of FUS in RNA processing and on ALS pathogenesis

Two genes with related functions in RNA biogenesis were recently reported in patients with familial ALS: the FUS/TLS gene at the ALS6 locus and the TARDBP/TDP-43 gene at the ALS10 locus [1, 2]. FUS has been implicated to function in several steps of gene expression, including transcription regulation [3], RNA splicing [4, 5], mRNA transport in neurons [6] and, interestingly, in microRNA (miRNA) processing [7]. The goal of this project is to identify the molecular mechanisms leading to the development of FUS mutations-associated ALS. Specifically, we want to test the hypothesis that these FUS mutations misregulate miRNA levels that in turn affect the expression of genes critical for motor neuron survival. In addition we want to test whether misregulation of the miRNA profile is a common feature in ALS. We have performed immunoprecipitations from total extracts of 293T cells expressing FLAG-tagged FUS to characterize its interactome by mass spectrometry. This proteomic study not only revealed a strong interaction of FUS with splicing factors, but shows that FUS might be involved in many, quite different pathways. To map which parts of the FUS protein contribute to the interaction with splicing factors, we have performed a set of experiments with a series of missense and deletion mutants. With this approach, we will not only gain information on the binding partners of FUS along with a map of the required domains for the interactions, but it will also help to unravel whether certain ALS-associated FUS mutations lead to a loss or gain of function due to gain or loss of interactors. Additionally, we have performed quantitative interactomics using SILAC to identify interactome differences of ALS-associated FUS mutants. To this end we have performed immunoprecipitations of total extract from 293T cells, stably transduced with constructs expressing wild-type FUS-FLAG as well as three different ALS-associated mutants (G156E, R244C, P525L). First results indicate striking differences in the interactome with certain RNA binding proteins. We are now validating these candidates in order to reveal the importance of these differential interactions in the context of ALS.

The H19/let-7 double-negative feedback loop contributes to glucose metabolism in muscle cells

The H19 lncRNA has been implicated in development and growth control and is associated with human genetic disorders and cancer. Acting as a molecular sponge, H19 inhibits microRNA (miRNA) let-7. Here we report that H19 is significantly decreased in muscle of human subjects with type-2 diabetes and insulin resistant rodents. This decrease leads to increased bioavailability of let-7, causing diminished expression of let-7 targets, which is recapitulated in vitro where H19 depletion results in impaired insulin signaling and decreased glucose uptake. Furthermore, acute hyperinsulinemia downregulates H19, a phenomenon that occurs through PI3K/AKT-dependent phosphorylation of the miRNA processing factor KSRP, which promotes biogenesis of let-7 and its mediated H19 destabilization. Our results reveal a previously undescribed double-negative feedback loop between sponge lncRNA and target miRNA that contributes to glucose regulation in muscle cells.

Transcriptional and translational regulators of gene expression in germ cell development

Germ cell development is a highly coordinated process driven, in part, by regulatory mechanisms that control gene expression. Not only transcription, but also translation, is under regulatory control to direct proper germ cell development. In this dissertation, I have focused on two regulators of germ cell development. One is the homeobox protein RHOX10, which has the potential to be both a transcriptional and translational regulator in mouse male germ cell development. The other is the RNA-binding protein, Hermes, which functions as a translational regulator in Xenopus laevis female germ cell development. ^ Rhox10 is a member of reproductive homeobox gene X-(linked (Rhox) gene cluster, of which expression is developmentally regulated in developing mouse testes. To identify the cell types and developmental stages in which Rhox10 might function, I characterized its temporal and spatial expression pattern in mouse embryonic, neonatal, and adult tissues. Among other things, this analysis revealed that both the level and the subcellular localization of RHOX10 are regulated during germ cell development. To understand the role of Rhox10 in germ cell development, I generated transgenic mice expressing an artificial microRNA (miRNA) targeting Rhox10. While this artificial miRNA robustly downregulated RHOX10 protein expression in vitro, it did not significantly reduce RHOX10 expression in vivo. So I next elected to knockdown RHOX10 levels in spermatogonial stem cells (SSCs), which I found highly express both Rhox10 mRNA and RHOX10 protein. Using a recently developed in vitro culture system for SSCs combined with a short-hairpin RNA (shRNA) approach, I strongly depleted RHOX10 expression in SSCs. These RHOX10-depleted cells exhibited a defect in the ability to form stem cell clusters in vitro. Expression profiling analysis revealed many genes regulated by Rhox10, including many meiotic genes, which could be downstream of Rhox10 in a molecular pathway that controls SSC differentiation. ^ RNA recognition motif (RRM) containing protein, Hermes is localized in germ plasm, where dormant mRNAs are also located, of Xenopus oocytes, which implicates its role in translational regulator. To understand the function of Hermes in oocyte meiosis, I used a morpholino oligonucleotide (MO) based knockdown approach. Microinjection of Hermes MO into fully grown oocytes, which are arrested in meiotic prophase, caused acceleration of oocytes reentry into meiosis (i.e., maturation) upon progesterone induction. Using a candidate approach, I identified at least three targets of Hermes: Ringo/Spy, Xcat2, and Mos. Ringo/Spy and Mos are known to have functions in oocyte maturation, while Ringo/Spy, Xcat2 mRNA are localized in the germ plasm of oocytes, which drives germ cell specification after fertilization. This led me to propose that Hermes functions in both oocyte maturation and germ cell development through its ability to regulate 3 crucial target mRNAs. ^

Activity-Regulated microRNAs: Modulators of Synaptic Growth at the Drosophila Neuromuscular Junction

It is well established that long-term changes in synaptic structure and function are mediated by rapid activity-dependent gene transcription and new protein synthesis. A growing body of evidence supports the involvement of the microRNA (miRNA) pathway in these processes. We have used the Drosophila neuromuscular junction (NMJ) as a model synapse to characterize activity-regulated miRNAs and their important mRNA targets. Here, we have identified five neuronal miRNAs (miRs-1, -8, -289, -314, and -958) that are significantly downregulated in response to neuronal activity. Furthermore we have discovered that neuronal misexpression of three of these miRNAs (miR-8, -289, and -958) is capable of suppressing new synaptic growth in response to activity suggesting that these miRNAs control the translation of biologically relevant target mRNAs. Putative targets of the activity-regulated miRNAs-8 and -289 are significantly enriched in clusters mapping to functional processes including axon development, pathfinding, and axon growth. We demonstrate that activity-regulated miR-8 regulates the 3'UTR of wingless, a presynaptic regulatory protein involved in the process of activity-dependent axon terminal growth. Additionally, we show that the 3'UTR of the protein tyrosine phosophatase leukocyte antengen related (lar), a protein required for axon guidance and synaptic growth, is regulated by activity-regulated miRNAs-8, -289, and -958 in vitro. Both wg and lar were identified as relevant putative targets for co-regulation based through our functional cluster analysis. One putative target of miR-289 is the Ca2+/calmodulin-dependent protein kinase II (CamKII). While CamKII is not predicted as a target for co-regulation by multiple activity-regulated miRNAs we identified it as an especially pertinent target for analysis in our system for two reasons. First, CamKII has an extremely well characterized role in postsynaptic plasticity, but its presynaptic role is less well characterized and bears further analysis. Second, local translation of CamKII mRNA is regulated in part by the miRNA pathway in an activity-dependent manner in dendrites. We find that the CamKII 3'UTR is regulated by miR-289 in-vitro and this regulation is alleviated by mutating the `seed region' of the miR-289 binding site within the CamKII 3'UTR. Furthermore, we demonstrate a requirement for local translation of CamKII in motoneurons in the process of activity-regulated axon terminal growth.

Neuropsychiatric disease (NPD) clustering in families with Autism Spectrum Disorder (ASD): genetic, epigenetic and environmental issues

Autism Spectrum Disorder (ASD) is a neurodevelopmental disorder characterized by deficits in social communication/interaction and by unusual repetitive and restricted behaviors and interests. ASD often co-occurs in the same families with other neuropsychiatric diseases (NPD), such as intellectual disability, schizophrenia, epilepsy, depression and attention deficit hyperactivity disorder. Genetic factors have an important role in ASD etiology. Multiple copy number variants (CNVs) and single nucleotide variants (SNVs) in candidate genes have been associated with an increased risk to develop ASD. Nevertheless, recent heritability estimates and the high genotypic and phenotypic heterogeneity characteristic of ASD indicate a role of environmental and epigenetic factors, such as long noncoding RNA (lncRNA) and microRNA (miRNA), as modulators of genetic expression and further clinical presentation. Both miRNA and lncRNA are functional RNA molecules that are transcribed from DNA but not translated into proteins, instead they act as powerful regulators of gene expression. While miRNA are small noncoding RNAs with 22-25 nucleotides in length that act at the post-transcriptional level of gene expression, the lncRNA are bigger molecules (>200 nucleotides in length) that are capped, spliced, and polyadenylated, similar to messenger RNA. Although few lncRNA were well characterized until date, there is a great evidence that they are implicated in several levels of gene expression (transcription/post-transcription/post-translation, organization of protein complexes, cell– cell signaling as well as recombination) as shown in figure 1.

Expression and cellular localization of microRNA-29b and RAX, an activator of the RNA-dependent protein kinase (PKR), in the retina of streptozotocin-induced diabetic rats

Purpose: The apoptosis of retinal neurons plays a critical role in the pathogenesis of diabetic retinopathy (DR), but the molecular mechanisms underlying this phenomenon remain unclear. The purpose of this study was to investigate the cellular localization and the expression of microRNA-29b (miR-29b) and its potential target PKR associated protein X (RAX), an activator of the pro-apoptotic RNA-dependent protein kinase (PKR) signaling pathway, in the retina of normal and diabetic rats. Methods: Retinas were obtained from normal and diabetic rats within 35 days after streptozotocin (STZ) injection. In silico analysis indicated that RAX is a potential target of miR-29b. The cellular localization of miR-29b and RAX was assessed by in situ hybridization and immunofluorescence, respectively. The expression levels of miR-29b and RAX mRNA were evaluated by quantitative reverse transcription PCR (qRT-PCR), and the expression of RAX protein was evaluated by western blot. A luciferase reporter assay and inhibition of endogenous RAX were performed to confirm whether RAX is a direct target of miR-29b as predicted by the in silico analysis. Results: We found that miR-29b and RAX are localized in the retinal ganglion cells (RGCs) and the cells of the inner nuclear layer (INL) of the retinas from normal and diabetic rats. Thus, the expression of miR-29b and RAX, as assessed in the retina by quantitative RT-PCR, reflects their expression in the RGCs and the cells of the INL. We also revealed that RAX protein is upregulated (more than twofold) at 3, 6, 16, and 22 days and downregulated (70%) at 35 days, whereas miR-29b is upregulated (more than threefold) at 28 and 35 days after STZ injection. We did not confirm the computational prediction that RAX is a direct target of miR-29b. Conclusions: Our results suggest that RAX expression may be indirectly regulated by miR-29b, and the upregulation of this miRNA at the early stage of STZ-induced diabetes may have a protective effect against the apoptosis of RGCs and cells of the INL by the pro-apoptotic RNA-dependent protein kinase (PKR) signaling pathway.

Ultra-Deep Sequencing Reveals the microRNA Expression Pattern of the Human Stomach

Background: While microRNAs (miRNAs) play important roles in tissue differentiation and in maintaining basal physiology, little is known about the miRNA expression levels in stomach tissue. Alterations in the miRNA profile can lead to cell deregulation, which can induce neoplasia. Methodology/Principal Findings: A small RNA library of stomach tissue was sequenced using high-throughput SOLiD sequencing technology. We obtained 261,274 quality reads with perfect matches to the human miRnome, and 42% of known miRNAs were identified. Digital Gene Expression profiling (DGE) was performed based on read abundance and showed that fifteen miRNAs were highly expressed in gastric tissue. Subsequently, the expression of these miRNAs was validated in 10 healthy individuals by RT-PCR showed a significant correlation of 83.97% (P<0.05). Six miRNAs showed a low variable pattern of expression (miR-29b, miR-29c, miR-19b, miR-31, miR-148a, miR-451) and could be considered part of the expression pattern of the healthy gastric tissue. Conclusions/Significance: This study aimed to validate normal miRNA profiles of human gastric tissue to establish a reference profile for healthy individuals. Determining the regulatory processes acting in the stomach will be important in the fight against gastric cancer, which is the second-leading cause of cancer mortality worldwide.

miRNA analysis of prostate cancer by quantitative real time PCR: Comparison between formalin-fixed paraffin embedded and fresh-frozen tissue

Objective: Micro RNA (miRNA) is a class of small noncoding RNA that plays a major role in the regulation of gene expression, which has been related to cancer behavior. The possibility of analyzing miRNA from the archives of pathology laboratories is exciting, as it allows for large retrospective studies. Formalin is the most common fixative used in the surgical pathology routine, and its promotion of nucleic acid degradation is well known. Our aim is to compare miRNA profiles from formalin-fixed paraffin embedded (FFPE) tissues with fresh-frozen prostate cancer tissues. Methods: The expression of 14 miRNAs was determined by quantitative real time polymerase chain reaction (qRT-PCR) in 5 paired fresh-frozen and FFPE tissues, which were representative of prostate carcinoma. Results: There was a very good correlation of the miRNA expression of miR-let7c and miR-32 between the fresh-frozen and FFPE tissues, with Pearson`s correlation coefficients of 0.927 (P = 0.023) and 0.960 (P = 0.010), respectively. For the remaining miRNAs, the correlation was good with Spearman correlation coefficient of 0.638 (P < 0.001). Conclusion: Analysis of miRNAs from routinely processed and stored FFPE prostate tissue is feasible for some miRNAs using qRT-PCR. Further studies should be conducted to confirm the reliability of using stock tissues for miRNA expression determination. (C) 2011 Elsevier Inc. All rights reserved.

Conserved microRNA editing in mammalian evolution, development and disease.

BACKGROUND: Mammalian microRNAs (miRNAs) are sometimes subject to adenosine-to-inosine RNA editing, which can lead to dramatic changes in miRNA target specificity or expression levels. However, although a few miRNAs are known to be edited at identical positions in human and mouse, the evolution of miRNA editing has not been investigated in detail. In this study, we identify conserved miRNA editing events in a range of mammalian and non-mammalian species. RESULTS: We demonstrate deep conservation of several site-specific miRNA editing events, including two that date back to the common ancestor of mammals and bony fishes some 450 million years ago. We also find evidence of a recent expansion of an edited miRNA family in placental mammals and show that editing of these miRNAs is associated with changes in target mRNA expression during primate development and aging. While global patterns of miRNA editing tend to be conserved across species, we observe substantial variation in editing frequencies depending on tissue, age and disease state: editing is more frequent in neural tissues compared to heart, kidney and testis; in older compared to younger individuals; and in samples from healthy tissues compared to tumors, which together suggests that miRNA editing might be associated with a reduced rate of cell proliferation. CONCLUSIONS: Our results show that site-specific miRNA editing is an evolutionarily conserved mechanism, which increases the functional diversity of mammalian miRNA transcriptomes. Furthermore, we find that although miRNA editing is rare compared to editing of long RNAs, miRNAs are greatly overrepresented among conserved editing targets.

Birth and expression evolution of mammalian microRNA genes.

MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression, yet their origins and functional evolution in mammals remain little understood due to the lack of appropriate comparative data. Using RNA sequencing, we have generated extensive and comparable miRNA data for five organs in six species that represent all main mammalian lineages and birds (the evolutionary outgroup) with the aim to unravel the evolution of mammalian miRNAs. Our analyses reveal an overall expansion of miRNA repertoires in mammals, with threefold accelerated birth rates of miRNA families in placentals and marsupials, facilitated by the de novo emergence of miRNAs in host gene introns. Generally, our analyses suggest a high rate of miRNA family turnover in mammals with many newly emerged miRNA families being lost soon after their formation. Selectively preserved mammalian miRNA families gradually evolved higher expression levels, as well as altered mature sequences and target gene repertoires, and were apparently mainly recruited to exert regulatory functions in nervous tissues. However, miRNAs that originated on the X chromosome evolved high expression levels and potentially diverse functions during spermatogenesis, including meiosis, through selectively driven duplication-divergence processes. Overall, our study thus provides detailed insights into the birth and evolution of mammalian miRNA genes and the associated selective forces.

A microRNA Signature Associated with Early Recurrence in Breast Cancer.

Recurrent breast cancer occurring after the initial treatment is associated with poor outcome. A bimodal relapse pattern after surgery for primary tumor has been described with peaks of early and late recurrence occurring at about 2 and 5 years, respectively. Although several clinical and pathological features have been used to discriminate between low- and high-risk patients, the identification of molecular biomarkers with prognostic value remains an unmet need in the current management of breast cancer. Using microarray-based technology, we have performed a microRNA expression analysis in 71 primary breast tumors from patients that either remained disease-free at 5 years post-surgery (group A) or developed early (group B) or late (group C) recurrence. Unsupervised hierarchical clustering of microRNA expression data segregated tumors in two groups, mainly corresponding to patients with early recurrence and those with no recurrence. Microarray data analysis and RT-qPCR validation led to the identification of a set of 5 microRNAs (the 5-miRNA signature) differentially expressed between these two groups: miR-149, miR-10a, miR-20b, miR-30a-3p and miR-342-5p. All five microRNAs were down-regulated in tumors from patients with early recurrence. We show here that the 5-miRNA signature defines a high-risk group of patients with shorter relapse-free survival and has predictive value to discriminate non-relapsing versus early-relapsing patients (AUC = 0.993, p-value<0.05). Network analysis based on miRNA-target interactions curated by public databases suggests that down-regulation of the 5-miRNA signature in the subset of early-relapsing tumors would result in an overall increased proliferative and angiogenic capacity. In summary, we have identified a set of recurrence-related microRNAs with potential prognostic value to identify patients who will likely develop metastasis early after primary breast surgery.

A neuron-specific deletion of the microRNA-processing enzyme DICER induces severe but transient obesity in mice.

MicroRNAs (miRNAs) are small, non-coding RNA molecules that regulate gene expression post-transcriptionally. MiRNAs are implicated in various biological processes associated with obesity, including adipocyte differentiation and lipid metabolism. We used a neuronal-specific inhibition of miRNA maturation in adult mice to study the consequences of miRNA loss on obesity development. Camk2a-CreERT2 (Cre+) and floxed Dicer (Dicerlox/lox) mice were crossed to generate tamoxifen-inducible conditional Dicer knockouts (cKO). Vehicle- and/or tamoxifen-injected Cre+;Dicerlox/lox and Cre+;Dicer+/+ served as controls. Four cohorts were used to a) measure body composition, b) follow food intake and body weight dynamics, c) evaluate basal metabolism and effects of food deprivation, and d) assess the brain transcriptome consequences of miRNA loss. cKO mice developed severe obesity and gained 18 g extra weight over the 5 weeks following tamoxifen injection, mainly due to increased fat mass. This phenotype was highly reproducible and observed in all 38 cKO mice recorded and in none of the controls, excluding possible effects of tamoxifen or the non-induced transgene. Development of obesity was concomitant with hyperphagia, increased food efficiency, and decreased activity. Surprisingly, after reaching maximum body weight, obese cKO mice spontaneously started losing weight as rapidly as it was gained. Weight loss was accompanied by lowered O2-consumption and respiratory-exchange ratio. Brain transcriptome analyses in obese mice identified several obesity-related pathways (e.g. leptin, somatostatin, and nemo-like kinase signaling), as well as genes involved in feeding and appetite (e.g. Pmch, Neurotensin) and in metabolism (e.g. Bmp4, Bmp7, Ptger1, Cox7a1). A gene cluster with anti-correlated expression in the cerebral cortex of post-obese compared to obese mice was enriched for synaptic plasticity pathways. While other studies have identified a role for miRNAs in obesity, we here present a unique model that allows for the study of processes involved in reversing obesity. Moreover, our study identified the cortex as a brain area important for body weight homeostasis.

The lincRNA HOTAIRM1, located in the HOXA genomic region, is expressed in acute myeloid leukemia, impacts prognosis in patients in the intermediate-risk cytogenetic category, and is associated with a distinctive microRNA signature.

Long non-coding RNAs (lncRNAs) are deregulated in several tumors, although their role in acute myeloid leukemia (AML) is mostly unknown.We have examined the expression of the lncRNA HOX antisense intergenic RNA myeloid 1 (HOTAIRM1) in 241 AML patients. We have correlated HOTAIRM1 expression with a miRNA expression profile. We have also analyzed the prognostic value of HOTAIRM1 expression in 215 intermediate-risk AML (IR-AML) patients.The lowest expression level was observed in acute promyelocytic leukemia (P < 0.001) and the highest in t(6;9) AML (P = 0.005). In 215 IR-AML patients, high HOTAIRM1 expression was independently associated with shorter overall survival (OR:2.04;P = 0.001), shorter leukemia-free survival (OR:2.56; P < 0.001) and a higher cumulative incidence of relapse (OR:1.67; P = 0.046). Moreover, HOTAIRM1 maintained its independent prognostic value within the favorable molecular subgroup (OR: 3.43; P = 0.009). Interestingly, HOTAIRM1 was overexpressed in NPM1-mutated AML (P < 0.001) and within this group retained its prognostic value (OR: 2.21; P = 0.01). Moreover, HOTAIRM1 expression was associated with a specific 33-microRNA signature that included miR-196b (P < 0.001). miR-196b is located in the HOX genomic region and has previously been reported to have an independent prognostic value in AML. miR-196b and HOTAIRM1 in combination as a prognostic factor can classify patients as high-, intermediate-, or low-risk (5-year OS: 24% vs 42% vs 70%; P = 0.004).Determination of HOTAIRM1 level at diagnosis provided relevant prognostic information in IR-AML and allowed refinement of risk stratification based on common molecular markers. The prognostic information provided by HOTAIRM1 was strengthened when combined with miR-196b expression. Furthermore, HOTAIRM1 correlated with a 33-miRNA signature.

MicroRNA-122 modulates the rhythmic expression profile of the circadian deadenylase Nocturnin in mouse liver.

Nocturnin is a circadian clock-regulated deadenylase thought to control mRNA expression post-transcriptionally through poly(A) tail removal. The expression of Nocturnin is robustly rhythmic in liver at both the mRNA and protein levels, and mice lacking Nocturnin are resistant to diet-induced obesity and hepatic steatosis. Here we report that Nocturnin expression is regulated by microRNA-122 (miR-122), a liver specific miRNA. We found that the 3'-untranslated region (3'-UTR) of Nocturnin mRNA harbors one putative recognition site for miR-122, and this site is conserved among mammals. Using a luciferase reporter construct with wild-type or mutant Nocturnin 3'-UTR sequence, we demonstrated that overexpression of miR-122 can down-regulate luciferase activity levels and that this effect is dependent on the presence of the putative miR-122 recognition site. Additionally, the use of an antisense oligonucleotide to knock down miR-122 in vivo resulted in significant up-regulation of both Nocturnin mRNA and protein expression in mouse liver during the night, resulting in Nocturnin rhythms with increased amplitude. Together, these data demonstrate that the normal rhythmic profile of Nocturnin expression in liver is shaped in part by miR-122. Previous studies have implicated Nocturnin and miR-122 as important post-transcriptional regulators of both lipid metabolism and circadian clock controlled gene expression in the liver. Therefore, the demonstration that miR-122 plays a role in regulating Nocturnin expression suggests that this may be an important intersection between hepatic metabolic and circadian control.

Design and evaluation of a panel os single-nucleotide polymorphisms in microRNA genomic regions for association studies in human disease

MicroRNAs (miRNA) are recognized posttranscriptional gene repressors involved in the control of almost every biological process. Allelic variants in these regions may be an important source of phenotypic diversity and contribute to disease susceptibility. We analyzed the genomic organization of 325 human miRNAs (release 7.1, miRBase) to construct a panel of 768 single-nucleotide polymorphisms (SNPs) covering approximately 1 Mb of genomic DNA, including 131 isolated miRNAs (40%) and 194 miRNAs arranged in 48 miRNA clusters, as well as their 5-kb flanking regions. Of these miRNAs, 37% were inside known protein-coding genes, which were significantly associated with biological functions regarding neurological, psychological or nutritional disorders. SNP coverage analysis revealed a lower SNP density in miRNAs compared with the average of the genome, with only 24 SNPs located in the 325 miRNAs studied. Further genotyping of 340 unrelated Spanish individuals showed that more than half of the SNPs in miRNAs were either rare or monomorphic, in agreement with the reported selective constraint on human miRNAs. A comparison of the minor allele frequencies between Spanish and HapMap population samples confirmed the applicability of this SNP panel to the study of complex disorders among the Spanish population, and revealed two miRNA regions, hsa-mir-26a-2 in the CTDSP2 gene and hsa-mir-128-1 in the R3HDM1 gene, showing geographical allelic frequency variation among the four HapMap populations, probably because of differences in natural selection. The designed miRNA SNP panel could help to identify still hidden links between miRNAs and human disease.