Gene variants of adhesion molecules predispose to MS: A case-control study
Citation Manager Formats
Make Comment
See Comments

Abstract
Objective To examine the effect of variants in genes encoding molecules that are implicated in leukocyte trafficking into the CNS on the development of MS.
Methods A total of 389 Greek MS cases and 336 controls were recruited by 3 MS centers in Cyprus and Greece. In total, 147 tagging single nucleotide polymorphisms across 9 genes encoding for P-selectin (SELP), integrins (ITGA4, ITGB1, and ITGB7), adhesion molecules (ICAM1, VCAM1, and MADCAM1), fibronectin 1 (FN1), and osteopontin (SPP1) were genotyped. The clinical end point of the study was diagnosis of MS according to the 2005 revised McDonald criteria. Permutation analysis was used for adjusting for multiple comparisons.
Results Overall, 21 variants across SELP, ITGA4, ITGB1, ICAM1, VCAM1, MADCAM1, FN1, and SSP1 genes were each associated with MS (pperm < 0.05). The most significant were rs3917779 and rs2076074 (SELP), rs6721763 (ITGA4), and rs1250258 (FN1), all with a permutation p value of less than 1e-004.
Conclusions The current study provides preliminary evidence that variants across genes encoding adhesion molecules, responsible for lymphocyte adhesion and trafficking within the CNS, are implicated in the risk of developing MS.
Glossary
- BBB=
- blood-brain barrier;
- CGAS=
- candidate gene association study;
- DMT=
- disease-modifying therapy;
- EDSS=
- Kurtzke Expanded Disability Status Scale;
- GWAS=
- genome-wide association study;
- LD=
- linkage disequilibrium;
- MAF=
- minor allele frequency;
- MSSS=
- Multiple Sclerosis Severity Score;
- OR=
- odds ratio;
- SELP=
- P-selectin;
- SNP=
- single nucleotide polymorphism
MS is a common complex immune-mediated disease of the CNS.1 Inflammation, demyelination, axonal loss, and neurodegeneration appear during the course of MS.2 Although the proximal pathogenic mechanisms that lead to MS development still remain largely unclear, there is evidence that hints toward an interplay between genetic background, environment factors, and epigenetic modifications conferring susceptibility to MS.3
A few genetic variants that confer susceptibility to MS have been identified through candidate gene association studies (CGASs), linkage and genome-wide linkage analyses, and genome-wide association studies (GWASs).4 Over 110 genetic risk factors of MS have been identified.5 However, most studies have yielded inconsistent results, with limited replication of the reported associations.6
Several molecules have been implicated in the adhesion and diapedesis of leukocytes via the blood-brain barrier (BBB), as well as their migration from the extracellular matrix to the target tissue.7 The protective effect of an antibody against α4 integrin (natalizumab) in MS supports the hypothesis that the ingress of leukocytes into the brain parenchyma is a crucial step.8 We have previously reported that variants of genes coding for proteins implicated in leukocyte trafficking across the BBB influence the severity of MS (namely rs6721763 [ITGA4] and rs6532040 [SPP1]) and have an effect on disease onset (rs1250249 [FN1]).9 Therefore, we also hypothesize that variants of these molecules may also predispose to MS development.
The present study examines tagging single nucleotide polymorphisms (SNPs) in P-selectin (SELP), integrins α4, β1, and β7, adhesion molecules ICAM1, VCAM1, and MADCAM1, fibronectin 1, and osteopontin genes for possible association with MS risk.
Methods
Study population
Our cohort consisted of Greek patients with MS and controls from 3 MS centers: the Cyprus Institute of Neurology and Genetics in Cyprus, the University Hospital of Larissa, Greece, and the AHEPA Hospital of Aristotle University in Thessaloniki, Greece, that have been previously described in detail.9 In short, a total of 389 MS cases and 336 controls have participated in the study. The diagnosis of MS was assessed according to the 2005 revised McDonald criteria.10 All participants were aged ≥18 years and had a disease duration of ≥5 years. We did not include primary progressive MS cases.
The same criteria for collection of clinical and demographic data were applied at each center. Participants' sex, current age, age at disease onset (defined as age of the first neurologic symptom suggestive of MS), use of disease-modifying therapies (DMTs), and duration of treatment were recorded. Using the Kurtzke Expanded Disability Status Scale (EDSS),11 the disability status was assessed in the absence of relapse for at least 6 months. The time from the first symptom suggestive of MS to the last EDSS assessment was considered as the disease duration (in years) for each patient. The progression index was applied to express the rate of disease progression. By means of the Multiple Sclerosis Severity Score (MSSS), disease severity was estimated.12
The unrelated healthy control group was matched for age and sex and consisted of healthy volunteers from the same ethnic regions.
Standard protocol approvals, registrations, and patient consents
The study received approval from the respective institutional ethical standards committee on human experimentation for any experiments using human patients. Written informed consent was obtained from all patients (or guardians of patients) participating in the study (consent for research).
DNA isolation, selection of tagging SNPs, and genotyping procedure
With the puregene DNA purification kit (Qiagen, Valencia, CA), we isolated DNA from peripheral blood samples. Variants across the following genes were genotyped: (1) SELP, (2) ITGA4, (3) ITGB1, (4) ITGB, (5) ICAM1, (6) VCAM1, (7) MADCAM1, (8) FN1, and (9) SPP1.
Tagging SNPs were identified on the basis of linkage disequilibrium (LD) blocks according to the HapMap project. More precisely, the selection was based on the CEU (Utah residents with Northern and Western European ancestry from the CEPH [Centre d'Etude du Polymorphisme Humain collection]) population at the HapMap Release 27 database (Phase II + III, Feb09, on National Center for Biotechnology Information B36 assembly, dbSNP b126). Pairwise r2 values ≥0.8 and minor allele frequency (MAF) >0.05 were the applied criteria. Thirty-two tagging SNPs in the SELP gene, 36 in the ITGA4 gene, 16 in ITGB1, 3 in ITGB7, 5 in ICAM1, 15 in VCAM1, 3 in MADCAM1, 38 in FN1, and 8 in SPP1, which add up to 156 SNPs, were retrieved. Studied tagging SNPs' characteristics and information regarding their genes are available in table e-1 (links.lww.com/NXG/A137).
TaqMan 5′-nuclease assays were used to genotype SNPs by acquiring specific assays for each SNP and then performing allelic discrimination assays on the ABI 7900HT Real-Time System (Applied Biosystems, Foster City, CA). Nine SNPs failed to be genotyped, whereas the genotypic rates were above 96%.
Statistical analysis
The exact test, with a p value of ≤0.05, was considered as indicative for deviation from the Hardy-Weinberg equilibrium in the control group for each SNP.13 Haploview software version 4.2 (broadinstitute.org/scientific-community/science/programs/medical-and-population-genetics/haploview/haploview) was applied in each gene in order for pairwise D′ and r2 between variants and LD blocks to be calculated.14,15 The power of the sample was calculated using CaTS software (csg.sph.umich.edu/abecasis/cats/gas_power_calculator/index.html).16 Quality control with a threshold lower than 0.05% was performed by calculating MAFs. Allelic association regression analysis was performed using PLINK 1.07 software (zzz.bwh.harvard.edu/plink/).17 p values were corrected for multiple testing by means of permutation analysis. A total of 10,000 permutations were performed to correct the p value. Value less than 0.05 was set as the statistically significant threshold.
Data availability statement
Anonymized data will be shared by request from any qualified investigator.
Results
Detailed characteristics of the recruited cohort have been previously extensively described.9 In brief, 251 patients (64.5%) were women. Two hundred seventy-one patients (69.7%) were under DMT. Two hundred thirty-four healthy controls (69.6%) were women. The mean age of the control group was 36.3 years. Detailed clinical and demographic data of controls and patients are presented in table 1.
Demographic and clinical characteristics of study participants
According to the power analysis, our study had a power of 80.0% to detect an association of a variant with a genetic relative risk of 1.729, assuming the multiplicative model, MAF equal to 5% (the lowest in MS cases for the rs3917714), type I error level of 0.05, in a sample of 336 controls and 389 MS cases.
The distribution of the variants that were analyzed in our cohort, in distinct LD blocks across genes, and the pairwise D′ and r2 values between them are demonstrated in figure e-1 (links.lww.com/NXG/A136).
According to allelic association regression analysis, significant associations (with pperm < 0.05) were found for 21 SNPs in total. Namely, 8 SELP polymorphisms (rs3917779, rs3917768, rs2142760, rs3917744, rs2076074, rs3917740, rs3917727, and rs3917709), 7 ITGA4 polymorphisms (rs12988934, rs11694175, rs17224277, rs155103, rs17225354, rs11689738, and rs6721763), 1 ITGB1 polymorphism (rs10763902), 1 ICAM1 polymorphism (rs281437), 1 VCAM1 polymorphism (rs3176878), 1 MADCAM1 polymorphism (rs12982646), 1 FN1 polymorphism (rs1250258), and 1 SSP1 polymorphism (rs6532040) were found to reach a statistical significant level. The most statistically significant pperm values were reached from rs3917779 and rs2076074 (SELP), rs6721763 (ITGA4), and rs1250258 (FN1), with a permutation p value equal to 9.999e-005. The odds ratios [ORs], 95 CIs, and the permutation p values after correction for multiple testing for statistically significant variants (stressed according to the permutation p value) are presented in table 2. The ORs, 95 CIs, p values, and the permutation p values after correction for multiple testing for each variant are presented in table e-2 (links.lww.com/NXG/A137).
Rs number, gene, OR, 95% CI, and permutation p value for the statistical significant variants resulted from allelic association analysis association with MS risk (MS cases vs healthy controls)
Discussion
The present research aimed at investigating a possible association between predisposition to MS and variants of genes that are implicated in leukocyte trafficking into the CNS. The 147 tagging SNPs across 9 genes were selected with the LD block method and further analyzed. Our analysis revealed associations between 21 tagging SNPs across SELP, ITGA4, ITGB1, ICAM1, VCAM1, MADCAM1, FN1, and SSP1 and MS. This result needs further replication in ethnically different and even larger groups.
Recently, a renewed effort to identify genetic variants conferring susceptibility to MS has been made. A few variants have emerged through CGASs and GWASs.5,6,18,–,20 However, CGASs and GWASs have yielded inconsistent findings, with limited replication of genetic variants.6 Previous CGASs have focused on pathways related to autoimmunity, immune function, myelin structure, and the human leukocyte antigen, which also yielded inconsistent results.21,22 It is possible that the complexity of MS pathophysiology,5 the interplay between environmental and genetic risk factors, epistasis between loci, and the lack of systematic exploration of rare variations in MS susceptibility may be responsible, to some degree, for the low reproducibility of MS loci in different studies.6,18
The SELP protein is stored in the alpha-granules of platelets and Weibel-Palade bodies of endothelial cells.23 During platelet activation and degranulation, SELP is redistributed to the plasma membrane and mediates the interaction between leukocytes and activated endothelial cells or platelets.23 Eight SELP polymorphisms (rs3917779, rs3917768, rs2142760, rs3917744, rs2076074, rs3917740, rs3917727, and rs3917709) have reached the statistically significant threshold in our study. The rs3917779, rs3917768, and rs2142760 are in block 2 according to LD values. This region of the SELP gene, that spans approximately 7.8 kb between the 2 extremes (rs3917786 and rs2142760, with a value of D′ equal to 0.986), appears to represent an important part across the SELP gene and may affect the expression of the gene. Every single SNP within this region may confer an additional risk to the development of MS.
The ITGA4 gene encodes the α4 subunit of the α4β1 integrin. This integrin is involved in the adhesion and migration procedures of leukocytes through the BBB.7 ITGA4 is already an approved treatment target in clinical practice in relapsing-remitting MS, with the monoclonal antibody natalizumab.24 A few variants (rs1143676, rs1449263, and rs6721763) in ITGA4 have been linked to MS risk.9,25,–,27 Of interest, the intronic rs6721763 (intron 17–18), which has been found to be associated with MS risk, has also been associated with disease severity, as evaluated with the MSSS.9
Osteopontin consists a secreted matrix protein, which is involved in processes regarding activation, proliferation, and migration of T cells.28 The intronic rs6532040 has been found to confer susceptibility to MS, based on our results. The rs6532040 of SPP1 has been previously reported to be a genetic modifier of MS course.9
Fibronectin 1 is an extracellular matrix glycoprotein. Its main function includes cytoskeleton organization, cell adhesion, and migration.29 Previously, the FN1 rs1250249 has been found to exert a dose-dependent effect on MS onset.9 However, this polymorphism does not seem to be associated with the risk of developing MS in our study. In contrast, the FN1 rs1250258 (pperm = 9.999e-005) is strongly associated with MS risk in our cohort. The rs1250258 is 12.9 kb away from the rs1250249, and in low LD (D′ = 0.048), suggesting that additional loci across FN1 may consist a risk factor for MS.
One polymorphism across each of the ITGB1, ICAM1, VCAM1, and MADCAM1 genes has been associated with MS in our study. Namely, the polymorphisms rs10763902, rs281437, rs3176878, and rs12982646, respectively, have reached pperm values of <0.05. ITGB1, ICAM1, VCAM1, and MADCAM1 are adhesion molecules expressed in the CNS extracellular space.7,9,30 They are mainly implicated in the regulation of the reactivation and migration of lymphocytes to myelinated axons.7,9,30
The clinically well-characterized MS group is one of the strengths of our study. Moreover, according to the number of participants, our sample was sufficiently powered to reveal common susceptibility variants of modest effect. In addition, the correction for multiple testing was assessed by means of permutation analysis, an unbiased method.
However, our study has certain limitations. First, it carries all the inherent limitations of studies with a case-control design. Moreover, our concluding remarks would have been more robust if our analysis was accompanied by supportive functional analyses, a validation in an independent sample or an in silico analysis.
In brief, our study implicates variants of genes encoding adhesion molecules, which are mainly involved in lymphocyte ingress and trafficking in the CNS, in the risk of developing MS. The investigation of additional genetic factors that predispose to MS may provide physicians with personalized tools for diagnosis or treatment, particularly in patients with clinically or radiologically isolated syndromes.31
Study funding
The study was supported by the Project YGEIA/BIOS/0609 (BIE/01), which is cofinanced by the European Regional Development Fund and the Republic of Cyprus through the Research Promotion Foundation.
Disclosure
Efthimios Dardiotis, Elena Panayiotou, Vasileios Siokas, and Athina-Maria Aloizou report no disclosures. Kyproula Christodoulou has received foundation/society research support from the Research Promotion Foundation (Cyprus), Telethon, the Norwegian Financial Mechanism Programme, the Innovative Medicines Initiative, Precisesads, the EU, and Biorise. Andreas Hadjisavvas and Marios Pantzaris report no disclosures. Nikolaos Grigoriadis has served on the scientific advisory boards of Biogen, Bayer, Celgene, Genesis Pharma, Merck Serono, Sanofi, ParadigMS, NavigatingMS, and Helani and has served as a consultant of Biogen, Bayer, Novartis, Celgene, Genesis Pharma, Merck Serono, and Sanofi. Georgios Hadjigeorgiou and Theodoros Kyriakides report no disclosures. Full disclosure form information provided by the authors is available with the full text of this article at Neurology.org/NG.
Appendix 1. Author contributions

Footnotes
Funding information and disclosures are provided at the end of the article. Full disclosure form information provided by the authors is available with the full text of this article at Neurology.org/NG.
The Article Processing Charge was funded by the authors.
- Received August 1, 2018.
- Accepted in final form November 5, 2018.
- Copyright © 2019 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology.
This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND), which permits downloading and sharing the work provided it is properly cited. The work cannot be changed in any way or used commercially without permission from the journal.
References
- 1.↵
- 2.↵
- 3.↵
- Thompson AJ,
- Baranzini SE,
- Geurts J,
- Hemmer B,
- Ciccarelli O
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- 9.↵
- Dardiotis E,
- Panayiotou E,
- Provatas A, et al
- 10.↵
- 11.↵
- 12.↵
- 13.↵
- 14.↵
- 15.↵
- Gabriel SB,
- Schaffner SF,
- Nguyen H, et al
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- Canto E,
- Oksenberg JR
- 20.↵
- Tizaoui K
- 21.↵
- Ramagopalan SV,
- Deluca GC,
- Degenhardt A,
- Ebers GC
- 22.↵
- 23.↵
- Edwards EE,
- Thomas SN
- 24.↵
- 25.↵
- Durmanova V,
- Shawkatova I,
- Javor J, et al
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- Fernández O,
- Delvecchio M,
- Edan G, et al
Letters: Rapid online correspondence
REQUIREMENTS
If you are uploading a letter concerning an article:
You must have updated your disclosures within six months: http://submit.neurology.org
Your co-authors must send a completed Publishing Agreement Form to Neurology Staff (not necessary for the lead/corresponding author as the form below will suffice) before you upload your comment.
If you are responding to a comment that was written about an article you originally authored:
You (and co-authors) do not need to fill out forms or check disclosures as author forms are still valid
and apply to letter.
Submission specifications:
- Submissions must be < 200 words with < 5 references. Reference 1 must be the article on which you are commenting.
- Submissions should not have more than 5 authors. (Exception: original author replies can include all original authors of the article)
- Submit only on articles published within 6 months of issue date.
- Do not be redundant. Read any comments already posted on the article prior to submission.
- Submitted comments are subject to editing and editor review prior to posting.
You May Also be Interested in
Dr. Jeffrey Allen and Dr. Nicholas Purcell
► Watch
Related Articles
- No related articles found.
Topics Discussed
Alert Me
Recommended articles
-
Article
Gene variants of adhesion molecules act as modifiers of disease severity in MSEfthimios Dardiotis, Elena Panayiotou, Antonios Provatas et al.Neurology: Neuroimmunology & Neuroinflammation, April 24, 2017 -
Articles
Methylprednisolone reduces adhesion molecules in blood and cerebrospinal fluid in patients with MSI. Elovaara, M. Lällä, E. Spåre et al.Neurology, December 01, 1998 -
Articles
Adhesion molecule expression in symptomatic and asymptomatic carotid stenosisK. Nuotio, P. J. Lindsberg, O. Carpén et al.Neurology, June 24, 2003 -
Views and Reviews
A neurologist’s guide to genome-wide association studiesS. A. Mullen, D. E. Crompton, P. W. Carney et al.Neurology, February 09, 2009