Risk alleles for systemic lupus erythematosus in a large case-control collection and associations with clinical subphenotypes.

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. Recent studies have greatly expanded the number of established SLE risk alleles, but the distribution of multiple risk alleles in cases versus controls and their relationship to subphenotypes have not been studied. We studied 22 SLE susceptibility polymorphisms with previous genome-wide evidence of association (p < 5 x 10⁻¹²8) in 1919 SLE cases from 9 independent Caucasian SLE case series and 4813 independent controls. The mean number of risk alleles in cases was 15.1 (SD 3.1) while the mean in controls was 13.1 (SD 2.8), with trend p = 4 x 10⁻8. We defined a genetic risk score (GRS) for SLE as the number of risk alleles with each weighted by the SLE risk odds ratio (OR). The OR for high-low GRS tertiles, adjusted for intra-European ancestry, sex, and parent study, was 4.4 (95% CI 3.8-5.1). We studied associations of individual SNPs and the GRS with clinical manifestations for the cases: age at diagnosis, the 11 American College of Rheumatology classification criteria, and double-stranded DNA antibody (anti-dsDNA) production. Six subphenotypes were significantly associated with the GRS, most notably anti-dsDNA (OR(high-low) = 2.36, p = 9e-9), the immunologic criterion (OR(high-low) = 2.23, p = 3e-7), and age at diagnosis (OR(high-low) = 1.45, p = 0.0060). Finally, we developed a subphenotype-specific GRS (sub-GRS) for each phenotype with more power to detect cumulative genetic associations. The sub-GRS was more strongly associated than any single SNP effect for 5 subphenotypes (the above plus hematologic disorder and oral ulcers), while single loci are more significantly associated with renal disease (HLA-DRB1, OR = 1.37, 95% CI 1.14-1.64) and arthritis (ITGAM, OR = 0.72, 95% CI 0.59-0.88). We did not observe significant associations for other subphenotypes, for individual loci or the sub-GRS. Thus our analysis categorizes SLE subphenotypes into three groups: those having cumulative, single, and no known genetic association with respect to the currently established SLE risk loci.

Differential genetic associations for systemic lupus erythematosus based on anti-dsDNA autoantibody production.

Systemic lupus erythematosus (SLE) is a clinically heterogeneous, systemic autoimmune disease characterized by autoantibody formation. Previously published genome-wide association studies (GWAS) have investigated SLE as a single phenotype. Therefore, we conducted a GWAS to identify genetic factors associated with anti-dsDNA autoantibody production, a SLE-related autoantibody with diagnostic and clinical importance. Using two independent datasets, over 400,000 single nucleotide polymorphisms (SNPs) were studied in a total of 1,717 SLE cases and 4,813 healthy controls. Anti-dsDNA autoantibody positive (anti-dsDNA +, n = 811) and anti-dsDNA autoantibody negative (anti-dsDNA -, n = 906) SLE cases were compared to healthy controls and to each other to identify SNPs associated specifically with these SLE subtypes. SNPs in the previously identified SLE susceptibility loci STAT4, IRF5, ITGAM, and the major histocompatibility complex were strongly associated with anti-dsDNA + SLE. Far fewer and weaker associations were observed for anti-dsDNA - SLE. For example, rs7574865 in STAT4 had an OR for anti-dsDNA + SLE of 1.77 (95% CI 1.57-1.99, p = 2.0E-20) compared to an OR for anti-dsDNA - SLE of 1.26 (95% CI 1.12-1.41, p = 2.4E-04), with p(heterogeneity)<0.0005. SNPs in the SLE susceptibility loci BANK1, KIAA1542, and UBE2L3 showed evidence of association with anti-dsDNA + SLE and were not associated with anti-dsDNA - SLE. In conclusion, we identified differential genetic associations with SLE based on anti-dsDNA autoantibody production. Many previously identified SLE susceptibility loci may confer disease risk through their role in autoantibody production and be more accurately described as autoantibody propensity loci. Lack of strong SNP associations may suggest that other types of genetic variation or non-genetic factors such as environmental exposures have a greater impact on susceptibility to anti-dsDNA - SLE.

Clinical and biological data integration for biomarker discovery.

Biomarkers hold promise for increasing success rates of clinical trials. Biomarker discovery requires searching for associations across a spectrum of data. The field of biomedical data integration has made strides in developing management and analysis tools for structured biological data, but best practices are still evolving for the integration of high-throughput data with less structured clinical data. Integrated repositories are needed to support data analysis, storage and access. We describe a data integration strategy that implements a clinical and biological database and a wiki interface. We integrated parameters across clinical trials and associated genetic, gene expression and protein data. We provide examples to illustrate the utility of data integration to explore disease heterogeneity and develop predictive biomarkers.

European population substructure is associated with mucocutaneous manifestations and autoantibody production in systemic lupus erythematosus.

OBJECTIVE: To determine whether genetic substructure in European-derived populations is associated with specific manifestations of systemic lupus erythematosus (SLE), including mucocutaneous phenotypes, autoantibody production, and renal disease. METHODS: SLE patients of European descent (n=1,754) from 8 case collections were genotyped for >1,400 ancestry informative markers that define a north-south gradient of European substructure. Using the Structure program, each SLE patient was characterized in terms of percent Northern (versus percent Southern) European ancestry based on these genetic markers. Nonparametric methods, including tests for trend, were used to identify associations between Northern European ancestry and specific SLE manifestations. RESULTS: In multivariate analyses, increasing levels of Northern European ancestry were significantly associated with photosensitivity (Ptrend=0.0021, odds ratio for highest quartile of Northern European ancestry versus lowest quartile [ORhigh-low] 1.64, 95% confidence interval [95% CI] 1.13-2.35) and discoid rash (Ptrend=0.014, ORhigh-low 1.93, 95% CI 0.98-3.83). In contrast, increasing levels of Northern European ancestry had a protective effect against the production of anticardiolipin autoantibodies (Ptrend=1.6x10(-4), ORhigh-low 0.46, 95% CI 0.30-0.69) and anti-double-stranded DNA autoantibodies (Ptrend=0.017, ORhigh-low 0.67, 95% CI 0.46-0.96). CONCLUSION: This study demonstrates that specific SLE manifestations vary according to Northern versus Southern European ancestry. Thus, genetic ancestry may contribute to the clinical heterogeneity and variation in disease outcomes among SLE patients of European descent. Moreover, these results suggest that genetic studies of SLE subphenotypes will need to carefully address issues of population substructure based on genetic ancestry.

Specificity of the STAT4 genetic association for severe disease manifestations of systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. A polymorphism in the STAT4 gene has recently been established as a risk factor for SLE, but the relationship with specific SLE subphenotypes has not been studied. We studied 137 SNPs in the STAT4 region genotyped in 4 independent SLE case series (total n = 1398) and 2560 healthy controls, along with clinical data for the cases. Using conditional testing, we confirmed the most significant STAT4 haplotype for SLE risk. We then studied a SNP marking this haplotype for association with specific SLE subphenotypes, including autoantibody production, nephritis, arthritis, mucocutaneous manifestations, and age at diagnosis. To prevent possible type-I errors from population stratification, we reanalyzed the data using a subset of subjects determined to be most homogeneous based on principal components analysis of genome-wide data. We confirmed that four SNPs in very high LD (r(2) = 0.94 to 0.99) were most strongly associated with SLE, and there was no compelling evidence for additional SLE risk loci in the STAT4 region. SNP rs7574865 marking this haplotype had a minor allele frequency (MAF) = 31.1% in SLE cases compared with 22.5% in controls (OR = 1.56, p = 10(-16)). This SNP was more strongly associated with SLE characterized by double-stranded DNA autoantibodies (MAF = 35.1%, OR = 1.86, p<10(-19)), nephritis (MAF = 34.3%, OR = 1.80, p<10(-11)), and age at diagnosis<30 years (MAF = 33.8%, OR = 1.77, p<10(-13)). An association with severe nephritis was even more striking (MAF = 39.2%, OR = 2.35, p<10(-4) in the homogeneous subset of subjects). In contrast, STAT4 was less strongly associated with oral ulcers, a manifestation associated with milder disease. We conclude that this common polymorphism of STAT4 contributes to the phenotypic heterogeneity of SLE, predisposing specifically to more severe disease.

An interferon signature in the peripheral blood of dermatomyositis patients is associated with disease activity.

Recent studies have shown increased expression of interferon (IFN)-regulated genes in the peripheral blood cells of patients with systemic lupus erythematosus. A similar interferon signature has been observed in affected muscle tissue from patients with dermatomyositis (DM), but it has not yet been determined if this signature extends to the peripheral blood in DM. We performed global gene expression profiling of peripheral blood cells from adult and juvenile DM patients and healthy controls. Several interesting groups of genes were differentially expressed in DM, including genes with immune function, and others that function in muscle or are involved in mitochondrial/oxidative phosphorylation. Investigation of type I IFN-regulated transcripts revealed a striking interferon signature present in most DM patients studied. Levels of type I IFN-regulated proteins were also elevated in DM serum samples. Furthermore, both the transcript and serum protein IFN signatures were associated with disease activity. These data suggest that the IFN signature may be a useful marker for DM disease activity, and that sampling peripheral blood may be a more practical alternative to muscle biopsy for measuring this signature.

Three functional variants of IFN regulatory factor 5 (IRF5) define risk and protective haplotypes for human lupus.

Systematic genome-wide studies to map genomic regions associated with human diseases are becoming more practical. Increasingly, efforts will be focused on the identification of the specific functional variants responsible for the disease. The challenges of identifying causal variants include the need for complete ascertainment of genetic variants and the need to consider the possibility of multiple causal alleles. We recently reported that risk of systemic lupus erythematosus (SLE) is strongly associated with a common SNP in IFN regulatory factor 5 (IRF5), and that this variant altered spicing in a way that might provide a functional explanation for the reproducible association to SLE risk. Here, by resequencing and genotyping in patients with SLE, we find evidence for three functional alleles of IRF5: the previously described exon 1B splice site variant, a 30-bp in-frame insertion/deletion variant of exon 6 that alters a proline-, glutamic acid-, serine- and threonine-rich domain region, and a variant in a conserved polyA+ signal sequence that alters the length of the 3' UTR and stability of IRF5 mRNAs. Haplotypes of these three variants define at least three distinct levels of risk to SLE. Understanding how combinations of variants influence IRF5 function may offer etiological and therapeutic insights in SLE; more generally, IRF5 and SLE illustrates how multiple common variants of the same gene can together influence risk of common disease.

Elevated serum levels of interferon-regulated chemokines are biomarkers for active human systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) is a serious systemic autoimmune disorder that affects multiple organ systems and is characterized by unpredictable flares of disease. Recent evidence indicates a role for type I interferon (IFN) in SLE pathogenesis; however, the downstream effects of IFN pathway activation are not well understood. Here we test the hypothesis that type I IFN-regulated proteins are present in the serum of SLE patients and correlate with disease activity. METHODS AND FINDINGS: We performed a comprehensive survey of the serologic proteome in human SLE and identified dysregulated levels of 30 cytokines, chemokines, growth factors, and soluble receptors. Particularly striking was the highly coordinated up-regulation of 12 inflammatory and/or homeostatic chemokines, molecules that direct the movement of leukocytes in the body. Most of the identified chemokines were inducible by type I IFN, and their levels correlated strongly with clinical and laboratory measures of disease activity. CONCLUSIONS: These data suggest that severely disrupted chemokine gradients may contribute to the systemic autoimmunity observed in human SLE. Furthermore, the levels of serum chemokines may serve as convenient biomarkers for disease activity in lupus.

A common haplotype of interferon regulatory factor 5 (IRF5) regulates splicing and expression and is associated with increased risk of systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease characterized by activation of the type I interferon (IFN) pathway. Here we convincingly replicate association of the IFN regulatory factor 5 (IRF5) rs2004640 T allele with SLE in four independent case-control cohorts (P = 4.4 x 10(-16)) and by family-based transmission disequilibrium test analysis (P = 0.0006). The rs2004640 T allele creates a 5' donor splice site in an alternate exon 1 of IRF5, allowing expression of several unique IRF5 isoforms. We also identify an independent cis-acting variant associated with elevated expression of IRF5 and linked to the exon 1B splice site. Haplotypes carrying the variant associated with elevated expression and lacking the exon 1B donor site do not confer risk of SLE. Thus, a common IRF5 haplotype driving elevated expression of multiple unique isoforms of IRF5 is an important genetic risk factor for SLE, establishing a causal role for type I IFN pathway genes in human autoimmunity.

Fine-mapping chromosome 20 in 230 systemic lupus erythematosus sib pair and multiplex families: evidence for genetic epistasis with chromosome 16q12.

The presence of systemic lupus erythematosus (SLE) susceptibility genes on chromosome 20 is suggested by the observation of genetic linkage in several independent SLE family collections. To further localize the genetic effects, we typed 59 microsatellites in the two best regions, as defined by genome screens. Genotypes were analyzed for statistical linkage and/or association with SLE, by use of a combination of nonparametric linkage methods, family-based tests of association (transmission/disequilibrium and pedigree disequilibrium tests), and haplotype-sharing statistics (haplotype runs test), in a set of 230 SLE pedigrees. Maximal evidence for linkage to SLE was to 20p12 (LOD = 2.84) and 20q13.1 (LOD = 1.64) in the white pedigrees. Subsetting families on the basis of evidence for linkage to 16q12 significantly improved the LOD scores at both chromosome 20 locations (20p12 LOD = 5.06 and 20q13 LOD = 3.65), consistent with epistasis. We then typed 162 single-nucleotide polymorphism markers across a 1.3-Mb candidate region on 20q13.1 and identified several SNPs that demonstrated significant evidence for association. These data provide additional support for linkage and association to 20p12 and 20q13.1 in SLE and further refine the intervals of interest. These data further suggest the possibility of epistatic relationships among loci within the 20q12, 20q13, and 16q12 regions in SLE families.

Progress towards understanding the genetic pathogenesis of systemic lupus erythematosus.

In order to better understand the genetic factors that initiate systemic lupus erythematosus (SLE), we are using both linkage and association approaches to identify susceptibility genes for the disease. Association studies have recently identified three HLA Class II haplotypes as well as a functional missense polymorphism in protein tyrosine phosphatase (PTP) PTPN22 as important risk alleles for SLE. Here, we will review these data, and explain how these findings contribute to an understanding of the genetic architecture of human SLE.

Genetic association of the R620W polymorphism of protein tyrosine phosphatase PTPN22 with human SLE.

We genotyped 525 independent North American white individuals with systemic lupus erythematosus (SLE) for the PTPN22 R620W polymorphism and compared the results with data generated from 1,961 white control individuals. The R620W SNP was associated with SLE (genotypic P=.00009), with estimated minor (T) allele frequencies of 12.67% in SLE cases and 8.64% in controls. A single copy of the T allele (W620) increases risk of SLE (odds ratio [OR]=1.37; 95% confidence interval [CI] 1.07-1.75), and two copies of the allele more than double this risk (OR=4.37; 95% CI 1.98-9.65). Together with recent evidence showing association of this SNP with type 1 diabetes and rheumatoid arthritis, these data provide compelling evidence that PTPN22 plays a fundamental role in regulating the immune system and the development of autoimmunity.

Interferon-inducible gene expression signature in peripheral blood cells of patients with severe lupus.

Systemic lupus erythematosus (SLE) is a complex, inflammatory autoimmune disease that affects multiple organ systems. We used global gene expression profiling of peripheral blood mononuclear cells to identify distinct patterns of gene expression that distinguish most SLE patients from healthy controls. Strikingly, about half of the patients studied showed dysregulated expression of genes in the IFN pathway. Furthermore, this IFN gene expression "signature" served as a marker for more severe disease involving the kidneys, hematopoetic cells, and/or the central nervous system. These results provide insights into the genetic pathways underlying SLE, and identify a subgroup of patients who may benefit from therapies targeting the IFN pathway.

Visualizing human leukocyte antigen class II risk haplotypes in human systemic lupus erythematosus.

Human leukocyte antigen (HLA) class I and class II alleles are implicated as genetic risk factors for many autoimmune diseases. However, the role of the HLA loci in human systemic lupus erythematosus (SLE) remains unclear. Using a dense map of polymorphic microsatellites across the HLA region in a large collection of families with SLE, we identified three distinct haplotypes that encompassed the class II region and exhibited transmission distortion. DRB1 and DQB1 typing of founders showed that the three haplotypes contained DRB1*1501/ DQB1*0602, DRB1*0801/ DQB1*0402, and DRB1*0301/DQB1*0201 alleles, respectively. By visualizing ancestral recombinants, we narrowed the disease-associated haplotypes containing DRB1*1501 and DRB1*0801 to an approximately 500-kb region. We conclude that HLA class II haplotypes containing DRB1 and DQB1 alleles are strong risk factors for human SLE.