A comprehensive analysis of shared loci between systemic lupus erythematosus (SLE) and sixteen autoimmune diseases reveals limited genetic overlap

In spite of the well-known clustering of multiple autoimmune disorders in families, analyses of specific shared genes and polymorphisms between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) have been limited. Therefore, we comprehensively tested autoimmune variants for association with SLE, aiming to identify pleiotropic genetic associations between these diseases. We compiled a list of 446 non–Major Histocompatibility Complex (MHC) variants identified in genome-wide association studies (GWAS) of populations of European ancestry across 17 ADs. We then tested these variants in our combined Caucasian SLE cohorts of 1,500 cases and 5,706 controls. We tested a subset of these polymorphisms in an independent Caucasian replication cohort of 2,085 SLE cases and 2,854 controls, allowing the computation of a meta-analysis between all cohorts. We have uncovered novel shared SLE loci that passed multiple comparisons adjustment, including the VTCN1 (rs12046117, P = 2.02×10⁻⁰⁶) region. We observed that the loci shared among the most ADs include IL23R, OLIG3/TNFAIP3, and IL2RA. Given the lack of a universal autoimmune risk locus outside of the MHC and variable specificities for different diseases, our data suggests partial pleiotropy among ADs. Hierarchical clustering of ADs suggested that the most genetically related ADs appear to be type 1 diabetes with rheumatoid arthritis and Crohn's disease with ulcerative colitis. These findings support a relatively distinct genetic susceptibility for SLE. For many of the shared GWAS autoimmune loci, we found no evidence for association with SLE, including IL23R. Also, several established SLE loci are apparently not associated with other ADs, including the ITGAM-ITGAX and TNFSF4 regions. This study represents the most comprehensive evaluation of shared autoimmune loci to date, supports a relatively distinct non–MHC genetic susceptibility for SLE, provides further evidence for previously and newly identified shared genes in SLE, and highlights the value of studies of potentially pleiotropic genes in autoimmune diseases.

It is well known that multiple autoimmune disorders cluster in families. However, all of the genetic variants that explain this clustering have not been discovered, and the specific genetic variants shared between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) are not known. In order to better understand the genetic factors that explain this predisposition to autoimmunity, we performed a comprehensive evaluation of shared autoimmune genetic variants. First we considered results from 17 ADs and compiled a list with 446 significant genetic variants from these studies. We identified some genetic variants extensively shared between ADs, as well as the ADs that share the most variants. The genetic overlap between SLE and other ADs was modest. Next we tested how important all the 446 genetic variants were in our collection with a minimum of 1,500 SLE patients. Among the most significant variants in SLE, the majority had already been identified in previous studies, but we also discovered variants in two important immune genes. In summary, our data identified diseases with common genetic risk factors and novel SLE effects, and this supports a relatively distinct genetic susceptibility for SLE. This study helps delineate the genetic architecture of ADs.

Association of PPP2CA polymorphisms with systemic lupus erythematosus susceptibility in multiple ethnic groups

OBJECTIVE

T cells from patients with systemic lupus erythematosus (SLE) express increased amounts of PP2Ac, which contributes to decreased production of interleukin-2 (IL-2). Because IL-2 is important in the regulation of several aspects of the immune response, it has been proposed that PP2Ac contributes to the expression of SLE. This study was designed to determine whether genetic variants of PPP2AC are linked to the expression of SLE and specific clinical manifestations and account for the increased expression of PP2Ac.

METHODS

We conducted a trans-ethnic study of 8,695 SLE cases and 7,308 controls of 4 different ancestries. Eighteen single-nucleotide polymorphisms (SNPs) across PPP2CA were genotyped using an Illumina custom array. PPP2CA expression in SLE and control T cells was analyzed by real-time polymerase chain reaction.

RESULTS

A 32-kb haplotype comprising multiple SNPs of PPP2CA showed significant association with SLE in Hispanic Americans, European Americans, and Asians, but not in African Americans. Conditional analyses revealed that SNP rs7704116 in intron 1 showed consistently strong association with SLE across Asian, European American, and Hispanic American populations (odds ratio 1.3 [95% confidence interval 1.14-1.31], meta-analysis P=3.8×10(-7)). In European Americans, the largest ethnic data set studied, the risk A allele of rs7704116 was associated with the presence of renal disease, anti-double-stranded DNA, and anti-RNP antibodies. PPP2CA expression was ∼2-fold higher in SLE patients carrying the rs7704116 AG genotype than those carrying the GG genotype (P=0.007).

CONCLUSION

Our data provide the first evidence of an association between PPP2CA polymorphisms and elevated PP2Ac transcript levels in T cells, which implicates a new molecular pathway for SLE susceptibility in European Americans, Hispanic Americans, and Asians.

Polymorphisms in MYH9 are associated with diabetic nephropathy in European Americans

BACKGROUND

Polymorphisms in the non-muscle myosin IIA gene (MYH9) are associated with focal segmental glomerulosclerosis (FSGS) and non-diabetic end-stage renal disease (ESRD) in African Americans and FSGS in European Americans. We tested for association of single nucleotide polymorphisms (SNPs) in MYH9 with T2DM-ESRD in European Americans; additionally, three APOL1 gene variants were evaluated.

METHODS

Fifteen MYH9 SNPs and two APOL1 SNPs plus a 6-bp deletion were genotyped in 1963 European Americans, 536 cases with T2DM-ESRD and 1427 non-nephropathy controls (467 with T2DM and 960 without diabetes).

RESULTS

Comparing T2DM-ESRD cases with the 467 T2DM non-nephropathy controls, single variant associations trending toward significance were detected with SNPs rs4821480, rs2032487 and rs4281481 comprising part of the major MYH9 E1 risk haplotype [P-values 0.053-0.055 recessive, odds ratio (OR) 6.08-6.14]. Comparing T2DM-ESRD cases to all 1427 non-nephropathy controls, we confirmed evidence of association in these three SNPs as well as in the fourth E1 SNP (rs3752462) (P-values 0.017-0.035, OR 1.41-3.72). APOL1 G1/G2 nephropathy risk variants were rare in individuals of European American heritage, present in 0.28% of chromosomes in T2DM-ESRD cases and 0.32% of controls.

CONCLUSIONS

MYH9 SNPs rs4821480, rs2032487, rs4281481 and rs3752462 are associated with T2DM-ESRD susceptibility in European Americans. The APOL1 risk variants are not present at appreciable frequency in this cohort with T2DM-ESRD. Therefore, polymorphisms in MYH9 appear to influence nephropathy risk in this sample.

Genetic analyses of interferon pathway-related genes reveal multiple new loci associated with systemic lupus erythematosus

OBJECTIVE

The overexpression of interferon (IFN)-inducible genes is a prominent feature of systemic lupus erythematosus (SLE); it serves as a marker for active and more severe disease, and is also observed in other autoimmune and inflammatory conditions. This study was undertaken to investigate the genetic variations responsible for sustained activation of IFN-responsive genes in SLE.

METHODS

We systematically evaluated association of SLE with a total of 1,754 IFN pathway-related genes, including IFN-inducible genes known to be differentially expressed in SLE patients and their direct regulators. We used a 3-stage study design in which 2 cohorts (total of 939 SLE cases and 3,398 controls) were analyzed independently and jointly for association with SLE, and the results were adjusted for the number of comparisons.

RESULTS

A total of 15,166 single-nucleotide polymorphisms (SNPs) passed all quality control filters; 305 of these SNPs demonstrated replicated association with SLE in both cohorts. Nine variants were further genotyped for confirmation in an average of 1,316 independent SLE cases and 3,215 independent controls. Association with SLE was confirmed for several genes, including those for the transmembrane receptor CD44 (CD44 [rs507230]; P = 3.98 × 10⁻¹²), the cytokine pleiotrophin (PTN [rs919581]; P = 5.38 × 10⁻⁴), the heat-shock protein DnaJ (DNAJA1 [rs10971259]; P = 6.31 × 10⁻³), and the nuclear import protein karyopherin α1 (KPNA [rs6810306]; P = 4.91 × 10⁻²).

CONCLUSION

This study expands the number of candidate genes that have been shown to be associated with SLE and highlights potential of pathway-based approaches for gene discovery. Identification of the causal alleles will help elucidate the molecular mechanisms responsible for activation of the IFN system in SLE.

Association of genetic variants in complement factor H and factor H-related genes with systemic lupus erythematosus susceptibility

Systemic lupus erythematosus (SLE), a complex polygenic autoimmune disease, is associated with increased complement activation. Variants of genes encoding complement regulator factor H (CFH) and five CFH-related proteins (CFHR1-CFHR5) within the chromosome 1q32 locus linked to SLE, have been associated with multiple human diseases and may contribute to dysregulated complement activation predisposing to SLE. We assessed 60 SNPs covering the CFH-CFHRs region for association with SLE in 15,864 case-control subjects derived from four ethnic groups. Significant allelic associations with SLE were detected in European Americans (EA) and African Americans (AA), which could be attributed to an intronic CFH SNP (rs6677604, in intron 11, P_meta = 6.6×10⁻⁸, OR = 1.18) and an intergenic SNP between CFHR1 and CFHR4 (rs16840639, P_meta = 2.9×10⁻⁷, OR = 1.17) rather than to previously identified disease-associated CFH exonic SNPs, including I62V, Y402H, A474A, and D936E. In addition, allelic association of rs6677604 with SLE was subsequently confirmed in Asians (AS). Haplotype analysis revealed that the underlying causal variant, tagged by rs6677604 and rs16840639, was localized to a ∼146 kb block extending from intron 9 of CFH to downstream of CFHR1. Within this block, the deletion of CFHR3 and CFHR1 (CFHR3-1Δ), a likely causal variant measured using multiplex ligation-dependent probe amplification, was tagged by rs6677604 in EA and AS and rs16840639 in AA, respectively. Deduced from genotypic associations of tag SNPs in EA, AA, and AS, homozygous deletion of CFHR3-1Δ (P_meta = 3.2×10⁻⁷, OR = 1.47) conferred a higher risk of SLE than heterozygous deletion (P_meta = 3.5×10⁻⁴, OR = 1.14). These results suggested that the CFHR3-1Δ deletion within the SLE-associated block, but not the previously described exonic SNPs of CFH, might contribute to the development of SLE in EA, AA, and AS, providing new insights into the role of complement regulators in the pathogenesis of SLE.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease, associated with increased complement activation. Previous studies have provided evidence for the presence of SLE susceptibility gene(s) in the chromosome 1q31-32 locus. Within 1q32, genes encoding complement regulator factor H (CFH) and five CFH-related proteins (CFHR1-CFHR5) may contribute to the development of SLE, because genetic variants of these genes impair complement regulation and predispose to various human diseases. In this study, we tested association of genetic variants in the region containing CFH and CFHRs with SLE. We identified genetic variants predisposing to SLE in European American, African American, and Asian populations, which might be attributed to the deletion of CFHR3 and CFHR1 genes but not previously identified disease-associated exonic variants of CFH. This study provides the first evidence for consistent association between CFH/CFHRs and SLE across multi-ancestral SLE datasets, providing new insights into the role of complement regulators in the pathogenesis of SLE.

The susceptibility loci juvenile idiopathic arthritis shares with other autoimmune diseases extend to PTPN2, COG6, and ANGPT1

OBJECTIVE

To test for associations between non-major histocompatibility complex susceptibility loci previously reported in autoimmune diseases and juvenile idiopathic arthritis (JIA).

METHODS

Published autoimmune disease genome-wide association studies were reviewed, and 519 single-nucleotide polymorphisms (SNPs) were selected for association testing. The initial cohort included 809 JIA cases and 3,535 controls of non-Hispanic, European ancestry. Of the SNPs, 257 were successfully genotyped, while 168 were imputed with quality. Based on findings in the initial cohort, replication was sought for 21 SNPs in a second cohort of 1,015 JIA cases and 1,569 controls collected in the US and Germany. For the initial cohort, tests for association were adjusted for potential confounding effects of population structure by including principal components derived from a genome-wide association study as covariates in logistic regression models. Odds ratios (ORs) and 95% confidence intervals were calculated.

RESULTS

Testing for association of previously reported autoimmune disease genetic associations in the initial cohort suggested associations with JIA in 13 distinct loci. Of these, 7 were validated in the replication cohort. Meta-analysis results for the replicating loci included PTPN22 (rs6679677 [OR 1.58, P = 1.98 × 10(-12) ], rs2476601 [OR 1.64, P = 1.90 × 10(-13) ], and rs2488457 [OR 1.32, P = 6.74 × 10(-8) ]), PTPN2 (rs1893217 [OR = 1.33, P = 1.60 × 10(-9) ] and rs7234029 [OR 1.35, P = 1.86 × 10(-10) ]), ADAD1-IL2-IL21 (rs17388568 [OR 1.24, P = 1.13 × 10(-6) ] and rs13143866 [OR 0.83, P = 1.95 × 10(-4) ]), STAT4 (rs3821236 [OR = 1.27, P = 2.36 × 10(-6) ] and rs7574865 [OR = 1.31, P = 2.21 × 10(-6) ]), C12orf30 (rs17696736 [OR = 1.19, P = 2.59 × 10(-5) ]), COG6 (rs7993214 [OR = 0.76, P = 1.10 × 10(-5) ]), and ANGPT1 (rs1010824 [OR = 0.79, P = 2.91 × 10(-4) ]). These polymorphisms have been reported in diseases such as rheumatoid arthritis, type 1 diabetes mellitus, Crohn's disease, and multiple sclerosis.

CONCLUSION

General susceptibility loci for autoimmunity are shared across diseases, including JIA, suggesting the potential for common therapeutic targets and mechanisms.

The non-muscle Myosin heavy chain 9 gene (MYH9) is not associated with lupus nephritis in African Americans

BACKGROUND

African Americans (AA) disproportionately develop lupus nephritis (LN) relative to European Americans and familial clustering supports causative genes. Since MYH9 underlies approximately 40% of end-stage renal disease (ESRD) in AA, we tested for genetic association with LN.

METHODS

Seven MYH9 single nucleotide polymorphisms (SNPs) and the E1 risk haplotype were tested for association with LN in three cohorts of AA.

RESULTS

A preliminary analysis revealed that the MYH9 E1 risk haplotype was associated with ESRD in 25 cases with presumed systemic lupus erythematosus (SLE)-associated ESRD, compared to 735 non-SLE controls (odds ratio 3.1; p = 0.010 recessive). Replication analyses were performed in 583 AA with SLE in the PROFILE cohort (318 with LN; 265 with SLE but without nephropathy) and 60 AA from the NIH (39 with LN; 21 with SLE but without nephropathy). Analysis of the NIH and larger PROFILE cohorts, as well as a combined analysis, did not support this association.

CONCLUSIONS

These results suggest that AA with ESRD and coincident SLE who were recruited from dialysis clinics more likely have kidney diseases in the MYH9-associated spectrum of focal segmental glomerulosclerosis. PROFILE and NIH participants, recruited from rheumatology practices, demonstrate that MYH9 does not contribute substantially to the development of LN in AA.

The acetyl-coenzyme A carboxylase beta (ACACB) gene is associated with nephropathy in Chinese patients with type 2 diabetes

BACKGROUND

A single-nucleotide polymorphism (SNP), rs2268388, in the acetyl-coenzyme A carboxylase beta (ACACB) gene is associated with susceptibility to type 2 diabetic nephropathy (T2DN) in Japanese and European-American populations. Whether this association also exists in Chinese patients is unclear. Attempts at replication in small Singaporean and Korean samples were not significant.

METHODS

Eight ACACB SNPs were genotyped in 595 subjects with type 2 diabetes mellitus born in Hong Kong or southern China, 295 with advanced T2DN and 300 with long-standing diabetes lacking nephropathy. Association analyses were focused primarily on SNP rs2268388 and secondarily on flanking SNPs and haplotypes.

RESULTS

Adjusting for age, gender and diabetes duration, ACACB SNP rs2268388 was significantly associated with advanced T2DN (odds ratio = 2.39; recessive model; P = 0.0129).

CONCLUSION

These results in the Chinese replicate the association between T2DN and rs2268388, as seen in Japanese and European Americans. The ACACB gene and attendant alterations in fatty acid oxidation may play important roles in susceptibility to T2DN. Targeting this pathway may provide novel treatment options for the prevention of diabetic nephropathy.

Is collapsing C1q nephropathy another MYH9-associated kidney disease? A case report

C1q nephropathy is a rare kidney disease that can present with nephrotic syndrome and typically has the histologic phenotype of either minimal change disease or focal segmental glomerulosclerosis (FSGS). Disagreement exists about whether it is a distinct immune complex-mediated glomerulopathy or it resides in the spectrum of FSGS-minimal change disease. Two African American patients with C1q nephropathy histologically presenting as the collapsing variant of FSGS (collapsing C1q nephropathy) and rapid loss of kidney function were genotyped for polymorphisms in the non-muscle myosin heavy chain 9 gene (MYH9). Both cases were homozygous for the MYH9 E1 risk haplotype, the variant strongly associated with idiopathic FSGS, collapsing FSGS in human immunodeficiency virus-associated nephropathy, and focal global glomerulosclerosis (historically attributed to hypertensive nephrosclerosis). Collapsing C1q nephropathy with rapid progression to end-stage renal disease appears to reside in the MYH9-associated disease spectrum.

Non-muscle myosin heavy chain 9 gene MYH9 associations in African Americans with clinically diagnosed type 2 diabetes mellitus-associated ESRD

BACKGROUND

Although MYH9 is strongly associated with biopsy-proven idiopathic and HIV-associated focal segmental glomerulosclerosis (FSGS) and clinically diagnosed 'hypertension-associated' end-stage renal disease (ESRD) in African Americans, its role in type 2 diabetes mellitus (T2DM)-associated ESRD is unclear.

METHODS

To assess whether MYH9 was associated with T2DM-ESRD, 751 African Americans with T2DM-ESRD, 227 with T2DM lacking nephropathy and 925 non-diabetic non-nephropathy controls were genotyped for 14 MYH9 SNPs. Association analyses used SNPGWA and Dandelion.

RESULTS

Comparing T2DM-ESRD cases with non-diabetic controls, single SNP associations were detected with 8 of 14 SNPs, gender- and admixture-adjusted P-values 0.047-0.005 [recessive model, odds ratio (OR) range 1.30-1.55]. The previously associated MYH9 E1 and L1 haplotypes were associated with T2DM-ESRD (E1: OR 1.27, 95% CI 1.04-1.56, P = 0.021 recessive and L1: OR 1.43, 95% CI 1.09-1.87, P = 0.009 dominant). Contrasting the 751 T2DM-ESRD cases with 227 T2DM non-nephropathy controls revealed that E1 haplotype SNPs rs4821480, rs2032487 and rs4821481 were associated with kidney failure (OR 1.38-1.40 recessive, all P < 0.048). Among E1 and L1 risk homozygotes, respectively, mean (SD) diabetes duration prior to renal replacement therapy was 16.6 (9.7) and 16.4 (10.0) years, and 65% had diabetic retinopathy.

CONCLUSIONS

Genetic dissection of T2DM-associated ESRD reveals that MYH9 underlies a portion of this clinically diagnosed disorder in African Americans. It is likely that a subset of African Americans with T2DM and coincident nephropathy have primary MYH9-related kidney disease (e.g. FSGS or global glomerulosclerosis), although renal biopsy studies need to be performed.