A robust pipeline with high replication rate for detection of somatic variants in the adaptive immune system as a source of common genetic variation in autoimmune disease

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.Peer reviewe

Multiple sclerosis risk variants alter expression of co-stimulatory genes in B cells.

The increasing evidence supporting a role for B cells in the pathogenesis of multiple sclerosis prompted us to investigate the influence of known susceptibility variants on the surface expression of co-stimulatory molecules in these cells. Using flow cytometry we measured surface expression of CD40 and CD86 in B cells from 68 patients and 162 healthy controls that were genotyped for the multiple sclerosis associated single nucleotide polymorphisms (SNPs) rs4810485, which maps within the CD40 gene, and rs9282641, which maps within the CD86 gene. We found that carrying the risk allele rs4810485*T lowered the cell-surface expression of CD40 in all tested B cell subtypes (in total B cells P ≤ 5.10 × 10-5 in patients and ≤4.09 × 10-6 in controls), while carrying the risk allele rs9282641*G increased the expression of CD86, with this effect primarily seen in the naïve B cell subset (P = 0.048 in patients and 5.38 × 10-5 in controls). In concordance with these results, analysis of RNA expression demonstrated that the risk allele rs4810485*T resulted in lower total CD40 expression (P = 0.057) but with an increased proportion of alternative splice-forms leading to decoy receptors (P = 4.00 × 10-7). Finally, we also observed that the risk allele rs4810485*T was associated with decreased levels of interleukin-10 (P = 0.020), which is considered to have an immunoregulatory function downstream of CD40. Given the importance of these co-stimulatory molecules in determining the immune reaction that appears in response to antigen our data suggest that B cells might have an important antigen presentation and immunoregulatory role in the pathogenesis of multiple sclerosis

CHIT1 at Diagnosis Reflects Long-Term Multiple Sclerosis Disease Activity

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.status: publishe

Genetic Variation in WNT9B Increases Relapse Hazard in Multiple Sclerosis

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 x 10(-10)). A pathway analysis identified an association of the pathway response to vitamin D with relapse hazard (p = 4.33 x 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. ANN NEUROL 202

A robust pipeline with high replication rate for detection of somatic variants in the adaptive immune system as a source of common genetic variation in autoimmune disease.

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 >55%. Validation in an independent dataset 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 nonsynonymous 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 CADD and 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.status: publishe

LncRNAs expression profile in peripheral blood mononuclear cells from multiple sclerosis patients

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.status: publishe

Multiple sclerosis risk variants alter expression of co-stimulatory genes in B cells

The increasing evidence supporting a role for B cells in the pathogenesis of multiple sclerosis prompted us to investigate the influence of known susceptibility variants on the surface expression of co-stimulatory molecules in these cells. Using flow cytometry we measured surface expression of CD40 and CD86 in B cells from 68 patients and 162 healthy controls that were genotyped for the multiple sclerosis associated single nucleotide polymorphisms (SNPs) rs4810485, which maps within the CD40 gene, and rs9282641, which maps within the CD86 gene. We found that carrying the risk allele rs4810485*T lowered the cell-surface expression of CD40 in all tested B cell subtypes (in total B cells P ≤ 5.10 × 10-5 in patients and ≤4.09 × 10-6 in controls), while carrying the risk allele rs9282641*G increased the expression of CD86, with this effect primarily seen in the naïve B cell subset (P = 0.048 in patients and 5.38 × 10-5 in controls). In concordance with these results, analysis of RNA expression demonstrated that the risk allele rs4810485*T resulted in lower total CD40 expression (P = 0.057) but with an increased proportion of alternative splice-forms leading to decoy receptors (P = 4.00 × 10-7). Finally, we also observed that the risk allele rs4810485*T was associated with decreased levels of interleukin-10 (P = 0.020), which is considered to have an immunoregulatory function downstream of CD40. Given the importance of these co-stimulatory molecules in determining the immune reaction that appears in response to antigen our data suggest that B cells might have an important antigen presentation and immunoregulatory role in the pathogenesis of multiple sclerosis.status: publishe

Genetic Architecture of Adaptive Immune System Identifies Key Immune Regulators

The 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.status: publishe

Genetic Architecture of Adaptive Immune System Identifies Key Immune Regulators

Summary: The 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. : Lagou et al. identify genetic factors explaining interindividual variation in composition of the adaptive immune system. Factors pinpoint key human immune regulators controlling B and T cell differentiation and levels of disease-relevant T helper and regulatory cells. These findings shed light on mechanisms of autoimmune disease and offer therapeutic perspectives. Keywords: adaptive immune system, immune phenotype, genetics, association, genome-wide association, autoimmunity, susceptibilit