Genome-wide significant association with seven novel multiple sclerosis risk loci.

OBJECTIVE: A recent large-scale study in multiple sclerosis (MS) using the ImmunoChip platform reported on 11 loci that showed suggestive genetic association with MS. Additional data in sufficiently sized and independent data sets are needed to assess whether these loci represent genuine MS risk factors. METHODS: The lead SNPs of all 11 loci were genotyped in 10 796 MS cases and 10 793 controls from Germany, Spain, France, the Netherlands, Austria and Russia, that were independent from the previously reported cohorts. Association analyses were performed using logistic regression based on an additive model. Summary effect size estimates were calculated using fixed-effect meta-analysis. RESULTS: Seven of the 11 tested SNPs showed significant association with MS susceptibility in the 21 589 individuals analysed here. Meta-analysis across our and previously published MS case-control data (total sample size n=101 683) revealed novel genome-wide significant association with MS susceptibility (p<5×10(-8)) for all seven variants. This included SNPs in or near LOC100506457 (rs1534422, p=4.03×10(-12)), CD28 (rs6435203, p=1.35×10(-9)), LPP (rs4686953, p=3.35×10(-8)), ETS1 (rs3809006, p=7.74×10(-9)), DLEU1 (rs806349, p=8.14×10(-12)), LPIN3 (rs6072343, p=7.16×10(-12)) and IFNGR2 (rs9808753, p=4.40×10(-10)). Cis expression quantitative locus effects were observed in silico for rs6435203 on CD28 and for rs9808753 on several immunologically relevant genes in the IFNGR2 locus. CONCLUSIONS: This study adds seven loci to the list of genuine MS genetic risk factors and further extends the list of established loci shared across autoimmune diseases.

HERV-W polymorphism in chromosome X is associated with multiple sclerosis risk and with differential expression of MSRV.

BACKGROUND: Multiple Sclerosis (MS) is an autoimmune demyelinating disease that occurs more frequently in women than in men. Multiple Sclerosis Associated Retrovirus (MSRV) is a member of HERV-W, a multicopy human endogenous retroviral family repeatedly implicated in MS pathogenesis. MSRV envelope protein is elevated in the serum of MS patients and induces inflammation and demyelination but, in spite of this pathogenic potential, its exact genomic origin and mechanism of generation are unknown. A possible link between the HERV-W copy on chromosome Xq22.3, that contains an almost complete open reading frame, and the gender differential prevalence in MS has been suggested. RESULTS: MSRV transcription levels were higher in MS patients than in controls (U-Mann-Whitney; p = 0.004). Also, they were associated with the clinical forms (Spearman; p = 0.0003) and with the Multiple Sclerosis Severity Score (MSSS) (Spearman; p = 0.016). By mapping a 3 kb region in Xq22.3, including the HERV-W locus, we identified three polymorphisms: rs6622139 (T/C), rs6622140 (G/A) and rs1290413 (G/A). After genotyping 3127 individuals (1669 patients and 1458 controls) from two different Spanish cohorts, we found that in women rs6622139 T/C was associated with MS susceptibility: [χ2; p = 0.004; OR (95% CI) = 0.50 (0.31-0.81)] and severity, since CC women presented lower MSSS scores than CT (U-Mann-Whitney; p = 0.039) or TT patients (U-Mann-Whitney; p = 0.031). Concordantly with the susceptibility conferred in women, rs6622139*T was associated with higher MSRV expression (U-Mann-Whitney; p = 0.003). CONCLUSIONS: Our present work supports the hypothesis of a direct involvement of HERV-W/MSRV in MS pathogenesis, identifying a genetic marker on chromosome X that could be one of the causes underlying the gender differences in MS.

MicroRNAs Associated with Disability Progression and Clinical Activity in Multiple Sclerosis Patients Treated with Glatiramer Acetate.

MicroRNAs (miRNAs) are promising biomarkers in multiple sclerosis (MS). This study aims to investigate the association between a preselected list of miRNAs in serum with therapeutic response to Glatiramer Acetate (GA) and with the clinical evolution of a cohort of relapsing-remitting MS (RRMS) patients. We conducted a longitudinal study for 5 years, with cut-off points at 2 and 5 years, including 26 RRMS patients treated with GA for at least 6 months. A total of 6 miRNAs from a previous study (miR-9.5p, miR-126.3p, mir-138.5p, miR-146a.5p, miR-200c.3p, and miR-223.3p) were selected for this analysis. Clinical relapse, MRI activity, confirmed disability progression (CDP), alone or in combination (No Evidence of Disease Activity-3) (NEDA-3), and Expanded Disability Status Scale (EDSS), were studied. After multivariate regression analysis, miR-9.5p was associated with EDSS progression at 2 years (ß = 0.23; 95% CI: 0.04-0.46; p = 0.047). Besides this, mean miR-138.5p values were lower in those patients with NEDA-3 at 2 years (p = 0.033), and miR-146a.5p and miR-126.3p were higher in patients with CDP progression at 2 years (p = 0.044 and p = 0.05 respectively. These results reinforce the use of microRNAs as potential biomarkers in multiple sclerosis. We will need more studies to corroborate these data and to better understand the role of microRNAs in the pathophysiology of this disease.

Automated quantitative trait locus analysis (AutoQTL).

Background: Quantitative Trait Locus (QTL) analysis and Genome-Wide Association Studies (GWAS) have the power to identify variants that capture significant levels of phenotypic variance in complex traits. However, effort and time are required to select the best methods and optimize parameters and pre-processing steps. Although machine learning approaches have been shown to greatly assist in optimization and data processing, applying them to QTL analysis and GWAS is challenging due to the complexity of large, heterogenous datasets. Here, we describe proof-of-concept for an automated machine learning approach, AutoQTL, with the ability to automate many complex decisions related to analysis of complex traits and generate diverse solutions to describe relationships that exist in genetic data. Results: Using a dataset of 18 putative QTL from a large-scale GWAS of body mass index in the laboratory rat, Rattus norvegicus , AutoQTL captures the phenotypic variance explained under a standard additive model while also providing evidence of non-additive effects including deviations from additivity and 2-way epistatic interactions from simulated data via multiple optimal solutions. Additionally, feature importance metrics provide different insights into the inheritance models and predictive power of multiple GWAS-derived putative QTL. Conclusions: This proof-of-concept illustrates that automated machine learning techniques can be applied to genetic data and has the potential to detect both additive and non-additive effects via various optimal solutions and feature importance metrics. In the future, we aim to expand AutoQTL to accommodate omics-level datasets with intelligent feature selection strategies.

Automated quantitative trait locus analysis (AutoQTL).

BACKGROUND: Quantitative Trait Locus (QTL) analysis and Genome-Wide Association Studies (GWAS) have the power to identify variants that capture significant levels of phenotypic variance in complex traits. However, effort and time are required to select the best methods and optimize parameters and pre-processing steps. Although machine learning approaches have been shown to greatly assist in optimization and data processing, applying them to QTL analysis and GWAS is challenging due to the complexity of large, heterogenous datasets. Here, we describe proof-of-concept for an automated machine learning approach, AutoQTL, with the ability to automate many complicated decisions related to analysis of complex traits and generate solutions to describe relationships that exist in genetic data. RESULTS: Using a publicly available dataset of 18 putative QTL from a large-scale GWAS of body mass index in the laboratory rat, Rattus norvegicus, AutoQTL captures the phenotypic variance explained under a standard additive model. AutoQTL also detects evidence of non-additive effects including deviations from additivity and 2-way epistatic interactions in simulated data via multiple optimal solutions. Additionally, feature importance metrics provide different insights into the inheritance models and predictive power of multiple GWAS-derived putative QTL. CONCLUSIONS: This proof-of-concept illustrates that automated machine learning techniques can complement standard approaches and have the potential to detect both additive and non-additive effects via various optimal solutions and feature importance metrics. In the future, we aim to expand AutoQTL to accommodate omics-level datasets with intelligent feature selection and feature engineering strategies.

Genome-wide association studies of human and rat BMI converge on synapse, epigenome, and hormone signaling networks.

A vexing observation in genome-wide association studies (GWASs) is that parallel analyses in different species may not identify orthologous genes. Here, we demonstrate that cross-species translation of GWASs can be greatly improved by an analysis of co-localization within molecular networks. Using body mass index (BMI) as an example, we show that the genes associated with BMI in humans lack significant agreement with those identified in rats. However, the networks interconnecting these genes show substantial overlap, highlighting common mechanisms including synaptic signaling, epigenetic modification, and hormonal regulation. Genetic perturbations within these networks cause abnormal BMI phenotypes in mice, too, supporting their broad conservation across mammals. Other mechanisms appear species specific, including carbohydrate biosynthesis (humans) and glycerolipid metabolism (rodents). Finally, network co-localization also identifies cross-species convergence for height/body length. This study advances a general paradigm for determining whether and how phenotypes measured in model species recapitulate human biology.

MHC2TA mRNA levels and human herpesvirus 6 in multiple sclerosis patients treated with interferon beta along two-year follow-up.

BACKGROUND: In previous studies we found that MHC2TA +1614 genotype frequency was very different when MS patients with and without human herpesvirus 6 (HHV-6) in serum samples were compared; a different clinical behavior was also described. The purpose of the study was: 1. To evaluate if MHC2TA expression in MS patients was influenced by interferon beta (IFN-beta) treatment. 2. To study MHC2TA expression in MS patients with and without minor allele C. 3. To analyze the relation between MHC2TA mRNA levels and HHV-6 active infection in MS patients. METHODS: Blood and serum samples of 154 MS patients were collected in five programmed visits: basal (prior to beginning IFN-beta treatment), six, twelve, eighteen and twenty-four months later. HHV-6 in serum and MHC2TA mRNA levels were evaluated by PCR and RT-PCR, respectively. Neutralizing antibodies (NAbs) against IFN-beta were analyzed by the cytopathic effect assay. RESULTS: We found that MHC2TA mRNA levels were significantly lower among MS patients with HHV-6 active infection at the basal visit (without treatment) than in those MS patients without HHV-6 active infection at the basal visit (p = 0.012); in all the positive samples we only found variant A. Furthermore, 58/99 (58.6%) MS patients without HHV-6 along the five programmed visits and an increase of MHC2TA expression after two-years of IFN-beta treatment were clinical responders vs. 5/21 (23.8%) among those MS patients with HHV-6 and a decrease of MHC2TA mRNA levels along the two-years with IFN-beta treatment (p = 0.004); no differences were found between patients with and without NAbs. CONCLUSIONS: MHC2TA mRNA levels could be decreased by the active replication of HHV-6; the absence of HHV-6 in serum and the increase of MHC2TA expression could be further studied as markers of good clinical response to IFN-beta treatment.

Epstein-Barr Virus and multiple sclerosis in a Spanish cohort: A two-years longitudinal study.

Objectives: 1. To analyze the prevalence and levels of anti-EBNA-1 and anti-VCA IgG antibodies of Epstein-Barr virus (EBV) in a Spanish cohort of multiple sclerosis (MS) patients and their interactions with other environmental and genetic risk factors. 2. To analyze the association of the evolution of these antibodies with the clinical response to different disease modifying therapies (DMTs) after two-years of follow-up. 3. To assess their possible correlation with the class II HLA alleles as well as with several SNPs identified in GWAS related to disease susceptibility. Materials and methods: We recruited 325 MS patients without DMT (serum samples were collected 1-3 months before starting a therapy) and 295 healthy controls (HC). For each patient we also collected serum samples 6, 12, 18 and 24 months after starting the DMT. EBNA-1 and VCA IgG titers were analyzed by ELISA; 25(OH)D levels were analyzed by immunoassay; HLA DRB1*15:01 allelic variant was analyzed by Taqman technology. Results: 1. 97.8% (318/325) vs. 87.1% (257/295) positives for EBNA-1 in MS patients and HC, respectively (p<0.0001; O.R. = 6.7); 99.7% (324/325) vs. 94.6% (279/295) for VCA in MS patients and HC, respectively (p=0.0001; O.R. = 18.6). All MS patients were positive for EBNA-1 and/or VCA IgG antibodies vs. 280/295 (94.9%) HC (p<0.0001). IgG titers were also significantly higher in MS patients than in HC. 2. We did not find any statistical correlation in the variation of the EBNA-1 and VCA IgG titers between baseline and 24 month visits with the number of relapses, progression, clinical response, NEDA-3 condition or therapeutic failure. 3. When we compared different epidemiological and clinical variables between those with genetic factors associated with lower EBNA-1 IgG titers and all other MS patients, we found MS started 3.5 years later among the first. Conclusions: These results confirm that MS occurs rarely in absence of EBV. An intriguing association between genetic burden and lower EBNA-1 IgG titers was associated with an earlier age of disease onset. Similar studies with B-cell-targeted therapies should be performed.

microRNA Expression and Its Association With Disability and Brain Atrophy in Multiple Sclerosis Patients Treated With Glatiramer Acetate.

Background: MicroRNAs are small non-coding RNA that regulate gene expression at a post-transcriptional level affecting several cellular processes including inflammation, neurodegeneration and remyelination. Different patterns of miRNAs expression have been demonstrated in multiple sclerosis compared to controls, as well as in different courses of the disease. For these reason they have been postulated as promising biomarkers candidates in multiple sclerosis. Objective: to correlate serum microRNAs profile expression with disability, cognitive functioning and brain volume in patients with remitting-relapsing multiple sclerosis. Methods: cross-sectional study in relapsing-remitting multiple sclerosis patients treated with glatiramer acetate. Disability was measured with Expanded Disability Status Scale (EDSS) and cognitive function was studied with Symbol Digit Modalities Test (SDMT). Brain volume was analyzed with automatic software NeuroQuant®. Results: We found an association between miR.146a.5p (rs:0.434, p=0.03) and miR.9.5p (rs:0.516, p=0.028) with EDSS; and miR-146a.5p (rs:-0.476, p=0.016) and miR-126.3p (rs:-0.528, p=0.007) with SDMT. Regarding to the brain volume, miR.9.5p correlated with thalamus (rs:-0.545, p=0.036); miR.200c.3p with pallidum (rs:-0.68, p=0.002) and cerebellum (rs:-0.472, p=0.048); miR-138.5p with amygdala (rs:0.73, p=0.016) and pallidum (rs:0.64, p=0.048); and miR-223.3p with caudate (rs:0.46, p=0.04). Conclusions: These data support the hypothesis of microRNA as potential biomarkers in this disease. More studies are needed to validate these results and to better understand the role of microRNAs in the pathogenesis, monitoring and therapeutic response of multiple sclerosis.

Genome-Wide Association Study on Three Behaviors Tested in an Open Field in Heterogeneous Stock Rats Identifies Multiple Loci Implicated in Psychiatric Disorders.

Many personality traits are influenced by genetic factors. Rodents models provide an efficient system for analyzing genetic contribution to these traits. Using 1,246 adolescent heterogeneous stock (HS) male and female rats, we conducted a genome-wide association study (GWAS) of behaviors measured in an open field, including locomotion, novel object interaction, and social interaction. We identified 30 genome-wide significant quantitative trait loci (QTL). Using multiple criteria, including the presence of high impact genomic variants and co-localization of cis-eQTL, we identified 17 candidate genes (Adarb2, Ankrd26, Cacna1c, Cacng4, Clock, Ctu2, Cyp26b1, Dnah9, Gda, Grxcr1, Eva1a, Fam114a1, Kcnj9, Mlf2, Rab27b, Sec11a, and Ube2h) for these traits. Many of these genes have been implicated by human GWAS of various psychiatric or drug abuse related traits. In addition, there are other candidate genes that likely represent novel findings that can be the catalyst for future molecular and genetic insights into human psychiatric diseases. Together, these findings provide strong support for the use of the HS population to study psychiatric disorders.

Genome-wide association study finds multiple loci associated with intraocular pressure in HS rats.

Elevated intraocular pressure (IOP) is influenced by environmental and genetic factors. Increased IOP is a major risk factor for most types of glaucoma, including primary open angle glaucoma (POAG). Investigating the genetic basis of IOP may lead to a better understanding of the molecular mechanisms of POAG. The goal of this study was to identify genetic loci involved in regulating IOP using outbred heterogeneous stock (HS) rats. HS rats are a multigenerational outbred population derived from eight inbred strains that have been fully sequenced. This population is ideal for a genome-wide association study (GWAS) owing to the accumulated recombinations among well-defined haplotypes, the relatively high allele frequencies, the accessibility to a large collection of tissue samples, and the large allelic effect size compared to human studies. Both male and female HS rats (N = 1,812) were used in the study. Genotyping-by-sequencing was used to obtain ∼3.5 million single nucleotide polymorphisms (SNP) from each individual. SNP heritability for IOP in HS rats was 0.32, which agrees with other studies. We performed a GWAS for the IOP phenotype using a linear mixed model and used permutation to determine a genome-wide significance threshold. We identified three genome-wide significant loci for IOP on chromosomes 1, 5, and 16. Next, we sequenced the mRNA of 51 whole eye samples to find cis-eQTLs to aid in identification of candidate genes. We report 5 candidate genes within those loci: Tyr, Ctsc, Plekhf2, Ndufaf6 and Angpt2. Tyr, Ndufaf6 and Angpt2 genes have been previously implicated by human GWAS of IOP-related conditions. Ctsc and Plekhf2 genes represent novel findings that may provide new insight into the molecular basis of IOP. This study highlights the efficacy of HS rats for investigating the genetics of elevated IOP and identifying potential candidate genes for future functional testing.

CD46 in a Spanish cohort of multiple sclerosis patients: genetics, mRNA expression and response to interferon-beta treatment.

BACKGROUND: In a prior study of our group we found an up-regulation of CD46 expression in a cohort of Spanish multiple sclerosis (MS) patients. OBJECTIVE: To evaluate the potential role of CD46 in the response to interferon-beta treatment in MS patients through the analysis of five tagging single nucleotide polymorphisms (SNPs) and measurement of mRNA. METHODS: A total of 406 MS patients and 513 control patients were analysed for five SNPs at the CD46 locus. Furthermore, 163 MS patients and 163 matched control patients were analysed by RT-PCR for the CD46 mRNA expression in three blood samples (basal, and at 6 and 12 months of interferon-beta treatment) collected in the course of a 1-year follow-up. RESULTS: Two genotypes of rs2724385 polymorphism (AT and TT) could be markers of response to interferon-beta therapy in MS patients (p=0.007 and p=0.006, respectively). Furthermore, the frequency of interferon-beta responders was 44.4% (32/72) in MS patients with an increased CD46 mRNA expression, vs. 65.9% (60/91) in patients with a decreased CD46 mRNA expression (p=0.006). CONCLUSION: The present study shows that CD46 could be associated with the response to interferon-beta therapy; however, the genetic results should be replicated in an independent cohort and further studies are needed to confirm the role of CD46.

TailTimer: A device for automating data collection in the rodent tail immersion assay.

The tail immersion assay is a widely used method for measuring acute thermal pain in a way which is quantifiable and reproducible. It is non-invasive and measures response to a stimulus that may be encountered by an animal in its natural environment. However, quantification of tail withdrawal latency relies on manual timing of tail flick using a stopwatch, and precise temperatures of the water at the time of measurement are most often not recorded. These two factors greatly reduce the reproducibility of tail immersion assay data and likely contribute to some of the discrepancies present among relevant literature. We designed a device, TailTimer, which uses a Raspberry Pi single-board computer, a digital temperature sensor, and two electrical wires, to automatically record tail withdrawal latency and water temperature. We programmed TailTimer to continuously display and record water temperature and to only permit the assay to be conducted when the water is within ± 0.25°C of the target temperature. Our software also records the identification of the animals using a radio frequency identification (RFID) system. We further adapted the RFID system to recognize several specific keys as user interface commands, allowing TailTimer to be operated via RFID fobs for increased usability. Data recorded using the TailTimer device showed a negative linear relationship between tail withdrawal latency and water temperature when tested between 47-50°C. We also observed a previously unreported, yet profound, effect of water mixing speed on latency. In one experiment using TailTimer, we observed significantly longer latencies following administration of oral oxycodone versus a distilled water control when measured after 15 mins or 1 h, but not after 4 h. TailTimer also detected significant strain differences in baseline latency. These findings valorize TailTimer in its sensitivity and reliability for measuring thermal pain thresholds.

Predictive factors and early biomarkers of response in multiple sclerosis patients treated with natalizumab.

There are an increasing number of treatments available for multiple sclerosis (MS). The early identification of optimal responders to individual treatments is important to achieve individualized therapy. With this aim, we performed a multicenter retrospective longitudinal study including 186 MS patients treated with natalizumab who were followed for 2 years. We analyzed the following variables at recruitment: sex, current age, age at disease onset, disease duration, EDSS, number of T2 and Gd + lesions, IgG and IgM oligoclonal bands, HLA class II (DR, DRB, DQA, DQB, and DRB1*15:01), IgG and IgM antibody titers against human herpesvirus 6 (HHV-6) and the antibody response to Epstein-Barr virus (EBV) through the measurement of the anti-EBNA-1 and anti-VCA IgG titers, in relation to clinical response (no relapses or disability progression), and to NEDA-3 (no evidence of disease activity in terms of clinical response and no changes in MRI scans either) after 2-years follow-up. Baseline EDSS score, baseline EBNA-1 IgG titers and percentage change of HHV6 IgG titers between baseline and 6 month visits were significantly different in clinical responders and in NEDA-3 status (all of them remained significant in the multivariate analysis). We identified three variables for the early identification of natalizumab optimal responders in a rapid and cost-effective approach.

Genome-Wide Association Study in 3,173 Outbred Rats Identifies Multiple Loci for Body Weight, Adiposity, and Fasting Glucose.

OBJECTIVE: Obesity is influenced by genetic and environmental factors. Despite the success of human genome-wide association studies, the specific genes that confer obesity remain largely unknown. The objective of this study was to use outbred rats to identify the genetic loci underlying obesity and related morphometric and metabolic traits. METHODS: This study measured obesity-relevant traits, including body weight, body length, BMI, fasting glucose, and retroperitoneal, epididymal, and parametrial fat pad weight in 3,173 male and female adult N/NIH heterogeneous stock (HS) rats across three institutions, providing data for the largest rat genome-wide association study to date. Genetic loci were identified using a linear mixed model to account for the complex family relationships of the HS and using covariates to account for differences among the three phenotyping centers. RESULTS: This study identified 32 independent loci, several of which contained only a single gene (e.g., Epha5, Nrg1, Klhl14) or obvious candidate genes (e.g., Adcy3, Prlhr). There were strong phenotypic and genetic correlations among obesity-related traits, and there was extensive pleiotropy at individual loci. CONCLUSIONS: This study demonstrates the utility of HS rats for investigating the genetics of obesity-related traits across institutions and identify several candidate genes for future functional testing.

Fingolimod Use for the Treatment of Multiple Sclerosis in a Clinical Practice Setting in Madrid.

OBJECTIVE: To assess the effectiveness and safety of fingolimod use in a Spanish clinical practice setting. METHODS: Retrospective study with multiple sclerosis patients who received at least 1 fingolimod dose between January 2004 and January 2015. Effectiveness and safety data were collected during the entire treatment of each patient. Analysis was performed for the total population and stratified according to prior treatment, sex, and age at treatment initiation. RESULTS: A total of 167 patients were included, 50.9% had prior immunomodulator use, 33.5% natalizumab use, and 15.6% were naive patients. The annual relapse rate (ARR) decreased for the total population at month 12 (62%) and month 24 (84%) (P < 0.0001, in both cases); for naive patients (P < 0.05) and patients with prior immunomodulator use (P < 0.0001); for patients with prior natalizumab use, the ARR kept low after treatment initiation (0.23). After 24 months, the proportion of relapse-free patients was 70% or greater and disability progression-free patients was 80% or greater. No significant differences were observed when the results were compared by prior treatment, sex, or age. Thirty-two patients (19.2%) reported adverse drug reactions and 9.6% discontinued: 4.8% due to adverse drug reactions and 4.8% for lack of effectiveness. CONCLUSIONS: The results support fingolimod use due to clinical effectiveness, tolerability, and ease of administration.

Analysis of Plasminogen Genetic Variants in Multiple Sclerosis Patients.

Multiple sclerosis (MS) is a prevalent neurological disease of complex etiology. Here, we describe the characterization of a multi-incident MS family that nominated a rare missense variant (p.G420D) in plasminogen (PLG) as a putative genetic risk factor for MS. Genotyping of PLG p.G420D (rs139071351) in 2160 MS patients, and 886 controls from Canada, identified 10 additional probands, two sporadic patients and one control with the variant. Segregation in families harboring the rs139071351 variant, identified p.G420D in 26 out of 30 family members diagnosed with MS, 14 unaffected parents, and 12 out of 30 family members not diagnosed with disease. Despite considerably reduced penetrance, linkage analysis supports cosegregation of PLG p.G420D and disease. Genotyping of PLG p.G420D in 14446 patients, and 8797 controls from Canada, France, Spain, Germany, Belgium, and Austria failed to identify significant association with disease (P = 0.117), despite an overall higher prevalence in patients (OR = 1.32; 95% CI = 0.93-1.87). To assess whether additional rare variants have an effect on MS risk, we sequenced PLG in 293 probands, and genotyped all rare variants in cases and controls. This analysis identified nine rare missense variants, and although three of them were exclusively observed in MS patients, segregation does not support pathogenicity. PLG is a plausible biological candidate for MS owing to its involvement in immune system response, blood-brain barrier permeability, and myelin degradation. Moreover, components of its activation cascade have been shown to present increased activity or expression in MS patients compared to controls; further studies are needed to clarify whether PLG is involved in MS susceptibility.

Anti-human herpesvirus 6A/B IgG correlates with relapses and progression in multiple sclerosis.

OBJECTIVE: To analyze the titers of the IgG and IgM antibodies against human herpesvirus 6A/B (HHV-6A/B) in multiple sclerosis (MS) patients treated with different disease modified therapies (DMTs) along two-years of follow-up. METHODS: We collected 2163 serum samples from 596 MS; for 301 MS patients a 2-years follow-up was performed. Serum samples of 337 healthy controls were also analyzed. Anti-HHV-6A/B IgG and IgM were analyzed by ELISA (Panbio). RESULTS: We found that 129/187 (69.0%) MS patients with a decrease of the anti-HHV-6A/B IgG titers after 2-years with DMTs were free of relapses and progression vs. 46/113 (40.7%) of MS patients with an increase of the anti-HHV-6A/B IgG titers (p = 0.0000015); the higher significance was found for natalizumab. Furthermore, we found that anti-HHV-6A/B IgG titers reached their highest value two weeks before the relapse (p = 0.0142), while the anti-HHV-6A/B IgM titers reached their highest value one month before the relapse (p = 0.0344). CONCLUSION: The measurement of the anti-HHV-6A/B IgG titers could be a good biomarker of clinical response to the different DMTs. The increase of the anti-HHV-6A/B IgG and IgM titers predicts the upcoming clinical relapses. However, further longitudinal studies are needed to validate these results.

Anti-JCV antibodies detection and JCV DNA levels in PBMC, serum and urine in a cohort of Spanish Multiple Sclerosis patients treated with natalizumab.

One of the most effective multiple sclerosis (MS) treatment is natalizumab. Nevertheless, it has been associated with an increased risk of progressive multifocal leukoencephalopathy (PML) caused by the JC virus (JCV). Our main objective was to assess the utility of testing JCV-DNA, apart from anti-JCV antibodies, to determine which natalizumab-treated MS patients has been previously in contact with the virus. For this purpose, 138 MS natalizumab/non-natalizumab treated patients participated in several studies. Cross-sectional study (CS): association of several epidemiological variables with anti-JCV antibodies and JCV-DNA levels in PBMC/serum/urine. First longitudinal study (A): evaluation of JCV-DNA prevalence in urine throughout the treatment. Second longitudinal study (B): simultaneous assessment of antibodies and viral DNA levels in PBMC/serum/urine at two time points. CS: The seropositivity rate for anti-JCV antibodies (62.3 %) and JCV prevalence in urine (59.4 %) were similar; although 26 % of our population was positive only using one of the two techniques. A: The viral prevalence in urine seemed to increase between the baseline visit and the others (Baseline-Visit/V18months, p = 0.006). B: Our rate of positive antibody seroconversion was 36 %. Nearly all patients with detectable JCV-DNA levels in PBMC excreted the virus intermittently in urine; while our PML case, positive in PBMC and serum samples 2 month before the PML, excreted JCV permanently. In conclusion, the determination of JCV DNA levels in urine could be complementary to anti-JCV antibodies for identifying MS patients who has been infected by the JCV. Further research would be necessary to understand the different JCV excretion profiles in urine.

The DNA copy number of human endogenous retrovirus-W (MSRV-type) is increased in multiple sclerosis patients and is influenced by gender and disease severity.

BACKGROUND: Multiple sclerosis is an autoimmune disease more prevalent in women than in men. Multiple Sclerosis Associated Retrovirus element (MSRV) is a member of type-W endogenous retrovirus family (HERV-W), known to be associated to MS. Most HERVs are unable to replicate but MSRV expression associated with reverse-transcriptase activity in MS would explain reported DNA copy number increase in MS patients. A potential link between HERV-W copies on chromosome X and gender differential prevalence has been suggested. The present study addresses MSRV-type DNA load in relation with the gender differences and clinical status in MS and healthy controls. RESULTS: 178 MS patients (62.9% women) and 124 controls (56.5% women) were included. MSRV env load (copies/pg of DNA) was analyzed by real time qPCR with specific primers and probe for its env gene, in DNA from peripheral blood mononuclear cells (PBMCs). MSRV load was more elevated in MS patients than in controls (p = 4.15e-7). MS women presented higher MSRV load than control women (p = 0.009) and MS men also had higher load than control men (p = 2.77e-6). Besides, women had higher levels than men, both among patients (p = 0.007) and controls (p = 1.24e-6). Concordantly, EDSS and MSSS scores were higher among female patients with an elevated MSRV load (p = 0.03 and p = 0.04, respectively). CONCLUSIONS: MSRV increases its copy number in PBMC of MS patients and particularly in women with high clinical scores. This may explain causes underlying the higher prevalence of MS in women. The association with the clinical severity calls for further investigations on MSRV load in PBMCs as a biomarker for MS.