Implication of genetic variants in primary microRNA processing sites in the risk of multiple sclerosis

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system with a well-established genetic contribution to susceptibility. Over 200 genetic regions have been linked to the inherited risk of developing MS, but the disease-causing variants and their functional effects at the molecular level are still largely unresolved. We hypothesised that MS-associated single-nucleotide polymorphisms (SNPs) affect the recognition and enzymatic cleavage of primary microRNAs (pri-miRNAs).

Our study focused on 11 pri-miRNAs (9 primate-specific) that are encoded in genetic risk loci for MS. The levels of mature miRNAs and potential isoforms (isomiRs) produced from those pri-miRNAs were measured in B cells obtained from the peripheral blood of 63 MS patients and 28 healthy controls. We tested for associations between SNP genotypes and miRNA expression in cis using quantitative trait locus (cis-miR-eQTL) analyses. Genetic effects on miRNA stem-loop processing efficiency were verified using luciferase reporter assays. Potential direct miRNA target genes were identified by transcriptome profiling and computational binding site assessment.

Mature miRNAs and isomiRs from hsa-mir-26a-2, hsa-mir-199a-1, hsa-mir-4304, hsa-mir-4423, hsa-mir-4464 and hsa-mir-4492 could be detected in all B-cell samples. When MS patient subgroups were compared with healthy controls, a significant differential expression was observed for miRNAs from the 5' and 3' strands of hsa-mir-26a-2 and hsa-mir-199a-1. The cis-miR-eQTL analyses and reporter assays pointed to a slightly more efficient Drosha-mediated processing of hsa-mir-199a-1 when the MS risk allele T of SNP rs1005039 is present. On the other hand, the MS risk allele A of SNP rs817478, which substitutes the first C in a CNNC sequence motif, was found to cause a markedly lower efficiency in the processing of hsa-mir-4423. Overexpression of hsa-mir-199a-1 inhibited the expression of 60 protein-coding genes, including IRAK2, MIF, TNFRSF12A and TRAF1. The only target gene identified for hsa-mir-4423 was TMEM47.

We found that MS-associated SNPs in sequence determinants of pri-miRNA processing can affect the expression of mature miRNAs. Our findings complement the existing literature on the dysregulation of miRNAs in MS. Further studies on the maturation and function of miRNAs in different cell types and tissues may help to gain a more detailed functional understanding of the genetic basis of MS.

This study was funded by the Rostock University Medical Center (FORUN program, grant: 889002), Sanofi Genzyme (grant: GZ-2016-11560) and Merck Serono GmbH (Darmstadt, Germany, an affiliate of Merck KGaA, CrossRef Funder ID: 10.13039/100009945, grant: 4501860307). NB was supported by the Stiftung der Deutschen Wirtschaft (sdw) and the FAZIT foundation. EP was supported by the Landesgraduiertenförderung Mecklenburg-Vorpommern.