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

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

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

OBJECTIVES: Systemic lupus erythematosus (SLE; OMIM 152700) is characterised by the production of antibodies to nuclear antigens. We previously identified variants in complement receptor 2 (CR2/CD21) that were associated with decreased risk of SLE. This study aimed to identify the causal variant for this association. METHODS: Genotyped and imputed genetic variants spanning CR2 were assessed for association with SLE in 15 750 case-control subjects from four ancestral groups. Allele-specific functional effects of associated variants were determined using quantitative real-time PCR, quantitative flow cytometry, electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP)-PCR. RESULTS: The strongest association signal was detected at rs1876453 in intron 1 of CR2 (pmeta=4.2×10(-4), OR 0.85), specifically when subjects were stratified based on the presence of dsDNA autoantibodies (case-control pmeta=7.6×10(-7), OR 0.71; case-only pmeta=1.9×10(-4), OR 0.75). Although allele-specific effects on B cell CR2 mRNA or protein levels were not identified, levels of complement receptor 1 (CR1/CD35) mRNA and protein were significantly higher on B cells of subjects harbouring the minor allele (p=0.0248 and p=0.0006, respectively). The minor allele altered the formation of several DNA protein complexes by EMSA, including one containing CCCTC-binding factor (CTCF), an effect that was confirmed by ChIP-PCR. CONCLUSIONS: These data suggest that rs1876453 in CR2 has long-range effects on gene regulation that decrease susceptibility to lupus. Since the minor allele at rs1876453 is preferentially associated with reduced risk of the highly specific dsDNA autoantibodies that are present in preclinical, active and severe lupus, understanding its mechanisms will have important therapeutic implications.

Lupus Risk Variant Increases pSTAT1 Binding and Decreases ETS1 Expression.

Genetic variants at chromosomal region 11q23.3, near the gene ETS1, have been associated with systemic lupus erythematosus (SLE), or lupus, in independent cohorts of Asian ancestry. Several recent studies have implicated ETS1 as a critical driver of immune cell function and differentiation, and mice deficient in ETS1 develop an SLE-like autoimmunity. We performed a fine-mapping study of 14,551 subjects from multi-ancestral cohorts by starting with genotyped variants and imputing to all common variants spanning ETS1. By constructing genetic models via frequentist and Bayesian association methods, we identified 16 variants that are statistically likely to be causal. We functionally assessed each of these variants on the basis of their likelihood of affecting transcription factor binding, miRNA binding, or chromatin state. Of the four variants that we experimentally examined, only rs6590330 differentially binds lysate from B cells. Using mass spectrometry, we found more binding of the transcription factor signal transducer and activator of transcription 1 (STAT1) to DNA near the risk allele of rs6590330 than near the non-risk allele. Immunoblot analysis and chromatin immunoprecipitation of pSTAT1 in B cells heterozygous for rs6590330 confirmed that the risk allele increased binding to the active form of STAT1. Analysis with expression quantitative trait loci indicated that the risk allele of rs6590330 is associated with decreased ETS1 expression in Han Chinese, but not other ancestral cohorts. We propose a model in which the risk allele of rs6590330 is associated with decreased ETS1 expression and increases SLE risk by enhancing the binding of pSTAT1.

Association of two independent functional risk haplotypes in TNIP1 with systemic lupus erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified >30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting NF-κB signaling. The present study was undertaken to investigate the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway. METHODS: We analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog TNIP2, in case-control populations of diverse ethnic origins. TNIP1, TNIP2, and TAX1BP1 were fine-mapped in a total of 8,372 SLE cases and 7,492 healthy controls from European-ancestry, African American, Hispanic, East Asian, and African American Gullah populations. Levels of TNIP1 messenger RNA (mRNA) and ABIN1 protein in Epstein-Barr virus-transformed human B cell lines were analyzed by quantitative reverse transcription-polymerase chain reaction and Western blotting, respectively. RESULTS: We found significant associations between SLE and genetic variants within TNIP1, but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified 2 independent risk haplotypes within TNIP1 in individuals of European ancestry that were also present in African American and Hispanic populations. Levels of TNIP1 mRNA and ABIN1 protein were reduced among subjects with these haplotypes, suggesting that they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression. CONCLUSION: Our results confirm the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis.

Association of genetic variants in complement factor H and factor H-related genes with systemic lupus erythematosus susceptibility.

Systemic lupus erythematosus (SLE), a complex polygenic autoimmune disease, is associated with increased complement activation. Variants of genes encoding complement regulator factor H (CFH) and five CFH-related proteins (CFHR1-CFHR5) within the chromosome 1q32 locus linked to SLE, have been associated with multiple human diseases and may contribute to dysregulated complement activation predisposing to SLE. We assessed 60 SNPs covering the CFH-CFHRs region for association with SLE in 15,864 case-control subjects derived from four ethnic groups. Significant allelic associations with SLE were detected in European Americans (EA) and African Americans (AA), which could be attributed to an intronic CFH SNP (rs6677604, in intron 11, P(meta) = 6.6×10(-8), OR = 1.18) and an intergenic SNP between CFHR1 and CFHR4 (rs16840639, P(meta) = 2.9×10(-7), OR = 1.17) rather than to previously identified disease-associated CFH exonic SNPs, including I62V, Y402H, A474A, and D936E. In addition, allelic association of rs6677604 with SLE was subsequently confirmed in Asians (AS). Haplotype analysis revealed that the underlying causal variant, tagged by rs6677604 and rs16840639, was localized to a ~146 kb block extending from intron 9 of CFH to downstream of CFHR1. Within this block, the deletion of CFHR3 and CFHR1 (CFHR3-1Δ), a likely causal variant measured using multiplex ligation-dependent probe amplification, was tagged by rs6677604 in EA and AS and rs16840639 in AA, respectively. Deduced from genotypic associations of tag SNPs in EA, AA, and AS, homozygous deletion of CFHR3-1Δ (P(meta) = 3.2×10(-7), OR = 1.47) conferred a higher risk of SLE than heterozygous deletion (P(meta) = 3.5×10(-4), OR = 1.14). These results suggested that the CFHR3-1Δ deletion within the SLE-associated block, but not the previously described exonic SNPs of CFH, might contribute to the development of SLE in EA, AA, and AS, providing new insights into the role of complement regulators in the pathogenesis of SLE.