Objectives Altered signaling in B-cells is a predominant feature of systemic lupus erythematosus (SLE). The genes BANK1 and BLK were recently described as associated with SLE. BANK1 codes for a B-cell-specific cytoplasmic protein involved in B-cell receptor signaling and BLK codes for an Src tyrosine kinase with important roles in B-cell development. To characterize the role of BANK1 and BLK in SLE, we performed a genetic interaction analysis hypothesizing that genetic interactions could reveal functional pathways relevant to disease pathogenesis. Methods We Used the method GPAT16 to analyze the gene-gene interactions of BANK1 and BLK. Confocal microscopy was used to investigate co-localization, and immunoprecipitation was used to verify the physical interaction of BANK1 and BLK. Results Epistatic interactions between BANK1 and BLK polymorphisms associated with SLE were observed in a discovery set of 279 patients and 515 controls from Northern Europe. A meta-analysis with 4399 European individuals confirmed the genetic interactions between BANK1 and BLK. As BANK1 was identified as a binding partner of the Src tyrosine kinase LYN, we tested the possibility that BANK1 and BLK could also show a protein-protein interaction. We demonstrated co-immunoprecipitation and co-localization of BLK and BANK1. In a Daudi cell line and primary naïve B-cells the endogenous binding was enhanced upon B-cell receptor stimulation using anti-IgM antibodies. Conclusions Here, we show a genetic interaction between BANK1 and BLK, and demonstrate that these molecules interact physically. Our results have important consequences for the understanding of SLE and other autoimmune diseases and identify a potential new signaling pathway.

A B S T R A C TBackground: The clinical and pathologic diversity of systemic lupus erythematosus (SLE) hinders diagnosis, management, and treatment development. This study addresses heterogeneity in SLE through comprehensive molecular phenotyping and machine learning clustering. Methods: Adult SLE patients (n = 198) provided plasma, serum, and RNA. Disease activity was scored by modified SELENA-SLEDAI. Twenty-nine co-expression module scores were calculated from microarray geneexpression data. Plasma soluble mediators (n = 23) and autoantibodies (n = 13) were assessed by multiplex bead-based assays and ELISAs. Patient clusters were identified by machine learning combining K-means clustering and random forest analysis of co-expression module scores and soluble mediators. Findings: SLEDAI scores correlated with interferon, plasma cell, and select cell cycle modules, and with circulating IFN-a, IP10, and IL-1a levels. Co-expression modules and soluble mediators differentiated seven clusters of SLE patients with unique molecular phenotypes. Inflammation and interferon modules were elevated in Clusters 1 (moderately) and 4 (strongly), with decreased T cell modules in Cluster 4. Monocyte, neutrophil, plasmablast, B cell, and T cell modules distinguished the remaining clusters. Active clinical features were similar across clusters. Clinical SLEDAI trended highest in Clusters 3 and 4, though Cluster 3 lacked strong interferon and inflammation signatures. Renal activity was more frequent in Cluster 4, and rare in Clusters 2, 5, and 7. Serology findings were lowest in Clusters 2 and 5. Musculoskeletal and mucocutaneous activity were common in all clusters. Interpretation: Molecular profiles distinguish SLE subsets that are not apparent from clinical information. Prospective longitudinal studies of these profiles may help improve prognostic evaluation, clinical trial design, and precision medicine approaches.

We targeted LYN, a src-tyosine kinase involved in B-cell activation, in case-control association studies using populations of European-American, African-American and Korean subjects. Our combined European-derived population, consisting of 2463 independent cases and 3131 unrelated controls, shows significant association with rs6983130 in a female-only analysis with 2254 cases and 2228 controls (P ¼ 1.1 Â 10 À4 , odds ratio (OR) ¼ 0.81 (95% confidence interval: 0.73-0.90)). This single nucleotide polymorphism (SNP) is located in the 5 0 untranslated region within the first intron near the transcription initiation site of LYN. In addition, SNPs upstream of the first exon also show weak and sporadic association in subsets of the total EuropeanAmerican population. Multivariate logistic regression analysis implicates rs6983130 as a protective factor for systemic lupus erythematosus (SLE) susceptibility when anti-dsDNA, anti-chromatin, anti-52 kDa Ro or anti-Sm autoantibody status were used as covariates. Subset analysis of the European-American female cases by American College of Rheumatology classification criteria shows a reduction in the risk of hematological disorder with rs6983130 compared with cases without hematological disorders (P ¼ 1.5 Â 10 À3 , OR ¼ 0.75 (95% CI: 0.62À0.89)). None of the 90 SNPs tested show significant association with SLE in the African American or Korean populations. These results support an association of LYN with European-derived individuals with SLE, especially within autoantibody or clinical subsets.

Efforts to identify lupus-associated causal variants in the FAM167A/BLK locus on 8p21 are hampered by highly associated noncausal variants. In this report, we used a trans-population mapping and sequencing strategy to identify a common variant (rs922483) in the proximal BLK promoter and a tri-allelic variant (rs1382568) in the upstream alternative BLK promoter as putative causal variants for association with systemic lupus erythematosus. The risk allele (T) at rs922483 reduced proximal promoter activity and modulated alternative promoter usage. Allelic differences at rs1382568 resulted in altered promoter activity in B progenitor cell lines. Thus, our results demonstrated that both lupus-associated functional variants contribute to the autoimmune disease association by modulating transcription of BLK in B cells and thus potentially altering immune responses.

Systemic lupus erythematosus (SLE) is an autoimmune disease with highly variable clinical presentation. Patients suffer from immunological abnormalities that target T-cell, B-cell and accessory cell functions. B cells are hyperactive in SLE patients. An adapter protein expressed in B cells called BANK1 (B-cell scaffold protein with ankyrin repeats) was reported in a previous study to be associated with SLE in a European population. The objective of this study was to assess the BANK1 genotypephenotype association in an independent replication sample. We genotyped 38 single nucleotide polymorphisms (SNPs) in BANK1 on 1892 European-derived SLE patients and 2652 European-derived controls. The strongest associations with SLE and BANK1 were at rs17266594 (corrected P-value ¼ 1.97 Â 10 À5 , odds ratio (OR) ¼ 1.22, 95% CI 1.12-1.34) and rs10516487 (corrected P-value ¼ 2.59 Â 10 À5 , OR ¼ 1.22, 95% CI 1.11-1.34). Our findings suggest that the association is explained by these two SNPs, confirming previous reports that these polymorphisms contribute to the risk of developing lupus. Analysis of patient subsets enriched for hematological, immunological and renal ACR criteria or the levels of autoantibodies, such as anti-RNP A and anti-SmRNP, uncovers additional BANK1 associations. Our results suggest that BANK1 polymorphisms alter immune system development and function to increase the risk for developing lupus.