miR2GO: comparative functional analysis for microRNAs

Investigation of the role of miRNA variants in neurodegenerative brain diseases

Introduction

miRNAs are small noncoding elements known to regulate different molecular processes, including developmental and executive functions in the brain. Dysregulation of miRNAs could contribute to brain neurodegeneration, as suggested by miRNA profiling studies of individuals suffering from neurodegenerative brain diseases (NBDs). Here, we report rare miRNA variants in patients with Alzheimer's dementia (AD) and frontotemporal dementia (FTD).

Methods

We initially used whole exome sequencing data in a subset of FTD patients (n = 209) from Flanders-Belgium. We then performed targeted resequencing of variant-harboring miRNAs in an additional subset of FTD patients (n = 126) and control individuals (n = 426). Lastly, we sequenced the MIR885 locus in a Flanders-Belgian AD cohort (n = 947) and a total number of n = 755 controls.

Results

WES identified rare seed variants in MIR656, MIR423, MIR122 and MIR885 in FTD patients. Most of these miRNAs bind to FTD-associated genes, implicated in different biological pathways. Additionally, some miRNA variants create novel binding sites for genes associated with FTD. Sequencing of the MIR885 locus in the AD cohort initially showed a significant enrichment of MIR885 variants in AD patients compared to controls (SKAT-O, p-value = 0.026). Genetic association was not maintained when we included sex and APOE status as covariates. Using the miRVaS prediction tool, variants rs897551430 and rs993255773 appeared to evoke significant structural changes in the primary miRNA. These variants are also predicted to strongly downregulate mature miR885 levels, in line with what is reported for MIR885 in the context of AD.

Discussion

Functional investigation of miRNAs/variants described in this study could propose novel miRNA-mediated molecular cascades in FTD and AD pathogenicity. Furthermore, we believe that the genetic evidence presented here suggests a role for MIR885 in molecular mechanisms involved in AD and warrants genetic follow-up in larger cohorts to explore this hypothesis.

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.

Functional characterization of human genomic variation linked to polygenic diseases

The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.

MiRLog and dbmiR: prioritization and functional annotation tools to study human microRNA sequence variants

The recent identification of noncoding variants with pathogenic effects suggests that these variations could underlie a significant number of undiagnosed cases. Several computational methods have been developed to predict the functional impact of noncoding variants, but they exhibit only partial concordance and are not integrated with functional annotation resources, making the interpretation of these variants still challenging. MicroRNAs (miRNAs) are small noncoding RNA molecules that act as fine regulators of gene expression and play crucial functions in several biological processes, such as cell proliferation and differentiation. An increasing number of studies demonstrate a significant impact of miRNA single nucleotide variants (SNVs) both in Mendelian diseases and complex traits. To predict the functional effect of miRNA SNVs, we implemented a new meta‐predictor, MiRLog, and we integrated it into a comprehensive database, dbmiR, which includes a precompiled list of all possible miRNA allelic SNVs, providing their biological annotations at nucleotide and miRNA levels. MiRLog and dbmiR were used to explore the genetic variability of miRNAs in 15,708 human genomes included in the gnomAD project, finding several ultra‐rare SNVs with a potentially deleterious effect on miRNA biogenesis and function representing putative contributors to human phenotypes.

MicroRNA Targeting

MicroRNAs (miRNAs) are small noncoding elements that play essential roles in the posttranscriptional regulation of biochemical processes. miRNAs recognize and target multiple mRNAs; therefore, investigating miRNA dysregulation is an indispensable strategy to understand pathological conditions and to design innovative drugs. Targeting miRNAs in diseases improve outcomes of several therapeutic strategies thus, this present study highlights miRNA targeting methods through experimental assays and bioinformatics tools. The first part of this review focuses on experimental miRNA targeting approaches for elucidating key biochemical pathways. A growing body of evidence about the miRNA world reveals the fact that it is not possible to uncover these molecules' structural and functional characteristics related to the biological processes with a deterministic approach. Instead, a systemic point of view is needed to truly understand the facts behind the natural complexity of interactions and regulations that miRNA regulations present. This task heavily depends both on computational and experimental capabilities. Fortunately, several miRNA bioinformatics tools catering to nonexperts are available as complementary wet-lab approaches. For this purpose, this work provides recent research and information about computational tools for miRNA targeting research.

Bioinformatics Accelerates the Major Tetrad: A Real Boost for the Pharmaceutical Industry

With advanced technology and its development, bioinformatics is one of the avant-garde fields that has managed to make amazing progress in the pharmaceutical-medical field by modeling the infrastructural dimensions of healthcare and integrating computing tools in drug innovation, facilitating prevention, detection/more accurate diagnosis, and treatment of disorders, while saving time and money. By association, bioinformatics and pharmacovigilance promoted both sample analyzes and interpretation of drug side effects, also focusing on drug discovery and development (DDD), in which systems biology, a personalized approach, and drug repositioning were considered together with translational medicine. The role of bioinformatics has been highlighted in DDD, proteomics, genetics, modeling, miRNA discovery and assessment, and clinical genome sequencing. The authors have collated significant data from the most known online databases and publishers, also narrowing the diversified applications, in order to target four major areas (tetrad): DDD, anti-microbial research, genomic sequencing, and miRNA research and its significance in the management of current pandemic context. Our analysis aims to provide optimal data in the field by stratification of the information related to the published data in key sectors and to capture the attention of researchers interested in bioinformatics, a field that has succeeded in advancing the healthcare paradigm by introducing developing techniques and multiple database platforms, addressed in the manuscript.

Adenosine-to-Inosine Editing of Vasoactive MicroRNAs Alters Their Targetome and Function in Ischemia

Adenosine-to-inosine (A-to-I) editing in the seed sequence of microRNAs can shift the microRNAs’ targetomes and thus their function. Using public RNA-sequencing data, we identified 35 vasoactive microRNAs that are A-to-I edited. We quantified A-to-I editing of the primary (pri-)microRNAs in vascular fibroblasts and endothelial cells. Nine pri-microRNAs were indeed edited, and editing consistently increased under ischemia. We determined mature microRNA editing for the highest expressed microRNAs, i.e., miR-376a-3p, miR-376c-3p, miR-381-3p, and miR-411-5p. All four mature microRNAs were edited in their seed sequence. We show that both ADAR1 and ADAR2 (adenosine deaminase acting on RNA 1 and RNA 2) can edit pri-microRNAs in a microRNA-specific manner. MicroRNA editing also increased under ischemia in vivo in a murine hindlimb ischemia model and ex vivo in human veins. For each edited microRNA, we confirmed a shift in targetome. Expression of the edited microRNA targetomes, not the wild-type targetomes, was downregulated under ischemia in vivo. Furthermore, microRNA editing enhanced angiogenesis in vitro and ex vivo. In conclusion, we show that microRNA A-to-I editing is a widespread phenomenon, induced by ischemia. Each editing event results in a novel microRNA with a unique targetome, leading to increased angiogenesis.

Computational annotation of miRNA transcription start sites

Motivation

MicroRNAs (miRNAs) are small noncoding RNAs that play important roles in gene regulation and phenotype development. The identification of miRNA transcription start sites (TSSs) is critical to understand the functional roles of miRNA genes and their transcriptional regulation. Unlike protein-coding genes, miRNA TSSs are not directly detectable from conventional RNA-Seq experiments due to miRNA-specific process of biogenesis. In the past decade, large-scale genome-wide TSS-Seq and transcription activation marker profiling data have become available, based on which, many computational methods have been developed. These methods have greatly advanced genome-wide miRNA TSS annotation.

Results

In this study, we summarized recent computational methods and their results on miRNA TSS annotation. We collected and performed a comparative analysis of miRNA TSS annotations from 14 representative studies. We further compiled a robust set of miRNA TSSs (RSmirT) that are supported by multiple studies. Integrative genomic and epigenomic data analysis on RSmirT revealed the genomic and epigenomic features of miRNA TSSs as well as their relations to protein-coding and long non-coding genes.

Contact

xiaoman@mail.ucf.edu, haihu@cs.ucf.edu

miRPathDB 2.0: a novel release of the miRNA Pathway Dictionary Database

Since the initial release of miRPathDB, tremendous progress has been made in the field of microRNA (miRNA) research. New miRNA reference databases have emerged, a vast amount of new miRNA candidates has been discovered and the number of experimentally validated target genes has increased considerably. Hence, the demand for a major upgrade of miRPathDB, including extended analysis functionality and intuitive visualizations of query results has emerged. Here, we present the novel release 2.0 of the miRNA Pathway Dictionary Database (miRPathDB) that is freely accessible at https://mpd.bioinf.uni-sb.de/. miRPathDB 2.0 comes with a ten-fold increase of pre-processed data. In total, the updated database provides putative associations between 27 452 (candidate) miRNAs, 28 352 targets and 16 833 pathways for Homo sapiens, as well as interactions of 1978 miRNAs, 24 898 targets and 6511 functional categories for Mus musculus. Additionally, we analyzed publications citing miRPathDB to identify common use-cases and further extensions. Based on this evaluation, we added new functionality for interactive visualizations and down-stream analyses of bulk queries. In summary, the updated version of miRPathDB, with its new custom-tailored features, is one of the most comprehensive and advanced resources for miRNAs and their target pathways.

Gene Ontology Semantic Similarity Analysis Using GOSemSim

The GOSemSim package, an R-based tool within the Bioconductor project, offers several methods based on information content and graph structure for measuring semantic similarity among GO terms, gene products and gene clusters. In this chapter, I illustrate the use of GOSemSim on a list of regulators in preimplantation embryos. A step-by-step analysis was provided as well as instructions on interpretation and visualization of the results. GOSemSim is open-source and is available from https://www.bioconductor.org/packages/GOSemSim .

MicroRNA-411 and Its 5'-IsomiR Have Distinct Targets and Functions and Are Differentially Regulated in the Vasculature under Ischemia

MicroRNAs are posttranscriptional regulators of gene expression. As microRNAs can target many genes simultaneously, microRNAs can regulate complex multifactorial processes, including post-ischemic neovascularization, a major recovery pathway in cardiovascular disease. MicroRNAs select their target mRNAs via full complementary binding with their seed sequence, i.e., nucleotides 2-8 from the 5' end of a microRNA. The exact sequence of a mature microRNA, and thus of its 5' and 3' ends, is determined by two sequential cleavage steps of microRNA precursors, Drosha/DGCR8 and Dicer. When these cleavage steps result in nucleotide switches at the 5' end, forming a so-called 5'-isomiR, this results in a shift in the mature microRNA's seed sequence. The role of 5'-isomiRs in cardiovascular diseases is still unknown. Here, we characterize the expression and function of the 5'-isomiR of miR-411 (ISO-miR-411). ISO-miR-411 is abundantly expressed in human primary vascular cells. ISO-miR-411 has a different "targetome" from WT-miR-411, with only minor overlap. The ISO-miR-411/WT-miR-411 ratio is downregulated under acute ischemia, both in cells and a murine ischemia model, but is upregulated instead in chronically ischemic human blood vessels. ISO-miR-411 negatively influences vascular cell migration, whereas WT-miR-411 does not. Our data demonstrate that isomiR formation is a functional pathway that is actively regulated during ischemia.

Computational Resources for Prediction and Analysis of Functional miRNA and Their Targetome

microRNAs are evolutionarily conserved, endogenously produced, noncoding RNAs (ncRNAs) of approximately 19-24 nucleotides (nts) in length known to exhibit gene silencing of complementary target sequence. Their deregulated expression is reported in various disease conditions and thus has therapeutic implications. In the last decade, various computational resources are published in this field. In this chapter, we have reviewed bioinformatics resources, i.e., miRNA-centered databases, algorithms, and tools to predict miRNA targets. First section has enlisted more than 75 databases, which mainly covers information regarding miRNA registries, targets, disease associations, differential expression, interactions with other noncoding RNAs, and all-in-one resources. In the algorithms section, we have compiled about 140 algorithms from eight subcategories, viz. for the prediction of precursor (pre-) and mature miRNAs. These algorithms are developed on various sequence, structure, and thermodynamic based features incorporated into different machine learning techniques (MLTs). In addition, computational identification of miRNAs from high-throughput next generation sequencing (NGS) data and their variants, viz. isomiRs, differential expression, miR-SNPs, and functional annotation, are discussed. Prediction and analysis of miRNAs and their associated targets are also evaluated under miR-targets section providing knowledge regarding novel miRNA targets and complex host-pathogen interactions. In conclusion, we have provided comprehensive review of in silico resources published in miRNA research to help scientific community be updated and choose the appropriate tool according to their needs.

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.

Virus-encoded miRNAs in Ebola virus disease

Ebola virus (EBOV) is a negative-strand RNA virus that replicates in the cytoplasm and causes an often-fatal hemorrhagic fever. EBOV, like other viruses, can reportedly encode its own microRNAs (miRNAs) to subvert host immune defenses. miRNAs are short noncoding RNAs that can regulate gene expression by hybridizing to multiple mRNAs, and viral miRNAs can enhance viral replication and infectivity by regulating host or viral genes. To date, only one EBOV miRNA has been examined in human infection. Here, we assayed mouse, rhesus macaque, cynomolgus macaque, and human samples infected with three EBOV variants for twelve computationally predicted viral miRNAs using RT-qPCR. Ten miRNAs aligned to EBOV variants and were detectable in the four species during disease with several viral miRNAs showing presymptomatic amplification in animal models. miRNA abundances in both the mouse and nonhuman primate models mirrored the human cohort, with miR-1-5p, miR-1-3p, and miR-T3-3p consistently at the highest levels. These striking similarities in the most abundant miRNAs during infection with different EBOV variants and hosts indicate that these miRNAs are potential valuable diagnostic markers and key effectors of EBOV pathogenesis.

Structure and functional impact of seed region variant in MIR-499 gene family in bronchial asthma

BACKGROUND

Small non-coding RNAs (microRNAs) have been evolved to master numerous cellular processes. Genetic variants within microRNA seed region might influence microRNA biogenesis and function. The study aimed at determining the role of microRNA-499 (MIR-499) gene family polymorphism as a marker for susceptibility and progression of bronchial asthma and to analyze the structural and functional impact of rs3746444 within the seed region.

METHODS

Genotyping for 192 participants (96 patients and 96 controls) in the discovery phase and 319 subjects (115 patients and 204 controls) in the replication phase was performed via Real Time-Polymerase Chain Reaction technology. Patients underwent the methacholine challenge test and biochemical analysis. Gene structural and functional analysis, target prediction, annotation clustering, and pathway enrichment analysis were executed. Predicted functional effect of rs37464443 SNP was analyzed.

RESULTS

miR-499 gene family is highly implicated in inflammation-related signaling pathways. Rs374644 (A > G) in MIR499A and MIR499B within the seed region could disrupt target genes and create new genes. The G variant was associated with high risk of developing asthma under all genetic association models (G versus A: OR = 3.27, 95% CI = 2.53-4.22; GG versus AA: OR = 9.52, 95% CI = 5.61-16.5; AG versus AA: OR = 2.13, 95% CI = 1.24-3.46; GG + AG versus AA: OR = 4.43, 95% CI = 2.88-6.82). GG genotype was associated with poor pre-bronchodilator FEV1 (p = 0.047) and the worst bronchodilator response after Salbutamol inhalation, represented in low peaked expiratory flow rate (p = 0.035).

CONCLUSIONS

miR-499 rs3746444 (A > G) polymorphism was associated with asthma susceptibility and bronchodilator response in Egyptian children and adolescents. Further functional analysis is warranted to develop more specific theranostic agents for selecting targeted therapy.

A compilation of Web-based research tools for miRNA analysis

Since the discovery of microRNAs (miRNAs), a class of noncoding RNAs that regulate the gene expression posttranscriptionally in sequence-specific manner, there has been a release of number of tools useful for both basic and advanced applications. This is because of the significance of miRNAs in many pathophysiological conditions including cancer. Numerous bioinformatics tools that have been developed for miRNA analysis have their utility for detection, expression, function, target prediction and many other related features. This review provides a comprehensive assessment of web-based tools for the miRNA analysis that does not require prior knowledge of any computing languages.

FARNA: knowledgebase of inferred functions of non-coding RNA transcripts

Non-coding RNA (ncRNA) genes play a major role in control of heterogeneous cellular behavior. Yet, their functions are largely uncharacterized. Current available databases lack in-depth information of ncRNA functions across spectrum of various cells/tissues. Here, we present FARNA, a knowledgebase of inferred functions of 10,289 human ncRNA transcripts (2,734 microRNA and 7,555 long ncRNA) in 119 tissues and 177 primary cells of human. Since transcription factors (TFs) and TF co-factors (TcoFs) are crucial components of regulatory machinery for activation of gene transcription, cellular processes and diseases in which TFs and TcoFs are involved suggest functions of the transcripts they regulate. In FARNA, functions of a transcript are inferred from TFs and TcoFs whose genes co-express with the transcript controlled by these TFs and TcoFs in a considered cell/tissue. Transcripts were annotated using statistically enriched GO terms, pathways and diseases across cells/tissues based on guilt-by-association principle. Expression profiles across cells/tissues based on Cap Analysis of Gene Expression (CAGE) are provided. FARNA, having the most comprehensive function annotation of considered ncRNAs across widest spectrum of human cells/tissues, has a potential to greatly contribute to our understanding of ncRNA roles and their regulatory mechanisms in human. FARNA can be accessed at: http://cbrc.kaust.edu.sa/farna

Systematic Prediction of the Impacts of Mutations in MicroRNA Seed Sequences

MicroRNAs are a class of small non-coding RNAs that are involved in many important biological processes and the dysfunction of microRNA has been associated with many diseases. The seed region of a microRNA is of crucial importance to its target recognition. Mutations in microRNA seed regions may disrupt the binding of microRNAs to their original target genes and make them bind to new target genes. Here we use a knowledge-based computational method to systematically predict the functional effects of all the possible single nucleotide mutations in human microRNA seed regions. The result provides a comprehensive reference for the functional assessment of the impacts of possible natural and artificial single nucleotide mutations in microRNA seed regions.

Computational and Experimental Identification of Tissue-Specific MicroRNA Targets

In this chapter we discuss computational methods for the prediction of microRNA (miRNA) targets. More specifically, we consider machine learning-based approaches and explain why these methods have been relatively unsuccessful in reducing the number of false positive predictions. Further we suggest approaches designed to improve their performance by considering tissue-specific target regulation. We argue that the miRNA targetome differs depending on the tissue type and introduce a novel algorithm that predicts miRNA targets specifically for colorectal cancer. We discuss features of miRNAs and target sites that affect target recognition, and how next-generation sequencing data can support the identification of novel miRNAs, differentially expressed miRNAs and their tissue-specific mRNA targets. In addition, we introduce some experimental approaches for the validation of miRNA targets as well as web-based resources sharing predicted and validated miRNA target interactions.

Precursor miR-499a Variant but not miR-196a2 is Associated with Rheumatoid Arthritis Susceptibility in an Egyptian Population

INTRODUCTION

Rheumatoid arthritis (RA) has a complex component induced by several genes that interact together with environmental and hormonal factors. We aimed to investigate the association of miR-196a2 rs11614913 (C/T) and miR-499a rs3746444 (A/G) polymorphisms and their combination with RA susceptibility and disease activity in an Egyptian population, and to evaluate their impact on methotrexate drug response and toxicity.

MATERIALS AND METHODS

Bioinformatics databases were searched to select potential micro RNA (miRNA)-messenger RNA (mRNA) interactions involved in RA pathogenesis. Ninety-five RA patients diagnosed according to the American College of Rheumatology and 200 healthy controls were genotyped using real-time polymerase chain reaction technology.

RESULTS

In overall and stratified analysis, miR-499a, but not miR-196a2, was associated with RA risk. Heterozygote carriers with rs3746444*A/G displayed protection against developing RA (p = 0.005) with an odds ratio of 0.2 (95 % confidence interval 0.17-0.62). The carriage of the combinations (miR499a*AG + miR196a2*CC) and (miR499a*AA + miR196a2*TT) were 3 and 7.5 times more likely to develop RA, respectively, while the combinations (miR499a*GG + miR196a2*CC), (miR499a*AG + miR196a2*TT) and (miR499a*AA + miR196a2*CT) show less susceptibility to have RA disease (all p < 0.05). rs3746444*AA genotype had a higher disease activity score (DAS28) [p = 0.023], tender joint count (TJC) (p = 0.007), and methotrexate-induced gastrointestinal toxicity (p = 0.043) compared with both AG/GG genotypes. rs11614913*C carriers were associated with higher DAS28 activity (p = 0.021). Homozygote male patients (CC and TT) had higher TJC (p = 0.046) and higher rheumatoid factor levels (p = 0.026), whereas, TT homozygote females had higher levels of ALT (p = 0.022).

CONCLUSIONS

Different genotypes of miR-499a rs3746444 single nucleotide polymorphisms (SNPs) are associated with RA risk, disease activity, and methotrexate toxicity in our population. In combination with specific miR-196a2 rs11614913 genotypes, this risk could increase or decrease according to the type of combination. Further functional analysis of the SNP and its impact on mRNA targets is required to confirm the relationship between genotype and phenotype.

A passenger strand variant in miR-196a2 contributes to asthma severity in children and adolescents: A preliminary study

There is emerging evidence to support the role of microRNAs in allergic airway diseases and inflammation. Genetic variants in microRNA genes might affect microRNA-mediated cell regulation. This preliminary study was designed to investigate the association of the microRNA-196a2 rs11614913 (C/T) polymorphism with susceptibility to asthma and clinical outcomes in children and adolescents. Genotyping of rs11614913 polymorphism was determined in 96 patients with bronchial asthma (6-18 years of age) and 96 unrelated controls, using real-time polymerase chain reaction technology. In-silico target prediction and network core analyses were performed. The asthmatics did not show significant differences in genotype distribution (p = 0.609) and allele frequencies (p = 0.428) compared with the controls. There were also no associations with disease duration, age at onset, asthma phenotype, asthma control, therapeutic level, airway hyper-responsiveness, or biochemical parameters in the blood. However, the CC genotype was associated with a more severe degree of asthma (p = 0. 023) and higher frequency of nocturnal asthma (p = 0.002). Carriers for CC were 17 times more likely to develop nocturnal asthma, and had a more than 2.5-fold increased risk for poor disease outcome compared with CT and TT individuals. In conclusion, microRNA-196a2 rs11614913 polymorphism might be associated with asthma severity in our sample of the Egyptian population. Further investigations in studies with a larger sample size and functional tests are needed to validate our findings and to explore the detailed biological mechanisms.

ImiRP: a computational approach to microRNA target site mutation

BACKGROUND

MicroRNAs (miRNAs) are small ~22 nucleotide non-coding RNAs that function as post-transcriptional regulators of messenger RNA (mRNA) through base-pairing to 6-8 nucleotide long target sites, usually located within the mRNA 3' untranslated region. A common approach to validate and probe microRNA-mRNA interactions is to mutate predicted target sites within the mRNA and determine whether it affects miRNA-mediated activity. The introduction of miRNA target site mutations, however, is potentially problematic as it may generate new, "illegitimate sites" target sites for other miRNAs, which may affect the experimental outcome. While it is possible to manually generate and check single miRNA target site mutations, this process can be time consuming, and becomes particularly onerous and error prone when multiple sites are to be mutated simultaneously. We have developed a modular Java-based system called ImiRP (Illegitimate miRNA Predictor) to solve this problem and to facilitate miRNA target site mutagenesis.

RESULTS

The ImiRP interface allows users to input a sequence of interest, specify the locations of multiple predicted target sites to mutate, and set parameters such as species, mutation strategy, and disallowed illegitimate target site types. As mutant sequences are generated, ImiRP utilizes the miRBase high confidence miRNA dataset to identify illegitimate target sites in each mutant sequence by comparing target site predictions between input and mutant sequences. ImiRP then assembles a final mutant sequence in which all specified target sites have been mutated.

CONCLUSIONS

ImiRP is a mutation generator program that enables selective disruption of specified miRNA target sites while ensuring predicted target sites for other miRNAs are not inadvertently created. ImiRP supports mutagenesis of single and multiple miRNA target sites within a given sequence, including sites that overlap. This software will be particularly useful for studies looking at microRNA cooperativity, where mutagenesis of multiple microRNA target sites may be desired. The software is available at imirp.org and is available open source for download through GitHub ( https://github.com/imirp ).

Bioinformatic tools for microRNA dissection

Recently, microRNAs (miRNAs) have emerged as important elements of gene regulatory networks. MiRNAs are endogenous single-stranded non-coding RNAs (~22-nt long) that regulate gene expression at the post-transcriptional level. Through pairing with mRNA, miRNAs can down-regulate gene expression by inhibiting translation or stimulating mRNA degradation. In some cases they can also up-regulate the expression of a target gene. MiRNAs influence a variety of cellular pathways that range from development to carcinogenesis. The involvement of miRNAs in several human diseases, particularly cancer, makes them potential diagnostic and prognostic biomarkers. Recent technological advances, especially high-throughput sequencing, have led to an exponential growth in the generation of miRNA-related data. A number of bioinformatic tools and databases have been devised to manage this growing body of data. We analyze 129 miRNA tools that are being used in diverse areas of miRNA research, to assist investigators in choosing the most appropriate tools for their needs.

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.