Melissa Sorosina

MicroRNA and mRNA expression profile screening in multiple sclerosis patients to unravel novel pathogenic steps and identify potential biomarkers

Identification of novel targets and biomarkers, such as microRNAs, is extremely helpful to understand the pathogenetic mechanisms in a disease like multiple sclerosis (MS). We tested the expression profile of 1145 microRNAs in peripheral blood mononuclear cells (PBMCs) of 19 MS patients and 14 controls, and we further explored their function by performing a whole-genome mRNA profiling in same subjects and using bioinformatic prediction tool. A total of 104 miRNAs have been identified as deregulated in MS patients; 2/10 which ranked highest (let-7g and miR-150) have been validated in a replication sample, leading to the identification of putative target genes.

A genome-wide association study in progressive multiple sclerosis

Background: The role played by genetic factors in influencing the clinical course of multiple sclerosis (MS) is not yet well established. Objective: We aimed to identify genetic variants associated with progressive MS (PrMS). Methods: We conducted a genome-wide association study (GWAS) in 197 patients with PrMS and 234 controls of Italian origin. We tested the top 20 single nucleotide polymorphisms (SNPs) with suggestive evidence of association (p-value<10−4) in two independent sets of primary progressive MS cases and controls. Results: We identified a risk-associated SNP in the HLA region in linkage disequilibrium (LD) with DRB1*1501 and DQB*0602 loci, with genome-wide significance (rs3129934T, p combined=6.7×10-16, OR=2.34, 95% CI=1.90–2.87), and a novel locus on chromosome 7q35 with suggestive evidence of association (rs996343G, p combined=2.4×10-5, OR=0.70, 95% CI=0.59–0.83) which maps within a human endogenous retroviral (HERV) element. The new locus did not have a ‘cis’ effect on RNA expression in lymphoblastic cell lines, but pathway analyses of ‘trans’ effects point to an expression regulation of genes involved in neurodegeneration, including glutamate metabolism (p<0.01) and axonal guidance signalling (p<0.02). Conclusions: We have confirmed the established association with the HLA region and, despite the low statistical power of the study, we found suggestive evidence for association with a novel locus on chromosome 7, with a putative regulatory role.

Genetic variants are major determinants of CSF antibody levels in multiple sclerosis

Immunological hallmarks of multiple sclerosis include the production of antibodies in the central nervous system, expressed as presence of oligoclonal bands and/or an increased immunoglobulin G index-the level of immunoglobulin G in the cerebrospinal fluid compared to serum. However, the underlying differences between oligoclonal band-positive and -negative patients with multiple sclerosis and reasons for variability in immunoglobulin G index are not known. To identify genetic factors influencing the variation in the antibody levels in the cerebrospinal fluid in multiple sclerosis, we have performed a genome-wide association screen in patients collected from nine countries for two traits, presence or absence of oligoclonal bands (n = 3026) and immunoglobulin G index levels (n = 938), followed by a replication in 3891 additional patients. We replicate previously suggested association signals for oligoclonal band status in the major histocompatibility complex region for the rs9271640*A-rs6457617*G haplotype, correlated with HLA-DRB1*1501, and rs34083746*G, correlated with HLA-DQA1*0301 (P comparing two haplotypes = 8.88 × 10(-16)). Furthermore, we identify a novel association signal of rs9807334, near the ELAC1/SMAD4 genes, for oligoclonal band status (P = 8.45 × 10(-7)). The previously reported association of the immunoglobulin heavy chain locus with immunoglobulin G index reaches strong evidence for association in this data set (P = 3.79 × 10(-37)). We identify two novel associations in the major histocompatibility complex region with immunoglobulin G index: the rs9271640*A-rs6457617*G haplotype (P = 1.59 × 10(-22)), shared with oligoclonal band status, and an additional independent effect of rs6457617*G (P = 3.68 × 10(-6)). Variants identified in this study account for up to 2-fold differences in the odds of being oligoclonal band positive and 7.75% of the variation in immunoglobulin G index. Both traits are associated with clinical features of disease such as female gender, age at onset and severity. This is the largest study population so far investigated for the genetic influence on antibody levels in the cerebrospinal fluid in multiple sclerosis, including 6950 patients. We confirm that genetic factors underlie these antibody levels and identify both the major histocompatibility complex and immunoglobulin heavy chain region as major determinants.

A pharmacogenetic study implicates SLC9a9 in multiple sclerosis disease activity

A proportion of multiple sclerosis (MS) patients experience disease activity despite treatment. The early identification of the most effective drug is critical to impact long‐term outcome and to move toward a personalized approach. The aim of the present study is to identify biomarkers for further clinical development and to yield insights into the pathophysiology of disease activity.

Multiple sclerosis risk loci and disease severity in 7,125 individuals from 10 studies

Objective: We investigated the association between 52 risk variants identified through genome-wide association studies and disease severity in multiple sclerosis (MS). Methods: Ten unique MS case data sets were analyzed. The Multiple Sclerosis Severity Score (MSSS) was calculated using the Expanded Disability Status Scale at study entry and disease duration. MSSS was considered as a continuous variable and as 2 dichotomous variables (median and extreme ends; MSSS of ≤5 vs >5 and MSSS of <2.5 vs ≥7.5, respectively). Single nucleotide polymorphisms (SNPs) were examined individually and as both combined weighted genetic risk score (wGRS) and unweighted genetic risk score (GRS) for association with disease severity. Random-effects meta-analyses were conducted and adjusted for cohort, sex, age at onset, and HLA-DRB1*15:01. Results: A total of 7,125 MS cases were analyzed. The wGRS and GRS were not strongly associated with disease severity after accounting for cohort, sex, age at onset, and HLA-DRB1*15:01. After restricting analyses to cases with disease duration ≥10 years, associations were null (p value ≥0.05). No SNP was associated with disease severity after adjusting for multiple testing. Conclusions: The largest meta-analysis of established MS genetic risk variants and disease severity, to date, was performed. Results suggest that the investigated MS genetic risk variants are not associated with MSSS, even after controlling for potential confounders. Further research in large cohorts is needed to identify genetic determinants of disease severity using sensitive clinical and MRI measures, which are critical to understanding disease mechanisms and guiding development of effective treatments.

A Strong Anti-Inflammatory Signature Revealed by Liver Transcription Profiling of Tmprss6−/− Mice

Control of systemic iron homeostasis is interconnected with the inflammatory response through the key iron regulator, the antimicrobial peptide hepcidin. We have previously shown that mice with iron deficiency anemia (IDA)-low hepcidin show a pro-inflammatory response that is blunted in iron deficient-high hepcidin Tmprss6 KO mice. The transcriptional response associated with chronic hepcidin overexpression due to genetic inactivation of Tmprss6 is unknown. By using whole genome transcription profiling of the liver and analysis of spleen immune-related genes we identified several functional pathways differentially expressed in Tmprss6 KO mice, compared to IDA animals and thus irrespective of the iron status. In the effort of defining genes potentially targets of Tmprss6 we analyzed liver gene expression changes according to the genotype and independently of treatment. Tmprss6 inactivation causes down-regulation of liver pathways connected to immune and inflammatory response as well as spleen genes related to macrophage activation and inflammatory cytokines production. The anti-inflammatory status of Tmprss6 KO animals was confirmed by the down-regulation of pathways related to immunity, stress response and intracellular signaling in both liver and spleen after LPS treatment. Opposite to Tmprss6 KO mice, Hfe−/− mice are characterized by iron overload with inappropriately low hepcidin levels. Liver expression profiling of Hfe−/− deficient versus iron loaded mice show the opposite expression of some of the genes modulated by the loss of Tmprss6. Altogether our results confirm the anti-inflammatory status of Tmprss6 KO mice and identify new potential target pathways/genes of Tmprss6.

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.