Genetic variants regulating immune cell levels in health and disease.

The complex network of specialized cells and molecules in the immune system has evolved to defend against pathogens, but inadvertent immune system attacks on "self" result in autoimmune disease. Both genetic regulation of immune cell levels and their relationships with autoimmunity are largely undetermined. Here, we report genetic contributions to quantitative levels of 95 cell types encompassing 272 immune traits, in a cohort of 1,629 individuals from four clustered Sardinian villages. We first estimated trait heritability, showing that it can be substantial, accounting for up to 87% of the variance (mean 41%). Next, by assessing ∼8.2 million variants that we identified and confirmed in an extended set of 2,870 individuals, 23 independent variants at 13 loci associated with at least one trait. Notably, variants at three loci (HLA, IL2RA, and SH2B3/ATXN2) overlap with known autoimmune disease associations. These results connect specific cellular phenotypes to specific genetic variants, helping to explicate their involvement in disease.

Neutralizing Antibody Responses to Viral Infections Are Linked to the Non-classical MHC Class II Gene H2-Ob.

Select humans and animals control persistent viral infections via adaptive immune responses that include production of neutralizing antibodies. The precise genetic basis for the control remains enigmatic. Here, we report positional cloning of the gene responsible for production of retrovirus-neutralizing antibodies in mice of the I/LnJ strain. It encodes the beta subunit of the non-classical major histocompatibility complex class II (MHC-II)-like molecule H2-O, a negative regulator of antigen presentation. The recessive and functionally null I/LnJ H2-Ob allele supported the production of virus-neutralizing antibodies independently of the classical MHC haplotype. Subsequent bioinformatics and functional analyses of the human H2-Ob homolog, HLA-DOB, revealed both loss- and gain-of-function alleles, which could affect the ability of their carriers to control infections with human hepatitis B (HBV) and C (HCV) viruses. Thus, understanding of the previously unappreciated role of H2-O (HLA-DO) in immunity to infections may suggest new approaches in achieving neutralizing immunity to viruses.

The MHC Class II Antigen-Processing and Presentation Pathway Is Dysregulated in Type 1 Diabetes.

Peptide loading of MHC class II (MHCII) molecules is facilitated by HLA-DM (DM), which catalyzes CLIP release, stabilizes empty MHCII, and edits the MHCII-bound peptide repertoire. HLA-DO (DO) binds to DM and modulates its activity, resulting in an altered set of peptides presented at the cell surface. MHCII-peptide presentation in individuals with type 1 diabetes (T1D) is abnormal, leading to a breakdown in tolerance; however, no direct measurement of the MHCII pathway activity in T1D patients has been performed. In this study, we measured MHCII Ag-processing pathway activity in humans by determining MHCII, MHCII-CLIP, DM, and DO levels by flow cytometry for peripheral blood B cells, dendritic cells, and monocytes from 99 T1D patients and 97 controls. Results showed that MHCII levels were similar for all three APC subsets. In contrast, MHCII-CLIP levels, independent of sex, age at blood draw, disease duration, and diagnosis age, were significantly increased for all three APCs, with B cells showing the largest increase (3.4-fold). DM and DO levels, which usually directly correlate with MHCII-CLIP levels, were unexpectedly identical in T1D patients and controls. Gene expression profiling on PBMC RNA showed that DMB mRNA was significantly elevated in T1D patients with residual C-peptide. This resulted in higher levels of DM protein in B cells and dendritic cells. DO levels were also increased, suggesting that the MHCII pathway maybe differentially regulated in individuals with residual C-peptide. Collectively, these studies show a dysregulation of the MHCII Ag-processing pathway in patients with T1D.

Relationships between phenotypic plasticity and epigenetic variation in two Caribbean Acropora corals.

The plastic ability for a range of phenotypes to be exhibited by the same genotype allows organisms to respond to environmental variation and may modulate fitness in novel environments. Differing capacities for phenotypic plasticity within a population, apparent as genotype by environment interactions (GxE), can therefore have both ecological and evolutionary implications. Epigenetic gene regulation alters gene function in response to environmental cues without changes to the underlying genetic sequence and likely mediates phenotypic variation. DNA methylation is currently the most well described epigenetic mechanism and is related to transcriptional homeostasis in invertebrates. However, evidence quantitatively linking variation in DNA methylation with that of phenotype is lacking in some taxa, including reef-building corals. In this study, spatial and seasonal environmental variation in Bonaire, Caribbean Netherlands was utilized to assess relationships between physiology and DNA methylation profiles within genetic clones across different genotypes of Acropora cervicornis and A. palmata corals. The physiology of both species was highly influenced by environmental variation compared to the effect of genotype. GxE effects on phenotype were only apparent in A. cervicornis. DNA methylation in both species differed between genotypes and seasons and epigenetic variation was significantly related to coral physiological metrics. Furthermore, plastic shifts in physiology across seasons were significantly positively correlated with shifts in DNA methylation profiles in both species. These results highlight the dynamic influence of environmental conditions and genetic constraints on the physiology of two important Caribbean coral species. Additionally, this study provides quantitative support for the role of epigenetic DNA methylation in mediating phenotypic plasticity in invertebrates.

Effect of Genetic Factors, Age and Sex on Levels of Circulating Extracellular Vesicles and Platelets.

Extracellular vesicles (EVs) mediate cell interactions in biological processes, such as receptor activation or molecule transfer. Estimates of variation by age and sex have been limited by small sample size, and no report has assessed the contribution of genetic factors to levels of EVs. Here, we evaluated blood levels of 25 EV and 3 platelet traits in 974 individuals (933 genotyped) and reported the first genome-wide association study (GWAS) on levels of these traits. EV levels all decreased with age, whereas the trend for their surface markers was more heterogeneous. Platelets and CD31dim platelet EVs significantly increased in females compared to males, although CD31 expression on both platelets and platelet EVs decreased in females. Levels of the other EV subsets were similar between sexes. GWAS revealed three statistically significant genetic signals associated with EV levels in the F10 and GBP1 genes and in the intergenic region between LRIG1 and KBTBD8. These add to a signal in the 3'UTR of RHOF associated with CD31 expression on platelets that was previously found to be associated with other platelet traits. These findings suggest that EV formation is not a simple, constant adjunct of metabolism but is under both age-related and genetic control that can be independent of the regulation of the levels of the cells from which the EVs derive.

Human Hepatitis B Viral Infection Outcomes Are Linked to Naturally Occurring Variants of HLA-DOA That Have Altered Function.

HLA molecules of the MHC class II (MHCII) bind and present pathogen-derived peptides for CD4 T cell activation. Peptide loading of MHCII in the endosomes of cells is controlled by the interplay of the nonclassical MHCII molecules, HLA-DM (DM) and HLA-DO (DO). DM catalyzes peptide loading, whereas DO, an MHCII substrate mimic, prevents DM from interacting with MHCII, resulting in an altered MHCII-peptide repertoire and increased MHCII-CLIP. Although the two genes encoding DO (DOA and DOB) are considered nonpolymorphic, there are rare natural variants. Our previous work identified DOB variants that altered DO function. In this study, we show that natural variation in the DOA gene also impacts DO function. Using the 1000 Genomes Project database, we show that ∼98% of individuals express the canonical DOA*0101 allele, and the remaining individuals mostly express DOA*0102, which we found was a gain-of-function allele. Analysis of 25 natural occurring DOα variants, which included the common alleles, identified three null variants and one variant with reduced and nine with increased ability to modulate DM activity. Unexpectedly, several of the variants produced reduced DO protein levels yet efficiently inhibited DM activity. Finally, analysis of associated single-nucleotide polymorphisms genetically linked the DOA*0102 common allele, a gain-of-function variant, with human hepatitis B viral persistence. In contrast, we found that the DOα F114L null allele was linked with viral clearance. Collectively, these studies show that natural variation occurring in the human DOA gene impacts DO function and can be linked to specific outcomes of viral infections.

PRF1 mutation alters immune system activation, inflammation, and risk of autoimmunity.

BACKGROUND: Defective alleles within the PRF1 gene, encoding the pore-forming protein perforin, in combination with environmental factors, cause familial type 2 hemophagocytic lymphohistiocytosis (FHL2), a rare, severe autosomal recessive childhood disorder characterized by massive release of cytokines-cytokine storm. OBJECTIVE: The aim of this study was to determine the function of hypomorph PRF1:p.A91V g.72360387 G > A on multiple sclerosis (MS) and type 1 diabetes (T1D). METHODS: We cross-compare the association data for PRF1:p.A91V mutation derived from GWAS on adult MS and pediatric T1D in Sardinians. The novel association with T1D was replicated in metanalysis in 12,584 cases and 17,692 controls from Sardinia, the United Kingdom, and Scotland. To dissect this mutation function, we searched through the coincident association immunophenotypes in additional set of general population Sardinians. RESULTS: We report that PRF1:p.A91V, is associated with increase of lymphocyte levels, especially within the cytotoxic memory T-cells, at general population level with reduced interleukin 7 receptor expression on these cells. The minor allele increased risk of MS, in 2903 cases and 2880 controls from Sardinia p = 2.06 × 10-4, odds ratio OR = 1.29, replicating a previous finding, whereas it protects from T1D p = 1.04 × 10-5, OR = 0.82. CONCLUSION: Our results indicate opposing contributions of the cytotoxic T-cell compartment to MS and T1D pathogenesis.

Overexpression of the Cytokine BAFF and Autoimmunity Risk.

BACKGROUND: Genomewide association studies of autoimmune diseases have mapped hundreds of susceptibility regions in the genome. However, only for a few association signals has the causal gene been identified, and for even fewer have the causal variant and underlying mechanism been defined. Coincident associations of DNA variants affecting both the risk of autoimmune disease and quantitative immune variables provide an informative route to explore disease mechanisms and drug-targetable pathways. METHODS: Using case-control samples from Sardinia, Italy, we performed a genomewide association study in multiple sclerosis followed by TNFSF13B locus-specific association testing in systemic lupus erythematosus (SLE). Extensive phenotyping of quantitative immune variables, sequence-based fine mapping, cross-population and cross-phenotype analyses, and gene-expression studies were used to identify the causal variant and elucidate its mechanism of action. Signatures of positive selection were also investigated. RESULTS: A variant in TNFSF13B, encoding the cytokine and drug target B-cell activating factor (BAFF), was associated with multiple sclerosis as well as SLE. The disease-risk allele was also associated with up-regulated humoral immunity through increased levels of soluble BAFF, B lymphocytes, and immunoglobulins. The causal variant was identified: an insertion-deletion variant, GCTGT→A (in which A is the risk allele), yielded a shorter transcript that escaped microRNA inhibition and increased production of soluble BAFF, which in turn up-regulated humoral immunity. Population genetic signatures indicated that this autoimmunity variant has been evolutionarily advantageous, most likely by augmenting resistance to malaria. CONCLUSIONS: A TNFSF13B variant was associated with multiple sclerosis and SLE, and its effects were clarified at the population, cellular, and molecular levels. (Funded by the Italian Foundation for Multiple Sclerosis and others.).

A genome-wide association scan on the levels of markers of inflammation in Sardinians reveals associations that underpin its complex regulation.

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals-5 of which were identified only with the custom arrays-and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10(-29)); for ESR, at the HBB (rs4910472, p = 2.31×10(-11)) and UCN119B/SPPL3 (rs11829037, p = 8.91×10(-10)) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10(-13)) and in CADM3 (rs3026968, p = 7.63×10(-13)); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10(-21)), near DARC (rs3845624, p = 1.43×10(-10)), UNC119B/SPPL3 (rs11829037, p = 1.50×10(-14)), and ICOSLG/AIRE (rs113459440, p = 1.54×10(-08)) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process.

Cross-sectional analysis of the humoral response after SARS-CoV-2 vaccination in Sardinian multiple sclerosis patients, a follow-up study.

Monitoring immune responses to SARS-CoV-2 vaccination and its clinical efficacy over time in Multiple Sclerosis (MS) patients treated with disease-modifying therapies (DMTs) help to establish the optimal strategies to ensure adequate COVID-19 protection without compromising disease control offered by DMTs. Following our previous observations on the humoral response one month after two doses of BNT162b2 vaccine (T1) in MS patients differently treated, here we present a cross-sectional and longitudinal follow-up analysis six months following vaccination (T2, n=662) and one month following the first booster (T3, n=185). Consistent with results at T1, humoral responses were decreased in MS patients treated with fingolimod and anti-CD20 therapies compared with untreated patients also at the time points considered here (T2 and T3). Interestingly, a strong upregulation one month after the booster was observed in patients under every DMTs analyzed, including those treated with fingolimod and anti-CD20 therapies. Although patients taking these latter therapies had a higher rate of COVID-19 infection five months after the first booster, only mild symptoms that did not require hospitalization were reported for all the DMTs analyzed here. Based on these findings we anticipate that additional vaccine booster shots will likely further improve immune responses and COVID-19 protection in MS patients treated with any DMT.

Effect of Different Disease-Modifying Therapies on Humoral Response to BNT162b2 Vaccine in Sardinian Multiple Sclerosis Patients.

Objectives: Vaccination against COVID-19 is highly recommended to patients affected by multiple sclerosis (MS); however, the impact of MS disease-modifying therapies (DMTs) on the immune response following vaccination has been only partially investigated. Here, we aimed to elucidate the effect of DMTs on the humoral immune response to mRNA-based anti-SARS-CoV-2 vaccines in MS patients. Methods: We obtained sera from 912 Sardinian MS patients and 63 healthy controls 30 days after the second dose of BNT162b2 vaccine and tested them for SARS-CoV-2 response using anti-Spike (S) protein-based serology. Previous SARS-CoV-2 infection was assessed by anti-Nucleocapsid (N) serology. Patients were either untreated or undergoing treatment with a total of 13 different DMTs. Differences between treatment groups comprised of at least 10 patients were assessed by generalized linear mixed-effects model. Demographic and clinical data and smoking status were analyzed as additional factors potentially influencing humoral immunity from COVID-19 vaccine. Results: MS patients treated with natalizumab, teriflunomide, azathioprine, fingolimod, ocrelizumab, and rituximab showed significantly lower humoral responses compared to untreated patients. We did not observe a statistically significant difference in response between patients treated with the other drugs (dimethyl fumarate, interferon, alemtuzumab and glatiramer acetate) and untreated patients. In addition, older age, male sex and active smoking were significantly associated with lower antibody titers against SARS-CoV-2. MS patients previously infected with SARS-CoV-2 had significantly higher humoral responses to vaccine than uninfected patients. Conclusion: Humoral response to BNT162b2 is significantly influenced by the specific DMTs followed by patients, as well as by other factors such as previous SARS-CoV-2 infection, age, sex, and smoking status. These results are important to inform targeted strategies to prevent clinically relevant COVID-19 in MS patients.

Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.

We report on the influence of ~22 million variants on 731 immune cell traits in a cohort of 3,757 Sardinians. We detected 122 significant (P < 1.28 × 10-11) independent association signals for 459 cell traits at 70 loci (53 of them novel) identifying several molecules and mechanisms involved in cell regulation. Furthermore, 53 signals at 36 loci overlapped with previously reported disease-associated signals, predominantly for autoimmune disorders, highlighting intermediate phenotypes in pathogenesis. Collectively, our findings illustrate complex genetic regulation of immune cells with highly selective effects on autoimmune disease risk at the cell-subtype level. These results identify drug-targetable pathways informing the design of more specific treatments for autoimmune diseases.

Author Correction: Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.