Cis-regulation of IRF5 expression is unable to fully account for systemic lupus erythematosus association: analysis of multiple experiments with lymphoblastoid cell lines

Coordination between innate immune cells, type I IFNs and IRF5 drives SLE pathogenesis

Systemic lupus erythematosus (SLE) is a complex autoimmune disease which affects multiple organs. The type I interferon (IFN) gene signature and circulating autoantibodies are hallmarks of SLE. Plasmacytoid dendritic cells (pDCs) are considered the main producers of type I IFN and production is modulated by multiple other immune cell types. In SLE, essentially every immune cell type is dysregulated and aberrant deregulation is thought to be due, in part, to direct or indirect exposure to IFN. Genetic variants within or around the transcription factor interferon regulatory factor 5 (IRF5) associate with SLE risk. Elevated IFNα activity was detected in the sera of SLE patients carrying IRF5 risk polymorphisms who were positive for either anti-RNA binding protein (anti-RBP) or anti-double-stranded DNA (anti-dsDNA) autoantibodies. Neutrophils are also an important source of type I IFNs and are found in abundance in human blood. Neutrophil extracellular traps (NETs) are considered a potential source of antigenic trigger in SLE that can lead to type I IFN gene induction, as well as increased autoantibody production. In this review, we will focus on immune cell types that produce type I IFNs and/or are affected by type I IFN in SLE. In addition, we will discuss potential inducers of endogenous type I IFN production in SLE. Last, we will postulate how the different immune cell populations may be affected by an IRF5-SLE risk haplotype.

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.

The IRF5-TNPO3 association with systemic lupus erythematosus has two components that other autoimmune disorders variably share

Exploiting genotyping, DNA sequencing, imputation and trans-ancestral mapping, we used Bayesian and frequentist approaches to model the IRF5-TNPO3 locus association, now implicated in two immunotherapies and seven autoimmune diseases. Specifically, in systemic lupus erythematosus (SLE), we resolved separate associations in the IRF5 promoter (all ancestries) and with an extended European haplotype. We captured 3230 IRF5-TNPO3 high-quality, common variants across 5 ethnicities in 8395 SLE cases and 7367 controls. The genetic effect from the IRF5 promoter can be explained by any one of four variants in 5.7 kb (P-valuemeta = 6 × 10(-49); OR = 1.38-1.97). The second genetic effect spanned an 85.5-kb, 24-variant haplotype that included the genes IRF5 and TNPO3 (P-valuesEU = 10(-27)-10(-32), OR = 1.7-1.81). Many variants at the IRF5 locus with previously assigned biological function are not members of either final credible set of potential causal variants identified herein. In addition to the known biologically functional variants, we demonstrated that the risk allele of rs4728142, a variant in the promoter among the lowest frequentist probability and highest Bayesian posterior probability, was correlated with IRF5 expression and differentially binds the transcription factor ZBTB3. Our analytical strategy provides a novel framework for future studies aimed at dissecting etiological genetic effects. Finally, both SLE elements of the statistical model appear to operate in Sjögren's syndrome and systemic sclerosis whereas only the IRF5-TNPO3 gene-spanning haplotype is associated with primary biliary cirrhosis, demonstrating the nuance of similarity and difference in autoimmune disease risk mechanisms at IRF5-TNPO3.

Genome-wide association study combined with biological context can reveal more disease-related SNPs altering microRNA target seed sites

Background

Emerging studies demonstrate that single nucleotide polymorphisms (SNPs) resided in the microRNA recognition element seed sites (MRESSs) in 3′UTR of mRNAs are putative biomarkers for human diseases and cancers. However, exhaustively experimental validation for the causality of MRESS SNPs is impractical. Therefore bioinformatics have been introduced to predict causal MRESS SNPs. Genome-wide association study (GWAS) provides a way to detect susceptibility of millions of SNPs simultaneously by taking linkage disequilibrium (LD) into account, but the multiple-testing corrections implemented to suppress false positive rate always sacrificed the sensitivity. In our study, we proposed a method to identify candidate causal MRESS SNPs from 12 GWAS datasets without performing multiple-testing corrections. Alternatively, we used biological context to ensure credibility of the selected SNPs.

Results

In 11 out of the 12 GWAS datasets, MRESS SNPs were over-represented in SNPs with p-value ≤ 0.05 (odds ratio (OR) ranged from 1.1 to 2.4). Moreover, host genes of susceptible MRESS SNPs in each of the 11 GWAS dataset shared biological context with reported causal genes. There were 286 MRESS SNPs identified by our method, while only 13 SNPs were identified by multiple-testing corrections with a given threshold of 1 × 10⁻⁵, which is a common cutoff used in GWAS. 27 out of the 286 candidate SNPs have been reported to be deleterious while only 2 out of 13 multiple-testing corrected SNPs were documented in PubMed. MicroRNA-mRNA interactions affected by the 286 candidate SNPs were likely to present negatively correlated expression. These SNPs introduced greater alternation of binding free energy than other MRESS SNPs, especially when grouping by haplotypes (4210 vs. 4105 cal/mol by mean, 9781 vs. 8521 cal/mol by mean, respectively).

Conclusions

MRESS SNPs are promising disease biomarkers in multiple GWAS datasets. The method of integrating GWAS p-value and biological context is stable and effective for selecting candidate causal MRESS SNPs, it reduces the loss of sensitivity compared to multiple-testing corrections. The 286 candidate causal MRESS SNPs provide researchers a credible source to initialize their design of experimental validations in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-669) contains supplementary material, which is available to authorized users.

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

The rs2004640 single nucleotide polymorphism and the CGGGG copy-number variant (rs77571059) are promoter polymorphisms within interferon regulatory factor 5 (IRF5). They have been implicated as susceptibility factors for several autoimmune diseases. IRF5 uses alternative promoter splicing, where any of 4 first exons begin the mRNA. The CGGGG indel is in exon 1A's promoter; the rs2004640 allele creates a splicing recognition site, enabling usage of exon 1B. This study aimed at characterizing alterations in IRF5 mRNA due to these polymorphisms. Cells with risk polymorphisms exhibited ~2-fold higher levels of IRF5 mRNA and protein, but demonstrated no change in mRNA stability. Quantitative PCR demonstrated decreased usage of exons 1C and 1D in cell lines with the risk polymorphisms. RNA folding analysis revealed a hairpin in exon 1B; mutational analysis showed that the hairpin shape decreased translation 5-fold. Although translation of mRNA that uses exon 1B is low due to a hairpin, increased IRF5 mRNA levels in individuals with the rs2004640 risk allele lead to higher overall protein expression. In addition, several new splice variants of IRF5 were sequenced. IRF5's promoter polymorphisms alter first exon usage and increase transcription levels. High levels of IRF5 may bias the immune system toward autoimmunity.

Four Promoters of IRF5 Respond Distinctly to Stimuli and are Affected by Autoimmune-Risk Polymorphisms

Introduction: Autoimmune diseases such as systemic lupus erythematosus, rheumatoid arthritis, and multiple sclerosis affect millions of people worldwide. Interferon regulatory factor 5 (IRF5) contains polymorphisms associated with these autoimmune diseases. Two of these functional polymorphisms are found upstream of the IRF5 gene. rs2004640, which is a single nucleotide polymorphism and the CGGGG insertion/deletion (indel) were studied. IRF5 uses four different promoters for its four first exons: 1A, 1B, 1C, and 1D. Each promoter was analyzed, including functional differences due to the autoimmune-risk polymorphisms.

Results: IRF5 promoters were analyzed using ChIP-Seq data (ENCODE database) and the FactorBook database to define transcription factor binding sites. To verify promoter activity, the promoters were cloned into luciferase plasmids. Each construct exhibited luciferase activity. Exons 1A and 1D contain putative PU.1 and NFkB binding sites. Imiquimod, a Toll-like receptor 7 (TLR7) ligand, was used to activate these transcription factors. IRF5 levels were doubled after imiquimod treatment (p < 0.001), with specific increases in the 1A promoter (2.2-fold, p = 0.03) and 1D promoter (2.8-fold, p = 0.03). A putative binding site for p53, which affects apoptosis, was found in the promoter for exon 1B. However, site-directed mutagenesis of the p53 site showed no effect in a reporter assay.

Conclusion: The IRF5 exon 1B promoter has been characterized, and the responses of each IRF5 promoter to TLR7 stimulation have been determined. Changes in promoter activity and gene expression are likely due to specific and distinct transcription factors that bind to each promoter. Since high expression of IRF5 contributes to the development of autoimmune disease, understanding the source of increased IRF5 levels is key to understanding autoimmune etiology.

Identification of three new cis-regulatory IRF5 polymorphisms: in vitro studies

Background

Polymorphisms in the interferon regulatory factor 5 (IRF5) gene are associated with susceptibility to systemic lupus erythematosus, rheumatoid arthritis and other diseases through independent risk and protective haplotypes. Several functional polymorphisms are already known, but they do not account for the protective haplotypes that are tagged by the minor allele of rs729302.

Methods

Polymorphisms in linkage disequilibrium (LD) with rs729302 or particularly associated with IRF5 expression were selected for functional screening, which involved electrophoretic mobility shift assays (EMSAs) and reporter gene assays.

Results

A total of 54 single-nucleotide polymorphisms in the 5' region of IRF5 were genotyped. Twenty-four of them were selected for functional screening because of their high LD with rs729302 or protective haplotypes. In addition, two polymorphisms were selected for their prominent association with IRF5 expression. Seven of these twenty-six polymorphisms showed reproducible allele differences in EMSA. The seven were subsequently analyzed in gene reporter assays, and three of them showed significant differences between their two alleles: rs729302, rs13245639 and rs11269962. Haplotypes including the cis-regulatory polymorphisms correlated very well with IRF5 mRNA expression in an analysis based on previous data.

Conclusion

We have found that three polymorphisms in LD with the protective haplotypes of IRF5 have differential allele effects in EMSA and in reporter gene assays. Identification of these cis-regulatory polymorphisms will allow more accurate analysis of transcriptional regulation of IRF5 expression, more powerful genetic association studies and deeper insight into the role of IRF5 in disease susceptibility.

RNA-Seq for enrichment and analysis of IRF5 transcript expression in SLE

Polymorphisms in the interferon regulatory factor 5 (IRF5) gene have been consistently replicated and shown to confer risk for or protection from the development of systemic lupus erythematosus (SLE). IRF5 expression is significantly upregulated in SLE patients and upregulation associates with IRF5-SLE risk haplotypes. IRF5 alternative splicing has also been shown to be elevated in SLE patients. Given that human IRF5 exists as multiple alternatively spliced transcripts with distinct function(s), it is important to determine whether the IRF5 transcript profile expressed in healthy donor immune cells is different from that expressed in SLE patients. Moreover, it is not currently known whether an IRF5-SLE risk haplotype defines the profile of IRF5 transcripts expressed. Using standard molecular cloning techniques, we identified and isolated 14 new differentially spliced IRF5 transcript variants from purified monocytes of healthy donors and SLE patients to generate an IRF5 variant transcriptome. Next-generation sequencing was then used to perform in-depth and quantitative analysis of full-length IRF5 transcript expression in primary immune cells of SLE patients and healthy donors by next-generation sequencing. Evidence for additional alternatively spliced transcripts was obtained from de novo junction discovery. Data from these studies support the overall complexity of IRF5 alternative splicing in SLE. Results from next-generation sequencing correlated with cloning and gave similar abundance rankings in SLE patients thus supporting the use of this new technology for in-depth single gene transcript profiling. Results from this study provide the first proof that 1) SLE patients express an IRF5 transcript signature that is distinct from healthy donors, 2) an IRF5-SLE risk haplotype defines the top four most abundant IRF5 transcripts expressed in SLE patients, and 3) an IRF5 transcript signature enables clustering of SLE patients with the H2 risk haplotype.

An integrative functional genomics approach for discovering biomarkers in schizophrenia

Schizophrenia (SZ) is a complex disorder resulting from both genetic and environmental causes with a lifetime prevalence world-wide of 1%; however, there are no specific, sensitive and validated biomarkers for SZ. A general unifying hypothesis has been put forward that disease-associated single nucleotide polymorphisms (SNPs) from genome-wide association study (GWAS) are more likely to be associated with gene expression quantitative trait loci (eQTL). We will describe this hypothesis and review primary methodology with refinements for testing this paradigmatic approach in SZ. We will describe biomarker studies of SZ and testing enrichment of SNPs that are associated both with eQTLs and existing GWAS of SZ. SZ-associated SNPs that overlap with eQTLs can be placed into gene-gene expression, protein-protein and protein-DNA interaction networks. Further, those networks can be tested by reducing/silencing the gene expression levels of critical nodes. We present pilot data to support these methods of investigation such as the use of eQTLs to annotate GWASs of SZ, which could be applied to the field of biomarker discovery. Those networks that have association with SNP markers, especially cis-regulated expression, might lead to a more clear understanding of important candidate genes that predispose to disease and alter expression. This method has general application to many complex disorders.