Epidemiology of seropositive myasthenia gravis in Sardinia: A population-based study in the district of Sassari.

INTRODUCTION/AIMS: The global incidence and prevalence of myasthenia gravis (MG) range between 6-31/million and 10-37/100,000, respectively. Sardinia is a high-risk region for different immune-mediated disorders, but the epidemiology of MG remains unclear. We determined the epidemiology of MG with acetylcholine receptor (AChR)-immunoglobulin G (IgG) and muscle-specific tyrosine kinase (MuSK)-IgG in the district of Sassari (North-Western Sardinia; population, 325,288). METHODS: From the laboratory of the University Hospital of Sassari (reference for AChR/MuSK-IgG testing in Sardinia since 1998) and the main neurology units in Sardinia, we retrospectively identified MG patients with (1) AChR-IgG and/or MuSK-IgG positivity by radioimmunoprecipitation assay; and (2) residency in the district of Sassari. Incidence (January 2010-December 2019) and prevalence (December 31, 2019) were calculated. RESULTS: A total of 202 patients were included (incident, 107; prevalent, 180). Antibody specificities were AChR (n = 187 [93%]) and MuSK (n = 15 [7%]). The crude MG incidence (95% confidence interval) was 32.6 (26.8-39.2)/million, while prevalence was 55.3 (47.7-63.9)/100,000. After age-standardization to the world population, incidence decreased to 18.4 (14.3-22.5)/million, while prevalence decreased to 31.6 (26.1-37.0)/100,000. Among incident cases, age strata (years) at MG onset were: <18 (2%), 18-49 (14%), 50-64 (21%), and ≥65 (63%). DISCUSSION: Sardinia is a high-risk region for MG, with a prevalence that exceeds the European threshold for rare disease. Identification of the environmental and genetic determinants of this risk may improve our understanding of disease pathophysiology.

Homozygosity Haplotype and Whole-Exome Sequencing Analysis to Identify Potentially Functional Rare Variants Involved in Multiple Sclerosis among Sardinian Families.

Multiple Sclerosis (MS) is a complex multifactorial autoimmune disease, whose sex- and age-adjusted prevalence in Sardinia (Italy) is among the highest worldwide. To date, 233 loci were associated with MS and almost 20% of risk heritability is attributable to common genetic variants, but many low-frequency and rare variants remain to be discovered. Here, we aimed to contribute to the understanding of the genetic basis of MS by investigating potentially functional rare variants. To this end, we analyzed thirteen multiplex Sardinian families with Immunochip genotyping data. For five families, Whole Exome Sequencing (WES) data were also available. Firstly, we performed a non-parametric Homozygosity Haplotype analysis for identifying the Region from Common Ancestor (RCA). Then, on these potential disease-linked RCA, we searched for the presence of rare variants shared by the affected individuals by analyzing WES data. We found: (i) a variant (43181034 T > G) in the splicing region on exon 27 of CUL9; (ii) a variant (50245517 A > C) in the splicing region on exon 16 of ATP9A; (iii) a non-synonymous variant (43223539 A > C), on exon 9 of TTBK1; (iv) a non-synonymous variant (42976917 A > C) on exon 9 of PPP2R5D; and v) a variant (109859349-109859354) in 3'UTR of MYO16.

An exome study of Parkinson's disease in Sardinia, a Mediterranean genetic isolate.

Parkinson's disease (PD) is a common neurodegenerative disorder of complex aetiology. Rare, highly penetrant PD-causing mutations and common risk factors of small effect size have been identified in several genes/loci. However, these mutations and risk factors only explain a fraction of the disease burden, suggesting that additional, substantial genetic determinants remain to be found. Genetically isolated populations offer advantages for dissecting the genetic architecture of complex disorders, such as PD. We performed exome sequencing in 100 unrelated PD patients from Sardinia, a genetic isolate. SNPs absent from dbSNP129 and 1000 Genomes, shared by at least five patients, and of functional effects were genotyped in an independent Sardinian case-control sample (n = 500). Variants associated with PD with nominal p value <0.05 and those with odds ratio (OR) ≥3 were validated by Sanger sequencing and typed in a replication sample of 2965 patients and 2678 controls from Italy, Spain, and Portugal. We identified novel moderately rare variants in several genes, including SCAPER, HYDIN, UBE2H, EZR, MMRN2 and OGFOD1 that were specifically present in PD patients or enriched among them, nominating these as novel candidate risk genes for PD, although no variants achieved genome-wide significance after Bonferroni correction. Our results suggest that the genetic bases of PD are highly heterogeneous, with implications for the design of future large-scale exome or whole-genome analyses of this disease.

Multiple Sclerosis Heritability Estimation on Sardinian Ascertained Extended Families Using Bayesian Liability Threshold Model.

Heritability studies represent an important tool to investigate the main sources of variability for complex diseases, whose etiology involves both genetics and environmental factors. In this paper, we aimed to estimate multiple sclerosis (MS) narrow-sense heritability (h2), on a liability scale, using extended families ascertained from affected probands sampled in the Sardinian province of Nuoro, Italy. We also investigated the sources of MS liability variability among shared environment effects, sex, and categorized year of birth (<1946, ≥1946). The latter can be considered a proxy for different early environmental exposures. To this aim, we implemented a Bayesian liability threshold model to obtain posterior distributions for the parameters of interest adjusting for ascertainment bias. Our analysis highlighted categorized year of birth as the main explanatory factor, explaining ~70% of MS liability variability (median value = 0.69, 95% CI: 0.64, 0.73), while h2 resulted near to 0% (median value = 0.03, 95% CI: 0.00, 0.09). By performing a year of birth-stratified analysis, we found a high h2 only in individuals born on/after 1946 (median value = 0.82, 95% CI: 0.68, 0.93), meaning that the genetic variability acquired a high explanatory role only when focusing on this subpopulation. Overall, the results obtained highlighted early environmental exposures, in the Sardinian population, as a meaningful factor involved in MS to be further investigated.

Heritability Estimation of Multiple Sclerosis Related Plasma Protein Levels in Sardinian Families with Immunochip Genotyping Data.

This work aimed at estimating narrow-sense heritability, defined as the proportion of the phenotypic variance explained by the sum of additive genetic effects, via Haseman-Elston regression for a subset of 56 plasma protein levels related to Multiple Sclerosis (MS). These were measured in 212 related individuals (with 69 MS cases and 143 healthy controls) obtained from 20 Sardinian families with MS history. Using pedigree information, we found seven statistically significant heritable plasma protein levels (after multiple testing correction), i.e., Gc (h2 = 0.77; 95%CI: 0.36, 1.00), Plat (h2 = 0.70; 95%CI: 0.27, 0.95), Anxa1 (h2 = 0.68; 95%CI: 0.27, 1.00), Sod1 (h2 = 0.58; 95%CI: 0.18, 0.96), Irf8 (h2 = 0.56; 95%CI: 0.19, 0.99), Ptger4 (h2 = 0.45; 95%CI: 0.10, 0.96), and Fadd (h2 = 0.41; 95%CI: 0.06, 0.84). A subsequent analysis was performed on these statistically significant heritable plasma protein levels employing Immunochip genotyping data obtained in 155 healthy controls (92 related and 63 unrelated); we found a meaningful proportion of heritable plasma protein levels' variability explained by a small set of SNPs. Overall, the results obtained, for these seven MS-related proteins, emphasized a high additive genetic variance component explaining plasma levels' variability.

Plasma Protein Levels Analysis in Multiple Sclerosis Sardinian Families Identified C9 and CYP24A1 as Candidate Biomarkers.

Here we investigate protein levels in 69 multiple sclerosis (MS) cases and 143 healthy controls (HC) from twenty Sardinian families to search for promising biomarkers in plasma. Using antibody suspension bead array technology, the plasma levels of 56 MS-related proteins were obtained. Differences between MS cases and HC were estimated using Linear Mixed Models or Linear Quantile Mixed Models. The proportion of proteins level variability, explained by a set of 119 MS-risk SNPs as to the literature, was also quantified. Higher plasma C9 and CYP24A1 levels were found in MS cases compared to HC (p < 0.05 after Holm multiple testing correction), with protein level differences estimated as, respectively, 0.53 (95% CI: 0.25, 0.81) and 0.42 (95% CI: 0.19, 0.65) times plasma level standard deviation measured in HC. Furthermore, C9 resulted in both statistically significantly higher relapsing-remitting MS (RRMS) and secondary-progressive MS (SPMS) compared to HC, with SPMS showing the highest differences. Instead, CYP24A1 was statistically significantly higher only in RRMS as compared to HC. Respectively, 26% (95% CI: 10%, 44%) and 16% (95% CI: 9%, 39%) of CYP24A1 and C9 plasma level variability was explained by known MS-risk SNPs. Our results highlight C9 and CYP24A1 as potential biomarkers in plasma for MS and allow us to gain insight into molecular disease mechanisms.

Broadening the phenotype of TARDBP mutations: the TARDBP Ala382Thr mutation and Parkinson's disease in Sardinia.

Mutations in the TARDBP gene are a cause of autosomal dominant amyotrophic lateral sclerosis (ALS) and of frontotemporal lobar degeneration (FTLD), but they have not been found so far in patients with Parkinson's disease (PD). A founder TARDBP mutation (p.Ala382Thr) was recently identified as the cause of ~30% of ALS cases in Sardinia, a Mediterranean genetic isolate. We studied 327 consecutive Sardinian patients with clinically diagnosed PD (88 familial, 239 sporadic) and 578 Sardinian controls. One family with FTLD and parkinsonism was also included. The p.Ala382Thr heterozygous mutation was detected in eight unrelated PD patients (2.5%). The three patients from the FTLD/parkinsonism family also carried this mutation. Within the control group, there were three heterozygous mutation carriers. During follow-up, one of these individuals developed motoneuron disease and another, a rapidly progressive dementia; the third remains healthy at the age of 79 but two close relatives developed motoneuron disease and dementia. The eight PD patients carrying the p.Ala382Thr mutation had all sporadic disease presentation. Their average onset age was 70.0 years (SD 9.4, range 51-79), which is later but not significantly different from that of the patients who did not carry this mutation. In conclusion, we expand the clinical spectrum associated with TARDBP mutations to FTLD with parkinsonism without motoneuron disease and to clinically definite PD. The TDP-43 protein might be directly involved in a broader neurodegenerative spectrum, including not only motoneuron disease and FTLD but also PD.

Investigating the Causal Effect of Brain Expression of CCL2, NFKB1, MAPK14, TNFRSF1A, CXCL10 Genes on Multiple Sclerosis: A Two-Sample Mendelian Randomization Approach.

Multiple Sclerosis (MS) exhibits considerable heterogeneity in phenotypic expression, course, prognosis and response to therapy. This suggests this disease involves multiple, as yet poorly understood, causal mechanisms. In this work we assessed the possible causal link between gene expression level of five selected genes related to the pro-inflammatory NF-κB signaling pathway (i.e., CCL2, NFKB1, MAPK14, TNFRSF1A, CXCL10) in ten different brain tissues (i.e., cerebellum, frontal cortex, hippocampus, medulla, occipital cortex, putamen, substantia nigra, thalamus, temporal cortex and intralobular white matter) and MS. We adopted a two-stage Mendelian Randomization (MR) approach for the estimation of the causal effects of interest, based on summary-level data from 20 multiplex Sardinian families and data provided by the United Kingdom Brain Expression Consortium (UKBEC). Through Radial-MR and Cochrane's Q statistics we identified and removed genetic variants which are most likely to be invalid instruments. To estimate the total causal effect, univariable MR was carried out separately for each gene and brain region. We used Inverse-Variance Weighted estimator (IVW) as main analysis and MR-Egger Regression estimator (MR-ER) and Weighted Median Estimator (WME) as sensitivity analysis. As these genes belong to the same pathway and thus they can be closely related, we also estimated their direct causal effects by applying IVW and MR-ER within a multivariable MR (MVMR) approach using set of genetic instruments specific and common (composite) to each multiple exposures represented by the expression of the candidate genes. Univariate MR analysis showed a significant positive total causal effect for CCL2 and NFKB1 respectively in medulla and cerebellum. MVMR showed a direct positive causal effect for NFKB1 and TNFRSF1A, and a direct negative causal effect for CCL2 in cerebellum; while in medulla we observed a direct positive causal effect for CCL2. Since in general we observed a different magnitude for the gene specific causal effect we hypothesize that in cerebellum and medulla the effect of each gene expression is direct but also mediated by the others. These results confirm the importance of the involvement of NF-κB signaling pathway in brain tissue for the development of the disease and improve our understanding in the pathogenesis of MS.

Association between protective and deleterious HLA alleles with multiple sclerosis in Central East Sardinia.

The human leukocyte antigen (HLA) complex on chromosome 6p21 has been unambiguously associated with multiple sclerosis (MS). The complex features of the HLA region, especially its high genic content, extreme polymorphism, and extensive linkage disequilibrium, has prevented to resolve the nature of HLA association in MS. We performed a family based association study on the isolated population of the Nuoro province (Sardinia) to clarify the role of HLA genes in MS. The main stage of our study involved an analysis of the ancestral haplotypes A2Cw7B58DR2DQ1 and A30Cw5B18DR3DQ2. On the basis of a multiplicative model, the effect of the first haplotype is protective with an odds ratio (OR) = 0.27 (95% confidence interval CI 0.13-0.57), while that of the second is deleterious, OR 1.78 (95% CI 1.26-2.50). We found both class I (A, Cw, B) and class II (DR, DQ) loci to have an effect on MS susceptibility, but we saw that they act independently from each other. We also performed an exploratory analysis on a set of 796 SNPs in the same HLA region. Our study supports the claim that Class I and Class II loci act independently on MS susceptibility and this has a biological explanation. Also, the analysis of SNPs suggests that there are other HLA genes involved in MS, but replication is needed. This opens up new perspective on the study of MS.

Association between the ACCN1 gene and multiple sclerosis in Central East Sardinia.

Multiple genome screens have been performed to identify regions in linkage or association with Multiple Sclerosis (MS, OMIM 126200), but little overlap has been found among them. This may be, in part, due to a low statistical power to detect small genetic effects and to genetic heterogeneity within and among the studied populations. Motivated by these considerations, we studied a very special population, namely that of Nuoro, Sardinia, Italy. This is an isolated, old, and genetically homogeneous population with high prevalence of MS. Our study sample includes both nuclear families and unrelated cases and controls. A multi-stage study design was adopted. In the first stage, microsatellites were typed in the 17q11.2 region, previously independently found to be in linkage with MS. One significant association was found at microsatellite D17S798. Next, a bioinformatic screening of the region surrounding this marker highlighted an interesting candidate MS susceptibility gene: the Amiloride-sensitive Cation Channel Neuronal 1 (ACCN1) gene. In the second stage of the study, we resequenced the exons and the 3' untranslated (UTR) region of ACCN1, and investigated the MS association of Single Nucleotide Polymorphisms (SNPs) identified in that region. For this purpose, we developed a method of analysis where complete, phase-solved, posterior-weighted haplotype assignments are imputed for each study individual from incomplete, multi-locus, genotyping data. The imputed assignments provide an input to a number of proposed procedures for testing association at a microsatellite level or of a sequence of SNPs. These include a Mantel-Haenszel type test based on expected frequencies of pseudocase/pseudocontrol haplotypes, as well as permutation based tests, including a combination of permutation and weighted logistic regression analysis. Application of these methods allowed us to find a significant association between MS and the SNP rs28936 located in the 3' UTR segment of ACCN1 with p = 0.0004 (p = 0.002, after adjusting for multiple testing). This result is in tune with several recent experimental findings which suggest that ACCN1 may play an important role in the pathogenesis of MS.