Genetic mapping across autoimmune diseases reveals shared associations and mechanisms

Autoimmune and inflammatory diseases are polygenic disorders of the immune system. Many genomic loci harbor risk alleles for several diseases, but the limited resolution of genetic mapping prevents determining whether the same allele is responsible, indicating a shared underlying mechanism. Here, using a collection of 129,058 cases and controls across 6 diseases, we show that ~40% of overlapping associations are due to the same allele. We improve fine-mapping resolution for shared alleles twofold by combining cases and controls across diseases, allowing us to identify more expression quantitative trait loci driven by the shared alleles. The patterns indicate widespread sharing of pathogenic mechanisms but not a single global autoimmune mechanism. Our approach can be applied to any set of traits and is particularly valuable as sample collections become depleted.

A review of genetic risk in systemic lupus erythematosus

INTRODUCTION

Systemic Lupus Erythematosus (SLE) is a complex multisystem autoimmune disease with a wide range of signs and symptoms in affected individuals. The utilization of genome-wide association study (GWAS) technology has led to an explosion in the number of genetic risk factors mapped for autoimmune diseases, including SLE.

AREAS COVERED

In this review, we summarize the more recent genetic risk loci mapped in SLE, which bring the total number of loci mapped to approximately 200. We review prioritization analyses of the associated variants and experimental validation of the putative causal variants. This includes the implementation of new bioinformatic techniques to align genomic and functional data and the use of transcriptomics with single-cell RNA-sequencing, CRISPR genome editing, and Massive Parallel Reporter Assays to analyze non-coding regulatory genetics.

EXPERT OPINION

Despite progress in identifying more genetic risk loci and variant-gene pairs for SLE, understanding its pathogenesis and applying findings clinically remains challenging. The polygenic risk score (PRS) has been used as an application of SLE genetics, but with limited performance in non-EUR populations. In the next few years, advancements in proteomics, post-translational modification estimation, and whole-genome sequencing will enhance disease understanding.

Genetics of SLE: does this explain susceptibility and severity across racial groups?

The prevalence and severity of systemic lupus erythematosus (SLE) have been found to vary across populations of different ancestries. This review explores whether these differences can be explained by the genetic aetiology of the condition. Large genetic studies suggest that populations of different ancestry share the same risk loci but the individual risk alleles are more common in some, leading to a higher prevalence, severity, and earlier onset of the condition. Despite many of the loci being shared across populations, some have been found to be ancestry-specific and these are hypothesised to have undergone differential selective pressure in recent human history. Additionally, the effectiveness of some of the drugs used in SLE has been found to vary across ancestries, which might affect the progression of the disease, but it is unclear whether these differences are pharmacogenetic. We concluded that to understand the full role of genetics in the risk, presentation and response to treatment of SLE, larger studies including individuals from a wider representation of ancestries will be required.

Comprehensive genetic and functional analyses of Fc gamma receptors influence on response to rituximab therapy for autoimmunity

BACKGROUND

Rituximab is widely used to treat autoimmunity but clinical response varies. Efficacy is determined by the efficiency of B-cell depletion, which may depend on various Fc gamma receptor (FcγR)-dependent mechanisms. Study of FcγR is challenging due to the complexity of the FCGR genetic locus. We sought to assess the effect of FCGR variants on clinical response, B-cell depletion and NK-cell-mediated killing in rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE).

METHODS

A longitudinal cohort study was conducted in 835 patients [RA = 573; SLE = 262]. Clinical outcome measures were two-component disease activity score in 28-joints (2C-DAS28CRP) for RA and British Isles Lupus Assessment Group (BILAG)-2004 major clinical response (MCR) for SLE at 6 months. B-cells were evaluated by highly-sensitive flow cytometry. Single nucleotide polymorphism and copy number variation for genes encoding five FcγRs were measured using multiplex ligation-dependent probe amplification. Ex vivo studies assessed NK-cell antibody-dependent cellular cytotoxicity (ADCC) and FcγR expression.

FINDINGS

In RA, carriage of FCGR3A-158V and increased FCGR3A-158V copies were associated with greater 2C-DAS28CRP response (adjusted for baseline 2C-DAS28CRP). In SLE, MCR was associated with increased FCGR3A-158V, OR 1.64 (95% CI 1.12-2.41) and FCGR2C-ORF OR 1.93 (95% CI 1.09-3.40) copies. 236/413 (57%) patients with B-cell data achieved complete depletion. Homozygosity for FCGR3A-158V and increased FCGR3A-158V copies were associated with complete depletion in combined analyses. FCGR3A genotype was associated with rituximab-induced ADCC, and increased NK-cell FcγRIIIa expression was associated with improved clinical response and depletion in vivo. Furthermore, disease status and concomitant therapies impacted both NK-cell FcγRIIIa expression and ADCC.

INTERPRETATION

FcγRIIIa is the major low affinity FcγR associated with rituximab response. Increased copies of the FCGR3A-158V allele (higher affinity for IgG1), influences clinical and biological responses to rituximab in autoimmunity. Enhancing FcγR-effector functions could improve the next generation of CD20-depleting therapies and genotyping may stratify patients for optimal treatment protocols.

FUNDING

Medical Research Council, National Institute for Health and Care Research, Versus Arthritis.

Age acquired skewed X chromosome inactivation is associated with adverse health outcomes in humans

BACKGROUND

Ageing is a heterogenous process characterised by cellular and molecular hallmarks, including changes to haematopoietic stem cells and is a primary risk factor for chronic diseases. X chromosome inactivation (XCI) randomly transcriptionally silences either the maternal or paternal X in each cell of 46, XX females to balance the gene expression with 46, XY males. Age acquired XCI-skew describes the preferential selection of cells across a tissue resulting in an imbalance of XCI, which is particularly prevalent in blood tissues of ageing females, and yet its clinical consequences are unknown.

METHODS

We assayed XCI in 1575 females from the TwinsUK population cohort using DNA extracted from whole blood. We employed prospective, cross-sectional, and intra-twin study designs to characterise the relationship of XCI-skew with molecular and cellular measures of ageing, cardiovascular disease risk, and cancer diagnosis.

RESULTS

We demonstrate that XCI-skew is independent of traditional markers of biological ageing and is associated with a haematopoietic bias towards the myeloid lineage. Using an atherosclerotic cardiovascular disease risk score, which captures traditional risk factors, XCI-skew is associated with an increased cardiovascular disease risk both cross-sectionally and within XCI-skew discordant twin pairs. In a prospective 10 year follow-up study, XCI-skew is predictive of future cancer incidence.

CONCLUSIONS

Our study demonstrates that age acquired XCI-skew captures changes to the haematopoietic stem cell population and has clinical potential as a unique biomarker of chronic disease risk.

FUNDING

KSS acknowledges funding from the Medical Research Council [MR/M004422/1 and MR/R023131/1]. JTB acknowledges funding from the ESRC [ES/N000404/1]. MM acknowledges funding from the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. TwinsUK is funded by the Wellcome Trust, Medical Research Council, European Union, Chronic Disease Research Foundation (CDRF), Zoe Global Ltd and the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London.

COVID-19 and systemic lupus erythematosus genetics: A balance between autoimmune disease risk and protection against infection

Genome wide association studies show there is a genetic component to severe COVID-19. We find evidence that the genome-wide genetic association signal with severe COVID-19 is correlated with that of systemic lupus erythematosus (SLE), having formally tested this using genetic correlation analysis by LD score regression. To identify the shared associated loci and gain insight into the shared genetic effects, using summary level data we performed meta-analyses, a local genetic correlation analysis and fine-mapping using stepwise regression and functional annotation. This identified multiple loci shared between the two traits, some of which exert opposing effects. The locus with most evidence of shared association is TYK2, a gene critical to the type I interferon pathway, where the local genetic correlation is negative. Another shared locus is CLEC1A, where the direction of effects is aligned, that encodes a lectin involved in cell signaling, and the anti-fungal immune response. Our analyses suggest that several loci with reciprocal effects between the two traits have a role in the defense response pathway, adding to the evidence that SLE risk alleles are protective against infection.

We observed a correlation between the genetic associations with severe COVID-19 and those with systemic lupus erythematosus (SLE, Lupus), and aimed to discover which genetic loci were shared by these diseases and what biological processes were involved. This resulted in the discovery of several genetic loci, some of which had alleles that were risk for both diseases and some of which were risk for severe COVID-19 yet protective for SLE. The locus with most evidence of shared association (TYK2) is involved in interferon production, a process that is important in response to viral infection and known to be dysregulated in SLE patients. Other shared associated loci contained genes also involved in the defense response and the immune system signaling. These results add to the growing evidence that there are alleles in the human genome that provide protection against viral infection yet are risk for autoimmune disease.

Immunodeficiency, autoimmunity, and increased risk of B cell malignancy in humans with TRAF3 mutations

Tumor necrosis factor receptor-associated factor 3 (TRAF3) is a central regulator of immunity. TRAF3 is often somatically mutated in B cell malignancies, but its role in human immunity is not defined. Here, in five unrelated families, we describe an immune dysregulation syndrome of recurrent bacterial infections, autoimmunity, systemic inflammation, B cell lymphoproliferation, and hypergammaglobulinemia. Affected individuals each had monoallelic mutations in TRAF3 that reduced TRAF3 expression. Immunophenotyping showed that patients' B cells were dysregulated, exhibiting increased nuclear factor-κB 2 activation, elevated mitochondrial respiration, and heightened inflammatory responses. Patients had mild CD4+ T cell lymphopenia, with a reduced proportion of naïve T cells but increased regulatory T cells and circulating T follicular helper cells. Guided by this clinical phenotype, targeted analyses demonstrated that common genetic variants, which also reduce TRAF3 expression, are associated with an increased risk of B cell malignancies, systemic lupus erythematosus, higher immunoglobulin levels, and bacterial infections in the wider population. Reduced TRAF3 conveys disease risks by driving B cell hyperactivity via intrinsic activation of multiple intracellular proinflammatory pathways and increased mitochondrial respiration, with a likely contribution from dysregulated T cell help. Thus, we define monogenic TRAF3 haploinsufficiency syndrome and demonstrate how common TRAF3 variants affect a range of human diseases.

OA13 Comprehensive genetic and functional analyses of Fc gamma receptors explain response to rituximab therapy for autoimmune rheumatic diseases

Background/Aims

Rituximab is widely used to treat rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) but clinical response varies. Efficacy is determined by the efficiency of depletion, which may depend on a variety of Fc gamma receptor (FcγR)-dependent mechanisms. Previous research was limited by complexity of the FCGR locus, not integrating copy number variation with functional SNP, and small sample size. Here, we aimed to assess the effect of the full range of FcγRs variants on depletion, clinical response and functional effect on NK-cell-mediated killing in two rheumatic diseases with a view to personalised B-cell depleting therapies.

Methods

A prospective longitudinal study was conducted in 873 patients (RA = 611; SLE=262) from four cohorts (BSRBR-RA, BILAG-BR, Leeds RA and Leeds SLE Biologics). For RA, the outcome measures were 3C-DAS28CRP and 2C-DAS28CRP at 6 (+/-3) months post-rituximab (adjusted for baseline DAS28). For SLE, major clinical response (MCR) was defined as improvement of active BILAG-2004 domains to grade C/better at 6 months. B-cell depletion was evaluated by highly-sensitive flow cytometry. Qualitative and quantitative polymorphisms for five major FcγRs were measured using a commercial multiplex ligation-dependent probe amplification. Median NK cell FcγRIIIa expression (CD3-CD56+CD16+) and NK-cell degranulation (CD107a) in the presence of rituximab-coated Daudi/Raji B-cell lines were assessed using flow cytometry.

Results

In RA, for FCGR3A, carriage of V allele (coefficient -0.25 [SE 0.11]; p = 0.02) and increased copies of V allele (-0.20 [0.09]; p = 0.02) were associated with better 2C-DAS28 response. Irrespective of FCGR3A genotype, increased gene copies were associated with a better response. In SLE, 177/262 (67.6%) achieved BILAG response (MCR=34.4%; Partial=33.2%). MCR was associated with increased copies of FCGR3A-158V allele, OR 1.64 (95% CI 1.12-2.41) and FCGR2C-ORF allele 1.93 (1.09-3.40). Of patients with B-cells data in the combined cohort, 236/413 (57%) achieved complete depletion post-rituximab. Only homozygosity for V allele and higher copies of FCGR3A V allele were associated with increased odds of depletion. Patients with complete depletion had higher NK cell FcγRIIIa expression at rituximab initiation than those with incomplete depletion (p = 0.04) and this higher expression was associated with improved EULAR response in RA. Moreover, for FCGR3A, degranulation activity was increased in V allele carriers vs FF genotype in the combined cohort; p = 0.02.

Conclusion

FcγRIIIa is the major low affinity FcγR and increased copies of the FCGR3A-158V allele, encoding the allotype with a higher affinity for IgG1, was associated with clinical and biological responses to rituximab in two autoimmune diseases. This was supported by functional data on NK cell-mediated cytotoxicity. In SLE, increased copies of the FCGR2C-ORF allele was also associated with improved response. Our findings indicate that enhancing FcγR-effector functions could improve the next generation of CD20-depleting therapies and genotyping could stratify patients for optimal treatment protocols.

Disclosure

M. Md Yusof: None. J. Robinson: None. V. Davies: None. D. Wild: None. M. Morgan: None. J. Taylor: None. Y. El-Sherbiny: None. D. Morris: None. L. Liu: None. A. Rawstron: None. M. Buch: None. D. Plant: None. H. Cordell: None. J. Isaacs: None. I. Bruce: None. P. Emery: Grants/research support; PE has received consultancy fees and funding for research from Roche within the last 3 years. A. Barton: None. T. Vyse: None. J. Barrett: None. E. Vital: Grants/research support; EMV has received consultancy fees and funding for research from Roche within the last 3 years. A. Morgan: None.

Genome-wide assessment of genetic risk for systemic lupus erythematosus and disease severity

Using three European and two Chinese genome-wide association studies (GWAS), we investigated the performance of genetic risk scores (GRSs) for predicting the susceptibility and severity of systemic lupus erythematosus (SLE), using renal disease as a proxy for severity. We used four GWASs to test the performance of GRS both cross validating within the European population and between European and Chinese populations. The performance of GRS in SLE risk prediction was evaluated by receiver operating characteristic (ROC) curves. We then analyzed the polygenic nature of SLE statistically. We also partitioned patients according to their age-of-onset and evaluated the predictability of GRS in disease severity in each age group. We found consistently that the best GRS in the prediction of SLE used SNPs associated at the level of P < 1e−05 in all GWAS data sets and that SNPs with P-values above 0.2 were inflated for SLE true positive signals. The GRS results in an area under the ROC curve ranging between 0.64 and 0.72, within European and between the European and Chinese populations. We further showed a significant positive correlation between a GRS and renal disease in two independent European GWAS (P_cohort1 = 2.44e−08; P_cohort2 = 0.00205) and a significant negative correlation with age of SLE onset (P_cohort1 = 1.76e−12; P_cohort2 = 0.00384). We found that the GRS performed better in the prediction of renal disease in the ‘later onset’ compared with the ‘earlier onset’ group. The GRS predicts SLE in both European and Chinese populations and correlates with poorer prognostic factors: young age-of-onset and lupus nephritis.

Identification of 38 novel loci for systemic lupus erythematosus and genetic heterogeneity between ancestral groups

Systemic lupus erythematosus (SLE), a worldwide autoimmune disease with high heritability, shows differences in prevalence, severity and age of onset among different ancestral groups. Previous genetic studies have focused more on European populations, which appear to be the least affected. Consequently, the genetic variations that underlie the commonalities, differences and treatment options in SLE among ancestral groups have not been well elucidated. To address this, we undertake a genome-wide association study, increasing the sample size of Chinese populations to the level of existing European studies. Thirty-eight novel SLE-associated loci and incomplete sharing of genetic architecture are identified. In addition to the human leukocyte antigen (HLA) region, nine disease loci show clear ancestral differences and implicate antibody production as a potential mechanism for differences in disease manifestation. Polygenic risk scores perform significantly better when trained on ancestry-matched data sets. These analyses help to reveal the genetic basis for disparities in SLE among ancestral groups.

Identification of susceptibility loci for Takayasu arteritis through a large multi-ancestral genome-wide association study

Takayasu arteritis is a rare inflammatory disease of large arteries. We performed a genetic study in Takayasu arteritis comprising 6,670 individuals (1,226 affected individuals) from five different populations. We discovered HLA risk factors and four non-HLA susceptibility loci in VPS8, SVEP1, CFL2, and chr13q21 and reinforced IL12B, PTK2B, and chr21q22 as robust susceptibility loci shared across ancestries. Functional analysis proposed plausible underlying disease mechanisms and pinpointed ETS2 as a potential causal gene for chr21q22 association. We also identified >60 candidate loci with suggestive association (p < 5 × 10^-5) and devised a genetic risk score for Takayasu arteritis. Takayasu arteritis was compared to hundreds of other traits, revealing the closest genetic relatedness to inflammatory bowel disease. Epigenetic patterns within risk loci suggest roles for monocytes and B cells in Takayasu arteritis. This work enhances understanding of the genetic basis and pathophysiology of Takayasu arteritis and provides clues for potential new therapeutic targets.

Nucleolin acts as the receptor for C1QTNF4 and supports C1QTNF4-mediated innate immunity modulation

The C1q and TNF related 4 (C1QTNF4) protein is a structurally unique member of the C1QTNF family, a family of secreted proteins that have structural homology with both complement C1q and the tumor necrosis factor superfamily. C1QTNF4 has been linked to the autoimmune disease systemic lupus erythematosus through genetic studies; however, its role in immunity and inflammation remains poorly defined and a cell surface receptor of C1QTNF4 has yet to be identified. Here we report identification of nucleolin as a cell surface receptor of C1QTNF4 using mass spectrometric analysis. Additionally, we present evidence that the interaction between C1QTNF4 and nucleolin is mediated by the second C1q-like domain of C1QTNF4 and the C terminus of nucleolin. We show that monocytes and B cells are target cells of C1QTNF4 and observe extensive binding to dead cells. Imaging flow cytometry experiments in monocytes show that C1QTNF4 becomes actively internalized upon cell binding. Our results suggest that nucleolin may serve as a docking molecule for C1QTNF4 and act in a context-dependent manner through coreceptors. Taken together, these findings further our understanding of C1QTNF4's function in the healthy immune system and how dysfunction may contribute to the development of systemic lupus erythematosus.

Trans-Ancestral Fine-Mapping and Epigenetic Annotation as Tools to Delineate Functionally Relevant Risk Alleles at IKZF1 and IKZF3 in Systemic Lupus Erythematosus

Background: Prioritizing tag-SNPs carried on extended risk haplotypes at susceptibility loci for common disease is a challenge. Methods: We utilized trans-ancestral exclusion mapping to reduce risk haplotypes at IKZF1 and IKZF3 identified in multiple ancestries from SLE GWAS and ImmunoChip datasets. We characterized functional annotation data across each risk haplotype from publicly available datasets including ENCODE, RoadMap Consortium, PC Hi-C data from 3D genome browser, NESDR NTR conditional eQTL database, GeneCards Genehancers and TF (transcription factor) binding sites from Haploregv4. Results: We refined the 60 kb associated haplotype upstream of IKZF1 to just 12 tag-SNPs tagging a 47.7 kb core risk haplotype. There was preferential enrichment of DNAse I hypersensitivity and H3K27ac modification across the 3′ end of the risk haplotype, with four tag-SNPs sharing allele-specific TF binding sites with promoter variants, which are eQTLs for IKZF1 in whole blood. At IKZF3, we refined a core risk haplotype of 101 kb (27 tag-SNPs) from an initial extended haplotype of 194 kb (282 tag-SNPs), which had widespread DNAse I hypersensitivity, H3K27ac modification and multiple allele-specific TF binding sites. Dimerization of Fox family TFs bound at the 3′ and promoter of IKZF3 may stabilize chromatin looping across the locus. Conclusions: We combined trans-ancestral exclusion mapping and epigenetic annotation to identify variants at both IKZF1 and IKZF3 with the highest likelihood of biological relevance. The approach will be of strong interest to other complex trait geneticists seeking to attribute biological relevance to risk alleles on extended risk haplotypes in their disease of interest.

Independent Replication on Genome-Wide Association Study Signals Identifies IRF3 as a Novel Locus for Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is a genetically complex autoimmune disease. Despite the significant progress made in identifying susceptibility genes for SLE, the genetic architecture of the disease is far from being understood. In this study, we set to replicate a number of suggestive association signals found in genome-wide association studies (GWASs) in additional independent cohorts. Replication studies were performed on Han Chinese cohorts from Hong Kong and Anhui, involving a total of 2,269 cases and 5,073 controls. We identified a missense variant in IRF3 (rs7251) reaching genome-wide significance through a joint analysis of GWAS and replication data (OR = 0.876, P = 4.40E-08). A significant correlation was observed between rs7251 and lupus nephritis (LN) by subphenotype stratification (OR = 0.785, P = 0.0128). IRF3 is a key molecule in type I interferon production upon nucleic acid antigen stimulations and may inhibit regulatory T cell differentiation. Further elucidation of the mechanism of this association could help us better understand the pathogenesis of SLE.

Complement genes contribute sex-biased vulnerability in diverse disorders

Many common illnesses differentially affect men and women for unknown reasons. The autoimmune diseases lupus and Sjögren’s syndrome affect nine times more women than men¹, whereas schizophrenia affects men more frequently and severely². All three illnesses have their strongest common genetic associations in the Major Histocompatibility Complex (MHC) locus, an association that in lupus and Sjögren’s syndrome has long been thought to arise from alleles of the human leukocyte antigen (HLA) genes at that locus^3–6. Here we show that the complement component 4 (C4) genes, which are also in the MHC locus and were recently found to increase risk for schizophrenia⁷, generate 7-fold variation in risk for lupus (95% CI: 5.88–8.61; p < 10⁻¹¹⁷ in total) and 16-fold variation in risk for Sjögren’s syndrome (95% CI: 8.59–30.89; p < 10⁻²³ in total) among individuals with common C4 genotypes, with C4A protecting more strongly than C4B in both illnesses. The same alleles that increase risk for schizophrenia greatly reduced risk for lupus and Sjögren’s syndrome. In all three illnesses, C4 alleles acted more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for lupus, 31-fold variation in risk for Sjögren’s syndrome, and 1.7-fold variation in schizophrenia risk among men (vs. 6-fold, 15-fold, and 1.26-fold among women respectively). At a protein level, both C4 and its effector C3 were present at greater levels in men than women in cerebrospinal fluid (p < 10⁻⁵ for both C4 and C3) and plasma^8,9 among adults ages 20–50, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help explain the larger effects of C4 alleles in men, women’s greater risk of SLE and Sjögren’s, and men’s greater vulnerability in schizophrenia. These results implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

Genetic overlap between autoimmune diseases and non-Hodgkin lymphoma subtypes

Epidemiologic studies show an increased risk of non-Hodgkin lymphoma (NHL) in patients with autoimmune disease (AD), due to a combination of shared environmental factors and/or genetic factors, or a causative cascade: chronic inflammation/antigen-stimulation in one disease leads to another. Here we assess shared genetic risk in genome-wide-association-studies (GWAS). Secondary analysis of GWAS of NHL subtypes (chronic lymphocytic leukemia, diffuse large B-cell lymphoma, follicular lymphoma, and marginal zone lymphoma) and ADs (rheumatoid arthritis, systemic lupus erythematosus, and multiple sclerosis). Shared genetic risk was assessed by (a) description of regional genetic of overlap, (b) polygenic risk score (PRS), (c)"diseasome", (d)meta-analysis. Descriptive analysis revealed few shared genetic factors between each AD and each NHL subtype. The PRS of ADs were not increased in NHL patients (nor vice versa). In the diseasome, NHLs shared more genetic etiology with ADs than solid cancers (p = .0041). A meta-analysis (combing AD with NHL) implicated genes of apoptosis and telomere length. This GWAS-based analysis four NHL subtypes and three ADs revealed few weakly-associated shared loci, explaining little total risk. This suggests common genetic variation, as assessed by GWAS in these sample sizes, may not be the primary explanation for the link between these ADs and NHLs.

Interferon inducible X-linked gene CXorf21 may contribute to sexual dimorphism in Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is an autoimmune disease, characterised by increased expression of type I interferon (IFN)-regulated genes and a striking sex imbalance towards females. Through combined genetic, in silico, in vitro, and ex vivo approaches, we define CXorf21, a gene of hitherto unknown function, which escapes X-chromosome inactivation, as a candidate underlying the Xp21.2 SLE association. We demonstrate that CXorf21 is an IFN-response gene and that the sexual dimorphism in expression is magnified by immunological challenge. Fine-mapping reveals a single haplotype as a potential causal cis-eQTL for CXorf21. We propose that expression is amplified through modification of promoter and 3'-UTR chromatin interactions. Finally, we show that the CXORF21 protein colocalises with TLR7, a pathway implicated in SLE pathogenesis. Our study reveals modulation in gene expression affected by the combination of two hallmarks of SLE: CXorf21 expression increases in a both an IFN-inducible and sex-specific manner.

Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans

Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A-C-B-DRB1-DQA-DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA- B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA-C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

Tissue-Restricted Adaptive Type 2 Immunity Is Orchestrated by Expression of the Costimulatory Molecule OX40L on Group 2 Innate Lymphoid Cells

The local regulation of type 2 immunity relies on dialog between the epithelium and the innate and adaptive immune cells. Here we found that alarmin-induced expression of the co-stimulatory molecule OX40L on group 2 innate lymphoid cells (ILC2s) provided tissue-restricted T cell co-stimulation that was indispensable for Th2 and regulatory T (Treg) cell responses in the lung and adipose tissue. Interleukin (IL)-33 administration resulted in organ-specific surface expression of OX40L on ILC2s and the concomitant expansion of Th2 and Treg cells, which was abolished upon deletion of OX40L on ILC2s (Il7raCre/+Tnfsf4fl/fl mice). Moreover, Il7raCre/+Tnfsf4fl/fl mice failed to mount effective Th2 and Treg cell responses and corresponding adaptive type 2 pulmonary inflammation arising from Nippostrongylus brasiliensis infection or allergen exposure. Thus, the increased expression of OX40L in response to IL-33 acts as a licensing signal in the orchestration of tissue-specific adaptive type 2 immunity, without which this response fails to establish.

Meta-analysis of GWAS on both Chinese and European populations identifies GPR173 as a novel X chromosome susceptibility gene for SLE

BACKGROUND

Systemic lupus erythematous (SLE) is a complex autoimmune disease with female predominance, particularly affecting those of childbearing age. We performed analysis of three genome-wide genotyping datasets of populations of both Chinese and European origin.

METHODS

This study involved 5695 cases and 10,357 controls in the discovery stage. The lead signal on chromosome X was followed by replication in three additional Asian cohorts, with 2300 cases and 4244 controls in total. Conditional analysis of the known associated loci on chromosome X was also performed to further explore independent signals.

RESULTS

Single-nucleotide polymorphism rs13440883 in GPR173 was found to be significantly associated with SLE (P_meta = 7.53 × 10^- 9, OR_meta= 1.16), whereas conditional analysis provided evidence of a potential independent signal in the L1CAM-IRAK1-MECP2 region in Asian populations (rs5987175 [LCA10]).

CONCLUSIONS

We identified a novel SLE susceptibility locus on the X chromosome. This finding emphasizes the importance of the X chromosome in disease pathogenesis and highlights the role of sex chromosomes in the female bias of SLE.

Identification of ST3AGL4, MFHAS1, CSNK2A2 and CD226 as loci associated with systemic lupus erythematosus (SLE) and evaluation of SLE genetics in drug repositioning

OBJECTIVES

Systemic lupus erythematosus (SLE) is a prototype autoimmune disease with a strong genetic component in its pathogenesis. Through genome-wide association studies (GWAS), we recently identified 10 novel loci associated with SLE and uncovered a number of suggestive loci requiring further validation. This study aimed to validate those loci in independent cohorts and evaluate the role of SLE genetics in drug repositioning.

METHODS

We conducted GWAS and replication studies involving 12 280 SLE cases and 18 828 controls, and performed fine-mapping analyses to identify likely causal variants within the newly identified loci. We further scanned drug target databases to evaluate the role of SLE genetics in drug repositioning.

RESULTS

We identified three novel loci that surpassed genome-wide significance, including ST3AGL4 (rs13238909, p_meta=4.40E-08), MFHAS1 (rs2428, p_meta=1.17E-08) and CSNK2A2 (rs2731783, p_meta=1.08E-09). We also confirmed the association of CD226 locus with SLE (rs763361, p_meta=2.45E-08). Fine-mapping and functional analyses indicated that the putative causal variants in CSNK2A2 locus reside in an enhancer and are associated with expression of CSNK2A2 in B-lymphocytes, suggesting a potential mechanism of association. In addition, we demonstrated that SLE risk genes were more likely to be interacting proteins with targets of approved SLE drugs (OR=2.41, p=1.50E-03) which supports the role of genetic studies to repurpose drugs approved for other diseases for the treatment of SLE.

CONCLUSION

This study identified three novel loci associated with SLE and demonstrated the role of SLE GWAS findings in drug repositioning.

De novo mutations implicate novel genes in systemic lupus erythematosus

The omnigenic model of complex disease stipulates that the majority of the heritability will be explained by the effects of common variation on genes in the periphery of core disease pathways. Rare variant associations, expected to explain far less of the heritability, may be enriched in core disease genes and thus will be instrumental in the understanding of complex disease pathogenesis and their potential therapeutic targets. Here, using complementary whole-exome sequencing, high-density imputation, and in vitro cellular assays, we identify candidate core genes in the pathogenesis of systemic lupus erythematosus (SLE). Using extreme-phenotype sampling, we sequenced the exomes of 30 SLE parent-affected-offspring trios and identified 14 genes with missense de novo mutations (DNM), none of which are within the >80 SLE susceptibility loci implicated through genome-wide association studies. In a follow-up cohort of 10, 995 individuals of matched European ancestry, we imputed genotype data to the density of the combined UK10K-1000 genomes Phase III reference panel across the 14 candidate genes. Gene-level analyses indicate three functional candidates: DNMT3A, PRKCD, and C1QTNF4. We identify a burden of rare variants across PRKCD associated with SLE risk (P = 0.0028), and across DNMT3A associated with two severe disease prognosis sub-phenotypes (P = 0.0005 and P = 0.0033). We further characterise the TNF-dependent functions of the third candidate gene C1QTNF4 on NF-κB activation and apoptosis, which are inhibited by the p.His198Gln DNM. Our results identify three novel genes in SLE susceptibility and support extreme-phenotype sampling and DNM gene discovery to aid the search for core disease genes implicated through rare variation.

Profiling RNA-Seq at multiple resolutions markedly increases the number of causal eQTLs in autoimmune disease

Genome-wide association studies have identified hundreds of risk loci for autoimmune disease, yet only a minority (~25%) share genetic effects with changes to gene expression (eQTLs) in immune cells. RNA-Seq based quantification at whole-gene resolution, where abundance is estimated by culminating expression of all transcripts or exons of the same gene, is likely to account for this observed lack of colocalisation as subtle isoform switches and expression variation in independent exons can be concealed. We performed integrative cis-eQTL analysis using association statistics from twenty autoimmune diseases (560 independent loci) and RNA-Seq data from 373 individuals of the Geuvadis cohort profiled at gene-, isoform-, exon-, junction-, and intron-level resolution in lymphoblastoid cell lines. After stringently testing for a shared causal variant using both the Joint Likelihood Mapping and Regulatory Trait Concordance frameworks, we found that gene-level quantification significantly underestimated the number of causal cis-eQTLs. Only 5.0-5.3% of loci were found to share a causal cis-eQTL at gene-level compared to 12.9-18.4% at exon-level and 9.6-10.5% at junction-level. More than a fifth of autoimmune loci shared an underlying causal variant in a single cell type by combining all five quantification types; a marked increase over current estimates of steady-state causal cis-eQTLs. Causal cis-eQTLs detected at different quantification types localised to discrete epigenetic annotations. We applied a linear mixed-effects model to distinguish cis-eQTLs modulating all expression elements of a gene from those where the signal is only evident in a subset of elements. Exon-level analysis detected disease-associated cis-eQTLs that subtly altered transcription globally across the target gene. We dissected in detail the genetic associations of systemic lupus erythematosus and functionally annotated the candidate genes. Many of the known and novel genes were concealed at gene-level (e.g. IKZF2, TYK2, LYST). Our findings are provided as a web resource.

X Chromosome Dose and Sex Bias in Autoimmune Diseases: Increased Prevalence of 47,XXX in Systemic Lupus Erythematosus and Sjögren's Syndrome

OBJECTIVE

More than 80% of autoimmune disease predominantly affects females, but the mechanism for this female bias is poorly understood. We suspected that an X chromosome dose effect accounts for this, and we undertook this study to test our hypothesis that trisomy X (47,XXX; occurring in ∼1 in 1,000 live female births) would be increased in patients with female-predominant diseases (systemic lupus erythematosus [SLE], primary Sjögren's syndrome [SS], primary biliary cirrhosis, and rheumatoid arthritis [RA]) compared to patients with diseases without female predominance (sarcoidosis) and compared to controls.

METHODS

All subjects in this study were female. We identified subjects with 47,XXX using aggregate data from single-nucleotide polymorphism arrays, and, when possible, we confirmed the presence of 47,XXX using fluorescence in situ hybridization or quantitative polymerase chain reaction.

RESULTS

We found 47,XXX in 7 of 2,826 SLE patients and in 3 of 1,033 SS patients, but in only 2 of 7,074 controls (odds ratio in the SLE and primary SS groups 8.78 [95% confidence interval 1.67-86.79], P = 0.003 and odds ratio 10.29 [95% confidence interval 1.18-123.47], P = 0.02, respectively). One in 404 women with SLE and 1 in 344 women with SS had 47,XXX. There was an excess of 47,XXX among SLE and SS patients.

CONCLUSION

The estimated prevalence of SLE and SS in women with 47,XXX was ∼2.5 and ∼2.9 times higher, respectively, than that in women with 46,XX and ∼25 and ∼41 times higher, respectively, than that in men with 46,XY. No statistically significant increase of 47,XXX was observed in other female-biased diseases (primary biliary cirrhosis or RA), supporting the idea of multiple pathways to sex bias in autoimmunity.

Identification of a Sjögren's syndrome susceptibility locus at OAS1 that influences isoform switching, protein expression, and responsiveness to type I interferons

Sjögren's syndrome (SS) is a common, autoimmune exocrinopathy distinguished by keratoconjunctivitis sicca and xerostomia. Patients frequently develop serious complications including lymphoma, pulmonary dysfunction, neuropathy, vasculitis, and debilitating fatigue. Dysregulation of type I interferon (IFN) pathway is a prominent feature of SS and is correlated with increased autoantibody titers and disease severity. To identify genetic determinants of IFN pathway dysregulation in SS, we performed cis-expression quantitative trait locus (eQTL) analyses focusing on differentially expressed type I IFN-inducible transcripts identified through a transcriptome profiling study. Multiple cis-eQTLs were associated with transcript levels of 2'-5'-oligoadenylate synthetase 1 (OAS1) peaking at rs10774671 (PeQTL = 6.05 × 10-14). Association of rs10774671 with SS susceptibility was identified and confirmed through meta-analysis of two independent cohorts (Pmeta = 2.59 × 10-9; odds ratio = 0.75; 95% confidence interval = 0.66-0.86). The risk allele of rs10774671 shifts splicing of OAS1 from production of the p46 isoform to multiple alternative transcripts, including p42, p48, and p44. We found that the isoforms were differentially expressed within each genotype in controls and patients with and without autoantibodies. Furthermore, our results showed that the three alternatively spliced isoforms lacked translational response to type I IFN stimulation. The p48 and p44 isoforms also had impaired protein expression governed by the 3' end of the transcripts. The SS risk allele of rs10774671 has been shown by others to be associated with reduced OAS1 enzymatic activity and ability to clear viral infections, as well as reduced responsiveness to IFN treatment. Our results establish OAS1 as a risk locus for SS and support a potential role for defective viral clearance due to altered IFN response as a genetic pathophysiological basis of this complex autoimmune disease.

Large-Scale Identification of Common Trait and Disease Variants Affecting Gene Expression

Genome-wide association studies (GWASs) have identified a multitude of genetic loci involved with traits and diseases. However, it is often unclear which genes are affected in such loci and whether the associated genetic variants lead to increased or decreased gene function. To mitigate this, we integrated associations of common genetic variants in 57 GWASs with 24 studies of expression quantitative trait loci (eQTLs) from a broad range of tissues by using a Mendelian randomization approach. We discovered a total of 3,484 instances of gene-trait-associated changes in expression at a false-discovery rate < 0.05. These genes were often not closest to the genetic variant and were primarily identified in eQTLs derived from pathophysiologically relevant tissues. For instance, genes with expression changes associated with lipid traits were mostly identified in the liver, and those associated with cardiovascular disease were identified in arterial tissue. The affected genes additionally point to biological processes implicated in the interrogated traits, such as the interleukin-27 pathway in rheumatoid arthritis. Further, comparing trait-associated gene expression changes across traits suggests that pleiotropy is a widespread phenomenon and points to specific instances of both agonistic and antagonistic pleiotropy. For instance, expression of SNX19 and ABCB9 is positively correlated with both the risk of schizophrenia and educational attainment. To facilitate interpretation, we provide this lexicon of how common trait-associated genetic variants alter gene expression in various tissues as the online database GWAS2Genes.

Dysregulated CD46 shedding interferes with Th1-contraction in systemic lupus erythematosus

IFN‐γ‐producing T helper 1 (Th1) cell responses mediate protection against infections but uncontrolled Th1 activity also contributes to a broad range of autoimmune diseases. Autocrine complement activation has recently emerged as key in the induction and contraction of human Th1 immunity: activation of the complement regulator CD46 and the C3aR expressed by CD4⁺ T cells via autocrine generated ligands C3b and C3a, respectively, are critical to IFN‐γ production. Further, CD46‐mediated signals also induce co‐expression of immunosuppressive IL‐10 in Th1 cells and transition into a (self)‐regulating and contracting phase. In consequence, C3 or CD46‐deficient patients suffer from recurrent infections while dysregulation of CD46 signaling contributes to Th1 hyperactivity in rheumatoid arthritis and multiple sclerosis. Here, we report a defect in CD46‐regulated Th1 contraction in patients with systemic lupus erythematosus (SLE). We observed that MMP‐9‐mediated increased shedding of soluble CD46 by Th1 cells was associated with this defect and that inhibition of MMP‐9 activity normalized release of soluble CD46 and restored Th1 contraction in patients’ T cells. These data may deliver the first mechanistic explanation for the increased serum CD46 levels observed in SLE patients and indicate that targeting CD46‐cleaving proteases could be a novel avenue to modulate Th1 responses.

Mapping eQTLs with RNA-seq reveals novel susceptibility genes, non-coding RNAs and alternative-splicing events in systemic lupus erythematosus

Studies attempting to functionally interpret complex-disease susceptibility loci by GWAS and eQTL integration have predominantly employed microarrays to quantify gene-expression. RNA-Seq has the potential to discover a more comprehensive set of eQTLs and illuminate the underlying molecular consequence. We examine the functional outcome of 39 variants associated with Systemic Lupus Erythematosus (SLE) through the integration of GWAS and eQTL data from the TwinsUK microarray and RNA-Seq cohort in lymphoblastoid cell lines. We use conditional analysis and a Bayesian colocalisation method to provide evidence of a shared causal-variant, then compare the ability of each quantification type to detect disease relevant eQTLs and eGenes. We discovered the greatest frequency of candidate-causal eQTLs using exon-level RNA-Seq, and identified novel SLE susceptibility genes (e.g. NADSYN1 and TCF7) that were concealed using microarrays, including four non-coding RNAs. Many of these eQTLs were found to influence the expression of several genes, supporting the notion that risk haplotypes may harbour multiple functional effects. Novel SLE associated splicing events were identified in the T-reg restricted transcription factor, IKZF2, and other candidate genes (e.g. WDFY4) through asQTL mapping using the Geuvadis cohort. We have significantly increased our understanding of the genetic control of gene-expression in SLE by maximising the leverage of RNA-Seq and performing integrative GWAS-eQTL analysis against gene, exon, and splice-junction quantifications. We conclude that to better understand the true functional consequence of regulatory variants, quantification by RNA-Seq should be performed at the exon-level as a minimum, and run in parallel with gene and splice-junction level quantification.

Multiple signals at the extended 8p23 locus are associated with susceptibility to systemic lupus erythematosus

BACKGROUND

A major systemic lupus erythematosus (SLE) susceptibility locus lies within a common inversion polymorphism region (encompassing 3.8 - 4.5  Mb) located at 8p23. Initially implicated genes included FAM167A-BLK and XKR6, of which BLK received major attention due to its known role in B-cell biology. Recently, additional SLE risk carried in non-inverted background was also reported.

OBJECTIVE AND METHODS

In this case -control study, we further investigated the 'extended' 8p23 locus (~ 4  Mb) where we observed multiple SLE signals and assessed these signals for their relation to the inversion affecting this region. The study involved a North American discovery data set (~ 1200  subjects) and a replication data set (> 10 000  subjects) comprising European-descent individuals.

RESULTS

Meta-analysis of 8p23 SNPs, with p < 0.05 in both data sets, identified 51 genome-wide significant SNPs (p < 5.0 × 10^-8). While most of these SNPs were related to previously implicated signals (XKR6-FAM167A-BLK subregion), our results also revealed two 'new' SLE signals, including SGK223-CLDN23-MFHAS1 (6.06 × 10^-9 ≤ meta p ≤ 4.88 × 10^-8) and CTSB (meta p = 4.87 × 10^-8) subregions that are located > 2 Mb upstream and ~ 0.3  Mb downstream from previously reported signals. Functional assessment of relevant SNPs indicated putative cis-effects on the expression of various genes at 8p23. Additional analyses in discovery sample, where the inversion genotypes were inferred, replicated the association of non-inverted status with SLE risk and suggested that a number of SLE risk alleles are predominantly carried in non-inverted background.

CONCLUSIONS

Our results implicate multiple (known+novel) SLE signals/genes at the extended 8p23 locus, beyond previously reported signals/genes, and suggest that this broad locus contributes to SLE risk through the effects of multiple genes/pathways.

A Genome-wide Association Study Identifies Risk Alleles in Plasminogen and P4HA2 Associated with Giant Cell Arteritis

Giant cell arteritis (GCA) is the most common form of vasculitis in individuals older than 50 years in Western countries. To shed light onto the genetic background influencing susceptibility for GCA, we performed a genome-wide association screening in a well-powered study cohort. After imputation, 1,844,133 genetic variants were analyzed in 2,134 case subjects and 9,125 unaffected individuals from ten independent populations of European ancestry. Our data confirmed HLA class II as the strongest associated region (independent signals: rs9268905, p = 1.94 × 10^-54, per-allele OR = 1.79; and rs9275592, p = 1.14 × 10^-40, OR = 2.08). Additionally, PLG and P4HA2 were identified as GCA risk genes at the genome-wide level of significance (rs4252134, p = 1.23 × 10^-10, OR = 1.28; and rs128738, p = 4.60 × 10^-9, OR = 1.32, respectively). Interestingly, we observed that the association peaks overlapped with different regulatory elements related to cell types and tissues involved in the pathophysiology of GCA. PLG and P4HA2 are involved in vascular remodelling and angiogenesis, suggesting a high relevance of these processes for the pathogenic mechanisms underlying this type of vasculitis.

Tartrate-Resistant Acid Phosphatase Deficiency in the Predisposition to Systemic Lupus Erythematosus

OBJECTIVE

Mutations in the ACP5 gene, which encodes tartrate-resistant acid phosphatase (TRAP), cause the immuno-osseous disorder spondyloenchondrodysplasia, which includes as disease features systemic lupus erythematosus (SLE) and a type I interferon (IFN) signature. Our aims were to identify TRAP substrates, determine the consequences of TRAP deficiency in immune cells, and assess whether ACP5 mutations are enriched in sporadic cases of SLE.

METHODS

Interaction between TRAP and its binding partners was tested by a yeast 2-hybrid screening, confocal microscopy, and immunoprecipitation/Western blotting. TRAP knockdown was performed using small interfering RNA. Phosphorylation of osteopontin (OPN) was analyzed by mass spectrometry. Nucleotide sequence analysis of ACP5 was performed by Sanger sequencing or next-generation sequencing.

RESULTS

TRAP and OPN colocalized and interacted in human macrophages and plasmacytoid dendritic cells (PDCs). TRAP dephosphorylated 3 serine residues on specific OPN peptides. TRAP knockdown resulted in increased OPN phosphorylation and increased nuclear translocation of IRF7 and P65, with resultant heightened expression of IFN-stimulated genes and IL6 and TNF following Toll-like receptor 9 stimulation. An excess of heterozygous ACP5 missense variants was observed in SLE compared to controls (P = 0.04), and transfection experiments revealed a significant reduction in TRAP activity in a number of variants.

CONCLUSION

Our findings indicate that TRAP and OPN colocalize and that OPN is a substrate for TRAP in human immune cells. TRAP deficiency in PDCs leads to increased IFNα production, providing at least a partial explanation for how ACP5 mutations cause lupus in the context of spondyloenchondrodysplasia. Detection of ACP5 missense variants in a lupus cohort suggests that impaired TRAP functioning may increase susceptibility to sporadic lupus.

Superresolution imaging of the cytoplasmic phosphatase PTPN22 links integrin-mediated T cell adhesion with autoimmunity

Integrins are heterodimeric transmembrane proteins that play a fundamental role in the migration of leukocytes to sites of infection or injury. We found that protein tyrosine phosphatase nonreceptor type 22 (PTPN22) inhibits signaling by the integrin lymphocyte function-associated antigen-1 (LFA-1) in effector T cells. PTPN22 colocalized with its substrates at the leading edge of cells migrating on surfaces coated with the LFA-1 ligand intercellular adhesion molecule-1 (ICAM-1). Knockout or knockdown of PTPN22 or expression of the autoimmune disease-associated PTPN22-R620W variant resulted in the enhanced phosphorylation of signaling molecules downstream of integrins. Superresolution imaging revealed that PTPN22-R620 (wild-type PTPN22) was present as large clusters in unstimulated T cells and that these disaggregated upon stimulation of LFA-1, enabling increased association of PTPN22 with its binding partners at the leading edge. The failure of PTPN22-R620W molecules to be retained at the leading edge led to increased LFA-1 clustering and integrin-mediated cell adhesion. Our data define a previously uncharacterized mechanism for fine-tuning integrin signaling in T cells, as well as a paradigm of autoimmunity in humans in which disease susceptibility is underpinned by inherited phosphatase mutations that perturb integrin function.

Genome-Wide Association Study in an Amerindian Ancestry Population Reveals Novel Systemic Lupus Erythematosus Risk Loci and the Role of European Admixture

OBJECTIVE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease with a strong genetic component. We undertook the present work to perform the first genome-wide association study on individuals from the Americas who are enriched for Native American heritage.

METHODS

We analyzed 3,710 individuals from the US and 4 countries of Latin America who were diagnosed as having SLE, and healthy controls. Samples were genotyped with HumanOmni1 BeadChip. Data on out-of-study controls genotyped with HumanOmni2.5 were also included. Statistical analyses were performed using SNPtest and SNPGWA. Data were adjusted for genomic control and false discovery rate. Imputation was performed using Impute2 and, for classic HLA alleles, HiBag. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated.

RESULTS

The IRF5-TNPO3 region showed the strongest association and largest OR for SLE (rs10488631: genomic control-adjusted P [Pgcadj ] = 2.61 × 10(-29), OR 2.12 [95% CI 1.88-2.39]), followed by HLA class II on the DQA2-DQB1 loci (rs9275572: Pgcadj  = 1.11 × 10(-16), OR 1.62 [95% CI 1.46-1.80] and rs9271366: Pgcadj  = 6.46 × 10(-12), OR 2.06 [95% CI 1.71-2.50]). Other known SLE loci found to be associated in this population were ITGAM, STAT4, TNIP1, NCF2, and IRAK1. We identified a novel locus on 10q24.33 (rs4917385: Pgcadj  = 1.39 × 10(-8)) with an expression quantitative trait locus (eQTL) effect (Peqtl  = 8.0 × 10(-37) at USMG5/miR1307), and several new suggestive loci. SLE risk loci previously identified in Europeans and Asians were corroborated. Local ancestry estimation showed that the HLA allele risk contribution is of European ancestral origin. Imputation of HLA alleles suggested that autochthonous Native American haplotypes provide protection against development of SLE.

CONCLUSION

Our results demonstrate that studying admixed populations provides new insights in the delineation of the genetic architecture that underlies autoimmune and complex diseases.

The complement receptor 3 (CD11b/CD18) agonist Leukadherin-1 suppresses human innate inflammatory signalling

Complement receptor 3 (CR3, CD11b/CD18) is a multi‐functional receptor expressed predominantly on myeloid and natural killer (NK) cells. The R77H variant of CD11b, encoded by the ITGAM rs1143679 polymorphism, is associated robustly with development of the autoimmune disease systemic lupus erythematosus (SLE) and impairs CR3 function, including its regulatory role on monocyte immune signalling. The role of CR3 in NK cell function is unknown. Leukadherin‐1 is a specific small‐molecule CR3 agonist that has shown therapeutic promise in animal models of vascular injury and inflammation. We show that Leukadherin‐1 pretreatment reduces secretion of interferon (IFN)‐γ, tumour necrosis factor (TNF) and macrophage inflammatory protein (MIP)‐1β by monokine‐stimulated NK cells. It was associated with a reduction in phosphorylated signal transducer and activator of transcription (pSTAT)‐5 following interleukin (IL)‐12 + IL‐15 stimulation (P < 0·02) and increased IL‐10 secretion following IL‐12 + IL‐18 stimulation (P < 0·001). Leukadherin‐1 pretreatment also reduces secretion of IL‐1β, IL‐6 and TNF by Toll‐like receptor (TLR)‐2 and TLR‐7/8‐stimulated monocytes (P < 0·01 for all). The R77H variant did not affect NK cell response to Leukadherin‐1 using ex‐vivo cells from homozygous donors; nor did the variant influence CR3 expression by these cell types, in contrast to a recent report. These data extend our understanding of CR3 biology by demonstrating that activation potently modifies innate immune inflammatory signalling, including a previously undocumented role in NK cell function. We discuss the potential relevance of this to the pathogenesis of SLE. Leukadherin‐1 appears to mediate its anti‐inflammatory effect irrespective of the SLE‐risk genotype of CR3, providing further evidence to support its evaluation of Leukadherin‐1 as a potential therapeutic for autoimmune disease.

Genome-wide association meta-analysis in Chinese and European individuals identifies ten new loci associated with systemic lupus erythematosus

Systemic lupus erythematosus (SLE; OMIM 152700) is a genetically complex autoimmune disease. Genome-wide association studies (GWASs) have identified more than 50 loci as robustly associated with the disease in single ancestries, but genome-wide transancestral studies have not been conducted. We combined three GWAS data sets from Chinese (1,659 cases and 3,398 controls) and European (4,036 cases and 6,959 controls) populations. A meta-analysis of these studies showed that over half of the published SLE genetic associations are present in both populations. A replication study in Chinese (3,043 cases and 5,074 controls) and European (2,643 cases and 9,032 controls) subjects found ten previously unreported SLE loci. Our study provides further evidence that the majority of genetic risk polymorphisms for SLE are contained within the same regions across both populations. Furthermore, a comparison of risk allele frequencies and genetic risk scores suggested that the increased prevalence of SLE in non-Europeans (including Asians) has a genetic basis.

Identification of a Systemic Lupus Erythematosus Risk Locus Spanning ATG16L2, FCHSD2, and P2RY2 in Koreans

OBJECTIVE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disorder whose etiology is incompletely understood, but likely involves environmental triggers in genetically susceptible individuals. Using an unbiased genome-wide association (GWA) scan and replication analysis, we sought to identify the genetic loci associated with SLE in a Korean population.

METHODS

A total of 1,174 SLE cases and 4,246 population controls from Korea were genotyped and analyzed with a GWA scan to identify single-nucleotide polymorphisms (SNPs) significantly associated with SLE, after strict quality control measures were applied. For select variants, replication of SLE risk loci was tested in an independent data set of 1,416 SLE cases and 1,145 population controls from Korea and China.

RESULTS

Eleven regions outside the HLA exceeded the genome-wide significance level (P = 5 × 10(-8) ). A novel SNP-SLE association was identified between FCHSD2 and P2RY2, peaking at rs11235667 (P = 1.03 × 10(-8) , odds ratio [OR] 0.59) on a 33-kb haplotype upstream of ATG16L2. In the independent replication data set, the SNP rs11235667 continued to show a significant association with SLE (replication meta-analysis P = 0.001, overall meta-analysis P = 6.67 × 10(-11) ; OR 0.63). Within the HLA region, the SNP-SLE association peaked in the class II region at rs116727542, with multiple independent effects observed in this region. Classic HLA allele imputation analysis identified HLA-DRB1*1501 and HLA-DQB1*0602, each highly correlated with one another, as most strongly associated with SLE. Ten previously established SLE risk loci were replicated: STAT1-STAT4, TNFSF4, TNFAIP3, IKZF1, HIP1, IRF5, BLK, WDFY4, ETS1, and IRAK1-MECP2. Of these loci, previously unreported, independent second risk effects of SNPs in TNFAIP3 and TNFSF4, as well as differences in the association with a putative causal variant in the WDFY4 region, were identified.

CONCLUSION

Further studies are needed to identify true SLE risk effects in other loci suggestive of a significant association, and to identify the causal variants in the regions of ATG16L2, FCHSD2, and P2RY2.

Decreased SMG7 expression associates with lupus-risk variants and elevated antinuclear antibody production

OBJECTIVES

Following up the systemic lupus erythematosus (SLE) genome-wide association studies (GWAS) identification of NMNAT2 at rs2022013, we fine-mapped its 150 kb flanking regions containing NMNAT2 and SMG7 in a 15 292 case-control multi-ancestry population and tested functions of identified variants.

METHODS

We performed genotyping using custom array, imputation by IMPUTE 2.1.2 and allele specific functions using quantitative real-time PCR and luciferase reporter transfections. SLE peripheral blood mononuclear cells (PBMCs) were cultured with small interfering RNAs to measure antinuclear antibody (ANA) and cyto/chemokine levels in supernatants using ELISA.

RESULTS

We confirmed association at NMNAT2 in European American (EA) and Amerindian/Hispanic ancestries, and identified independent signal at SMG7 tagged by rs2702178 in EA only (p=2.4×10^-8, OR=1.23 (95% CI 1.14 to 1.32)). In complete linkage disequilibrium with rs2702178, rs2275675 in the promoter region robustly associated with SMG7 mRNA levels in multiple expression quantitative trait locus (eQTL) datasets. Its risk allele was dose-dependently associated with decreased SMG7 mRNA levels in PBMCs of 86 patients with SLE and 119 controls (p=1.1×10^-3 and 6.8×10^-8, respectively) and conferred reduced transcription activity in transfected HEK-293 (human embryonic kidney cell line) and Raji cells (p=0.0035 and 0.0037, respectively). As a critical component in the nonsense-mediated mRNA decay pathway, SMG7 could regulate autoantigens including ribonucleoprotein (RNP) and Smith (Sm). We showed SMG7 mRNA levels in PBMCs correlated inversely with ANA titres of patients with SLE (r=-0.31, p=0.01), and SMG7 knockdown increased levels of ANA IgG and chemokine (C-C motif) ligand 19 in SLE PBMCs (p=2.0×10^-5 and 2.0×10^-4, respectively).

CONCLUSION

We confirmed NMNAT2 and identified independent SMG7 association with SLE. The inverse relationship between levels of the risk allele-associated SMG7 mRNAs and ANA suggested the novel contribution of mRNA surveillance pathway to SLE pathogenesis.

Identification of a New Susceptibility Locus for Systemic Lupus Erythematosus on Chromosome 12 in Individuals of European Ancestry

OBJECTIVE

Genome-wide association studies (GWAS) in individuals of European ancestry identified a number of systemic lupus erythematosus (SLE) susceptibility loci using earlier versions of high-density genotyping platforms. Followup studies on suggestive GWAS regions using larger samples and more markers identified additional SLE loci in subjects of European descent. This multistage study was undertaken to identify novel SLE loci.

METHODS

In stage 1, we conducted a new GWAS of SLE in a North American case-control sample of subjects of European ancestry (n = 1,166) genotyped on Affymetrix Genome-Wide Human SNP Array 6.0. In stage 2, we further investigated top new suggestive GWAS hits by in silico evaluation and meta-analysis using an additional data set of subjects of European descent (>2,500 individuals), followed by replication of top meta-analysis findings in another data set of subjects of European descent (>10,000 individuals) in stage 3.

RESULTS

As expected, our GWAS revealed the most significant associations at the major histocompatibility complex locus (6p21), which easily surpassed the genome-wide significance threshold (P < 5 × 10(-8)). Several other SLE signals/loci previously implicated in Caucasians and/or Asians were also confirmed in the stage 1 discovery sample, and the strongest signals were observed at 2q32/STAT4 (P = 3.6 × 10(-7)) and at 8p23/BLK (P = 8.1 × 10(-6)). Stage 2 meta-analyses identified a new genome-wide significant SLE locus at 12q12 (meta P = 3.1 × 10(-8)), which was replicated in stage 3.

CONCLUSION

Our multistage study identified and replicated a new SLE locus that warrants further followup in additional studies. Publicly available databases suggest that this newly identified SLE signal falls within a functionally relevant genomic region and near biologically important genes.

Genetic associations of leptin-related polymorphisms with systemic lupus erythematosus

Leptin is abnormally elevated in the plasma of patients with systemic lupus erythematosus (SLE), where it is thought to promote and/or sustain pro-inflammatory responses. Whether this association could reflect an increased genetic susceptibility to develop SLE is not known, and studies of genetic associations with leptin-related polymorphisms in SLE patients have been so far inconclusive. Here we genotyped DNA samples from 15,706 SLE patients and healthy matched controls from four different ancestral groups, to correlate polymorphisms of genes of the leptin pathway to risk for SLE. It was found that although several SNPs showed weak associations, those associations did not remain significant after correction for multiple testing. These data do not support associations between defined leptin-related polymorphisms and increased susceptibility to develop SLE.

Reduced Fluorescence versus Forward Scatter Time-of-Flight and Increased Peak versus Integral Fluorescence Ratios Indicate Receptor Clustering in Flow Cytometry

Clustering of surface receptors is often required to initiate signal transduction, receptor internalization, and cellular activation. To study the kinetics of clustering, we developed an economic high-throughput method using flow cytometry. The quantification of receptor clustering by flow cytometry is based on the following two observations: first, the fluorescence signal length (FL time-of-flight [ToF]) decreases relative to the forward scatter signal length (FSc-ToF), and second, the peak FL (FL-peak) increases relative to the integral FL (FL-integral) upon clustering of FL-labeled surface receptors. Receptor macroclustering can therefore be quantified using the ratios FL-ToF/FSc-ToF (method ToF) or FL-peak/FL-integral (method Peak). We have used these methods to analyze clustering of two immune receptors known to undergo different conformational and oligomeric states: the BCR and the complement receptor 3 (CR3), on murine splenocytes, purified B cells, and human neutrophils. Engagement of both the BCR and CR3, on immortalized as well as primary murine B cells and human neutrophil, respectively, resulted in decreased FL-ToF/FSc-ToF and increased FL-peak/FL-integral ratios. Manipulation of the actin-myosin cytoskeleton altered BCR clustering which could be measured using the established parameters. To confirm clustering of CR3 on neutrophils, we applied imaging flow cytometry. Because receptor engagement is as a biological process dependent on cell viability, energy metabolism, and temperature, receptor clustering can only be quantified by gating on viable cells under physiological conditions. In summary, with this novel method, receptor clustering on nonadherent cells can easily be monitored by high-throughput conventional flow cytometry.

Identification of Susceptibility Loci in IL6, RPS9/LILRB3, and an Intergenic Locus on Chromosome 21q22 in Takayasu Arteritis in a Genome-Wide Association Study

OBJECTIVE

Takayasu arteritis is a rare large vessel vasculitis with incompletely understood etiology. This study was undertaken to perform the first unbiased genome-wide association analysis of Takayasu arteritis.

METHODS

Two independent cohorts of patients with Takayasu arteritis from Turkey and North America were included in our study. The Turkish cohort consisted of 559 patients and 489 controls, and the North American cohort consisted of 134 patients and 1,047 controls of European ancestry. Genotyping was performed using the Omni1-Quad and Omni2.5 genotyping arrays. Genotyping data were subjected to rigorous quality control measures and subsequently analyzed to discover genetic susceptibility loci for Takayasu arteritis.

RESULTS

We identified genetic susceptibility loci for Takayasu arteritis with a genome-wide level of significance in IL6 (rs2069837) (odds ratio [OR] 2.07, P = 6.70 × 10(-9)), RPS9/LILRB3 (rs11666543) (OR 1.65, P = 2.34 × 10(-8)), and an intergenic locus on chromosome 21q22 (rs2836878) (OR 1.79, P = 3.62 × 10(-10)). The genetic susceptibility locus in RPS9/LILRB3 lies within the leukocyte receptor complex gene cluster on chromosome 19q13.4, and the disease risk variant in this locus correlates with reduced expression of multiple genes including the inhibitory leukocyte immunoglobulin-like receptor gene LILRB3 (P = 2.29 × 10(-8)). In addition, we identified candidate susceptibility genes with suggestive levels of association (P < 1 × 10(-5)) with Takayasu arteritis, including PCSK5, LILRA3, PPM1G/NRBP1, and PTK2B.

CONCLUSION

Our findings indicate novel genetic susceptibility loci for Takayasu arteritis and uncover potentially important aspects of the pathophysiology of this form of vasculitis.

Lupus risk variants in the PXK locus alter B-cell receptor internalization

Genome wide association studies have identified variants in PXK that confer risk for humoral autoimmune diseases, including systemic lupus erythematosus (SLE or lupus), rheumatoid arthritis and more recently systemic sclerosis. While PXK is involved in trafficking of epidermal growth factor Receptor (EGFR) in COS-7 cells, mechanisms linking PXK to lupus pathophysiology have remained undefined. In an effort to uncover the mechanism at this locus that increases lupus-risk, we undertook a fine-mapping analysis in a large multi-ancestral study of lupus patients and controls. We define a large (257kb) common haplotype marking a single causal variant that confers lupus risk detected only in European ancestral populations and spans the promoter through the 3′ UTR of PXK. The strongest association was found at rs6445972 with P < 4.62 × 10⁻¹⁰, OR 0.81 (0.75–0.86). Using stepwise logistic regression analysis, we demonstrate that one signal drives the genetic association in the region. Bayesian analysis confirms our results, identifying a 95% credible set consisting of 172 variants spanning 202 kb. Functionally, we found that PXK operates on the B-cell antigen receptor (BCR); we confirmed that PXK influenced the rate of BCR internalization. Furthermore, we demonstrate that individuals carrying the risk haplotype exhibited a decreased rate of BCR internalization, a process known to impact B cell survival and cell fate. Taken together, these data define a new candidate mechanism for the genetic association of variants around PXK with lupus risk and highlight the regulation of intracellular trafficking as a genetically regulated pathway mediating human autoimmunity.

Defective removal of ribonucleotides from DNA promotes systemic autoimmunity

Genome integrity is continuously challenged by the DNA damage that arises during normal cell metabolism. Biallelic mutations in the genes encoding the genome surveillance enzyme ribonuclease H2 (RNase H2) cause Aicardi-Goutières syndrome (AGS), a pediatric disorder that shares features with the autoimmune disease systemic lupus erythematosus (SLE). Here we determined that heterozygous parents of AGS patients exhibit an intermediate autoimmune phenotype and demonstrated a genetic association between rare RNASEH2 sequence variants and SLE. Evaluation of patient cells revealed that SLE- and AGS-associated mutations impair RNase H2 function and result in accumulation of ribonucleotides in genomic DNA. The ensuing chronic low level of DNA damage triggered a DNA damage response characterized by constitutive p53 phosphorylation and senescence. Patient fibroblasts exhibited constitutive upregulation of IFN-stimulated genes and an enhanced type I IFN response to the immunostimulatory nucleic acid polyinosinic:polycytidylic acid and UV light irradiation, linking RNase H2 deficiency to potentiation of innate immune signaling. Moreover, UV-induced cyclobutane pyrimidine dimer formation was markedly enhanced in ribonucleotide-containing DNA, providing a mechanism for photosensitivity in RNase H2-associated SLE. Collectively, our findings implicate RNase H2 in the pathogenesis of SLE and suggest a role of DNA damage-associated pathways in the initiation of autoimmunity.

Preferential association of a functional variant in complement receptor 2 with antibodies to double-stranded DNA

Objectives

Systemic lupus erythematosus (SLE; OMIM 152700) is characterised by the production of antibodies to nuclear antigens. We previously identified variants in complement receptor 2 (CR2/CD21) that were associated with decreased risk of SLE. This study aimed to identify the causal variant for this association.

Methods

Genotyped and imputed genetic variants spanning CR2 were assessed for association with SLE in 15 750 case-control subjects from four ancestral groups. Allele-specific functional effects of associated variants were determined using quantitative real-time PCR, quantitative flow cytometry, electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP)-PCR.

Results

The strongest association signal was detected at rs1876453 in intron 1 of CR2 (p_meta=4.2×10⁻⁴, OR 0.85), specifically when subjects were stratified based on the presence of dsDNA autoantibodies (case-control p_meta=7.6×10⁻⁷, OR 0.71; case-only p_meta=1.9×10⁻⁴, OR 0.75). Although allele-specific effects on B cell CR2 mRNA or protein levels were not identified, levels of complement receptor 1 (CR1/CD35) mRNA and protein were significantly higher on B cells of subjects harbouring the minor allele (p=0.0248 and p=0.0006, respectively). The minor allele altered the formation of several DNA protein complexes by EMSA, including one containing CCCTC-binding factor (CTCF), an effect that was confirmed by ChIP-PCR.

Conclusions

These data suggest that rs1876453 in CR2 has long-range effects on gene regulation that decrease susceptibility to lupus. Since the minor allele at rs1876453 is preferentially associated with reduced risk of the highly specific dsDNA autoantibodies that are present in preclinical, active and severe lupus, understanding its mechanisms will have important therapeutic implications.