Ferdinando Clarelli

The burden of multiple sclerosis variants in continental Italians and Sardinians

Background: Recent studies identified > 100 non-HLA (human leukocyte antigen) multiple sclerosis (MS) susceptibility variants in Northern European populations, but their role in Southern Europeans is largely unexplored. Objective: We aimed to investigate the cumulative impact of those variants in two Mediterranean populations: Continental Italians and Sardinians. Methods: We calculated four weighted Genetic Risk Scores (wGRS), using up to 102 non-HLA MS risk variants and 5 HLA MS susceptibility markers in 1691 patients and 2194 controls from continental Italy; and 2861 patients and 3034 controls from Sardinia. We then assessed the differences between populations using Nagelkerke’s R2 and the area under the Receiver Operating Characteristic (ROC) curves. Results: As expected, the genetic burden (mean wGRS value) was significantly higher in MS patients than in controls, in both populations. Of note, the burden was significantly higher in Sardinians. Conversely, the proportion of variability explained and the predictive power were significantly higher in continental Italians. Notably, within the Sardinian patients, we also observed a significantly higher burden of non-HLA variants in individuals who do not carry HLA risk alleles. Conclusions: The observed differences in MS genetic burden between the two Mediterranean populations highlight the need for more genetic studies in South Europeans, to further expand the knowledge of MS genetics.

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.

Association between DPP6 polymorphism and the risk of progressive multiple sclerosis in Northern and Southern Europeans

BACKGROUND In this study, we investigated the role of the dipeptidyl-peptidase-6 (DPP6) gene in the etiopathogenesis of progressive forms of multiple sclerosis (PrMS). This gene emerged as a candidate gene in a genome-wide association study (GWAS) performed in an Italian sample of PrMS and controls in which two SNPs located in the gene (rs6956703 and rs11767658) showed evidence of association (nominal p-value<10(-4)) (Martinelli-Boneschi et al.) [18]. Moreover, the gene is highly expressed in the central nervous system, and it has been found to be associated with sporadic cases of amyotrophic lateral sclerosis which shares some feature with PrMS. METHODS We genotyped 19 SNPs selected using a direct and tagging approach in 244 Italian PrMS and 225 controls, and we measured the expression levels of the gene in 13 PrMS cases and 25 controls. RESULTS Five out of 19 SNPs were found to be associated with the disease (adjusted p<0.05), and they have been tested in an independent sample of 179 primary progressive MS and 198 controls from Northern Europe. None of the SNPs was replicated, but combined analysis confirmed the presence of association for rs2046748 (p=2.5×10(-3),OR=1.82, 95%CI=1.24-2.69). CONCLUSIONS These results, inflated by the limited sample size determined by the rarity of this condition, suggest a possible role of this gene in the susceptibility to PrMS, at least in Southern Europeans. Moreover, DPP6 was over-expressed in PrMS patients compared to controls.

Next Generation Sequencing of Pooled Samples: Guideline for Variants’ Filtering

Sequencing large number of individuals, which is often needed for population genetics studies, is still economically challenging despite falling costs of Next Generation Sequencing (NGS). Pool-seq is an alternative cost- and time-effective option in which DNA from several individuals is pooled for sequencing. However, pooling of DNA creates new problems and challenges for accurate variant call and allele frequency (AF) estimation. In particular, sequencing errors confound with the alleles present at low frequency in the pools possibly giving rise to false positive variants. We sequenced 996 individuals in 83 pools (12 individuals/pool) in a targeted re-sequencing experiment. We show that Pool-seq AFs are robust and reliable by comparing them with public variant databases and in-house SNP-genotyping data of individual subjects of pools. Furthermore, we propose a simple filtering guideline for the removal of spurious variants based on the Kolmogorov-Smirnov statistical test. We experimentally validated our filters by comparing Pool-seq to individual sequencing data showing that the filters remove most of the false variants while retaining majority of true variants. The proposed guideline is fairly generic in nature and could be easily applied in other Pool-seq experiments.

Pharmacogenetic study of long-term response to interferon-β treatment in multiple sclerosis.

The aim of the study is the identification of genetic factors that influence the long-term response to interferon-β (IFNβ) (4-year follow-up). We performed a genome-wide association study in 337 IFNβ-treated Italian multiple sclerosis patients at the extreme of treatment response, and we meta-analyzed association effects, integrating results with pathway analysis, gene-expression profiling of IFNβ-stimulated peripheral blood mononuclear cells from 20 healthy controls (HC) and expression quantitative locus (eQTL) analyses. From meta-analysis, 43 markers were associated at P<10−4, and two of them (rs7298096 and rs4726460) pointed to two genes, NINJ2 and TBXAS1, that were significantly downregulated after IFNβ stimulation in HC (P=3.1 × 10−9 and 5.6 × 10−10). We also observed an eQTL effect for the allele associated with favorable treatment response (rs4726460A); moreover, TBXAS1 appeared downregulated upon IFNβ administration (β=−0.39; P=0.02). Finally, we found an enrichment of pathways related to inflammatory processes and presynaptic membrane, the latter with involvement of genes related to glutamatergic system (GRM3 and GRIK2), confirming its potential role in the response to IFNβ.

Impact of MS genetic loci on familial aggregation, clinical phenotype, and disease prediction.

Objective: To investigate the role of known multiple sclerosis (MS)-associated genetic variants in MS familial aggregation, clinical expression, and accuracy of disease prediction in sporadic and familial cases. Methods: A total of 1,443 consecutive patients were screened for MS and familial autoimmune history in a hospital-based Italian cohort. Among them, 461 sporadic and 93 familial probands were genotyped for 107 MS-associated polymorphisms. Their effect sizes were combined to calculate the weighted genetic risk score (wGRS). Results: Family history of MS was reported by 17.2% of probands, and 33.8% reported a familial autoimmune disorder, with autoimmune thyroiditis and psoriasis being the most frequent. No difference in wGRS was observed between sporadic and familial MS cases. In contrast, a lower wGRS was observed in probands with greater familial aggregation (>1 first-degree relative or >2 relatives with MS) (p = 0.03). Also, female probands of familial cases with greater familial aggregation had a lower wGRS than sporadic cases (p = 0.0009) and male probands of familial cases (p = 0.04). An inverse correlation between wGRS and age at onset was observed (p = 0.05). The predictive performance of the genetic model including all known MS variants was modest but greater in sporadic vs familial cases (area under the curve = 0.63 and 0.57). Conclusions: Additional variants outside the known MS-associated loci, rare variants, and/or environmental factors may explain disease occurrence within families; in females, hormonal and epigenetic factors probably have a predominant role in explaining familial aggregation. The inclusion of these additional factors in future versions of aggregated genetic measures could improve their predictive ability.

Inverse correlation of genetic risk score with age at onset in bout-onset and progressive-onset multiple sclerosis.

We correlated the weighted genetic risk score measured using 107 established susceptibility variants for multiple sclerosis (MS) with the age at onset in bout-onset (BOMS, n=906) and progressive-onset MS Italian patients (PrMS) (n=544). We observed an opposite relationship in the two disease courses: a higher weighted genetic risk score was associated with an earlier age at onset in BOMS (rho= −0.1; p=5 × 10−3) and a later age at onset in PrMS cases (rho=0.07; p=0.15) (p of difference of regression=1.4 × 10−2). These findings suggest that established MS risk variants anticipate the onset of the inflammatory phase, while they have no impact on, or even delay, the onset of the progressive phase.

Genetic burden of common variants in progressive and bout-onset multiple sclerosis

Background: The contribution of genetic variants underlying the susceptibility to different clinical courses of multiple sclerosis (MS) is still unclear. Objective: The aim of the study is to evaluate and compare the proportion of liability explained by common SNPs and the genetic burden of MS-associated SNPs in progressive onset (PrMS) and bout-onset (BOMS) cases. Methods: We estimated the proportion of variance in disease liability explained by 296,391 autosomal SNPs in cohorts of Italian PrMS and BOMS patients using the genome-wide complex trait analysis (GCTA) tool, and we calculated a weighted genetic risk score (wGRS) based on the known MS-associated loci. Results: Our results identified that common SNPs explain a greater proportion of phenotypic variance in BOMS (36.5%±10.1%) than PrMS (20.8%±6.0%) cases, and a trend of decrease was observed when testing primary progressive (PPMS) without brain MRI inflammatory activity (p = 7.9 × 10−3). Similarly, the wGRS and the variance explained by MS-associated SNPs were higher in BOMS than PPMS in males (wGRS: 6.63 vs 6.51, p = 0.04; explained variance: 4.8%±1.5% vs 1.7%±0.6%; p = 0.05). Conclusions: Our results suggest that the liability of disease is better captured by common genetic variants in BOMS than PrMS cases. The absence of inflammatory activity and male gender further raise the difference between clinical courses.

Response to interferon-beta treatment in multiple sclerosis patients: A genome-wide association study.

Up to 50% of multiple sclerosis (MS) patients do not respond to interferon-beta (IFN-β) treatment and determination of response requires lengthy clinical follow-up of up to 2 years. Response predictive genetic markers would significantly improve disease management. We aimed to identify IFN-β treatment response genetic marker(s) by performing a two-stage genome-wide association study (GWAS). The GWAS was carried out using data from 151 Australian MS patients from the ANZgene/WTCCC2 MS susceptibility GWAS (responder (R)=51, intermediate responders=24 and non-responders (NR)=76). Of the single-nucleotide polymorphisms (SNP) that were validated in an independent group of 479 IFN-β-treated MS patients from Australia, Spain and Italy (R=273 and NR=206), eight showed evidence of association with treatment response. Among the replicated associations, the strongest was observed for FHIT (Fragile Histidine Triad; combined P-value 6.74 × 10−6) and followed by variants in GAPVD1 (GTPase activating protein and VPS9 domains 1; combined P-value 5.83 × 10−5) and near ZNF697 (combined P-value 8.15 × 10−5).

CHRNA7 Gene and Response to Cholinesterase Inhibitors in an Italian Cohort of Alzheimer's Disease Patients.

Previous studies suggest that genetic variants in CHRNA7, which encodes for the major subunit of the acetylcholine receptor (α7-nAChR), are associated with the clinical response to cholinesterase inhibitors (ChEI) in Alzheimers disease (AD) patients. We sought to replicate the association of two SNPs in the CHRNA7 gene, rs6494223 and rs8024987, with response to ChEI treatment in an Italian cohort of 169 AD patients, further extending the study to gene-level analysis. None of the tested variants was associated with clinical response. However, rs6494223 showed a consistent effect direction (OR = 1.4; p = 0.17), which after meta-analysis with previous study yielded a significant result (OR = 1.57, p = 0.02, I2 = 0%).