Genome wide association studies have identified variants in PXK that confer risk for humoral autoimmune diseases, including systemic lupus erythematosus (SLE or lupus), rheumatoid arthritis and more recently systemic sclerosis. While PXK is involved in trafficking of epidermal growth factor Receptor (EGFR) in COS-7 cells, mechanisms linking PXK to lupus pathophysiology have remained undefined. In an effort to uncover the mechanism at this locus that increases lupus-risk, we undertook a fine-mapping analysis in a large multi-ancestral study of lupus patients and controls. We define a large (257kb) common haplotype marking a single causal variant that confers lupus risk detected only in European ancestral populations and spans the promoter through the 3′ UTR of PXK. The strongest association was found at rs6445972 with P < 4.62 × 10−10, OR 0.81 (0.75–0.86). Using stepwise logistic regression analysis, we demonstrate that one signal drives the genetic association in the region. Bayesian analysis confirms our results, identifying a 95% credible set consisting of 172 variants spanning 202 kb. Functionally, we found that PXK operates on the B-cell antigen receptor (BCR); we confirmed that PXK influenced the rate of BCR internalization. Furthermore, we demonstrate that individuals carrying the risk haplotype exhibited a decreased rate of BCR internalization, a process known to impact B cell survival and cell fate. Taken together, these data define a new candidate mechanism for the genetic association of variants around PXK with lupus risk and highlight the regulation of intracellular trafficking as a genetically regulated pathway mediating human autoimmunity.
lupus; PXK; fine-mapping; B cells; BCR
Multiple MHC loci encoding human leukocyte antigens (HLA) have allelic variants unequivocally associated with differential immune control of HIV-1 infection. Fine mapping based on single nucleotide polymorphisms (SNPs) in the extended MHC (xMHC) region is expected to reveal causal or novel factors and to justify a search for functional mechanisms. We have tested the utility of a custom fine-mapping platform (the ImmunoChip) for 172 HIV-1 seroconverters (SCs) and 449 seroprevalent individuals (SPs) from Lusaka, Zambia, with a focus on more than 6,400 informative xMHC SNPs. When conditioned on HLA and non-genetic factors previously associated with HIV-1 viral load (VL) in the study cohort, penalized approaches (HyperLasso models) identified an intergenic SNP (rs3094626 between RPP21 and HLA-E) and an intronic SNP (rs3134931 in NOTCH4) as novel correlates of early set-point VL in SCs. The minor allele of rs2857114 (downstream from HLA-DOB) was an unfavorable factor in SPs. Joint models based on demographic features, HLA alleles and the newly identified SNP variants could explain 29% and 15% of VL variance in SCs and SPs, respectively. These findings and bioinformatics strongly suggest that both classic and non-classic MHC genes deserve further investigation, especially in Africans with relatively short haplotype blocks.
HIV-1; HLA; human MHC; SNP; viral load
Methotrexate (MTX) has emerged as first-line therapy for early moderate to severe rheumatoid arthritis (RA), but individual variation in treatment response remains unexplained. We tested the associations between 863 known pharmacogenetic variants and MTX response in 471 TEAR Trial participants with early RA. Efficacy and toxicity were modeled using multiple regression, adjusted for demographic and clinical covariates. Penalized regression models were used to test joint associations of markers and/or covariates with the outcomes. The strongest genetic associations with efficacy were in CHST11 (five markers with P <0.003), encoding carbohydrate (chondroitin 4) sulfotransferase 11. Top markers associated with MTX toxicity were in the cytochrome p450 genes CYP20A1 and CYP39A1, solute carrier genes SLC22A2 and SLC7A7, and the mitochondrial aldehyde dehydrogenase gene ALDH2. The selected markers explained a consistently higher proportion of variation in toxicity than efficacy. These findings could inform future development of personalized therapeutic approaches.
Methotrexate; rheumatoid arthritis; pharmacogenetics
To characterize the clinical features of familial lupus, and determine its influence on damage accrual and survival using data from LUMINA, a longitudinal multiethnic US cohort.
Familial lupus was defined as patients with a first degree relative with SLE. Relative risks were estimated by logistic regression; odds ratios (OR) and their 95% confidence intervals (CI) were the measure of association for familial lupus. Hazard Ratios (HR) were calculated using Cox proportional hazard adjusted for potential confounders for damage and survival.
Thirty-two of 644 patients had familial and 612 had sporadic lupus; both groups were of comparable age (~ 36 years). Familial lupus patients were in decreasing order of frequency siblings, parents and children. In multivariable analyses, mucosal ulcers (OR=1.92, 95% CI 0.65–5.70), mitral valve prolapse (OR=1.74, 95% CI 0.50–6.10), cerebrovascular disease (OR=4.18, 95% CI 0.98–17.76) and oral contraceptive use (ever/never; OR=2.51, 95% CI 0.88–7.19) were more likely in familial lupus but a history of low platelet count (<150,000/mm3; OR=0.31, 95% CI 0.08–1.17) and pulmonary disease activity (OR=0.39, 95% CI 0.14–1.20) were less likely. However, none of these associations reached statistical significance. Familial lupus was not significantly associated with a shorter time to either damage accrual or death (HR=0.77, 95% CI 0.37–1.59, p = 0.4746 and HR=0.20, 95% CI 0.03–1.47, p = 0.2020, respectively).
Although some clinical differences were observed in patients with familial and sporadic lupus, familial lupus was not associated with a significantly greater disease burden (damage, survival) than sporadic lupus.
familial lupus; lupus; sporadic lupus; LUMINA; multiethnic cohort
Lupus nephritis (LN) is a severe manifestation of systemic lupus erythematosus (SLE) that exhibits familial aggregation and may progress to end-stage renal disease (ESRD). LN is more prevalent among African Americans than among European Americans. This study was undertaken to investigate the hypothesis that the apolipoprotein L1 gene (APOL1) nephropathy risk alleles G1/G2, common in African Americans and rare in European Americans, contribute to the ethnic disparity in risk.
APOL1 G1 and G2 nephropathy alleles were genotyped in 855 African American SLE patients with LN-ESRD (cases) and 534 African American SLE patients without nephropathy (controls) and tested for association under a recessive genetic model, by logistic regression.
Ninety percent of the SLE patients were female. The mean ± SD age at SLE diagnosis was significantly lower in LN-ESRD cases than in SLE non-nephropathy controls (27.3 ± 10.9 years versus 39.5 ± 12.2 years). The mean ± SD time from SLE diagnosis to development of LN-ESRD in cases was 7.3 ± 7.2 years. The G1/G2 risk alleles were strongly associated with SLE-ESRD, with 25% of cases and 12% of controls having 2 nephropathy alleles (odds ratio [OR] 2.57, recessive model P = 1.49 × 10−9), and after adjustment for age, sex, and ancestry admixture (OR 2.72, P = 6.23 × 10−6). The age-, sex-, and admixture-adjusted population attributable risk for ESRD among patients with G1/G2 polymorphisms was 0.26, compared to 0.003 among European American patients. The mean time from SLE diagnosis to ESRD development was ~2 years earlier among individuals with APOL1 risk genotypes (P = 0.01).
APOL1 G1/G2 alleles strongly impact the risk of LN-ESRD in African Americans, as well as the time to progression to ESRD. The high frequency of these alleles in African Americans with near absence in European Americans explains an important proportion of the increased risk of LN-ESRD in African Americans.
The Disease Activity Score based on 28 joints (DAS28) has been increasingly used in clinical practice and research studies of rheumatoid arthritis (RA). Studies have reported discordance between DAS28 based on erythrocyte sedimentation rate (ESR) versus C-reactive protein (CRP) in RA patients. However such comparison is lacking in African-Americans with RA.
This analysis included participants from the Consortium for the Longitudinal Evaluation of African Americans with Early Rheumatoid Arthritis (CLEAR) Registry which enrolls self-declared African-Americans with RA. Using tender and swollen joint counts separate ESR-based and CRP-based DAS28 scores (DAS28-ESR3 and DAS28-CRP3) were calculated, as were DAS28-ESR4 and DAS28-CRP4, which included the patient’s assessment of disease activity. The scores were compared using paired t-test, simple agreement and kappa, correlation coefficient and Bland-Altman plots.
Of the 233 included participants, 85% were women, mean age at enrollment was 52.6 years, and median disease duration at enrollment was 21 months. Mean DAS28-ESR3 was significantly higher than DAS28-CRP3 (4.8 vs. 3.9; p<0.001). Similarly, mean DAS28-ESR4 was significantly higher than DAS28-CRP4 (4.7 vs. 3.9; p<0.001). ESR-based DAS28 remained higher than CRP-based DAS28 even when stratified by age, sex, and disease duration. Overall agreement was not high between DAS28-ESR3 and DAS28-CRP3 (50%) or between DAS28-ESR4 and DAS28-CRP4 (59%). DAS28-CRP3 underestimated disease activity in 47% of the participants relative to DAS28-ESR3 and DAS28-CRP4 in 40% of the participants relative to DAS28-ESR4.
There was significant discordance between the ESR-based and CRP-based DAS28 which could impact clinical treatment decisions in African-Americans with RA.
DAS28; Rheumatoid Arthritis; African-Americans
The increased risk of thrombosis in systemic lupus erythematosus (SLE) may be partially explained by interrelated genetic pathways for thrombosis and SLE. In a case-control analysis, we investigated whether 33 established and novel single nucleotide polymorphisms (SNP) in 20 genes involved in hemostasis pathways that have been associated with deep venous thrombosis in the general population were risk factors for SLE development among Asians.
Patients in the discovery cohort were enrolled in one of two North American SLE cohorts. Patients in the replication cohort were enrolled in one of four Asian or two North American cohorts. SLE cases met American College of Rheumatology classification criteria. We first genotyped 263 Asian SLE and 357 healthy Asian control individuals for 33 SNPs using Luminex multiplex technology in the discovery phase, and then used Taqman and Immunochip assays to examine 5 SNPs in up to an additional 1496 cases and 993 controls in the Replication phase. SLE patients were compared to healthy controls for association with minor alleles in allelic models. Principal components analysis was used to control for intra-Asian ancestry in an analysis of the replication cohort.
Two genetic variants in the gene VKORC1, rs9934438 and rs9923231, were highly significant in both the discovery and replication cohorts: OR(disc) = 2.45 (p=2×10−9), OR(rep) = 1.53 (p=5×10−6) and OR(disc) = 2.40 (p=6×10−9), OR(rep) = 1.53 (p=5×10−6), respectively. These associations were significant in the replication cohort after adjustment for intra-Asian ancestry: rs9934438 OR(adj) = 1.34 (p=0.0029) and rs9923231 OR(adj) = 1.34 (p=0.0032).
Genetic variants in VKORC1, involved in vitamin K reduction and associated with DVT, are associated with SLE development in Asians. These results suggest intersecting genetic pathways for the development of SLE and thrombosis.
systemic lupus erythematosus; single nucleotide polymorphisms; genetic risk factors
The vast initial diversity of the antibody repertoire is generated centrally by means of a complex series of V(D)J gene rearrangement events, variation in the site of gene segment joining, and TdT catalyzed N-region addition. Although the diversity is great, close inspection has revealed distinct and unique characteristics in the antibody repertoires expressed by different B cell developmental subsets. In order to illustrate our approach to repertoire analysis, we present an in-depth comparison of V(D)J gene usage, hydrophobicity, length, DH reading frame, and amino acid usage between heavy chain repertoires expressed by immature, transitional, mature, memory IgD+, memory IgD−, and plasmacytes isolated from the blood of a single individual. Our results support the view that in both human and mouse, the H chain repertoires expressed by individual, developmental B cell subsets appear to differ in sequence content. Sequencing of unsorted B cells from the blood is thus likely to yield an incomplete or compressed view of what is actually happening in the immune response of the individual. Our findings support the view that studies designed to correlate repertoire expression with diseases of immune function will likely require deep sequencing of B cells sorted by subset.
human antibody repertoire; CDR-H3; B cells subsets
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 more than 30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting the NF-κB signaling. In order to better understand the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway, we analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog, TNIP2, in case-control sets of diverse ethnic origins.
We fine-mapped TNIP1, TNIP2, and TAX1BP1 in a total of 8372 SLE cases and 7492 healthy controls from European-ancestry, African-American, Hispanic, East Asian, and African-American Gullah populations. Levels of TNIP1 mRNA and ABIN1 protein were analyzed using quantitative RT-PCR and Western blotting, respectively, in EBV-transformed human B cell lines.
We found significant associations between genetic variants within TNIP1 and SLE but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified two independent risk haplotypes within TNIP1 in individuals of European-ancestry that were also present in African-American and Hispanic populations. These risk haplotypes produced lower levels of TNIP1 mRNA and ABIN1 protein suggesting they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression.
Our results confirmed 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.
Amerindian-Europeans, Asians and African-Americans have an excess morbidity from SLE and higher prevalence of lupus nephritis than Caucasians. The aim of this study was to analyze the relationship between genetic ancestry and socio-demographic characteristics and clinical features in a large cohort of Amerindian-European SLE patients.
A total of 2116 SLE patients of Amerindian-European origin and 4001 SLE patients of European descent with clinical data were used in the study. Genotyping of 253 continental ancestry informative markers was performed on the Illumina platform. The STRUCTURE and ADMIXTURE software were used to determine genetic ancestry of each individual. Correlation between ancestry and socio-demographic and clinical data were analyzed using logistic regression.
The average Amerindian genetic ancestry of 2116 SLE patients was 40.7%. There was an increased risk of having renal involvement (P<0.0001, OR= 3.50 95%CI 2.63-4.63) and an early age of onset with the presence of Amerindian genetic ancestry (P<0.0001). Amerindian ancestry protected against photosensitivity (P<0.0001, OR= 0.58 95%CI 0.44-0.76), oral ulcers (P<0.0001, OR= 0.55 95%CI 0.42-0.72), and serositis (P<0.0001, OR= 0.56 95%CI 0.41-0.75) after adjustment by age, gender and age of onset. However, gender and age of onset had stronger effects on malar rash, discoid rash, arthritis and neurological involvement than genetic ancestry.
In general, genetic Amerindian ancestry correlates with lower socio-demographic status and increases the risk for developing renal involvement and SLE at an earlier age of onset.
Immunoregulatory cytokine interleukin-10 (IL-10) is elevated in sera from patients with systemic lupus erythematosus (SLE) correlating with disease activity. The established association of IL10 with SLE and other autoimmune diseases led us to fine map causal variant(s) and to explore underlying mechanisms. We assessed 19 tag SNPs, covering the IL10 gene cluster including IL19, IL20 and IL24, for association with SLE in 15,533 case and control subjects from four ancestries. The previously reported IL10 variant, rs3024505 located at 1 kb downstream of IL10, exhibited the strongest association signal and was confirmed for association with SLE in European American (EA) (P = 2.7×10−8, OR = 1.30), but not in non-EA ancestries. SNP imputation conducted in EA dataset identified three additional SLE-associated SNPs tagged by rs3024505 (rs3122605, rs3024493 and rs3024495 located at 9.2 kb upstream, intron 3 and 4 of IL10, respectively), and SLE-risk alleles of these SNPs were dose-dependently associated with elevated levels of IL10 mRNA in PBMCs and circulating IL-10 protein in SLE patients and controls. Using nuclear extracts of peripheral blood cells from SLE patients for electrophoretic mobility shift assays, we identified specific binding of transcription factor Elk-1 to oligodeoxynucleotides containing the risk (G) allele of rs3122605, suggesting rs3122605 as the most likely causal variant regulating IL10 expression. Elk-1 is known to be activated by phosphorylation and nuclear localization to induce transcription. Of interest, phosphorylated Elk-1 (p-Elk-1) detected only in nuclear extracts of SLE PBMCs appeared to increase with disease activity. Co-expression levels of p-Elk-1 and IL-10 were elevated in SLE T, B cells and monocytes, associated with increased disease activity in SLE B cells, and were best downregulated by ERK inhibitor. Taken together, our data suggest that preferential binding of activated Elk-1 to the IL10 rs3122605-G allele upregulates IL10 expression and confers increased risk for SLE in European Americans.
Systemic lupus erythematosus (SLE), a debilitating autoimmune disease characterized by the production of pathogenic autoantibodies, has a strong genetic basis. Variants of the IL10 gene, which encodes cytokine interleukin-10 (IL-10) with known function of promoting B cell hyperactivity and autoantibody production, are associated with SLE and other autoimmune diseases, and serum IL-10 levels are elevated in SLE patients correlating with increased disease activity. In this study, to discover SLE-predisposing causal variant(s), we assessed variants within the genomic region containing IL10 and its gene family member IL19, IL20 and IL24 for association with SLE in case and control subjects from diverse ancestries. We identified SLE-associated SNP rs3122605 located at 9.2 kb upstream of IL10 as the most likely causal variant in subjects of European ancestry. The SLE-risk allele of rs3122605 was dose-dependently associated with elevated IL10 expression at both mRNA and protein levels in peripheral blood samples from SLE patients and controls, which could be explained, at least in part, by its preferential binding to Elk-1, a transcription factor activated in B cells during active disease of SLE patients. Elk-1-mediated IL-10 overexpression could be downregulated by inhibiting activation of mitogen-activated protein kinases, suggesting a potential therapeutic target for SLE.
Methotrexate (MTX) is one of the most commonly prescribed and most effective drugs for the treatment of rheumatoid arthritis (RA). Given the partial response of many patients and the side effect profile of the drug, there is considerable interest in identification of biomarkers to guide MTX therapy in RA. Upon entering cells, MTX is polyglutamated. Measuring methotrexate polyglutamates (MTX PGs) levels in circulating red blood cells (RBC) has been proposed as an objective measure that can help to optimize MTX therapy in RA. There is conflicting data with regard to the clinical utility of measurement of MTX PGs measurements as a predictor of the efficacy or toxicity of low-dose MTX effects in RA. Should large, randomized clinical trials of this assay show consistent, reproducible, long-term correlations between MTX PG levels and efficacy and toxicity, this test could become a prominent tool for clinicians to optimize the use of MTX in RA.
Rheumatoid Arthritis; Methotrexate; Polyglutamates
Protein tyrosine phosphatase non-receptor type 22 (PTPN22) is a negative regulator of T-cell activation associated with several autoimmune diseases, including systemic lupus erythematosus (SLE). Missense rs2476601 is associated with SLE in individuals with European ancestry. Since the rs2476601 risk allele frequency differs dramatically across ethnicities, we assessed robustness of PTPN22 association with SLE and its clinical sub-phenotypes across four ethnically diverse populations. Ten SNPs were genotyped in 8220 SLE cases and 7369 controls from in European-Americans (EA), African-Americans (AA), Asians (AS), and Hispanics (HS). We performed imputation-based association followed by conditional analysis to identify independent associations. Significantly associated SNPs were tested for association with SLE clinical sub-phenotypes, including autoantibody profiles. Multiple testing was accounted for by using false discovery rate. We successfully imputed and tested allelic association for 107 SNPs within the PTPN22 region and detected evidence of ethnic-specific associations from EA and HS. In EA, the strongest association was at rs2476601 (P = 4.7×10−9, OR = 1.40 (95% CI = 1.25–1.56)). Independent association with rs1217414 was also observed in EA, and both SNPs are correlated with increased European ancestry. For HS imputed intronic SNP, rs3765598, predicted to be a cis-eQTL, was associated (P = 0.007, OR = 0.79 and 95% CI = 0.67–0.94). No significant associations were observed in AA or AS. Case-only analysis using lupus-related clinical criteria revealed differences between EA SLE patients positive for moderate to high titers of IgG anti-cardiolipin (aCL IgG >20) versus negative aCL IgG at rs2476601 (P = 0.012, OR = 1.65). Association was reinforced when these cases were compared to controls (P = 2.7×10−5, OR = 2.11). Our results validate that rs2476601 is the most significantly associated SNP in individuals with European ancestry. Additionally, rs1217414 and rs3765598 may be associated with SLE. Further studies are required to confirm the involvement of rs2476601 with aCL IgG.
Little is known about the genetic etiology of systemic lupus erythematosus (SLE) in individuals of African ancestry, despite its higher prevalence and greater disease severity. Overproduction of nitric oxide (NO) and reactive oxygen species are implicated in the pathogenesis and severity of SLE, making NO synthases and other reactive intermediate related genes biological candidates for disease susceptibility. This study analyzed variation in reactive intermediate genes for association with SLE in two populations with African ancestry.
A total of 244 SNPs from 53 regions were analyzed in non-Gullah African Americans (AA; 1432 cases and 1687 controls) and the genetically more homogeneous Gullah of the Sea Islands of South Carolina (133 cases and 112 controls) and. Single-marker, haplotype, and two-locus interaction tests were computed for these populations.
The glutathione reductase gene GSR (rs2253409, P=0.0014, OR [95% CI]=1.26 [1.09–1.44]) was the most significant single-SNP association in AA. In the Gullah, the NADH dehydrogenase NDUFS4 (rs381575, P=0.0065, OR [95%CI]=2.10 [1.23–3.59]) and nitric oxide synthase gene NOS1 (rs561712, P=0.0072, OR [95%CI]=0.62 [0.44–0.88]) were most strongly associated with SLE. When both populations were analyzed together, GSR remained the most significant effect (rs2253409, P=0.00072, OR [95%CI]=1.26 [1.10–1.44]). Haplotype and two-locus interaction analyses also uncovered different loci in each population.
These results suggest distinct patterns of association with SLE in African-derived populations; specific loci may be more strongly associated within select population groups.
systemic lupus erythematosus; African Americans; genetic association studies; oxygen compounds; single nucleotide polymorphism
We previously established an 80 kb haplotype upstream of TNFSF4 as a susceptibility locus in the autoimmune disease SLE. SLE-associated alleles at this locus are associated with inflammatory disorders, including atherosclerosis and ischaemic stroke. In Europeans, the TNFSF4 causal variants have remained elusive due to strong linkage disequilibrium exhibited by alleles spanning the region. Using a trans-ancestral approach to fine-map the locus, utilising 17,900 SLE and control subjects including Amerindian/Hispanics (1348 cases, 717 controls), African-Americans (AA) (1529, 2048) and better powered cohorts of Europeans and East Asians, we find strong association of risk alleles in all ethnicities; the AA association replicates in African-American Gullah (152,122). The best evidence of association comes from two adjacent markers: rs2205960-T (P = 1.71×10−34, OR = 1.43[1.26–1.60]) and rs1234317-T (P = 1.16×10−28, OR = 1.38[1.24–1.54]). Inference of fine-scale recombination rates for all populations tested finds the 80 kb risk and non-risk haplotypes in all except African-Americans. In this population the decay of recombination equates to an 11 kb risk haplotype, anchored in the 5′ region proximal to TNFSF4 and tagged by rs2205960-T after 1000 Genomes phase 1 (v3) imputation. Conditional regression analyses delineate the 5′ risk signal to rs2205960-T and the independent non-risk signal to rs1234314-C. Our case-only and SLE-control cohorts demonstrate robust association of rs2205960-T with autoantibody production. The rs2205960-T is predicted to form part of a decameric motif which binds NF-κBp65 with increased affinity compared to rs2205960-G. ChIP-seq data also indicate NF-κB interaction with the DNA sequence at this position in LCL cells. Our research suggests association of rs2205960-T with SLE across multiple groups and an independent non-risk signal at rs1234314-C. rs2205960-T is associated with autoantibody production and lymphopenia. Our data confirm a global signal at TNFSF4 and a role for the expressed product at multiple stages of lymphocyte dysregulation during SLE pathogenesis. We confirm the validity of trans-ancestral mapping in a complex trait.
Systemic lupus erythematosus (SLE/lupus) is a complex disease in which the body's immune cells cause inflammation in one or more systems to cause the associated morbidity. Hormones, the environment and genes are all causal contributors to SLE and over the past several years the genetic component of SLE has been firmly established. Several genes which are regulators of the immune system are associated with disease risk. We have established one of these, the tumour-necrosis family superfamily member 4 (TNFSF4) gene, as a lupus susceptibility gene in Northern Europeans. A major obstacle in pinpointing the marker(s) at TNFSF4 which best explain the risk of SLE has been the strong correlation (linkage disequilibrium, LD) between adjacent markers across the TNFSF4 region in this population. To address this, we have typed polymorphisms in several populations in addition to the European groups. The mixed ancestry of these populations gives a different LD pattern than that found in Europeans, presenting a method of pinpointing the section of the TNFSF4 region which results in SLE susceptibility. The Non-European populations have allowed identification of a polymorphism likely to regulate expression of TNFSF4 to increase susceptibility to SLE.
Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease with significant immune system aberrations resulting from complex heritable genetics as well as environmental factors. TRAF6 is a candidate gene for SLE, which has a major role in several signaling pathways that are important for immunity and organ development.
Fifteen single-nucleotide polymorphisms (SNPs), across TRAF6 were evaluated in 7,490 SLE and 6,780 control subjects from different ancestries. Population-based case-control association analyses and meta-analyses were performed. P values, false discovery rate q values, and odds ratios with 95% confidence intervals were calculated.
Evidence of associations in multiple SNPs was detected. The best overall p values were obtained for SNPs rs5030437 and rs4755453 (p=7.85×10−5 and p=4.73×10−5, respectively) without significant heterogeneity among populations (p=0.67 and p=0.50 in Q-statistic). In addition, rs540386 previously reported to be associated with RA was found to be in LD with these two SNPs (r2= 0.95) and demonstrated evidence of association with SLE in the same direction (meta-analysis p=9.15×10−4, OR=0.89, 95%CI=0.83–0.95). Thrombocytopenia improved the overall results in different populations (meta-analysis p=1.99×10−6, OR=0.57, 95%CI=0.45–0.72, for rs5030470). Finally evidence of family based association in 34 African-American pedigrees with the presence of thrombocytopenia were detected in one available SNP rs5030437 with Z score magnitude of 2.28 (p=0.02) under a dominant model.
Our data indicate the presence of association of TRAF6 with SLE in agreement with the previous report of association with RA. These data provide further support for the involvement of TRAF6 in the pathogenesis of autoimmunity.
TRAF6; polymorphism; systemic lupus erythematosus
To determine the clinical manifestations and disease damage associated with discoid rash in a large multiethnic systemic lupus erythematosus (SLE) cohort.
SLE patients (per ACR criteria), age ≥ 16 years, disease duration ≤ 10 years at enrollment, and defined ethnicity (African American, Hispanic or Caucasian), from a longitudinal cohort were studied. Socioeconomic-demographic features, clinical manifestations and disease damage [as per the Systemic Lupus International Collaborating Clinics Damage Index (SDI)] were determined. The association of DLE with clinical manifestations and disease damage was examined using multivariable logistic regression.
A total of 2,228 SLE patients were studied. The mean (standard deviation, SD) age at diagnosis was 34.3 (12.8) years and the mean (SD) disease duration was 7.9 (6.0) years; 91.8% were women. Discoid lupus was observed in 393 (17.6%) of patients with SLE. In the multivariable analysis, patients with discoid lupus were more likely to be smokers and of African-American ethnicity, and to have malar rash, photosensitivity, oral ulcers, leukopenia and vasculitis. DLE patients were less likely to be of Hispanic (from Texas) ethnicity, and to have arthritis, end-stage renal disease (ESRD), and antinuclear, anti-dsDNA and anti-phospholipid antibodies. Patients with DLE had more damage accrual, particularly chronic seizures, scarring alopecia, scarring of the skin, and skin ulcers.
In this cohort of SLE patients, discoid lupus was associated with several clinical features including serious manifestations such as vasculitis and chronic seizures.
discoid rash; systemic lupus erythematosus; disease damage
The Xq28 region containing IRAK1 and MECP2 has been identified as a risk locus for systemic lupus erythematosus (SLE) in previous genetic association studies. However, due to the strong linkage disequilibrium between IRAK1 and MECP2, it remains unclear which gene is affected by the underlying causal variant(s) conferring risk of SLE.
We fine-mapped ≥136 SNPs in a ~227kb region on Xq28, containing IRAK1, MECP2 and 7 adjacent genes (L1CAM, AVPR2, ARHGAP4, NAA10, RENBP, HCFC1 and TMEM187), for association with SLE in 15,783 case-control subjects derived from 4 different ancestral groups.
Multiple SNPs showed strong association with SLE in European Americans, Asians and Hispanics at P<5×10−8 with consistent association in subjects with African ancestry. Of these, 6 SNPs located in the TMEM187-IRAK1-MECP2 region captured the underlying causal variant(s) residing in a common risk haplotype shared by all 4 ancestral groups. Among them, rs1059702 best explained the Xq28 association signals in conditional testings and exhibited the strongest P value in trans-ancestral meta-analysis (Pmeta=1.3×10−27, OR=1.43), and thus was considered to be the most-likely causal variant. The risk allele of rs1059702 results in the amino acid substitution S196F in IRAK1 and had previously been shown to increase NF-κB activity in vitro. We also found that the homozygous risk genotype of rs1059702 was associated with lower mRNA levels of MECP2, but not IRAK1, in SLE patients (P=0.0012) and healthy controls (P=0.0064).
These data suggest contributions of both IRAK1 and MECP2 to SLE susceptibility.
Systemic Lupus Erythematosus; Gene Polymorphism; Xq28; IRAK1; MECP2
We previously reported that the G allele of rs3853839 at 3′untranslated region (UTR) of Toll-like receptor 7 (TLR7) was associated with elevated transcript expression and increased risk for systemic lupus erythematosus (SLE) in 9,274 Eastern Asians [P = 6.5×10−10, odds ratio (OR) (95%CI) = 1.27 (1.17–1.36)]. Here, we conducted trans-ancestral fine-mapping in 13,339 subjects including European Americans, African Americans, and Amerindian/Hispanics and confirmed rs3853839 as the only variant within the TLR7-TLR8 region exhibiting consistent and independent association with SLE (Pmeta = 7.5×10−11, OR = 1.24 [1.18–1.34]). The risk G allele was associated with significantly increased levels of TLR7 mRNA and protein in peripheral blood mononuclear cells (PBMCs) and elevated luciferase activity of reporter gene in transfected cells. TLR7 3′UTR sequence bearing the non-risk C allele of rs3853839 matches a predicted binding site of microRNA-3148 (miR-3148), suggesting that this microRNA may regulate TLR7 expression. Indeed, miR-3148 levels were inversely correlated with TLR7 transcript levels in PBMCs from SLE patients and controls (R2 = 0.255, P = 0.001). Overexpression of miR-3148 in HEK-293 cells led to significant dose-dependent decrease in luciferase activity for construct driven by TLR7 3′UTR segment bearing the C allele (P = 0.0003). Compared with the G-allele construct, the C-allele construct showed greater than two-fold reduction of luciferase activity in the presence of miR-3148. Reduced modulation by miR-3148 conferred slower degradation of the risk G-allele containing TLR7 transcripts, resulting in elevated levels of gene products. These data establish rs3853839 of TLR7 as a shared risk variant of SLE in 22,613 subjects of Asian, EA, AA, and Amerindian/Hispanic ancestries (Pmeta = 2.0×10−19, OR = 1.25 [1.20–1.32]), which confers allelic effect on transcript turnover via differential binding to the epigenetic factor miR-3148.
Systemic lupus erythematosus (SLE) is a debilitating autoimmune disease contributed to by excessive innate immune activation involving toll-like receptors (TLRs, particularly TLR7/8/9) and type I interferon (IFN) signaling pathways. TLR7 responds against RNA–containing nuclear antigens and activates IFN-α pathway, playing a pivotal role in the development of SLE. While a genomic duplication of Tlr7 promotes lupus-like disease in the Y-linked autoimmune accelerator (Yaa) murine model, the lack of common copy number variations at TLR7 in humans led us to identify a functional single nucleotide polymorphism (SNP), rs3853839 at 3′ UTR of the TLR7 gene, associated with SLE susceptibility in Eastern Asians. In this study, we fine-mapped the TLR7-TLR8 region and confirmed rs3853839 exhibiting the strongest association with SLE in European Americans, African Americans, and Amerindian/Hispanics. Individuals carrying the risk G allele of rs3853839 exhibited increased TLR7 expression at the both mRNA and protein level and decreased transcript degradation. MicroRNA-3148 (miR-3148) downregulated the expression of non-risk allele (C) containing transcripts preferentially, suggesting a likely mechanism for increased TLR7 levels in risk-allele carriers. This trans-ancestral mapping provides evidence for the global association with SLE risk at rs3853839, which resides in a microRNA–gene regulatory site affecting TLR7 expression.
Systemic lupus erythematosus (SLE) is an inflammatory autoimmune disease with a strong genetic component. African-Americans (AA) are at increased risk of SLE, but the genetic basis of this risk is largely unknown. To identify causal variants in SLE loci in AA, we performed admixture mapping followed by fine mapping in AA and European-Americans (EA). Through genome-wide admixture mapping in AA, we identified a strong SLE susceptibility locus at 2q22–24 (LOD = 6.28), and the admixture signal is associated with the European ancestry (ancestry risk ratio ∼1.5). Large-scale genotypic analysis on 19,726 individuals of African and European ancestry revealed three independently associated variants in the IFIH1 gene: an intronic variant, rs13023380 [Pmeta = 5.20×10−14; odds ratio, 95% confidence interval = 0.82 (0.78–0.87)], and two missense variants, rs1990760 (Ala946Thr) [Pmeta = 3.08×10−7; 0.88 (0.84–0.93)] and rs10930046 (Arg460His) [Pdom = 1.16×10−8; 0.70 (0.62–0.79)]. Both missense variants produced dramatic phenotypic changes in apoptosis and inflammation-related gene expression. We experimentally validated function of the intronic SNP by DNA electrophoresis, protein identification, and in vitro protein binding assays. DNA carrying the intronic risk allele rs13023380 showed reduced binding efficiency to a cellular protein complex including nucleolin and lupus autoantigen Ku70/80, and showed reduced transcriptional activity in vivo. Thus, in SLE patients, genetic susceptibility could create a biochemical imbalance that dysregulates nucleolin, Ku70/80, or other nucleic acid regulatory proteins. This could promote antibody hypermutation and auto-antibody generation, further destabilizing the cellular network. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.
African-Americans (AA) are at increased risk of systemic lupus erythematosus (SLE), but the genetic basis of this risk increase is largely unknown. We used admixture mapping to localize disease-causing genetic variants that differ in frequency across populations. This approach is advantageous for localizing susceptibility genes in recently admixed populations like AA. Our genome-wide admixture scan identified seven admixture signals, and we followed the best signal at 2q22–24 with fine-mapping, imputation-based association analysis and experimental validation. We identified two independent coding variants and a non-coding variant within the IFIH1 gene associated with SLE. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.
Several confirmed genetic susceptibility loci for lupus have been described. To date, no clear evidence for genetic epistasis is established in lupus. We test for gene-gene interactions in a number of known lupus susceptibility loci.
Eighteen SNPs tagging independent and confirmed lupus susceptibility loci were genotyped in a set of 4,248 lupus patients and 3,818 normal healthy controls of European descent. Epistasis was tested using a 2-step approach utilizing both parametric and non-parametric methods. The false discovery rate (FDR) method was used to correct for multiple testing.
We detected and confirmed gene-gene interactions between the HLA region and CTLA4, IRF5, and ITGAM, and between PDCD1 and IL21 in lupus patients. The most significant interaction detected by parametric analysis was between rs3131379 in the HLA region and rs231775 in CTLA4 (Interaction odds ratio=1.19, z-score= 3.95, P= 7.8×10−5 (FDR≤0.05), PMDR= 5.9×10−45). Importantly, our data suggest that in lupus patients the presence of the HLA lupus-risk alleles in rs1270942 and rs3131379 increases the odds of also carrying the lupus-risk allele in IRF5 (rs2070197) by 17% and 16%, respectively (P= 0.0028 and 0.0047).
We provide evidence for gene-gene epistasis in systemic lupus erythematosus. These findings support a role for genetic interaction contributing to the complexity of lupus heritability.
Systemic lupus erythematosus (SLE) is an autoimmune disease with diverse clinical manifestations characterized by the development of pathogenic autoantibodies manifesting in inflammation of target organs such as the kidneys, skin and joints. Genome-wide association studies have identified genetic variants in the UBE2L3 region that are associated with SLE in subjects of European and Asian ancestry. UBE2L3 encodes an ubiquitin-conjugating enzyme, UBCH7, involved in cell proliferation and immune function. In this study, we sought to further characterize the genetic association in the region of UBE2L3 and use molecular methods to determine the functional effect of the risk haplotype. We identified significant associations between variants in the region of UBE2L3 and SLE in individuals of European and Asian ancestry that exceeded a Bonferroni corrected threshold (P < 1 × 10−4). A single risk haplotype was observed in all associated populations. Individuals harboring the risk haplotype display a significant increase in both UBE2L3 mRNA expression (P = 0.0004) and UBCH7 protein expression (P = 0.0068). The results suggest that variants carried on the SLE associated UBE2L3 risk haplotype influence autoimmunity by modulating UBCH7 expression.
Systemic Lupus Erythematosus; UBE2L3; Multi Ethnic Association Study; UBCH7 Expression
Candidate gene and genome-wide association studies have identified several disease susceptibility loci in lupus patients. These studies have been largely performed in European-derived and Asian lupus patients. In this study, we examine if some of these same susceptibility loci increase lupus risk in African-American individuals.
Single nucleotide polymorphisms tagging 15 independent lupus susceptibility loci were genotyped in a set of 1,724 lupus patients and 2,024 normal 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.
We provide the first evidence for genetic association between lupus and five susceptibility loci in African-American patients (C8orf13-BLK, BANK1, TNFSF4, KIAA1542 andCTLA4; P values= 8.0 × 10−6, 1.9 × 10−5, 5.7 × 10−5, 0.00099, 0.0045, respectively). Further, we confirm the genetic association between lupus and five additional lupus susceptibility loci (ITGAM, MSH5, CFB, STAT4, and FCGR2A; P values= 7.5 × 10−11, 5.2 × 10−8, 8.7 × 10−7, 0.0058, and 0.0070, respectively), and provide evidence for a genome-wide significance for the association between ITGAM and MSH5 (HLA region) for the first time in African-American lupus patients.
These findings provide evidence for 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.
T cells from patients with SLE express increased amounts of PP2Ac which contribute to decreased production of IL-2. Because IL-2 is important in the regulation of several aspects of the immune response, it has been proposed that PP2Ac contributes to the expression of SLE. This study was designed to determine whether genetic variants of PPP2AC are linked to the expression of SLE and specific clinical manifestations and account for the increased expression of PP2Ac.
We conducted a trans-ethnic study consisting of 8,695 SLE cases and 7,308 controls from four different ancestries. Eighteen single-nucleotide polymorphisms (SNPs) across the PPP2CA were genotyped using an Illumina custom array. PPP2CA expression in SLE and control T cells was analyzed by real-time PCR.
A 32-kb haplotype comprised of multiple SNPs of PPP2CA showed significant association with SLE in Hispanic Americans (HA), European Americans (EA) and Asians but not in African-Americans (AA). Conditional analyses revealed that SNP rs7704116 in intron 1 showed consistently strong association with SLE across Asian, EA and HA populations (pmeta=3.8×10−7, OR=1.3[1.14–1.31]). In EA, the largest ethnic dataset, the risk A allele of rs7704116 was associated with the presence of renal disease, anti-dsDNA and anti-RNP antibodies. PPP2CA expression was approximately 2-fold higher in SLE patients carrying the rs7704116 AG genotype than those carrying GG genotype (p = 0.008).
Our data provide the first evidence for an association between PPP2CA polymorphisms and elevated PP2Ac transcript levels in T cells, which implicates a new molecular pathway for SLE susceptibility in EA, HA and Asians.
The overexpression of interferon (IFN)-inducible genes is a prominent feature of SLE, serves as a marker for active and more severe disease, and is also observed in other autoimmune and inflammatory conditions. The genetic variations responsible for sustained activation of IFN responsive genes are unknown.
We systematically evaluated association of SLE with a total of 1,754 IFN-pathway related genes, including IFN-inducible genes known to be differentially expressed in SLE patients and their direct regulators. We performed a three-stage design where two cohorts (total n=939 SLE cases, 3,398 controls) were analyzed independently and jointly for association with SLE, and the results were adjusted for the number of comparisons.
A total of 16,137 SNPs passed all quality control filters of which 316 demonstrated replicated association with SLE in both cohorts. Nine variants were further genotyped for confirmation in an average of 1,316 independent SLE cases and 3,215 independent controls. Association with SLE was confirmed for several genes, including the transmembrane receptor CD44 (rs507230, P = 3.98×10−12), cytokine pleiotrophin (PTN) (rs919581, P = 5.38×10−04), the heat-shock DNAJA1 (rs10971259, P = 6.31×10−03), and the nuclear import protein karyopherin alpha 1 (KPNA1) (rs6810306, P = 4.91×10−02).
This study expands the number of candidate genes associated with SLE and highlights the potential of pathway-based approaches for gene discovery. Identification of the causal alleles will help elucidate the molecular mechanisms responsible for activation of the IFN system in SLE.