Human immune diversity: from evolution to modernity.

The extreme diversity of the human immune system, forged and maintained throughout evolutionary history, provides a potent defense against opportunistic pathogens. At the same time, this immune variation is the substrate upon which a plethora of immune-associated diseases develop. Genetic analysis suggests that thousands of individually weak loci together drive up to half of the observed immune variation. Intense selection maintains this genetic diversity, even selecting for the introgressed Neanderthal or Denisovan alleles that have reintroduced variation lost during the out-of-Africa migration. Variations in age, sex, diet, environmental exposure, and microbiome each potentially explain the residual variation, with proof-of-concept studies demonstrating both plausible mechanisms and correlative associations. The confounding interaction of many of these variables currently makes it difficult to assign definitive contributions. Here, we review the current state of play in the field, identify the key unknowns in the causality of immune variation, and identify the multidisciplinary pathways toward an improved understanding.

The cellular composition of the human immune system is shaped by age and cohabitation.

Detailed population-level description of the human immune system has recently become achievable. We used a 'systems-level' approach to establish a resource of cellular immune profiles of 670 healthy individuals. We report a high level of interindividual variation, with low longitudinal variation, at the level of cellular subset composition of the immune system. Despite the profound effects of antigen exposure on individual antigen-specific clones, the cellular subset structure proved highly elastic, with transient vaccination-induced changes followed by a return to the individual's unique baseline. Notably, the largest influence on immunological variation identified was cohabitation, with 50% less immunological variation between individuals who share an environment (as parents) than between people in the wider population. These results identify local environmental conditions as a key factor in shaping the human immune system.

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 × 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. ANN NEUROL 2021;89:884-894.

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.

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.

IFN-γ stimulates CpG-induced IL-10 production in B cells via p38 and JNK signalling pathways.

The production of IL-10, a potent immunosuppressive cytokine, must be strictly regulated to ensure a balanced immune response. IFN-γ, a key cytokine in multiple immune processes and pathologies, is known as an inhibitor of IL-10 production by monocytes and macrophages, but also has some regulatory functions. In the present study, we explored the role of IFN-γ on Toll-like receptor (TLR)-induced IL-10 production in murine peritoneal and spleen cells and in human peripheral blood mononuclear cells. IFN-γ inhibited IL-10 production induced by TLR2, TLR3, TLR4 and TLR7/8 agonists, but stimulated IL-10 production when cells were triggered with CpG oligodeoxynucleotides, a specific TLR9 agonist. The stimulatory effect of IFN-γ on TLR9-induced IL-10 was restricted to B cells. In line with the increased IL-10, B cells stimulated with CpG and IFN-γ profoundly inhibited CD4 T cell proliferation. Further research into the mechanisms involved, revealed that the mitogen-activated protein kinases p38 and JNK are essential players in this stimulatory effect, and that the phosphatase MKP1 - an inhibitor of p38 and JNK activity - is downregulated after combined stimulation with IFN-γ and CpG. Our data may represent a novel immunoregulatory role of IFN-γ in B cells after triggering of TLR9, by stimulating IL-10 production.

Machine learning identifies an immunological pattern associated with multiple juvenile idiopathic arthritis subtypes.

OBJECTIVES: Juvenile idiopathic arthritis (JIA) is the most common class of childhood rheumatic diseases, with distinct disease subsets that may have diverging pathophysiological origins. Both adaptive and innate immune processes have been proposed as primary drivers, which may account for the observed clinical heterogeneity, but few high-depth studies have been performed. METHODS: Here we profiled the adaptive immune system of 85 patients with JIA and 43 age-matched controls with indepth flow cytometry and machine learning approaches. RESULTS: Immune profiling identified immunological changes in patients with JIA. This immune signature was shared across a broad spectrum of childhood inflammatory diseases. The immune signature was identified in clinically distinct subsets of JIA, but was accentuated in patients with systemic JIA and those patients with active disease. Despite the extensive overlap in the immunological spectrum exhibited by healthy children and patients with JIA, machine learning analysis of the data set proved capable of discriminating patients with JIA from healthy controls with ~90% accuracy. CONCLUSIONS: These results pave the way for large-scale immune phenotyping longitudinal studies of JIA. The ability to discriminate between patients with JIA and healthy individuals provides proof of principle for the use of machine learning to identify immune signatures that are predictive to treatment response group.

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.

Genetic ablation of IP3 receptor 2 increases cytokines and decreases survival of SOD1G93A mice.

Amyotrophic lateral sclerosis (ALS) is a devastating progressive neurodegenerative disease characterized by the selective death of motor neurons. Disease pathophysiology is complex and not yet fully understood. Higher gene expression of the inositol 1,4,5-trisphosphate receptor 2 gene (ITPR2), encoding the IP3 receptor 2 (IP3R2), was detected in sporadic ALS patients. Here, we demonstrate that IP3R2 gene expression was also increased in spinal cords of ALS mice. Moreover, an increase of IP3R2 expression was observed in other models of chronic and acute neurodegeneration. Upregulation of IP3R2 gene expression could be induced by lipopolysaccharide (LPS) in murine astrocytes, murine macrophages and human fibroblasts indicating that it may be a compensatory response to inflammation. Preventing this response by genetic deletion of ITPR2 from SOD1G93A mice had a dose-dependent effect on disease duration, resulting in a significantly shorter lifespan of these mice. In addition, the absence of IP3R2 led to increased innate immunity, which may contribute to the decreased survival of the SOD1G93A mice. Besides systemic inflammation, IP3R2 knockout mice also had increased IFNγ, IL-6 and IL1α expression. Altogether, our data indicate that IP3R2 protects against the negative effects of inflammation, suggesting that the increase in IP3R2 expression in ALS patients is a protective response.

Genetic basis for relapse rate in multiple sclerosis: Association with LRP2 genetic variation.

BACKGROUND: In contrast to successes for multiple sclerosis (MS) susceptibility, the genetic basis for clinical heterogeneity remains largely unresolved. OBJECTIVES: We investigate the first reported genetic association with relapse rate. METHODS: We genotyped variant rs12988804 in LRP2 in a homogeneous study population of 527 Belgian MS patients with 970 documented relapses. RESULTS: The rs12988804*T allele is associated with a 1.16-fold increased hazard rate for a relapse occurring ( P = 0.0078) and a higher baseline relapse rate prior to immunomodulatory treatment ( P = 0.044). CONCLUSION: Variant rs12988804 in LRP2, the first example of a genome-wide significant association with relapse rate in MS, is replicated in an independent study.

Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis.

Multiple sclerosis is a common disease of the central nervous system in which the interplay between inflammatory and neurodegenerative processes typically results in intermittent neurological disturbance followed by progressive accumulation of disability. Epidemiological studies have shown that genetic factors are primarily responsible for the substantially increased frequency of the disease seen in the relatives of affected individuals, and systematic attempts to identify linkage in multiplex families have confirmed that variation within the major histocompatibility complex (MHC) exerts the greatest individual effect on risk. Modestly powered genome-wide association studies (GWAS) have enabled more than 20 additional risk loci to be identified and have shown that multiple variants exerting modest individual effects have a key role in disease susceptibility. Most of the genetic architecture underlying susceptibility to the disease remains to be defined and is anticipated to require the analysis of sample sizes that are beyond the numbers currently available to individual research groups. In a collaborative GWAS involving 9,772 cases of European descent collected by 23 research groups working in 15 different countries, we have replicated almost all of the previously suggested associations and identified at least a further 29 novel susceptibility loci. Within the MHC we have refined the identity of the HLA-DRB1 risk alleles and confirmed that variation in the HLA-A gene underlies the independent protective effect attributable to the class I region. Immunologically relevant genes are significantly overrepresented among those mapping close to the identified loci and particularly implicate T-helper-cell differentiation in the pathogenesis of multiple sclerosis.

Power estimation for non-standardized multisite studies.

A concern for researchers planning multisite studies is that scanner and T1-weighted sequence-related biases on regional volumes could overshadow true effects, especially for studies with a heterogeneous set of scanners and sequences. Current approaches attempt to harmonize data by standardizing hardware, pulse sequences, and protocols, or by calibrating across sites using phantom-based corrections to ensure the same raw image intensities. We propose to avoid harmonization and phantom-based correction entirely. We hypothesized that the bias of estimated regional volumes is scaled between sites due to the contrast and gradient distortion differences between scanners and sequences. Given this assumption, we provide a new statistical framework and derive a power equation to define inclusion criteria for a set of sites based on the variability of their scaling factors. We estimated the scaling factors of 20 scanners with heterogeneous hardware and sequence parameters by scanning a single set of 12 subjects at sites across the United States and Europe. Regional volumes and their scaling factors were estimated for each site using Freesurfer's segmentation algorithm and ordinary least squares, respectively. The scaling factors were validated by comparing the theoretical and simulated power curves, performing a leave-one-out calibration of regional volumes, and evaluating the absolute agreement of all regional volumes between sites before and after calibration. Using our derived power equation, we were able to define the conditions under which harmonization is not necessary to achieve 80% power. This approach can inform choice of processing pipelines and outcome metrics for multisite studies based on scaling factor variability across sites, enabling collaboration between clinical and research institutions.

No evidence for shared genetic basis of common variants in multiple sclerosis and amyotrophic lateral sclerosis.

Genome-wide association studies have been successful in identifying common variants that influence the susceptibility to complex diseases. From these studies, it has emerged that there is substantial overlap in susceptibility loci between diseases. In line with those findings, we hypothesized that shared genetic pathways may exist between multiple sclerosis (MS) and amyotrophic lateral sclerosis (ALS). While both diseases may have inflammatory and neurodegenerative features, epidemiological studies have indicated an increased co-occurrence within individuals and families. To this purpose, we combined genome-wide data from 4088 MS patients, 3762 ALS patients and 12 030 healthy control individuals in whom 5 440 446 single-nucleotide polymorphisms (SNPs) were successfully genotyped or imputed. We tested these SNPs for the excess association shared between MS and ALS and also explored whether polygenic models of SNPs below genome-wide significance could explain some of the observed trait variance between diseases. Genome-wide association meta-analysis of SNPs as well as polygenic analyses fails to provide evidence in favor of an overlap in genetic susceptibility between MS and ALS. Hence, our findings do not support a shared genetic background of common risk variants in MS and ALS.

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

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

Interleukin 7 receptor alpha chain (IL7R) shows allelic and functional association with multiple sclerosis.

Multiple sclerosis is a demyelinating neurodegenerative disease with a strong genetic component. Previous genetic risk studies have failed to identify consistently linked regions or genes outside of the major histocompatibility complex on chromosome 6p. We describe allelic association of a polymorphism in the gene encoding the interleukin 7 receptor alpha chain (IL7R) as a significant risk factor for multiple sclerosis in four independent family-based or case-control data sets (overall P = 2.9 x 10(-7)). Further, the likely causal SNP, rs6897932, located within the alternatively spliced exon 6 of IL7R, has a functional effect on gene expression. The SNP influences the amount of soluble and membrane-bound isoforms of the protein by putatively disrupting an exonic splicing silencer.

Novel COL4A1 mutations cause cerebral small vessel disease by haploinsufficiency.

Mutations in COL4A1 have been identified in families with hereditary small vessel disease of the brain presumably due to a dominant-negative mechanism. Here, we report on two novel mutations in COL4A1 in two families with porencephaly, intracerebral hemorrhage and severe white matter disease caused by haploinsufficiency. Two families with various clinical presentations of cerebral microangiopathy and autosomal dominant inheritance were examined. Clinical, neuroradiological and genetic investigations were performed. Electron microscopy of the skin was also performed. In one of the families, sequence analysis revealed a one base deletion, c.2085del, leading to a frameshift and a premature stopcodon, p.(Gly696fs). In the other family, a splice site mutation was identified, c.2194-1G>A, which most likely leads to skipping of an exon with a frameshift and premature termination as a result. In fibroblasts of affected individuals from both the families, nonsense-mediated decay (NMD) of the mutant COL4A1 messenger RNAs (mRNAs) and a clear reduction of COL4A1 protein expression were demonstrated, indicating haploinsufficiency of COL4A1. Moreover, thickening of the capillary basement membrane in the skin was documented, similar to reports in patients with COL4A1 missense mutations. These findings suggest haploinsufficiency, a different mechanism from the commonly assumed dominant-negative effect, for COL4A1 mutations as a cause of (antenatal) intracerebral hemorrhage and white matter disease.

Burden of risk variants correlates with phenotype of multiple sclerosis.

BACKGROUND: More than 100 common variants underlying multiple sclerosis (MS) susceptibility have been identified, but their effect on disease phenotype is still largely unknown. OBJECTIVE: The objective of this paper is to assess whether the cumulative genetic risk score of currently known susceptibility variants affects clinical presentation. METHODS: A cumulative genetic risk score was based on four human leukocyte antigen (HLA) and 106 non-HLA risk loci genotyped or imputed in 842 Belgian MS patients and 321 controls. Non-parametric analyses were applied. RESULTS: An increased genetic risk is observed for MS patients, including subsets such as oligoclonal band-negative and primary progressive MS patients, compared to controls. Within the patient group, a stronger association between HLA risk variants and the presence of oligoclonal bands, an increased immunoglobulin G (IgG) index and female gender was apparent. Results suggest an association between a higher accumulation of non-HLA risk variants and increased relapse rate as well as shorter relapse-free intervals after disease onset. CONCLUSION: MS patients display a significantly increased genetic risk compared to controls, irrespective of disease course or presence of oligoclonal bands. Whereas the cumulative burden of non-HLA risk variants appears to be reflected in the relapses of MS patients, the HLA region influences intrathecal IgG levels.

Whole-genome sequencing reveals a coding non-pathogenic variant tagging a non-coding pathogenic hexanucleotide repeat expansion in C9orf72 as cause of amyotrophic lateral sclerosis.

Motor neuron degeneration in amyotrophic lateral sclerosis (ALS) has a familial cause in 10% of patients. Despite significant advances in the genetics of the disease, many families remain unexplained. We performed whole-genome sequencing in five family members from a pedigree with autosomal-dominant classical ALS. A family-based elimination approach was used to identify novel coding variants segregating with the disease. This list of variants was effectively shortened by genotyping these variants in 2 additional unaffected family members and 1500 unrelated population-specific controls. A novel rare coding variant in SPAG8 on chromosome 9p13.3 segregated with the disease and was not observed in controls. Mutations in SPAG8 were not encountered in 34 other unexplained ALS pedigrees, including 1 with linkage to chromosome 9p13.2-23.3. The shared haplotype containing the SPAG8 variant in this small pedigree was 22.7 Mb and overlapped with the core 9p21 linkage locus for ALS and frontotemporal dementia. Based on differences in coverage depth of known variable tandem repeat regions between affected and non-affected family members, the shared haplotype was found to contain an expanded hexanucleotide (GGGGCC)(n) repeat in C9orf72 in the affected members. Our results demonstrate that rare coding variants identified by whole-genome sequencing can tag a shared haplotype containing a non-coding pathogenic mutation and that changes in coverage depth can be used to reveal tandem repeat expansions. It also confirms (GGGGCC)n repeat expansions in C9orf72 as a cause of familial ALS.