Currently, it remains difficult to identify which single nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWAS) are functional and how various functional SNPs (fSNPs) interact and contribute to disease susceptibility. GWAS have identified a CD40 locus that is associated with rheumatoid arthritis (RA). We previously used two techniques developed in our laboratory, single nucleotide polymorphism-next-generation sequencing (SNP-seq) and flanking restriction enhanced DNA pulldown-mass spectrometry (FREP-MS), to determine that the RA risk gene RBPJ regulates CD40 expression via a fSNP at the RA-associated CD40 locus. In the present work, by applying the same approach, we report the identification of six proteins that regulate RBPJ expression via binding to two fSNPs on the RA-associated RBPJ locus. Using these findings, together with the published data, we constructed an RA-associated signal transduction and transcriptional regulation network (STTRN) that functionally connects multiple RA-associated risk genes via transcriptional regulation networks (TRNs) linked by CD40-induced nuclear factor kappa B (NF-kB) signaling. Remarkably, this STTRN provides insight into the potential mechanism of action for the histone deacetylase inhibitor givinostat, an approved therapy for systemic juvenile idiopathic arthritis. Thus, the generation of disease-associated STTRNs based on post-GWAS functional studies is demonstrated as a novel and effective approach to apply GWAS for mechanistic studies and target identification.

Previously, using FREP-MS, we identified a protein complex including eight proteins that specifically bind to the functional SNP (fSNP) rs6032664 at a CD40 locus associated with autoimmune diseases. Among these eight proteins, four are ribosomal proteins RPL26, RPL4, RPL8, and RPS9 that normally make up the ribosomal subunits involved in the cellular process of protein translation. So far, no publication has shown these ribosomal proteins function as transcriptional regulators. In this work, we demonstrate that four ribosomal proteins: RPL26, RPL4, RPL8, and RPS9 are bona fide CD40 transcriptional regulators via binding to rs6032664. In addition, we show that suppression of CD40 expression by RPL26 RNAi knockdown inactivates NF-κB p65 by dephosphorylation via NF-κB signaling pathway in fibroblast-like synoviocytes (FLS), which further reduces the transcription of disease-associated risk genes such as STAT4, CD86, TRAF1 and ICAM1 as the direct targets of NF-κB p65. Based on these findings, a disease-associated risk gene transcriptional regulation network (TRN) is generated, in which decreased expression of, at least, RPL26 results in the downregulation of risk genes: STAT4, CD86, TRAF1 and ICAM1, as well as the two proinflammatory cytokines: IL1β and IL6 via CD40-induced NF-κB signaling. We believe that further characterization of this disease-associated TRN in the CD40-induced NF-κB signaling by identifying both the upstream and downstream regulators will potentially enable us to identify the best targets for drug development.