Association Study of Two Polymorphisms in Vitamin D Pathway with Multiple Sclerosis in the Moroccan Population

Background: Hypovitaminosis D is reported through the literature to be involved in autoimmune diseases such as multiple sclerosis (MS). In the last decade, numerous studies have investigated the association of single nucleotide polymorphisms (SNPs) with MS, including rs2248359 (CYP24A1) and rs703842 (CYP27B1) that are involved in vitamin D metabolic pathway. However, results were conflicting, probably due to ethnic differences between the studied populations. In this context, the present study aimed to analyze the association between these two SNPs and MS within the Moroccan population. Methods: rs2248359 and rs703842 were genotyped in 113 patients and 146 healthy controls. To assess their association with the disease risk, we compared the genotypic and allelic frequencies between the study groups. We also explored their possible influence on certain clinical features (age at onset, type, disability status and severity score) and with vitamin D3 serum level (DSL) by comparing mean values of these variables between the different genotypes. Results: No statistically significant differences in the distribution of both SNPs were found between patients and controls. A trend has emerged concerning the minor G allele of rs703842 which appears to have a protective effect on developing MS, but this result remained slightly below significance. Also, the two polymorphisms had no impact on the clinical features tested and the DSL. Conclusion: There is no convincing evidence that rs2248359 and rs703842 are associated with MS risk, its clinical features or vitamin D level in MorocHow to cite this paper: Skalli, A., Ait Ben Haddou, E., Tazi-Ahnini, R., Mpandzou, G.A., Tibar, H., Bouslam, N., Benomar, A., Hajjout, K., Yahyaoui, M., El Fahime, E. and Bouhouche, A. (2017) Association Study of Two Polymorphisms in Vitamin D Pathway with Multiple Sclerosis in the Moroccan Population. World Journal of Neuroscience, 7, 363-375. https://doi.org/10.4236/wjns.2017.74028 Received: September 11, 2017 Accepted: October 21, 2017 Published: October 24, 2017 Copyright © 2017 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access


Introduction
Multiple sclerosis (MS) is a chronic autoimmune disease of the central nervous system and the most common non-traumatic cause of neurologic disability in young and middle-aged adults [1].It is characterized by a complex pathogenesis in which demyelinating lesions, axonal degeneration and gradual accumulation of sclerotic plaques are mainly responsible for this handicap [2].Both environmental risk factors and genetic susceptibility contribute to MS etiology [3] [4].
Among the known environmental risk factors, low concentration of vitamin D in blood circulation seems to have a significant role on MS pathogenesis, probably due to the intervention of vitamin D in several pathways of the immune system [5].Indeed, besides its fundamental role in calcium homeostasis and bone metabolism, vitamin D also has numerous additional functions, especially anti-inflammatory and immunomodulatory actions which could be involved in autoimmune diseases like MS.These include inhibition of B and T cells proliferation and differentiation, regulation of antigen presenting cells differentiation and conversion of cytokine environment from a pro-inflammatory state (Th1/Th17-cell-mediated) to a rather anti-inflammatory one (Th2-cell-mediated).
Furthermore, calcitriol (1.25(OH)₂D), which is the active metabolite of vitamin D, seems to have potential actions within central nervous system cells, such as neuronal functioning, neuroprotection and myelination.Resulting from this, the vitamin D status likely contributes to the risk of developing MS, and could also influence the main clinical variables of MS patients after the disease onset [6] [7].
Concerning the genetic susceptibility of MS, haplotypes in the major histocompatibility complex (MHC), particularly HLA DRB1 * 1501-DQB1 * 0602 (HLA-DR15), are the risk factors exerting the strongest effect [4].But over the last decade, numerous genome-wide association studies (GWASs) have reported more than 50 non-MHC loci involved in the MS predisposition.It was therefore suggested that the MS susceptibility profile consisted of many risk alleles, each contributing a relatively small effect to the overall risk [8]- [13].Among these, single nucleotide polymorphisms (SNPs) of vitamin D metabolizing genes, including rs2248359 and rs703842, have been found to be implicated in MS pathogenicity in Caucasian populations [12] [13].rs703842 is located in the 5' region of the CYP27B1 gene (12q13.1-q13.3)which encodes a 1-α-hydroxylase enzyme (cytochrome P450, family 27, subfamily B, peptide 1) that catalyzes the conversion of 25(OH)D to 1.25(OH)₂D3, the physiologically active form of vitamin D. Similarly, rs2248359 is located in the 5' region of the CYP24A1 gene (20q13) which encodes a 24-hydroxylase enzyme (cytochrome P450, family 24, subfamily A, peptide 1) that initiates deactivation and degradation of both 25(OH)D and 1.25(OH)₂D3.Because of their location in genomic regions predicted to be promoters for their neighboring genes, both SNPs might change the expression of these genes, and thus contribute to tissue availability of active vitamin D3.Thereby, they might modify the vitamin D influence on immune function, which could in turn affect susceptibility to MS [14] [15] [16].
Many candidate gene studies carried out in different countries have attempted to replicate the association between the two SNPs and MS but results are conflicting, possibly due to differences in ethnic origins of the concerned participants.
The present study aimed to analyze the genetic association of rs2248359 and rs703842 with MS disease in Moroccan population, and to explore their possible interactions with certain clinical features (age at onset, MS type, disability status and severity score) as well as with the vitamin D3 serum level (DSL).

Blood Collection and DNA Extraction
For both patients and controls, blood was drawn into a 6 mL EDTA-containing vacutainer tubes.After lysis of red cells, leucocyte pellets obtained were stored at −25˚C.DNA extraction was carried out using Isolate II genomic DNA Kit from Bioline, and the DNA quantity and quality was determined by the Nanodrop 8000 Spectrophotometer from Thermo Scientific.The DNA samples were finally stored at −20˚C until use.

SNP Genotyping
Two SNPs, rs2248359 and rs703842 respectively located near CYP24A1 and CYP27B1 genes, were selected for their role in the vitamin D metabolism and their potential association with MS based on previously published genome-wide and candidate-gene studies.
Genotyping of rs2248359 (C > T) was performed with the PCR-RFLP method.
The following primers were designed with Primer3 v4.0 software and used for PCR: forward 5'-AGTTAGGAAATGCGCCTTGAG-3' and reverse
Descriptive statistics were used to assess the demographic variables, the disease characteristics and the DSL.Prior to association tests, deviation from Hardy-Weinberg equilibrium (HWE) was examined in the control group using a chi-square (χ 2 ) test.Associations of rs2248359 and rs703842 with MS risk were assessed by comparing genotypic and allelic frequencies between patients and controls.Concerning genotypic frequencies, subjects were stratified according to three genetic models: additive, dominant (common allele homozygotes versus heterozygotes and minor allele homozygotes) and recessive (common allele homozygotes and heterozygotes versus minor allele homozygotes).Differences in the genotypic and allelic distribution between study groups were determined using the χ 2 test for statistical significance.Odds ratios (ORs) with 95% confidence intervals (CIs) were estimated for dominant, recessive and allelic models.
Four clinical characteristics (age at onset, MS type, EDSS score and MSSS) and DSL were also tested for association with both polymorphisms.Differences in mean values of these continuous variables between the different genotypes were evaluated using the univariate analysis of variance (ANOVA); except for MS type (categorical variable) which frequency analysis in regards to genotypes was performed using the contingency coefficient (C).In all cases, p values less than 0.05 (two-sided) were considered statistically significant.

Results
One hundred thirteen patients with MS and 146 unrelated healthy controls were enrolled.Detailed characteristics of the study population have been previously described (in press) and are summarized in Table 1.
Genotype frequencies for rs2248359 and rs703842 were in HWE among controls (p = 0.26 and p = 0.63 respectively).The genotypic and allelic distribution of both polymorphisms in MS cases and controls are presented in  statistically significant differences in the proportions of either polymorphism were found between the two study groups (p = 0.50 for rs2248359 and p = 0.18 for rs703842).However, regarding the dominant genetic and the allelic models , there was a trend for a higher frequency of the minor G allele of rs703842 in controls than patients, but these results remained slightly below significance (p = 0.07 and p = 0.08 respectively).These results suggest that both SNPs were not associated with MS and thus did not influence the risk of developing MS disease among Moroccans.

Discussion
Hypovitaminosis D has been proposed as a possible risk factor for MS over 40   [13].This cohort was then supplemented by about 15,000 patients and 25,000 controls [32].Both SNPs (rs2248359 and rs703842) have been reported to be associated with MS in the first study and after enlargement of the cohort.Against this background, and in an effort to investigate the association reported by the IMSGC/WTCCC2 study [13] within the Moroccan population, we conducted a preliminary assessment of rs2248359 and rs703842 effects on MS risk, some clinical features (age at onset, MS type, EDSS score and MSSS) and DSL.
Our results showed no association between the tested SNPs and MS status.However, since the p value was very close to the significance threshold when considering the dominant genetic and the allelic models, it should be mentioned that an interesting trend emerges concerning rs703842 (CYP27B1).In fact, the minor G allele carriers (one or two copies) are more likely to not develop MS than non-carriers.This result suggests that the G allele might have a protective effect on developing MS in our population.Results from previous studies diverge again.Regarding rs2248359 (CYP24A1), besides the IMSGC/WTCCC2 study [13], only one other found an association between this polymorphism and MS risk.It was conducted exclusively on African Americans subjects (803 MS cases and 1516 controls) in order to determine if MS-associated loci discovered in European populations also influence susceptibility in African Americans [33].However, two other studies found no association between rs2248359 and MS status within German [34] and Chinese populations [35].The same applies to rs703842 (CYP27B1), as in the IMSGC/WTCCC2 study [13], it has been identified as a confirmed MS susceptibility locus in another GWAS conducted by the Australia and New Zealand Multiple Sclerosis Genetics Consortium (ANZgene) on 3874 MS cases and 5723 controls of European ancestry [12].This result has been replicated by two studies which include in their samples subsets of the IMSGC or ANZgene cohorts [36] [37] and two other studies within Swedish [38] and Chinese populations [35].Negative results concerning such a link between rs703842 and MS status have also been reported within Canadian [39] and American [40] populations.Furthermore, we found no evidence that any of the tested SNPs modified the clinical phenotypes of age at onset, MS type, EDSS score or MSSS in our patients.The same finding was reported by an Australian study for rs703842 (CYP27B1) [41].We also did not detect an effect of the two SNPs on DSL, although they were expected to modulate it.This was consistent with an Anglo-American study showing no association between rs703842 (CYP27B1) and DSL [37], but contrasted with a Canadian one in which DSL differed significantly between the different CYP27B1 genotypes [39].
The divergence of results obtained from each study is likely due to ethnic differences between the studied populations, but environmental or behavioral factors related to vitamin D could also influence the expected association between the two SNPs and MS risk [37].Interestingly, similar studies on association of vitamin D receptor (VDR) polymorphisms with MS or with type 1 diabetes mel-World Journal of Neuroscience litus (another disease that has been associated with DSL through the literature) suggested that an association might only be penetrant in populations with an adequate vitamin D exposure [42] [43] [44].Hypothetically in a similar way, the T allele of CYP24A1 and the G allele of CYP27B1 that yield protection for developing MS might not affect disease status in a population generally deficient in vitamin D like ours.Indeed, about 98% of our control subjects had hypovitaminosis D (in press), which could explain why we did not find an association.Almost all studies related to CYP24A1 and CYP27B1 mentioned above had not measured DSLs, nor even taken into account determinants of circulating vitamin D (sun exposure, shadow-seeking behavior, skin pigmentation, use of sunscreen, wearing concealing clothing, dietary intake, multivitamin or cod-liver oil sup-plementation…).Therefore, it is possible that those which reported an association involved populations with high exposure to vitamin D.
There are certain limitations to the present study.Indeed, because of the relatively small sample size, we were underpowered to detect modest effects of SNPs on DSL or MS.Therefore, our findings only provide evidence against the existence of strong effects.Then, the case-control approach we used can lead to imperfect matching between cases and controls.Potential associations may be consequently under or over-estimated.A family-based design would have been more appropriate for this kind of study, but unfortunately, relatives were often not available for included patients.Otherwise, we did a single DSL measure for each participant, at the time of recruitment.Thus, DSLs obtained may not reflect the long-term DSL status.Moreover, due to the inclusion of prevalent MS cases, we could not always measure DSLs of patients at disease onset.This would have allowed us to better evaluate the SNPs effects on DSL in patients, before some disease symptoms change their behavioral habits regarding sun exposure.

Conclusion
In conclusion, there was no convincing evidence that rs2248359 (CYP24A1) and rs703842 (CYP27B1) were associated with MS risk, some clinical features (age at onset, MS type, EDSS score and MSSS) or DSL within the Moroccan population.These findings require further confirmation from large independent studies; and it would also be interesting to carry out similar studies in other populations ethnically close to ours (Algeria, Libya...) to see if there are similarities.
5'-GGATCAGGTTGAAAGGATTCG-3'.The 326 bp PCR-generated amplicons were systematically checked by electrophoresis on 2% agarose gel.Restriction enzyme SacII (recognition site: CCGCGG) from New England Biolabs served for the amplicons' digestion.SacII recognizes the C allele of the CYP24A1 gene and cut to generate two fragments of 100 bp and 226 bp which segregate after electrophoresis on 3% agarose gel.Three randomly selected samples of each genotype (namely 7% of the total sample) were subjected to regenotyping through DNA sequencing.Concordance rate for these quality-control samples was 100%.Rs703842 (A > G) has been genotyped using PCR followed by DNA sequencing.As previous, primers used for amplification were designed with Primer3 v4.0 software: forward 5'-GGGTAGTCATGAACACCAGT-3' and reverse 5'-GGGGAAGAAAGTTCTACGTAAT-3'.PCR products were purified with the ExoSAP-IT kit from Affymetrix, then sequenced using the reverse PCR primer as sequencing primer and the ABI Prism BigDye® v3.1 Cycle Sequence Kit on an ABI 3130xl Genetic Analyzer, both from Applied Biosystems, according to the manufacturer's protocol.Raw data were analyzed using Seqscape v2.5 software, also from Applied Biosystems.Synthesizing of primers and sequencing were both performed at in the Functional Genomics Platform (PGF) of the National Center for Scientific and Technical Research (CNRST, Rabat, Morocco).

MS type. The disease course was divided into two subcategories, PPMS and RRMS-SPMS since SPMS patients were first RRMS ones. The age at onset of the disease corresponds to self-reported patient's age at the apparition of first symptoms. To assess the disease severity, the EDSS score and the disease duration were used to calculate the Multiple Sclerosis Severity
A. Skalli et al.DOI: 10.4236/wjns.2017.74028366 World Journal of Neuroscience

Table 2 . No Table 1 .
Characteristics of MS patients and controls.

Table 2 .
Comparison of genotypic and allelic distribution of rs2248359 (CYP24A1) and rs703842 (CYP27B1) polymorphisms between MS patients and controls.

Table 4 .
Comparison of vitamin D3 serum levels between the different genotypes of rs2248359 (CYP24A1) and rs703842 (CYP27B1) polymorphisms in patients and controls.IMSGC/WTCCC2) in 15 different countries, on nearly 10,000 MS patients and 17,000 controls, almost all from European ancestry [19]n and standard deviation; DSL: vitamin D3 serum level; MS: Multiple Sclerosis; CYP24A1 genotypes (CC = wild type; CT = heterozygote; CC = mutant type); CYP27B1 genotypes (AA = wild type; AG = heterozygote; GG = mutant type).yearsago,initially to explain the observation of a latitude gradient in MS prevalence[18].Since then, several studies have focused on the subject to see if low DSLs were associated with a high risk of developing MS, but results obtained are conflicting.Indeed, some case-control studies reported no difference in DSLs between MS patients and controls[19][20] [21] [22] [23], while others have ob-MS severity markers; Further research is needed in this direction.The most important to date was a collaborative GWAS conducted by the International Multiple Sclerosis Genetics Consortium and the Wellcome Trust Case Control Con-A.Skalli et al.DOI: 10.4236/wjns.2017.74028370 World Journal of Neuroscience sortium 2 (