MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge.

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis.

Inflammatory and neurodegenerative serum protein biomarkers increase sensitivity to detect clinical and radiographic disease activity in multiple sclerosis.

The multifaceted nature of multiple sclerosis requires quantitative biomarkers that can provide insights related to diverse physiological pathways. To this end, proteomic analysis of deeply-phenotyped serum samples, biological pathway modeling, and network analysis were performed to elucidate inflammatory and neurodegenerative processes, identifying sensitive biomarkers of multiple sclerosis disease activity. Here, we evaluated the concentrations of > 1400 serum proteins in 630 samples from three multiple sclerosis cohorts for association with clinical and radiographic new disease activity. Twenty proteins were associated with increased clinical and radiographic multiple sclerosis disease activity for inclusion in a custom assay panel. Serum neurofilament light chain showed the strongest univariate correlation with gadolinium lesion activity, clinical relapse status, and annualized relapse rate. Multivariate modeling outperformed univariate for all endpoints. A comprehensive biomarker panel including the twenty proteins identified in this study could serve to characterize disease activity for a patient with multiple sclerosis.

High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis.

Background: The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods. Methods: We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls. Results: We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk. Results: The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

Inflammatory and neurodegenerative serum protein biomarkers increase sensitivity to detect disease activity in multiple sclerosis.

Background/Objectives: Serum proteomic analysis of deeply-phenotyped samples, biological pathway modeling and network analysis were performed to elucidate the inflammatory and neurodegenerative processes of multiple sclerosis (MS) and identify sensitive biomarkers of MS disease activity (DA). Methods: Over 1100 serum proteins were evaluated in >600 samples from three MS cohorts to identify biomarkers of clinical and radiographic (gadolinium-enhancing lesions) new MS DA. Protein levels were analyzed and associated with presence of gadolinium-enhancing lesions, clinical relapse status (CRS), and annualized relapse rate (ARR) to create a custom assay panel. Results: Twenty proteins were associated with increased clinical and radiographic MS DA. Serum neurofilament light chain (NfL) showed the strongest univariate correlation with radiographic and clinical DA measures. Multivariate modeling significantly outperformed univariate NfL to predict gadolinium lesion activity, CRS and ARR. Discussion: These findings provide insight regarding correlations between inflammatory and neurodegenerative biomarkers and clinical and radiographic MS DA. Funding: Octave Bioscience, Inc (Menlo Park, CA).

Specific hypomethylation programs underpin B cell activation in early multiple sclerosis.

Epigenetic changes have been consistently detected in different cell types in multiple sclerosis (MS). However, their contribution to MS pathogenesis remains poorly understood partly because of sample heterogeneity and limited coverage of array-based methods. To fill this gap, we conducted a comprehensive analysis of genome-wide DNA methylation patterns in four peripheral immune cell populations isolated from 29 MS patients at clinical disease onset and 24 healthy controls. We show that B cells from new-onset untreated MS cases display more significant methylation changes than other disease-implicated immune cell types, consisting of a global DNA hypomethylation signature. Importantly, 4,933 MS-associated differentially methylated regions in B cells were identified, and this epigenetic signature underlies specific genetic programs involved in B cell differentiation and activation. Integration of the methylome to changes in gene expression and susceptibility-associated regions further indicates that hypomethylated regions are significantly associated with the up-regulation of cell activation transcriptional programs. Altogether, these findings implicate aberrant B cell function in MS etiology.

High Resolution Haplotype Analyses of Classical HLA Genes in Families With Multiple Sclerosis Highlights the Role of HLA-DP Alleles in Disease Susceptibility.

Multiple sclerosis (MS) susceptibility shows strong genetic associations with HLA alleles and haplotypes. We genotyped 11 HLA genes in 477 non-Hispanic European MS patients and their 954 unaffected parents using a validated next-generation sequencing (NGS) methodology. HLA haplotypes were assigned unequivocally by tracing HLA allele transmissions. We explored HLA haplotype/allele associations with MS using the genotypic transmission disequilibrium test (gTDT) and multiallelic TDT (mTDT). We also conducted a case-control (CC) study with all patients and 2029 healthy unrelated ethnically matched controls. We performed separate analyses of 54 extended multi-case families by reviewing transmission of haplotype blocks. The haplotype fragment including DRB5*01:01:01~DRB1*15:01:01:01 was significantly associated with predisposition (gTDT: p < 2.20e-16; mTDT: p =1.61e-07; CC: p < 2.22e-16) as reported previously. A second risk allele, DPB1*104:01 (gTDT: p = 3.69e-03; mTDT: p = 2.99e-03; CC: p = 1.00e-02), independent from the haplotype bearing DRB1*15:01 was newly identified. The allele DRB1*01:01:01 showed significant protection (gTDT: p = 8.68e-06; mTDT: p = 4.50e-03; CC: p = 1.96e-06). Two DQB1 alleles, DQB1*03:01 (gTDT: p = 2.86e-03; mTDT: p = 5.56e-02; CC: p = 4.08e-05) and DQB1*03:03 (gTDT: p = 1.17e-02; mTDT: p = 1.16e-02; CC: p = 1.21e-02), defined at two-field level also showed protective effects. The HLA class I block, A*02:01:01:01~C*03:04:01:01~B*40:01:02 (gTDT: p = 5.86e-03; mTDT: p = 3.65e-02; CC: p = 9.69e-03) and the alleles B*27:05 (gTDT: p = 6.28e-04; mTDT: p = 2.15e-03; CC: p = 1.47e-02) and B*38:01 (gTDT: p = 3.20e-03; mTDT: p = 6.14e-03; CC: p = 1.70e-02) showed moderately protective effects independently from each other and from the class II associated factors. By comparing statistical significance of 11 HLA loci and 19 haplotype segments with both untruncated and two-field allele names, we precisely mapped MS candidate alleles/haplotypes while eliminating false signals resulting from 'hitchhiking' alleles. We assessed genetic burden for the HLA allele/haplotype identified in this study. This family-based study including the highest-resolution of HLA alleles proved to be powerful and efficient for precise identification of HLA genotypes associated with both, susceptibility and protection to development of MS.

Next-Generation Sequencing Identifies Extended HLA Class I and II Haplotypes Associated With Early-Onset and Late-Onset Myasthenia Gravis in Italian, Norwegian, and Swedish Populations.

Genetic susceptibility to myasthenia gravis (MG) associates with specific HLA alleles and haplotypes at the class I and II regions in various populations. Previous studies have only examined alleles at a limited number of HLA loci that defined only broad serotypes or alleles defined at the protein sequence level. Consequently, genetic variants in noncoding and untranslated HLA gene segments have not been fully explored but could also be important determinants for MG. To gain further insight into the role of HLA in MG, we applied next-generation sequencing to analyze sequence variation at eleven HLA genes in early-onset (EO) and late-onset (LO) non-thymomatous MG patients positive for the acetylcholine receptor (AChR) antibodies and ethnically matched controls from Italy, Norway, and Sweden. For all three populations, alleles and haplotype blocks present on the ancestral haplotype AH8.1 were associated with risk in AChR-EOMG patients. HLA-B*08:01:01:01 was the dominant risk allele in Italians (OR = 3.28, P = 1.83E-05), Norwegians (OR = 3.52, P = 4.41E-16), and in Swedes HLA-B*08:01 was the primary risk allele (OR = 4.24, P <2.2E-16). Protective alleles and haplotype blocks were identified on the HLA-DRB7, and HLA-DRB13.1 class II haplotypes in Italians and Norwegians, whereas in Swedes HLA-DRB7 exhibited the main protective effect. For AChR-LOMG patients, the HLA-DRB15.1 haplotype and associated alleles were significantly associated with susceptibility in all groups. The HLA-DR13-HLA-DR-HLA-DQ haplotype was associated with protection in all AChR-LOMG groups. This study has confirmed and extended previous findings that the immunogenetic predisposition profiles for EOMG and LOMG are distinct. In addition, the results are consistent with a role for non-coding HLA genetic variants in the pathogenesis of MG.

High-Resolution Characterization of KIR Genes in a Large North American Cohort Reveals Novel Details of Structural and Sequence Diversity.

The KIR (killer-cell immunoglobulin-like receptor) region is characterized by structural variation and high sequence similarity among genes, imposing technical difficulties for analysis. We undertook the most comprehensive study to date of KIR genetic diversity in a large population sample, applying next-generation sequencing in 2,130 United States European-descendant individuals. Data were analyzed using our custom bioinformatics pipeline specifically designed to address technical obstacles in determining KIR genotypes. Precise gene copy number determination allowed us to identify a set of uncommon gene-content KIR haplotypes accounting for 5.2% of structural variation. In this cohort, KIR2DL4 is the framework gene that most varies in copy number (6.5% of all individuals). We identified phased high-resolution alleles in large multi-locus insertions and also likely founder haplotypes from which they were deleted. Additionally, we observed 250 alleles at 5-digit resolution, of which 90 have frequencies ≥1%. We found sequence patterns that were consistent with the presence of novel alleles in 398 (18.7%) individuals and contextualized multiple orphan dbSNPs within the KIR complex. We also identified a novel KIR2DL1 variant, Pro151Arg, and demonstrated by molecular dynamics that this substitution is predicted to affect interaction with HLA-C. No previous studies have fully explored the full range of structural and sequence variation of KIR as we present here. We demonstrate that pairing high-throughput sequencing with state-of-art computational tools in a large cohort permits exploration of all aspects of KIR variation including determination of population-level haplotype diversity, improving understanding of the KIR system, and providing an important reference for future studies.

Cell type-specific transcriptomics identifies neddylation as a novel therapeutic target in multiple sclerosis.

Multiple sclerosis is an autoimmune disease of the CNS in which both genetic and environmental factors are involved. Genome-wide association studies revealed more than 200 risk loci, most of which harbour genes primarily expressed in immune cells. However, whether genetic differences are translated into cell-specific gene expression profiles and to what extent these are altered in patients with multiple sclerosis are still open questions in the field. To assess cell type-specific gene expression in a large cohort of patients with multiple sclerosis, we sequenced the whole transcriptome of fluorescence-activated cell sorted T cells (CD4+ and CD8+) and CD14+ monocytes from treatment-naive patients with multiple sclerosis (n = 106) and healthy subjects (n = 22). We identified 479 differentially expressed genes in CD4+ T cells, 435 in monocytes, and 54 in CD8+ T cells. Importantly, in CD4+ T cells, we discovered upregulated transcripts from the NAE1 gene, a critical subunit of the NEDD8 activating enzyme, which activates the neddylation pathway, a post-translational modification analogous to ubiquitination. Finally, we demonstrated that inhibition of NEDD8 activating enzyme using the specific inhibitor pevonedistat (MLN4924) significantly ameliorated disease severity in murine experimental autoimmune encephalomyelitis. Our findings provide novel insights into multiple sclerosis-associated gene regulation unravelling neddylation as a crucial pathway in multiple sclerosis pathogenesis with implications for the development of tailored disease-modifying agents.

A splice acceptor variant in HLA-DRA affects the conformation and cellular localization of the class II DR alpha-chain.

Class II human leucocyte antigen (HLA) proteins are involved in the immune response by presenting pathogen-derived peptides to CD4+ T lymphocytes. At the molecular level, they are constituted by α/ß-heterodimers on the surface of professional antigen-presenting cells. Here, we report that the acceptor variant (rs8084) in the HLA-DRA gene mediates the transcription of an alternative version of the α-chain lacking 25 amino acids in its extracellular domain. Molecular dynamics simulations suggest this isoform undergoes structural refolding which in turn affects its stability and cellular trafficking. The short HLA-DRA isoform cannot reach the cell surface, although it is still able to bind the corresponding ß-chain. Conversely, it remains entrapped within the endoplasmic reticulum where it is targeted for degradation. Furthermore, we demonstrate that the short isoform can be transported to the cell membrane via interactions with the peptide-binding site of canonical HLA heterodimers. Altogether, our findings indicate that short HLA-DRA functions as a novel intact antigen for class II HLA molecules.

Ataxin-1 regulates B cell function and the severity of autoimmune experimental encephalomyelitis.

Ataxin-1 (ATXN1) is a ubiquitous polyglutamine protein expressed primarily in the nucleus where it binds chromatin and functions as a transcriptional repressor. Mutant forms of ataxin-1 containing expanded glutamine stretches cause the movement disorder spinocerebellar ataxia type 1 (SCA1) through a toxic gain-of-function mechanism in the cerebellum. Conversely, ATXN1 loss-of-function is implicated in cancer development and Alzheimer's disease (AD) pathogenesis. ATXN1 was recently nominated as a susceptibility locus for multiple sclerosis (MS). Here, we show that Atxn1-null mice develop a more severe experimental autoimmune encephalomyelitis (EAE) course compared to wildtype mice. The aggravated phenotype is mediated by increased T helper type 1 (Th1) cell polarization, which in turn results from the dysregulation of B cell activity. Ataxin-1 ablation in B cells leads to aberrant expression of key costimulatory molecules involved in proinflammatory T cell differentiation, including cluster of differentiation (CD)44 and CD80. In addition, comprehensive phosphoflow cytometry and transcriptional profiling link the exaggerated proliferation of ataxin-1 deficient B cells to the activation of extracellular signal-regulated kinase (ERK) and signal transducer and activator of transcription (STAT) pathways. Lastly, selective deletion of the physiological binding partner capicua (CIC) demonstrates the importance of ATXN1 native interactions for correct B cell functioning. Altogether, we report a immunomodulatory role for ataxin-1 and provide a functional description of the ATXN1 locus genetic association with MS risk.

Killer Cell Immunoglobulin-like Receptor Variants Are Associated with Protection from Symptoms Associated with More Severe Course in Parkinson Disease.

Immune dysfunction plays a role in the development of Parkinson disease (PD). NK cells regulate immune functions and are modulated by killer cell immunoglobulin-like receptors (KIR). KIR are expressed on the surface of NK cells and interact with HLA class I ligands on the surface of all nucleated cells. We investigated KIR-allelic polymorphism to interrogate the role of NK cells in PD. We sequenced KIR genes from 1314 PD patients and 1978 controls using next-generation methods and identified KIR genotypes using custom bioinformatics. We examined associations of KIR with PD susceptibility and disease features, including age at disease onset and clinical symptoms. We identified two KIR3DL1 alleles encoding highly expressed inhibitory receptors associated with protection from PD clinical features in the presence of their cognate ligand: KIR3DL1*015/HLA-Bw4 from rigidity (p c = 0.02, odds ratio [OR] = 0.39, 95% confidence interval [CI] 0.23-0.69) and KIR3DL1*002/HLA-Bw4i from gait difficulties (p c = 0.05, OR = 0.62, 95% CI 0.44-0.88), as well as composite symptoms associated with more severe disease. We also developed a KIR3DL1/HLA interaction strength metric and found that weak KIR3DL1/HLA interactions were associated with rigidity (pc = 0.05, OR = 9.73, 95% CI 2.13-172.5). Highly expressed KIR3DL1 variants protect against more debilitating symptoms of PD, strongly implying a role of NK cells in PD progression and manifestation.

Telomere Length Is Associated with Disability Progression in Multiple Sclerosis.

OBJECTIVE: To assess whether biological aging as measured by leukocyte telomere length (LTL) is associated with clinical disability and brain volume loss in multiple sclerosis (MS). METHODS: Adults with MS/clinically isolated syndrome in the University of California, San Francisco EPIC cohort study were included. LTL was measured on DNA samples by quantitative polymerase chain reaction and expressed as telomere to somatic DNA (T/S) ratio. Expanded Disability Status Scale (EDSS) and 3-dimensional T1-weighted brain magnetic resonance imaging were performed at baseline and follow-up. Associations of baseline LTL with cross-sectional and longitudinal outcomes were assessed using simple and mixed effects linear regression models. A subset (n = 46) had LTL measured over time, and we assessed the association of LTL change with EDSS change with mixed effects models. RESULTS: Included were 356 women and 160 men (mean age = 43 years, median disease duration = 6 years, median EDSS = 1.5 [range = 0-7], mean T/S ratio = 0.97 [standard deviation = 0.18]). In baseline analyses adjusted for age, disease duration, and sex, for every 0.2 lower LTL, EDSS was 0.27 higher (95% confidence interval [CI] = 0.13-0.42, p < 0.001) and brain volume was 7.4mm3 lower (95% CI = 0.10-14.7, p = 0.047). In longitudinal adjusted analyses, those with lower baseline LTL had higher EDSS and lower brain volumes over time. In adjusted analysis of the subset, LTL change was associated with EDSS change over 10 years; for every 0.2 LTL decrease, EDSS was 0.34 higher (95% CI = 0.08-0.61, p = 0.012). INTERPRETATION: Shorter telomere length was associated with disability independent of chronological age, suggesting that biological aging may contribute to neurological injury in MS. Targeting aging-related mechanisms is a potential therapeutic strategy against MS progression. ANN NEUROL 2019;86:671-682.

Association Between Serum Neurofilament Light Chain Levels and Long-term Disease Course Among Patients With Multiple Sclerosis Followed up for 12 Years.

IMPORTANCE: Blood sample-based biomarkers that are associated with clinically meaningful outcomes for patients with multiple sclerosis (MS) have not been developed. OBJECTIVE: To evaluate the potential of serum neurofilament light chain (sNFL) measurements as a biomarker of disease activity and progression in a longitudinal MS data set. DESIGN, SETTING, AND PARTICIPANTS: Single-center, ongoing, prospective observational cohort study of 607 patients with MS from the longitudinal EPIC (Expression, Proteomics, Imaging, Clinical) study at the University of California, San Francisco from July 1, 2004, through August 31, 2017. Clinical evaluations and sample collection were performed annually for 5 years, then at different time points for up to 12 years, with a median follow-up duration of 10 (interquartile range, 7-11) years. Serum NFL levels were measured using a sensitive single molecule array platform and compared with clinical and magnetic resonance imaging variables with the use of univariable and multivariable analyses. MAIN OUTCOMES AND MEASURES: The main outcomes were disability progression defined as clinically significant worsening on the Expanded Disability Status Scale (EDSS) score and brain fraction atrophy. RESULTS: Mean (SD) age of the 607 study participants at study entry was 42.5 (9.8) years; 423 (69.7%) were women; and all participants were of non-Hispanic European descent. Of 3911 samples sequentially collected, 3904 passed quality control for quantification of sNFL. Baseline sNFL levels showed significant associations with EDSS score (ß, 1.080; 95% CI, 1.047-1.114; P < .001), MS subtype (ß, 1.478; 95% CI, 1.279-1.707; P < .001), and treatment status (ß, 1.120; 95% CI, 1.007-1.245; P = .04). A significant interaction between EDSS worsening and change in levels of sNFL over time was found (ß, 1.015; 95% CI, 1.007-1.023; P < .001). Baseline sNFL levels alone were associated with approximately 11.6% of the variance in brain fraction atrophy at year 10. In a multivariable analysis that considered sex, age, and disease duration, baseline sNFL levels were associated with 18.0% of the variance in brain fraction atrophy at year 10. After 5 years' follow-up, active treatment was associated with lower levels of sNFL, with high-potency treatments associated with the greater decreases in sNFL levels compared with platform therapies (high-potency vs untreated: ß, 0.946; 95% CI, 0.915-0.976; P < .001; high-potency vs platform: ß, 0.972; 95% CI, 0.948-0.998; P = .04). CONCLUSIONS AND RELEVANCE: This study found that statistically significant associations of sNFL with relevant clinical and neuroimaging outcomes in MS were confirmed and extended, supporting the potential of sNFL as an objective surrogate of ongoing MS disease activity. In this data set of patients with MS who received early treatment, the prognostic power of sNFL for relapse activity and long-term disability progression was limited. Further prospective studies are necessary to assess the assay's utility for decision-making in individual patients.

Next-generation sequencing reveals new information about HLA allele and haplotype diversity in a large European American population.

The human leukocyte antigen (HLA) genes are extremely polymorphic and are useful molecular markers to make inferences about human population history. However, the accuracy of the estimation of genetic diversity at HLA loci very much depends on the technology used to characterize HLA alleles; high-resolution genotyping of long-range HLA gene products improves the assessment of HLA population diversity as well as other population parameters compared to lower resolution typing methods. In this study we examined allelic and haplotype HLA diversity in a large healthy European American population sourced from the UCSF-DNA bank. A high-resolution next-generation sequencing method was applied to define non-ambiguous 3- and 4-field alleles at the HLA-A, HLA-C, HLA-B, HLA-DRB1, HLA-DRB3/4/5, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1 loci in samples provided by 2248 unrelated individuals. A number of population parameters were examined including balancing selection and various measurements of linkage disequilibrium were calculated. There were no detectable deviations from Hardy-Weinberg proportions at HLA-A, HLA-DRB1, HLA-DQA1 and HLA-DQB1. For the remaining loci moderate and significant deviations were detected at HLA-C, HLA-B, HLA-DRB3/4/5, HLA-DPA1 and HLA-DPB1 loci mostly from population substructures. Unique 4-field associations were observed among alleles at 2 loci and haplotypes extending large intervals that were not apparent in results obtained using testing methodologies with limited sequence coverage and phasing. The high diversity at HLA-DPA1 results from detection of intron variants of otherwise well conserved protein sequences. It may be speculated that divergence in exon sequences may be negatively selected. Our data provides a valuable reference source for future population studies that may allow for precise fine mapping of coding and non-coding sequences determining disease susceptibility and allo-immunogenicity.

miRNA contributions to pediatric-onset multiple sclerosis inferred from GWAS.

OBJECTIVE: Onset of multiple sclerosis (MS) occurs in childhood for approximately 5% of cases (pediatric MS, or ped-MS). Epigenetic influences are strongly implicated in MS pathogenesis in adults, including the contribution from microRNAs (miRNAs), small noncoding RNAs that affect gene expression by binding target gene mRNAs. Few studies have specifically examined miRNAs in ped-MS, but individuals developing MS at an early age may carry a relatively high burden of genetic risk factors, and miRNA dysregulation may therefore play a larger role in the development of ped-MS than in adult-onset MS. This study aimed to look for evidence of miRNA involvement in ped-MS pathogenesis. METHODS: GWAS results from 486 ped-MS cases and 1362 controls from the U.S. Pediatric MS Network and Kaiser Permanente Northern California membership were investigated for miRNA-specific signals. First, enrichment of miRNA-target gene network signals was evaluated using MIGWAS software. Second, SNPs in miRNA genes and in target gene binding sites (miR-SNPs) were tested for association with ped-MS, and pathway analysis was performed on associated target genes. RESULTS: MIGWAS analysis showed that miRNA-target gene signals were enriched in GWAS (P = 0.038) and identified 39 candidate biomarker miRNA-target gene pairs, including immune and neuronal signaling genes. The miR-SNP analysis implicated dysregulation of miRNA binding to target genes in five pathways, mainly involved in immune signaling. INTERPRETATION: Evidence from GWAS suggests that miRNAs play a role in ped-MS pathogenesis by affecting immune signaling and other pathways. Candidate biomarker miRNA-target gene pairs should be further studied for diagnostic, prognostic, and/or therapeutic utility.

Sex-specific Tau methylation patterns and synaptic transcriptional alterations are associated with neural vulnerability during chronic neuroinflammation.

The molecular events underlying the transition from initial inflammatory flares to the progressive phase of multiple sclerosis (MS) remain poorly understood. Here, we report that the microtubule-associated protein (MAP) Tau exerts a gender-specific protective function on disease progression in the MS model experimental autoimmune encephalomyelitis (EAE). A detailed investigation of the autoimmune response in Tau-deficient mice excluded a strong immunoregulatory role for Tau, suggesting that its beneficial effects are presumably exerted within the central nervous system (CNS). Spinal cord transcriptomic data show increased synaptic dysfunctions and alterations in the NF-kB activation pathway upon EAE in Tau-deficient mice as compared to wildtype animals. We also performed the first comprehensive characterization of Tau post-translational modifications (PTMs) in the nervous system upon EAE. We report that the methylation levels of the conserved lysine residue K306 are significantly decreased in the chronic phase of the disease. By combining biochemical assays and molecular dynamics (MD) simulations, we demonstrate that methylation at K306 decreases the affinity of Tau for the microtubule network. Thus, the down-regulation of this PTM might represent a homeostatic response to enhance axonal stability against an autoimmune CNS insult. The results, altogether, position Tau as key mediator between the inflammatory processes and neurodegeneration that seems to unify many CNS diseases.

A specific amino acid motif of HLA-DRB1 mediates risk and interacts with smoking history in Parkinson's disease.

Parkinson's disease (PD) is a neurodegenerative disease in which genetic risk has been mapped to HLA, but precise allelic associations have been difficult to infer due to limitations in genotyping methodology. Mapping PD risk at highest possible resolution, we performed sequencing of 11 HLA genes in 1,597 PD cases and 1,606 controls. We found that susceptibility to PD can be explained by a specific combination of amino acids at positions 70-74 on the HLA-DRB1 molecule. Previously identified as the primary risk factor in rheumatoid arthritis and referred to as the "shared epitope" (SE), the residues Q/R-K/R-R-A-A at positions 70-74 in combination with valine at position 11 (11-V) is highly protective in PD, while risk is attributable to the identical epitope in the absence of 11-V. Notably, these effects are modified by history of cigarette smoking, with a strong protective effect mediated by a positive history of smoking in combination with the SE and 11-V (P = 10-4; odds ratio, 0.51; 95% confidence interval, 0.36-0.72) and risk attributable to never smoking in combination with the SE without 11-V (P = 0.01; odds ratio, 1.51; 95% confidence interval, 1.08-2.12). The association of specific combinations of amino acids that participate in critical peptide-binding pockets of the HLA class II molecule implicates antigen presentation in PD pathogenesis and provides further support for genetic control of neuroinflammation in disease. The interaction of HLA-DRB1 with smoking history in disease predisposition, along with predicted patterns of peptide binding to HLA, provide a molecular model that explains the unique epidemiology of smoking in PD.

Deconstruction of HLA-DRB1*04:01:01 and HLA-DRB1*15:01:01 class II haplotypes using next-generation sequencing in European-Americans with multiple sclerosis.

BACKGROUND: The association between HLA-DRB1*15:01 with multiple sclerosis (MS) susceptibility is well established, but the contribution of the tightly associated HLA-DRB5*01:01 allele has not yet been completely ascertained. Similarly, the effects of HLA-DRB1*04:01 alleles and haplotypes, defined at the full-gene resolution level with MS risk remains to be elucidated. OBJECTIVES: To characterize the molecular architecture of class II HLA-DR15 and HLA-DR4 haplotypes associated with MS. METHODS: Next-generation sequencing was used to determine HLA-DQB1, HLA-DQA1, and HLA-DRB1/4/5 alleles in 1403 unrelated European-American patients and 1425 healthy unrelated controls. Effect sizes of HLA alleles and haplotypes on MS risk were measured by odds ratio (OR) with 95% confidence intervals. RESULTS: HLA-DRB1*15:01:01:01SG (OR = 3.20, p < 2.2E-16), HLA-DRB5*01:01:01 (OR = 2.96, p < 2.2E-16), and HLA-DRB5*01:01:01v1_STR1 (OR = 8.18, p = 4.3E-05) alleles all occurred at significantly higher frequencies in MS patients compared to controls. The most significant predis-posing haplotypes were HLA-DQB1*06:02:01~ HLA-DQA1*01:02:01:01SG~HLA-DRB1*15:01:01:01SG~HLA-DRB5*01:01:01 and HLA-DQB1*06:02:01~HLA-DQA1*01:02:01:01SG~HLA-DRB1*15:01:01:01SG~HLA-DRB5*01:01:01v1_STR1 (OR = 3.19, p < 2.2E-16; OR = 9.30, p = 9.7E-05, respectively). Analyses of the HLA-DRB1*04 cohort in the absence of HLA-DRB1*15:01 haplotypes revealed that the HLA-DQB1*03:01:01:01~HLA-DQA1*03:03:01:01~HLA-DRB1*04:01:01:01SG~HLA-DRB4*01:03:01:01 haplotype was protective (OR = 0.64, p = 0.028), whereas the HLA-DQB1*03:02:01~HLA-DQA1*03:01:01~HLA-DRB1*04:01:01:01SG~HLA-DRB4*01:03:01:01 haplotype was associated with MS susceptibility (OR = 1.66, p = 4.9E-03). CONCLUSION: HLA-DR15 haplotypes, including genomic variants of HLA-DRB5, and HLA-DR4 haplotypes affect MS risk.

Genome sequencing uncovers phenocopies in primary progressive multiple sclerosis.

OBJECTIVE: Primary progressive multiple sclerosis (PPMS) causes accumulation of neurological disability from disease onset without clinical attacks typical of relapsing multiple sclerosis (RMS). However, whether genetic variation influences the disease course remains unclear. We aimed to determine whether mutations causative of neurological disorders that share features with multiple sclerosis (MS) contribute to risk for developing PPMS. METHODS: We examined whole-genome sequencing (WGS) data from 38 PPMS and 81 healthy subjects of European ancestry. We selected pathogenic variants exclusively found in PPMS patients that cause monogenic neurological disorders and performed two rounds of replication genotyping in 746 PPMS, 3,049 RMS, and 1,000 healthy subjects. To refine our findings, we examined the burden of rare, potentially pathogenic mutations in 41 genes that cause hereditary spastic paraplegias (HSPs) in PPMS (n = 314), secondary progressive multiple sclerosis (SPMS; n = 587), RMS (n = 2,248), and healthy subjects (n = 987) genotyped using the MS replication chip. RESULTS: WGS and replication studies identified three pathogenic variants in PPMS patients that cause neurological disorders sharing features with MS: KIF5A p.Ala361Val in spastic paraplegia 10; MLC1 p.Pro92Ser in megalencephalic leukodystrophy with subcortical cysts, and REEP1 c.606 + 43G>T in Spastic Paraplegia 31. Moreover, we detected a significant enrichment of HSP-related mutations in PPMS patients compared to controls (risk ratio [RR] = 1.95; 95% confidence interval [CI], 1.27-2.98; p = 0.002), as well as in SPMS patients compared to controls (RR = 1.57; 95% CI, 1.18-2.10; p = 0.002). Importantly, this enrichment was not detected in RMS. INTERPRETATION: This study provides evidence to support the hypothesis that rare Mendelian genetic variants contribute to the risk for developing progressive forms of MS. Ann Neurol 2018;83:51-63.