Background: Systemic lupus erythematosus (SLE) is an autoimmune inflammatory rheumatic disease. SLE susceptibility is affected by multiple genetic elements, environmental factors, and their interactions. We aimed in this study to statistically and functionally characterize a gene-gene interaction (epistasis) recently documented to affect SLE risk. Methods: Two single-nucleotide polymorphisms, rs2230926 in TNFAIP3 (A20) gene and rs131654 in UBE2L3 (UBCH7) gene, were genotyped in all 3525 Korean participants, and their SLE risk association and epistasis were statistically analyzed by calculating odds ratio (OR), 95% confidence interval (CI), and P values in genotype comparisons between 1318 SLE patients and 2207 healthy controls. Furthermore, their effects on gene functions were assessed by comparatively examining separate and combined effects of TNFAIP3 and UBE2L3 knockdowns on NF-κB transcription factor activity in human cells. Results: SLE susceptibility is associated with TNFAIP3 rs2230926 (OR = 1.9, 95% CI 1.6-2.4, P = 8.6 × 10 −11) and UBE2L3 rs131654 (OR = 1.2, 95% CI 1.1-1.4, P = 1.1 × 10 −4) in a Korean population of this study. Their risk-associated alleles synergistically elevate SLE susceptibility in both multivariate logistic regression analysis (OR interaction = 1.6, P = 0.0028) and genotype-stratified analysis (OR interaction = 2.4), confirming the synergistic TNFAIP3-UBE2L3 interaction in SLE risk. Additionally, the SLE-susceptible alleles confer decreased TNFAIP3 expression (P = 1.1 × 10 −6 , n = 610) and increased UBE2L3 expression (P = 9.5 × 10 −11 , n = 475), respectively, in B cell analysis of the International HapMap Project individuals with adjustment for ethnicity. Furthermore, when compared with TNFAIP3 non-knockdown and UBE2L3 knockdown in human HeLa cells, TNFAIP3 knockdown and UBE2L3 non-knockdown synergistically increase three cytokines, CCL2, CXCL8 (IL8), and IL6, all regulated by NF-κB in the human TNFR signaling pathway.