Whole Exome Sequencing of Patients from Multicase Families with Systemic Lupus Erythematosus Identifies Multiple Rare Variants.

In an effort to identify rare alleles associated with SLE, we have performed whole exome sequencing of the most distantly related affected individuals from two large Icelandic multicase SLE families followed by Ta targeted genotyping of additional relatives. We identified multiple rare likely pathogenic variants in nineteen genes co-segregating with the disease through multiple generations. Gene co-expression and protein-protein interaction analysis identified a network of highly connected genes comprising several loci previously implicated in autoimmune diseases. These genes were significantly enriched for immune system development, lymphocyte activation, DNA repair, and V(D)J gene recombination GO-categories. Furthermore, we found evidence of aggregate association and enrichment of rare variants at the FAM71E1/EMC10 locus in an independent set of 4,254 European SLE-cases and 4,349 controls. Our study presents evidence supporting that multiple rare likely pathogenic variants, in newly identified genes involved in known disease pathogenic pathways, segregate with SLE at the familial and population level.

PXK locus in systemic lupus erythematosus: fine mapping and functional analysis reveals novel susceptibility gene ABHD6.

OBJECTIVES: To perform fine mapping of the PXK locus associated with systemic lupus erythematosus (SLE) and study functional effects that lead to susceptibility to the disease. METHODS: Linkage disequilibrium (LD) mapping was conducted by using 1251 SNPs (single nucleotide polymorphism) covering a 862 kb genomic region on 3p14.3 comprising the PXK locus in 1467 SLE patients and 2377 controls of European origin. Tag SNPs and genotypes imputed with IMPUTE2 were tested for association by using SNPTEST and PLINK. The expression QTLs data included three independent datasets for lymphoblastoid cells of European donors: HapMap3, MuTHER and the cross-platform eQTL catalogue. Correlation analysis of eQTLs was performed using Vassarstats. Alternative splicing for the PXK gene was analysed on mRNA from PBMCs. RESULTS: Fine mapping revealed long-range LD (>200 kb) extended over the ABHD6, RPP14, PXK, and PDHB genes on 3p14.3. The highly correlated variants tagged an SLE-associated haplotype that was less frequent in the patients compared with the controls (OR=0.89, p=0.00684). A robust correlation between the association with SLE and enhanced expression of ABHD6 gene was revealed, while neither expression, nor splicing alterations associated with SLE susceptibility were detected for PXK. The SNP allele frequencies as well as eQTL pattern analysed in the CEU and CHB HapMap3 populations indicate that the SLE association and the effect on ABHD6 expression are specific to Europeans. CONCLUSIONS: These results confirm the genetic association of the locus 3p14.3 with SLE in Europeans and point to the ABHD6 and not PXK, as the major susceptibility gene in the region. We suggest a pathogenic mechanism mediated by the upregulation of ABHD6 in individuals carrying the SLE-risk variants.

Fine mapping and conditional analysis identify a new mutation in the autoimmunity susceptibility gene BLK that leads to reduced half-life of the BLK protein.

OBJECTIVES: To perform fine mapping of the autoimmunity susceptibility gene BLK and identify functional variants involved in systemic lupus erythematosus (SLE). METHODS: Genotyping of 1163 European SLE patients and 1482 controls and imputation were performed covering the BLK gene with 158 single-nucleotide polymorphisms. Logistic regression analysis was done using PLINK and conditional analyses using GENABEL's test score. Transfections of BLK constructs on HEK293 cells containing the novel mutation or the wild type form were analysed for their effect on protein half-life using a protein stability assay, cycloheximide and western blot. CHiP-qPCR for detection of nuclear factor κ B (NFkB) binding. RESULTS: Fine mapping of BLK identified two independent genetic effects with functional consequences: one represented by two tightly linked associated haplotype blocks significantly enriched for NFκB-binding sites and numerous putative regulatory variants whose risk alleles correlated with low BLK mRNA levels. Binding of NFkBp50 and p65 to an associated 1.2 Kb haplotype segment was confirmed. A second independent genetic effect was represented by an Ala71Thr, low-frequency missense substitution with an OR=2.31 (95% CI 1.38 to 3.86). The 71Thr decreased BLK protein half-life. CONCLUSIONS: These results show that rare and common regulatory variants in BLK are involved in disease susceptibility and both, albeit independently, lead to reduced levels of BLK protein.

Genetic and physical interaction of the B-cell systemic lupus erythematosus-associated genes BANK1 and BLK.

OBJECTIVES: Altered signalling 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 signalling and BLK codes for an Src tyrosine kinase with important roles in B-cell development. To characterise the role of BANK1 and BLK in SLE, a genetic interaction analysis was performed hypothesising that genetic interactions could reveal functional pathways relevant to disease pathogenesis. METHODS: The GPAT16 method was used to analyse the gene-gene interactions of BANK1 and BLK. Confocal microscopy was used to investigate co-localisation, 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, the possibility that BANK1 and BLK could also show a protein-protein interaction was tested. The co-immunoprecipitation and co-localisation of BLK and BANK1 were demonstrated. In a Daudi cell line and primary naive B cells endogenous binding was enhanced upon B-cell receptor stimulation using anti-IgM antibodies. CONCLUSION: This study shows a genetic interaction between BANK1 and BLK, and demonstrates that these molecules interact physically. The results have important consequences for the understanding of SLE and other autoimmune diseases and identify a potential new signalling pathway.

A 3'-untranslated region variant is associated with impaired expression of CD226 in T and natural killer T cells and is associated with susceptibility to systemic lupus erythematosus.

OBJECTIVE: Costimulatory receptor CD226 plays an important role in T cell activation, differentiation, and cytotoxicity. This study was undertaken to investigate the genetic association of CD226 with susceptibility to systemic lupus erythematosus (SLE) and to assess the functional implications of this association. METHODS: Twelve tag single-nucleotide polymorphisms (SNPs) in CD226 were typed in 1,163 SLE patients and 1,482 healthy control subjects from Europe or of European ancestry. Analyses of association were performed by single-marker Cochran-Mantel-Haenszel meta-analysis, followed by haplotype analysis. Gene expression was analyzed by quantitative real-time polymerase chain reaction analyses of RNA from peripheral blood mononuclear cells, and by fluorescence-activated cell sorter analysis. To study the functional impact of the associated variants, luciferase reporter constructs containing different portions of the 3'-untranslated region (3'-UTR) of the gene were prepared and used in transfection experiments. RESULTS: A 3-variant haplotype, rs763361;rs34794968;rs727088 (ATC), in the last exon of CD226 was associated with SLE (P = 1.3 × 10(-4) , odds ratio 1.24, 95% confidence interval 1.11-1.38). This risk haplotype correlated with low CD226 transcript expression and low CD226 protein levels on the surface of CD4+ and CD8+ T cells and natural killer T (NKT) cells. NK cells expressed high levels of CD226, but this expression was independent of the haplotype. Reporter assays with deletion constructs indicated that only the presence of rs727088 could account for the differences in the levels of luciferase transcripts. CONCLUSION: This study identified an association of CD226 with SLE in individuals of European ancestry. These data support the importance of the 3'-UTR SNP rs727088 in the regulation of CD226 transcription both in T cells and in NKT cells.

Genetic associations in type I interferon related pathways with autoimmunity.

Type I interferons play an outstanding role in innate and adaptive immunity by enhancing functions of dendritic cells, inducing differentiation of monocytes, promoting immunoglobulin class switching in B cells and stimulating effector functions of T cells. The increased production of IFNα/ß by plasmacytoid dendritic cells could be responsible for not only efficient antiviral defence, but it also may be a pathological factor in the development of various autoimmune disorders. The first evidence of a genetic link between type I interferons and autoimmune diseases was the observation that elevated IFNα activity is frequently detected in the sera of patients with systemic lupus erythematosus, and that this trait shows high heritability and familial aggregation in their first-degree healthy relatives. To date, a number of genes involved in interferon signalling have been associated with various autoimmune diseases. Patients with systemic lupus erythematosus, Sjögren's syndrome, dermatomyositis, psoriasis, and a fraction of patients with rheumatoid arthritis display a specific expression pattern of interferon-dependent genes in their leukocytes, termed the interferon signature. Here, in an attempt to understand the role of type I interferons in the pathogenesis of autoimmunity, we review the recent advances in the genetics of autoimmune diseases focusing on the association of genes involved in type I interferon pathways.

Promoter insertion/deletion in the IRF5 gene is highly associated with susceptibility to systemic lupus erythematosus in distinct populations, but exerts a modest effect on gene expression in peripheral blood mononuclear cells.

OBJECTIVE: We examined the genetic association of the promoter insertion/deletion (indel) in IRF5 gene with systemic lupus erythematosus (SLE) in distinct populations and assessed its role in gene expression. METHODS: Four IRF5 polymorphisms were genotyped in 1488 SLE patients and 1466 controls. Gene expression was analyzed by quantitative real-time PCR using RNA from peripheral blood mononuclear cells (PBMC). RESULTS: The promoter indel and rs2070197 had independent genetic effects, which accounted for the association of rs2004640 and rs10954213. Gene expression analysis revealed that rs10954213 exerted the greatest influence on IRF5 transcript levels. CONCLUSION: We corroborated the association of the promoter indel with SLE in 5 different populations and revealed that rs10954213 is the main single-nucleotide polymorphism responsible for altered IRF5 expression in PBMC.

Association of STAT4 and BLK, but not BANK1 or IRF5, with primary antiphospholipid syndrome.

OBJECTIVE: Primary antiphospholipid syndrome (APS) is formally classified by the presence of antiphospholipid antibodies, recurrent thrombosis, and/or pregnancy morbidity in the absence of any underlying full-blown systemic autoimmune disease. However, systemic manifestations in patients with primary APS have been recently reported, as has the presence of serologic markers in common with systemic lupus erythematosus (SLE). In spite of similarities between the 2 diseases, only a minority of cases of primary APS evolve into full-blown SLE, even after a long followup period. The aim of this study was to investigate whether the analysis of SLE susceptibility genes may provide at least a partial explanation for such a discrepancy. METHODS: One hundred thirty-three patients with primary APS classified according to the Sydney criteria and 468 healthy control subjects from the same geographic area were recruited. We genotyped 3 single-nucleotide polymorphisms (SNPs) in IRF5 (rs2004640, rs2070197, and rs10954213), 4 SNPs in STAT4 (rs1467199, rs3821236, rs3024866, and rs7574865), 2 SNPs in BANK1 (rs10516487 and rs3733197), and 1 SNP in BLK (rs2736340). RESULTS: STAT4 and BLK displayed a strong genetic association with primary APS (for rs7574865, odds ratio [OR] 2.19, P=5.17x10(-7); for rs2736340, OR 2.06, P=1.78x10(-6)), while a weak association with IRF5 and no association with BANK1 were observed. CONCLUSION: The presence of a strong genetic association with only a few SLE susceptibility genes and the absence of a more complex gene association may contribute to the lack of cases of full-blown SLE developing in patients with primary APS, in spite of the clinical and serologic similarities between SLE and primary APS.

Kallikrein genes are associated with lupus and glomerular basement membrane-specific antibody-induced nephritis in mice and humans.

Immune-mediated nephritis contributes to disease in systemic lupus erythematosus, Goodpasture syndrome (caused by antibodies specific for glomerular basement membrane [anti-GBM antibodies]), and spontaneous lupus nephritis. Inbred mouse strains differ in susceptibility to anti-GBM antibody-induced and spontaneous lupus nephritis. This study sought to clarify the genetic and molecular factors that maybe responsible for enhanced immune-mediated renal disease in these models. When the kidneys of 3 mouse strains sensitive to anti-GBM antibody-induced nephritis were compared with those of 2 control strains using microarray analysis, one-fifth of the underexpressed genes belonged to the kallikrein gene family,which encodes serine esterases. Mouse strains that upregulated renal and urinary kallikreins exhibited less evidence of disease. Antagonizing the kallikrein pathway augmented disease, while agonists dampened the severity of anti-GBM antibody-induced nephritis. In addition, nephritis-sensitive mouse strains had kallikrein haplotypes that were distinct from those of control strains, including several regulatory polymorphisms,some of which were associated with functional consequences. Indeed, increased susceptibility to anti-GBM antibody-induced nephritis and spontaneous lupus nephritis was achieved by breeding mice with a genetic interval harboring the kallikrein genes onto a disease-resistant background. Finally, both human SLE and spontaneous lupus nephritis were found to be associated with kallikrein genes, particularly KLK1 and the KLK3 promoter, when DNA SNPs from independent cohorts of SLE patients and controls were compared. Collectively, these studies suggest that kallikreins are protective disease-associated genes in anti-GBM antibody-induced nephritis and lupus.

Replication of the TNFSF4 (OX40L) promoter region association with systemic lupus erythematosus.

The tumor necrosis factor ligand superfamily member 4 gene (TNFSF4) encodes the OX40 ligand (OX40L), a costimulatory molecule involved in T-cell activation. A recent study demonstrated the association of TNFSF4 haplotypes located in the upstream region with risk for or protection from systemic lupus erythematosus (SLE). To replicate this association, five single nucleotide polymorphisms (SNPs) tagging the previously associated haplotypes and passing the proper quality-control filters were tested in 1312 cases and 1801 controls from Germany, Italy, Spain and Argentina. The association of TNFSF4 with SLE was replicated in all the sets except Spain. There was a unique risk haplotype tagged by the minor alleles of the SNPs rs1234317 (pooled odds ratio (OR)=1.39, P=0.0009) and rs12039904 (pooled OR=1.38, P=0.0012). We did not observe association to a single protective marker (rs844644) or haplotype as the first study reported; instead, we observed different protective haplotypes, all carrying the major alleles of both SNPs rs1234317 and rs12039904. Association analysis conditioning on the haplotypic background confirmed that these two SNPs explain the entire haplotype effect. This first replication study confirms the association of genetic variation in the upstream region of TNFSF4 with susceptibility to SLE.

Association of STAT4 with rheumatoid arthritis: a replication study in three European populations.

OBJECTIVE: This study was undertaken to investigate the previously reported association of the STAT4 polymorphism rs7574865 with rheumatoid arthritis (RA) in 3 different European populations from Spain, Sweden, and The Netherlands, comprising a total of 2,072 patients and 2,474 controls. METHODS: Three different cohorts were included in the study: 923 RA patients and 1,296 healthy controls from Spain, 273 RA patients and 285 healthy controls from Sweden, and 876 RA patients and 893 healthy controls from The Netherlands. DNA from patients and controls was obtained from peripheral blood. Samples were genotyped for the STAT4 single-nucleotide polymorphism rs7574865 using a TaqMan 5'-allele discrimination assay. The chi-square test was performed to compare allele and genotype distributions. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated. RESULTS: We observed a significantly increased frequency of the minor T allele in RA patients compared with healthy controls in the Spanish population (24.8% versus 20.8%; P = 0.001, OR 1.26 [95% CI 1.09-1.45]). This association was confirmed in both the Swedish population (P = 0.03, OR 1.35 [95% CI 1.03-1.77]) and the Dutch population (P = 0.03, OR 1.45 [95% CI 1.21-1.73]). The overall P value for all 3 populations was 9.79 x 10(-6) (OR 1.25 [95% CI 1.13-1.37]). No association between rs7574865 and the presence of rheumatoid factor or anti-cyclic citrullinated peptide autoantibodies was observed. A meta-analysis of all published STAT4 associations revealed an OR of 1.25 (95% CI 1.19-1.33) (P = 1 x 10(-5)). CONCLUSION: Our findings indicate an association between the STAT4 polymorphism rs7574865 and RA in 3 different populations, from Spain, Sweden, and The Netherlands, thereby confirming previous data.

Meta-analysis of HLA-DRB1 polymorphism in Latin American patients with rheumatoid arthritis.

OBJECTIVES: To estimate the common effect size of HLA-DRB1 alleles on rheumatoid arthritis (RA) susceptibility across Latin America populations through a meta-analysis combining the results of published data. METHODS: Case-control studies on HLA-DRB1 association with RA in Latin America were searched up to October 2006. Genotype frequencies were extracted according to both shared epitope (SE) and HLA-DR4 positive or negative alleles. The effect summary odds ratio (OR) and 95% confidence intervals was obtained. Heterogeneity and publication bias were assessed. RESULTS: Eight studies containing 684 cases and 1015 controls were included. Under the random effects model, the common OR was 3.28 (1.93, 5.60) (p<0.0001) and 3.54 (2.47, 5.05) (p=4.22 x 10(-12)) for HLA-DR4 and SE, respectively. There was no evidence of publication bias according to Funnel plot and Egger's regression test (p=0,445 for DR4 and p=0,464 for SE meta-analysis). Significant heterogeneity was observed for HLA-DR4 (I2=81.06%, Q=36.96, p=0.000005) but not for the SE meta-analysis. CONCLUSIONS: HLA-DR4 and SE positive HLA-DRB1 alleles (mainly HLA-DRB10404) are associated with RA in Latin Americans. Heterogeneity is expected owing to the diverse degree of admixture between the examined populations. Our findings support the HLA as a major susceptibility locus for RA and validate the SE hypothesis in Latin America.