In silico prioritisation of microRNA-associated common variants in multiple sclerosis

Abstract

Abstract Background Genome-wide association studies (GWAS) have highlighted over 200 autosomal variants associated with multiple sclerosis (MS). However, variants in non-coding regions such as those encoding microRNAs have not been explored thoroughly, despite strong evidence of microRNA dysregulation in MS patients and model organisms. This study explores the effect of microRNA-associated variants in MS, through the largest publicly available GWAS, which involved 47,429 MS cases and 68,374 controls. Methods We identified SNPs within the coordinates of microRNAs, ± 5-kb microRNA flanking regions and predicted 3′UTR target-binding sites using miRBase v22, TargetScan 7.0 RNA22 v2.0 and dbSNP v151. We established the subset of microRNA-associated SNPs which were tested in the summary statistics of the largest MS GWAS by intersecting these datasets. Next, we prioritised those microRNA-associated SNPs which are among known MS susceptibility SNPs, are in strong linkage disequilibrium with the former or meet a microRNA-specific Bonferroni-corrected threshold. Finally, we predicted the effects of those prioritised SNPs on their microRNAs and 3′UTR target-binding sites using TargetScan v7.0, miRVaS and ADmiRE. Results We have identified 30 candidate microRNA-associated variants which meet at least one of our prioritisation criteria. Among these, we highlighted one microRNA variant rs1414273 (MIR548AC) and four 3′UTR microRNA-binding site variants within SLC2A4RG (rs6742), CD27 (rs1059501), MMEL1 (rs881640) and BCL2L13 (rs2587100). We determined changes to the predicted microRNA stability and binding site recognition of these microRNA and target sites. Conclusions We have systematically examined the functional, structural and regulatory effects of candidate MS variants among microRNAs and 3′UTR targets. This analysis allowed us to identify candidate microRNA-associated MS SNPs and highlights the value of prioritising non-coding RNA variation in GWAS. These candidate SNPs could influence microRNA regulation in MS patients. Our study is the first thorough investigation of both microRNA and 3′UTR target-binding site variation in multiple sclerosis using GWAS summary statistics.