Abstract
Multiple sclerosis (MS) is a chronic debilitating neurological condition with a wide range of phenotype variability. A complex interplay of genetic and environmental factors contributes to disease onset and progression in MS patients. Vitamin D deficiency is a known susceptibility factor for MS, however the underlying mechanism of vitamin D-gene interactions in MS etiology is still poorly understood. Vitamin D receptor super-enhancers (VSEs) are enriched in MS risk variants and may modulate these environment-gene interactions. mRNA expression in total of 64 patients with contrasting MS severity was quantified in select genes. First, RNA-seq was performed on a discovery cohort (10 mild, 10 severe MS phenotype) and ten genes regulated by VSEs that have been linked to MS risk were analyzed. Four candidates showed a significant positive association (GRINA, PLEC, PARP10, and LRG1) in the discovery cohort and were then quantified using digital droplet PCR (ddPCR) in a validation cohort (33 mild, 11 severe MS phenotype). A significant differential expression persisted in the validation cohort for three of the VSE-MS genes: GRINA (p = 0.0138), LRG1 (p = 0.0157), and PLEC (p = 0.0391). In summary, genes regulated by VSE regions that contain known MS risk variants were shown to have differential expression based on disease severity (p<0.05). The findings implicate a role for vitamin D super-enhancers in modulating disease activity. In addition, expression levels may have some utility as prognostic biomarkers in the future.
Introduction
Multiple sclerosis (MS) is a chronic and debilitating autoimmune condition of the central nervous system (1). A complex interplay of genetic and environmental factors contributes to disease onset and progression in MS patients (2). Vitamin D deficiency is a known susceptibility factor for MS, however the underlying mechanism of vitamin D gene-environment interactions in MS etiology is still poorly understood (3).
There are several lines of evidence supporting the role of vitamin D in MS susceptibility, including higher MS prevalence in regions with lower annual UVB exposure (4–7), and migrant studies showing changes to MS risk based on latitude changes (8–10). Increased actinic sun damage, and reporting more time outdoors during childhood/adolescence, have both been associated with reduced risk of MS (11–13). Large prospective studies measuring vitamin D status prior to MS onset showed that participants with low serum vitamin D levels were at greater risk of developing MS (14–16). Clinical studies have reported that vitamin D supplementation in MS patients lowered disease relapse risk (17–21), while lower vitamin D status at MS onset is associated with worse cognitive function and neuronal integrity (22). MS genetic risk loci, identified by GWAS, are enriched with binding sites for the vitamin D receptor (VDR) (23), and several MS risk variants have been localized to VDR binding sites (24).
Recent in silico evidence suggests that vitamin D receptor (VDR) super-enhancers (VSE) may be influenced by nearby MS risk SNPs, thereby impacting the activity of these VSEs in relation to disease (25). Vitamin D exerts its biological effects when the active hormone [1,25-dihydroxyvitamin D, 1,25(OH)2D] binds VDR, which then acts as a transcription factor upon binding to VDR response elements in target genes (26). These VDR target genes are often associated with clusters of VDR binding sites that form super-enhancers, regulatory genomic regions consisting of dense clusters of enhancers (27). VSEs contain clusters of VDR binding sites, and loop to hundreds of promoter regions, upregulating vitamin D target genes (27).
Lu et al. (25) identified specific genes near VSE regions that overlap with MS risk variants. We investigated the expression of ten VSE-MS-linked genes in a cohort of 64 persons with MS representing the two extremes of disease activity, to assess whether expression differences were associated with mild vs. severe MS phenotype.
Methods
Ethics statement
The Conjoint Health Research Ethics Board (REB17-1193) at the University of Calgary provided approval for this study. Written informed consent was obtained from all participants.
Study population and design
Selection criteria and experimental methodology have been described previously (28). Briefly, study participants (n = 64) representing a range of phenotypic severity were selected from a large prospective cohort consisting of 2,831 participants recruited for an ongoing prospective cohort study at the University of Calgary MS Clinic in Calgary, Alberta. Patients were identified according to their age-related MS severity (ARMSS) scores over three or more clinical visits. 600 patients with the lowest ARMSS scores, and 600 patients with the highest ARMSS scores were designated as mild and severe phenotypes, respectively. From this subset, 64 patients were recruited for this pilot project: 43 patients with a mild phenotype (mean ARMSS score 1.34, range 0.06–3.08) and 21 patients with a severe phenotype (mean ARMSS score 6.90, range: 5.73–9.94). Clinical parameters were assessed to ensure relative homogeneity between the two groups, including age of MS onset, age at recruitment, sex, and use of disease modifying therapies. Mann-Whitney-Wilcoxon test for ordinal variables and Chi-square test for nominal variables were used to assess statistical significance.
A two-phase discovery and validation approach was employed. Candidate MS-VSE genes were identified in a discovery cohort consisting of 20 patients (10 mild, 10 severe phenotype) in which transcriptome sequencing (RNA-seq) was conducted. Genes demonstrating significant differential expression between mild and severe phenotypic subgroups were then retained and quantified in a validation cohort using digital droplet PCR (ddPCR).
RNA isolation and quality control
The Paxgene RNA Purification Kit was utilized to isolate total RNA from whole blood samples of recruited participants. Sample integrity was verified by fluorimetry using an RNA-specific dye (Qubit fluorimeter), and RIN analysis with an Agilent TapeStation 2200 (RIN score range 6.4–8.4).
High throughput RNA-sequencing
Low Sample protocol was followed using the #20020596 TruSeq Stranded Total RNA Library Preparation Kit (H/M/R) on 500 ng of each RNA sample. RiboZero magnetic beads were used to remove rRNA and keep the Illumina TruSeq i7 indices in the RNA samples. The fragments present in the RNA samples were enriched via 15 cycles of PCR amplification. The resulting libraries were validated by TapeStation analysis and qPCR using the Kapa qPCR Library Quant Kit for Illumina. Following validation, the libraries were then pooled and sequenced on an Illumina NextSeq 500 sequencer. The pool of 20 samples underwent single-end sequencing, which involved three consecutive high-output NextSeq V2 sequencing runs, consisting of 75 cycles for each run. Paired-end sequencing was also conducted on the 20-sample pool using a single 150 cycle (2 × 75 bp) high-output NextSeq V2 sequencing run. The single-end sequencing and paired-sequencing yielded an average of 80 million clusters PF per sample and 26 million clusters PF per samples, respectively.
Bioinformatics analysis
Transcriptomic data was processed by aligning counts to the EnsEMBL GRCh38.p12 genome reference assembly utilizing FeatureCounts. Differential expression between mild and severe phenotype patients was then conducted using the DeSeq2 R Package. Ten candidate genes were selected based on previous findings. Firstly, five genes regulated by VSEs which contain known MS risk variants were previously identified (25) and thus selected for inclusion: UBASH3B, IRF8, PLEC, PARP10, GRINA. The second set of genes overlap with MS risk variants and were under significant regulation by 1,25(OH)D via VSEs, with high expression levels (25): DENND6B, USP2, ASAP2, SEMA6B, LRG1. The targets that showed a statistically significant difference (Log2 fold change, p < 0.05, Wald test) were then retained for analysis in the validation cohort.
Digital droplet PCR (ddPCR) validation analysis
Validation analysis were performed using digital droplet PCR due to its increased sensitivity and absolute quantification method compared to conventional PCR techniques. ddPCR custom primers were designed using the Ensembl database and Primer3 software, which generated sequences for forward and reverse primers (Supplementary Table 1). The Bio-Rad iScript cDNA Synthesis Kit and random hexamers were used to perform reverse transcription on the remaining 44 RNA samples. Based on the Bio-Rad QX200 ddPCR system protocol, custom primers designed for each candidate and Bio-Rad ddPCR EvaGreen Supermix was used. The expression of each candidate was normalized to the expression of HPRT1 and B2M given their relatively consistent expression profiles. The normalized counts for each candidate were then analyzed using a one-tailed t-test in GraphPad Prism 8.4.1 to determine if a significant difference in expression existed between the mild and severe MS phenotypes.
Results
Clinical characteristics of the study participants are listed in Table 1. There were large differences in ARMSS scores between the mild and severe phenotype groups, as expected due to study design (Table 1). Across all study participants, the mean ARMSS scores were 1.34 and 6.90, for the mild and severe phenotypes, respectively. Other variables, including the age at onset, age at recruitment, and sex (% female) were not significantly different between mild and severe phenotypes (Table 1). The severe phenotype groups in both cohorts consisted of only female patients (Table 1). The number of participants in each group receiving at least one of the MS disease modifying therapies, as well as the types of treatments, did not differ significantly between the groups; 58.1% in the mild group, 61.9% in the severe group (data not shown).
Table 1
| Discovery cohort | |||
|---|---|---|---|
| Mild phenotype (n = 10) | Severe phenotype (n = 10) | p-value | |
| Mean age (range) | 53.4 (32.4–67.4) | 52.1 (39.6–59.5) | 0.7133 |
| Mean age at MS Onset (range) | 35.6 (20.9–46.2) | 28.7 (14.8–46.1) | 0.7649 |
| Sex (% female) | 70% | 100% | 0.2105 |
| Mean ARMSS score (range) | 1.177 (0.060–2.194) | 7.255 (6.187–9.401) | < 0.0001 |
| Validation cohort | |||
| Mild phenotype (n = 33) | Severe phenotype (n = 11) | p-value | |
| Mean age (range) | 57.8 (45.3–70.1) | 54.6 (34.4–67.1) | 0.4480 |
| Mean age at MS Onset (range) | 33.9 (16.6–56.2) | 32.2 (14.6-47.3) | 0.09692 |
| Sex (% female) | 60.6% | 100% | 0.0189 |
| Mean ARMSS score (range) | 1.864 (0.225–3.017) | 6.518 (5.804–7.891) | <0.0001 |
Clinical characteristics of included study participants.
In the discovery cohort, whole transcriptome sequencing was conducted for the ten candidate MS-linked-VSE genes selected based on prior evidence (Figure 1). Mean expression levels were analyzed for significant differences in the mild vs. severe phenotype groups (≥ 2-fold-change, p < 0.05). Four of the candidate MS-linked-VSE genes GRINA, PLEC, PARP10, and LRG1 were significantly upregulated in the severe phenotype MS subgroup (Table 2), and thus were evaluated in the validation cohort (n = 44) using custom-designed ddPCR assays. This analysis demonstrated three genes that were differentially expressed between MS phenotypes including GRINA, PLEC, and LRG1 (Figure 2).
Figure 1
Table 2
| Discovery cohort | Validation cohort | |||||||
|---|---|---|---|---|---|---|---|---|
| MS-VSE candidates | Mean counts [Standard deviation] | Mean counts [standard deviation] | ||||||
| Mild | Severe | log2FC | p-value, differential transcript levels | Mild | Severe | log2FC | p-value, normalized to HPRT1 and B2M | |
| GRINA | 832 [268] | 1,240 [368] | 0.57 | 0.0034* | 0.85 [0.60] | 1.45 [1.11] | 0.77 | 0.0138* |
| PLEC | 2,885 [1,019] | 3,608 [1,606] | 0.32 | 0.0006* | 0.90 [0.44] | 1.29 [1.01] | 0.52 | 0.0381* |
| PARP10 | 717 [255] | 1,071 [652] | 0.58 | 0.0028* | 1.00 [0.70] | 0.99 [0.86] | −0.01 | 0.4925 |
| LRG1 | 529 [140] | 780 [348] | 0.78 | 0.0063* | 0.84 [0.54] | 1.47 [1.36] | 0.81 | 0.0157* |
Differential gene expression between mild and severe MS phenotypes in the discovery and validation cohorts; mean counts were determined by RNA-seq (number of sequenced fragments) and by ddPCR (number of gene copies, normalized to HPRT1 and B2M expression), respectively.
The asterisks (*) indicate statistically significant results (p < 0.05).
Figure 2
Discussion
This pilot study compared expression of MS-VSE genes in mild vs. severe MS phenotypes, in order to elucidate a potential role for vitamin D in MS etiology, as well as identifying potential future biomarkers for differentiating these phenotype subgroups. Vitamin D super-enhancers in MS is a relatively new area of research and this is the first study, to our knowledge, to directly investigate the putative role of VSE-linked genes associated with MS risk SNPs in MS phenotypic severity.
Vitamin D deficiency has long been implicated in MS susceptibility, yet the exact vitamin D gene-environment or epigenetic interactions in MS etiology remain unclear. A better understanding may at best help in prevention or modulation of disease severity, or at least create awareness of the importance of vitamin D sufficiency. Active vitamin D, 1,25(OH)D exerts its biological effects by binding to the VDR and altering gene expression. Vitamin D super-enhancers contain multiple VDR binding sites, which then loop to transcription start site of target gene(s) to stimulate transcription (29). VSEs are believed to have more pronounced effects on gene expression than traditional enhancers (and many are signal-inducible via 1,25(OH)D exposure (26, 27). Recent in silico evidence identified specific genes with VSE regions overlapping with MS risk variants as well as VSE genes overlapping with MS risk variants (25). We compared the expression of ten of these identified candidates in our cohort. It is hypothesized that MS risk SNPs within VSEs, or genes associated with VSEs may alter the VSE activity to impact disease risk and activity (30).
RNA-sequencing and ddPCR were used to assess differential expression of select candidates. Three VSE genes (PLEC, GRINA, LRG1) that were previously identified as linked to MS variants (25) showed consistent differences in expression between the mild vs. severe phenotype groups in both the discovery and confirmation cohorts. It should be noted that for neurodegenerative diseases, even small fold-changes in expression levels can relate to disease development and activity (24, 31, 32). Interestingly, a significant increase in transcription levels was associated with these VSE-genes in the severe phenotype. The findings implicate a role for vitamin D super-enhancers in modulating disease via environment-gene interactions. In addition, expression levels may have some utility as prognostic biomarkers in the future.
The molecular mechanisms involved in VSE regulation of gene expression and its impact on MS has not been studied. GRINA encodes a glutamate ionotropic receptor N-methyl-D-aspartate (NMDA)-type subunit associated protein 1, expressed on the postsynaptic membrane of neuronal synapses (33). The receptor is an excitatory glutamate-gated ion channel found throughout the body (34). Excessive glutamate at the sites of demyelination and overstimulation of glutamate receptors leading to neuronal death have been reported (35, 36). The PLEC gene encodes the cytoskeleton plectin protein, which is found in many tissues, including nervous tissue (37, 38). The degradation of cerebral tissue integrity and elasticity in MS patients, especially in glial cells (39), may be affected. LRG1 encodes the leucine-rich alpha-2-glycoprotein, which has recently been reported to promote angiogenesis in the brain (40, 41).
While it is hypothesized that VSE-related MS risk variants may lead to changes in the epigenomic landscape involved in MS pathogenesis, we cannot rule out that the converse is true—it may be that disease severity affects VDR related changes in chromatin accessibility and thus regulation of gene expression. It may also be a combination of both effects. Studying expression in individuals with first de-myelinating events, and/or longitudinally, would be needed to better support any potential cause-and-effect relationship.
Given the small sample size in this pilot study, follow up in larger phenotype cohorts or case-control designs, will be useful to validate a role for VSE genes in MS risk and severity, as well as to explore a potential role as clinically useful biomarkers. In particular, a comparison of mild cases vs. healthy controls and severe cases vs. healthy controls would be warranted in this validation. Furthermore, the severe phenotype group in the validation cohort contained exclusively females, and obtaining larger sample sizes with more males are warranted to confirm the findings. While the cohorts do not differ in the baseline demographic characteristics, we acknowledge the limitation on the results of not having more precise 1:1 propensity score matching between the individuals in the discovery and validation cohorts.
Another potential limitation of this study involves the use of whole blood, which has both advantages and disadvantages, as we have discussed previously (28). Many recent publications studying MS biomarkers have used whole blood (42). Other blood sample types (serum, plasma) require more onerous standardization to avoid variability in sample contents introduced (43). We believe that whole blood remains a reasonable choice for gene expression studies in MS because the collection is easily standardized, and there is precedent from numerous prior studies. Recent publications continue to demonstrate the important role of peripheral blood cells in the pathogenesis of MS (44, 45). Nonetheless, future studies extracting transcripts from different cell types could help determine if the results are more generalizable, given that the genomic binding pattern of VDR can be variable between tissues (29). While our findings support a link between VSE genes and MS risk variants, the mechanisms involved have not been studied. Extension to murine EAE models could explore functional effects.
The study of vitamin D super-enhancers, and their effect on human health, is a relatively new and growing area of study. There are very limited publications in this area for the field of multiple sclerosis. The candidate genes investigated may be important players in the interaction between the environment (vitamin D) and genetic risk associated with multiple sclerosis.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics statement
The Conjoint Health Research Ethics Board (REB17-1193) at the University of Calgary provided approval for this study. Written informed consent was obtained from all participants.
Author contributions
SO: conceived study and design, lab work, data analysis, data interpretation, manuscript writing, and manuscript editing. AS: lab work, data collection, data analysis, and manuscript writing. MG: data collection, data analysis, and manuscript editing. KM: technical assistance and experiments. LM: cohort collection and study design. AK: cohort collection, study design, and data analysis. GP: cohort collection, study design, data interpretation, and manuscript revision. All authors reviewed and approved the final manuscript version.
Funding
The research was supported by: Mount Royal University Faculty of Science and Technology New Faculty Funding (No. 101683); Mount Royal University Faculty of Science and Technology (No. 102207) and Alberta MS Network Summer Studentship funding (AS); Markin Studentship (MG); Infrastructure used was supported by a Canada Foundation for Innovation grant to Gerald Pfeffer; Operating funds for sequencing was supported from an Alberta MS Collaboration grant to Gerald Pfeffer. Sequencing was performed at the Center for Health Genomics and Informatics (CHGI), which is supported by the University of Calgary.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2022.1064008/full#supplementary-material
References
1.
FreedmanMSDevonshireVDuquettePGiacominiPSGiulianiFLevinMCet al. Treatment optimization in multiple sclerosis: Canadian MS working group recommendations. Can J Neurol Sci. (2020) 47:437–55. 10.1017/cjn.2020.66
2.
WaubantELucasRMowryEGravesJOlssonTAlfredssonLet al. Environmental and genetic risk factors for MS: an integrated review. Ann Clin Transl Neurol. (2019) 6:1905–22. 10.1002/acn3.50862
3.
Berlanga-TaylorAJDisantoGEbersGCRamagopalanSV. Vitamin D–gene interactions in multiple sclerosis. J Neurol Sci. (2011) 311:32–6. 10.1016/j.jns.2011.08.041
4.
AchesonEDBachrachCA. The distribution of multiple sclerosis in U. S veterans by birthplace. Am J Epidemiol. (1960) 72:88–99. 10.1093/oxfordjournals.aje.a120137
5.
OrtonSMWaldLConfavreuxCVukusicSKrohnJPRamagopalanSVet al. Association of UV radiation with multiple sclerosis prevalence and sex ratio in France. Neurology. (2011) 76:425–31. 10.1212/WNL.0b013e31820a0a9f
6.
SimpsonSBlizzardLOtahalPVan der MeiITaylorB. Latitude is significantly associated with the prevalence of multiple sclerosis: a meta-analysis. J Neurol Neurosurg Psychiatry. (2011) 82:1132–41. 10.1136/jnnp.2011.240432
7.
JelinekGADe LiveraAMMarckCHBrownCRNeateSLTaylorKLet al. Lifestyle, medication and socio-demographic determinants of mental and physical health-related quality of life in people with multiple sclerosis. BMC Neurol. (2016) 16:1–9. 10.1186/s12883-016-0763-4
8.
DeanGKurtzkeJF. On the risk of multiple sclerosis according to age at immigration to South Africa. BMJ. (1971) 3:725–9. 10.1136/bmj.3.5777.725
9.
HammondSR. The age-range of risk of developing multiple sclerosis: evidence from a migrant population in Australia. Brain. (2000) 123:968–74. 10.1093/brain/123.5.968
10.
NielsenNMCornGFrischMStenagerEKoch-HenriksenNWohlfahrtJet al. Multiple sclerosis among first- and second-generation immigrants in Denmark: a population-based cohort study. Brain. (2019) 142:1587–97. 10.1093/brain/awz088
11.
van der MeiIAFPonsonbyALDwyerTBlizzardLSimmonsRTaylorBVet al. Past exposure to sun, skin phenotype, and risk of multiple sclerosis: case-control study. BMJ. (2003) 327:1–6. 10.1136/bmj.327.7410.316
12.
TremlettHZhuFAscherioAMungerKL. Sun exposure over the life course and associations with multiple sclerosis. Neurology. (2018) 90:e1191–9. 10.1212/WNL.0000000000005257
13.
MagalhaesSPugliattiMRiiseTMyhrKMCiampiABjornevikKet al. Shedding light on the link between early life sun exposure and risk of multiple sclerosis: results from the EnvIMS Study. Int J Epidemiol. (2019) 48:1073–82. 10.1093/ije/dyy269
14.
MungerKLLevinLIHollisBWHowardNSAscherioA. Serum 25-hydroxyvitamin D levels and risk of multiple sclerosis. JAMA. (2006) 296:2832–8. 10.1001/jama.296.23.2832
15.
NielsenNMMungerKLKoch-HenriksenNHougaardDMMagyariMJørgensenKTet al. Neonatal vitamin D status and risk of multiple sclerosis: a population-based case-control study. Neurology. (2017) 88:44–51. 10.1212/WNL.0000000000003454
16.
MungerKLHongellKÅivoJSoilu-HänninenMSurcelHMAscherioA. 25-hydroxyvitamin D deficiency and risk of MS among women in the finnish maternity cohort. Neurology. (2017) 89:1578–83. 10.1212/WNL.0000000000004489
17.
RuniaTFHopWCJRijke YBdeBuljevacDHintzenRQ. Lower serum vitamin D levels are associated with a higher relapse risk in multiple sclerosis. Neurology. (2012) 79:261–6. 10.1212/WNL.0b013e31825fdec7
18.
MowryEMKruppLBMilazzoMChabasDStroberJBBelmanALet al. Vitamin D status is associated with relapse rate in pediatric-onset multiple sclerosis. Ann Neurol. (2010) 67:618–24. 10.1002/ana.21972
19.
SimpsonSTaylorBBlizzardLPonsonbyALPittasFTremlettHet al. Higher 25-hydroxyvitamin D is associated with lower relapse risk in multiple sclerosis. Ann Neurol. (2010) 68:193–203. 10.1002/ana.22043
20.
BerezowskaMCoeSDawesH. Effectiveness of vitamin D supplementation in the management of multiple sclerosis: a systematic review. Int J Mol Sci. (2019) 20:1–19. 10.3390/ijms20061301
21.
MicleaAMicleaMPistorMHoepnerAChanAHoepnerR. Vitamin D supplementation differentially affects seasonal multiple sclerosis disease activity. Brain Behav. (2017) 7:1–6. 10.1002/brb3.761
22.
CorteseMMungerKLMartínez-LapiscinaEHBarroCEdanGFreedmanMSet al. Vitamin D, smoking, EBV, and long-term cognitive performance in MS: 11-year follow-up of BENEFIT. Neurology. (2020) 94:e1950–60. 10.1212/WNL.0000000000009371
23.
DisantoGSandveGKBerlanga-TaylorAJRagneddaGMorahanJMWatsonCTet al. Vitamin D receptor binding, chromatin states and association with multiple sclerosis. Hum Mol Genet. (2012) 21:3575–86. 10.1093/hmg/dds189
24.
RamagopalanSVHegerABerlangaAJMaugeriNJLincolnMRBurrellAet al. A ChIP-seq defined genome-wide map of vitamin D receptor binding: associations with disease and evolution. Genome Res. (2010) 20:1352–60. 10.1101/gr.107920.110
25.
LuMMcComishBJBurdonKPTaylorBVKörnerH. The association between vitamin D and multiple sclerosis risk: 1,25(OH)2D3 induces super-enhancers bound by VDR. Front Immunol. (2019) 10:1–16. 10.3389/fimmu.2019.00488
26.
LuCLYeihDFHouYCJow GM LiZYLiuWCet al. The emerging role of nutritional vitamin D in secondary hyperparathyroidism in CKD. Nutrients. (2018) 10:1–19. 10.3390/nu10121890
27.
HanelAMalmbergHRCarlbergC. Genome-wide effects of chromatin on vitamin D signaling. J Mol Endocrinol. (2020) 64:R45–56. 10.1530/JME-19-0246
28.
GuptaMMartensKMetzLMde KoningAJPfefferG. Long noncoding RNAs associated with phenotypic severity in multiple sclerosis. Mult Scler Relat Disord. (2019) 36:1–6. 10.1016/j.msard.2019.101407
29.
CarlbergCVelleuerE. Vitamin D and the risk for cancer: a molecular analysis. Biochem Pharmacol. (2022) 196:1–10. 10.1016/j.bcp.2021.114735
30.
YamagataKNakayamadaSTanakaY. Critical roles of super-enhancers in the pathogenesis of autoimmune diseases. Inflamm Regen. (2020) 40:2–9. 10.1186/s41232-020-00124-9
31.
Cooper-KnockJKirbyJFerraiuoloLHeathPRRattrayMShawPJ. Gene expression profiling in human neurodegenerative disease. Nat Rev Neurol. (2012) 8:518–30. 10.1038/nrneurol.2012.156
32.
SonntagKC. miRNAs and deregulated gene expression networks in neurodegeneration. Brain Res. (2010) 1338C:48–57. 10.1016/j.brainres.2010.03.106
33.
ChenKYangLNLaiCLiuDZhuLQ. Role of GRINA/NMDARA1 in the central nervous system diseases. Curr Neuropharmacol. (2020) 18:861–7. 10.2174/1570159X18666200303104235
34.
BlankeMLVanDongenAMJ. Activation mechanisms of the NMDA receptor. In: Van Dongen AM, ed Biology of the NMDA Receptor. Boca Raton (FL): CRC Press (2009). p. 283–313. 10.1201/9781420044157.ch13
35.
FrigoMCogoMGFuscoMLGardinettiMFrigeniB. Glutamate and multiple sclerosis. Curr Med Chem. (2012) 19:1295–9. 10.2174/092986712799462559
36.
LeviteM. Glutamate, T cells and multiple sclerosis. J Neural Transm. (2017) 124:775–98. 10.1007/s00702-016-1661-z
37.
KetemaMSecadesPKreftMNahidiazarLJanssenHJalinkKet al. The rod domain is not essential for the function of plectin in maintaining tissue integrity. Mol Biol Cell. (2015) 26:2402–17. 10.1091/mbc.E15-01-0043
38.
GostyńskaKBNijenhuisMLemminkHPasHHPasmooijAMGLangKKet al. Mutation in exon 1a of PLEC, leading to disruption of plectin isoform 1a, causes autosomal-recessive skin-only epidermolysis bullosa simplex. Hum Mol Genet. (2015) 24:3155–62. 10.1093/hmg/ddv066
39.
WuerfelJPaulFBeierbachBHamhaberUKlattDPapazoglouSet al. MR-elastography reveals degradation of tissue integrity in multiple sclerosis. Neuroimage. (2010) 49:2520–5. 10.1016/j.neuroimage.2009.06.018
40.
WangXAbrahamSMcKenzieJAGJeffsNSwireMTripathiVBet al. LRG1 promotes angiogenesis by modulating endothelial TGF-β signalling. Nature. (2013) 499:306–11. 10.1038/nature12345
41.
MengHSongYZhuJLiuQLuPYeNet al. LRG1 promotes angiogenesis through upregulating the TGF-β1 pathway in ischemic rat brain. Mol Med Rep. (2016) 14:5535–43. 10.3892/mmr.2016.5925
42.
NocitiVSantoroM. What do we know about the role of lncRNAs in multiple sclerosis?Neural Regen Res. (2021) 16:1715–22. 10.4103/1673-5374.306061
43.
TuckMKChanDWChiaDGodwinAKGrizzleWEKruegerKEet al. Standard operating procedures for serum and plasma collection: early detection research network consensus statement standard operating procedure integration working group. J Proteome Res. (2009) 8:113–7. 10.1021/pr800545q
44.
AkaishiTMisuTFujiharaKNakayaNNakamuraTKogureMet al. White blood cell count profiles in multiple sclerosis during attacks before the initiation of acute and chronic treatments. Sci Rep. (2021) 11:1–9. 10.1038/s41598-021-01942-8
45.
van LangelaarJRijversLSmoldersJvan LuijnMM. B and T cells driving multiple sclerosis: identity, mechanisms and potential triggers. Front Immunol. (2020) 11:1–12. 10.3389/fimmu.2020.00760
Summary
Keywords
multiple sclerosis, vitamin D deficiency, super-enhancers, vitamin D receptor (VDR), vitamin D
Citation
Orton SM, Sangha A, Gupta M, Martens K, Metz LM, de Koning APJ and Pfeffer G (2022) Expression of risk genes linked to vitamin D receptor super-enhancer regions and their association with phenotype severity in multiple sclerosis. Front. Neurol. 13:1064008. doi: 10.3389/fneur.2022.1064008
Received
07 October 2022
Accepted
05 December 2022
Published
28 December 2022
Volume
13 - 2022
Edited by
Davide Cossu, University of Sassari, Italy
Reviewed by
Gloria Dalla Costa, San Raffaele Scientific Institute (IRCCS), Italy; Mohammad Reza Asadi, Tabriz University of Medical Sciences, Iran
Updates
Copyright
© 2022 Orton, Sangha, Gupta, Martens, Metz, de Koning and Pfeffer.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sarah M. Orton ✉sorton@mtroyal.ca
†These authors have contributed equally to this work and share first authorship
This article was submitted to Multiple Sclerosis and Neuroimmunology, a section of the journal Frontiers in Neurology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.