PolymiRTS Database: linking polymorphisms in microRNA target sites with complex traits

Insights into Online microRNA Bioinformatics Tools

MicroRNAs (miRNAs) are members of the small non-coding RNA family regulating gene expression at the post-transcriptional level. MiRNAs have been found to have critical roles in various biological and pathological processes. Research in this field has significantly progressed, with increased recognition of the importance of miRNA regulation. As a result of the vast data and information available regarding miRNAs, numerous online tools have emerged to address various biological questions related to their function and influence across essential cellular processes. This review includes a brief introduction to available resources for an investigation covering aspects such as miRNA sequences, target prediction/validation, miRNAs associated with disease, pathway analysis and genetic variants within miRNAs.

Application of the targeted sequencing approach reveals the single nucleotide polymorphism (SNP) repertoire in microRNA genes in the pig genome

MicroRNAs (miRNAs) are recognized as gene expression regulators, indirectly orchestrating a plethora of biological processes. Single nucleotide polymorphism (SNP), one of the most common genetic variations in the genome, is established to affect miRNA functioning and influence complex traits and diseases. SNPs in miRNAs have also been associated with important production traits in livestock. Thus, the aim of our study was to reveal the SNP variability of miRNA genes in the genome of the pig, which is a significant farm animal and large-mammal human model. To this end, we applied the targeted sequencing approach, enabling deep sequencing of specified genomic regions. As a result, 73 SNPs localized in 50 distinct pre-miRNAs were identified. In silico analysis revealed that many of the identified SNPs influenced the structure and energy of the hairpin precursors. Moreover, SNPs localized in the seed regions were shown to alter targeted genes and, as a result, enrich different biological pathways. The obtained results corroborate a significant impact of SNPs on the miRNA processing and broaden the state of knowledge in the field of animal genomics. We also report the targeted sequencing approach to be a promising alternative for the whole genome sequencing in miRNA genes focused studies.

Computational analysis of Cyclin D1 gene SNPs and association with breast cancer

CCND1 encodes for Cyclin D1 protein and single-nucleotide polymorphisms (SNPs) can modulate its activity. In the present study, the impact of CCND1 SNPs on structure and/or function of Cyclin D1 protein using in silico tools was investigated. Our analysis revealed only one splice site SNP (c.1988+5G<A) can effect CCND1 function. Subsequently, 78 out of 169 missense variants were predicted as pathogenic by Polyphen2, SIFT, PROVEAN, SNPs&GO, and PANTHER, and 4/78 missense SNPs were further evaluated because these four SNPs were found to be reside in highly conserved region of Cyclin D1. However, they did not show any major impact on tertiary structure and domain of Cyclin D1 but overall R15S and A190S has displayed a significant diseased phenotype and an altered molecular mechanism predicted by MutPred, FATHMM, SNPeffect, SNAP2, and PredictSNP. Consistently, A190S, R179L, and R15S may also cause a decrease in stability of Cyclin D1 anticipated by I-Mutant, HOPE and SNP effect. Furthermore, the Kaplan–Meier plotter has explained that high expression of CCND1 is associated with less survival rate of breast cancer patients. Altogether our study suggests that c.1988+5G<A, R15S, R179L, and A190S SNPs could directly or indirectly destabilize Cyclin D1.

FKBP-CaN-NFAT pathway polymorphisms selected by in silico biological function prediction are associated with tacrolimus efficacy in renal transplant patients

AIM

The aim of the present study was to investigate the potential effects of genetic variations in the FKBP-CaN-NFAT pathway on clinical events associated with tacrolimus efficacy in Chinese renal transplant patients.

METHODS

One hundred and forty Chinese renal transplant patients of Han ethnicity with over five years of follow-up were enrolled in our study. A pool of single nucleotide polymorphisms (SNPs) (1284 SNPs) was extracted from the Ensembl database according to chromosomal regions of the candidate genes. Next, 109 SNPs were screened out from this pool using multiple bioinformatics tools for subsequent genotyping using the MALDI-TOF-MS method. The associations of these candidate SNPs with acute rejection, nephrotoxicity, pneumonia and post-transplant estimated glomerular filtration rate (eGFR) were explored.

RESULTS

Fourty-four SNPs were found to be associated with tacrolimus-related clinical drug response. Specifically, eight SNPs were associated with the incidence of biopsy-proven acute rejection, four SNPs were associated with the rate of nephrotoxicity, 16 SNPs were correlated with the onset of pneumonia, and 26 SNPs were found to significantly influence post-transplant eGFR trend. An elaborate scoring system was implemented to prioritize the validation of these potentially causal SNPs. In particular, NFATC2 rs150348438 (G>T) performed well during integrative scoring (P_total=23.8) and was significantly associated with the occurrence of pneumonia (P = 0.0035, HR=0.91, 95% CI=0.85-0.97) and post-transplant eGFR levels (P = 0.000003).

CONCLUSIONS

NFATC2 rs150348438, rs6013219, rs1052653, and NFATC1 rs754093, ranking high in scoring, significantly affected the post-transplant eGFR and the incidence of pneumonia, acute rejection, and nephrotoxicity in renal transplant patients taking tacrolimus. Those SNPs may alter the expression and regulation of FKBP-CaN-NFAT pathway by influencing transcription regulation, mature mRNA degradation and RNA splicing, or protein coding. Critical SNPs of high ranking may serve as PD-associated pharmacogenetic biomarkers indicating individual response variability of TAC, and thus aid the clinical management of renal transplant patients.

A review of databases predicting the effects of SNPs in miRNA genes or miRNA-binding sites

Modern precision medicine comprises the knowledge and understanding of individual differences in the genomic sequence of patients to provide tailor-made treatments. Regularly, such variants are considered in coding regions only, and their effects are predicted based on their impact on the amino acid sequence of expressed proteins. However, assessing the effects of variants in noncoding elements, in particular microRNAs (miRNAs) and their binding sites, is important as well, as a single miRNA can influence the expression patterns of many genes at the same time. To analyze the effects of variants in miRNAs and their target sites, several databases storing variant impact predictions have been published. In this review, we will compare the core functionalities and features of these databases and discuss the importance of up-to-date data resources in the context of web applications. Finally, we will outline some recommendations for future developments in the field.

Genome-Wide Analysis of MicroRNA-related Single Nucleotide Polymorphisms (SNPs) in Mouse Genome

MicroRNAs are widely referred to as gene expression regulators for different diseases. The integration between single nucleotide polymorphisms (SNPs) and miRNAs has been associated with both human and animal diseases. In order to gain new insights on the effects of SNPs on miRNA and their related sequences, we steadily characterized a whole mouse genome miRNA related SNPs, analyzed their effects on the miRNA structural stability and target alteration. In this study, we collected 73643859 SNPs across the mouse genome, analyzed 1187 pre-miRNAs and 2027 mature miRNAs. Upon mapping the SNPs, 1700 of them were identified in 702 pre-miRNAs and 609 SNPs in mature miRNAs. We also discovered that SNP densities of the pre-miRNA and mature miRNAs are lower than the adjacent flanking regions. Also the flanking regions far away from miRNAs appeared to have higher SNP density. In addition, we also found that transitions were more frequent than transversions in miRNAs. Notably, 841 SNPs could change their corresponding miRNA's secondary structure from stable to unstable. We also performed target gain and loss analysis of 163 miRNAs and our results showed that few miRNAs remained unchanged and many miRNAs from wild mice gained target site. These results outline the first case of SNP variations in the mouse whole genome scale. Those miRNAs with changes in structure or target could be of interest for further studies.

Towards precision medicine: interrogating the human genome to identify drug pathways associated with potentially functional, population-differentiated polymorphisms

Drug response variations amongst different individuals/populations are influenced by several factors including allele frequency differences of single nucleotide polymorphisms (SNPs) that functionally affect drug-response genes. Here, we aim to identify drugs that potentially exhibit population differences in response using SNP data mining and analytics. Ninety-one pairwise-comparisons of >22,000,000 SNPs from the 1000 Genomes Project, across 14 different populations, were performed to identify ‘population-differentiated’ SNPs (pdSNPs). Potentially-functional pdSNPs (pf-pdSNPs) were then selected, mapped into genes, and integrated with drug–gene databases to identify ‘population-differentiated’ drugs enriched with genes carrying pf-pdSNPs. 1191 clinically-approved drugs were found to be significantly enriched (Z > 2.58) with genes carrying SNPs that were differentiated in one or more population-pair comparisons. Thirteen drugs were found to be enriched with such differentiated genes across all 91 population-pairs. Notably, 82% of drugs, which were previously reported in the literature to exhibit population differences in response were also found by this method to contain a significant enrichment of population specific differentiated SNPs. Furthermore, drugs with genetic testing labels, or those suspected to cause adverse reactions, contained a significantly larger number (P < 0.01) of population-pairs with enriched pf-pdSNPs compared with those without these labels. This pioneering effort at harnessing big-data pharmacogenomics to identify ‘population differentiated’ drugs could help to facilitate data-driven decision-making for a more personalized medicine.

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.

Architecture of polymorphisms in the human genome reveals functionally important and positively selected variants in immune response and drug transporter genes

Background

Genetic polymorphisms can contribute to phenotypic differences amongst individuals, including disease risk and drug response. Characterization of genetic polymorphisms that modulate gene expression and/or protein function may facilitate the identification of the causal variants. Here, we present the architecture of genetic polymorphisms in the human genome focusing on those predicted to be potentially functional/under natural selection and the pathways that they reside.

Results

In the human genome, polymorphisms that directly affect protein sequences and potentially affect function are the most constrained variants with the lowest single-nucleotide variant (SNV) density, least population differentiation and most significant enrichment of rare alleles. SNVs which potentially alter various regulatory sites, e.g. splicing regulatory elements, are also generally under negative selection.

Interestingly, genes that regulate the expression of transcription/splicing factors and histones are conserved as a higher proportion of these genes is non-polymorphic, contain ultra-conserved elements (UCEs) and/or has no non-synonymous SNVs (nsSNVs)/coding INDELs. On the other hand, major histocompatibility complex (MHC) genes are the most polymorphic with SNVs potentially affecting the binding of transcription/splicing factors and microRNAs (miRNA) exhibiting recent positive selection (RPS). The drug transporter genes carry the most number of potentially deleterious nsSNVs and exhibit signatures of RPS and/or population differentiation. These observations suggest that genes that interact with the environment are highly polymorphic and targeted by RPS.

Conclusions

In conclusion, selective constraints are observed in coding regions, master regulator genes, and potentially functional SNVs. In contrast, genes that modulate response to the environment are highly polymorphic and under positive selection.

Electronic supplementary material

The online version of this article (10.1186/s40246-018-0175-1) contains supplementary material, which is available to authorized users.

Exploring the selective vulnerability in Alzheimer disease using tissue specific variant analysis

The selective vulnerability of distinct regions of the brain is a critical factor in neurodegenerative disorders. In Alzheimer's disease (AD), neurons in hippocampus situated in medial temporal lobe are immensely damaged. Identifying tissue-specific variants is essential in order to perceive the selective vulnerability in AD. In current work, we aligned mRNA-seq data with HG19/HG38 genomic assembly and identified specific variations present in temporal, frontal and other lobes of the AD using sequence alignment map tools. We compared the results with the genome-wide association and gene expression quantitative trait loci studies of the various neurological disorders. We also distinguished variants and epitranscriptomic modifications through the RNA-modification database and evaluated the variant effect in the coding/UTR regions. In addition, we developed genetic and functional interaction networks to understand the relationship between predicted vulnerable variations and differentially expressed genes. We found that genes involved in gliogenesis, intermediate filament organization are altered in the temporal lobe. Oxidative phosphorylation, and calcium ion homeostasis are modified in the frontal lobe, and protein degradation, apoptotic signaling are altered in other lobes. From this study, we propose that disruption of glial cell structural integrity, defective gliogenesis, and failure in glia-neuron communication are the primary factors for selective vulnerability.

Lack of Evidence of the Role of APOA5 3'UTR Polymorphisms in Iranian Children and Adolescents with Metabolic Syndrome

BACKGROUND

Metabolic syndrome (MetS) is a complex and multifactorial disorder characterized by insulin resistance, dyslipidaemia, hyperglycemia, abdominal obesity, and elevated blood pressure. The apolipoprotein A5 (APOA5) gene variants have been reported to correlate with two major components of MetS, including low levels of high density lipoprotein cholesterol (HDL-C) and high levels of triglyceride. In the present study, we explored the associations between five single nucleotide polymorphisms (SNPs) of APOA5 gene and the MetS risk.

METHODS

In a case-control design, 120 Iranian children and adolescents with/without MetS were genotyped by polymerase chain reaction-sequencing for these SNPs. Then, we investigated the association of SNPs, individually or in haplotype constructs, with MetS risk.

RESULTS

The rs34089864 variant and H1 haplotype (harboring the two major alleles of rs619054 and rs34089864) were associated with HDL-C levels. However, there was no significant association between different haplotypes/individual SNPs and MetS risk.

CONCLUSION

These results presented no association of APOA5 3'UTR SNPs with MetS. Further studies, including other polymorphisms, are required to investigate the involvement of APOA5 gene in the genetic susceptibility to MetS in the pediatric age group.

Genetic variants in 3'-UTRs of MTHFR in the pregnancies complicated with preeclampsia and bioinformatics analysis

Preeclampsia (PE) as a pregnancy-specific disorder is the major cause of mortality and morbidity of mothers and fetuses. This study attempts to investigate the possible association between the 2572C>A (rs4846049) and 4869C>G (rs1537514) polymorphisms in the 3'- untranslated region of MTHFR gene and the risk of PE. A total of 198 patients diagnosed with PE and 171 unrelated, age matched healthy pregnant women, were recruited for this case-control study. The MTHFR 2572C>A and 4869C>G genotyping was performed by the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. The CG genotype of MTHFR 4869C>G was associated with decreased risk of PE, and this genotype was found to be a protective factor for PE susceptibility. There was no significant difference in the genotypes of MTHFR 2572C>A polymorphism between PE patients and control group. The frequency of combined AC/CG genotypes of MTHFR 2572C>A and 4869C>G polymorphisms were less frequent in PE patients and were associated with a lower risk of PE. The C-G and A-G haplotypes of MTHFR 2572C>A and 4869C>G polymorphisms were significantly lower in PE patients. In conclusion, the CG genotype of MTHFR 4869C>G polymorphism was associated with a lower risk of PE. No association was found between MTHFR 2572C>A polymorphism and PE.

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.

miRNAs target databases: developmental methods and target identification techniques with functional annotations

PURPOSE

microRNA (miRNA) regulates diverse biological mechanisms and metabolisms in plants and animals. Thus, the discoveries of miRNA has revolutionized the life sciences and medical research.The miRNA represses and cleaves the targeted mRNA by binding perfect or near perfect or imperfect complementary base pairs by RNA-induced silencing complex (RISC) formation during biogenesis process. One miRNA interacts with one or more mRNA genes and vice versa, hence takes part in causing various diseases. In this paper, the different microRNA target databases and their functional annotations developed by various researchers have been reviewed. The concurrent research review aims at comprehending the significance of miRNA and presenting the existing status of annotated miRNA target resources built by researchers henceforth discovering the knowledge for diagnosis and prognosis.

METHODS AND RESULTS

This review discusses the applications and developmental methodologies for constructing target database as well as the utility of user interface design. An integrated architecture is drawn and a graphically comparative study of present status of miRNA targets in diverse diseases and various biological processes is performed. These databases comprise of information such as miRNA target-associated disease, transcription factor binding sites (TFBSs) in miRNA genomic locations, polymorphism in miRNA target, A-to-I edited target, Gene Ontology (GO), genome annotations, KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, target expression analysis, TF-miRNA and miRNA-mRNA interaction networks, drugs-targets interactions, etc.

CONCLUSION

miRNA target databases contain diverse experimentally and computationally predicted target through various algorithms. The comparison of various miRNA target database has been performed on various parameters. The computationally predicted target databases suffer from false positive information as there is no common theory for prediction of miRNA targets. The review conclusion emphasizes the need of more intelligent computational improvement for the miRNA target identification, their functional annotations and datasbase development.

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.

Genetic variants in microRNAs and their binding sites within gene 3'UTRs associate with susceptibility to age-related macular degeneration

Age-related macular degeneration (AMD), the leading cause of blindness in the elderly, is a complex disease that results from multiple genetic and environmental factors. MicroRNAs (miRNAs) are small noncoding RNAs that post-transcriptionally regulate target mRNAs and are frequently implicated in human diseases. Here, we investigated the association of genetic variants in miRNAs and miRNA-binding sites within gene 3'-untranslated regions (3'UTRs) with AMD using data from the largest AMD genome-wide association study. First, we identified three variants in miRNAs significantly associated with AMD. These include rs2168518:G>A in the miR-4513 seed sequence, rs41292412:C>T in pre-miR-122/miR-3591, and rs4351242:C>T in the terminal-loop of pre-miR-3135b. We demonstrated that these variants reduce expression levels of the mature miRNAs in vitro and pointed the target genes that may mediate downstream effects of these miRNAs in AMD. Second, we identified 54 variants (in 31 genes) in miRNA-binding sites associated with AMD. Based on stringent prioritization criteria, we highlighted the variants that are more likely to have an impact on the miRNA-target interactions. Further, we selected rs4151672:C>T within the CFB 3'UTR and experimentally showed that while miR-210-5p downregulates expression of CFB, the variant decreases miR-210-5p-mediated repression of CFB. Together, our findings support the notion that miRNAs may play a role in AMD.

Identification of coagulation gene 3'UTR variants that are potentially regulated by microRNAs

MicroRNAs have been recognized as critical regulators of gene expression and might affect the risk of venous thrombosis. We aimed to identify 3' untranslated region (UTR) variants in coagulation genes that influence coagulation factor levels and venous thrombosis risk. The 3'UTR of coagulation genes were sequenced in subjects with extremely high or low plasma levels of these factors in two case-control studies. In total, 28 variants were identified. Five single nucleotide polymorphisms (SNPs) were predominantly present in one extreme level group (F2 rs1799963, F8 rs1050705 and F11 rs4253429, rs4253430 and rs1062547). Additional to F2 rs1799963, F8 rs1050705 (in men) and F11 rs4253430 were associated with an increased risk of venous thrombosis albeit confidence intervals were wide. The three F11 SNPs were in high linkage disequilibrium with functional variants rs2289252 and rs2036914. Rs1062547 and rs4253430 were associated with a significant increase of plasma FXI activity in heterozygotes and homozygotes in wild-type controls. In silico prediction revealed that these SNPs might disturb the binding sites of miR-544 and miR-513a-3p. Only miR-544 provoked a significant decrease of the luciferase activity that was not observed with a rs4253430 mutated vector. In conclusion, these results reinforce that microRNAs are candidates to play a role in haemostasis and complex disorders, such as thrombosis.

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.

Functional variomics and network perturbation: connecting genotype to phenotype in cancer

Proteins interact with other macromolecules in complex cellular networks for signal transduction and biological function. In cancer, genetic aberrations have been traditionally thought to disrupt the entire gene function. It has been increasingly appreciated that each mutation of a gene could have a subtle but unique effect on protein function or network rewiring, contributing to diverse phenotypic consequences across cancer patient populations. In this Review, we discuss the current understanding of cancer genetic variants, including the broad spectrum of mutation classes and the wide range of mechanistic effects on gene function in the context of signalling networks. We highlight recent advances in computational and experimental strategies to study the diverse functional and phenotypic consequences of mutations at the base-pair resolution. Such information is crucial to understanding the complex pleiotropic effect of cancer genes and provides a possible link between genotype and phenotype in cancer.

Identification of a New Single-nucleotide Polymorphism within the Apolipoprotein A5 Gene, Which is Associated with Metabolic Syndrome

Background:

Metabolic syndrome (MetS) is a common disorder which is a constellation of clinical features including abdominal obesity, increased level of serum triglycerides (TGs) and decrease of serum high-density lipoprotein-cholesterol (HDL-C), elevated blood pressure, and glucose intolerance. The apolipoprotein A5 (APOA5) is involved in lipid metabolism, influencing the level of plasma TG and HDL-C. In the present study, we aimed to investigate the associations between four INDEL variants of APOA5 gene and the MetS risk.

Materials and Methods:

In this case–control study, we genotyped 116 Iranian children and adolescents with/without MetS by using Sanger sequencing method for these INDELs. Then, we explored the association of INDELs with MetS risk and their clinical components by logistic regression and one-way analysis of variance analyses.

Results:

We identified a novel insertion polymorphism, c. *282–283 insAG/c. *282–283 insG variant, which appears among case and control groups. rs72525532 showed a significant difference for TG levels between various genotype groups. In addition, there were significant associations between newly identified single-nucleotide polymorphism (SNP) and rs72525532 with MetS risk.

Conclusions:

These results show that rs72525532 and the newly identified SNP may influence the susceptibility of the individuals to MetS.

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

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

Genome-wide identification of microRNA-related variants associated with risk of Alzheimer's disease

MicroRNAs (miRNAs) serve as key post-transcriptional regulators of gene expression. Genetic variation in miRNAs and miRNA-binding sites may affect miRNA function and contribute to disease risk. Here, we investigated the extent to which variants within miRNA-related sequences could constitute a part of the functional variants involved in developing Alzheimer’s disease (AD), using the largest available genome-wide association study of AD. First, among 237 variants in miRNAs, we found rs2291418 in the miR-1229 precursor to be significantly associated with AD (p-value = 6.8 × 10⁻⁵, OR = 1.2). Our in-silico analysis and in-vitro miRNA expression experiments demonstrated that the variant’s mutant allele enhances the production of miR-1229-3p. Next, we found miR-1229-3p target genes that are associated with AD and might mediate the miRNA function. We demonstrated that miR-1229-3p directly controls the expression of its top AD-associated target gene (SORL1) using luciferase reporter assays. Additionally, we showed that miR-1229-3p and SORL1 are both expressed in the human brain. Second, among 42,855 variants in miRNA-binding sites, we identified 10 variants (in the 3′ UTR of 9 genes) that are significantly associated with AD, including rs6857 that increases the miR-320e-mediated regulation of PVRL2. Collectively, this study shows that miRNA-related variants are associated with AD and suggests miRNA-dependent regulation of several AD genes.

Tumor-promoting function of single nucleotide polymorphism rs1836724 (C3388T) alters multiple potential legitimate microRNA binding sites at the 3'-untranslated region of ErbB4 in breast cancer

ErbB4 can act as either a tumor-suppressor gene or an oncogene in breast cancer. Multiple genetic factors including single nucleotide polymorphisms (SNPs) affect gene expression patterns. Multiple 3'-untranslated region (3'-UTR) SNPs reside within the target binding site of microRNAs, which can strengthen or weaken binding to target genes. The present study aimed to predict potential 3'‑UTR variants of ErbB4 that alter the target binding site of microRNAs (miRNAs) and to clarify the association of the potential variant with the risk of developing breast cancer. In silico prediction was performed to identify potential functional SNPs within miRNA target binding sites in the 3'‑UTR of ErbB4. Thus, 146 patients and controls were genotyped using restriction fragment length polymorphism-polymerase chain reaction. In addition to the Cochran-Armitage test for trend, allele and genotype frequency differences were determined to investigate the association between rs1836724 and the susceptibility to breast cancer. Bioinformatics analysis identified rs1836724 to be a polymorphism in the seed region of four miRNA binding sites (hsa-miR335-5p, hsa-miR-28-5p, has‑miR‑708‑5p and has‑miR‑665), which may participate in the development of breast cancer. Logistic regression data indicated that the T allele of the polymorphism [OR (95% CI)=1.72 (1.056‑2.808), P=0.029] is associated with the risk of breast cancer. Using bioinformatics tools, a correlation was indicated between the presence of the T allele and a reduction in ErbB4 RNA silencing based on miRNA interaction. Furthermore, case subgroup data analysis revealed an association between the C/T genotype and an ER positive phenotype [OR (95% CI)=6.00 (1.082‑33.274), P=0.028] compared with the T/T genotype. ErbB4 and estrogen receptor 1 (ESR1) are regulated by identical miRNAs thus there may be a competition for binding sites. Due to this pattern, if the interaction between miRNAs with one gene is reduced, it may be consistent with the increase in interaction with another one. Therefore, more interaction with rs1836724 C variant within ErbB4 may be associated with higher expression of ESR1 (ER‑positive phenotype). miRNAs interact with ErbB4 mRNA more frequently when it carries C allele at the rs1836724 position compared with the T carriers. Therefore, the identical miRNA interacts with ESR1 less frequently when ErbB4 mRNA has a C allele. Therefore, ESR1 expression may be higher when ErbB4 mRNA has a C allele.

MiRNA-Related Genetic Variations Associated with Radiotherapy-Induced Toxicities in Patients with Locally Advanced Non-Small Cell Lung Cancer

Severe radiation-induced toxicities limit treatment efficacy and compromise outcomes of lung cancer. We aimed to identify microRNA-related genetic variations as biomarkers for the prediction of radiotherapy-induced acute toxicities. We genotyped 233 SNPs (161 in microRNA binding site and 72 in processing gene) and analyzed their associations with pneumonitis and esophagitis in 167 stage III NSCLC patients received definitive radiation therapy. Sixteen and 11 SNPs were associated with esophagitis and pneumonitis, respectively. After multiple comparison correction, RPS6KB2:rs10274, SMO:rs1061280, SMO:rs1061285 remained significantly associated with esophagitis, while processing gene DGCR8:rs720014, DGCR8:rs3757, DGCR8:rs1633445 remained significantly associated with pneumonitis. Patients with the AA genotype of RPS6KB2:rs10274 had an 81% reduced risk of developing esophagitis (OR: 0.19, 95% CI: 0.07–0.51, p = 0.001, q = 0.06). Patients with the AG+GG genotype of SMO:rs1061280 had an 81% reduced risk of developing esophagitis (OR: 0.19, 95% CI: 0.07–0.53, p = 0.001, q = 0.06). Patients with the GG+GA genotype of DGCR8:rs720014 had a 3.54-fold increased risk of pneumonitis (OR: 3.54, 95% CI: 1.65–7.61, p <0.05, q <0.1). Significantly cumulative effects of the top SNPs were observed for both toxicities (P-trend <0.001). Using bioinformatics tools, we found that the genotype of rs10274 was associated with altered expression of the RPS6KB2 gene. Gene-based analysis showed DGCR8 (p = 0.010) and GEMIN4 (p = 0.039) were the top genes associated with the risk of developing pneumonitis. Our results provide strong evidence that microRNA-related genetic variations contribute to the development of radiotherapy-induced acute esophagitis and pneumonitis and could thus serve as biomarkers to help accurately predict radiotherapy-induced toxicity in NSCLC patients.

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.

Altered Gene Expression Associated with microRNA Binding Site Polymorphisms

Allele-specific gene expression associated with genetic variation in regulatory regions can play an important role in the development of complex traits. We hypothesized that polymorphisms in microRNA (miRNA) response elements (MRE-SNPs) that either disrupt a miRNA binding site or create a new miRNA binding site can affect the allele-specific expression of target genes. By integrating public expression quantitative trait locus (eQTL) data, miRNA binding site predictions, small RNA sequencing, and Argonaute crosslinking immunoprecipitation (AGO-CLIP) datasets, we identified genetic variants that can affect gene expression by modulating miRNA binding efficiency. We also identified MRE-SNPs located in regions associated with complex traits, indicating possible causative mechanisms associated with these loci. The results of this study expand the current understanding of gene expression regulation and help to interpret the mechanisms underlying eQTL effects.

Micro-RNA Binding Site Polymorphisms in the WFS1 Gene Are Risk Factors of Diabetes Mellitus

The absolute or relative lack of insulin is the key factor in the pathogenesis of diabetes mellitus. Although the connection between loss of function mutations of the WFS1 gene and DIDMOAD-syndrome including diabetes mellitus underpins the significance of wolframin in the pathogenesis, exact role of WFS1 polymorphic variants in the development of type 1 and type 2 diabetes has not been discovered yet. In this analysis, 787 patients with diabetes and 900 healthy people participated. Genotyping of the 7 WFS1 SNPs was carried out by TaqMan assays. Association study was performed by χ2-test in combination with correction for multiple testing. For functional analysis, the entire 3' UTR of the WFS1 gene was subcloned in a pMIR-Report plasmid and relative luciferase activities were determined. Linkage disequilibrium analysis showed a generally high LD within the investigated region, however the rs1046322 locus was not in LD with the other SNPs. The two miR-SNPs, rs1046322 and rs9457 showed significant association with T1DM and T2DM, respectively. Haplotype analysis also confirmed the association between the 3' UTR loci and both disease types. In vitro experiments showed that miR-185 reduces the amount of the resulting protein, and rs9457 miRSNP significantly influences the rate of reduction in a luciferase reporter assay. Genetic variants of the WFS1 gene might contribute to the genetic risk of T1DM and T2DM. Furthermore demonstrating the effect of rs9457 in binding of miR-185, we suggest that the optimal level of wolframin protein, potentially influenced by miR-regulation, is crucial in normal beta cell function.

A genetic variant in the seed region of miR-4513 shows pleiotropic effects on lipid and glucose homeostasis, blood pressure, and coronary artery disease

MicroRNAs (miRNA) play a crucial role in the regulation of diverse biological processes by post-transcriptional modulation of gene expression. Genetic polymorphisms in miRNA-related genes can potentially contribute to a wide range of phenotypes. The effect of such variants on cardiometabolic diseases has not yet been defined. We systematically investigated the association of genetic variants in the seed region of miRNAs with cardiometabolic phenotypes, using the thus far largest genome-wide association studies on 17 cardiometabolic traits/diseases. We found that rs2168518:G>A, a seed region variant of miR-4513, associates with fasting glucose, low-density lipoprotein-cholesterol, total cholesterol, systolic and diastolic blood pressure, and risk of coronary artery disease. We experimentally showed that miR-4513 expression is significantly reduced in the presence of the rs2168518 mutant allele. We sought to identify miR-4513 target genes that may mediate these associations and revealed five genes (PCSK1, BNC2, MTMR3, ANK3, and GOSR2) through which these effects might be taking place. Using luciferase reporter assays, we validated GOSR2 as a target of miR-4513 and further demonstrated that the miRNA-mediated regulation of this gene is changed by rs2168518. Our findings indicate a pleiotropic effect of miR-4513 on cardiometabolic phenotypes and may improve our understanding of the pathophysiology of cardiometabolic diseases.

Identification of a novel FGFRL1 MicroRNA target site polymorphism for bone mineral density in meta-analyses of genome-wide association studies

MicroRNAs (miRNAs) are critical post-transcriptional regulators. Based on a previous genome-wide association (GWA) scan, we conducted a polymorphism in microRNA target sites (poly-miRTS)-centric multistage meta-analysis for lumbar spine (LS)-, total hip (HIP)- and femoral neck (FN)-bone mineral density (BMD). In stage I, 41 102 poly-miRTSs were meta-analyzed in seven cohorts with a genome-wide significance (GWS) α = 0.05/41 102 = 1.22 × 10(-6). By applying α = 5 × 10(-5) (suggestive significance), 11 poly-miRTSs were selected, with FGFRL1 rs4647940 and PRR5 rs3213550 as top signals for FN-BMD (P = 7.67 × 10(-6) and 1.58 × 10(-5)) in gender-combined sample. In stage II in silico replication (two cohorts), FGFRL1 rs4647940 was the only signal marginally replicated for FN-BMD (P = 5.08 × 10(-3)) at α = 0.10/11 = 9.09 × 10(-3). PRR5 rs3213550 was also selected based on biological significance. In stage III de novo genotyping replication (two cohorts), FGFRL1 rs4647940 was the only signal significantly replicated for FN-BMD (P = 7.55 × 10(-6)) at α = 0.05/2 = 0.025 in gender-combined sample. Aggregating three stages, FGFRL1 rs4647940 was the single stage I-discovered and stages II- and III-replicated signal attaining GWS for FN-BMD (P = 8.87 × 10(-12)). Dual-luciferase reporter assays demonstrated that FGFRL1 3' untranslated region harboring rs4647940 appears to be hsa-miR-140-5p's target site. In a zebrafish microinjection experiment, dre-miR-140-5p is shown to exert a dramatic impact on craniofacial skeleton formation. Taken together, we provided functional evidence for a novel FGFRL1 poly-miRTS rs4647940 in a previously known 4p16.3 locus, and experimental and clinical genetics studies have shown both FGFRL1 and hsa-miR-140-5p are important for bone formation.

Genetic Variations in MicroRNA-Binding Sites Affect MicroRNA-Mediated Regulation of Several Genes Associated With Cardio-metabolic Phenotypes

BACKGROUND

Genome-wide association studies enabled us to discover a large number of variants and genomic loci contributing to cardiovascular and metabolic disorders. However, because the vast majority of the identified variants are thought to merely be proxies for other functional variants, the causal mechanisms remain to be elucidated. We hypothesized that the part of the functional variants involved in deregulating cardiometabolic genes is located in microRNA (miRNA)-binding sites.

METHODS AND RESULTS

Using the largest genome-wide association studies available on glycemic indices, lipid traits, anthropometric measures, blood pressure, coronary artery diseases, and type 2 diabetes mellitus, we identified 11,067 variants that are associated with cardiometabolic phenotypes. Of these, 230 variants are located within miRNA-binding sites in the 3'-untranslated region of 155 cardiometabolic genes. Thirty-seven of 230 variants were found to fulfill our predefined criteria for being functional in their genomic loci. Ten variants were subsequently selected for experimental validation based on genome-wide association studies results, expression quantitative trait loci (eQTL) analyses, and coexpression of their host genes and regulatory miRNAs in relevant tissues. Luciferase reporter assays revealed an allele-specific regulation of genes hosting the variants by miRNAs. These cotransfection experiments showed that rs174545 (FADS1:miR-181a-2), rs1059611 (LPL:miR-136), rs13702 (LPL:miR-410), rs1046875 (FN3KRP:miR-34a), rs7956 (MKRN2:miR-154), rs3217992 (CDKN2B:miR-138-2-3p), and rs11735092 (HSD17B13:miR-375) decrease or abrogate miRNA-dependent regulation of the genes. Conversely, 2 variants, rs6857 (PVRL2:miR-320e) and rs907091 (IKZF3:miR-326), were shown to enhance the activity of miRNAs on their host genes.

CONCLUSIONS

We provide evidence for a model in which polymorphisms in miRNA-binding sites can both positively and negatively affect miRNA-mediated regulation of cardiometabolic genes.

Genetic variants in miRNA processing genes and pre-miRNAs are associated with the risk of chronic lymphocytic leukemia

Genome wide association studies (GWAS) have identified several low-penetrance susceptibility alleles in chronic lymphocytic leukemia (CLL). Nevertheless, these studies scarcely study regions that are implicated in non-coding molecules such as microRNAs (miRNAs). Abnormalities in miRNAs, as altered expression patterns and mutations, have been described in CLL, suggesting their implication in the development of the disease. Genetic variations in miRNAs can affect levels of miRNA expression if present in pre-miRNAs and in miRNA biogenesis genes or alter miRNA function if present in both target mRNA and miRNA sequences. Therefore, the present study aimed to evaluate whether polymorphisms in pre-miRNAs, and/or miRNA processing genes contribute to predisposition for CLL. A total of 91 SNPs in 107 CLL patients and 350 cancer-free controls were successfully analyzed using TaqMan Open Array technology. We found nine statistically significant associations with CLL risk after FDR correction, seven in miRNA processing genes (rs3805500 and rs6877842 in DROSHA, rs1057035 in DICER1, rs17676986 in SND1, rs9611280 in TNRC6B, rs784567 in TRBP and rs11866002 in CNOT1) and two in pre-miRNAs (rs11614913 in miR196a2 and rs2114358 in miR1206). These findings suggest that polymorphisms in genes involved in miRNAs biogenesis pathway as well as in pre-miRNAs contribute to the risk of CLL. Large-scale studies are needed to validate the current findings.

miR2GO: comparative functional analysis for microRNAs

miR2GO is a web-based platform for comparative analyses of human miRNA functions. It includes two programs: miRmut2GO and miRpair2GO. miRmut2GO implements a knowledge-based method to assess the functional effects of genetic and somatic mutations in microRNA seed regions. The functional effects of a mutation are analysed by semantic comparison of enriched gene ontology (GO) annotations of the target gene sets for the wild-type and mutated alleles. miRpair2GO compares the functions of two different miRNAs based on the enriched functional annotations of their target gene sets.

AVAILABILITY AND IMPLEMENTATION

The miR2GO web server is available at http://compbio.uthsc.edu/miR2GO.

The association of common polymorphisms in miR-196a2 with waist to hip ratio and miR-1908 with serum lipid and glucose

OBJECTIVE

MicroRNAs (miRNAs) have been implicated in the regulation of cardiometabolic disorders. Given the crucial role of miRNAs in gene expression, genetic variation within miRNA genes is expected to affect miRNA function and substantially contribute to disease risk.

METHODS

2,320 variants in miRNA-encoding sequences were systematically retrieved, and their associations with 17 cardiometabolic traits/diseases were investigated, using genome-wide association studies (GWAS) on glycemic indices, anthropometric measures, lipid traits, blood pressure, coronary artery disease, and type 2 diabetes. Next, target genes of the identified miRNAs that may mediate their effect on the phenotypes were examined. Furthermore, trans- expression quantitative trait loci analysis and luciferase reporter assay to provide functional evidence for our findings were performed.

RESULTS

rs11614913:C/T in miR-196a2 was associated with waist to hip ratio (P-value=1.7 × 10(-5) , β = 0.023). Two target genes, SFMBT1 and HOXC8, which may mediate this association were identfied, and they were shown experimentally as direct targets of miR-196a2. Moreover, rs174561:C/T in miR-1908 was found to be associated with total cholesterol (P-value=6.5 × 10(-16) , β=0.044), LDL-cholesterol (P-value=4.3 × 10(-18) , β=0.049), HDL-cholesterol (P-value=1.7 × 10(-6) , β=0.026), triglyceride (P-value=7.8 × 10(-14) , β=0.038), and fasting glucose (P-value=4.3 × 10(-10) , β=0.02). In addition, a number of miR-1908 target genes were highlighted as potential mediators.

CONCLUSIONS

The results indicated miRNA-dependent regulation of fat distribution by miR-196a2 and of lipid metabolism by miR-1908.

Identification of a functional SNP in the 3'UTR of CXCR2 that is associated with reduced risk of lung cancer

Global changes in gene expression accompany the development of cancer. Thus, inherited variants in miRNA-binding sites are likely candidates for conferring inherited susceptibility. Using an in silico approach, we compiled a comprehensive list of SNPs predicted to modulate miRNA binding in genes from several key lung cancer pathways. We then investigated whether these SNPs were associated with lung cancer risk in two independent populations. In general, SNPs in miRNA-binding sites are rare. However, some allelic variation was observed. We found that rs1126579 in CXCR2 was associated with a reduced risk of lung cancer in both European American [ORTT vs. CC 0.56 (0.37-0.88); P = 0.008] and Japanese [ORTT vs. CC 0.62 (0.38-1.00); P = 0.049] populations. Furthermore, we found that the SNP disrupted a novel binding site for miR-516a-3p, led to a moderate increase in CXCR2 mRNA and protein expression, and increased MAPK signaling. Moreover, analysis of rs1126579 with serum levels of IL8, its endogenous ligand, supported an interaction whereby rs1126579-T and high serum IL8 conferred synergistic protection from lung cancer. Our findings demonstrate a function for a 3'UTR SNP in modulating CXCR2 expression, signaling, and susceptibility to lung cancer.

MicroRNA-binding site polymorphisms in hematological malignancies

Dysregulation of microRNA networks has been implicated in hematological malignancies. One of the reasons for disturbed miRNA-mediated regulation are polymorphisms in miRNA-binding sites (miRSNPs), which alter the strength of miRNA interaction with target transcripts. In the recent years the first findings of miRSNPs associated with risk and prognosis in hematological malignancies have been reported. From the studies described in this review miRSNPs not only emerge as novel markers of risk and prognosis but can also lead to better understanding of the role of miRNAs in regulating gene expression in health and disease.

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.

Computational Methods for MicroRNA Target Prediction

MicroRNAs (miRNAs) have been identified as one of the most important molecules that regulate gene expression in various organisms. miRNAs are short, 21–23 nucleotide-long, single stranded RNA molecules that bind to 3' untranslated regions (3' UTRs) of their target mRNAs. In general, they silence the expression of their target genes via degradation of the mRNA or by translational repression. The expression of miRNAs, on the other hand, also varies in different tissues based on their functions. It is significantly important to predict the targets of miRNAs by computational approaches to understand their effects on the regulation of gene expression. Various computational methods have been generated for miRNA target prediction but the resulting lists of candidate target genes from different algorithms often do not overlap. It is crucial to adjust the bioinformatics tools for more accurate predictions as it is equally important to validate the predicted target genes experimentally.

A functional SNP catalog of overlapping miRNA-binding sites in genes implicated in prion disease and other neurodegenerative disorders

The involvement of SNPs in miRNA target sites remains poorly investigated in neurodegenerative disease. In addition to associations with disease risk, such genetic variations can also provide novel insight into mechanistic pathways that may be responsible for disease etiology and/or pathobiology. To identify SNPs associated specifically with degenerating neurons, we restricted our analysis to genes that are dysregulated in CA1 hippocampal neurons of mice during early, preclinical phase of Prion disease. The 125 genes chosen are also implicated in other numerous degenerative and neurological diseases and disorders and are therefore likely to be of fundamental importance. We predicted those SNPs that could increase, decrease, or have neutral effects on miRNA binding. This group of genes was more likely to possess DNA variants than were genes chosen at random. Furthermore, many of the SNPs are common within the human population, and could contribute to the growing awareness that miRNAs and associated SNPs could account for detrimental neurological states. Interestingly, SNPs that overlapped miRNA-binding sites in the 3'-UTR of GABA-receptor subunit coding genes were particularly enriched. Moreover, we demonstrated that SNP rs9291296 would strengthen miR-26a-5p binding to a highly conserved site in the 3'-UTR of gamma-aminobutyric acid receptor subunit alpha-4.

microRNAs: role in leukemia and their computational perspective

MicroRNAs (miRNAs) belong to the family of noncoding RNAs (ncRNAs) and had gained importance due to its role in complex biochemical pathways. Changes in the expression of protein coding genes are the major cause of leukemia. Role of miRNAs as tumor suppressors has provided a new insight in the field of leukemia research. Particularly, the miRNAs mediated gene regulation involves the modulation of multiple mRNAs and cooperative action of different miRNAs to regulate a particular gene expression. This highly complex array of regulatory pathway network indicates the great possibility in analyzing and identifying novel findings. Owing to the conventional, slow experimental identification process of miRNAs and their targets, the last decade has witnessed the development of a large amount of computational approaches to deal with the complex interrelations present within biological systems. This article describes the various roles played by miRNAs in regulating leukemia and the role of computational approaches in exploring new possibilities.

Genetic variation in the 3'-UTR of CYP1A2, CYP2B6, CYP2D6, CYP3A4, NR1I2, and UGT2B7: potential effects on regulation by microRNA and pharmacogenomics relevance

Introduction: Pharmacogenomics research has concentrated on variation in genes coding for drug metabolizing enzymes, transporters and nuclear receptors. However, variation affecting microRNA could also play a role in drug response. This project set out to investigate potential microRNA target sites in 11 genes and the extent of variation in the 3′-UTR of six selected genes; CYP1A2, CYP2B6, CYP2D6, CYP3A4, NR1I2, and UGT2B7.

Methods: Fifteen microRNA target prediction algorithms were used to identify microRNAs predicted to regulate 11 genes. The 3′-UTR of the 6 genes which topped the list of potential microRNA targets was sequenced in 30 black South Africans. In addition, genetic variants within these genes were investigated for interference with mRNA-microRNA interactions. Potential effects of observed variants were determined using in silico prediction tools.

Results: The 11 genes coding for DMEs, transporters and nuclear receptors were predicted to be targets of microRNAs with CYP2B6, NR1I2 (PXR), CYP3A4, and CYP1A2, interacting with the most microRNAs. The majority of identified genetic variants were predicted to interfere with microRNA regulation. For example, the variant, rs1054190C in NR1I2 was predicted to result in the presence of a binding site for the microRNA miR-1250-5p, while the variant rs1054191G was predicted to result in the absence of a recognition site for miR-371b-3p, miR-4258 and miR-4707-3p. Fifteen of the seventeen, novel variants occurred within microRNA target sequences.

Conclusion: The 3′-UTR harbors variation that is likely to influence regulation of specific genes by microRNA. In silico prediction followed by functional validation could aid in decoding the contribution of variation in the 3′-UTR, to some unexplained heritability that affects drug response. Understanding the specific role of each microRNA may lead to identification of markers for targeted therapy and therefore improve personalized drug treatment.

Polymorphisms in microRNA target sites modulate risk of lymphoblastic and myeloid leukemias and affect microRNA binding

BACKGROUND

MicroRNA dysregulation is a common event in leukemia. Polymorphisms in microRNA-binding sites (miRSNPs) in target genes may alter the strength of microRNA interaction with target transcripts thereby affecting protein levels. In this study we aimed at identifying miRSNPs associated with leukemia risk and assessing impact of these miRSNPs on miRNA binding to target transcripts.

METHODS

We analyzed with specialized algorithms the 3' untranslated regions of 137 leukemia-associated genes and identified 111 putative miRSNPs, of which 10 were chosen for further investigation. We genotyped patients with acute myeloid leukemia (AML, n = 87), chronic myeloid leukemia (CML, n = 140), childhood acute lymphoblastic leukemia (ALL, n = 101) and healthy controls (n = 471). Association between SNPs and leukemia risk was calculated by estimating odds ratios in the multivariate logistic regression analysis. For miRSNPs that were associated with leukemia risk we performed luciferase reporter assays to examine whether they influence miRNA binding.

RESULTS

Here we show that variant alleles of TLX1_rs2742038 and ETV6_rs1573613 were associated with increased risk of childhood ALL (OR (95% CI) = 3.97 (1.43-11.02) and 1.9 (1.16-3.11), respectively), while PML_rs9479 was associated with decreased ALL risk (OR = 0.55 (0.36-0.86). In adult myeloid leukemias we found significant associations between the variant allele of PML_rs9479 and decreased AML risk (OR = 0.61 (0.38-0.97), and between variant alleles of IRF8_ rs10514611 and ARHGAP26_rs187729 and increased CML risk (OR = 2.4 (1.12-5.15) and 1.63 (1.07-2.47), respectively). Moreover, we observed a significant trend for an increasing ALL and CML risk with the growing number of risk genotypes with OR = 13.91 (4.38-44.11) for carriers of ≥3 risk genotypes in ALL and OR = 4.9 (1.27-18.85) for carriers of 2 risk genotypes in CML. Luciferase reporter assays revealed that the C allele of ARHGAP26_rs187729 creates an illegitimate binding site for miR-18a-3p, while the A allele of PML_rs9479 enhances binding of miR-510-5p and the C allele of ETV6_rs1573613 weakens binding of miR-34c-5p and miR-449b-5p.

CONCLUSIONS

Our study implicates that microRNA-binding site polymorphisms modulate leukemia risk by interfering with the miRNA-mediated regulation. Our findings underscore the significance of variability in 3' untranslated regions in leukemia.

In silico mapping of polymorphic miRNA-mRNA interactions in autoimmune thyroid diseases

MicroRNAs (miRNAs) are important regulators of gene expression and translation. The genetic variants altering miRNA targets have been associated with many diseases. Here we systematically mapped the human genetic polymorphisms that may affect miRNA-mRNA interactions in the autoimmune thyroid disease (AITD) pathway. We also mapped the polymorphic miRNA target sites in the genes that have been linked to AITDs or other thyroid-related diseases/phenotypes in genome-wide association studies (GWAS). These genetic polymorphisms may potentially contribute to the pathogenesis of AITDs and other thyroid diseases. The polymorphic miRNA-mRNA interactions we mapped in the AITD pathway and the GWAS-informed thyroid disease loci may provide insights into the possible miRNA-mediated molecular mechanisms through which genetic variants assert their influences on thyroid diseases and phenotypes.