SummaryThe immune system is highly diverse, but characterization of its genetic architecture has lagged behind the vast progress made by genome-wide association studies (GWASs) of emergent diseases. Our GWAS for 54 functionally relevant phenotypes of the adaptive immune system in 489 healthy individuals identifies eight genome-wide significant associations explaining 6%–20% of variance. Coding and splicing variants in PTPRC and COMMD10 are involved in memory T cell differentiation. Genetic variation controlling disease-relevant T helper cell subsets includes RICTOR and STON2 associated with Th2 and Th17, respectively, and the interferon-lambda locus controlling regulatory T cell proliferation. Early and memory B cell differentiation stages are associated with variation in LARP1B and SP4. Finally, the latrophilin family member ADGRL2 correlates with baseline pro-inflammatory interleukin-6 levels. Suggestive associations reveal mechanisms of autoimmune disease associations, in particular related to pro-inflammatory cytokine production. Pinpointing these key human immune regulators offers attractive therapeutic perspectives.

Increasing evidence implicates B cells in the pathogenesis of multiple sclerosis. Smets et al. report that the CD40 multiple sclerosis risk allele lowers CD40 expression, whereas the CD86 risk allele increases CD86 expression. B cells may have an important antigen presentation and immunoregulatory role in multiple sclerosis.

Objective Evidence for a role of microglia in the pathogenesis of multiple sclerosis (MS) is growing. We investigated association of microglial markers at time of diagnostic lumbar puncture (LP) with different aspects of disease activity (relapses, disability, magnetic resonance imaging parameters) up to 6 years later in a cohort of 143 patients. Methods In cerebrospinal fluid (CSF), we measured 3 macrophage and microglia‐related proteins, chitotriosidase (CHIT1), chitinase‐3–like protein 1 (CHI3L1 or YKL‐40), and soluble triggering receptor expressed on myeloid cells 2 (sTREM2), as well as a marker of neuronal damage, neurofilament light chain (NfL), using enzyme‐linked immunosorbent assay and electrochemiluminescence. We investigated the same microglia‐related markers in publicly available RNA expression data from postmortem brain tissue. Results CHIT1 levels at diagnostic LP correlated with 2 aspects of long‐term disease activity after correction for multiple testing. First, CHIT1 increased with reduced tissue integrity in lesions at a median 3 years later (p = 9.6E‐04). Second, CHIT1 reflected disease severity at a median 5 years later (p = 1.2E‐04). Together with known clinical covariates, CHIT1 levels explained 12% and 27% of variance in these 2 measures, respectively, and were able to distinguish slow and fast disability progression (area under the curve = 85%). CHIT1 was the best discriminator of chronic active versus chronic inactive lesions and the only marker correlated with NfL (r = 0.3, p = 0.0019). Associations with disease activity were, however, independent of NfL. Interpretation CHIT1 CSF levels measured during the diagnostic LP reflect microglial activation early on in MS and can be considered a valuable prognostic biomarker for future disease activity. ANN NEUROL 2020;87:633–645

LncRNA PCR arrays containing 90 common LncRNAs were used to screen lncRNA expression levels in PBMC from a discovery population of patients with MS. Data from discovery and replications cohorts showed a generalized dysregulation of lncRNA levels in MS patients compared with controls. MALAT1, MEG9, NRON, ANRIL, TUG1, XIST, SOX2OT, GOMAFU, HULC, BACE-1AS were significantly downregulated in MS patients in comparison with controls. Therefore, we performed a validation analysis in an independent cohort of Belgian origin. In this study, NRON and TUG1 downregulations in MS patients compared with controls were confirmed (p ≤ .05 and p ≤ .0001 respectively), whereas considering the other lncRNAs, the statistical threshold was not reached. LncRNAs profiling could thus represent a new challenge in the research of easy detectable biomarkers of disease susceptibility and progression.

The role of somatic variants in diseases beyond cancer is increasingly being recognized, with potential roles in autoinflammatory and autoimmune diseases. However, as mutation rates and allele fractions are lower, studies in these diseases are substantially less tolerant of false positives, and bio-informatics algorithms require high replication rates. We developed a pipeline combining two variant callers, MuTect2 and VarScan2, with technical filtering and prioritization. Our pipeline detects somatic variants with allele fractions as low as 0.5% and achieves a replication rate of >55%. Validation in an independent data set demonstrates excellent performance (sensitivity > 57%, specificity > 98%, replication rate > 80%). We applied this pipeline to the autoimmune disease multiple sclerosis (MS) as a proof-of-principle. We demonstrate that 60% of MS patients carry 2–10 exonic somatic variants in their peripheral blood T and B cells, with the vast majority (80%) occurring in T cells and variants persisting over time. Synonymous variants significantly co-occur with non-synonymous variants. Systematic characterization indicates somatic variants are enriched for being novel or very rare in public databases of germline variants and trend towards being more damaging and conserved, as reflected by higher phred-scaled combined annotation-dependent depletion (CADD) and genomic evolutionary rate profiling (GERP) scores. Our pipeline and proof-of-principle now warrant further investigation of common somatic genetic variation on top of inherited genetic variation in the context of autoimmune disease, where it may offer subtle survival advantages to immune cells and contribute to the capacity of these cells to participate in the autoimmune reaction.

Objective: Many multiple sclerosis (MS) genetic susceptibility variants have been identified, but understanding disease heterogeneity remains a key challenge. Relapses are a core feature of MS and a common primary outcome of clinical trials, with prevention of relapses benefiting patients immediately and potentially limiting long-term disability accrual. We aim to identify genetic variation associated with relapse hazard in MS by analyzing the largest study population to date. Methods: We performed a genomewide association study (GWAS) in a discovery cohort and investigated the genomewide significant variants in a replication cohort. Combining both cohorts, we captured a total of 2,231 relapses occurring before the start of any immunomodulatory treatment in 991 patients. For assessing time to relapse, we applied a survival analysis utilizing Cox proportional hazards models. We also investigated the association between MS genetic risk scores and relapse hazard and performed a gene ontology pathway analysis. Results: The low-frequency genetic variant rs11871306 within WNT9B reached genomewide significance in predicting relapse hazard and replicated (meta-analysis hazard ratio (HR) = 2.15, 95% confidence interval (CI) = 1.70-2.78, p = 2.07 × 10 −10 ). A pathway analysis identified an association of the pathway "response to vitamin D" with relapse hazard (p = 4.33 × 10 −6 ). The MS genetic risk scores, however, were not associated with relapse hazard. Interpretation: Genetic factors underlying disease heterogeneity differ from variants associated with MS susceptibility. Our findings imply that genetic variation within the Wnt signaling and vitamin D pathways contributes to differences in relapse occurrence. The present study highlights these cross-talking pathways as potential modulators of MS disease activity.

Variant rs12988804 in LRP2, the first example of a genome-wide significant association with relapse rate in MS, is replicated in an independent study.

Although fingolimod and interferon-β are two mechanistically different multiple sclerosis (MS) treatments, they both induce B cell activating factor (BAFF) and shift the B cell pool towards a regulatory phenotype. However, whether there is a shared mechanism between both treatments in how they influence the B cell compartment remains elusive. In this study, we collected a cross-sectional study population of 112 MS patients (41 untreated, 42 interferon-β, 29 fingolimod) and determined B cell subsets, cell-surface and RNA expression of BAFF-receptor (BAFF-R) and transmembrane activator and cyclophilin ligand interactor (TACI) as well as plasma and/or RNA levels of BAFF, BAFF splice forms and interleukin-10 (IL-10) and -35 (IL-35). We added an in vitro B cell culture with four stimulus conditions (Medium, CpG, BAFF and CpG+BAFF) for untreated and interferon-β treated patients including measurement of intracellular IL-10 levels. Our flow experiments showed that interferon-β and fingolimod induced BAFF protein and mRNA expression (P ≤ 3.15 x 10-4) without disproportional change in the antagonizing splice form. Protein BAFF correlated with an increase in transitional B cells (P = 5.70 x 10-6), decrease in switched B cells (P = 3.29 x 10-4), and reduction in B cell-surface BAFF-R expression (P = 2.70 x 10-10), both on TACI-positive and -negative cells. TACI and BAFF-R RNA levels remained unaltered. RNA, plasma and in vitro experiments demonstrated that BAFF was not associated with increased IL-10 and IL-35 levels. In conclusion, treatment-induced BAFF correlates with a shift towards transitional B cells which are enriched for cells with an immunoregulatory function. However, BAFF does not directly influence the expression of the immunoregulatory cytokines IL-10 and IL-35. Furthermore, the post-translational mechanism of BAFF-induced BAFF-R cell surface loss was TACI-independent. These observations put the failure of pharmaceutical anti-BAFF strategies in perspective and provide insights for targeted B cell therapies.