PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Genes Immun. Author manuscript; available in PMC Jan 1, 2013.
Published in final edited form as:
PMCID: PMC3411915
NIHMSID: NIHMS372304
A functional haplotype of UBE2L3 confers risk for Systemic Lupus Erythematosus
Shaofeng Wang,1,29 Indra Adrianto,1,29 Graham B. Wiley,1,29 Christopher J. Lessard,1,2 Jennifer A. Kelly,1 Adam J. Adler,1 Stuart B. Glenn,1 Adrienne H. Williams,3 Julie T. Ziegler,3 Mary E. Comeau,3 Miranda C. Marion,3 Benjamin E. Wakeland,4 Chaoying Liang,4 Kenneth M. Kaufman,5,6 Joel M. Guthridge,1 Marta E. Alarcón-Riquelme,1,7 on behalf of the BIOLUPUS and GENLES Networks, Graciela S. Alarcón,8 Juan-Manuel Anaya,9 Sang-Cheol Bae,10 Jae-Hoon Kim,10 Young Bin Joo,10 Susan A. Boackle,11 Elizabeth E. Brown,12 Michelle A. Petri,13 Rosalind Ramsey-Goldman,14 John D. Reveille,15 Luis M. Vilá,16 Lindsey A. Criswell,17 Jeffrey C. Edberg,8 Barry I. Freedman,18 Gary S. Gilkeson,19 Chaim O. Jacob,20 Judith A. James,1,2 Diane L. Kamen,19 Robert P. Kimberly,8 Javier Martin,21 Joan T. Merrill,22 Timothy B. Niewold,23 Bernardo A. Pons-Estel,24 R. Hal Scofield,1,2 Anne M. Stevens,25,26 Betty P. Tsao,27 Timothy J. Vyse,28 Carl D. Langefeld,3 John B. Harley,5,6 Edward K. Wakeland,4 Kathy L. Moser,1,2 Courtney G. Montgomery,1,30 and Patrick M. Gaffney1,30
1Arthritis and Clinical Immunology Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK
2College of Medicine, University of Oklahoma Health Sciences Center, Oklahoma City, OK
3Department of Biostatistical Sciences, Wake Forest School of Medicine, Winston-Salem, NC
4Department of Immunology, University of Texas Southwestern Medical Center, Dallas, TX
5Division of Rheumatology, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH
6US Department of Veterans Affairs Medical Center, Cincinnati, OH
7Centro de Genómica e Investigaciones Oncológicas (GENyO), Pfizer-Universidad de Granada-Junta de Andalucía Granada, Spain
8Department of Medicine, University of Alabama at Birmingham, Birmingham, AL
9Center for Autoimmune Diseases Research (CREA), Universidad del Rosario, Bogotá, Colombia
10Department of Rheumatology, Hanyang University Hospital for Rheumatic Diseases, Seoul, Republic of Korea
11Division of Rheumatology, University of Colorado Denver, Aurora, CO
12Department of Epidemiology, University of Alabama at Birmingham, Birmingham, AL
13Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD
14Division of Rheumatology, Northwestern University Feinberg School of Medicine, Chicago, IL
15Rheumatology and Clinical Immunogenetics, University of Texas Health Science Center at Houston, Houston, TX
16Department of Medicine, Division of Rheumatology, University of Puerto Rico Medical Sciences Campus, San Juan, Puerto Rico
17Rosalind Russell Medical Research Center for Arthritis, University of California San Francisco, San Francisco, CA
18Department of Internal Medicine/Nephrology, Wake Forest School of Medicine, Winston-Salem, NC
19Division of Rheumatology, Medical University of South Carolina, Charleston, SC
20Department of Medicine, University of Southern California, Los Angeles, CA
21Instituto de Parasitología y Biomedicina “López-Neyra”, Consejo Superior de Investigaciones Científicas (CSIC), Granada, Spain
22Clinical Pharmacology Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK
23Section of Rheumatology and Gwen Knapp Center for Lupus and Immunology Research, University of Chicago, Chicago, IL
24Sanatorio Parque, Rosario, Argentina
25Division of Rheumatology, Department of Pediatrics, University of Washington, Seattle, WA
26Center for Immunity and Immunotherapies, Seattle Children's Research Institute, Seattle, WA
27Division of Rheumatology, Department of Medicine, University of California Los Angeles, Los Angeles, CA
28Divisions of Genetics and Molecular Medicine and Immunology, Infection and Inflammatory Diseases, King’s College London, London, UK
Corresponding Author: Patrick M. Gaffney, M.D., Arthritis and Clinical Immunology Research Program, Oklahoma Medical Research Foundation, 825 N.E. 13th Street, MS #57, Oklahoma City, OK 73104, Phone: 405-271-2572, Fax: 405-271-2536, gaffneyp/at/omrf.org
29These authors contributed equally to this work.
30These authors jointly directed this work.
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.
Keywords: Systemic Lupus Erythematosus, UBE2L3, Multi Ethnic Association Study, UBCH7 Expression
Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by self-reactive antibodies that form immune complexes leading to systemic inflammation and organ failure. SLE susceptibility is strongly influenced by both genetic and environmental factors. Recent candidate gene and genome-wide association studies (GWAS) have identified more than 30 susceptibility loci for SLE 18. Single nucleotide polymorphisms (SNPs) in the region of UBE2L3, which encodes the ubiquitin conjugating enzyme, UBCH7, demonstrate association with SLE in multiple independent SLE cohorts of European and African American ancestry 9,10 and correlate most significantly with patients developing anti-dsDNA antibodies 11. Variants in the region of UBE2L3 have also been reported to be associated with several other autoimmune disorders such as Crohn’s disease 12, 13, celiac disease 14 and rheumatoid arthritis 9, 14. Gene expression studies suggest that variants in the vicinity of UBE2L3 regulate UBE2L3 expression, thus providing a potential mechanism by which UBE2L3 influences susceptibility to autoimmune diseases 13.
Post-translational ubiquitination of proteins is an important process in eukaryotes that is responsible for the degradation of short-lived and abnormal cytosolic proteins and the regulation of cellular signaling pathways 15. Three classes of enzymes, ubiquitin-activating enzymes (E1s), ubiquitin-conjugating enzymes (E2s) and ubiquitin-protein ligases (E3s) constitute the system by which ubiquitin is transferred to target proteins. UBE2L3, located at chromosome 22q11.2, is a member of the E2 ubiquitin conjugating enzyme family and has been demonstrated to participate in the ubiquitination of p53 16, c-Fos, and the NF-κB precursor p105 in vitro 17, 18. Recent studies have further revealed that UBE2L3 is involved in cell proliferation 19.
In order to more thoroughly evaluate the UBE2L3 locus in SLE, we fine mapped and imputed SNPs in five diverse ethnic populations using a custom genotyping array, publicly available datasets of human variation and a targeted resequencing dataset enriched for subjects with SLE risk haplotypes. We identified a single 67kb risk haplotype associated with SLE and characterized the effect of the risk haplotype on gene expression by using quantitative-PCR and Western blotting. Our data demonstrate that both UBE2L3 mRNA transcripts and UBCH7 protein expression is increased by variants carried on the SLE risk haplotype, suggesting a mechanism by which variants in the region of UBE2L3 influence the pathogenesis of SLE.
Genome-wide association studies have identified genetic association with variants in the vicinity of UBE2L3 and multiple autoimmune diseases. In an effort to identify the causal variants responsible for association with SLE, we genotyped 57 SNPs in and around UBE2L3 along with 347 ancestry-informative markers (AIMs) in 8922 independent SLE cases and 8077 independent controls across five ethnic populations (Table 1, Supplementary Figure 1, Supplementary Tables 1, 2 and 3). After applying a series of quality control filters, 55 genotyped SNPs and 262 AIMs were available for further analyses. To enrich our dataset for additional untyped SNPs, we imputed a minimum of 285 SNPs from the 1000 Genomes Project. Single-marker logistic regression analyses, adjusting for gender and global ancestry estimates, revealed significant associations between multiple SNPs and SLE surpassing a Bonferroni corrected P < 1 × 10−4. In individuals of European-ancestry the strongest signal was observed at rs131658 (P = 6.50 × 10−7, odds ratio [OR] = 1.24, 95% confidence interval [CI] = 1.14–1.35, Figure 1A). In the Asian population the strongest signal occurred at rs5754177 (P = 1.98 × 10−6, OR = 1.33, 95% CI = 1.18–1.50, Figure 1B). We also observed weaker evidence of association not exceeding the Bonferroni corrected threshold in other populations with the optimal signals at rs11089629 for African Americans (P = 1.23 × 10−3, OR = 1.18, 95% CI = 1.07–1.30, Figure 1C), rs390408 for Hispanics (P = 2.89 × 10−3, OR = 1.23, 95% CI = 1.07–1.42, Figure 1D) and rs11705317 for Gullah (P = 1.74 × 10−2, OR = 0.27, 95%CI = 0.09–0.79, Figure 1E). When all populations were combined in meta-analysis, rs7444 produced the most significant association (Pcombined = 2.21 × 10−14, Supplementary Table 4) with no evidence of heterogeneity (the Cochran’s Q test P = 0.672 and the inconsistency index I2 = 0%, see Methods).
Table 1
Table 1
Samples available for analysis following quality control adjustments.
Figure 1
Figure 1
SNPs in and around the UBE2L3 region associated with SLE. (A) European-ancestry, (B) Asian, (C) African American, (D) Hispanic and (E) African-American Gullah populations. The dashed line in each panel signifies the Bonferroni corrected level of significance (more ...)
To capture novel variants enriched on the UBE2L3 risk haplotype that were not genotyped or imputed with the 1000 Genomes Project reference panel, we resequenced 174 subjects of European-ancestry enriched for SLE risk haplotypes including UBE2L3. The phased haplotypes of these sequenced individuals were then imputed into the European-ancestry dataset. This procedure added 5 novel variants (3 SNPs and 2 deletion/insertion polymorphisms [DIPs]) that were not present in dbSNP 132 (Supplementary Table 5). Among these five variants, a single base insertion located in the 3’ UTR of UBE2L3 demonstrated significant association with SLE (P = 2.56 × 10−6, OR = 1.23, 95% CI = 1.13–1.33) and is in strong linkage disequilibrium with the most significant SNP in European-ancestry (rs131658, r2 = 0.99).
To determine if differences in the linkage disequilibrium patterns across populations (trans-population mapping) could help define a minimal risk segment, we performed haplotype analysis using the thirty-four variants with P < 1 × 10−4 defined in subjects of European-ancestry (Table 2). In the European-ancestry population we observed a single 67 kb risk haplotype (P = 1.17 × 10−7) spanning the UBE2L3 region (haplotype H2, Figure 2). Similarly, a single risk haplotype harboring the majority of alleles in the EA risk haplotype was also present in Asian (haplotype H2, Supplementary Figure 2A), African American (haplotypes H2, Supplementary Figure 2B), and Hispanic populations (haplotype H2, Supplementary Figure 2C). Strong linkage disequilibrium was observed on the risk haplotype in all four populations and limited the utility of trans-population mapping or conditional analysis to further isolate a minimal risk segment. These results suggest that a single risk effect common to these populations may be responsible for the association with SLE.
Table 2
Table 2
SNPs in the region of UBE2L3 associated with SLE.
Figure 2
Figure 2
Analyses of 34 associated SNPs present on UBE2L3 region in European-ancestry population. Top: UBE2L3 haplotype association analysis with haplotype frequencies > 5%. Alleles in white boxes represent the major alleles and those in gray boxes represent (more ...)
Previous studies have demonstrated that variants in the region of UBE2L3 influence UBE2L3 transcript expression 13, therefore, we evaluated whether the SLE associated risk haplotype produced a similar molecular phenotype. To evaluate UBE2L3 mRNA and UBCH7 protein expression, quantitative real-time PCR and western blotting was performed in an independent set of EBV-transformed B cell lines under resting conditions. Cell lines were selected based on whether they contained 0, 1, or 2 copies of the UBE2L3 risk haplotype as defined by the rs7444-C risk allele. Concordant with other published studies, we observed increased UBE2L3 mRNA expression and increased expression of UBCH7 protein as a function of the number of copies of the risk haplotype (P = 0.0004 and P = 0.0068, respectively (one-way ANOVA), Figures 3A, 3B and Supplementary Figure 3).
Figure 3
Figure 3
Effect of the risk haplotype on UBE2L3 (A) mRNA and UBCH7 (B) protein expression. On the X-axis, the three different genotypes for SNP rs7444 are displayed corresponding to homozygote of risk haplotype (C/C), heterozygote (C/T), and homozygote of non-risk (more ...)
In this study, we observed significant associations between variants in UBE2L3 and SLE in individuals of European, Asian, and African-American ancestry. Weaker association evidence was also observed in the Hispanic and Gullah populations due in part to the smaller samples sizes of these two groups (Table 1). Risk variants were carried on a 67 kb risk haplotype tagged by the proxy SNP rs7444, in all populations demonstrating association with SLE. Since the variants in this haplotype block were highly correlated across the different populations, we were unable to further narrow this SLE associated DNA segment using conditional analyses or trans-population mapping.
In line with data published in Crohn’s Disease, we observed higher levels of UBE2L3 mRNA and UBCH7 protein expression in EBV cell lines carrying the risk haplotype. This suggests that similar molecular mechanisms in UBE2L3 that influence susceptibility to autoimmunity are shared between SLE and CD. The precise mechanism by which causal variants on the UBE2L3 risk haplotype influence expression of UBE2L3 is not yet defined but we hypothesize that this could be due to the modification of mRNA stability and/or modification of the binding affinity of transcription factors to the UBE2L3 promoter. Further studies geared toward identification of the causal variant(s), which underlies the effect on gene and protein expression are required.
Ubiquitination is a critical post-translational protein modification for regulation of NF-κB signaling 20, however, little is known about how UBCH7 mediated ubiquitination might impact NF-κB signaling. In a cell free system, Orian et al. demonstrated that the NF-κB precursor protein, p105, was a substrate for UBCH7 mediated ubiquitination. At rest, p105, encoded by the gene, NF-κB1, undergoes constitutive proteosomal processing to yield the NF-κB subunit, p50 18. Unprocessed p105 functions as an inhibitor of NF-κB by retaining p50 homodimers in the cytoplasm using ankyrin repeats located in the C-terminal portion of the protein 21. Following cellular activation, p105 is phosphorylated and undergoes complete proteosomal degradation, allowing bound p50 homodimers to translocate to the nucleus. It is possible that UBCH7 mediated ubiquitination of p105 may result in increased proteosomal processing and/or degradation of p105, resulting in increased levels of free p50 homodimers.
UBCH7 has been demonstrated to function with the HECT (homologous to the E6-associated protein carboxy terminus) family E3 ubiquitin ligase, ITCH, in in vitro ubiquitination assays 22, 23. ITCH participates in regulation of NF-κB along with RNF11, TAX1BP1 and A20 as part of a protein complex known at the ubiquitin-editing complex 24. Recent data demonstrates that UBCH7 is restricted to HECT and RBR (RING-in-between-RING) type E3 ligases which underscores the possibility that UBCH7 and ITCH could function together to ubiquitinate substrate proteins, however, to our knowledge, a physical interaction between ITCH and UBCH7 has not yet been demonstrated in vivo.
In summary, our data support a role for variants in the UBE2L3 locus in the predisposition to SLE in multiple ethnic populations. The UBE2L3 locus demonstrates low haplotype diversity with a single risk haplotype associated with SLE. This risk haplotype carries causal variants that result in increased expression of UBE2L3 transcripts and UBCH7 protein. Future work will now focus on the isolation and characterization of the variants that result in this expression phenotype and on the role of UBCH7 function in immune cell signaling.
Subjects
In this study, the following independent case and control subjects were collected, respectively: African-American (1,569/1,893), Asian (1,328/1,348), European-ancestry (4,248/3,818), African-American Gullah (155/131) and Hispanic enriched for the Amerindian-European admixture (1,622/887) populations (Supplementary Table 1). SLE cases were determined by meeting at least four of the eleven 1997 ACR revised criteria for SLE. Case and control samples were obtained from multiple sites with the Institutional Review Board (IRB) approval from each institution and processed at the Oklahoma Medical Research Foundation (OMRF) under the OMRF IRB.
Genotyping and Quality Control
The Illumina iSelect platform at OMRF was employed to genotype 57 SNPs and 347 genome-wide ancestry-informative markers (AIMs) 25, 26. SNP quality control (QC) measures included well-defined cluster scatter plots, a call rate >90%, a minor allele frequency >0.001 and Hardy-Weinberg proportion test p-value in controls >0.001 for inclusion. For the AIMs, we removed AIMs with low call rates (<90%), low minor allele frequencies (<0.001), and that are in LD with each other (r2>0.2). We did not perform the Hardy-Weinberg proportion test for the AIM QC to avoid AIMs being inadvertently dropped due to monomophic states in one of the ethnic groups. Principal components 27 calculated using R and global ancestry estimated using ADMIXMAP 28, 29 (with ancestral allele frequencies from African, European, American, Indian, and East Asian population) were utilized to pinpoint population outliers (Supplementary Figure 1) and to adjust the logistic regression models for controlling population structure in our association analyses. A total of 1,135 samples were removed because they were duplicates (the proportion of alleles shared identity by descent (IBD) >0.4), sample heterozygosity outliers (>5 standard deviation from the mean), population outliers, low call rate (<90%), or gender discrepancies between reported gender and genetic data (Supplementary Table 3). The final dataset, following quality control exclusions comprised 55 SNPs and 262 AIMs and 15,864 samples (Table 1).
Association Analyses
Single marker association analyses were calculated using the logistic regression function in PLINK v1.07 30 under the additive model adjusting for gender and global ancestry estimates (African, European, and East Asian). Meta-analyses to combine p-values from different populations were performed using a weighted Z-score METAL 31. We used both the Cochran's Q test statistic and I2 index to test for heterogeneity in the meta-analysis. The Cochran’s Q test calculates the weighted sum of the squared deviations between individual study effects and the overall effect across studies 32 whereas the I2 index measures the degree or percentage of inconsistency across studies due to heterogeneity rather than by chance 33. LD between variants was estimated and probable haplotypes were calculated using Haploview 4.2 34 followed by haplotypic association for all haplotypes formed by the associated markers across the various populations. Constructing haplotypes using all variants yielded the same haplotypes as the analysis using only the associated SNPs (results not shown).
Imputation
IMPUTE2 software 35 was used to impute SNPs from 20.21 Mb to 20.34 Mb on chromosome 22 with genotype data as the source of observed genotypes and the 1000 Genomes Project from Phase I interim release (June 2011) for 1,094 individuals from Africa, Asia, Europe, and the Americas (Supplementary Table 6) as reference genotypes. Imputation using the sequence data from our European-ancestry samples along with the 1000 Genome Project haplotypes was also performed. IMPUTE2 calculates posterior probabilities for the three possible genotypes (i.e. AA, AB, and BB). These probabilities were converted to the most possible genotypes with a threshold of 0.8. Imputed SNPs with the information measure less than 0.4 were excluded.
Resequencing, Variant Detection and Quality Control
We resequenced 74 SLE cases and 100 controls of European-ancestry then included the sequenced haplotypes into the genotype imputation. For each sample 3–5 ug of whole genomic DNA were sheared and prepared using an Illumina Paired-End Genomic DNA Sample Prep Kit. The SureSelect Target Enrichment System was used to enrich targeted regions of interest from each sample by utilizing a custom designed bait pool. Resequencing was performed on an Illumina GAIIx platform using standard procedures with minimum average fold coverage of 25X. Illumina Pipeline software v.1.7 was the used to process post sequence data.
Duplicate reads were excluded using a custom script followed by alignment to the human reference genome build hg18 using BWA alignment software version 0.5.9 36. Realignment of reads around insertion/deletion sites and problematic areas, base quality score recalibration, and variation detection were processed using the Genome Analysis Tool Kit (GATK) software suite version 1.0 37, 38. Variants clustered within 10 base pairs were filtered out, as well as any variant with a quality score less than 30, a quality by depth score less than 5, inclusion within a homopolymer run of 5 or more bases, or a strand bias score of greater than −0.1. The program Beagle version 3.3 39 was utilized to determine variant phase. PLINK and IMPUTE2 format files were created using the vcftools software suite version 0.1.3 40.
In order to assess the quality of the sequence data, the sequence-based variant calls were compared with common SNPs previously genotyped with the Illumina iSelect platform. More than 99% concordance was observed suggesting high quality of our sequence data. Samples with more than 5% of variants inconsistent with genotype calls required a manual inspection of the assembled contig sequence to determine the sequence quality using the Integrative Genomics Viewer (IGV) program 41. The assembled contig sequence of each novel variant identified by our sequencing was also inspected using IGV.
Cell Culture
EBV-transformed B cell lines were requested from the Lupus Family Registry and Repository (LFRR) at OMRF with IRB approval. All cell lines in this study were EA samples and were stratified by rs7444 genotype, which is a proxy of the UBE2L3 risk haplotype. Cell lines are either homozygous (carry two copies) of non-risk haplotype, heterozygous (one copy of the risk haplotype and one copy of non risk haplotype), or homozygous (carry two copies) of risk haplotype. Cell lines were cultured in RPMI 1640 supplemented with 10% fetal bovine serum, penicillin, streptomycin, L-glutamine, and 55µM beta-mercaptoethanol. Equal numbers of cells were harvested under basal culture condition in log-phase growth.
RNA Isolation and Quantitative RT-PCR
Total RNA was isolated using the Trizol total RNA isolation reagent (Invitrogen Inc., Carlsbad, CA). The concentrations of total RNA were determined by using nanodrop, and were diluted with 20ng/µL of MS2-RNA (Hoffmann-La Roche, Inc., Nutley NJ) to a final concentration of 0.5µg/µL. Total RNA was treated with DNase and cDNA was synthesized using the iScript cDNA Synthesis Kits purchased from Bio-Rad Laboratories, Inc., Hercules, CA. Quantitative PCR was carried out using the SYBR Green method to determine the mRNA expression of UBE2L3. A pair of primers was designed and synthesized: sense, 5'-TTAGTGCCGAAAACTGGAAGC-3'; anti-sense, 5'-ATTCACCAGTGCTATGAGGGAC-3'. The PCR product corresponds to 346bp-416bp of UBE2L3 mRNA. Human HMBS gene was used in quantitative RT-PCR as a reference. The RT2 qPCR Primer Assay-SYBR Green Human HMBS Kit was purchased from SABiosciences Inc., Frederick, MD. mRNA expression of UBE2L3 was normalized to HMBS.
UBE2L3 Protein Expression
EBV-transformed B cells were harvested and lysed in Whole Cell Extraction Buffer (25mM Tris, 1% Triton X-100, 150mM NaCl, 1mM EDTA and protease inhibitors). Concentrations of protein in each cell line were determined using Quick Start Bradford Protein Assay Kits and were adjusted to a final protein concentration of 2mg/mL. Anti-UBE2L3 and Anti-GAPDH antibodies were purchased from Cell Signaling Technology, Inc., Danvers, MA, and were used to detect protein expression of UBCH7 and GAPDH, respectively. ECL Plus Western Blotting Detection System was purchased from GE Heathcare, Inc., Amersham, UK. The intensity of each band was analyzed using Image J (NIH) software. Protein expression of UBCH7 was normalized to GAPDH.
Supplementary Material
Supplementary Info
ACKNOWLEDGMENTS
We would like to thank all individuals including the SLE patients and the controls that participated in this study. We are grateful to the research assistants, coordinators and physicians that helped in the recruitment of participants. We would like to express our gratitude the following individuals for contributing samples genotyped in this study: S. D'Alfonso (Italy), R. Scorza (Italy), P. Junker and H. Laustrup (Denmark), M. Bijl (Holland), E. Endreffy (Hungary), C. Vasconcelos and B.M. da Silva (Portugal), A. Suarez and C. Gutierrez (Spain), I. Rúa-Figueroa (Spain) and C. Garcilazo (Argentina). For the Asociación Andaluza de Enfermedades Autoimmunes (AADEA) collaboration: N. Ortego-Centeno (Spain), J. Jimenez-Alonso (Spain), E. de Ramon (Spain) and J. Sanchez-Roman (Spain). For the collaboration on Hispanic populations enriched for Amerindian-European admixture: M. Cardiel (Mexico), I.G. de la Torre (Mexico), M. Maradiaga (Mexico), J.F. Moctezuma (Mexico), E. Acevedo (Peru), C. Castel and M. Busajm (Argentina), and J. Musuruana (Argentina). Other participants from the Argentine Collaborative Group are: H.R. Scherbarth, P.C. Marino, E.L. Motta, S. Gamron, C. Drenkard, E. Menso, A. Allievi, G.A. Tate, J.L. Presas, S.A. Palatnik, M. Abdala, M. Bearzotti, A. Alvarellos, F. Caeiro, A. Bertoli, S. Paira, S. Roverano, C.E. Graf, E. Bertero, C. Guillerón, S. Grimaudo, J. Manni, L.J. Catoggio, E.R. Soriano, C.D. Santos, C. Prigione, F.A. Ramos, S.M. Navarro, G.A. Berbotto, M. Jorfen, E.J. Romero, M.A. Garcia, J.C. Marcos, A.I. Marcos, C.E. Perandones, A. Eimon and C.G. Battagliotti.
We would like to thank P.S. Ramos and S. Frank for their assistance in genotyping, quality control analyses and clinical data management and the staff of the Lupus Family Registry and Repository (LFRR) for collecting and maintaining SLE samples. Support for this work was obtained from the US National Institutes of Health grants: R01 AI063274, R01 AR056360, P20 GM103456 (P.M.G.); R01 AR043274 (K.L.M.); N01 AR62277 (K.L.M. and J.B.H.); R37 24717, R01 AR042460, P01 AI083194, R01 DE018209 (J.B.H.); P01 AR49084 (R.P.K., J.C.E. and E.E.B); 5UL1 RR025777 (J.C.E.); R01 AR33062 (R.P.K.); P30 AR48311 (E.E.B); K08 AI083790, LRP AI071651, UL1 RR024999 (T.B.N.); R01 CA141700, RC1 AR058621 (M.E.A.R.); R01 AR051545-01A2, ULI RR025014-02 (A.M.S.); U19 AI082714, RC1 AR058554, P30 RR031152, P30 AR053483, N01 AI50026 (J.A.J. and J.M.G.); R21 AI070304 (S.A.B.); R01 AR43814 (B.P.T.); P60 AR053308, M01 RR-00079 (L.A.C.); R01 AR043727, UL1 RR025005 (M.A.P.), K24 AR002138, P60 2 AR30692, P01 AR49084, UL1 RR025741 (R.R.G.), 1U54 RR23417-01 (J.D.R.), R01 AR043727, UL1 RR025005 (M.A.P.), P60 AR049459 and UL1 RR029882 (D.L.K.). Additional support was obtained from the Alliance for Lupus Research (K.L.M.); Merit Award from the US Department of Veterans Affairs (J.B.H. and G.S.G.); the Swedish Research Council for Medicine, Gustaf Vth-80th Jubilee Fund and Swedish Association Against Rheumatism, Instituto de Salud Carlos III, Oklahoma Center for Advancement of Science and Technology (OCAST) HR09-106 (M.E.A.R.); the European Science Foundation funds the BIOLUPUS network (M.E.A.R. coordinator); the Barrett Scholarship Fund Oklahoma Medical Research Foundation (OMRF) (C.J.L.); the Korea Healthcare Technology Research and Development Project, Ministry for Health and Welfare, Republic of Korea (A111218-11-GM01, S.C.B.); Lupus Research Institute (T.B.N.); The Alliance for Lupus Research (T.B.N., L.A.C. and C.O.J.); the Arthritis National Research Foundation Eng Tan Scholar Award (T.B.N.); Arthritis Foundation (P.M.G. and A.M.S.); the Lupus Foundation of Minnesota (P.M.G. and K.L.M.); the Wellcome Trust (T.J.V.); Arthritis Research UK (T.J.V.); Kirkland Scholar Award (L.A.C.); Wake Forest University Health Sciences Center for Public Health Genomics (C.D.L.); and the Federico Wilhelm Agricola Foundation Research grant (B.A.P.E.) . The work reported on in this publication has been in part financially supported by the ESF, in the framework of the Research Networking Programme European Science Foundation - The Identification of Novel Genes and Biomarkers for Systemic Lupus Erythematosus (BIOLUPUS)” 07-RNP-083.
Footnotes
AUTHOR CONTRIBUTIONS
P.M.G., C.G.M., K.L.M., C.J.L., J.A.K., K.M.K., C.D.L., and J.B.H. selected SNPs and were responsible for the study design. J.M.G., M.E.A.R., G.S.A., J.M.A., S.C.B., S.A.B., E.E.B., M.A.P., R.R.G., J.D.R., L.M.V., L.A.C., J.C.E., B.I.F., G.S.G., C.O.J., J.A.J., D.L.K., R.P.K., J.M., J.T.M., T.B.N., B.A.P.E., A.M.S., B.P.T., L.M.V., T.J.V., J.B.H., K.L.M., and P.M.G. assisted in the collection and characterization of the SLE cases and controls. A.A., K.M.K. and P.M.G. performed the genotyping. S.B.G., A.W., J.Z., M.E.C., M.M., J.A.K., K.M.K. and C.D.L. performed genotyping quality control. I.A., S.W., G.W., C.J.L., C.G.M. and P.M.G. performed association analyses and imputation. G.W., B.E.W., C.L., E.K.W., S.W. and P.M.G. performed sequencing. G.W., I.A., S.W., C.J.L., S.B.G., C.G.M. and P.M.G. performed sequencing data analysis. S.W. and P.M.G. performed functional studies. S.W., I.A., G.W., C.G.M., and P.M.G. prepared the manuscript and all authors approved the final draft.
COMPETING FINANCIAL INTERESTS
The authors declare no competing financial interests.
1. Kozyrev SV, Abelson AK, Wojcik J, Zaghlool A, Linga Reddy MV, Sanchez E, et al. Functional variants in the B-cell gene BANK1 are associated with systemic lupus erythematosus. Nat Genet. 2008;40(2):211–216. [PubMed]
2. Gateva V, Sandling JK, Hom G, Taylor KE, Chung SA, Sun X, et al. A large-scale replication study identifies TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 as risk loci for systemic lupus erythematosus. Nat Genet. 2009;41(11):1228–1233. [PMC free article] [PubMed]
3. Han JW, Zheng HF, Cui Y, Sun LD, Ye DQ, Hu Z, et al. Genome-wide association study in a Chinese Han population identifies nine new susceptibility loci for systemic lupus erythematosus. Nat Genet. 2009;41(11):1234–1237. [PubMed]
4. Harley JB, Alarcon-Riquelme ME, Criswell LA, Jacob CO, Kimberly RP, Moser KL, et al. Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM, PXK, KIAA1542 and other loci. Nat Genet. 2008;40(2):204–210. [PMC free article] [PubMed]
5. Hom G, Graham RR, Modrek B, Taylor KE, Ortmann W, Garnier S, et al. Association of systemic lupus erythematosus with C8orf13-BLK and ITGAM-ITGAX. N Engl J Med. 2008;358(9):900–909. [PubMed]
6. Deng Y, Tsao BP. Genetic susceptibility to systemic lupus erythematosus in the genomic era. Nat Rev Rheumatol. 2010;6(12):683–692. [PMC free article] [PubMed]
7. Graham RR, Cotsapas C, Davies L, Hackett R, Lessard CJ, Leon JM, et al. Genetic variants near TNFAIP3 on 6q23 are associated with systemic lupus erythematosus. Nat Genet. 2008;40(9):1059–1061. [PMC free article] [PubMed]
8. Adrianto I, Wen F, Templeton A, Wiley G, King JB, Lessard CJ, et al. Association of a functional variant downstream of TNFAIP3 with systemic lupus erythematosus. Nat Genet. 2011;43(3):253–258. [PMC free article] [PubMed]
9. Orozco G, Eyre S, Hinks A, Bowes J, Morgan AW, Wilson AG, et al. Study of the common genetic background for rheumatoid arthritis and systemic lupus erythematosus. Annals of the rheumatic diseases. 2011;70(3):463–468. [PMC free article] [PubMed]
10. Agik S, Franek BS, Kumar AA, Kumabe M, Utset TO, Mikolaitis RA, et al. The autoimmune disease risk allele of UBE2L3 in African American patients with systemic lupus erythematosus: a recessive effect upon subphenotypes. J Rheumatol. 2012;39(1):73–78. [PMC free article] [PubMed]
11. Chung SA, Taylor KE, Graham RR, Nititham J, Lee AT, Ortmann WA, et al. Differential genetic associations for systemic lupus erythematosus based on anti-dsDNA autoantibody production. PLoS genetics. 2011;7(3):e1001323. [PMC free article] [PubMed]
12. Franke A, McGovern DP, Barrett JC, Wang K, Radford-Smith GL, Ahmad T, et al. Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci. Nature genetics. 2010;42(12):1118–1125. [PMC free article] [PubMed]
13. Fransen K, Visschedijk MC, van Sommeren S, Fu JY, Franke L, Festen EA, et al. Analysis of SNPs with an effect on gene expression identifies UBE2L3 and BCL3 as potential new risk genes for Crohn's disease. Hum Mol Genet. 2010;19(17):3482–3488. [PubMed]
14. Zhernakova A, Stahl EA, Trynka G, Raychaudhuri S, Festen EA, Franke L, et al. Meta-analysis of genome-wide association studies in celiac disease and rheumatoid arthritis identifies fourteen non-HLA shared loci. PLoS genetics. 2011;7(2):e1002004. [PMC free article] [PubMed]
15. Ciechanover A. The ubiquitin-proteasome pathway: on protein death and cell life. The EMBO journal. 1998;17(24):7151–7160. [PubMed]
16. Ciechanover A, Shkedy D, Oren M, Bercovich B. Degradation of the tumor suppressor protein p53 by the ubiquitin-mediated proteolytic system requires a novel species of ubiquitin-carrier protein, E2. The Journal of biological chemistry. 1994;269(13):9582–9589. [PubMed]
17. Moynihan TP, Ardley HC, Leek JP, Thompson J, Brindle NS, Markham AF, et al. Characterization of a human ubiquitin-conjugating enzyme gene UBE2L3. Mamm Genome. 1996;7(7):520–525. [PubMed]
18. Orian A, Whiteside S, Israel A, Stancovski I, Schwartz AL, Ciechanover A. Ubiquitin-mediated processing of NF-kappa B transcriptional activator precursor p105. Reconstitution of a cell-free system and identification of the ubiquitin-carrier protein, E2, and a novel ubiquitin-protein ligase, E3, involved in conjugation. The Journal of biological chemistry. 1995;270(37):21707–21714. [PubMed]
19. Whitcomb EA, Taylor A. Ubiquitin control of S phase: a new role for the ubiquitin conjugating enzyme, UbcH7. Cell Div. 2009;4:17. [PMC free article] [PubMed]
20. Wertz IE, Dixit VM. Signaling to NF-kappaB: regulation by ubiquitination. Cold Spring Harb Perspect Biol. 2010;2(3):a003350. [PMC free article] [PubMed]
21. Hayden MS, Ghosh S. Shared principles in NF-kappaB signaling. Cell. 2008;132(3):344–362. [PubMed]
22. Raimondo D, Giorgetti A, Bernassola F, Melino G, Tramontano A. Modelling and molecular dynamics of the interaction between the E3 ubiquitin ligase Itch and the E2 UbcH7. Biochem Pharmacol. 2008;76(11):1620–1627. [PubMed]
23. Rossi M, De Laurenzi V, Munarriz E, Green DR, Liu YC, Vousden KH, et al. The ubiquitin-protein ligase Itch regulates p73 stability. The EMBO journal. 2005;24(4):836–848. [PubMed]
24. Markson G, Kiel C, Hyde R, Brown S, Charalabous P, Bremm A, et al. Analysis of the human E2 ubiquitin conjugating enzyme protein interaction network. Genome research. 2009;19(10):1905–1911. [PubMed]
25. Smith MW, Patterson N, Lautenberger JA, Truelove AL, McDonald GJ, Waliszewska A, et al. A high-density admixture map for disease gene discovery in african americans. Am J Hum Genet. 2004;74(5):1001–1013. [PubMed]
26. Halder I, Shriver M, Thomas M, Fernandez JR, Frudakis T. A panel of ancestry informative markers for estimating individual biogeographical ancestry and admixture from four continents: utility and applications. Hum Mutat. 2008;29(5):648–658. [PubMed]
27. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat.Genet. 2006;38(8):904–909. [PubMed]
28. Hoggart CJ, Parra EJ, Shriver MD, Bonilla C, Kittles RA, Clayton DG, et al. Control of confounding of genetic associations in stratified populations. Am J Hum Genet. 2003;72(6):1492–1504. [PubMed]
29. Hoggart CJ, Shriver MD, Kittles RA, Clayton DG, McKeigue PM. Design and analysis of admixture mapping studies. Am J Hum Genet. 2004;74(5):965–978. [PubMed]
30. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–575. [PubMed]
31. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26(17):2190–2191. [PMC free article] [PubMed]
32. Cochran WG. The Combination of Estimates from Different Experiments. Biometrics. 1954;10(1):101–129.
33. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003;327(7414):557–560. [PMC free article] [PubMed]
34. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–265. [PubMed]
35. Howie BN, Donnelly P, Marchini J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. 2009;5(6):e1000529. [PMC free article] [PubMed]
36. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–1760. [PMC free article] [PubMed]
37. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–1303. [PubMed]
38. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43(5):491–498. [PMC free article] [PubMed]
39. Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007;81(5):1084–1097. [PubMed]
40. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–2158. [PMC free article] [PubMed]
41. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nature biotechnology. 2011;29(1):24–26. [PMC free article] [PubMed]