dbSMR: a novel resource of genome-wide SNPs affecting microRNA mediated regulation

Genetic association of Nuclear factor of activated T cells' 3'UTR and structural polymorphisms with susceptibility to generalized vitiligo in Gujarat population

Generalized vitiligo(GV) is an autoimmune skin depigmenting disease characterized by loss of functional melanocytes. Nuclear factor of activated T cells(NFATs) play a key role in regulatory T cells' (Tregs) activation and function. Our previous studies have highlighted the role of reduced NFATs expression and activity in impaired Tregs suppressive capacity, leading to GV pathogenesis. 3'UTR region and structural single nucleotide polymorphisms(SNPs) could lead to reduced NFAT expression and activity. Therefore, we studied the association of NFATs 3'UTR [NFATC2 rs4811198(T>G) &NFATC4 rs11848279(A>G)] and structural [NFATC1 rs754093(T>G) &NFATC2 rs12479626(T>C)] SNPs in 427 GV patients and 415 controls from Gujarat population by Polymerase chain reaction-restriction fragment length polymorphism(PCR-RFLP). Additionally, we carried out genotype-phenotype correlation and in silico analysis to assess the effect of NFATs SNPs on NFATs expression and structure. NFATC2 rs4811198(T>G) 3' UTR &NFATC2 rs12479626(T>C) structural SNPs were significantly associated with GV(p<0.0001). Interestingly, for NFATC2 rs4811198(T>G) SNP, there was a significant difference in the TT vs GG genotypes' frequencies (p=0.0034; Table 2), and for NFATC2 rs12479626(T>C) SNP there was a significant difference between TT vs TC and CC genotypes' frequencies(p<0.0001 & p=0.0002) between GV patients and controls. Furthermore, Odds ratio suggested that the susceptible alleles for NFATC2 rs4811198(T>G) &NFATC2 rs12479626(T>C) SNPs increased the risk of GV by 1.38 & 3.043 fold. However, the NFAT 3' UTR [NFATC2 rs4811198 (T>G)] and structural [NFATC1 rs754093(T>G)] SNPs were not significantly associated with GV. Interestingly, the genotype-phenotype correlation suggested that the susceptible 'G' allele of NFATC2 rs4811198(T>G) &NFATC4 rs11848279(A>G) 3' UTR SNPs lead to reduced NFATC2 and NFATC4 expression(p<0.0001). Furthermore, in silico analysis suggested that hsa-miR-3183 & hsa-miR-6720-3p miRNAs specifically bound to 'G' allele of NFATC2 rs4811198 SNP and has-miR-4652-3p miRNA specifically bound to 'G' allele of NFATC4 rs11848279 SNP. Overall, our study suggests that NFATC2 rs4811198(T>G) 3' UTR &NFATC2 rs12479626(T>C) structural SNPs may be associated with GV susceptibility in Gujarat population. Moreover, the susceptible alleles for the 3' UTR SNPs could lead to reduced NFATs levels, which may further possibly, affect the Treg suppressive function leading to GV.

Identification and analysis of RNA structural disruptions induced by single nucleotide variants using Riprap and RiboSNitchDB

RNA conformational alteration has significant impacts on cellular processes and phenotypic variations. An emerging genetic factor of RNA conformational alteration is a new class of single nucleotide variant (SNV) named riboSNitch. RiboSNitches have been demonstrated to be involved in many genetic diseases. However, identifying riboSNitches is notably difficult as the signals of RNA structural disruption are often subtle. Here, we introduce a novel computational framework–RIboSNitch Predictor based on Robust Analysis of Pairing probabilities (Riprap). Riprap identifies structurally disrupted regions around any given SNVs based on robust analysis of local structural configurations between wild-type and mutant RNA sequences. Compared to previous approaches, Riprap shows higher accuracy when assessed on hundreds of known riboSNitches captured by various experimental RNA structure probing methods including the parallel analysis of RNA structure (PARS) and the selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE). Further, Riprap detects the experimentally validated riboSNitch that regulates human catechol-O-methyltransferase haplotypes and outputs structurally disrupted regions precisely at base resolution. Riprap provides a new approach to interpreting disease-related genetic variants. In addition, we construct a database (RiboSNitchDB) that includes the annotation and visualization of all presented riboSNitches in this study as well as 24 629 predicted riboSNitches from human expression quantitative trait loci.

Critical Analysis of Genome-Wide Association Studies: Triple Negative Breast Cancer Quae Exempli Causa

Genome-wide association studies (GWAS) are useful in assessing and analyzing either differences or variations in DNA sequences across the human genome to detect genetic risk factors of diseases prevalent within a target population under study. The ultimate goal of GWAS is to predict either disease risk or disease progression by identifying genetic risk factors. These risk factors will define the biological basis of disease susceptibility for the purposes of developing innovative, preventative, and therapeutic strategies. As single nucleotide polymorphisms (SNPs) are often used in GWAS, their relevance for triple negative breast cancer (TNBC) will be assessed in this review. Furthermore, as there are different levels and patterns of linkage disequilibrium (LD) present within different human subpopulations, a plausible strategy to evaluate known SNPs associated with incidence of breast cancer in ethnically different patient cohorts will be presented and discussed. Additionally, a description of GWAS for TNBC will be presented, involving various identified SNPs correlated with miRNA sites to determine their efficacies on either prognosis or progression of TNBC in patients. Although GWAS have identified multiple common breast cancer susceptibility variants that individually would result in minor risks, it is their combined effects that would likely result in major risks. Thus, one approach to quantify synergistic effects of such common variants is to utilize polygenic risk scores. Therefore, studies utilizing predictive risk scores (PRSs) based on known breast cancer susceptibility SNPs will be evaluated. Such PRSs are potentially useful in improving stratification for screening, particularly when combining family history, other risk factors, and risk prediction models. In conclusion, although interpretation of the results from GWAS remains a challenge, the use of SNPs associated with TNBC may elucidate and better contextualize these studies.

Blood-based analysis of type-2 diabetes mellitus susceptibility genes identifies specific transcript variants with deregulated expression and association with disease risk

Despite significant progress by genome-wide association studies, the ability of genetic variants to conduce to the prediction or prognosis of type-2 diabetes (T2D) is weak. Expression analysis of the corresponding genes may suggest possible links between single-nucleotide polymorphisms and T2D phenotype and/or risk. Herein, we investigated the expression patterns of 24 T2D-susceptibility genes, and their individual transcript variants (tv), in peripheral blood of T2D patients and controls (CTs), applying RNA-seq and real-time qPCR methodologies, and explore possible associations with disease features. Our data revealed the deregulation of certain transcripts in T2D patients. Among them, the down-regulation of CAPN10 tv3 was confirmed as an independent predictor for T2D. In patients, increased expression of CDK5 tv2, CDKN2A tv3 or THADA tv5 correlated positively with serum insulin levels, of CDK5 tv1 positively with % HbA1c levels, while in controls, elevated levels of TSPAN8 were associated positively with the presence of T2D family history. Herein, a T2D-specific expression profile of specific transcripts of disease-susceptibility genes is for the first time described in human peripheral blood. Large-scale studies are needed to evaluate the potential of these molecules to serve as disease biomarkers.

Noncoding RNAs Databases: Current Status and Trends

One of the most important resources for researchers of noncoding RNAs is the information available in public databases spread over the internet. However, the effective exploration of this data can represent a daunting task, given the large amount of databases available and the variety of stored data. This chapter describes a classification of databases based on information source, type of RNA, source organisms, data formats, and the mechanisms for information retrieval, detailing the relevance of each of these classifications and its usability by researchers. This classification is used to update a 2012 review, indexing now more than 229 public databases. This review will include an assessment of the new trends for ncRNA research based on the information that is being offered by the databases. Additionally, we will expand the previous analysis focusing on the usability and application of these databases in pathogen and disease research. Finally, this chapter will analyze how currently available database schemas can help the development of new and improved web resources.

regQTLs: Single nucleotide polymorphisms that modulate microRNA regulation of gene expression in tumors

Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with trait diversity and disease susceptibility, yet their functional properties often remain unclear. It has been hypothesized that SNPs in microRNA binding sites may disrupt gene regulation by microRNAs (miRNAs), short non-coding RNAs that bind to mRNA and downregulate the target gene. While several studies have predicted the location of SNPs in miRNA binding sites, to date there has been no comprehensive analysis of their impact on miRNA regulation. Here we investigate the functional properties of genetic variants and their effects on miRNA regulation of gene expression in cancer. Our analysis is motivated by the hypothesis that distinct alleles may cause differential binding (from miRNAs to mRNAs or from transcription factors to DNA) and change the expression of genes. We previously identified pathways—systems of genes conferring specific cell functions—that are dysregulated by miRNAs in cancer, by comparing miRNA–pathway associations between healthy and tumor tissue. We draw on these results as a starting point to assess whether SNPs on dysregulated pathways are responsible for miRNA dysregulation of individual genes in tumors. Using an integrative regression analysis that incorporates miRNA expression, mRNA expression, and SNP genotype data, we identify functional SNPs that we term “regulatory QTLs (regQTLs)”: loci whose alleles impact the regulation of genes by miRNAs. We apply the method to breast, liver, lung, and prostate cancer data from The Cancer Genome Atlas, and provide a tool to explore the findings.

Genomics studies have identified single nucleotide polymorphisms (SNPs) associated with trait diversity and disease susceptibility, yet the mechanism of action of many genetic variants remains unclear. MicroRNAs (miRNAs) are a class of small non-coding RNA molecules that base-pair coding mRNAs to regulate gene transcription. We hypothesize that SNP variants may affect the ability of miRNAs to bind their target genes, thus influencing gene regulation. To identify these “regulatory QTLs” (regQTLs), we integrate miRNA expression, mRNA expression, and SNP data to identify miRNAs that are associated with pathway dysregulation in tumors, and assess whether SNPs on these pathways are responsible for disrupted miRNA-gene regulation. This data-driven approach enables the discovery of SNPs whose alleles impact gene regulation by miRNAs, with functional consequences for tumor biology. We detail the method, apply it to data from The Cancer Genome Atlas, and provide a tool to explore the findings.

miRNA Mediated Noise Making of 3'UTR Mutations in Cancer

Somatic mutations in 3'-untranslated regions (3'UTR) do not alter amino acids and are considered to be silent in cancers. We found that such mutations can promote tumor progression by altering microRNA (miRNA) targeting efficiency and consequently affecting miRNA⁻mRNA interactions. We identified 67,159 somatic mutations located in the 3'UTRs of messenger RNAs (mRNAs) which can alter miRNA⁻mRNA interactions (functional somatic mutations, funcMutations), and 69.3% of these funcMutations (the degree of energy change > 12 kcal/mol) were identified to significantly promote loss of miRNA-mRNA binding. By integrating mRNA expression profiles of 21 cancer types, we found that the expression of target genes was positively correlated with the loss of absolute affinity level and negatively correlated with the gain of absolute affinity level. Functional enrichment analysis revealed that genes carrying funcMutations were significantly enriched in the MAPK and WNT signaling pathways, and analysis of regulatory modules identified eighteen miRNA modules involved with similar cellular functions. Our findings elucidate a complex relationship between miRNA, mRNA, and mutations, and suggest that 3'UTR mutations may play an important role in tumor development.

Genetic variants in mRNA untranslated regions

Genome Wide Association Studies (GWAS) have mapped thousands of genetic variants associated with complex disease risk and regulating quantitative traits, thus exploiting an unprecedented high-resolution genetic characterization of the human genome. A small fraction (3.7%) of the identified associations is located in untranslated regions (UTRs), and the molecular mechanism has been elucidated for few of them. Genetic variations at UTRs may modify regulatory elements affecting the interaction of the UTRs with proteins and microRNAs. The overall functional consequences include modulation of messenger RNA (mRNA) transcription, secondary structure, stability, localization, translation, and access to regulators like microRNAs (miRNAs) and RNA-binding proteins (RBPs). Alterations of these regulatory mechanisms are known to modify molecular pathways and cellular processes, potentially leading to disease processes. Here, we analyze some examples of genetic risk variants mapping in the UTR regulatory elements. We describe a recently identified genetic variant localized in the 3'UTR of the TNFSF13B gene, associated with autoimmunity risk and responsible of an increased stability and translation of TNFSF13B mRNA. We discuss how the correct use and interpretation of public GWAS repositories could lead to a better understanding of etiopathogenetic mechanisms and the generation of robust biological hypothesis as starting point for further functional studies. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Evolution and Genomics > Computational Analyses of RNA RNA in Disease and Development > RNA in Disease.

A GDF15 3' UTR variant, rs1054564, results in allele-specific translational repression of GDF15 by hsa-miR-1233-3p

Growth differentiation factor 15 (GDF15) is a strong predictor of cardiovascular events and mortality in individuals with or without cardiovascular diseases. Single nucleotide polymorphisms (SNPs) in microRNA (miRNA) target sites, also known as miRSNPs, are known to enhance or weaken miRNA-mRNA interactions and have been linked to diseases such as cardiovascular disease and cancer. In this study, we aimed to elucidate the functional significance of the miRSNP rs1054564 in regulating GDF15 levels. Two rs1054564-containing binding sites for hsa-miR-873-5p and hsa-miR-1233-3p were identified in the 3′ untranslated region (UTR) of the GDF15 transcript using bioinformatics tools. Their activities were further characterized by in vitro reporter assays. Bioinformatics prediction suggested that miRNA binding sites harboring the rs1054564-G allele had lower free energies than those with the C allele and therefore were better targets with higher affinities for both hsa-miR-873-5p and hsa-miR-1233-3p. Reporter assays showed that luciferase activity was significantly decreased by rs1054564-G-containing 3′ UTRs for both miRNAs (P < 0.05) and was restored by miRNA inhibitors. Comparing the fold suppression of the two miRNAs, only that of hsa-miR-1233-3p showed significant changes between the rs1054564-G- and C-containing 3′ UTRs (P = 0.034). In addition, western blots showed that transfection of both miRNA mimics significantly decreased endogenous GDF15 expression in a melanoma cell line (P < 0.05). Taken together, our findings demonstrate that GDF15 is a target of hsa-miR-873-5p and hsa-miR-1233-3p and that the rs1054564-C allele partially abolishes hsa-miR-1233-3p-mediated translational suppression of GDF15. These results suggest that rs1054564 confers allele-specific translational repression of GDF15 via hsa-miR-1233-3p. Our work thus provides biological insight into the previously reported clinical association between rs1054564 and plasma GDF15 levels.

An Integrating Approach for Genome-Wide Screening of MicroRNA Polymorphisms Mediated Drug Response Alterations

MicroRNAs (miRNAs) are a class of evolutionarily conserved small noncoding RNAs, ~22 nt in length, and found in diverse organisms and play important roles in the regulation of mRNA translation and degradation. It was shown that miRNAs were involved in many key biological processes through regulating the expression of targets. Genetic polymorphisms in miRNA target sites may alter miRNA regulation and therefore result in the alterations of the drug targets. Recent studies have demonstrated that SNPs in miRNA target sites can affect drug efficiency. However, there are still a large number of specific genetic variants related to drug efficiency that are yet to be discovered. We integrated large scale of genetic variations, drug targets, gene interaction networks, biological pathways, and seeds region of miRNA to identify miRNA polymorphisms affecting drug response. In addition, harnessing the abundant high quality biological network/pathways, we evaluated the cascade distribution of tarSNP impacts. We showed that the predictions can uncover most of the known experimentally supported cases as well as provide informative candidates complementary to existing methods/tools. Although there are several existing databases predicting the gain or loss of targeting function of miRNA mediated by SNPs, such as PolymiRTS, miRNASNP, MicroSNiPer, and MirSNP, none of them evaluated the influences of tarSNPs on drug response alterations. We developed a user-friendly online database of this approach named Mir2Drug.

Identification of a functional polymorphism affecting microRNA binding in the susceptibility locus 1q25.3 for colorectal cancer

Genome-wide association studies (GWASs) have identified dozens of susceptibility loci for colorectal cancer (CRC). However, most of them lack functional genetic variants and clear biological mechanisms. MicroRNAs (miRNAs) are small noncoding RNAs involved in a variety of physiological and tumorigenic processes. Here we hypothesized that single nucleotide polymorphisms (SNPs) that affect miRNAs biogenesis and binding, could contribute to CRC risk in the Chinese population. To locate miRNA-related SNPs in established GWAS loci, we initially screened out five candidate SNPs using a systematic bioinformatics analysis. Then, we performed a two-stage case-control study consisting of 2347 cases and 3390 controls, and found a positive polymorphism rs1062044, which presented consistently significant associations with CRC in both stages, and with an odds ratio (OR) = 1.32 (95% confidence interval (95%CI) = 1.18-1.49, P = 3.43E-06) under the dominant model in the combined study. Further luciferase reporter gene assays indicated that the variant G allele obviously improved the specific binding between miR-423-5p and the gene LAMC1. These findings suggested that the functional SNP rs1062044 at 1q25.3 might be a genetic modifier for the occurrence and development of CRC.

Personalized genetics of the cholinergic blockade of neuroinflammation

Acetylcholine signaling is essential for cognitive functioning and blocks inflammation. To maintain homeostasis, cholinergic signaling is subjected to multi‐leveled and bidirectional regulation by both proteins and non‐coding microRNAs (‘CholinomiRs’). CholinomiRs coordinate the cognitive and inflammatory aspects of cholinergic signaling by targeting major cholinergic transcripts including the acetylcholine hydrolyzing enzyme acetylcholinesterase (AChE). Notably, AChE inhibitors are the only currently approved line of treatment for Alzheimer's disease patients. Since cholinergic signaling blocks neuroinflammation which is inherent to Alzheimer's disease, genomic changes modifying AChE's properties and its susceptibility to inhibitors and/or to CholinomiRs regulation may affect the levels and properties of inflammasome components such as NLRP3. This calls for genomic‐based medicine approaches based on genotyping of both coding and non‐coding single nucleotide polymorphisms (SNPs) in the genes involved in cholinergic signaling. An example is a SNP in a recognition element for the primate‐specific microRNA‐608 within the 3′ untranslated region of the AChE transcript. Carriers of the minor allele of that SNP present massively elevated brain AChE levels, increased trait anxiety and inflammation, accompanied by perturbed CholinomiR‐608 regulatory networks and elevated prefrontal activity under exposure to stressful insults. Several additional SNPs in the AChE and other cholinergic genes await further studies, and might likewise involve different CholinomiRs and pathways including those modulating the initiation and progression of neurodegenerative diseases. CholinomiRs regulation of the cholinergic system thus merits in‐depth interrogation and is likely to lead to personalized medicine approaches for achieving better homeostasis in health and disease.

This is an article for the special issue XVth International Symposium on Cholinergic Mechanisms.

CRY1 and CRY2 genetic variants in seasonality: A longitudinal and cross-sectional study

Cryptochromes are key components of the circadian clocks that generate and maintain seasonal variations. The aim of our study was to analyze the associations of CRY1 and CRY2 genetic variants with the problematicity of seasonal variations, and whether the problematicity of seasonal variations changed during the follow-up of 11 years. Altogether 21 CRY1 and 16 CRY2 single-nucleotide polymorphisms (SNPs) were genotyped and analyzed in 5910 individuals from a Finnish nationwide population-based sample who had filled in the self-report on the seasonal variations in mood and behavior in the year 2000. In the year 2011, 3356 of these individuals filled in the same self-report on the seasonal variations in mood and behavior. Regression models were used to test whether any of the SNPs associated with the problematicity of seasonal variations or with a change in the problematicity from 2000 to 2011. In the longitudinal analysis, CRY2 SNP rs61884508 was protective from worsening of problematicity of seasonal variations. In the cross-sectional analysis, CRY2 SNP rs72902437 showed evidence of association with problematicity of seasonal variations, as did SNP rs1554338 (in the MAPK8IP1 and downstream of CRY2).

The regulatory roles of non-coding RNAs in nerve injury and regeneration

Non-coding RNAs (ncRNAs), especially microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), have attracted much attention since their regulatory roles in diverse cell processes were recognized. Emerging studies demonstrate that many ncRNAs are differentially expressed after injury to the nervous system, significantly affecting nerve regeneration. In this review, we compile the miRNAs and lncRNAs that have been reported to be dysregulated following a variety of central and peripheral nerve injuries, including acquired brain injury, spinal cord injury, and peripheral nerve injury. We also list investigations on how these miRNAs and lncRNAs exert the regulatory actions in neurodegenerative and neuroregenerative processes through different mechanisms involving their interaction with target coding genes. We believe that comprehension of the expression profiles and the possible functions of ncRNAs during the processes of nerve injury and regeneration will help understand the molecular mechanisms responsible for post-nerve-injury changes, and may contribute to the potential use of ncRNAs as a diagnostic marker and therapeutic target for nerve injury.

The Human Myotrophin Variant Attenuates MicroRNA-Let-7 Binding Ability but Not Risk of Left Ventricular Hypertrophy in Human Essential Hypertension

Myotrophin, known as a myocardial hypertrophy-inducing factor, is responsible for the initiation of cardiac hypertrophy that transits to heart failure. MicroRNAs are small noncoding RNAs that down-regulate posttranscriptional expression of target molecules. We investigated the role of variants of the microRNA-binding site in myotrophin in affecting its expression and any association with cardiac hypertrophy. Bioinformatics demonstrated that variant rs17168525 was identified to be located in the let-7/miR-98-binding site of myotrophin. We further experimentally test to effects of the identified variant on myotrophin translation using luciferase reporter assay and Western blotting. We found that the C allele of rs17168525 suppressed myotrophin translation by facilitating let-7c binding, but not the T allele. Let-7c overexpression caused a significant decrease in the level of myotrophin protein. Next, we investigated the association of the variant with cardiac hypertrophy in 1614 hypertensive patients, including 552 with left ventricular hypertrophy and 1062 without left ventricular hypertrophy, as well as 591 healthy control subjects from a Han Chinese population. No significant association between the variant rs17168525 and left ventricular hypertrophy in hypertensive patients in a Han Chinese population (P>0.05). In conclusion, our experimental results provide evidence that the T allele of rs17168525 in the 3'-UTR of myotrophin might influence the level of myotrophin protein by interfering with let-7/miR-98 binding.

An update of miRNASNP database for better SNP selection by GWAS data, miRNA expression and online tools

MicroRNAs (miRNAs) are key regulators of gene expression involved in a broad range of biological processes. MiRNASNP aims to provide single nucleotide polymorphisms (SNPs) in miRNAs and genes that may impact miRNA biogenesis and/or miRNA target binding. Advanced miRNA research provided abundant data about miRNA expression, validated targets and related phenotypic variants. In miRNASNP v2.0, we have updated our previous database with several new data and features, including: (i) expression level and expression correlation of miRNAs and target genes in different tissues, (ii) linking SNPs to the results of genome-wide association studies, (iii) integrating experimentally validated miRNA:mRNA interactions, (iv) adding multiple filters to prioritize functional SNPs. In addition, as a supplement of the database, we have set up three flexible online tools to analyse the influence of novel variants on miRNA:mRNA binding. A new nice web interface was designed for miRNASNP v2.0 allowing users to browse, search and download. We aim to maintain the miRNASNP as a solid resource for function, genetics and disease studies of miRNA-related SNPs. Database URL: http://bioinfo.life. hust.edu.cn/miRNASNP2/.

Genetics of type 2 diabetes-pitfalls and possibilities

Type 2 diabetes (T2D) is a complex disease that is caused by a complex interplay between genetic, epigenetic and environmental factors. While the major environmental factors, diet and activity level, are well known, identification of the genetic factors has been a challenge. However, recent years have seen an explosion of genetic variants in risk and protection of T2D due to the technical development that has allowed genome-wide association studies and next-generation sequencing. Today, more than 120 variants have been convincingly replicated for association with T2D and many more with diabetes-related traits. Still, these variants only explain a small proportion of the total heritability of T2D. In this review, we address the possibilities to elucidate the genetic landscape of T2D as well as discuss pitfalls with current strategies to identify the elusive unknown heritability including the possibility that our definition of diabetes and its subgroups is imprecise and thereby makes the identification of genetic causes difficult.

Mature miRNAs form secondary structure, which suggests their function beyond RISC

The generally accepted model of the miRNA-guided RNA down-regulation suggests that mature miRNA targets mRNA in a nucleotide sequence-specific manner. However, we have shown that the nucleotide sequence of miRNA is not the only determinant of miRNA specificity. Using specific nucleases, T1, V1 and S1 as well as NMR, UV/Vis and CD spectroscopies, we found that miR-21, miR-93 and miR-296 can adopt hairpin and/or homoduplex structures. The secondary structure of those miRNAs in solution is a function of RNA concentration and ionic conditions. Additionally, we have shown that a formation of miRNA hairpin is facilitated by cellular environment.Looking for functional consequences of this observation, we have perceived that structure of these miRNAs resemble RNA aptamers, short oligonucleotides forming a stable 3D structures with a high affinity and specificity for their targets. We compared structures of anti-tenascin C (anti-Tn-C) aptamers, which inhibit brain tumor glioblastoma multiforme (GBM, WHO IV) and selected miRNA. A strong overexpression of miR-21, miR-93 as well Tn-C in GBM may imply some connections between them. The structural similarity of these miRNA hairpins and anti-Tn-C aptamers indicates that miRNAs may function also beyond RISC and are even more sophisticated regulators, that it was previously expected. We think that the knowledge of the miRNA structure may give a new insight into miRNA-dependent gene regulation mechanism and be a step forward in the understanding their function and involvement in cancerogenesis. This may improve design process of anti-miRNA therapeutics.

CRY1, CRY2 and PRKCDBP genetic variants in metabolic syndrome

The circadian clock affects metabolic cycles, and there is a link between circadian clock genes and metabolic syndrome. Therefore, we wanted to investigate whether variants of the core circadian clock genes, cryptochrome circadian clocks 1 and 2 (CRY1 and CRY2), or those of protein kinase C, delta binding protein (PRKCDBP), which regulate the interactions and abundance of dimers of the period and cryptochrome proteins, are associated with metabolic syndrome or its components. The association of 48 single-nucleotide polymorphisms (SNPs) from CRY1, CRY2 and PRKCDBP genes with metabolic disorder or its components was analyzed in a sample of 5910 individuals. Genotyping was performed using the Sequenom MassARRAY system. SNPs and haplotypes were analyzed using linear or logistic regression with additive models controlling for age and sex. Continuous phenotypes were permuted 10,000 times. False discovery rate q-values were calculated to correct for multiple testing. Overall, CRY1 and CRY2 variants showed nominal association with the metabolic syndrome components, hypertension and triglyceride levels, and one CRY2 variant had an association with metabolic syndrome, although none of these associations yielded significant q-values. However, the haplotype analysis of these variants supported the association of CRY1 with arterial hypertension and elevated blood pressure. Further studies are warranted regarding the role of CRY1 in cardiovascular diseases.

CardioGxE, a catalog of gene-environment interactions for cardiometabolic traits

Background

Genetic understanding of complex traits has developed immensely over the past decade but remains hampered by incomplete descriptions of contribution to phenotypic variance. Gene-environment (GxE) interactions are one of these contributors and in the guise of diet and physical activity are important modulators of cardiometabolic phenotypes and ensuing diseases.

Results

We mined the scientific literature to collect GxE interactions from 386 publications for blood lipids, glycemic traits, obesity anthropometrics, vascular measures, inflammation and metabolic syndrome, and introduce CardioGxE, a gene-environment interaction resource. We then analyzed the genes and SNPs supporting cardiometabolic GxEs in order to demonstrate utility of GxE SNPs and to discern characteristics of these important genetic variants. We were able to draw many observations from our extensive analysis of GxEs. 1) The CardioGxE SNPs showed little overlap with variants identified by main effect GWAS, indicating the importance of environmental interactions with genetic factors on cardiometabolic traits. 2) These GxE SNPs were enriched in adaptation to climatic and geographical features, with implications on energy homeostasis and response to physical activity. 3) Comparison to gene networks responding to plasma cholesterol-lowering or regression of atherosclerotic plaques showed that GxE genes have a greater role in those responses, particularly through high-energy diets and fat intake, than do GWAS-identified genes for the same traits. Other aspects of the CardioGxE dataset were explored.

Conclusions

Overall, we demonstrate that SNPs supporting cardiometabolic GxE interactions often exhibit transcriptional effects or are under positive selection. Still, not all such SNPs can be assigned potential functional or regulatory roles often because data are lacking in specific cell types or from treatments that approximate the environmental factor of the GxE. With research on metabolic related complex disease risk embarking on genome-wide GxE interaction tests, CardioGxE will be a useful resource.

Current trend of annotating single nucleotide variation in humans--A case study on SNVrap

As high throughput methods, such as whole genome genotyping arrays, whole exome sequencing (WES) and whole genome sequencing (WGS), have detected huge amounts of genetic variants associated with human diseases, function annotation of these variants is an indispensable step in understanding disease etiology. Large-scale functional genomics projects, such as The ENCODE Project and Roadmap Epigenomics Project, provide genome-wide profiling of functional elements across different human cell types and tissues. With the urgent demands for identification of disease-causal variants, comprehensive and easy-to-use annotation tool is highly in demand. Here we review and discuss current progress and trend of the variant annotation field. Furthermore, we introduce a comprehensive web portal for annotating human genetic variants. We use gene-based features and the latest functional genomics datasets to annotate single nucleotide variation (SNVs) in human, at whole genome scale. We further apply several function prediction algorithms to annotate SNVs that might affect different biological processes, including transcriptional gene regulation, alternative splicing, post-transcriptional regulation, translation and post-translational modifications. The SNVrap web portal is freely available at http://jjwanglab.org/snvrap.

Genetic variability of microRNA regulome in human

MicroRNAs are currently being extensively studied due to their important role as post-transcriptional regulators. During miRNA biogenesis, precursors undergo two cleavage steps performed by Drosha-DGCR8 (Microprocessor) cleaving of pri-miRNA to produce pre-miRNA and Dicer-mediated cleaving to create mature miRNA. Genetic variants within human miRNA regulome have been shown to influence miRNA expression, target interaction and to affect the phenotype. In this study, we reviewed the literature, existing bioinformatics tools and catalogs associated with polymorphic miRNA regulome, and organized them into four categories: (1) polymorphisms located within miRNA genes (miR-SNPs), (2) transcription factor-binding sites/miRNA regulatory regions (miR-rSNPs), (3) miRNA target sites (miR-TS-SNPs), and 4. miRNA silencing machinery (miR-SM-SNPs). Since the miR-SM-SNPs have not been systematically studied yet, we have collected polymorphisms associated with miRNA silencing machinery. We have developed two catalogs containing genetic variability within: (1) genes encoding three main catalytic components of the silencing machinery, DROSHA, DGCR8, and DICER1; (2) miRNA genes itself, overlapping Drosha and Dicer cleavage sites. The developed resource of polymorphisms is available online (http://www.integratomics-time.com/miRNA-regulome) and will be useful for further functional studies and development of biomarkers associated with diseases and phenotypic traits.

Identifying a polymorphic 'switch' that influences miRNAs' regulation of a myasthenia gravis risk pathway

The significant roles of genetic variants in myasthenia gravis (MG) pathogenesis have been demonstrated in many studies, and recently it has been revealed that aberrant level/regulation of microRNAs (miRNAs) might contribute to the initiation and progression of MG. However, the dysfunction of miRNA associated with single nucleotide polymorphisms (miRSNPs) has not been well investigated in MG. In this study, we created a contemporary catalog of 89 MG risk genes via manual literature-mining. Based on this risk gene catalog, we obtained 18 MG risk pathways. Furthermore, we identified 93 miRNAs that target MG risk pathways and revealed the miRSNPs ‘switches’ in miRNA regulation in the MG risk pathways by integrating the database information of miRSNPs. We also constructed a miRNA-mediated SNP switching pathway network (MSSPN) to intuitively analyze miRNA regulation of MG risk pathways and the relationship of the polymorphism ‘switch’ with these changes in regulation. Moreover, we carried out in-depth dissection on the correlation between hsa05200 (pathway in cancer) and MG development, and elaborated the significance of 4 high-risk genes. By network analysis and literature mining, we proposed a potential mechanism of miRSNPs→gene→pathway effects on MG pathogenesis, especially for rs28457673 (miR-15/16/195/424/497 family)→IGF1R→hsa05200 (pathway in cancer). Therefore, our studies have revealed a functional role for genetic modulators in MG pathogenesis at a systemic level, which could be informative for further miRNA and miRSNPs studies in MG.

Genetics of diabetes--are we missing the genes or the disease?

Diabetes is a group of metabolic diseases characterized by hyperglycemia resulting from defects in insulin secretion, insulin action, or both. The chronic hyperglycemia of diabetes is associated with long-term damage, dysfunction, and failure of different organs, especially the eyes, kidneys, nerves, heart, and blood vessels. Several pathogenic processes are involved in the development of diabetes. These range from autoimmune destruction of the beta-cells of the pancreas with consequent insulin deficiency to abnormalities that result in resistance to insulin action (American Diabetes Association, 2011). The vast majority of cases of diabetes fall into two broad categories. In type 1 diabetes (T1D), the cause is an absolute deficiency of insulin secretion, whereas in type 2 diabetes (T2D), the cause is a combination of resistance to insulin action and an inadequate compensatory insulin secretory response. However, the subdivision into two main categories represents a simplification of the real situation, and research during the recent years has shown that the disease is much more heterogeneous than a simple subdivision into two major subtypes assumes. Worldwide prevalence figures estimate that there are 280 million diabetic patients in 2011 and more than 500 million in 2030 (http://www.diabetesatlas.org/). In Europe, about 6-8% of the population suffer from diabetes, of them about 90% has T2D and 10% T1D, thereby making T2D to the fastest increasing disease in Europe and worldwide. This epidemic has been ascribed to a collision between the genes and the environment. While our knowledge about the genes is clearly better for T1D than for T2D given the strong contribution of variation in the HLA region to the risk of T1D, the opposite is the case for T2D, where our knowledge about the environmental triggers (obesity, lack of exercise) is much better than the understanding of the underlying genetic causes. This lack of knowledge about the underlying genetic causes of diabetes is often referred to as missing heritability (Manolio et al., 2009) which exceeds 80% for T2D but less than 25% for T1D. In the following review, we will discuss potential sources of this missing heritability which also includes the possibility that our definition of diabetes and its subgroups is imprecise and thereby making the identification of genetic causes difficult.

MicroRNA and noncoding RNA-related data sources

Noncoding RNAs (ncRNAs) are ribonucleic acids capable of controlling different genetic and metabolic functions. These molecules have been recently organized into different classes, and among them microRNAs (miRNAs) are extensively studied. MicroRNAs are short oligomers mainly involved in posttranscriptional gene silencing. The specific research field, focused on structural and functional characterization of microRNAs, is commonly called mirnomics. The exploitation of the interest in microRNAs has stimulated the organization of several databases that are often integrated with analytical tools in order to predict microRNA targets, or to find those miRNAs capable to inhibit the expression of a specific protein. This work attempts to provide an overview of accessible information about microRNAs and other noncoding RNAs that has been gathered in curated databases.

Identification of microRNAs regulated by isothiocyanates and association of polymorphisms inside their target sites with risk of sporadic colorectal cancer

Sporadic colorectal cancer (CRC) is a typical multifactorial disease. Isothiocyanates (ITC) have been recently shown to inhibit development of CRC in many experimental models. MicroRNAs (miRNAs) are short noncoding RNAs that posttranscriptionally regulate gene expression through binding to 3' untranslated regions (3'UTR) of target mRNAs. MiRNAs are regulated by natural agents, ITCs included. In our study, using global expression profiling based on TaqMan Low-Density Arrays, we identified 3 common miRNAs (miR-155, miR-23b, miR-27b) regulated by ITCs (sulforaphane, iberin) in colonic epithelial cell lines NCM460 and NCM356. In silico predictions allowed us to find 9 relevant single nucleotide polymorphisms (SNPs) localized within the 3'UTRs of genes (AGTR1, TNFAIP2, PRKCB, HSPA9, RABGAP1, DICER1, ADAM19, VWA5A, and SIRT5) targeted by these ITC-related miRNAs. Finally, we observed that homozygous CC genotype of DICER1, rs1057035, was significantly associated with decreased risk of CRC (odds ratio = 0.49; 95% confidence interval: 0.25-0.95, P = 0.036) when compared to TT homozygote genotype; also, the C allele tended to have a protective effect (P = 0.072). This study showed that miRNAs could be involved in chemoprotective effects of natural agents; their function alteration through SNPs in their binding sites and flanking regions presents a new class of CRC risk factors.

RNAsnp: efficient detection of local RNA secondary structure changes induced by SNPs

Structural characteristics are essential for the functioning of many noncoding RNAs and cis-regulatory elements of mRNAs. SNPs may disrupt these structures, interfere with their molecular function, and hence cause a phenotypic effect. RNA folding algorithms can provide detailed insights into structural effects of SNPs. The global measures employed so far suffer from limited accuracy of folding programs on large RNAs and are computationally too demanding for genome-wide applications. Here, we present a strategy that focuses on the local regions of maximal structural change between mutant and wild-type. These local regions are approximated in a “screening mode” that is intended for genome-wide applications. Furthermore, localized regions are identified as those with maximal discrepancy. The mutation effects are quantified in terms of empirical P values. To this end, the RNAsnp software uses extensive precomputed tables of the distribution of SNP effects as function of length and GC content. RNAsnp thus achieves both a noise reduction and speed-up of several orders of magnitude over shuffling-based approaches. On a data set comprising 501 SNPs associated with human-inherited diseases, we predict 54 to have significant local structural effect in the untranslated region of mRNAs. RNAsnp is available at http://rth.dk/resources/rnasnp.

CRY2 genetic variants associate with dysthymia

People with mood disorders often have disruptions in their circadian rhythms. Recent molecular genetics has linked circadian clock genes to mood disorders. Our objective was to study two core circadian clock genes, CRY1 and CRY2 as well as TTC1 that interacts with CRY2, in relation to depressive and anxiety disorders. Of these three genes, 48 single-nucleotide polymorphisms (SNPs) whose selection was based on the linkage disequilibrium and potential functionality were genotyped in 5910 individuals from a nationwide population-based sample. The diagnoses of major depressive disorder, dysthymia and anxiety disorders were assessed with a structured interview (M-CIDI). In addition, the participants filled in self-report questionnaires on depressive and anxiety symptoms. Logistic and linear regression models were used to analyze the associations of the SNPs with the phenotypes. Four CRY2 genetic variants (rs10838524, rs7121611, rs7945565, rs1401419) associated significantly with dysthymia (false discovery rate q<0.05). This finding together with earlier CRY2 associations with winter depression and with bipolar type 1 disorder supports the view that CRY2 gene has a role in mood disorders.

Prioritization of genetic variants in the microRNA regulome as functional candidates in genome-wide association studies

Comprehensive analyses of results from genome-wide association studies (GWAS) have demonstrated that complex disease/trait-associated loci are enriched in gene regulatory regions of the genome. The search for causal regulatory variation has focused primarily on transcriptional elements, such as promoters and enhancers. microRNAs (miRNAs) are now widely appreciated as critical posttranscriptional regulators of gene expression and are thought to impart stability to biological systems. Naturally occurring genetic variation in the miRNA regulome is likely an important contributor to phenotypic variation in the human population. However, the extent to which polymorphic miRNA-mediated gene regulation underlies GWAS signals remains unclear. In this study, we have developed the most comprehensive bioinformatic analysis pipeline to date for cataloging and prioritizing variants in the miRNA regulome as functional candidates in GWAS. We highlight specific findings, including a variant in the promoter of the miRNA let-7 that may contribute to human height variation. We also provide a discussion of how our approach can be expanded in the future. Overall, we believe that the results of this study will be valuable for researchers interested in determining whether GWAS signals implicate the miRNA regulome in their disease/trait of interest.

Emerging roles for miRNAs in the development, diagnosis, and treatment of diabetic nephropathy

Although the causes of diabetic nephropathy are not yet fully known, emerging evidence suggests a role for epigenetic factors in the development of the disease. In particular, microRNAs (miRNAs) are becoming recognized as important mediators of biological processes relevant to diabetic nephropathy. Until recently, investigations of miRNAs in the development of diabetic nephropathy have remained relatively limited; however, the number of reports identifying potential new candidates and mechanisms of impact is presently expanding at a rapid pace. This review seeks to summarize these recent findings, focusing on new candidates and/or novel mechanisms, including the intersection between genetic variation and miRNA function in modulating disease expression, emerging in the field. We also review the latest advances in the diagnostic and therapeutic potential of miRNAs in the treatment of diabetic nephropathy.

A fast and accurate SNP detection algorithm for next-generation sequencing data

Various methods have been developed for calling single-nucleotide polymorphisms from next-generation sequencing data. However, for satisfactory performance, most of these methods require expensive high-depth sequencing. Here, we propose a fast and accurate single-nucleotide polymorphism detection program that uses a binomial distribution-based algorithm and a mutation probability. We extensively assess this program on normal and cancer next-generation sequencing data from The Cancer Genome Atlas project and pooled data from the 1,000 Genomes Project. We also compare the performance of several state-of-the-art programs for single-nucleotide polymorphism calling and evaluate their pros and cons. We demonstrate that our program is a fast and highly accurate single-nucleotide polymorphism detection method, particularly when the sequence depth is low. The program can finish single-nucleotide polymorphism calling within four hours for 10-fold human genome next-generation sequencing data (30 gigabases) on a standard desktop computer.

Single Nucleotide Polymorphisms Associated with MicroRNA Regulation

Since the discovery of microRNA (miRNA), the polymorphisms that affect miRNA regulation had been extensively investigated by many independent studies. Recently, researchers utilized bioinformatics and statistical approaches for genome-wide analysis on the human polymorphisms that reside in the miRNA genes, targets, and/or genes involved in miRNA processing. In this review, we will give an overview about the important findings of these studies from three perspectives: architecture of the polymorphisms within miRNAs or their targets, potential functional consequences of the polymorphisms on miRNA processing or targeting, and the associations of the polymorphisms with miRNA or target gene expression. The results of the previous studies demonstrated the signatures of natural selections on the miRNA genes and their targets, and proposed a collection of potentially functional, expression-associated, and/or positively selected polymorphisms that are promising for further investigations. In the meantime, a few useful resources about the polymorphic miRNA regulation have been developed and the different features of these databases were discussed in this review. Though recent research had benefited from these comprehensive studies and resources, there were still gaps in our knowledge about the polymorphisms involved in miRNA regulation, and future investigations were expected to address these questions.

A global map for dissecting phenotypic variants in human lincRNAs

Large intergenic noncoding RNAs (lincRNAs) are emerging as key factors of multiple cellular processes. Cumulative evidence has linked lincRNA polymorphisms to diverse diseases. However, the global properties of lincRNA polymorphisms and their implications for human disease remain largely unknown. Here we performed a systematic analysis of naturally occurring variants in human lincRNAs, with a particular focus on lincRNA polymorphism as novel risk factor of disease etiology. We found that lincRNAs exhibited a relatively low level of polymorphisms, and low single-nucleotide polymorphism (SNP) density lincRNAs might have a broad range of functions. We also found that some polymorphisms in evolutionarily conserved regions of lincRNAs had significant effects on predicted RNA secondary structures, indicating their potential contribution to diseases. We mapped currently available phenotype-associated SNPs to lincRNAs and found that lincRNAs were associated with a wide range of human diseases. Some lincRNAs could be responsible for particular diseases. Our results provided not only a global perspective on genetic variants in human lincRNAs but also novel insights into the function and etiology of lincRNA. All the data in this study can be accessed and retrieved freely via a web server at http://bioinfo.hrbmu.edu.cn/lincPoly.

Genome-wide analysis of human SNPs at long intergenic noncoding RNAs

Long intergenic noncoding RNAs (lincRNAs) represent a large portion of the noncoding genes in mammals and other eukaryotes but remains among the least well-understood of genetic factors to date. Here, we systematically analyzed the human SNPs of lincRNAs at a genome level. We found a significantly lower SNP density in lincRNA regions than both their upstream and downstream flanking regions. Functional regions show lower SNP density than other regions in lincRNAs. We revealed that lincRNAs with higher expression levels and broader expression spectrum have significantly lower SNP density. Moreover, we identified lincRNAs that are under recent positive selection and revealed that these lincRNAs show distinct SNP density, expression level, and tissue specificity. Importantly, we identified a genetic variant (rs7990916:T>C) under recent positive selection at a brain-specific lincRNA that significantly affects the structure of normal brain. Analysis of brain magnetic resonance images showed that individuals with CC genotype have significant bigger regional gray matter volume than individuals with TT genotype. Moreover, the genotype of this SNP shows different distribution in normal elders, mild cognitive impairment, and Alzheimer disease subjects, suggesting that this lincRNA may have a role in physiology and pathophysiology of human brain.

miRNA regulatory variation in human evolution

Recent advancements have revealed a complex post-transcriptional regulatory network in humans involving miRNAs. However, the contribution of miRNAs to human evolution, especially interindividual variation associated with miRNAs, is only beginning to be studied. In this article, we illustrate the extent of variation in miRNA-mediated post-transcriptional regulation in humans. Based on evidence from recent studies, we argue that the evolution of post-transcriptional control may be adaptive, and that it not only complements the primary transcriptional regulation by transcription factors (TFs), but also diversifies gene expression phenotypes, thereby generating genetic novelty on which natural selection subsequently acts. Given that current evolutionary analyses and genotype-phenotype mapping are primarily focused on protein-coding genes and TF-mediated regulations, comprehensive examination of post-transcriptional variations should be included in future studies to add a new dimension to understanding of human phenotypic evolution.

Genetic variations in key microRNA processing genes and risk of head and neck cancer: a case-control study in Chinese population

MicroRNAs (miRNAs) have been reported to play a key role in oncogenesis. Genetic variations in miRNA processing genes and miRNA binding sites may affect the biogenesis of miRNA and the miRNA-mRNA interactions, hence promoting tumorigenesis. In the present study, we hypothesized that potentially functional polymorphisms in miRNA processing genes may contribute to head and neck cancer (HNC) susceptibility. To test this hypothesis, we genotyped three SNPs at miRNA binding sites of miRNA processing genes (rs1057035 in 3'UTR of DICER, rs3803012 in 3'UTR of RAN and rs10773771 in 3'UTR of HIWI) with a case-control study including 397 HNC cases and 900 controls matched by age and sex in Chinese. Although none of three SNPs was significantly associated with overall risk of HNC, rs1057035 in 3'UTR of DICER was associated with a significantly decreased risk of oral cancer (TC/CC vs. TT: adjusted OR  =  0.65, 95% CI  =  0.46-0.92). Furthermore, luciferase activity assay showed that rs1057035 variant C allele led to significantly lower expression levels as compared to the T allele, which may be due to the relatively high inhibition of hsa-miR-574-3p on DICER mRNA. These findings indicated that rs1057035 located at 3'UTR of DICER may contribute to the risk of oral cancer by affecting the binding of miRNAs to DICER. Large-scale and well-designed studies are warranted to validate our findings.

Single nucleotide polymorphisms in genes that are common targets of luteotropin and luteolysin in primate corpus luteum: computational exploration

Luteal insufficiency affects fertility and hence study of mechanisms that regulate corpus luteum (CL) function is of prime importance to overcome infertility problems. Exploration of human genome sequence has helped to study the frequency of single nucleotide polymorphisms (SNPs). Clinical benefits of screening SNPs in infertility are being recognized well in recent times. Examining SNPs in genes associated with maintenance and regression of CL may help to understand unexplained luteal insufficiency and related infertility. Publicly available microarray gene expression databases reveal the global gene expression patterns in primate CL during the different functional state. We intend to explore computationally the deleterious SNPs of human genes reported to be common targets of luteolysin and luteotropin in primate CL. Different computational algorithms were used to dissect out the functional significance of SNPs in the luteinizing hormone sensitive genes. The results raise the possibility that screening for SNPs might be integrated to evaluate luteal insufficiency associated with human female infertility for future studies.

Mutations and binding sites of human transcription factors

Mutations in any genome may lead to phenotype characteristics that determine ability of an individual to cope with adaptation to environmental challenges. In studies of human biology, among the most interesting ones are phenotype characteristics that determine responses to drug treatments, response to infections, or predisposition to specific inherited diseases. Most of the research in this field has been focused on the studies of mutation effects on the final gene products, peptides, and their alterations. Considerably less attention was given to the mutations that may affect regulatory mechanism(s) of gene expression, although these may also affect the phenotype characteristics. In this study we make a pilot analysis of mutations observed in the regulatory regions of 24,667 human RefSeq genes. Our study reveals that out of eight studied mutation types, "insertions" are the only one that in a statistically significant manner alters predicted transcription factor binding sites (TFBSs). We also find that 25 families of TFBSs have been altered by mutations in a statistically significant manner in the promoter regions we considered. Moreover, we find that the related transcription factors are, for example, prominent in processes related to intracellular signaling; cell fate; morphogenesis of organs and epithelium; development of urogenital system, epithelium, and tube; neuron fate commitment. Our study highlights the significance of studying mutations within the genes regulatory regions and opens way for further detailed investigations on this topic, particularly on the downstream affected pathways.

Modeling SNP mediated differential targeting of homologous 3'UTR by microRNA

We had previously proposed that the post-transcriptional regulation through microRNA as a mechanism for incomplete penetrance and variable expressivity, leads to lack of correlation between genotype and phenotype. Here we report the validation of miRNA-target interactions we predicted earlier and demonstrate the regulation of endogenous JAG1 by hsa-miR-214 and hsa-miR-124, and TGFBR2 by hsa-miR-34b*, through luciferase activity of reporter constructs and also the expression levels of the endogenous genes. Using these targets, we have modeled the diploid state for miRNA target site with heterozygosity for the SNP and demonstrate the differential targeting of an otherwise identical 3'UTR. We show that SNP rs8708 (A > G) at the target site of hsa-miR-214 can relieve the repression while an SNP rs11466532 (C > T) enhances the repression of reporter expression by hsa-miR-34b*. We discuss the results in the light of its implications in the context of penetrance of dominant mutations in miRNA targeted genes, using JAG1 as an example. These observations imply that disease causing mutations in JAG1 linked to the SNP rs8708G will be poorly targeted by hsa-miR-214 when present against a normal allele of JAG1 with rs8708A and will show penetrance of JAG1 mutations as Alagille syndrome, while mutant JAG1 linked to rs8708A against rs8708G on the normal allele will show either no disease or much attenuated symptoms and hence exhibit incomplete penetrance.

A miR-1207-5p binding site polymorphism abolishes regulation of HBEGF and is associated with disease severity in CFHR5 nephropathy

Heparin binding epidermal growth factor (HBEGF) is expressed in podocytes and was shown to play a role in glomerular physiology. MicroRNA binding sites on the 3′UTR of HBEGF were predicted using miRWalk algorithm and followed by DNA sequencing in 103 patients diagnosed with mild or severe glomerulopathy. A single nucleotide polymorphism, miRSNP C1936T (rs13385), was identified at the 3′UTR of HBEGF that corresponds to the second base of the hsa-miR-1207-5p seed region. When AB8/13 undifferentiated podocytes were transfected with miRNA mimics of hsa-miR-1207-5p, the HBEGF protein levels were reduced by about 50%. A DNA fragment containing the miRSNP allele-1936C was cloned into the pMIR-Report Luciferase vector and co-transfected with miRNA mimics of hsa-miR-1207-5p into AB8/13 podocytes. In agreement with western blot data, this resulted in reduced luciferase expression demonstrating the ability of hsa-miR-1207-5p to directly regulate HBEGF expression. On the contrary, in the presence of the miRSNP 1936T allele, this regulation was abolished. Collectively, these results demonstrate that variant 1936T of this miRSNP prevents hsa-miR-1207-5p from down-regulating HBEGF in podocytes. We hypothesized that this variant has a functional role as a genetic modifier. To this end, we showed that in a cohort of 78 patients diagnosed with CFHR5 nephropathy (also known as C3-glomerulopathy), inheritance of miRSNP 1936T allele was significantly increased in the group demonstrating progression to chronic renal failure on long follow-up. No similar association was detected in a cohort of patients with thin basement membrane nephropathy. This is the first report associating a miRSNP as genetic modifier to a monogenic renal disorder.

Genome-wide identification of SNPs in microRNA genes and the SNP effects on microRNA target binding and biogenesis

MicroRNAs (miRNAs) are studied as key regulators of gene expression involved in different diseases. Several single nucleotide polymorphisms (SNPs) in miRNA genes or target sites (miRNA-related SNPs) have been proved to be associated with human diseases by affecting the miRNA-mediated regulatory function. To systematically analyze miRNA-related SNPs and their effects, we performed a genome-wide scan for SNPs in human pre-miRNAs, miRNA flanking regions, target sites, and designed a pipeline to predict the effects of them on miRNA-target interaction. As a result, we identified 48 SNPs in human miRNA seed regions and thousands of SNPs in 3' untranslated regions with the potential to either disturb or create miRNA-target interactions. Furthermore, we experimentally confirmed seven loss-of-function SNPs and one gain-of-function SNP by luciferase assay. This is the first case of experimental validation of an SNP in an miRNA creating a novel miRNA target binding. All useful data were complied into miRNASNP, a user-friendly free online database (http://www.bioguo.org/miRNASNP/). These data will be a useful resource for studying miRNA function, identifying disease-associated miRNAs, and further personalized medicine.

Cataloguing functionally relevant polymorphisms in gene DNA ligase I: a computational approach

A computational approach for identifying functionally relevant SNPs in gene LIG1 has been proposed. LIG1 is a crucial gene which is involved in excision repair pathways and mutations in this gene may lead to increase sensitivity towards DNA damaging agents. A total of 792 SNPs were reported to be associated with gene LIG1 in dbSNP. Different web server namely SIFT, PolyPhen, CUPSAT, FASTSNP, MAPPER and dbSMR were used to identify potentially functional SNPs in gene LIG1. SIFT, PolyPhen and CUPSAT servers predicted eleven nsSNPs to be intolerant, thirteen nsSNP to be damaging and two nsSNPs have the potential to destabilize protein structure. The nsSNP rs11666150 was predicted to be damaging by all three servers and its mutant structure showed significant increase in overall energy. FASTSNP predicted twenty SNPs to be present in splicing modifier binding sites while rSNP module from MAPPER server predicted nine SNPs to influence the binding of transcription factors. The results from the study may provide vital clues in establishing affect of polymorphism on phenotype and in elucidating drug response.

dPORE-miRNA: polymorphic regulation of microRNA genes

BACKGROUND

MicroRNAs (miRNAs) are short non-coding RNA molecules that act as post-transcriptional regulators and affect the regulation of protein-coding genes. Mostly transcribed by PolII, miRNA genes are regulated at the transcriptional level similarly to protein-coding genes. In this study we focus on human miRNAs. These miRNAs are involved in a variety of pathways and can affect many diseases. Our interest is on possible deregulation of the transcription initiation of the miRNA encoding genes, which is facilitated by variations in the genomic sequence of transcriptional control regions (promoters).

METHODOLOGY

Our aim is to provide an online resource to facilitate the investigation of the potential effects of single nucleotide polymorphisms (SNPs) on miRNA gene regulation. We analyzed SNPs overlapped with predicted transcription factor binding sites (TFBSs) in promoters of miRNA genes. We also accounted for the creation of novel TFBSs due to polymorphisms not present in the reference genome. The resulting changes in the original TFBSs and potential creation of new TFBSs were incorporated into the Dragon Database of Polymorphic Regulation of miRNA genes (dPORE-miRNA).

CONCLUSIONS

The dPORE-miRNA database enables researchers to explore potential effects of SNPs on the regulation of miRNAs. dPORE-miRNA can be interrogated with regards to: a/miRNAs (their targets, or involvement in diseases, or biological pathways), b/SNPs, or c/transcription factors. dPORE-miRNA can be accessed at http://cbrc.kaust.edu.sa/dpore and http://apps.sanbi.ac.za/dpore/. Its use is free for academic and non-profit users.

Let-7/miR-98 regulate Fas and Fas-mediated apoptosis

Fas is ubiquitously expressed on a variety of cells and triggers apoptosis, which have critical roles in the immune system. MicroRNAs (miRNAs) have been recently identified as regulators that modulate target gene expression and are involved in diverse biological processes, such as cell proliferation and apoptosis. This study was undertaken to investigate the contribution of miRNA in the regulation of Fas expression and Fas-mediated apoptosis. Bioinformatics analysis indicated that Fas was a potential target of let-7/miR-98 family. Indeed ectopic expression of let-7/miR-98 reduced, whereas knockdown of endogenous let-7/miR-98 increased the expression of Fas at both mRNA and protein levels. Let-7/miR-98 was verified to target Fas 3' untranslated region directly by site-directed gene mutagenesis and reporter gene assay. More importantly, introduction of let-7/miR-98 could decrease the sensitivity to Fas-induced apoptosis. Furthermore, let-7/miR-98 expression was reduced in activation-induced cell death process, accompanied by increased expression of Fas. In conclusion, our study first demonstrated that let-7/miR-98 regulated Fas expression and the sensitivity of Fas-mediated apoptosis.

Cardiovascular Disease, Single Nucleotide Polymorphisms; and the Renin Angiotensin System: Is There a MicroRNA Connection?

Essential hypertension is a complex disorder, caused by the interplay between many genetic variants, gene-gene interactions, and environmental factors. Given that the renin-angiotensin system (RAS) plays an important role in blood pressure (BP) control, cardiovascular regulation, and cardiovascular remodeling, special attention has been devoted to the investigation of single-nucleotide polymorphisms (SNP) harbored in RAS genes that may be associated with hypertension and cardiovascular disease. MicroRNAs (miRNAs) are a family of small, ∼21-nucleotide long, and nonprotein-coding RNAs that recognize target mRNAs through partial complementary elements in the 3′-untranslated region (3′-UTR) of mRNAs and inhibit gene expression by targeting mRNAs for translational repression or destabilization. Since miRNA SNPs (miRSNPs) can create, destroy, or modify miRNA binding sites, this review focuses on the hypothesis that transcribed target SNPs harbored in RAS mRNAs, that alter miRNA gene regulation and consequently protein expression, may contribute to cardiovascular disease susceptibility.

A functional variant in the 3'-UTR of angiopoietin-1 might reduce stroke risk by interfering with the binding efficiency of microRNA 211

Angiopoietin-1 is a vascular strengthening factor during vascular development and a protective factor for pathological vascular inflammation and leakage. Brain vascular leaking and inflammation are two important pathological processes of stroke; therefore, we hypothesized that variants of the microRNA-binding site in angiopoietin-1 would affect its expression and confer a risk of stroke. To test our hypothesis, a predicted microRNA-binding site was found in the 3'-UTR of angiopoietin-1 using bioinformatics; variant rs2507800 was identified to be located in the miR-211-binding site of angiopoietin-1. Secondly, the effects of the identified variant on angiopoietin-1 translation were assessed using a luciferase reporter assay and ELISA. We found that the A allele of rs2507800 suppressed angiopoietin-1 translation by facilitating miR-211 binding, but not the T allele. Subjects carrying the TT genotype had higher plasma angiopoietin-1 levels than those with the A allele. Finally, the association of the variant with stroke was tested in 438 stroke patients and 890 controls, and replicated in an independent population of 1791 stroke patients and 1843 controls. The TT genotype resulted in a significant reduction in overall stroke risk {OR, 0.51 [95% confidence interval (CI), 0.36-0.74], P = 0.0003}, ischemic stroke [OR, 0.56 (95% CI, 0.36-0.85), P = 0.007] and hemorrhagic stroke [OR, 0.46 (95% CI, 0.26-0.80), P = 0.007]. These results were confirmed in an independent study. Our results provide evidence that the TT genotype (rs2507800) in the 3'-UTR of angiopoietin-1 might reduce the risk of stroke by interfering with miR-211 binding.

microRNAs in diseases: from candidate to modifier genes

Until recently, the search for genetic factors predisposing or causing Mendelian diseases focused almost exclusively on protein coding sequences. As essential components of the regulatory system of gene expression, microRNAs (miRNAs) hold great promises into elucidating a number of inherited diseases. The herein review focuses on the genetic variations, whether copy number variation (CNV) or single nucleotide polymorphism (SNP), alternatively at the levels of the miRNA gene itself and of its target genes. We consider miRNA as the candidate gene, or the regulator of a disease-causing gene, or the modifier gene. The best paradigms of the field are presented in both monogenic diseases and complex traits. The computational tools, which are essential into identifying miRNAs and characterizing miRNA targets, are overviewed.