Previous work has demonstrated that northern and southern European ancestries are associated with specific systemic lupus erythematosus (SLE) manifestations. Here, 1855 SLE cases of European descent were genotyped for 4965 single nucleotide polymorphisms and principal components analysis of genotype information was used to define population substructure. The first principal component (PC1) distinguished northern from southern European ancestry, PC2 differentiated eastern from western European ancestry, and PC3 delineated Ashkenazi Jewish ancestry. Compared to northern European ancestry, southern European ancestry was associated with autoantibody production (OR=1.40, 95% CI 1.07-1.83) and renal involvement (OR 1.41, 95% CI 1.06-1.87), and was protective for discoid rash (OR=0.51, 95% CI 0.32-0.82) and photosensitivity (OR=0.74, 95% CI 0.56-0.97). Both serositis (OR=1.46, 95% CI 1.12-1.89) and autoantibody production (OR=1.38, 95% CI 1.06-1.80) were associated with Western compared to Eastern European ancestry. Ashkenazi Jewish ancestry was protective against neurologic manifestations of SLE (OR=0.62, 95% CI 0.40-0.94). Homogeneous clusters of cases defined by multiple PCs demonstrated stronger phenotypic associations. Genetic ancestry may contribute to the development of SLE endophenotypes and should be accounted for in genetic studies of disease characteristics.
Systemic lupus erythematosus; epidemiology; population substructure; genetics
Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. Recent studies have greatly expanded the number of established SLE risk alleles, but the distribution of multiple risk alleles in cases versus controls and their relationship to subphenotypes have not been studied. We studied 22 SLE susceptibility polymorphisms with previous genome-wide evidence of association (p<5×10−8) in 1919 SLE cases from 9 independent Caucasian SLE case series and 4813 independent controls. The mean number of risk alleles in cases was 15.1 (SD 3.1) while the mean in controls was 13.1 (SD 2.8), with trend p = 4×10−128. We defined a genetic risk score (GRS) for SLE as the number of risk alleles with each weighted by the SLE risk odds ratio (OR). The OR for high-low GRS tertiles, adjusted for intra-European ancestry, sex, and parent study, was 4.4 (95% CI 3.8–5.1). We studied associations of individual SNPs and the GRS with clinical manifestations for the cases: age at diagnosis, the 11 American College of Rheumatology classification criteria, and double-stranded DNA antibody (anti-dsDNA) production. Six subphenotypes were significantly associated with the GRS, most notably anti-dsDNA (ORhigh-low = 2.36, p = 9e−9), the immunologic criterion (ORhigh-low = 2.23, p = 3e−7), and age at diagnosis (ORhigh-low = 1.45, p = 0.0060). Finally, we developed a subphenotype-specific GRS (sub-GRS) for each phenotype with more power to detect cumulative genetic associations. The sub-GRS was more strongly associated than any single SNP effect for 5 subphenotypes (the above plus hematologic disorder and oral ulcers), while single loci are more significantly associated with renal disease (HLA-DRB1, OR = 1.37, 95% CI 1.14–1.64) and arthritis (ITGAM, OR = 0.72, 95% CI 0.59–0.88). We did not observe significant associations for other subphenotypes, for individual loci or the sub-GRS. Thus our analysis categorizes SLE subphenotypes into three groups: those having cumulative, single, and no known genetic association with respect to the currently established SLE risk loci.
Systemic lupus erythematosus is a chronic disabling autoimmune disease, most commonly striking women in their thirties or forties. It can cause a wide variety of clinical manifestations, including kidney disease, arthritis, and skin disorders. Prognosis varies greatly depending on these clinical features, with kidney disease and related characteristics leading to greater morbidity and mortality. It is also complex genetically; while lupus runs in families, genes increase one's risk for lupus but do not fully determine the outcome. The interactions of multiple genes and/or interactions between genes and environmental factors may cause lupus, but the causes and disease pathways of this very heterogeneous disease are not well understood. By examining relationships between the presence of multiple lupus risk genes, lupus susceptibility, and clinical manifestations, we hope to better understand how lupus is triggered and by what biological pathways it progresses. We show in this work that certain clinical manifestations of lupus are highly associated with cumulative genetic variations, i.e. multiple risk alleles, while others are associated with a single variation or none at all.
Systemic Lupus Erythematosus (SLE) is an autoimmune disease with a very varied spectrum of clinical manifestations that could be partly determined by genetic factors. We aimed to determine the relationship between prevalence of 11 clinical features and age of disease onset with European population genetic substructure. Data from 1413 patients of European ancestry recruited in nine countries was tested for association with genotypes of top ancestry informative markers. This analysis was done with logistic regression between phenotypes and genotypes or principal components extracted from them. We used a genetic additive model and adjusted for gender and disease duration. Three clinical features showed association with ancestry informative markers: autoantibody production defined as immunologic disorder (P = 6.8×10−4), oral ulcers (P = 6.9×10−4) and photosensitivity (P = 0.002). Immunologic disorder was associated with genotypes more common in Southern European ancestries, whereas the opposite trend was observed for photosensitivity. Oral ulcers were specifically more common in patients of Spanish and Portuguese self-reported ancestry. These results should be taken into account in future research and suggest new hypotheses and possible underlying mechanisms to be investigated. A first hypothesis linking photosensitivity with variation in skin pigmentation is suggested.
High serum interferon α (IFNα) activity is a heritable risk factor for systemic lupus erythematosus (SLE). Auto-antibodies found in SLE form immune complexes which can stimulate IFNα production by activating endosomal Toll-like receptors and interferon regulatory factors (IRFs), including IRF5. Genetic variation in IRF5 is associated with SLE susceptibility; however, it is unclear how IRF5 functional genetic elements contribute to human disease.
1034 patients with SLE and 989 controls of European ancestry, 555 patients with SLE and 679 controls of African–American ancestry, and 73 patients with SLE of South African ancestry were genotyped at IRF5 polymorphisms, which define major haplotypes. Serum IFNα activity was measured using a functional assay.
In European ancestry subjects, anti-double-stranded DNA (dsDNA) and anti-Ro antibodies were each associated with different haplotypes characterised by a different combination of functional genetic elements (OR > 2.56, p >003C; 1.9×10−14 for both). These IRF5 haplotype-auto-antibody associations strongly predicted higher serum IFNα in patients with SLE and explained > 70% of the genetic risk of SLE due to IRF5. In African–American patients with SLE a similar relationship between serology and IFNα was observed, although the previously described European ancestry-risk haplotype was present at admixture proportions in African–American subjects and absent in African patients with SLE.
The authors define a novel risk haplotype of IRF5 that is associated with anti-dsDNA antibodies and show that risk of SLE due to IRF5 genotype is largely dependent upon particular auto-antibodies. This suggests that auto-antibodies are directly pathogenic in human SLE, resulting in increased IFNα in cooperation with particular combinations of IRF5 functional genetic elements.
SLE is a systemic autoimmune disorder affecting multiple organ systems including the skin, musculoskeletal, renal and haematopoietic systems. Humoral autoimmunity is a hallmark of SLE, and patients frequently have circulating auto-antibodies directed against dsDNA, as well as RNA binding proteins (RBP). Anti-RBP autoantibodies include antibodies which recognize Ro, La, Smith (anti-Sm), and ribonucleoprotein (anti-nRNP), collectively referred to as anti-retinol-binding protein). Anti-retinol-binding protein and anti-dsDNA auto-antibodies are rare in the healthy population.1 These auto-antibodies can be present in sera for years preceding the onset of clinical SLE illness2 and are likely pathogenic in SLE.34
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.
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.
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.
European population genetic substructure was examined in a diverse set of >1,000 individuals of European descent, each genotyped with >300 K SNPs. Both STRUCTURE and principal component analyses (PCA) showed the largest division/principal component (PC) differentiated northern from southern European ancestry. A second PC further separated Italian, Spanish, and Greek individuals from those of Ashkenazi Jewish ancestry as well as distinguishing among northern European populations. In separate analyses of northern European participants other substructure relationships were discerned showing a west to east gradient. Application of this substructure information was critical in examining a real dataset in whole genome association (WGA) analyses for rheumatoid arthritis in European Americans to reduce false positive signals. In addition, two sets of European substructure ancestry informative markers (ESAIMs) were identified that provide substantial substructure information. The results provide further insight into European population genetic substructure and show that this information can be used for improving error rates in association testing of candidate genes and in replication studies of WGA scans.
Ancestry differences corresponding to ethnic groups may be important in determining disease risk factors and optimizing treatment. Our study further defines ancestry relationship among different European ethnic groups by examining over 300 thousand variations in DNA, in over 2,000 individuals. This study allowed a clearer ascertainment of differences that could not be discerned in smaller studies using more limited numbers of DNA variations. We show clear differences among European American participants of different self-identified ethnic affiliation. The analyses showed multiple components of variation. The components showing the largest variations generally corresponded to the grandparental country or region of origin within Europe. We also show the importance of applying this information in determining genetic risk factors for complex diseases. Moreover, the results have enabled a better selection of smaller numbers of DNA variations that can be used in future disease studies to identify more homogenous participant groups and minimize false positive and false negative results in assessing genetic risk factors for disease.
Systemic lupus erythematosus (SLE) is a chronic multisystem genetically complex autoimmune disease characterised by the production of autoantibodies to nuclear and cellular antigens, tissue inflammation and organ damage. Genome-wide association studies have shown that variants within the major histocompatibility complex (MHC) region on chromosome 6 confer the greatest genetic risk for SLE in European and Chinese populations. However, the causal variants remain elusive due to tight linkage disequilibrium across disease-associated MHC haplotypes, the highly polymorphic nature of many MHC genes and the heterogeneity of the SLE phenotype.
A high-density case-control single nucleotide polymorphism (SNP) study of the MHC region was undertaken in SLE cohorts of Spanish and Filipino ancestry using a custom Illumina chip in order to fine-map association signals in these haplotypically diverse populations. In addition, comparative analyses were performed between these two datasets and a northern European UK SLE cohort. A total of 1433 cases and 1458 matched controls were examined.
Using this transancestral SNP mapping approach, novel independent loci were identified within the MHC region in UK, Spanish and Filipino patients with SLE with some evidence of interaction. These loci include HLA-DPB1, HLA-G and MSH5 which are independent of each other and HLA-DRB1 alleles. Furthermore, the established SLE-associated HLA-DRB1*15 signal was refined to an interval encompassing HLA-DRB1 and HLA-DQA1. Increased frequencies of MHC region risk alleles and haplotypes were found in the Filipino population compared with Europeans, suggesting that the greater disease burden in non-European SLE may be due in part to this phenomenon.
These data highlight the usefulness of mapping disease susceptibility loci using a transancestral approach, particularly in a region as complex as the MHC, and offer a springboard for further fine-mapping, resequencing and transcriptomic analysis.
African Americans, East Asians, and Hispanics with systemic lupus erythematosus (SLE) are more likely to develop renal disease than SLE patients of European descent. We investigated whether European genetic ancestry protects against the development of lupus nephritis and explored genetic and socioeconomic factors that might explain this effect.
This was a cross-sectional study of 1906 adults with SLE. Participants were genotyped for 126 single nucleotide polymorphisms (SNPs) informative for ancestry. A subset of participants was also genotyped for 80 SNPs in 14 candidate genes for renal disease in SLE. We used logistic regression to test the association between European ancestry and renal disease. Analyses adjusted for continental ancestries, socioeconomic status, and candidate genes.
Participants (n=1906) had on average 62.4% European, 15.8% African, 11.5% East Asian, 6.5% Amerindian, and 3.8% South Asian ancestry. Among participants, 34% (n=656) had renal disease. A 10% increase in European ancestry was associated with a 15% reduction in the odds of having renal disease after adjustment for disease duration and sex (OR 0.85, 95% CI 0.82-0.87, p=1.9 × 10−30). Adjusting for other genetic ancestries, measures of socioeconomic status, or SNPs in genes most associated with renal disease (IRF5 (rs4728142), BLK (rs2736340), STAT4 (rs3024912), ITGAM (rs9937837) and HLA-DRB1*0301 and DRB1*1501, p<0.05) did not substantively alter this relationship.
European ancestry is protective against the development of renal disease in SLE, an effect independent of other genetic ancestries, common risk alleles, and socioeconomic status.
Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. A polymorphism in the STAT4 gene has recently been established as a risk factor for SLE, but the relationship with specific SLE subphenotypes has not been studied. We studied 137 SNPs in the STAT4 region genotyped in 4 independent SLE case series (total n = 1398) and 2560 healthy controls, along with clinical data for the cases. Using conditional testing, we confirmed the most significant STAT4 haplotype for SLE risk. We then studied a SNP marking this haplotype for association with specific SLE subphenotypes, including autoantibody production, nephritis, arthritis, mucocutaneous manifestations, and age at diagnosis. To prevent possible type-I errors from population stratification, we reanalyzed the data using a subset of subjects determined to be most homogeneous based on principal components analysis of genome-wide data. We confirmed that four SNPs in very high LD (r2 = 0.94 to 0.99) were most strongly associated with SLE, and there was no compelling evidence for additional SLE risk loci in the STAT4 region. SNP rs7574865 marking this haplotype had a minor allele frequency (MAF) = 31.1% in SLE cases compared with 22.5% in controls (OR = 1.56, p = 10−16). This SNP was more strongly associated with SLE characterized by double-stranded DNA autoantibodies (MAF = 35.1%, OR = 1.86, p<10−19), nephritis (MAF = 34.3%, OR = 1.80, p<10−11), and age at diagnosis<30 years (MAF = 33.8%, OR = 1.77, p<10−13). An association with severe nephritis was even more striking (MAF = 39.2%, OR = 2.35, p<10−4 in the homogeneous subset of subjects). In contrast, STAT4 was less strongly associated with oral ulcers, a manifestation associated with milder disease. We conclude that this common polymorphism of STAT4 contributes to the phenotypic heterogeneity of SLE, predisposing specifically to more severe disease.
Systemic lupus erythematosus is a chronic disabling autoimmune disease, most commonly striking women in their thirties or forties. It can cause a wide variety of clinical manifestations, including kidney disease, arthritis, and skin disorders. Prognosis varies greatly depending on these clinical features, with kidney disease and related characteristics leading to greater morbidity and mortality. It is also complex genetically; while lupus runs in families, genes increase one’s risk for lupus but do not fully determine the outcome. It is thought that the interactions of multiple genes and/or interactions between genes and environmental factors may cause lupus, but the causes and disease pathways of this very heterogeneous disease are not well understood. By examining relationships between subtypes of lupus and specific genes, we hope to better understand how lupus is triggered and by what biological pathways it progresses. We show in this work that the STAT4 gene, very recently identified as a lupus risk gene, predisposes specifically to severe manifestations of lupus, including kidney disease.
Systemic lupus erythematosus (SLE) is a multifactorial disorder characterized by the presence of autoantibodies. We and others have implicated free radical mediated peroxidative damage in the pathogenesis of SLE. Since harmful free radical products are formed during this oxidative process, including 4-hydroxy 2-nonenol (4-HNE) and malondialdehyde (MDA), we hypothesized that specific HNE-protein adducts would be present in SLE red blood cell (RBC) membranes. Catalase is located on chromosome 11p13 where linkage analysis has revealed a marker in the same region of the genome among families with thrombocytopenia, a clinical manifestation associated with severe lupus in SLE affected pedigrees. Moreover, SLE afflicts African-Americans three times more frequently than their European-American counterparts. Hence we investigated the effects of a genetic polymorphism of catalase on risk and severity of SLE in 48 pedigrees with African American ancestry.
Tryptic digestion followed by matrix assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOFMS) analysis was used to identify the protein modified by HNE, following Coomassie staining to visualize the bands on the acrylamide gels. Genotyping analysis for the C → T, -262 bp polymorphism in the promoter region of catalase was performed by PCR-RFLP and direct PCR-sequencing. We used a "pedigree disequilibrium test" for the family based association analysis, implemented in the PDT program to analyze the genotyping results.
We found two proteins to be HNE-modified, migrating around 80 and 50 kD respectively. Tryptic digestion followed by matrix assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOFMS) analysis of the Coomassie stained 80 kD band revealed that the target of HNE modification was catalase, a protein shown to associate with RBC membrane proteins. All the test statistics carried out on the genotyping analysis for the C → T, -262 bp polymorphism in the promoter region of catalase were non-significant (p > 0.05) in our data, which suggested that this SNP is not associated with SLE.
Our results indicate that catalase is one of the proteins modified due to oxidative stress. However, catalase may not be a susceptibility gene for SLE. Nonetheless, catalase is oxidatively modified among SLE patients. This suggests a possible role between oxidative modification of catalase and its affects on enzymatic activity in SLE. An oxidatively modified catalase could be one of the reasons for lower enzymatic activity among SLE subjects, which in turn could favor the accumulation of deleterious hydrogen peroxide. Furthermore, HNE-products are potential neoantigens and could be involved in the pathogenesis of SLE. Decrease in catalase activity could affect the oxidant-antioxidant balance. Chronic disturbance of this balance in patients with SLE may work favorably for the premature onset of atherogenesis with severe vascular effect.
Systemic Lupus Erythematosus (SLE) disproportionately affects minorities, such as Hispanic-Americans. Prevalence of SLE is 3–5 times higher in Hispanic Americans (HA) than European derived populations, and have more active disease at the time of diagnosis, with more serious organ system involvement. HA is an admixed population, it is possible that there is an effect of admixture on the relative risk of disease. This admixture can create substantial increase of linkage disequilibrium (LD) in both magnitude and range, which can provide a unique opportunity for admixture mapping. Main objectives of this study are to (a) estimate hidden population structure in HA individuals; (b) estimate individual ancestry proportions and its impact on SLE risk; (c) assess impact of admixture on ITGAM association, a recently identified SLE susceptibility gene; and (d) estimate power of admixture mapping in HA. Our dataset contained 1,125 individuals, of whom 884 (657 SLE cases and 227 controls) were self classified as HA. Using 107 unlinked ancestry informative markers (AIMs) we estimated hidden population structure and individual ancestry in HA. Out of 5,671 possible pair-wise LD, 54% were statistically significant, indicating recent population admixture. The best fitted model for HA was a four population model with average ancestry of European (48%), American-Indian (40%), African (8%) and a fourth population (4%) with unknown ancestry. We also identified significant higher risk associated with American-Indian ancestry (OR=4.84, P=0.0001, 95%CI=2.14—10.95) on overall SLE. We showed that ITGAM is associated as a risk factor for SLE (OR= 2.06, P=8.74×10−5, 95%CI=1.44–2.97). This association is not affected by population substructure or admixture. We have demonstrated that HA have great potential and are an 3 appropriate population for admixture mapping. As expected, the case-only design is more powerful than case-control design, for any given admixture proportion or ancestry risk ratio.
SLE; Association; Hispanics; Admixture mapping; Hispanic-American; Population structure
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.
Systemic Lupus Erythematosus (SLE) shows a spectrum of clinical manifestations that complicate its diagnosis, treatment and research. This variability is likely related with environmental exposures and genetic factors among which known SLE susceptibility loci are prime candidates. The first published analyses seem to indicate that this is the case for some of them, but results are still inconclusive and we aimed to further explore this question.
European SLE patients, 1444, recruited at 17 centres from 10 countries were analyzed. Genotypes for 26 SLE associated SNPs were compared between patients with and without each of 11 clinical features: ten of the American College of Rheumatology (ACR) classification criteria (except ANAs) and age of disease onset. These analyses were adjusted for centre of recruitment, top ancestry informative markers, gender and time of follow-up. Overlap of samples with previous studies was excluded for assessing replication.
There were three new associations: the SNPs in XKR6 and in FAM167A-BLK were associated with lupus nephritis (OR = 0.76 and 1.30, Pcorr = 0.007 and 0.03, respectively) and the SNP of MECP2, which is in chromosome X, with earlier age of disease onset in men. The previously reported association of STAT4 with early age of disease onset was replicated. Some other results were suggestive of the presence of additional associations. Together, the association signals provided support to some previous findings and to the characterization of lupus nephritis, autoantibodies and age of disease onset as the clinical features more associated with SLE loci.
Some of the SLE loci shape the disease phenotype in addition to increase susceptibility to SLE. This influence is more prominent for some clinical features than for others. However, results are only partially consistent between studies and subphenotype specific GWAS are needed to unravel their genetic component.
Systemic lupus erythematosus (SLE) is a clinically heterogeneous, systemic autoimmune disease characterized by autoantibody formation. Previously published genome-wide association studies (GWAS) have investigated SLE as a single phenotype. Therefore, we conducted a GWAS to identify genetic factors associated with anti–dsDNA autoantibody production, a SLE–related autoantibody with diagnostic and clinical importance. Using two independent datasets, over 400,000 single nucleotide polymorphisms (SNPs) were studied in a total of 1,717 SLE cases and 4,813 healthy controls. Anti–dsDNA autoantibody positive (anti–dsDNA +, n = 811) and anti–dsDNA autoantibody negative (anti–dsDNA –, n = 906) SLE cases were compared to healthy controls and to each other to identify SNPs associated specifically with these SLE subtypes. SNPs in the previously identified SLE susceptibility loci STAT4, IRF5, ITGAM, and the major histocompatibility complex were strongly associated with anti–dsDNA + SLE. Far fewer and weaker associations were observed for anti–dsDNA – SLE. For example, rs7574865 in STAT4 had an OR for anti–dsDNA + SLE of 1.77 (95% CI 1.57–1.99, p = 2.0E-20) compared to an OR for anti–dsDNA – SLE of 1.26 (95% CI 1.12–1.41, p = 2.4E-04), with pheterogeneity<0.0005. SNPs in the SLE susceptibility loci BANK1, KIAA1542, and UBE2L3 showed evidence of association with anti–dsDNA + SLE and were not associated with anti–dsDNA – SLE. In conclusion, we identified differential genetic associations with SLE based on anti–dsDNA autoantibody production. Many previously identified SLE susceptibility loci may confer disease risk through their role in autoantibody production and be more accurately described as autoantibody propensity loci. Lack of strong SNP associations may suggest that other types of genetic variation or non-genetic factors such as environmental exposures have a greater impact on susceptibility to anti–dsDNA – SLE.
Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that can involve virtually any organ system. SLE patients produce antibodies that bind to their own cells and proteins (autoantibodies) which can cause irreversible organ damage. One particular SLE–related autoantibody directed at double-stranded DNA (anti–dsDNA) is associated with kidney involvement and more severe disease. Previous genome-wide association studies (GWAS) in SLE have studied SLE itself, not particular SLE manifestations. Therefore, we conducted this GWAS of anti–dsDNA autoantibody production to identify genetic associations with this clinically important autoantibody. We found that many previously identified SLE–associated genes are more strongly associated with anti–dsDNA autoantibody production than SLE itself, and they may be more accurately described as autoantibody propensity genes. No strong genetic associations were observed for SLE patients who do not produce anti–dsDNA autoantibodies, suggesting that other factors may have more influence in developing this type of SLE. Further investigation of these autoantibody propensity genes may lead to greater insight into the causes of autoantibody production and organ damage in SLE.
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.
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
We aimed to investigate whether the effect size of the systemic lupus erythematosus (SLE) risk alleles varies across European subpopulations.
European SLE patients (n = 1,742) and ethnically matched healthy controls (n = 2,101) were recruited at 17 centres from 10 different countries. Only individuals with self-reported ancestry from the country of origin were included. In addition, participants were genotyped for top ancestry informative markers and for 25 SLE associated SNPs. The results were used to compare effect sizes between the Central Eureopan and Southern European subgroups.
Twenty of the 25 SNPs showed independent association with SLE, These SNPs showed a significant bias to larger effect sizes in the Southern subgroup, with 15/20 showing this trend (P = 0.019) and a larger mean odds ratio of the 20 SNPs (1.46 vs. 1.34, P = 0.02) as well as a larger difference in the number of risk alleles (2.06 vs. 1.63, P = 0.027) between SLE patients and controls than for Central Europeans. This bias was reflected in a very significant difference in the cumulative genetic risk score (4.31 vs. 3.48, P = 1.8 × 10-32). Effect size bias was accompanied by a lower number of SLE risk alleles in the Southern subjects, both patients and controls, the difference being more marked between the controls (P = 1.1 × 10-8) than between the Southern and Central European patients (P = 0.016). Seven of these SNPs showed significant allele frequency clines.
Our findings showed a bias to larger effect sizes of SLE loci in the Southern Europeans relative to the Central Europeans together with clines of SLE risk allele frequencies. These results indicate the need to study risk allele clines and the implications of the polygenic model of inheritance in SLE.
Systemic lupus erythematosus (SLE) is a highly heterogeneous disorder, characterized by differences in autoantibody profile, serum cytokines, and clinical manifestations. SLE-associated autoantibodies and high serum interferon alpha (IFN-α) are important heritable phenotypes in SLE which are correlated with each other, and play a role in disease pathogenesis. These two heritable risk factors are shared between ancestral backgrounds. The aim of the study was to detect genetic factors associated with autoantibody profiles and serum IFN-α in SLE.
We undertook a case-case genome-wide association study of SLE patients stratified by ancestry and extremes of phenotype in serology and serum IFN-α. Single nucleotide polymorphisms (SNPs) in seven loci were selected for follow-up in a large independent cohort of 538 SLE patients and 522 controls using a multi-step screening approach based on novel metrics and expert database review. The seven loci were: leucine-rich repeat containing 20 (LRRC20); protein phosphatase 1 H (PPM1H); lysophosphatidic acid receptor 1 (LPAR1); ankyrin repeat and sterile alpha motif domain 1A (ANKS1A); protein tyrosine phosphatase, receptor type M (PTPRM); ephrin A5 (EFNA5); and V-set and immunoglobulin domain containing 2 (VSIG2).
SNPs in the LRRC20, PPM1H, LPAR1, ANKS1A, and VSIG2 loci each demonstrated strong association with a particular serologic profile (all odds ratios > 2.2 and P < 3.5 × 10-4). Each of these serologic profiles was associated with increased serum IFN-α. SNPs in both PTPRM and LRRC20 were associated with increased serum IFN-α independent of serologic profile (P = 2.2 × 10-6 and P = 2.6 × 10-3 respectively). None of the SNPs were strongly associated with SLE in case-control analysis, suggesting that the major impact of these variants will be upon subphenotypes in SLE.
This study demonstrates the power of using serologic and cytokine subphenotypes to elucidate genetic factors involved in complex autoimmune disease. The distinct associations observed emphasize the heterogeneity of molecular pathogenesis in SLE, and the need for stratification by subphenotypes in genetic studies. We hypothesize that these genetic variants play a role in disease manifestations and severity in SLE.
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.
To analyze if genetically determined Amerindian ancestry predicts the increased presence of risk alleles of known susceptibility genes for systemic lupus erythematosus.
Single nucleotide polymorphisms within 16 confirmed genetic susceptibility loci for SLE were genotyped in a set of 804 Mestizo lupus patients and 667 Mestizo normal healthy controls. In addition, 347 admixture informative markers were genotyped. Individual ancestry proportions were determined using STRUCTURE. Association analysis was performed using PLINK, and correlation of the presence of risk alleles with ancestry was done using linear regression.
A meta-analysis of the genetic association of the 16 SNPs across populations showed that TNFSF4, STAT4, PDCD1, ITGAM, and IRF5 were associated with lupus in a Hispanic-Mestizo cohort enriched for European and Amerindian ancestry. In addition, two SNPs within the MHC region, previously associated in a genome-wide association study in Europeans, were also associated in Mestizos. Using linear regression we predict an average increase of 2.34 risk alleles when comparing a lupus patient with 100% Amerindian ancestry to an SLE patient with 0% American Indian Ancestry (p<0.0001). SLE patients with 43% more Amerindian ancestry are predicted to carry one additional risk allele.
Amerindian ancestry increased the number of risk alleles for lupus.
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.
UBE2L3 is associated with susceptibility to systemic lupus erythematosus (SLE) and rheumatoid arthritis in European ancestry populations, and this locus has not been investigated fully in non-European populations. We studied the UBE2L3 risk allele for association with SLE, interferon-α (IFN-α), and autoantibodies in a predominantly African American SLE cohort.
We studied 395 patients with SLE and 344 controls. The UBE2L3 rs5754217 polymorphism was genotyped using Taqman primer-probe sets, and IFN-α was measured using a reporter cell assay.
The UBE2L3 rs5754217 T allele was strongly enriched in African American patients with anti-La antibodies as compared to controls, and a recessive model was the best fit for this association (OR 2.55, p = 0.0061). Serum IFN-α also demonstrated a recessive association with the rs5754217 genotype in African American patients, and the TT/anti-La-positive patients formed a significantly high IFN-α subgroup (p = 0.0040). Similar nonstatistically significant patterns of association were observed in the European American patients with SLE. Case-control analysis did not show large allele frequency differences, supporting the idea that this allele is most strongly associated with anti-La-positive patients.
This pattern of recessive influence within a subgroup of patients may explain why this allele does not produce a strong signal in standard case-control studies, and subphenotypes should be included in future studies of UBE2L3. The interaction we observed between UBE2L3 genotype and autoantibodies upon serum IFN-α suggests a biological role for this locus in patients with SLE in vivo.
SYSTEMIC LUPUS ERYTHEMATOSUS; GENETICS INTERFERON-α; AUTOANTIBODIES; UBE2L3 GENOTYPE
To evaluate FCGR3B copy number variation (CNV) in African and European populations and to determine if FCGR3B copy number is associated with SLE and SLE nephritis risk in Afro-Caribbeans, adjusting for African genetic ancestry.
We estimated FCGR3B to determine if there were ethnic variations in CNV (unrelated unadmixed Europeans and Africans). We then examined CNV at FCGR3B in relation to SLE and SLE nephritis within a case–control collection of 134 cases of SLE (37 with SLE nephritis) and 589 population controls of mainly Afro-Caribbean descent resident in Trinidad.
We found a significant difference in copy number FCGR3B distribution between unadmixed African and European UK cohorts, with 27 (29%) vs 3 (5%) for those with low (0 or 1) copy FCGR3B, respectively, P = 0.002. In a Trinidadian SLE case–control study, low FCGR3B CNV was associated with SLE risk 1.7 (95% CI 1.1, 2.8), P = 0.02, which remained after adjustment for African genetic ancestry; odds ratios (ORs) 1.7 (95% CI 1.0, 2.8), P = 0.04.
Our studies suggest that FCGR3B low copy number is associated with SLE risk in Afro-Caribbean populations independently of CNV due to African ancestry.
SLE; Admixture; FCGR3B; Copy number variation; African