Analysis of genes, pathways and networks involved in disease severity and age at onset in primary-progressive multiple sclerosis.

BACKGROUND: The role of genetic factors in influencing the clinical expression of multiple sclerosis (MS) is unclear. OBJECTIVE: The objective of this paper is to identify genes, pathways and networks implicated in age at onset (AAO) and severity, measured using the Multiple Sclerosis Severity Score (MSSS), of primary-progressive MS (PPMS). METHODS: We conducted a genome-wide association study (GWAS) of 470 PPMS patients of Italian origin:. Allelic association of 296,589 SNPs with AAO and MSSS was calculated. Pathway and network analyses were also conducted using different tools. RESULTS: No single association signal exceeded genome-wide significance in AAO and MSSS analyses. Nominally associated genes to AAO and MSSS were enriched in both traits for 10 pathways, including: "oxidative phosphorylation" (FDRAAO=9*10(-4); FDRMSSS=3.0*10(-2)), "citrate (TCA) cycle" (FDRAAO=1.6*10(-2); FDRMSSS=3.2*10(-3)), and "B cell receptor signaling" (FDRAAO=3.1*10(-2); FDRMSSS=2.2*10(-3)). In addition, an enrichment of "chemokine signaling pathway" (FDR=9*10(-4)) for AAO and of "leukocyte transendothelial migration" (FDR=2.4*10(-3)) for MSSS trait was observed, among others. Network analysis revealed that p53 and CREB1 were central hubs for AAO and MSSS traits, respectively. CONCLUSIONS: Despite the fact that no major effect signals emerged in the present GWAS, our data suggest that genetic variants acting in the context of oxidative stress and immune dysfunction could modulate the onset and severity of PPMS.

Burden of rare coding variants in an Italian cohort of familial multiple sclerosis.

BACKGROUND: Multiple Sclerosis (MS) is a chronic inflammatory and neurodegenerative demyelinating disease of the central nervous system. It is a complex and heterogeneous disease caused by a combination of genetic and environmental factors, and it can cluster in families. OBJECTIVE: to evaluate at gene-level the aggregate contribution of predicted damaging low-frequency and rare variants to MS risk in multiplex families. METHODS: We performed whole exome sequencing (WES) in 28 multiplex MS families with at least 3 MS cases (81 affected and 42 unaffected relatives) and 38 unrelated healthy controls. A gene-based burden test was then performed, focusing on two sets of candidate genes: i) literature-driven selection and ii) data-driven selection. RESULTS: We identified 11 genes enriched with predicted damaging low-frequency and rare variants in MS compared to healthy individuals. Among them, UBR2 and DST were the two genes with the strongest enrichment (p = 5 × 10-4 and 3 × 10-4, respectively); interestingly enough the association signal in UBR2 is driven by rs62414610, which was present in 25% of analysed families. CONCLUSION: Despite limitations, this is one of the first studies evaluating the aggregate contribution of predicted damaging low-frequency and rare variants in MS families using WES data. A replication effort in independent cohorts is warranted to validate our findings and to evaluate the role of identified genes in MS pathogenesis.

HLA-class I markers and multiple sclerosis susceptibility in the Italian population.

Previous studies reported an association with multiple sclerosis (MS) of distinct HLA-class I markers, namely HLA-A*02, HLA-Cw*05 and MOG-142L. In this work, we tested the association with MS of A*02 and Cw*05 in 1273 Italian MS patients and 1075 matched controls, which were previously analyzed for MOG-142, and explored the relationship among these three markers in modulating MS risk. HLA-A*02 conferred a statistically robust MS protection (odds ratio, OR=0.61; 95% confidence intervals, CI=0.51-0.72, P<10(-9)), which was independent of DRB1*15 and of any other DRB1* allele and remained similar after accounting for the other two analyzed class I markers. Conversely, the protective effect we previously observed for MOG-142L was secondary to its linkage disequilibrium with A*02. Cw*05 was not associated considering the whole sample, but its presence significantly enhanced the protection in the HLA-A*02-positive group, independently of DRB1: the OR conferred by A*02 in Cw*05-positive individuals (0.22, 95% CI=0.13-0.38) was significantly lower than in Cw*05-negative individuals (0.69, 95% CI=0.58-0.83) with a significant (P=4.94 x 10(-5)) multiplicative interaction between the two markers. In the absence of A*02, Cw*05 behaved as a risk factor, particularly in combination with DRB1*03 (OR=3.89, P=0.0006), indicating that Cw*05 might be a marker of protective or risk haplotypes, respectively.

The osteopontin gene +1239A/C single nucleotide polymorphism is associated with type 1 diabetes mellitus in the Italian population.

Secreted phosphoprotein 1, also known as Osteopontin (Opn), is a proinflammatory cytokine involved in the TH1 response and is highly expressed in the islets and pancreatic lymph nodes of non-obese diabetic mice before the onset of diabetes. In humans, typing of the +1239A/C single nucleotide polymorphism (SNP) in the 3UTR of the Opn gene (SPP1) showed that +1239C carriers displayed higher Opn serum levels than +1239A homozygotes and a higher risk of developing autoimmune/lymphoproliferative syndrome, multiple sclerosis, and systemic lupus erythematosus. The aim of this work is to evaluate whether +1239A/C is also associated with type 1 diabetes mellitus (T1DM). We typed +1239A/C in an initial cohort of 184 T1DM patients and 361 controls, and confirmed our data in a second cohort of 513 patients and 857 controls. In both cohorts, +1239C carriers displayed a significantly higher risk of T1DM than +1239A homozygotes (combined cohorts: OR=1.63, 95 percent CI: 1.34-1.97). Clinical analysis did not detect any differences between patients carrying or not +1239C in terms of gender distribution and age at T1DM diagnosis. These data suggest that SPP1 variants marked by +1239C are associated with T1DM development in the Italian population. The predisposing effect may depend on its effect on Opn levels.

A sequence variation in the MOG gene is involved in multiple sclerosis susceptibility in Italy.

Several studies suggest that the histocompatibility complex (HLA) class I region harbours genes modulating multiple sclerosis (MS) susceptibility independently from the effect of class II alleles. A candidate gene in this region is MOG, encoding the myelin oligodendrocyte glycoprotein. A significant association with the missense variation V142L (rs2857766) was previously reported in a small sample of 50 Italian MS patients. We confirmed this result in two independent Italian sample sets consisting of 878 MS patients and 890 matched controls (P=6.6 x 10(-4)) and 246 trio families (P=1.5 x 10(-3)). The comparison of genotype frequencies suggested a dominant-protective effect of L142. In the combined sample sets L142 conferred an odds ratio (OR)=0.70 (95% confidence interval (CI): 0.60-0.82) that remained similar after accounting for HLA-DRB1(*)15 carrier status. The association with MOG V142L was still significant after conditioning for all DRB1 alleles (P=0.035). Eleven additional single nucleotide polymorphisms in the MOG gene (namely -1077T/C, -910T/C, -875A/G, -93T/C, S5S, Indel L22, V145I, +814C/T, +900A/G, +1024A/T, +1059C/T), two microsatellites in the MOG 5' flanking (MOGCA) and 3' untranslated (MOGTAAA) regions and four microsatellites in the HLA-class I region, from HLA-B to HFE, (namely MIB, D6S265, D6S1683 and D6S2239) were tested by transmission disequilibrium test in 199 trio families. None of these polymorphisms or of their haplotypic combinations showed a significant transmission distortion, in the absence of V142L. In conclusion, MOG V142L, or an untested variant in tight-linkage disequilibrium with it, is an independent MS susceptibility-modulating factor in the HLA class I region.

Variations of the perforin gene in patients with multiple sclerosis.

Perforin is involved in cell-mediated cytotoxicity and mutations of its gene (PRF1) cause familial hemophagocytic lymphohistiocytosis (FLH2). PRF1 sequencing in 190 patients with multiple sclerosis and 268 controls detected two FLH2-associated variations (A91V, N252S) in both groups and six novel mutations (C999T, G1065A, G1428A, A1620G, G719A, C1069T) in patients. All together, carriers of these variations were more frequent in patients than in controls (phenotype frequency: 17 vs 9%, P=0.0166; odds ratio (OR)=2.06, 95% confidence interval (CI): 1.13-3.77). Although A91V was the most frequent variation and displayed a trend of association with multiple sclerosis (MS) in the first population of patients and controls (frequency of the 91V allele: 0.076 vs 0.043, P=0.044), we used it as a marker to confirm PRF1 involvement in MS and assessed its frequency in a second population of 966 patients and 1520 controls. Frequency of the 91V allele was significantly higher in patients than in controls also in the second population (0.075 vs 0.058%, P=0.019). In the combined cohorts of 1156 patients and 1788 controls, presence of the 91V allele in single or double dose conferred an OR=1.38 (95% CI=1.10-1.74). These data suggest that A91V and possibly other perforin variations indicate susceptibility to MS.

Opposed independent effects and epistasis in the complex association of IRF5 to SLE.

Genetic variation in the interferon regulatory factor 5 (IRF5) gene affects systemic lupus erythematosus (SLE) susceptibility. However, association is complex and incompletely defined. We obtained fourteen European sample collections with a total of 1383 SLE patients and 1614 controls to better define the role of the different IRF5 variants. Eleven polymorphisms were studied, including nine tag single nucleotide polymorphisms (SNPs) and two extra functional polymorphisms. Two tag SNPs showed independent and opposed associations: susceptibility (rs10488631, P<10(-17)) and protection (rs729302, P<10(-6)). Haplotype analyses showed that the susceptibility haplotype, identified by the minor allele of rs10488631, can be due to epistasis between three IRF5 functional polymorphisms. These polymorphisms determine increased mRNA expression, a splice variant with a different exon 1 and a longer proline-rich region in exon 6. This result is striking as none of the three polymorphisms had an independent effect on their own. Protection was independent of these polymorphisms and seemed to reside in the 5' side of the gene. In conclusion, our results help to understand the role of the IRF5 locus in SLE susceptibility by clearly separating protection from susceptibility as caused by independent polymorphisms. In addition, we have found evidence for epistasis between known functional polymorphisms for the susceptibility effect.

Two single-nucleotide polymorphisms in the 5' and 3' ends of the osteopontin gene contribute to susceptibility to systemic lupus erythematosus.

OBJECTIVE: To test the association of osteopontin (OPN) polymorphisms with systemic lupus erythematosus (SLE). METHODS: The coding 5' and 3' flanking regions of the OPN gene were scanned for polymorphisms by denaturing high-performance liquid chromatography. A case-control association study was performed in 394 Italian SLE patients and 479 matched controls. OPN serum levels were determined by enzyme-linked immunosorbent assay in 40 patients and 124 controls, and the mean levels were compared between the different OPN genotypes. RESULTS: Among the 13 detected single-nucleotide polymorphisms (SNPs), alleles -156G (frequency 0.714 versus 0.651; P = 0.006, corrected P [P(corr)] = 0.036) and +1239C (0.377 versus 0.297; P = 0.00094, P(corr) = 0.0056) were significantly increased in the SLE patients compared with the controls. The presence of the associated allele in single or double dose conferred an odds ratio (OR) of 2.35 (95% confidence interval [95% CI] 1.38-4.02) for SNP -156 and an OR of 1.57 (95% CI 1.16-2.13) for SNP +1239. These effects were independent of each other, i.e., not a consequence of linkage disequilibrium between the 2 alleles. The risk associated with a double dose of susceptibility alleles at both SNPs was 3.8-fold higher (95% CI 2.0-7.4) relative to the complete absence of susceptibility alleles. With regard to individual clinical and immunologic features, a significant association was seen between lymphadenopathy and -156 genotypes (overall P = 0.0011, P(corr) = 0.046). A significantly increased OPN serum level was detected in healthy individuals carrying +1239C (P = 0.002), which is indicative of an association between the SLE susceptibility allele and OPN levels. CONCLUSION: These data suggest the independent effect of a promoter (-156) and a 3'-untranslated region (+1239) SNP in SLE susceptibility. We can speculate that these sequence variants (or others in perfect linkage disequilibrium) create a predisposition to high production of OPN, and that this in turn may confer susceptibility to SLE.

Association tests with systemic lupus erythematosus (SLE) of IL10 markers indicate a direct involvement of a CA repeat in the 5' regulatory region.

Many lines of evidence suggest that IL10 is a strong candidate gene for systemic lupus erythematosus (SLE) susceptibility. In our previously reported study an allele (IL10.G-140bp) of the microsatellite IL10.G located at position -1100 was significantly increased in Italian SLE patients in comparison with controls. Starting from this observation, we tested if sequence variations in the vicinity of IL10.G were more strongly associated with SLE. We performed a comprehensive association study including 26 SNPs (of which four were newly identified in the present study by DHPLC analysis) spanning 8.5 Kb of the 5' flanking and the transcribed region of the IL10 gene. The association study was performed by the DNA pool method on an extended panel of Italian patients (205) and controls (631). Haplotypic associations were studied by individual typing of seven selected markers surrounding IL10.G. Gene, genotype and haplotype frequencies were not significantly different in patients and controls. Thus the IL10.G microsatellite remains to date the only IL10 marker associated with SLE in our population. A meta-analysis of all published results indicates a possible direct role of the IL10.G repeat number in SLE susceptibility.