To confirm and define the genetic association of STAT4 and systemic lupus erythematosus (SLE), investigate the possibility of correlations with differential splicing and/or expression levels, and genetic interaction with IRF5.

30 tag SNPs were genotyped in an independent set of Spanish cases and controls. SNPs surviving correction for multiple tests were genotyped in five new sets of cases and controls for replication. STAT4 cDNA was analysed by 5′-RACE PCR and sequencing. Expression levels were measured by quantitative PCR.

In the fine mapping, four SNPs were significant after correction for multiple testing, with rs3821236 and rs3024866 as the strongest signals, followed by the previously associated rs7574865, and by rs1467199. Association was replicated in all cohorts. After conditional regression analyses, two major independent signals, represented by SNPs rs3821236 and rs7574865, remained significant across the sets. These SNPs belong to separate haplotype blocks. High levels of STAT4 expression correlated with SNPs rs3821236, rs3024866 (both in the same haplotype block) and rs7574865 but not with other SNPs. Transcription of alternative tissue-specific exons 1, indicating the presence of tissue-specific promoters of potential importance in the expression of STAT4, was also detected. No interaction with associated SNPs of IRF5 was observed using regression analysis.

These data confirm STAT4 as a susceptibility gene for SLE and suggest the presence of at least two functional variants affecting levels of STAT4. The results also indicate that the genes STAT4 and IRF5 act additively to increase the risk for SLE.