Elevated levels of IL-12/IL-23p40 in Nova Scotia Duck Tolling Retrievers with autoimmune disease and lymphoma.

The Nova Scotia Duck Tolling Retriever (NSDTR) is predisposed to immune mediated rheumatic disease (IMRD), steroid-responsive meningitis-arteritis (SRMA) and certain forms of cancer. Cytokines are the main regulators of the immune system. Interleukin 2 is a cytokine involved in activation of T regulatory cells, playing a role in central tolerance and tumor immunity. Interleukin 12 and interleukin 23 share the same subunit, p40, and are both pro-inflammatory cytokines. The aim of this study was to compare levels of IL-2 in healthy NSDTRs to those with cancer or autoimmune disease and to compare levels of IL-12/IL-23p40 in healthy NSDTRs and beagles versus NSDTRs with cancer or autoimmune disease. 62 dogs were included in the analysis of IL-12/IL-23p40; healthy NSDTRs (n = 16), healthy beagles (n = 16), NSDTRs autoimmune (n = 18) and NDSTRs lymphoma/mastocytoma (n = 12) and 68 dogs for IL-2; healthy (n = 20), autoimmune (n = 36) and lymphoma/mastocytoma/adenocarcinoma (n = 12). NSDTRs with autoimmune disease had higher levels of IL-12/IL-23p40 compared to healthy dogs (p = 0.008). NSDTRs with lymphoma also had higher levels of IL-12/IL-23p40 compared to healthy NSDTRs (p = 0.002). There was no difference in levels of IL-2 between healthy and diseased NSDTR. Statistical analysis was performed using Bonferroni corrections for multiple testing. These findings can contribute to the knowledge of autoimmune disease and cancer in dogs.

A survey of ficolin-3 activity in Systemic Lupus Erythematosus reveals a link to hematological disease manifestations and autoantibody profile.

The complement system plays a central role in the pathogenesis of Systemic Lupus Erythematosus (SLE), but most studies have focused on the classical pathway. Ficolin-3 is the main initiator of the lectin pathway of complement in humans, but its role in systemic autoimmune disease has not been conclusively determined. Here, we combined biochemical and genetic approaches to assess the contribution of ficolin-3 to SLE risk and disease manifestations. Ficolin-3 activity was measured by a functional assay in serum or plasma samples from Swedish SLE patients (n = 786) and controls matched for age and sex (n = 566). Genetic variants in an extended 300 kb genomic region spanning the FCN3 locus were analyzed for their association with ficolin-3 activity and SLE manifestations in a Swedish multicenter cohort (n = 985). Patients with ficolin-3 activity in the highest tertile showed a strong enrichment in an SLE cluster defined by anti-Sm/DNA/nucleosome antibodies (OR 3.0, p < 0.001) and had increased rates of hematological disease (OR 1.4, p = 0.078) and lymphopenia (OR = 1.6, p = 0.039). Genetic variants associated with low ficolin-3 activity mapped to an extended haplotype in high linkage disequilibrium upstream of the FCN3 gene. Patients carrying the lead genetic variant associated with low ficolin-3 activity had a lower frequency of hematological disease (OR 0.67, p = 0.018) and lymphopenia (OR 0.63, p = 0.031) and fewer autoantibodies (p = 0.0019). Loss-of-function variants in the FCN3 gene were not associated with SLE, but four (0.5 %) SLE patients developed acquired ficolin-3 deficiency where ficolin-3 activity in serum was depleted following diagnosis of SLE. Taken together, our results provide genetic and biochemical evidence that implicate the lectin pathway in hematological SLE manifestations. We also identify lectin pathway activation through ficolin-3 as a factor that contributes to the autoantibody response in SLE.

Using evolutionary constraint to define novel candidate driver genes in medulloblastoma.

Current knowledge of cancer genomics remains biased against noncoding mutations. To systematically search for regulatory noncoding mutations, we assessed mutations in conserved positions in the genome under the assumption that these are more likely to be functional than mutations in positions with low conservation. To this end, we use whole-genome sequencing data from the International Cancer Genome Consortium and combined it with evolutionary constraint inferred from 240 mammals, to identify genes enriched in noncoding constraint mutations (NCCMs), mutations likely to be regulatory in nature. We compare medulloblastoma (MB), which is malignant, to pilocytic astrocytoma (PA), a primarily benign tumor, and find highly different NCCM frequencies between the two, in agreement with the fact that malignant cancers tend to have more mutations. In PA, a high NCCM frequency only affects the BRAF locus, which is the most commonly mutated gene in PA. In contrast, in MB, >500 genes have high levels of NCCMs. Intriguingly, several loci with NCCMs in MB are associated with different ages of onset, such as the HOXB cluster in young MB patients. In adult patients, NCCMs occurred in, e.g., the WASF-2/AHDC1/FGR locus. One of these NCCMs led to increased expression of the SRC kinase FGR and augmented responsiveness of MB cells to dasatinib, a SRC kinase inhibitor. Our analysis thus points to different molecular pathways in different patient groups. These newly identified putative candidate driver mutations may aid in patient stratification in MB and could be valuable for future selection of personalized treatment options.

Evolutionary constraint and innovation across hundreds of placental mammals.

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

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.

Strong Association of Combined Genetic Deficiencies in the Classical Complement Pathway With Risk of Systemic Lupus Erythematosus and Primary Sjögren's Syndrome.

OBJECTIVE: Complete genetic deficiency of the complement component C2 is a strong risk factor for monogenic systemic lupus erythematosus (SLE), but whether heterozygous C2 deficiency adds to the risk of SLE or primary Sjögren's syndrome (SS) has not been studied systematically. This study was undertaken to investigate potential associations of heterozygous C2 deficiency and C4 copy number variation with clinical manifestations in patients with SLE and patients with primary SS. METHODS: The presence of the common 28-bp C2 deletion rs9332736 and C4 copy number variation was examined in Scandinavian patients who had received a diagnosis of SLE (n = 958) or primary SS (n = 911) and in 2,262 healthy controls through the use of DNA sequencing. The concentration of complement proteins in plasma and classical complement function were analyzed in a subgroup of SLE patients. RESULTS: Heterozygous C2 deficiency-when present in combination with a low C4A copy number-substantially increased the risk of SLE (odds ratio [OR] 10.2 [95% confidence interval (95% CI) 3.5-37.0]) and the risk of primary SS (OR 13.0 [95% CI 4.5-48.4]) when compared to individuals with 2 C4A copies and normal C2. For patients heterozygous for rs9332736 with 1 C4A copy, the median age at diagnosis was 7 years earlier in patients with SLE and 12 years earlier in patients with primary SS when compared to patients with normal C2. Reduced C2 levels in plasma (P = 2 × 10-9 ) and impaired function of the classical complement pathway (P = 0.03) were detected in SLE patients with heterozygous C2 deficiency. Finally, in a primary SS patient homozygous for C2 deficiency, we observed low levels of anti-Scl-70, which suggests a risk of developing systemic sclerosis or potential overlap between primary SS and other systemic autoimmune diseases. CONCLUSION: We demonstrate that a genetic pattern involving partial deficiencies of C2 and C4A in the classical complement pathway is a strong risk factor for SLE and for primary SS. Our results emphasize the central role of the complement system in the pathogenesis of both SLE and primary SS.

Complement C4 Copy Number Variation is Linked to SSA/Ro and SSB/La Autoantibodies in Systemic Inflammatory Autoimmune Diseases.

OBJECTIVE: Copy number variation of the C4 complement components, C4A and C4B, has been associated with systemic inflammatory autoimmune diseases. This study was undertaken to investigate whether C4 copy number variation is connected to the autoimmune repertoire in systemic lupus erythematosus (SLE), primary Sjögren's syndrome (SS), or myositis. METHODS: Using targeted DNA sequencing, we determined the copy number and genetic variants of C4 in 2,290 well-characterized Scandinavian patients with SLE, primary SS, or myositis and 1,251 healthy controls. RESULTS: A prominent relationship was observed between C4A copy number and the presence of SSA/SSB autoantibodies, which was shared between the 3 diseases. The strongest association was detected in patients with autoantibodies against both SSA and SSB and 0 C4A copies when compared to healthy controls (odds ratio [OR] 18.0 [95% confidence interval (95% CI) 10.2-33.3]), whereas a weaker association was seen in patients without SSA/SSB autoantibodies (OR 3.1 [95% CI 1.7-5.5]). The copy number of C4 correlated positively with C4 plasma levels. Further, a common loss-of-function variant in C4A leading to reduced plasma C4 was more prevalent in SLE patients with a low copy number of C4A. Functionally, we showed that absence of C4A reduced the individuals' capacity to deposit C4b on immune complexes. CONCLUSION: We show that a low C4A copy number is more strongly associated with the autoantibody repertoire than with the clinically defined disease entities. These findings may have implications for understanding the etiopathogenetic mechanisms of systemic inflammatory autoimmune diseases and for patient stratification when taking the genetic profile into account.

Contribution of Rare Genetic Variation to Disease Susceptibility in a Large Scandinavian Myositis Cohort.

OBJECTIVE: Idiopathic inflammatory myopathies (IIMs) are a heterogeneous group of complex autoimmune conditions characterized by inflammation in skeletal muscle and extramuscular compartments, and interferon (IFN) system activation. We undertook this study to examine the contribution of genetic variation to disease susceptibility and to identify novel avenues for research in IIMs. METHODS: Targeted DNA sequencing was used to mine coding and potentially regulatory single nucleotide variants from ~1,900 immune-related genes in a Scandinavian case-control cohort of 454 IIM patients and 1,024 healthy controls. Gene-based aggregate testing, together with rare variant- and gene-level enrichment analyses, was implemented to explore genotype-phenotype relations. RESULTS: Gene-based aggregate tests of all variants, including rare variants, identified IFI35 as a potential genetic risk locus for IIMs, suggesting a genetic signature of type I IFN pathway activation. Functional annotation of the IFI35 locus highlighted a regulatory network linked to the skeletal muscle-specific gene PTGES3L, as a potential candidate for IIM pathogenesis. Aggregate genetic associations with AGER and PSMB8 in the major histocompatibility complex locus were detected in the antisynthetase syndrome subgroup, which also showed a less marked genetic signature of the type I IFN pathway. Enrichment analyses indicated a burden of synonymous and noncoding rare variants in IIM patients, suggesting increased disease predisposition associated with these classes of rare variants. CONCLUSION: Our study suggests the contribution of rare genetic variation to disease susceptibility in IIM and specific patient subgroups, and pinpoints genetic associations consistent with previous findings by gene expression profiling. These features highlight genetic profiles that are potentially relevant to disease pathogenesis.

Identification and functional characterization of a novel susceptibility locus for small vessel vasculitis with MPO-ANCA.

OBJECTIVE: To identify and characterize genetic loci associated with the risk of developing ANCA-associated vasculitides (AAV). METHODS: Genetic association analyses were performed after Illumina sequencing of 1853 genes and subsequent replication with genotyping of selected single nucleotide polymorphisms in a total cohort of 1110 Scandinavian cases with granulomatosis with polyangiitis or microscopic polyangiitis, and 1589 controls. A novel AAV-associated single nucleotide polymorphism was analysed for allele-specific effects on gene expression using luciferase reporter assay. RESULTS: PR3-ANCA+ AAV was significantly associated with two independent loci in the HLA-DPB1/HLA-DPA1 region [rs1042335, P = 6.3 × 10-61, odds ratio (OR) 0.10; rs9277341, P = 1.5 × 10-44, OR 0.22] and with rs28929474 in the SERPINA1 gene (P = 2.7 × 10-10, OR 2.9). MPO-ANCA+ AAV was significantly associated with the HLA-DQB1/HLA-DQA2 locus (rs9274619, P = 5.4 × 10-25, OR 3.7) and with a rare variant in the BACH2 gene (rs78275221, P = 7.9 × 10-7, OR 3.0), the latter a novel susceptibility locus for MPO-ANCA+ granulomatosis with polyangiitis/microscopic polyangiitis. The rs78275221-A risk allele reduced luciferase gene expression in endothelial cells, specifically, as compared with the non-risk allele. CONCLUSION: We identified a novel susceptibility locus for MPO-ANCA+ AAV and propose that the associated variant is of mechanistic importance, exerting a regulatory function on gene expression in specific cell types.

Hierarchical Structure of Protein Sequence.

Most non-communicable diseases are associated with dysfunction of proteins or protein complexes. The relationship between sequence and structure has been analyzed for a long time, and the analysis of the sequences organization in domains and motifs remains an actual research area. Here, we propose a mathematical method for revealing the hierarchical organization of protein sequences. The method is based on the pentapeptide as a unit of protein sequences. Employing the frequency of occurrence of pentapeptides in sequences of natural proteins and a special mathematical approach, this method revealed a hierarchical structure in the protein sequence. The method was applied to 24,647 non-homologous protein sequences with sizes ranging from 50 to 400 residues from the NRDB90 database. Statistical analysis of the branching points of the graphs revealed 11 characteristic values of y (the width of the inscribed function), showing the relationship of these multiple fragments of the sequences. Several examples illustrate how fragments of the protein spatial structure correspond to the elements of the hierarchical structure of the protein sequence. This methodology provides a promising basis for a mathematically-based classification of the elements of the spatial organization of proteins. Elements of the hierarchical structure of different levels of the hierarchy can be used to solve biotechnological and medical problems.

Molecular pathways in patients with systemic lupus erythematosus revealed by gene-centred DNA sequencing.

OBJECTIVES: Systemic lupus erythematosus (SLE) is an autoimmune disease with extensive heterogeneity in disease presentation between patients, which is likely due to an underlying molecular diversity. Here, we aimed at elucidating the genetic aetiology of SLE from the immunity pathway level to the single variant level, and stratify patients with SLE into distinguishable molecular subgroups, which could inform treatment choices in SLE. METHODS: We undertook a pathway-centred approach, using sequencing of immunological pathway genes. Altogether 1832 candidate genes were analysed in 958 Swedish patients with SLE and 1026 healthy individuals. Aggregate and single variant association testing was performed, and we generated pathway polygenic risk scores (PRS). RESULTS: We identified two main independent pathways involved in SLE susceptibility: T lymphocyte differentiation and innate immunity, characterised by HLA and interferon, respectively. Pathway PRS defined pathways in individual patients, who on average were positive for seven pathways. We found that SLE organ damage was more pronounced in patients positive for the T or B cell receptor signalling pathways. Further, pathway PRS-based clustering allowed stratification of patients into four groups with different risk score profiles. Studying sets of genes with priors for involvement in SLE, we observed an aggregate common variant contribution to SLE at genes previously reported for monogenic SLE as well as at interferonopathy genes. CONCLUSIONS: Our results show that pathway risk scores have the potential to stratify patients with SLE beyond clinical manifestations into molecular subsets, which may have implications for clinical follow-up and therapy selection.

A rare regulatory variant in the MEF2D gene affects gene regulation and splicing and is associated with a SLE sub-phenotype in Swedish cohorts.

Systemic lupus erythematosus (SLE) is an autoimmune disorder with heterogeneous clinical presentation and complex etiology involving the interplay between genetic, epigenetic, environmental and hormonal factors. Many common SNPs identified by genome wide-association studies (GWAS) explain only a small part of the disease heritability suggesting the contribution from rare genetic variants, undetectable in GWAS, and complex epistatic interactions. Using targeted re-sequencing of coding and conserved regulatory regions within and around 215 candidate genes selected on the basis of their known role in autoimmunity and genes associated with canine immune-mediated diseases, we identified a rare regulatory variant rs200395694:G > T located in intron 4 of the MEF2D gene encoding the myocyte-specific enhancer factor 2D transcription factor and associated with SLE in Swedish cohorts (504 SLE patients and 839 healthy controls, p = 0.014, CI = 1.1-10). Fisher's exact test revealed an association between the genetic variant and a triad of disease manifestations including Raynaud, anti-U1-ribonucleoprotein (anti-RNP), and anti-Smith (anti-Sm) antibodies (p = 0.00037) among the patients. The DNA-binding activity of the allele was further studied by EMSA, reporter assays, and minigenes. The region has properties of an active cell-specific enhancer, differentially affected by the alleles of rs200395694:G > T. In addition, the risk allele exerts an inhibitory effect on the splicing of the alternative tissue-specific isoform, and thus may modify the target gene set regulated by this isoform. These findings emphasize the potential of dissecting traits of complex diseases and correlating them with rare risk alleles with strong biological effects.

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.

Variants within the SP110 nuclear body protein modify risk of canine degenerative myelopathy.

Canine degenerative myelopathy (DM) is a naturally occurring neurodegenerative disease with similarities to some forms of amyotrophic lateral sclerosis (ALS). Most dogs that develop DM are homozygous for a common superoxide dismutase 1 gene (SOD1) mutation. However, not all dogs homozygous for this mutation develop disease. We performed a genome-wide association analysis in the Pembroke Welsh Corgi (PWC) breed comparing DM-affected and -unaffected dogs homozygous for the SOD1 mutation. The analysis revealed a modifier locus on canine chromosome 25. A haplotype within the SP110 nuclear body protein (SP110) was present in 40% of affected compared with 4% of unaffected dogs (P = 1.5 × 10(-5)), and was associated with increased probability of developing DM (P = 4.8 × 10(-6)) and earlier onset of disease (P = 1.7 × 10(-5)). SP110 is a nuclear body protein involved in the regulation of gene transcription. Our findings suggest that variations in SP110-mediated gene transcription may underlie, at least in part, the variability in risk for developing DM among PWCs that are homozygous for the disease-related SOD1 mutation. Further studies are warranted to clarify the effect of this modifier across dog breeds.

Multiple Changes of Gene Expression and Function Reveal Genomic and Phenotypic Complexity in SLE-like Disease.

The complexity of clinical manifestations commonly observed in autoimmune disorders poses a major challenge to genetic studies of such diseases. Systemic lupus erythematosus (SLE) affects humans as well as other mammals, and is characterized by the presence of antinuclear antibodies (ANA) in patients' sera and multiple disparate clinical features. Here we present evidence that particular sub-phenotypes of canine SLE-related disease, based on homogenous (ANA(H)) and speckled ANA (ANA(S)) staining pattern, and also steroid-responsive meningitis-arteritis (SRMA) are associated with different but overlapping sets of genes. In addition to association to certain MHC alleles and haplotypes, we identified 11 genes (WFDC3, HOMER2, VRK1, PTPN3, WHAMM, BANK1, AP3B2, DAPP1, LAMTOR3, DDIT4L and PPP3CA) located on five chromosomes that contain multiple risk haplotypes correlated with gene expression and disease sub-phenotypes in an intricate manner. Intriguingly, the association of BANK1 with both human and canine SLE appears to lead to similar changes in gene expression levels in both species. Our results suggest that molecular definition may help unravel the mechanisms of different clinical features common between and specific to various autoimmune disease phenotypes in dogs and humans.

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.