COVID-19 and systemic lupus erythematosus genetics: A balance between autoimmune disease risk and protection against infection

Genome wide association studies show there is a genetic component to severe COVID-19. We find evidence that the genome-wide genetic association signal with severe COVID-19 is correlated with that of systemic lupus erythematosus (SLE), having formally tested this using genetic correlation analysis by LD score regression. To identify the shared associated loci and gain insight into the shared genetic effects, using summary level data we performed meta-analyses, a local genetic correlation analysis and fine-mapping using stepwise regression and functional annotation. This identified multiple loci shared between the two traits, some of which exert opposing effects. The locus with most evidence of shared association is TYK2, a gene critical to the type I interferon pathway, where the local genetic correlation is negative. Another shared locus is CLEC1A, where the direction of effects is aligned, that encodes a lectin involved in cell signaling, and the anti-fungal immune response. Our analyses suggest that several loci with reciprocal effects between the two traits have a role in the defense response pathway, adding to the evidence that SLE risk alleles are protective against infection.

We observed a correlation between the genetic associations with severe COVID-19 and those with systemic lupus erythematosus (SLE, Lupus), and aimed to discover which genetic loci were shared by these diseases and what biological processes were involved. This resulted in the discovery of several genetic loci, some of which had alleles that were risk for both diseases and some of which were risk for severe COVID-19 yet protective for SLE. The locus with most evidence of shared association (TYK2) is involved in interferon production, a process that is important in response to viral infection and known to be dysregulated in SLE patients. Other shared associated loci contained genes also involved in the defense response and the immune system signaling. These results add to the growing evidence that there are alleles in the human genome that provide protection against viral infection yet are risk for autoimmune disease.

Genome-wide assessment of genetic risk for systemic lupus erythematosus and disease severity

Using three European and two Chinese genome-wide association studies (GWAS), we investigated the performance of genetic risk scores (GRSs) for predicting the susceptibility and severity of systemic lupus erythematosus (SLE), using renal disease as a proxy for severity. We used four GWASs to test the performance of GRS both cross validating within the European population and between European and Chinese populations. The performance of GRS in SLE risk prediction was evaluated by receiver operating characteristic (ROC) curves. We then analyzed the polygenic nature of SLE statistically. We also partitioned patients according to their age-of-onset and evaluated the predictability of GRS in disease severity in each age group. We found consistently that the best GRS in the prediction of SLE used SNPs associated at the level of P < 1e−05 in all GWAS data sets and that SNPs with P-values above 0.2 were inflated for SLE true positive signals. The GRS results in an area under the ROC curve ranging between 0.64 and 0.72, within European and between the European and Chinese populations. We further showed a significant positive correlation between a GRS and renal disease in two independent European GWAS (P_cohort1 = 2.44e−08; P_cohort2 = 0.00205) and a significant negative correlation with age of SLE onset (P_cohort1 = 1.76e−12; P_cohort2 = 0.00384). We found that the GRS performed better in the prediction of renal disease in the ‘later onset’ compared with the ‘earlier onset’ group. The GRS predicts SLE in both European and Chinese populations and correlates with poorer prognostic factors: young age-of-onset and lupus nephritis.

Nucleolin acts as the receptor for C1QTNF4 and supports C1QTNF4-mediated innate immunity modulation

The C1q and TNF related 4 (C1QTNF4) protein is a structurally unique member of the C1QTNF family, a family of secreted proteins that have structural homology with both complement C1q and the tumor necrosis factor superfamily. C1QTNF4 has been linked to the autoimmune disease systemic lupus erythematosus through genetic studies; however, its role in immunity and inflammation remains poorly defined and a cell surface receptor of C1QTNF4 has yet to be identified. Here we report identification of nucleolin as a cell surface receptor of C1QTNF4 using mass spectrometric analysis. Additionally, we present evidence that the interaction between C1QTNF4 and nucleolin is mediated by the second C1q-like domain of C1QTNF4 and the C terminus of nucleolin. We show that monocytes and B cells are target cells of C1QTNF4 and observe extensive binding to dead cells. Imaging flow cytometry experiments in monocytes show that C1QTNF4 becomes actively internalized upon cell binding. Our results suggest that nucleolin may serve as a docking molecule for C1QTNF4 and act in a context-dependent manner through coreceptors. Taken together, these findings further our understanding of C1QTNF4's function in the healthy immune system and how dysfunction may contribute to the development of systemic lupus erythematosus.

Trans-Ancestral Fine-Mapping and Epigenetic Annotation as Tools to Delineate Functionally Relevant Risk Alleles at IKZF1 and IKZF3 in Systemic Lupus Erythematosus

Background: Prioritizing tag-SNPs carried on extended risk haplotypes at susceptibility loci for common disease is a challenge. Methods: We utilized trans-ancestral exclusion mapping to reduce risk haplotypes at IKZF1 and IKZF3 identified in multiple ancestries from SLE GWAS and ImmunoChip datasets. We characterized functional annotation data across each risk haplotype from publicly available datasets including ENCODE, RoadMap Consortium, PC Hi-C data from 3D genome browser, NESDR NTR conditional eQTL database, GeneCards Genehancers and TF (transcription factor) binding sites from Haploregv4. Results: We refined the 60 kb associated haplotype upstream of IKZF1 to just 12 tag-SNPs tagging a 47.7 kb core risk haplotype. There was preferential enrichment of DNAse I hypersensitivity and H3K27ac modification across the 3′ end of the risk haplotype, with four tag-SNPs sharing allele-specific TF binding sites with promoter variants, which are eQTLs for IKZF1 in whole blood. At IKZF3, we refined a core risk haplotype of 101 kb (27 tag-SNPs) from an initial extended haplotype of 194 kb (282 tag-SNPs), which had widespread DNAse I hypersensitivity, H3K27ac modification and multiple allele-specific TF binding sites. Dimerization of Fox family TFs bound at the 3′ and promoter of IKZF3 may stabilize chromatin looping across the locus. Conclusions: We combined trans-ancestral exclusion mapping and epigenetic annotation to identify variants at both IKZF1 and IKZF3 with the highest likelihood of biological relevance. The approach will be of strong interest to other complex trait geneticists seeking to attribute biological relevance to risk alleles on extended risk haplotypes in their disease of interest.

Interferon inducible X-linked gene CXorf21 may contribute to sexual dimorphism in Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is an autoimmune disease, characterised by increased expression of type I interferon (IFN)-regulated genes and a striking sex imbalance towards females. Through combined genetic, in silico, in vitro, and ex vivo approaches, we define CXorf21, a gene of hitherto unknown function, which escapes X-chromosome inactivation, as a candidate underlying the Xp21.2 SLE association. We demonstrate that CXorf21 is an IFN-response gene and that the sexual dimorphism in expression is magnified by immunological challenge. Fine-mapping reveals a single haplotype as a potential causal cis-eQTL for CXorf21. We propose that expression is amplified through modification of promoter and 3'-UTR chromatin interactions. Finally, we show that the CXORF21 protein colocalises with TLR7, a pathway implicated in SLE pathogenesis. Our study reveals modulation in gene expression affected by the combination of two hallmarks of SLE: CXorf21 expression increases in a both an IFN-inducible and sex-specific manner.

Profiling RNA-Seq at multiple resolutions markedly increases the number of causal eQTLs in autoimmune disease

Genome-wide association studies have identified hundreds of risk loci for autoimmune disease, yet only a minority (~25%) share genetic effects with changes to gene expression (eQTLs) in immune cells. RNA-Seq based quantification at whole-gene resolution, where abundance is estimated by culminating expression of all transcripts or exons of the same gene, is likely to account for this observed lack of colocalisation as subtle isoform switches and expression variation in independent exons can be concealed. We performed integrative cis-eQTL analysis using association statistics from twenty autoimmune diseases (560 independent loci) and RNA-Seq data from 373 individuals of the Geuvadis cohort profiled at gene-, isoform-, exon-, junction-, and intron-level resolution in lymphoblastoid cell lines. After stringently testing for a shared causal variant using both the Joint Likelihood Mapping and Regulatory Trait Concordance frameworks, we found that gene-level quantification significantly underestimated the number of causal cis-eQTLs. Only 5.0-5.3% of loci were found to share a causal cis-eQTL at gene-level compared to 12.9-18.4% at exon-level and 9.6-10.5% at junction-level. More than a fifth of autoimmune loci shared an underlying causal variant in a single cell type by combining all five quantification types; a marked increase over current estimates of steady-state causal cis-eQTLs. Causal cis-eQTLs detected at different quantification types localised to discrete epigenetic annotations. We applied a linear mixed-effects model to distinguish cis-eQTLs modulating all expression elements of a gene from those where the signal is only evident in a subset of elements. Exon-level analysis detected disease-associated cis-eQTLs that subtly altered transcription globally across the target gene. We dissected in detail the genetic associations of systemic lupus erythematosus and functionally annotated the candidate genes. Many of the known and novel genes were concealed at gene-level (e.g. IKZF2, TYK2, LYST). Our findings are provided as a web resource.

B cell OX40L supports T follicular helper cell development and contributes to SLE pathogenesis

Objectives

TNFSF4 (encodes OX40L) is a susceptibility locus for systemic lupus erythematosus (SLE). Risk alleles increase TNFSF4 expression in cell lines, but the mechanism linking this effect to disease is unclear, and the OX40L-expressing cell types mediating the risk are not clearly established. Blockade of OX40L has been demonstrated to reduce disease severity in several models of autoimmunity, but not in SLE. We sought to investigate its potential therapeutic role in lupus.

Methods

We used a conditional knockout mouse system to investigate the function of OX40L on B and T lymphocytes in systemic autoimmunity.

Results

Physiologically, OX40L on both B and T cells contributed to the humoral immune response, but B cell OX40L supported the secondary humoral response and antibody affinity maturation. Our data also indicated that loss of B cell OX40L impeded the generation of splenic T follicular helper cells. We further show that in two models of SLE—a spontaneous congenic model and the H2-IA^bm12 graft-versus-host-induced model—loss of B cell OX40L ameliorates the autoimmune phenotype. This improvement was, in each case, accompanied by a decline in T follicular helper cell numbers. Importantly, the germline knockout did not exhibit a markedly different phenotype from the B cell knockout in these models.

Conclusions

These findings contribute to a model in which genetically determined increased OX40L expression promotes human SLE by several mechanisms, contingent on its cellular expression. The improvement in pathology in two models of systemic autoimmunity indicates that OX40L is an excellent therapeutic target in SLE.

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.

Mapping eQTLs with RNA-seq reveals novel susceptibility genes, non-coding RNAs and alternative-splicing events in systemic lupus erythematosus

Studies attempting to functionally interpret complex-disease susceptibility loci by GWAS and eQTL integration have predominantly employed microarrays to quantify gene-expression. RNA-Seq has the potential to discover a more comprehensive set of eQTLs and illuminate the underlying molecular consequence. We examine the functional outcome of 39 variants associated with Systemic Lupus Erythematosus (SLE) through the integration of GWAS and eQTL data from the TwinsUK microarray and RNA-Seq cohort in lymphoblastoid cell lines. We use conditional analysis and a Bayesian colocalisation method to provide evidence of a shared causal-variant, then compare the ability of each quantification type to detect disease relevant eQTLs and eGenes. We discovered the greatest frequency of candidate-causal eQTLs using exon-level RNA-Seq, and identified novel SLE susceptibility genes (e.g. NADSYN1 and TCF7) that were concealed using microarrays, including four non-coding RNAs. Many of these eQTLs were found to influence the expression of several genes, supporting the notion that risk haplotypes may harbour multiple functional effects. Novel SLE associated splicing events were identified in the T-reg restricted transcription factor, IKZF2, and other candidate genes (e.g. WDFY4) through asQTL mapping using the Geuvadis cohort. We have significantly increased our understanding of the genetic control of gene-expression in SLE by maximising the leverage of RNA-Seq and performing integrative GWAS-eQTL analysis against gene, exon, and splice-junction quantifications. We conclude that to better understand the true functional consequence of regulatory variants, quantification by RNA-Seq should be performed at the exon-level as a minimum, and run in parallel with gene and splice-junction level quantification.

Genome-wide association meta-analysis in Chinese and European individuals identifies ten new loci associated with systemic lupus erythematosus

Systemic lupus erythematosus (SLE; OMIM 152700) is a genetically complex autoimmune disease. Genome-wide association studies (GWASs) have identified more than 50 loci as robustly associated with the disease in single ancestries, but genome-wide transancestral studies have not been conducted. We combined three GWAS data sets from Chinese (1,659 cases and 3,398 controls) and European (4,036 cases and 6,959 controls) populations. A meta-analysis of these studies showed that over half of the published SLE genetic associations are present in both populations. A replication study in Chinese (3,043 cases and 5,074 controls) and European (2,643 cases and 9,032 controls) subjects found ten previously unreported SLE loci. Our study provides further evidence that the majority of genetic risk polymorphisms for SLE are contained within the same regions across both populations. Furthermore, a comparison of risk allele frequencies and genetic risk scores suggested that the increased prevalence of SLE in non-Europeans (including Asians) has a genetic basis.

Defective removal of ribonucleotides from DNA promotes systemic autoimmunity

Genome integrity is continuously challenged by the DNA damage that arises during normal cell metabolism. Biallelic mutations in the genes encoding the genome surveillance enzyme ribonuclease H2 (RNase H2) cause Aicardi-Goutières syndrome (AGS), a pediatric disorder that shares features with the autoimmune disease systemic lupus erythematosus (SLE). Here we determined that heterozygous parents of AGS patients exhibit an intermediate autoimmune phenotype and demonstrated a genetic association between rare RNASEH2 sequence variants and SLE. Evaluation of patient cells revealed that SLE- and AGS-associated mutations impair RNase H2 function and result in accumulation of ribonucleotides in genomic DNA. The ensuing chronic low level of DNA damage triggered a DNA damage response characterized by constitutive p53 phosphorylation and senescence. Patient fibroblasts exhibited constitutive upregulation of IFN-stimulated genes and an enhanced type I IFN response to the immunostimulatory nucleic acid polyinosinic:polycytidylic acid and UV light irradiation, linking RNase H2 deficiency to potentiation of innate immune signaling. Moreover, UV-induced cyclobutane pyrimidine dimer formation was markedly enhanced in ribonucleotide-containing DNA, providing a mechanism for photosensitivity in RNase H2-associated SLE. Collectively, our findings implicate RNase H2 in the pathogenesis of SLE and suggest a role of DNA damage-associated pathways in the initiation of autoimmunity.

Variants at multiple loci implicated in both innate and adaptive immune responses are associated with Sjögren's syndrome

Sjögren’s syndrome is a common autoimmune disease (~0.7% of European Americans) typically presenting as keratoconjunctivitis sicca and xerostomia. In addition to strong association within the HLA region at 6p21 (P_meta=7.65×10⁻¹¹⁴), we establish associations with IRF5-TNPO3 (P_meta=2.73×10⁻¹⁹), STAT4 (P_meta=6.80×10⁻¹⁵), IL12A (P_meta =1.17×10⁻¹⁰), FAM167A-BLK (P_meta=4.97×10⁻¹⁰), DDX6-CXCR5 (P_meta=1.10×10⁻⁸), and TNIP1 (P_meta=3.30×10⁻⁸). Suggestive associations with P_meta<5×10⁻⁵ were observed with 29 regions including TNFAIP3, PTTG1, PRDM1, DGKQ, FCGR2A, IRAK1BP1, ITSN2, and PHIP amongst others. These results highlight the importance of genes involved in both innate and adaptive immunity in Sjögren’s syndrome.

CSK regulatory polymorphism is associated with systemic lupus erythematosus and influences B-cell signaling and activation

C-src tyrosine kinase, Csk, physically interacts with the intracellular phosphatase Lyp (PTPN22) and can modify the activation state of downstream Src kinases, such as Lyn, in lymphocytes. We identified an association of Csk with systemic lupus erythematosus (SLE) and refined its location to an intronic polymorphism rs34933034 (OR 1.32, p = 1.04 × 10⁻⁹). The risk allele is associated with increased CSK expression and augments inhibitory phosphorylation of Lyn. In carriers of the risk allele, B cell receptor (BCR)-mediated activation of mature B cells, as well as plasma IgM, are increased. Moreover, the fraction of transitional B cells is doubled in the cord blood of carriers of the risk allele compared to non-risk haplotypes due to an expansion of the late transitional cells, a stage targeted by selection mechanisms. This suggests that the Lyp-Csk complex increases susceptibility to lupus at multiple maturation and activation points of B cells.

The genetics of lupus: a functional perspective

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic component and is characterized by chronic inflammation and the production of anti-nuclear auto-antibodies. In the era of genome-wide association studies (GWASs), elucidating the genetic factors present in SLE has been a very successful endeavor; 28 confirmed disease susceptibility loci have been mapped. In this review, we summarize the current understanding of the genetics of lupus and focus on the strongest associated risk loci found to date (P <1.0 × 10−8). Although these loci account for less than 10% of the genetic heritability and therefore do not account for the bulk of the disease heritability, they do implicate important pathways, which contribute to SLE pathogenesis. Consequently, the main focus of the review is to outline the genetic variants in the known associated loci and then to explore the potential functional consequences of the associated variants. We also highlight the genetic overlap of these loci with other autoimmune diseases, which indicates common pathogenic mechanisms. The importance of developing functional assays will be discussed and each of them will be instrumental in furthering our understanding of these associated variants and loci. Finally, we indicate that performing a larger SLE GWAS and applying a more targeted set of methods, such as the ImmunoChip and next generation sequencing methodology, are important for identifying additional loci and enhancing our understanding of the pathogenesis of SLE.

Dense mapping of IL18 shows no association in SLE

Systemic lupus erythematosus (SLE) is an autoimmune disease which behaves as a complex genetic trait. At least 20 SLE risk susceptibility loci have been mapped using both candidate gene and genome-wide association strategies. The gene encoding the pro-inflammatory cytokine, IL18, has been reported as a candidate gene showing an association with SLE. This pleiotropic cytokine is expressed in a range of immune cells and has been shown to induce interferon-γ and tumour necrosis factor-α. Serum interleukin-18 has been reported to be elevated in patients with SLE. Here we aimed to densely map single nucleotide polymorphisms (SNPs) across IL18 to investigate the association across this locus. We genotyped 36 across IL18 by Illumina bead express in 372 UK SLE trios. We also genotyped these SNPs in a further 508 non-trio UK cases and were able to accurately impute a dense marker set across IL18 in WTCCC2 controls with a total of 258 SNPs. To improve the study's power, we also imputed a total of 158 SNPs across the IL18 locus using data from an SLE genome-wide association study and performed association testing. In total, we analysed 1818 cases and 10 770 controls in this study. Our large well-powered study (98% to detect odds ratio = 1.5, with respect to rs360719) showed that no individual SNP or haplotype was associated with SLE in any of the cohorts studied. We conclude that we were unable to replicate the SLE association with rs360719 located upstream of IL18. No evidence for association with any other common variant at IL18 with SLE was found.

Genome-wide association studies in systemic lupus erythematosus: a perspective

Genome-wide association studies (GWAS) have been shown to be a powerful way of identifying novel susceptibility genes in systemic lupus erythematosus (SLE), as demonstrated by a series of publications in the past year. Lupus has been a late-comer to the GWAS community, being preceded by success stories for the GWAS approach in other autoimmune diseases, including type I diabetes, ankylosing spondylitis, rheumatoid arthritis, Crohn's disease and ulcerative colitis. The paper by Suarez-Gestal and colleagues seeks to exploit the wealth of data available from a total of four GWAS in SLE, three in European-American populations and one in a Swedish population. The authors describe replication of ten lupus susceptibility alleles in a Spanish SLE case-control study.

Association of the co-stimulator OX40L with systemic lupus erythematosus

The archetypal systemic autoimmune disease systemic lupus erythematosus (SLE) has incompletely understood pathogenesis, although evidence suggests a strong genetic component. Unlike organ-specific autoimmune diseases such as type 1 diabetes, the genetics of lupus are not as dominated by the effect of a single locus. Undoubtedly, the major histocompatibility complex is the greatest and most consistent genetic risk factor in SLE susceptibility; however, recent candidate gene and whole genome association (WGA) studies have identified several other genes that are likely to advance our understanding of this complex disease. One of these, the TNF superfamily member OX40L, interacts with its unique receptor OX40, to maintain T cell memory by providing a late-stage co-stimulatory signal to sustain the survival of activated T cells. The precise immunological consequences are yet to be determined; however, signalling through OX40-OX40L is bidirectional and the reverse signalling pathway via OX40L may quantitatively enhance B cell proliferation to augment the B cell hyperactivity found in SLE. Like OX40L, several genes recently identified in WGA studies are components of B cell pathways. Collectively, these genes will help us to unravel the mechanisms by which aberrant B cell signalling results in lupus pathogenesis.

FCGR3B copy number variation is associated with susceptibility to systemic, but not organ-specific, autoimmunity

Naturally occurring variation in gene copy number is increasingly recognized as a heritable source of susceptibility to genetically complex diseases. Here we report strong association between FCGR3B copy number and risk of systemic lupus erythematosus (P = 2.7 x 10(-8)), microscopic polyangiitis (P = 2.9 x 10(-4)) and Wegener's granulomatosis in two independent cohorts from the UK (P = 3 x 10(-3)) and France (P = 1.1 x 10(-4)). We did not observe this association in the organ-specific Graves' disease or Addison's disease. Our findings suggest that low FCGR3B copy number, and in particular complete FCGR3B deficiency, has a key role in the development of systemic autoimmunity.

Three functional variants of IFN regulatory factor 5 (IRF5) define risk and protective haplotypes for human lupus

Systematic genome-wide studies to map genomic regions associated with human diseases are becoming more practical. Increasingly, efforts will be focused on the identification of the specific functional variants responsible for the disease. The challenges of identifying causal variants include the need for complete ascertainment of genetic variants and the need to consider the possibility of multiple causal alleles. We recently reported that risk of systemic lupus erythematosus (SLE) is strongly associated with a common SNP in IFN regulatory factor 5 (IRF5), and that this variant altered spicing in a way that might provide a functional explanation for the reproducible association to SLE risk. Here, by resequencing and genotyping in patients with SLE, we find evidence for three functional alleles of IRF5: the previously described exon 1B splice site variant, a 30-bp in-frame insertion/deletion variant of exon 6 that alters a proline-, glutamic acid-, serine- and threonine-rich domain region, and a variant in a conserved polyA+ signal sequence that alters the length of the 3' UTR and stability of IRF5 mRNAs. Haplotypes of these three variants define at least three distinct levels of risk to SLE. Understanding how combinations of variants influence IRF5 function may offer etiological and therapeutic insights in SLE; more generally, IRF5 and SLE illustrates how multiple common variants of the same gene can together influence risk of common disease.

Association of IRF5 in UK SLE families identifies a variant involved in polyadenylation

Results from two studies have implicated the interferon regulatory gene IRF5 as a susceptibility gene in systemic lupus erythematosus (SLE). In this study, we conducted a family-based association analysis in 380 UK SLE nuclear families. Using a higher density of markers than has hitherto been screened, we show that there is association with two SNPs in the first intron, rs2004640 (P = 3.4 x 10(-4)) and rs3807306 (P = 4.9 x 10(-4)), and the association extends into the 3'-untranslated region (UTR). There is a single haplotype block encompassing IRF5 and we show for the first time that the gene comprises two over-transmitted haplotypes and a single under-transmitted haplotype. The strongest association is with a TCTAACT haplotype (T:U = 1.92, P = 5.8 x 10(-5)), which carries all the over-transmitted alleles from this study. Haplotypes carrying the T alleles of rs2004640 and rs2280714 and the A allele of rs10954213 are over-transmitted in SLE families. The TAT haplotype shows a dose-dependent relationship with mRNA expression. A differential expression pattern was seen between two expression probes located each side of rs10954213 in the 3'-UTR. rs10954213 shows the strongest association with RNA expression levels (P = 1 x 10(-14)). The A allele of rs10954213 creates a functional polyadenylation site and the A genotype correlates with increased expression of a transcript variant containing a shorter 3'-UTR. Expression levels of transcript variants with the shorter or longer 3'-UTRs are inversely correlated. Our data support a new mechanism by which an IRF5 polymorphism controls the expression of alternate transcript variants which may have different effects on interferon signalling.

Association analysis of the R620W polymorphism of protein tyrosine phosphatase PTPN22 in systemic lupus erythematosus families: increased T allele frequency in systemic lupus erythematosus patients with autoimmune thyroid disease

OBJECTIVE

Recent case-control studies show associations of the minor T allele (of the C1858T single-nucleotide polymorphism corresponding to the R620W amino acid substitution) of PTPN22 with multiple autoimmune diseases, including systemic lupus erythematosus (SLE). We performed family-based association studies of this polymorphism in 4 independent cohorts containing SLE patients and their parents and/or other family members.

METHODS

A total of 2,689 individuals from 902 independent Caucasian families with SLE were genotyped using polymerase chain reaction pyrosequencing (cohorts 1 and 2) and the Sequenom MassArray system (cohorts 3 and 4). The transmission disequilibrium test (TDT) and the pedigree disequilibrium test (PDT) were conducted to assess the evidence of association.

RESULTS

The 1858 C > T allele frequencies of the parents showed no deviation from Hardy-Weinberg equilibrium within each cohort. No evidence of preferential transmission of the T allele from heterozygous parents to their affected offspring was observed in each of the 4 cohorts or in the combined sample. Consistent with the TDT result, the PDT analysis revealed no significant association between the T allele and SLE. In 54 of the 661 SLE patients (cohorts 1 and 3) with documented autoimmune thyroid disease, the T allele frequency was higher than in individuals with SLE alone (16.7% versus 8.5%; P = 0.008, odds ratio 2.16 [95% confidence interval 1.25-3.72]).

CONCLUSION

The R620W polymorphism of the PTPN22 gene is not a major risk allele for SLE susceptibility in our sample of Caucasian individuals from northern America, the UK, or Finland, but it appears to be a risk factor for the concurrent autoimmune diseases of autoimmune thyroid disease and SLE.

Polymorphism at the C-reactive protein locus influences gene expression and predisposes to systemic lupus erythematosus

Relative deficiency of pentraxin proteins is implicated in the pathogenesis of systemic lupus erythematosus. The C-reactive protein (CRP) response is defective in patients with acute flares of disease, and mice with targeted deletions of the serum amyloid P component gene (Sap) develop a lupus-like illness. In humans, the genes for CRP (CRP) and SAP (APCS) map to 1q23.2 within an interval linked with SLE. We have investigated the candidate genes CRP and APCS in two cohorts totalling 586 UK simplex SLE families. The inheritance of an intronic dinucleotide repeat and seven single nucleotide polymorphisms in the CRP and APCS genes was examined by application of family-based tests of association and linkage. Basal levels of CRP were influenced independently by two polymorphisms at the CRP locus, CRP 2 and CRP 4. Furthermore, the latter polymorphism was linked/associated with SLE and antinuclear autoantibody production. Thus, the polymorphism associated with reduced basal CRP was also associated with the development of SLE. These data support the hypothesis that defective disposal of potentially immunogenic material is a contributory factor in lupus pathogenesis. The identification of polymorphisms that determine basal CRP levels has implications in ischaemic heart disease, where CRP level is an important predictor of risk.