Variants in BANK1 are associated with lupus nephritis of European ancestry.

The genetic background of lupus nephritis (LN) has not been completely elucidated. We performed a case-only study of 2886 SLE patients, including 947 (33%) with LN. Renal biopsies were available from 396 patients. The discovery cohort (Sweden, n = 1091) and replication cohort 1 (US, n = 962) were genotyped on the Immunochip and replication cohort 2 (Denmark/Norway, n = 833) on a custom array. Patients with LN, proliferative nephritis, or LN with end-stage renal disease were compared with SLE without nephritis. Six loci were associated with LN (p < 1 × 10-4, NFKBIA, CACNA1S, ITGA1, BANK1, OR2Y, and ACER3) in the discovery cohort. Variants in BANK1 showed the strongest association with LN in replication cohort 1 (p = 9.5 × 10-4) and proliferative nephritis in a meta-analysis of discovery and replication cohort 1. There was a weak association between BANK1 and LN in replication cohort 2 (p = 0.052), and in the meta-analysis of all three cohorts the association was strengthened (p = 2.2 × 10-7). DNA methylation data in 180 LN patients demonstrated methylation quantitative trait loci (meQTL) effects between a CpG site and BANK1 variants. To conclude, we describe genetic variations in BANK1 associated with LN and evidence for genetic regulation of DNA methylation within the BANK1 locus. This indicates a role for BANK1 in LN pathogenesis.

Role of the IL-12/IL-35 balance in patients with Sjögren syndrome.

BACKGROUND: An interferon signature is involved in the pathogenesis of primary Sjögren syndrome (pSS), but whether the signature is type 1 or type 2 remains controversial. Mouse models and genetic studies suggest the involvement of TH1 and type 2 interferon pathways. Likewise, polymorphisms of the IL-12A gene (IL12A), which encodes for IL-12p35, have been associated with pSS. The IL-12p35 subunit is shared by 2 heterodimers: IL-12 and IL-35. OBJECTIVE: We sought to confirm genetic association of the IL12A polymorphism and pSS and elucidate involvement of the IL-12/IL-35 balance in patients with pSS by using functional studies. METHODS: The genetic study involved 673 patients with pSS from 2 French pSS cohorts and 585 healthy French control subjects. Functional studies were performed on sorted monocytes, irrespective of whether they were stimulated. IL12A mRNA expression and IL-12 and IL-35 protein levels were assessed by using quantitative RT-PCR and ELISA and a multiplex kit for IL-35 and IL-12, respectively. RESULTS: We confirmed association of the IL12A rs485497 polymorphism and pSS and found an increased serum protein level of IL-12p70 in patients with pSS carrying the risk allele (P = .016). Serum levels of IL-12p70 were greater in patients than control subjects (P = .0001), especially in patients with more active disease (P = .05); conversely, IL-35 levels were decreased in patients (P = .0001), especially in patients with more active disease (P = .05). In blood cellular subsets both IL12p35 and EBV-induced gene protein 3 (EBI3) mRNAs were detected only in B cells, with a trend toward a lower level among patients with pSS. CONCLUSION: Our findings emphasize involvement of the IL-12/IL-35 balance in the pathogenesis of pSS. Serum IL-35 levels were associated with low disease activity, in contrast with serum IL-12p70 levels, which were associated with more active disease.

Candidate long-range regulatory sites acting on the IL17 pathway genes TRAF3IP2 and IL17RA are associated with psoriasis.

BACKGROUND: Drug-mediated disruption of IL17A, IL17F and IL17RA proteins is effective in psoriasis. However, disruption of the IL17 pathway by functional mutations has so far only been shown to affect risk in IL23R and TRAF3IP2. It is unclear whether this is due to rarity of disruptive mutations. OBJECTIVE: (a) To delineate the prevalence of mutations in key IL17 pathway genes and (b) to identify candidate regulatory sites acting on IL23R, IL17A, IL17RA and TRAF3IP2 from a distance. METHODS: Extraction of mutation frequencies from ExAc data, evolutionary sequence alignment; mapping of long-range interacting (LRI) enhancers; and genetic association testing in a novel psoriasis cohort. RESULTS: The prevalence of disruptive mutations in genes such as IL17RA is sufficient to have been detectable by existing data sets. Therefore, lack of their association with psoriasis indicates that genetic risk primarily resides in variants acting from a distance. We identify two LRI enhancer sites, regulating IL17RA and TRAF3IP2, respectively. The TRAF3IP2 regulator localizes to the TRAF3IP2 antisense promoter, suggesting feedback regulation. Both LRI sites are associated with psoriasis in a novel Scottish psoriasis cohort and the TRAF3IP2-LRI at rs71562294 replicates in the WTCCC cohort. CONCLUSION: Genetic risk for psoriasis may be encoded at LRI sites regulating IL17 pathway genes from a distance.

Multiple polymorphisms in the TNFAIP3 region are independently associated with systemic lupus erythematosus.

The TNFAIP3 (tumor necrosis factor alpha-induced protein 3) gene encodes a ubiquitin editing enzyme, A20, that restricts NF-kappaB-dependent signaling and prevents inflammation. We show that three independent SNPs in the TNFAIP3 region (rs13192841, rs2230926 and rs6922466) are associated with systemic lupus erythematosus (SLE) among individuals of European ancestry. These findings provide critical links between A20 and the etiology of SLE.

Unraveling multiple MHC gene associations with systemic lupus erythematosus: model choice indicates a role for HLA alleles and non-HLA genes in Europeans.

We have performed a meta-analysis of the major-histocompatibility-complex (MHC) region in systemic lupus erythematosus (SLE) to determine the association with both SNPs and classical human-leukocyte-antigen (HLA) alleles. More specifically, we combined results from six studies and well-known out-of-study control data sets, providing us with 3,701 independent SLE cases and 12,110 independent controls of European ancestry. This study used genotypes for 7,199 SNPs within the MHC region and for classical HLA alleles (typed and imputed). Our results from conditional analysis and model choice with the use of the Bayesian information criterion show that the best model for SLE association includes both classical loci (HLA-DRB1(∗)03:01, HLA-DRB1(∗)08:01, and HLA-DQA1(∗)01:02) and two SNPs, rs8192591 (in class III and upstream of NOTCH4) and rs2246618 (MICB in class I). Our approach was to perform a stepwise search from multiple baseline models deduced from a priori evidence on HLA-DRB1 lupus-associated alleles, a stepwise regression on SNPs alone, and a stepwise regression on HLA alleles. With this approach, we were able to identify a model that was an overwhelmingly better fit to the data than one identified by simple stepwise regression either on SNPs alone (Bayes factor [BF] > 50) or on classical HLA alleles alone (BF > 1,000).

Germline and somatic genetic variations of TNFAIP3 in lymphoma complicating primary Sjogren's syndrome.

Several autoimmune diseases, including primary Sjögren's syndrome (pSS), are associated with an increased risk for lymphoma. Polymorphisms of TNFAIP3, which encodes the A20 protein that plays a key role in controlling nuclear factor κB activation, have been associated with several autoimmune diseases. Somatic mutations of TNFAIP3 have been observed in the mucosa-associated lymphoid tissue lymphoma subtype frequently associated with pSS. We studied germline and somatic abnormalities of TNFAIP3 in 574 patients with pSS, including 25 with lymphoma. Nineteen additional patients with pSS and lymphoma were available for exome sequence analysis. Functional abnormalities of A20 were assessed by gene reporter assays. The rs2230926 exonic variant was associated with an increased risk for pSS complicated by lymphoma (odds ratio, 3.36 [95% confidence interval, 1.34-8.42], and odds ratio, 3.26 [95% confidence interval, 1.31-8.12], vs controls and pSS patients without lymphoma, respectively; P = .011). Twelve (60%) of the 20 patients with paired germline and lymphoma TNFAIP3 sequence data had functional abnormalities of A20: 6 in germline DNA, 5 in lymphoma DNA, and 1 in both. The frequency was even higher (77%) among pSS patients with mucosa-associated lymphoid tissue lymphoma. Some of these variants showed impaired control of nuclear factor κB activation. These results support a key role for germline and somatic variations of A20 in the transformation between autoimmunity and lymphoma.

Differentiation of follicular helper T cells by salivary gland epithelial cells in primary Sjögren's syndrome.

Follicular helper T cells (Tfh), which play a pivotal role in B cell activation and differentiation in lymphoid structures, secrete IL-21 whose augmented secretion is a hallmark of several autoimmune diseases. To decipher the cellular and molecular interactions occurring in salivary glands of patients suffering from primary Sjögren's syndrome (pSS), we investigated whether salivary gland epithelial cells (SGECs) were capable to induce Tfh differentiation. Co-cultures of naïve CD4(+) T cells and SGECs from both patients with pSS and controls were performed. Here, we report that IL-6 and ICOSL expression by SGECs contributes to naïve CD4(+) T differentiation into Tfh cells, as evidenced by their acquisition of a specific phenotype, characterized by Bcl-6, ICOS and CXCR5 expression and IL-21 secretion, but also but by their main functional feature: the capacity to enhance B lymphocytes survival. We demonstrated an increase of serum IL-21 with systemic activity. Finally, we analyzed the potential occurrence of a genetic association between IL-21 or IL-21R gene polymorphisms and pSS or elevated IL-21 secretion. This study, which demonstrates a direct induction of Tfh differentiation by SGECs, emphasizes a yet unknown pathogenic role of SGECs and suggests that Tfh and IL-21 might be relevant biomarkers and/or therapeutic targets in primary Sjögren's syndrome.

Brief Report: Single-nucleotide polymorphisms in VKORC1 are risk factors for systemic lupus erythematosus in Asians.

OBJECTIVE: The increased risk of thrombosis in systemic lupus erythematosus (SLE) may be partially explained by interrelated genetic pathways for thrombosis and SLE. The present study was undertaken to investigate whether 33 established and novel single-nucleotide polymorphisms (SNPs) in 20 genes involved in hemostasis pathways that have been associated with deep venous thrombosis (DVT) in the general population are risk factors for SLE among Asian subjects. METHODS: Patients in the discovery cohort were enrolled in 1 of 2 North American SLE cohorts. Patients in the replication cohort were enrolled in 1 of 4 Asian or 2 North American cohorts. We first genotyped 263 Asian patients with SLE and 357 healthy Asian control subjects for 33 SNPs in the discovery phase, and then genotyped 5 SNPs in up to an additional 1,496 patients and 993 controls in the replication phase. Patients were compared to controls for bivariate association with minor alleles. Principal components analysis was used to control for intra-Asian ancestry in the replication cohort. RESULTS: Two genetic variants in the gene VKORC1 were highly significant in both the discovery and replication cohorts: rs9934438 (in the discovery cohort, odds ratio [OR] 2.45, P=2×10(-9); in the replication cohort, OR 1.54, P=4×10(-6)) and rs9923231 (in the discovery cohort, OR 2.40, P=6×10(-9); in the replication cohort, OR 1.53, P=5×10(-6)). These associations were significant in the replication cohort after adjustment for intra-Asian ancestry: for rs9934438, OR 1.34, P=0.0029; for rs9923231, OR 1.34, P=0.0032. CONCLUSION: Genetic variants in VKORC1, which are involved in vitamin K reduction and associated with DVT, correlate with SLE development in Asian subjects. These results suggest that there may be intersecting genetic pathways for the development of SLE and thrombosis.

Differential genetic associations for systemic lupus erythematosus based on anti-dsDNA autoantibody production.

Systemic lupus erythematosus (SLE) is a clinically heterogeneous, systemic autoimmune disease characterized by autoantibody formation. Previously published genome-wide association studies (GWAS) have investigated SLE as a single phenotype. Therefore, we conducted a GWAS to identify genetic factors associated with anti-dsDNA autoantibody production, a SLE-related autoantibody with diagnostic and clinical importance. Using two independent datasets, over 400,000 single nucleotide polymorphisms (SNPs) were studied in a total of 1,717 SLE cases and 4,813 healthy controls. Anti-dsDNA autoantibody positive (anti-dsDNA +, n = 811) and anti-dsDNA autoantibody negative (anti-dsDNA -, n = 906) SLE cases were compared to healthy controls and to each other to identify SNPs associated specifically with these SLE subtypes. SNPs in the previously identified SLE susceptibility loci STAT4, IRF5, ITGAM, and the major histocompatibility complex were strongly associated with anti-dsDNA + SLE. Far fewer and weaker associations were observed for anti-dsDNA - SLE. For example, rs7574865 in STAT4 had an OR for anti-dsDNA + SLE of 1.77 (95% CI 1.57-1.99, p = 2.0E-20) compared to an OR for anti-dsDNA - SLE of 1.26 (95% CI 1.12-1.41, p = 2.4E-04), with p(heterogeneity)<0.0005. SNPs in the SLE susceptibility loci BANK1, KIAA1542, and UBE2L3 showed evidence of association with anti-dsDNA + SLE and were not associated with anti-dsDNA - SLE. In conclusion, we identified differential genetic associations with SLE based on anti-dsDNA autoantibody production. Many previously identified SLE susceptibility loci may confer disease risk through their role in autoantibody production and be more accurately described as autoantibody propensity loci. Lack of strong SNP associations may suggest that other types of genetic variation or non-genetic factors such as environmental exposures have a greater impact on susceptibility to anti-dsDNA - SLE.

Screening of Environmental Chemicals to Characterize Exposures in Participants With Systemic Lupus Erythematosus.

OBJECTIVE: There is a need to characterize exposures associated with the pathogenesis of systemic lupus erythematosus (SLE). In this pilot study, we explore a hypothesis-free approach that can measure thousands of exogenous chemicals in blood ("exposome") in patients with SLE and unaffected controls. METHODS: This cross-sectional study analyzed a cohort of patients with prevalent SLE (n = 285) and controls (n = 106). Plasma was analyzed by liquid chromatography-quadrupole time-of-flight mass spectrometry (LC-QTOF/MS). Mass spectrometry features present in at least 25% of all samples were selected for association analysis (n = 2,737). Features were matched to potential chemicals using available databases. Association analysis of abundances of features with SLE status was performed, adjusting for age and sex. We also explored features associated with SLE phenotypes, sociodemographic factors, and current medication use. RESULTS: We found 30 features significantly associated with SLE status (Bonferroni P < 0.05). Of these, seven matched chemical names based on databases. These seven features included phthalate metabolites, a formetanate metabolite, and eugenol. The abundance of acid pesticides differed between patients with SLE and controls (Bonferroni P < 0.05). Two unmatched features were associated with a history of lupus nephritis, and one with anti-double-stranded DNA antibody production (Bonferroni P < 0.05). Seventeen features varied by self-reported race and ethnicity, including a polyfluoroalkyl substance (analysis of variance P < 1.69 × 10-5). Eleven features correlated with antimalarials, 6 with mycophenolate mofetil, and 29 with prednisone use. CONCLUSION: This proof-of-concept study demonstrates that LC-QTOF/MS is a powerful tool that agnostically detects circulating exogenous compounds. These analyses can generate hypotheses of disease-related exposures for future prospective, longitudinal studies.

Local Ancestry at the Major Histocompatibility Complex Region is Not a Major Contributor to Disease Heterogeneity in a Multiethnic Lupus Cohort.

OBJECTIVE: Systemic lupus erythematosus (SLE) is an autoimmune disease resulting in debilitating clinical manifestations that vary in severity by race and ethnicity with a disproportionate burden in African American, Mestizo, and Asian populations compared with populations of European descent. Differences in global and local genetic ancestry may shed light on the underlying mechanisms contributing to these disparities, including increased prevalence of lupus nephritis, younger age of symptom onset, and presence of autoantibodies. METHODS: A total of 1,139 European, African American, and Mestizos patients with SLE were genotyped using the Affymetrix LAT1 World array. Global ancestry proportions were estimated using ADMIXTURE, and local ancestry was estimated using RFMIXv2.0. We investigated associations between lupus nephritis, age at onset, and autoantibody status with both global and local ancestry proportions within the Major Histocompatibility Complex region. RESULTS: Our results showed small effect sizes that did not meet the threshold for statistical significance for global or local ancestry proportions in either African American or Mestizo patients with SLE who presented with the clinical manifestations of interest compared with those who did not. CONCLUSION: These findings suggest that local genetic ancestry within the Major Histocompatibility Complex region is not a major contributor to these SLE manifestations among patients with SLE from admixed populations.

Single-cell RNA-seq analysis reveals cell subsets and gene signatures associated with Rheumatoid Arthritis Disease Activity.

Rheumatoid arthritis (RA) management lean toward achieving remission or low-disease activity. In this study, we conducted single-cell RNA sequencing (scRNAseq) of peripheral blood mononuclear cells (PBMCs) from 36 individuals (18 RA patients and 18 matched controls, accounting for age, sex, race, and ethnicity), to identify disease-relevant cell subsets and cell type-specific signatures associated with disease activity. Our analysis revealed 18 distinct PBMC subsets, including an IFITM3 overexpressing Interferon-activated (IFN-activated) monocyte subset. We observed an increase in CD4+ T effector memory cells in patients with moderate to high disease activity (DAS28-CRP ≥ 3.2), and a decrease in non-classical monocytes in patients with low disease activity or remission (DAS28-CRP < 3.2). Pseudobulk analysis by cell type identified 168 differentially expressed genes between RA and matched controls, with a downregulation of pro-inflammatory genes in the gamma-delta T cells subset, alteration of genes associated with RA predisposition in the IFN-activated subset, and non-classical monocytes. Additionally, we identified a gene signature associated with moderate-high disease activity, characterized by upregulation of pro-inflammatory genes such as TNF, JUN, EGR1, IFIT2, MAFB, G0S2, and downregulation of genes including HLA-DQB1, HLA-DRB5, TNFSF13B. Notably, cell-cell communication analysis revealed an upregulation of signaling pathways, including VISTA, in both moderate-high and remission-low disease activity contexts. Our findings provide valuable insights into the systemic cellular and molecular mechanisms underlying RA disease activity.

Identification of Cell-Specific Differential DNA Methylation Associated With Methotrexate Treatment Response in Rheumatoid Arthritis.

OBJECTIVE: We undertook this study to estimate changes in cell-specific DNA methylation (DNAm) associated with methotrexate (MTX) response using whole blood samples collected from rheumatoid arthritis (RA) patients before and after initiation of MTX treatment. METHODS: Patients included in this study were from the Rheumatoid Arthritis Medication Study (n = 66) and the University of California San Francisco Rheumatoid Arthritis study (n = 11). All patients met the American College of Rheumatology RA classification criteria. Blood samples were collected at baseline and following treatment. Disease Activity Scores in 28 joints using the C-reactive protein level were collected at baseline and after 3-6 months of treatment with MTX. Methylation profiles were generated using the Illumina Infinium HumanMethylation450 and MethylationEPIC v1.0 BeadChip arrays using DNA from whole blood. MTX response was defined using the EULAR response criteria (responders showed good/moderate response; nonresponders showed no response). Differentially methylated positions were identified using the Limma software package and Tensor Composition Analysis, which is a method for identifying cell-specific differential DNAm at the CpG level from tissue-level ("bulk") data. Differentially methylated regions were identified using Comb-p software. RESULTS: We found evidence of differential global methylation between treatment response groups. Further, we found patterns of cell-specific differential global methylation associated with MTX response. After correction for multiple testing, 1 differentially methylated position was associated with differential DNAm between responders and nonresponders at baseline in CD4+ T cells, CD8+ T cells, and natural killer cells. Thirty-nine cell-specific differentially methylated regions associated with MTX treatment response were identified. There were no significant findings in analyses of whole blood samples. CONCLUSION: We identified cell-specific changes in DNAm that were associated with MTX treatment response in RA patients. Future studies of DNAm and MTX treatment response should include measurements of DNAm from sorted cells.

GWAS meta-analysis of psoriasis identifies new susceptibility alleles impacting disease mechanisms and therapeutic targets.

Psoriasis is a common, debilitating immune-mediated skin disease. Genetic studies have identified biological mechanisms of psoriasis risk, including those targeted by effective therapies. However, the genetic liability to psoriasis is not fully explained by variation at robustly identified risk loci. To move towards a saturation map of psoriasis susceptibility we meta-analysed 18 GWAS comprising 36,466 cases and 458,078 controls and identified 109 distinct psoriasis susceptibility loci, including 45 that have not been previously reported. These include susceptibility variants at loci in which the therapeutic targets IL17RA and AHR are encoded, and deleterious coding variants supporting potential new drug targets (including in STAP2, CPVL and POU2F3). We conducted a transcriptome-wide association study to identify regulatory effects of psoriasis susceptibility variants and cross-referenced these against single cell expression profiles in psoriasis-affected skin, highlighting roles for the transcriptional regulation of haematopoietic cell development and epigenetic modulation of interferon signalling in psoriasis pathobiology.

High-density SNP screening of the major histocompatibility complex in systemic lupus erythematosus demonstrates strong evidence for independent susceptibility regions.

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99 x 10(-16)). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53 x 10(-12)), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80 x 10(-13). Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE-associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.

Dynamics of Methylation of CpG Sites Associated With Systemic Lupus Erythematosus Subtypes in a Longitudinal Cohort.

OBJECTIVE: Findings from cross-sectional studies have revealed associations between DNA methylation and systemic lupus erythematosus (SLE) outcomes. This study was undertaken to investigate the dynamics of DNA methylation by examining participants from an SLE longitudinal cohort using samples collected at 2 time points. METHODS: A total of 101 participants from the California Lupus Epidemiology Study were included in our analysis. DNA was extracted from blood samples collected at the time of enrolment in the cohort and samples collected after 2 years and was analyzed using Illumina EPIC BeadChip kit. Paired t-tests were used to identify genome-wide changes which included 256 CpG sites previously found to be associated with SLE subtypes. Linear mixed models were developed to understand the relationship between DNA methylation and disease activity, medication use, and sample cell-type proportions, adjusted for age, sex, and genetic principal components. RESULTS: The majority of CpGs that were previously determined to be associated with SLE subtypes remained stable over 2 years (185 CpGs [72.3%]; t-test false discovery rate >0.05). Compared to background genome-wide methylation, there was an enrichment of SLE subtype-associated CpGs that changed over time (27.7% versus 0.34%). Changes in cell-type proportions were associated with changes at 67 CpGs (P < 2.70 × 10-5 ), and 15 CpGs had at least 1 significant association with immunosuppressant use. CONCLUSION: In this longitudinal SLE cohort, we identified a subset of SLE subtype-associated CpGs that remained stable over time and may be useful as biomarkers of disease subtypes. Another subset of SLE subtype-associated CpGs changed at a higher proportion compared to the genome-wide methylome. Additional studies are needed to understand the etiology and impact of these changes on methylation of SLE-associated CpGs.

A Human MSH6 Germline Variant Associated With Systemic Lupus Erythematosus Induces Lupus-like Disease in Mice.

OBJECTIVE: To determine if single-nucleotide polymorphisms (SNPs) in DNA repair genes are enriched in individuals with systemic lupus erythematosus (SLE) and if they are sufficient to confer a disease phenotype in a mouse model. METHODS: Human exome chip data of 2499 patients with SLE and 1230 healthy controls were analyzed to determine if variants in 10 different mismatch repair genes (MSH4, EXO1, MSH2, MSH6, MLH1, MSH3, POLH, PMS2, ML3, and APEX2) were enriched in individuals with SLE. A mouse model of the MSH6 SNP, which was found to be enriched in individuals with SLE, was created using CRISPR/Cas9 gene targeting. Wildtype mice and mice heterozygous and homozygous for the MSH6 variant were then monitored for 2 years for the development of autoimmune phenotypes, including the presence of high levels of antinuclear antibodies (ANA). Additionally, somatic hypermutation frequencies and spectra of the intronic region downstream of the VH J558-rearranged JH4 immunoglobulin gene was characterized from Peyer's patches. RESULTS: Based on the human exome chip data, the MSH6 variant (rs63750897, p.Ser503Cys) is enriched among patients with SLE versus controls after we corrected for ancestry (odds ratio = 8.39, P = 0.0398). Mice homozygous for the MSH6 variant (Msh6S502C/S502C ) harbor significantly increased levels of ANA. Additionally, the Msh6S502C/S502C mice display a significant increase in the infiltration of CD68+ cells (a marker for monocytes and macrophages) into the lung alveolar space as well as apoptotic cells. Furthermore, characterization of somatic hypermutation in these mice reveals an increase in the DNA polymerase η mutational signature. CONCLUSION: An MSH6 mutation that is enriched in humans diagnosed with lupus was identified. Mice harboring this Msh6 mutation develop increased autoantibodies and an inflammatory lung disease. These results suggest that the human MSH6 variant is linked to the development of SLE.

Mer-tyrosine kinase: a novel susceptibility gene for SLE related end-stage renal disease.

OBJECTIVE: Lupus nephritis (LN) is a common and severe manifestation of SLE. The genetic risk for nephritis and progression to end-stage renal disease (ESRD) in patients with LN remains unclear. Herein, we aimed to identify novel genetic associations with LN, focusing on subphenotypes and ESRD. METHODS: We analysed genomic data on 958 patients with SLE (discovery cohort: LN=338) with targeted sequencing data from 1832 immunological pathway genes. We used an independent multiethnic cohort comprising 1226 patients with SLE (LN=603) as a replication dataset. Detailed functional annotation and functional epigenomic enrichment analyses were applied to predict functional effects of the candidate variants. RESULTS: A genetic variant (rs56097910) within the MERTK gene was associated with ESRD in both cohorts, meta-analysis OR=5.4 (2.8 to 10.6); p=1.0×10-6. We observed decreased methylation levels in peripheral blood cells from SLE patients with ESRD, compared with patients without renal SLE (p=2.7×10-4), at one CpG site (cg16333401) in close vicinity to the transcription start site of MERTK and located in a DNAse hypersensitivity region in T and B cells. Rs56097910 is linked to altered MERTK expression in kidney tissue in public eQTL databases. Two loci were replicated for association with proliferative LN: PRDM1 (rs6924535, pmeta=1.6×10-5, OR=0.58) and APOA1BP (NAXE) (rs942960, pmeta=1.2×10-5, OR=2.64). CONCLUSION: We identified a novel genetic risk locus, MERTK, associated with SLE-ESRD using the data from two large SLE cohorts. Through DNA methylation analysis and functional annotation, we showed that the risk could be mediated through regulation of gene expression. Our results suggest that variants in the MERTK gene are important for the risk of developing SLE-ESRD and suggest a role for PRDM1 and APOA1BP in proliferative LN.

P2RY8 variants in lupus patients uncover a role for the receptor in immunological tolerance.

B cell self-tolerance is maintained through multiple checkpoints, including restraints on intracellular signaling and cell trafficking. P2RY8 is a receptor with established roles in germinal center (GC) B cell migration inhibition and growth regulation. Somatic P2RY8 variants are common in GC-derived B cell lymphomas. Here, we identify germline novel or rare P2RY8 missense variants in lupus kindreds or the related antiphospholipid syndrome, including a "de novo" variant in a child with severe nephritis. All variants decreased protein expression, F-actin abundance, and GPCR-RhoA signaling, and those with stronger effects increased AKT and ERK activity and cell migration. Remarkably, P2RY8 was reduced in B cell subsets from some SLE patients lacking P2RY8 gene variants. Low P2RY8 correlated with lupus nephritis and increased age-associated B cells and plasma cells. By contrast, P2RY8 overexpression in cells and mice restrained plasma cell development and reinforced negative selection of DNA-reactive developing B cells. These findings uncover a role of P2RY8 in immunological tolerance and lupus pathogenesis.

Rheumatoid arthritis risk allele PTPRC is also associated with response to anti-tumor necrosis factor alpha therapy.

OBJECTIVE: Anti-tumor necrosis factor alpha (anti-TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). The aim of the present study was to test established RA genetic risk factors to determine whether the same alleles also influence the response to anti-TNF therapy. METHODS: A total of 1,283 RA patients receiving etanercept, infliximab, or adalimumab therapy were studied from among an international collaborative consortium of 9 different RA cohorts. The primary end point compared RA patients with a good treatment response according to the European League Against Rheumatism (EULAR) response criteria (n = 505) with RA patients considered to be nonresponders (n = 316). The secondary end point was the change from baseline in the level of disease activity according to the Disease Activity Score in 28 joints (triangle upDAS28). Clinical factors such as age, sex, and concomitant medications were tested as possible correlates of treatment response. Thirty-one single-nucleotide polymorphisms (SNPs) associated with the risk of RA were genotyped and tested for any association with treatment response, using univariate and multivariate logistic regression models. RESULTS: Of the 31 RA-associated risk alleles, a SNP at the PTPRC (also known as CD45) gene locus (rs10919563) was associated with the primary end point, a EULAR good response versus no response (odds ratio [OR] 0.55, P = 0.0001 in the multivariate model). Similar results were obtained using the secondary end point, the triangle upDAS28 (P = 0.0002). There was suggestive evidence of a stronger association in autoantibody-positive patients with RA (OR 0.55, 95% confidence interval [95% CI] 0.39-0.76) as compared with autoantibody-negative patients (OR 0.90, 95% CI 0.41-1.99). CONCLUSION: Statistically significant associations were observed between the response to anti-TNF therapy and an RA risk allele at the PTPRC gene locus. Additional studies will be required to replicate this finding in additional patient collections.