LINGO1 rs9652490 and rs11856808 polymorphisms are not associated with risk for multiple sclerosis

Background Some recent experimental data suggest a possible role of LINGO-1 in the pathogenesis of multiple sclerosis (MS). In an attempt to identify genetic biomarkers related to MS susceptibility, we genotyped two common SNPs in the LINGO1 gene which have been associated to other neurological conditions, in patients with MS and in healthy subjects. These SNPs are linked to several SNPs within the LINGO1 gene, especially in individuals of Oriental or Caucasian descent. Methods We analyzed the allelic and genotype frequency of two LINGO1 variants (rs9652490 and rs11856808) in 293 patients with MS and 318 healthy controls, using KASPar assays. Results LINGO1 rs9652490 and rs11856808 allelic and genotype frequencies did not differ significantly between MS patients and controls. The minor allele frequencies for rs9652490 were 0.171 (95% CI = 0.140-0.201) and 0.167 (95% CI = 0.138-0.196 for cases and controls respectively (p = 0.853). For rs11856808 the minor allele frequencies were 0.317 (95% CI = 0.280-0.355) and 0.310 (95% CI = 0.274-0.346) for cases and controls, respectively (p = 0.773). Allele and genotype frequencies were unrelated with the age of onset of MS, gender, and clinical course of MS. In addition, haplotype analyses did not reveal any putative risk related to haplotypes. Conclusions These results suggest that LINGO1 rs9652490 and rs11856808 polymorphisms are not related with risk for MS. This study adds to other published evidence indicating that, to date, the LINGO1 SNPs studied here could be useful risk biomarkers of developing essential tremor, but not other movement disorders.


Background
Multiple sclerosis (MS) is a chronic inflammatory demyelinating disorder with axonal degeneration affecting the Central Nervous system, which shows three major evolutive phenotypes: relapsing-remitting, primary progressive and secondary progressive. The etiology of MS is unknown, but it is probably multifactorial, with an interplay of genetic, ethnic, geographical and environmental factors (infectious or chemical) [1][2][3][4][5]. It has been proposed that MS is an autoimmune disorder with susceptibility influenced, if not determined, by a relatively small number of genes [1]. Findings from studies on seasonality in MS patients' birth, disease onset and exacerbations, as well as apparent temporal trends in incidence and gender ratio support an influential effect of viruses, metabolic and lifestyle factors on MS risk. Epstein-Barr virus, vitamin D status, and smoking are factors that may explain such epidemiological patterns [4].
A haplotype within the major histocompatibility region is the major risk factor for MS. But despite clear evidence for a genetic component additional risk, specific gene variants were not identified until the recent advent of genome-wide association studies (GWAS). Until 2010, 11 GWAS have been conducted on MS, and, together with follow-up studies, these GWAS have confirmed 16 loci with genome-wide significance [6,7]. Many of these common risk variants are located at, or near to, genes with central immunological functions (such as interleukin 2 and 7 receptors, CD58, CD6, CD40, TNFRSF1A and others) and the majority are associated with other autoimmune diseases [6,7]. A further report of the International Multiple Sclerosis Genetics Consortium identified at least 50 loci related with the risk for MS [8].
Although the underlying molecular mechanisms for the axonal degeneration are unknown, the degree of inflammatory demyelination correlates with the extent of axonal damage. This suggests an involvement of the proinflammatory mediators in inducing axonal degeneration [9]. However, the alternative possibility that axonal regeneration should be severely impaired in MS lesions could be suggested, since an accumulation of glial scar and neurite growth inhibitors provide a non-permissive environment for re-growth of damaged axons [10].
LINGO1 (leucine rich repeat and Ig domain containing Nogo receptor interacting protein-1) has a possible role in the pathogenesis of MS. LINGO1 is a transmembrane protein expressed in neural cells which inhibits the differentiation of oligodendrocyte precursor cells into mature oligodendrocytes, as well as myelination and remyelination [11,12]. LINGO1 comprises 12 leucine rich repeats followed by an immunoglobulin (Ig) domain and a short cytoplasmic tail. It is encoded by the LINGO1 gene (OMIM 609791, Gene Identity 84894) located in the chromosome 15q24.3 [13,14]. In neurons, LINGO1 simultaneously interacts with the Nogo-66 receptor (NgR) and p75 NTR or TROY to form a receptor complex that binds the structurally diverse associated glycoprotein and oligodendrocyte myelin glycoprotein, resulting in the restriction of axonal elongation via activation of the small GTPase RhoA [14][15][16]. Two LINGO1 variants designated as rs9652490 and rs11856808 have been claimed to be associated in case-control GWAS with other neurological conditions such as essential tremor [17,18] and Parkinson's disease [18,19]. Further studies confirmed the association with essential tremor, but discarded a major association with Parkinson's disease [20][21][22][23][24]. These single nucleotide polymorphisms (SNPs) are, according to HapMap, tag-SNPs for the following SNPs located within the LINGO1 gene: rs907400, rs8029432, rs1877294, rs7165679, rs9920101 and rs9920127, as well as nine additional SNPs in the 3' flanking region of the gene. Figure 1 shows that the linkage between the two SNPs analyzed in this study and the six SNPs located within the gene differ, depending on ethnicity.
In an attempt to identify additional factors involved in MS susceptibility, we genotyped the SNPs rs9652490 and rs11856808 in the LINGO1 gene, in patients with MS and in healthy subjects. Although LINGO1 polymorphisms were not significantly associated with the risk and hence are not mentioned among the possible susceptibility genes in GWAS studies, the possible role of LINGO1 in the pathogenesis of MS suggests that the LINGO1 gene should be a candidate gene for modifying MS risk.

Patients and controls
We recruited 293 unrelated Caucasian Spanish patients who fulfilled McDonald's criteria for definite MS [25], with no other previous neurological diseases. Recruiting sources were the following: the "Multiple Sclerosis Association of Madrid"; n = 165 cases), the Health Areas of the Hospital La-Mancha-Centro (Alcázar de San Juan, Ciudad Real; n = 65 cases), and University Hospitals "Doce de Octubre" (Madrid, n = 32 cases), and "Príncipe de Asturias" (Alcalá de Henares, Madrid; n = 31 cases). The control group was composed of 318 healthy unrelated Caucasian Spanish individuals gender and agematched with the patients (97 men, 221 women; mean age 43.76 ± 12.4 years). These patients participated in previous genetic studies [26][27][28]. The control individuals were students or professors from the University of Extremadura, Badajoz, Spain (n = 150), and the healthy spouses of patients with neurological disorders who came from different regions of Spain to the Department of Neurology, Clínica Universitaria de Navarra, University of Navarra School of Medicine, Pamplona, Spain (n = 168). All the participants were included in the study after giving written informed consent. Table 1 summarizes the characteristics of the individuals included in the study. The protocol was approved by the Ethics Committees of the University Hospitals "Príncipe de Asturias" and "Infanta Cristina" (Badajoz) and collaborating centres. The study was conducted according to the principles expressed in the declaration of Helsinki.

Genotyping of LINGO1 rs9652490 and rs18856808 polymorphisms
Genomic DNA was obtained from peripheral leukocytes and purified according to standard procedures. Two polymorphisms of LINGO1 gene, rs9652490 A/G and rs11856808 C/T, were genotyped using KASPar assays according to the manufacturer's protocol (www.kbioscience. co.uk). PCR was performed on a 96-well Tetrad 2 Peltier Thermal Cycler (BIO-RAD, Hercules, CA). PCR KASPar's protocol was performed as following: a denaturation step of 10 min, twenty-eight cycles of 15 sec denaturing at 94°C, annealing of 20 sec at 57°C, and extension of 30 sec at 72°C. PCR were followed by a final extension step of 5 min at 72°C. Genotype calling was performed in an allelic discrimination analysis module of the 7300   [31]. We used the default model for recombination rate variation with 1000 iterations, 500 burn-in iterations and a thinning interval of 1 as described elsewhere [32].

Results
The frequencies of LINGO1 rs9652490 and rs11856808 genotypes and alleles in patients with MS did not differ from those of controls ( Table 2). The genotype and allele frequencies in MS patients and healthy subjects were in Hardy-Weinberg's equilibrium. Both SNPs were in linkage disequilibrium with LOD values equal to 48.67 and 36.17 among patients and controls subjects, respectively (r-squared = 0.46 and 0.48, respectively). Mean age at onset of MS did not differ significantly between patients carrying LINGO1 rs9652490 A/A (mean ± SD = 33.0 ± 11.5 years), A/G (mean ± SD = 32.4 ± 9.8 years) and G/G (mean ± SD = 38.5 ± 24.7 years); (p = 0.818 for the comparison of carriers vs. non-carriers of variant alleles), and between patients with genotypes LINGO1 rs11856808 C/C (mean ± SD = 32.5 ± 11.5 years), C/T (mean ± SD = 32.6 ± 10,6 years) and T/T (mean ± SD = 35.6 ± 10.9 years); (p = 0.835 for the comparison of carriers vs. non-carriers of variant alleles). The distribution of rs9652490 and rs11856808 allelic and genotype frequencies were not influenced by gender ( Table 2). The distribution of the LINGO1 rs9652490 and rs11856808 genotype and allelic frequencies did not differ between each MS phenotype and controls (Table 3) or in the severity scores: expanded disability status scale or progression index. Haplotype analyses indicated that the commonest rs9652490-rs11856808 haplotype was

Discussion
The possible role of LINGO1 in the pathogenesis of MS makes it reasonable to analyse the possible relationship between LINGO1 polymorphisms and the risk of MS. In the present study, we found no significant differences in allele genotypes, or haplotypes frequencies for the rs9652490 and rs11856808 polymorphisms when comparing patients with MS and healthy control subjects. Nor were these polymorphisms related with the age at onset of MS or with the evolutive type of MS. The findings obtained, though negative, are novel and represent an incremental advance in the knowledge of the clinical implications of the LINGO1 gene polymorphism. Some experimental data suggest a possible role of LINGO1 in the pathogenesis of MS: (a) Nogo-A expression has been found to be enhanced in surviving oligodendrocytes, while NgR has been found to be up-regulated in reactive astrocytes and macrophages/microglia in chronic active demyelinating lesions of MS [33], (b) TROY has been found to be up-regulated, whereas LINGO1 expression has been found to be reduced, in MS brains [34], (c) LINGO1 knockout mice have shown earlier onset of myelination of CNS axons than the wild-type, and greater resistance to the development of myelin oligodendrocyte glycoprotein (MOG)-induced experimental autoimmune encephalomyelitis (EAE) [35]; (d) Treatment with antibody antagonists against LINGO1 function induces functional recovery and increases integrity of axons in MOG-induced EAE [35], and promotes oligodendrocyte precursor cell differentiation and remyelination in different experimental models of demyelination and remyelination [36].
The LINGO1 SNPs analyzed in this study have been studied as putative risk biomarkers for other movement disorders. A recent meta-analysis, which included 3,972 essential tremor (ET) patients and 20,714 controls for the LINGO1 rs9652490 polymorphism, and 2,076 ET patients and 18,792 controls for the rs11856808 polymorphism, concluded that the rs11856808 polymorphism was related with increased risk for both total and familial ET, whereas the rs9652490 polymorphism was related with increased risk for familial ET [21]. With regard to Parkinson's disease, another recent metaanalysis including 5,541 patients and 5,647 controls for the rs9652490 polymorphism and 3,276 patients and 3,371 controls for the rs11856808 polymorphism concluded that these polymorphisms could not be considered as major risk factors for susceptibility to PD [22]. Several SNPs have been described within the LINGO1 gene and in the 3' flanking region (Figure 1). It is to be noted that most genetic association studies on LINGO1 focused on the same SNPs which were analyzed in the present study [17][18][19][21][22][23][24][37][38][39][40][41][42][43][44][45][46][47][48][49][50]. Although several LIN GO1 nonsynonymous SNPs have been described, namely rs113329801, rs201732477, rs112205560, rs150289554, rs113096707, rs188738703, rs200688402, rs200463885, rs9855, rs201517725, rs184237450, rs77436810, rs1931 00227, rs111741384, rs202233236, rs199976207, rs199 628078, rs201438433, rs140914739, rs200528664, rs111 605415 and rs202223502, none of these SNPs show a minor allele frequency over 0.0005 (see the website http://www.ncbi.nlm.nih.gov/projects/SNP/snp_ref.cgi? showRare=on&chooseRs=all&go=Go&locusId=84894). Therefore, the identification of one heterozygous individual out of 1000 individuals studied is to be expected. This precluded the analysis of these nonsynonymous SNPs as putative risk factors. Nevertheless, the two SNPs studied are linked with other SNPs within the coding region of the LINGO1 gene. Figure 1  The present study has some limitations. First, the size of analyzed cohorts may not be sufficient for strict conclusions about LINGO1 role in MS. As was shown in a previous publication about the role of LINGO1 in essential tremor risk, individual studies of small number of patients gave very contradictory results [21]. Second, although the sample size is adequate to detect an OR as small as 1.5, a more modest association would not be detected. Third, because EDSS or progression index are not completely adequate measures of disease severity, the negative association observed in this study does not rule out a putative association with disease severity. Moreover, because the cohort study included MS patients with different degrees of severity, it is not adequate for the investigation of the influence of LINGO1 genotypes on the disability or severity of MS (the ideal study for this purpose should include genotyping of patients with a recent diagnosis of MS with similar followup periods).

Conclusions
In summary, taking in account the limitations of the present study, our results suggest that rs9652490 and rs11856808 genotype and allelic variants are not related with the risk for MS in Caucasian Spanish people.

Competing interests
The authors declare that they have no competing interests.
Authors' contributions EGM participated in the conception and design of the study, acquisition of data, analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript, administrative, technical, and material support, supervision, and obtaining funding. OLB participated in acquisition of data, analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript, administrative, technical, and material support. CM participated in acquisition of data, analysis and interpretation of data, critical revision of the manuscript, administrative, technical, and material support. PP participated in the conception and design of the study, acquisition of data, analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript, administrative, technical, and material support, supervision, and obtaining funding. JBL participated in acquisition of data, and critical revision of the manuscript. JM participated in acquisition of data, and critical revision of the manuscript. PC participated in acquisition of data, and critical revision of the manuscript. MDS participated in acquisition of data, and critical revision of the manuscript. DP participated in acquisition of data, and critical revision of the manuscript. LTF participated in acquisition of data, and critical revision of the manuscript. HAN participated in acquisition of data, analysis and interpretation of data, critical revision of the manuscript, administrative, technical, and material support. LAP participated in acquisition of data, and critical revision of the manuscript. DT participated in acquisition of data, and critical revision of the manuscript. EL participated in acquisition of data, and critical revision of the manuscript. JFPN participated in acquisition of data, and critical revision of the manuscript. JAGA participated in the conception and design of the study, acquisition of data, analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript, administrative, technical, and material support, supervision, and obtaining funding. FJJJ participated in the conception and design of the study, acquisition of data, analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript, administrative, technical, and material support, and supervision. All authors read and approved the final manuscript.