RESUMEN
To investigate the molecular mechanisms through which Epstein-Barr virus (EBV) may contribute to Systemic Lupus Erythematosus (SLE) pathogenesis, we interrogated SLE genetic risk loci for signatures of EBV infection. We first compared the gene expression profile of SLE risk genes across 459 different cell/tissue types. EBV-infected B cells (LCLs) had the strongest representation of highly expressed SLE risk genes. By determining an SLE risk allele effect on gene expression (expression quantitative trait loci, eQTL) in LCLs and 16 other immune cell types, we identified 79 SLE risk locus:gene pairs putatively interacting with EBV infection. A total of 10 SLE risk genes from this list (CD40, LYST, JAZF1, IRF5, BLK, IKZF2, IL12RB2, FAM167A, PTPRC and SLC15A) were targeted by the EBV transcription factor, EBNA2, differentially expressed between LCLs and B cells, and the majority were also associated with EBV DNA copy number, and expression level of EBV encoded genes. Our final gene network model based on these genes is suggestive of a nexus involving SLE risk loci and EBV latency III and B cell proliferation signalling pathways. Collectively, our findings provide further evidence to support the interaction between SLE risk loci and EBV infection that is in part mediated by EBNA2. This interplay may increase the tendency towards EBV lytic switching dependent on the presence of SLE risk alleles. These results support further investigation into targeting EBV as a therapeutic strategy for SLE.