Molecular effects of autoimmune-risk promoter polymorphisms on expression, exon choice, and translational efficiency of interferon regulatory factor 5

Characterization of Rheumatoid Arthritis Risk-Associated SNPs and Identification of Novel Therapeutic Sites Using an In-Silico Approach

Simple Summary

Rheumatoid arthritis (RA) is a complex disease resulting from multiple genetic and environmental pathogenic factors. The genetic factors include single-nucleotide polymorphisms (SNPs), which have been reported to be associated with RA, but their specific role in the pathogenesis of RA remains unexplained. This study explains the potential role of RA risk-associated SNPs in its pathogenesis in order to provide a basis for understanding the genetic complexity of RA. Several roles of these SNPs are described in this study, and may also aid in the design of a therapeutic strategy for RA. Furthermore, novel potential therapeutic sites have also been researched, resulting in the identification of three novel therapeutic targets. The therapeutic strategies for the treatment of RA include inflammatory pathway-targeting drugs, which alleviate inflammation in joints. There is always a need for novel therapeutic targets that can play a role in alleviating inflammation in autoimmune diseases including RA. Therefore, these novel therapeutic sites are very important, and further experimental studies are required.

Abstract

Single-nucleotide polymorphisms (SNPs) are reported to be associated with many diseases, including autoimmune diseases. In rheumatoid arthritis (RA), about 152 SNPs are reported to account for ~15% of its heritability. These SNPs may result in the alteration of gene expression and may also affect the stability of mRNA, resulting in diseased protein. Therefore, in order to predict the underlying mechanism of these SNPs and identify novel therapeutic sites for the treatment of RA, several bioinformatics tools were used. The damaging effect of 23 non-synonymous SNPs on proteins using different tools suggested four SNPs, including rs2476601 in PTPN22, rs5029941 and rs2230926 in TNFAIP3, and rs34536443 in TYK2, to be the most damaging. In total, 42 of 76 RA-associated intronic SNPs were predicted to create or abolish potential splice sites. Moreover, the analysis of 11 RA-associated UTR SNPs indicated that only one SNP, rs1128334, located in 3′UTR of ETS1, caused functional pattern changes in BRD-BOX. For the identification of novel therapeutics sites to treat RA, extensive gene–gene interaction network interactive pathways were established, with the identification of 13 potential target sites for the development of RA drugs, including three novel target genes. The anticipated effect of these findings on RA pathogenesis may be further validated in both in vivo and in vitro studies.

IRF5-mediated immune responses and its implications in immunological disorders

Transcription factors are gene regulators that activate or repress target genes. One family of the transcription factors that have been extensively studied for their crucial role in regulating gene network in the immune system is the interferon regulatory factors (IRFs). IRFs possess a novel turn-helix turn motif that recognizes a specific DNA consensus found in the promoters of many genes that are involved in immune responses. IRF5, a member of IRFs has recently gained much attention for its role in regulating inflammatory responses and autoimmune diseases. Here, we discuss the role of IRF5 in regulating immune cells functions and how the dysregulation of IRF5 contributes to the pathogenesis of immune disorders. We also review the latest findings of potential IRF5 inhibitors that modulate IRF5 activity in the effort of developing therapeutic approaches for treating inflammatory disorders.

Lineage-Specific Functionality of an Interferon Regulatory Factor 5 Lupus Risk Haplotype: Lack of B Cell Intrinsic Effects

Interferon regulatory factor 5 (IRF5) is widely recognized as a risk locus for systemic lupus erythematosus (SLE). Risk gene and IRF5 activation is triggered through toll-like receptor signaling. In myeloid cells, this leads to production of type I interferon and inflammatory cytokines, with enhanced production in cells of individuals harboring IRF5 risk alleles. Mouse models have also demonstrated the importance of IRF5 in B cell function, particularly plasma cell differentiation and isotype switching. Here, we evaluated the major SLE risk haplotype of IRF5 on the functional attributes of freshly isolated B cells from human subjects who do not have evidence of SLE or other forms of autoimmunity. We took this approach to avoid the complications of studying genotype-phenotype relationships in B cells that have been chronically exposed to an inflammatory disease environment before isolation. We focused on B cell endophenotypes that included gene expression, antibody secretion, class switching, and apoptotic susceptibility. We performed IRF5 overexpression studies, genetic reporter assays and electro-mobility shift assays on B and myeloid cell lines. Somewhat surprisingly, the results of our analyses indicate that IRF5 risk genotypes do not have a B cell intrinsic effect on these B cell functions. By contrast, we confirmed that the IRF5 risk and non-risk haplotypes exert differential effects in myeloid cells, including an increased susceptibility to apoptosis conferred by the risk haplotype. We also demonstrated an increased binding of the transcription factor specificity protein 1 to an insertion/deletion present in the risk haplotype. Our findings raise the specter that genetic risk alleles can have complex and unexpected lineage-specific effects, and these must be carefully considered when guiding or developing therapies based on understanding disease risk haplotypes.

Mangiferin inhibits macrophage classical activation via downregulating interferon regulatory factor 5 expression

Mangiferin is a natural polyphenol and the predominant effective component of Mangifera indica Linn. leaves. For hundreds of years, Mangifera indica Linn. leaf has been used as an ingredient in numerous traditional Chinese medicine preparations for the treatment of bronchitis. However, the pharmacological mechanism of mangiferin in the treatment of bronchitis remains to be elucidated. Macrophage classical activation is important role in the process of bronchial airway inflammation, and interferon regulatory factor 5 (IRF5) has been identified as a key regulatory factor for macrophage classical activation. The present study used the THP-1 human monocyte cell line to investigate whether mangiferin inhibits macrophage classical activation via suppressing IRF5 expression in vitro. THP-1 cells were differentiated to macrophages by phorbol 12-myristate 13-acetate. Macrophages were polarized to M1 macrophages following stimulation with lipopolysaccharide (LPS)/interferon-γ (IFN-γ). Flow cytometric analysis was conducted to detect the M1 macrophages. Reverse transcription-quantitative polymerase chain reaction was used to investigate cellular IRF5 gene expression. Levels of proinflammatory cytokines and IRF5 were assessed following cell culture and cellular homogenization using enzyme-linked immunosorbent assay. IRF5 protein and nuclei co-localization was performed in macrophages with laser scanning confocal microscope immunofluorescence analysis. The results of the present study demonstrated that mangiferin significantly inhibits LPS/IFN-γ stimulation-induced classical activation of macrophages in vitro and markedly decreases proinflammatory cytokine release. In addition, cellular IRF5 expression was markedly downregulated. These results suggest that the inhibitory effect of mangiferin on classical activation of macrophages may be exerted via downregulation of cellular IRF5 expression levels.

Discover hidden splicing variations by mapping personal transcriptomes to personal genomes

RNA-seq has become a popular technology for studying genetic variation of pre-mRNA alternative splicing. Commonly used RNA-seq aligners rely on the consensus splice site dinucleotide motifs to map reads across splice junctions. Consequently, genomic variants that create novel splice site dinucleotides may produce splice junction RNA-seq reads that cannot be mapped to the reference genome. We developed and evaluated an approach to identify ‘hidden’ splicing variations in personal transcriptomes, by mapping personal RNA-seq data to personal genomes. Computational analysis and experimental validation indicate that this approach identifies personal specific splice junctions at a low false positive rate. Applying this approach to an RNA-seq data set of 75 individuals, we identified 506 personal specific splice junctions, among which 437 were novel splice junctions not documented in current human transcript annotations. 94 splice junctions had splice site SNPs associated with GWAS signals of human traits and diseases. These involve genes whose splicing variations have been implicated in diseases (such as OAS1), as well as novel associations between alternative splicing and diseases (such as ICA1). Collectively, our work demonstrates that the personal genome approach to RNA-seq read alignment enables the discovery of a large but previously unknown catalog of splicing variations in human populations.