Variation in the upstream region of P-Selectin (SELP) is a risk factor for SLE

Systemic lupus erythematosus (SLE) is a complex autoimmune disease. Genome-wide linkage studies implicated a region containing the adhesion molecule P-Selectin. This family-based study revealed two regions of association within P-Selectin. The strongest signal, from a 21.4-kb risk haplotype, stretched from the promoter into the first two consensus repeat (CR) regions (P=8 × 10−4), with a second association from a 14.6-kb protective haplotype covering CR 2–9 (P=0.0198). The risk haplotype is tagged by the rare C allele of rs3753306, which disrupts the binding site of the trans-activating transcription factor HNF-1. One other variant (rs3917687) on the risk haplotype was significant after permutation (P10000<1 × 10−5), replicated in independent pseudo case-control analysis and was significant by meta-analysis (P=4.37 × 10−6). A third associated variant on the risk haplotype (rs3917657) replicated in 306 US SLE families and was significant in a joint UK-SLE data set after permutation. The protective haplotype is tagged by rs6133 (a non-synonymous variant in CR8 (P=9.00 × 10−4), which also shows association in the pseudo case-control analysis (P=1.09 × 10−3) and may contribute to another signal in P-Selectin. We propose that polymorphism in the upstream region may reduce expression of P-Selectin, the mechanism by which this promotes autoimmunity is unknown, although it may reduce the production of regulatory T cells

Association of NCF2, IKZF1, IRF8, IFIH1, and TYK2 with Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is a complex trait characterised by the production of a range of auto-antibodies and a diverse set of clinical phenotypes. Currently, ∼8% of the genetic contribution to SLE in Europeans is known, following publication of several moderate-sized genome-wide (GW) association studies, which identified loci with a strong effect (OR>1.3). In order to identify additional genes contributing to SLE susceptibility, we conducted a replication study in a UK dataset (870 cases, 5,551 controls) of 23 variants that showed moderate-risk for lupus in previous studies. Association analysis in the UK dataset and subsequent meta-analysis with the published data identified five SLE susceptibility genes reaching genome-wide levels of significance (Pcomb<5×10−8): NCF2 (Pcomb = 2.87×10−11), IKZF1 (Pcomb = 2.33×10−9), IRF8 (Pcomb = 1.24×10−8), IFIH1 (Pcomb = 1.63×10−8), and TYK2 (Pcomb = 3.88×10−8). Each of the five new loci identified here can be mapped into interferon signalling pathways, which are known to play a key role in the pathogenesis of SLE. These results increase the number of established susceptibility genes for lupus to ∼30 and validate the importance of using large datasets to confirm associations of loci which moderately increase the risk for disease

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

Haematopoietic stem cell-derived immune cells have reduced X chromosome inactivation skewing in systemic lupus erythematosus

Objectives: Systemic lupus erythematosus (SLE) shows a marked female bias in prevalence. X chromosome inactivation (XCI) is the mechanism which randomly silences one X chromosome to equalise gene expression between 46, XX females and 46, XY males. Though XCI is expected to result in a random pattern of mosaicism across tissues, some females display a significantly skewed ratio in immune cells, termed XCI-skew. We tested whether XCI was abnormal in females with SLE and hence contributes to sexual dimorphism. Methods: We assayed XCI in whole blood DNA in 181 female SLE cases, 796 female healthy controls and 10 twin pairs discordant for SLE. Using regression modelling and intra-twin comparisons, we assessed the effect of SLE on XCI and combined clinical, cellular and genetic data via a polygenic score to explore underlying mechanisms. Results: Accommodating the powerful confounder of age, XCI-skew was reduced in females with SLE compared with controls (p=1.3×10 -5), with the greatest effect seen in those with more severe disease. Applying an XCI threshold of &gt;80%, we observed XCI-skew in 6.6% of SLE cases compared with 22% of controls. This difference was not explained by differential white cell counts, medication or genetic susceptibility to SLE. Instead, XCI-skew correlated with a biomarker for type I interferon-regulated gene expression. Conclusions: These results refute current views on XCI-skew in autoimmunity and suggest, in lupus, XCI patterns of immune cells reflect the impact of disease state, specifically interferon signalling, on the haematopoietic stem cells from which they derive.</p

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 (Pcohort1 = 2.44e-08; Pcohort2 = 0.00205) and a significant negative correlation with age of SLE onset (Pcohort1 = 1.76e-12; Pcohort2 = 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.China Scholarship Council (CSC) & National Science Foundation of China grant 8180163

Loci associated with N-glycosylation of human IgG are not associated with rheumatoid arthritis: a Mendelian randomisation study

Objectives: A recent study identified 16 genetic variants associated with N-glycosylation of human IgG. Several of the genomic regions where these single nucleotide polymorphisms (SNPs) reside have also been associated with autoimmune disease (AID) susceptibility, suggesting there may be pleiotropy (genetic sharing) between loci controlling both N-glycosylation and AIDs. We investigated this by testing variants associated with levels of IgG N-glycosylation for association with rheumatoid arthritis (RA) susceptibility using a Mendelian randomisation study, and testing a subset of these variants in a less well-powered study of treatment response and severity. Methods: SNPs showing association with IgG N-glycosylation were analysed for association with RA susceptibility in 14 361 RA cases and 43 923 controls. Five SNPs were tested for association with response to anti-tumour necrosis factor (TNF) therapy in 1081 RA patient samples and for association with radiological disease severity in 342 patients. Results: Only one SNP (rs9296009) associated with N-glycosylation showed an association (p=6.92×10–266) with RA susceptibility, although this was due to linkage disequilibrium with causal human leukocyte antigen (HLA) variants. Four regions of the genome harboured SNPs associated with both traits (shared loci); although statistical analysis indicated that the associations observed for the two traits are independent. No SNPs showed association with response to anti-TNF therapy. One SNP rs12342831 was modestly associated with Larsen score (p=0.05). Conclusions: In a large, well-powered cohort of RA patients, we show SNPs driving levels of N-glycosylation have no association with RA susceptibility, indicating colocalisation of associated SNPs are not necessarily indicative of a shared genetic background or a role for glycosylation in disease susceptibility

Identification of novel genetic susceptibility loci in African American lupus patients in a candidate gene association study

Objective Candidate gene and genome‐wide association studies have identified several disease susceptibility loci in lupus patients. These studies have largely been performed in lupus patients who are Asian or of European ancestry. This study was undertaken to examine whether some of these same susceptibility loci increase lupus risk in African American individuals. Methods Single‐nucleotide polymorphisms tagging 15 independent lupus susceptibility loci were genotyped in a set of 1,724 lupus patients and 2,024 healthy controls of African American descent. The loci examined included PTPN22 , FCGR2A , TNFSF4 , STAT4 , CTLA4 , PDCD1 , PXK , BANK1 , MSH5 (HLA region), CFB (HLA region), C8orf13‐BLK region, MBL2 , KIAA1542 , ITGAM , and MECP2 / IRAK1 . Results We found the first evidence of genetic association between lupus in African American patients and 5 susceptibility loci ( C8orf13‐BLK , BANK1 , TNFSF4 , KIAA1542 , and CTLA4 ; P = 8.0 × 10 −6 , P = 1.9 × 10 −5 , P = 5.7 × 10 −5 , P = 0.00099, and P = 0.0045, respectively). Further, we confirmed the genetic association between lupus and 5 additional lupus susceptibility loci ( ITGAM , MSH5 , CFB , STAT4 , and FCGR2A ; P = 7.5 × 10 −11 , P = 5.2 × 10 −8 , P = 8.7 × 10 −7 , P = 0.0058, and P = 0.0070, respectively), and provided evidence, for the first time, of genome‐wide significance for the association between lupus in African American patients and ITGAM and MSH5 (HLA region). Conclusion These findings provide evidence of novel genetic susceptibility loci for lupus in African Americans and demonstrate that the majority of lupus susceptibility loci examined confer lupus risk across multiple ethnicities.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/88030/1/30563_ftp.pd

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