Eight-Fold Increased COVID-19 Mortality in Autosomal Dominant Tubulointerstitial Kidney Disease due to MUC1 Mutations: An Observational Study.

Background: MUC1 and UMOD pathogenic variants cause autosomal dominant tubulointerstitial kidney disease (ADTKD). MUC1 is expressed in kidney, nasal mucosa and respiratory tract, while UMOD is expressed only in kidney. Due to haplo-insufficiency ADTKD- MUC1 patients produce approximately 50% of normal mucin-1. Methods: To determine whether decreased mucin-1 production was associated with an increased COVID-19 risk, we sent a survey to members of an ADTKD registry in September 2021, after the initial, severe wave of COVID-19. We linked results to previously obtained ADTKD genotype and plasma CA15-3 (mucin-1) levels and created a longitudinal registry of COVID-19 related deaths. Results: Surveys were emailed to 637 individuals, with responses from 89 ADTKD- MUC1 and 132 ADTKD- UMOD individuals. 19/83 (23%) ADTKD- MUC1 survey respondents reported a prior COVID-19 infection vs. 14/125 (11%) ADTKD- UMOD respondents (odds ratio (OR) 2.35 (95%CI 1.60-3.11, P = 0.0260). Including additional familial cases reported from survey respondents, 10/41 (24%) ADTKD- MUC1 individuals died of COVID-19 vs. 1/30 (3%) with ADTKD- UMOD , with OR 9.21 (95%CI 1.22-69.32), P = 0.03. The mean plasma mucin-1 level prior to infection in 14 infected and 27 uninfected ADTKD- MUC1 individuals was 7.06±4.12 vs. 10.21±4.02 U/mL ( P = 0.035). Over three years duration, our longitudinal registry identified 19 COVID-19 deaths in 360 ADTKD- MUC1 individuals (5%) vs. 3 deaths in 478 ADTKD- UMOD individuals (0.6%) ( P = 0.0007). Multivariate logistic regression revealed the following odds ratios (95% confidence interval) for COVID-19 deaths: ADTKD- MUC1 8.4 (2.9-29.5), kidney transplant 5.5 (1.6-9.1), body mass index (kg/m 2 ) 1.1 (1.0-1.2), age (y) 1.04 (1.0-1.1). Conclusions: Individuals with ADTKD- MUC1 are at an eight-fold increased risk of COVID-19 mortality vs. ADTKD- UMOD individuals. Haplo-insufficient production of mucin-1 may be responsible.

Maternal health and pregnancy outcomes in autosomal dominant tubulointerstitial kidney disease.

Introduction: Autosomal dominant tubulointerstitial kidney disease (ADTKD) is an increasingly recognized cause of chronic kidney disease. ADTKD pregnancy outcomes have not previously been described. Methods: A cross-sectional survey was sent to women from ADTKD families. Results: Information was obtained from 85 afffected women (164 term pregnancies) and 23 controls (50 pregnancies). Only 16.5% of genetically affected women knew they had ADTKD during pregnancy. Eighteen percent of ADTKD mothers had hypertension during pregnancy versus 12% in controls (p = 0.54) and >40% in comparative studies of chronic kidney disease in pregnancy. Eleven percent of births of ADTKD mothers were <37 weeks versus 0 in controls (p < 0.0001). Cesarean section occurred in 19% of pregnancies in affected women versus 38% of unaffected individuals (p = 0.06). Only 12% of babies required a neonatal intensive care unit stay. Conclusions: ADTKD pregnancies had lower rates of hypertension during pregnancy versus other forms of chronic kidney disease, which may have contributed to good maternal and fetal outcomes.

Genome-Wide Association Study Identifies Loci for Liver Enzyme Concentrations in Mexican Americans: The GUARDIAN Consortium.

OBJECTIVE: Populations of Mexican American ancestry are at an increased risk for nonalcoholic fatty liver disease. The objective of this study was to determine whether loci in known and novel genes were associated with variation in aspartate aminotransferase (AST) (n = 3,644), alanine aminotransferase (ALT) (n = 3,595), and gamma-glutamyl transferase (GGT) (n = 1,577) levels by conducting the first genome-wide association study (GWAS) of liver enzymes, which commonly measure liver function, in individuals of Mexican American ancestry. METHODS: Levels of AST, ALT, and GGT were determined by enzymatic colorimetric assays. A multi-cohort GWAS of individuals of Mexican American ancestry was performed. Single-nucleotide polymorphisms (SNP) were tested for association with liver outcomes by multivariable linear regression using an additive genetic model. Association analyses were conducted separately in each cohort, followed by a nonparametric meta-analysis. RESULTS: In the PNPLA3 gene, rs4823173 (P = 3.44 × 10-10 ), rs2896019 (P = 7.29 × 10-9 ), and rs2281135 (P = 8.73 × 10-9 ) were significantly associated with AST levels. Although not genome-wide significant, these same SNPs were the top hits for ALT (P = 7.12 × 10-8 , P = 1.98 × 10-7 , and P = 1.81 × 10-7 , respectively). The strong correlation (r2 = 1.0) for these SNPs indicated a single hit in the PNPLA3 gene. No genome-wide significant associations were found for GGT. CONCLUSIONS: PNPLA3, a locus previously identified with ALT, AST, and nonalcoholic fatty liver disease in European and Japanese GWAS, is also associated with liver enzymes in populations of Mexican American ancestry.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus.

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

Preferential association of a functional variant in complement receptor 2 with antibodies to double-stranded DNA.

OBJECTIVES: Systemic lupus erythematosus (SLE; OMIM 152700) is characterised by the production of antibodies to nuclear antigens. We previously identified variants in complement receptor 2 (CR2/CD21) that were associated with decreased risk of SLE. This study aimed to identify the causal variant for this association. METHODS: Genotyped and imputed genetic variants spanning CR2 were assessed for association with SLE in 15 750 case-control subjects from four ancestral groups. Allele-specific functional effects of associated variants were determined using quantitative real-time PCR, quantitative flow cytometry, electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP)-PCR. RESULTS: The strongest association signal was detected at rs1876453 in intron 1 of CR2 (pmeta=4.2×10(-4), OR 0.85), specifically when subjects were stratified based on the presence of dsDNA autoantibodies (case-control pmeta=7.6×10(-7), OR 0.71; case-only pmeta=1.9×10(-4), OR 0.75). Although allele-specific effects on B cell CR2 mRNA or protein levels were not identified, levels of complement receptor 1 (CR1/CD35) mRNA and protein were significantly higher on B cells of subjects harbouring the minor allele (p=0.0248 and p=0.0006, respectively). The minor allele altered the formation of several DNA protein complexes by EMSA, including one containing CCCTC-binding factor (CTCF), an effect that was confirmed by ChIP-PCR. CONCLUSIONS: These data suggest that rs1876453 in CR2 has long-range effects on gene regulation that decrease susceptibility to lupus. Since the minor allele at rs1876453 is preferentially associated with reduced risk of the highly specific dsDNA autoantibodies that are present in preclinical, active and severe lupus, understanding its mechanisms will have important therapeutic implications.

Lupus Risk Variant Increases pSTAT1 Binding and Decreases ETS1 Expression.

Genetic variants at chromosomal region 11q23.3, near the gene ETS1, have been associated with systemic lupus erythematosus (SLE), or lupus, in independent cohorts of Asian ancestry. Several recent studies have implicated ETS1 as a critical driver of immune cell function and differentiation, and mice deficient in ETS1 develop an SLE-like autoimmunity. We performed a fine-mapping study of 14,551 subjects from multi-ancestral cohorts by starting with genotyped variants and imputing to all common variants spanning ETS1. By constructing genetic models via frequentist and Bayesian association methods, we identified 16 variants that are statistically likely to be causal. We functionally assessed each of these variants on the basis of their likelihood of affecting transcription factor binding, miRNA binding, or chromatin state. Of the four variants that we experimentally examined, only rs6590330 differentially binds lysate from B cells. Using mass spectrometry, we found more binding of the transcription factor signal transducer and activator of transcription 1 (STAT1) to DNA near the risk allele of rs6590330 than near the non-risk allele. Immunoblot analysis and chromatin immunoprecipitation of pSTAT1 in B cells heterozygous for rs6590330 confirmed that the risk allele increased binding to the active form of STAT1. Analysis with expression quantitative trait loci indicated that the risk allele of rs6590330 is associated with decreased ETS1 expression in Han Chinese, but not other ancestral cohorts. We propose a model in which the risk allele of rs6590330 is associated with decreased ETS1 expression and increases SLE risk by enhancing the binding of pSTAT1.

Genetic Variants Associated With Quantitative Glucose Homeostasis Traits Translate to Type 2 Diabetes in Mexican Americans: The GUARDIAN (Genetics Underlying Diabetes in Hispanics) Consortium.

Insulin sensitivity, insulin secretion, insulin clearance, and glucose effectiveness exhibit strong genetic components, although few studies have examined their genetic architecture or influence on type 2 diabetes (T2D) risk. We hypothesized that loci affecting variation in these quantitative traits influence T2D. We completed a multicohort genome-wide association study to search for loci influencing T2D-related quantitative traits in 4,176 Mexican Americans. Quantitative traits were measured by the frequently sampled intravenous glucose tolerance test (four cohorts) or euglycemic clamp (three cohorts), and random-effects models were used to test the association between loci and quantitative traits, adjusting for age, sex, and admixture proportions (Discovery). Analysis revealed a significant (P < 5.00 × 10(-8)) association at 11q14.3 (MTNR1B) with acute insulin response. Loci with P < 0.0001 among the quantitative traits were examined for translation to T2D risk in 6,463 T2D case and 9,232 control subjects of Mexican ancestry (Translation). Nonparametric meta-analysis of the Discovery and Translation cohorts identified significant associations at 6p24 (SLC35B3/TFAP2A) with glucose effectiveness/T2D, 11p15 (KCNQ1) with disposition index/T2D, and 6p22 (CDKAL1) and 11q14 (MTNR1B) with acute insulin response/T2D. These results suggest that T2D and insulin secretion and sensitivity have both shared and distinct genetic factors, potentially delineating genomic components of these quantitative traits that drive the risk for T2D.

Lupus nephritis susceptibility loci in women with systemic lupus erythematosus.

Lupus nephritis is a manifestation of SLE resulting from glomerular immune complex deposition and inflammation. Lupus nephritis demonstrates familial aggregation and accounts for significant morbidity and mortality. We completed a meta-analysis of three genome-wide association studies of SLE to identify lupus nephritis-predisposing loci. Through genotyping and imputation, >1.6 million markers were assessed in 2000 unrelated women of European descent with SLE (588 patients with lupus nephritis and 1412 patients with lupus without nephritis). Tests of association were computed using logistic regression adjusting for population substructure. The strongest evidence for association was observed outside the MHC and included markers localized to 4q11-q13 (PDGFRA, GSX2; P=4.5×10(-7)), 16p12 (SLC5A11; P=5.1×10(-7)), 6p22 (ID4; P=7.4×10(-7)), and 8q24.12 (HAS2, SNTB1; P=1.1×10(-6)). Both HLA-DR2 and HLA-DR3, two well established lupus susceptibility loci, showed evidence of association with lupus nephritis (P=0.06 and P=3.7×10(-5), respectively). Within the class I region, rs9263871 (C6orf15-HCG22) had the strongest evidence of association with lupus nephritis independent of HLA-DR2 and HLA-DR3 (P=8.5×10(-6)). Consistent with a functional role in lupus nephritis, intra-renal mRNA levels of PDGFRA and associated pathway members showed significant enrichment in patients with lupus nephritis (n=32) compared with controls (n=15). Results from this large-scale genome-wide investigation of lupus nephritis provide evidence of multiple biologically relevant lupus nephritis susceptibility loci.

End-stage renal disease in African Americans with lupus nephritis is associated with APOL1.

OBJECTIVE: Lupus nephritis (LN) is a severe manifestation of systemic lupus erythematosus (SLE) that exhibits familial aggregation and may progress to end-stage renal disease (ESRD). LN is more prevalent among African Americans than among European Americans. This study was undertaken to investigate the hypothesis that the apolipoprotein L1 gene (APOL1) nephropathy risk alleles G1/G2, common in African Americans and rare in European Americans, contribute to the ethnic disparity in risk. METHODS: APOL1 G1 and G2 nephropathy alleles were genotyped in 855 African American SLE patients with LN-ESRD (cases) and 534 African American SLE patients without nephropathy (controls) and tested for association under a recessive genetic model, by logistic regression. RESULTS: Ninety percent of the SLE patients were female. The mean ± SD age at SLE diagnosis was significantly lower in LN-ESRD cases than in SLE non-nephropathy controls (27.3 ± 10.9 years versus 39.5 ± 12.2 years). The mean ± SD time from SLE diagnosis to development of LN-ESRD in cases was 7.3 ± 7.2 years. The G1/G2 risk alleles were strongly associated with SLE-ESRD, with 25% of cases and 12% of controls having 2 nephropathy alleles (odds ratio [OR] 2.57, recessive model P = 1.49 × 10(-9)), and after adjustment for age, sex, and ancestry admixture (OR 2.72, P = 6.23 × 10(-6)). The age-, sex-, and admixture-adjusted population attributable risk for ESRD among patients with G1/G2 polymorphisms was 0.26, compared to 0.003 among European American patients. The mean time from SLE diagnosis to ESRD development was ∼2 years earlier among individuals with APOL1 risk genotypes (P = 0.01). CONCLUSION: APOL1 G1/G2 alleles strongly impact the risk of LN-ESRD in African Americans, as well as the time to progression to ESRD. The high frequency of these alleles in African Americans with near absence in European Americans explains an important proportion of the increased risk of LN-ESRD in African Americans.

Lupus risk variants in the PXK locus alter B-cell receptor internalization.

Genome wide association studies have identified variants in PXK that confer risk for humoral autoimmune diseases, including systemic lupus erythematosus (SLE or lupus), rheumatoid arthritis and more recently systemic sclerosis. While PXK is involved in trafficking of epidermal growth factor Receptor (EGFR) in COS-7 cells, mechanisms linking PXK to lupus pathophysiology have remained undefined. In an effort to uncover the mechanism at this locus that increases lupus-risk, we undertook a fine-mapping analysis in a large multi-ancestral study of lupus patients and controls. We define a large (257kb) common haplotype marking a single causal variant that confers lupus risk detected only in European ancestral populations and spans the promoter through the 3' UTR of PXK. The strongest association was found at rs6445972 with P < 4.62 × 10(-10), OR 0.81 (0.75-0.86). Using stepwise logistic regression analysis, we demonstrate that one signal drives the genetic association in the region. Bayesian analysis confirms our results, identifying a 95% credible set consisting of 172 variants spanning 202 kb. Functionally, we found that PXK operates on the B-cell antigen receptor (BCR); we confirmed that PXK influenced the rate of BCR internalization. Furthermore, we demonstrate that individuals carrying the risk haplotype exhibited a decreased rate of BCR internalization, a process known to impact B cell survival and cell fate. Taken together, these data define a new candidate mechanism for the genetic association of variants around PXK with lupus risk and highlight the regulation of intracellular trafficking as a genetically regulated pathway mediating human autoimmunity.

Variable association of reactive intermediate genes with systemic lupus erythematosus in populations with different African ancestry.

OBJECTIVE: Little is known about the genetic etiology of systemic lupus erythematosus (SLE) in individuals of African ancestry, despite its higher prevalence and greater disease severity. Overproduction of nitric oxide (NO) and reactive oxygen species are implicated in the pathogenesis and severity of SLE, making NO synthases and other reactive intermediate-related genes biological candidates for disease susceptibility. We analyzed variation in reactive intermediate genes for association with SLE in 2 populations with African ancestry. METHODS: A total of 244 single-nucleotide polymorphisms (SNP) from 53 regions were analyzed in non-Gullah African Americans (AA; 1432 cases and 1687 controls) and the genetically more homogeneous Gullah of the Sea Islands of South Carolina (133 cases and 112 controls). Single-marker, haplotype, and 2-locus interaction tests were computed for these populations. RESULTS: The glutathione reductase gene GSR (rs2253409; p = 0.0014, OR 1.26, 95% CI 1.09-1.44) was the most significant single SNP association in AA. In the Gullah, the NADH dehydrogenase NDUFS4 (rs381575; p = 0.0065, OR 2.10, 95% CI 1.23-3.59) and NO synthase gene NOS1 (rs561712; p = 0.0072, OR 0.62, 95% CI 0.44-0.88) were most strongly associated with SLE. When both populations were analyzed together, GSR remained the most significant effect (rs2253409; p = 0.00072, OR 1.26, 95% CI 1.10-1.44). Haplotype and 2-locus interaction analyses also uncovered different loci in each population. CONCLUSION: These results suggest distinct patterns of association with SLE in African-derived populations; specific loci may be more strongly associated within select population groups.

Fine mapping of Xq28: both MECP2 and IRAK1 contribute to risk for systemic lupus erythematosus in multiple ancestral groups.

OBJECTIVES: The Xq28 region containing IRAK1 and MECP2 has been identified as a risk locus for systemic lupus erythematosus (SLE) in previous genetic association studies. However, due to the strong linkage disequilibrium between IRAK1 and MECP2, it remains unclear which gene is affected by the underlying causal variant(s) conferring risk of SLE. METHODS: We fine-mapped ≥136 SNPs in a ∼227 kb region on Xq28, containing IRAK1, MECP2 and seven adjacent genes (L1CAM, AVPR2, ARHGAP4, NAA10, RENBP, HCFC1 and TMEM187), for association with SLE in 15 783 case-control subjects derived from four different ancestral groups. RESULTS: Multiple SNPs showed strong association with SLE in European Americans, Asians and Hispanics at p<5×10(-8) with consistent association in subjects with African ancestry. Of these, six SNPs located in the TMEM187-IRAK1-MECP2 region captured the underlying causal variant(s) residing in a common risk haplotype shared by all four ancestral groups. Among them, rs1059702 best explained the Xq28 association signals in conditional testings and exhibited the strongest p value in transancestral meta-analysis (p(meta )= 1.3×10(-27), OR=1.43), and thus was considered to be the most likely causal variant. The risk allele of rs1059702 results in the amino acid substitution S196F in IRAK1 and had previously been shown to increase NF-κB activity in vitro. We also found that the homozygous risk genotype of rs1059702 was associated with lower mRNA levels of MECP2, but not IRAK1, in SLE patients (p=0.0012) and healthy controls (p=0.0064). CONCLUSIONS: These data suggest contributions of both IRAK1 and MECP2 to SLE susceptibility.

Association of two independent functional risk haplotypes in TNIP1 with systemic lupus erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified >30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting NF-κB signaling. The present study was undertaken to investigate the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway. METHODS: We analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog TNIP2, in case-control populations of diverse ethnic origins. TNIP1, TNIP2, and TAX1BP1 were fine-mapped in a total of 8,372 SLE cases and 7,492 healthy controls from European-ancestry, African American, Hispanic, East Asian, and African American Gullah populations. Levels of TNIP1 messenger RNA (mRNA) and ABIN1 protein in Epstein-Barr virus-transformed human B cell lines were analyzed by quantitative reverse transcription-polymerase chain reaction and Western blotting, respectively. RESULTS: We found significant associations between SLE and genetic variants within TNIP1, but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified 2 independent risk haplotypes within TNIP1 in individuals of European ancestry that were also present in African American and Hispanic populations. Levels of TNIP1 mRNA and ABIN1 protein were reduced among subjects with these haplotypes, suggesting that they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression. CONCLUSION: Our results confirm the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis.

Identification of IRF8, TMEM39A, and IKZF3-ZPBP2 as susceptibility loci for systemic lupus erythematosus in a large-scale multiracial replication study.

Systemic lupus erythematosus (SLE) is a chronic heterogeneous autoimmune disorder characterized by the loss of tolerance to self-antigens and dysregulated interferon responses. The etiology of SLE is complex, involving both heritable and environmental factors. Candidate-gene studies and genome-wide association (GWA) scans have been successful in identifying new loci that contribute to disease susceptibility; however, much of the heritable risk has yet to be identified. In this study, we sought to replicate 1,580 variants showing suggestive association with SLE in a previously published GWA scan of European Americans; we tested a multiethnic population consisting of 7,998 SLE cases and 7,492 controls of European, African American, Asian, Hispanic, Gullah, and Amerindian ancestry to find association with the disease. Several genes relevant to immunological pathways showed association with SLE. Three loci exceeded the genome-wide significance threshold: interferon regulatory factor 8 (IRF8; rs11644034; p(meta-Euro) = 2.08 × 10(-10)), transmembrane protein 39A (TMEM39A; rs1132200; p(meta-all) = 8.62 × 10(-9)), and 17q21 (rs1453560; p(meta-all) = 3.48 × 10(-10)) between IKAROS family of zinc finger 3 (AIOLOS; IKZF3) and zona pellucida binding protein 2 (ZPBP2). Fine mapping, resequencing, imputation, and haplotype analysis of IRF8 indicated that three independent effects tagged by rs8046526, rs450443, and rs4843869, respectively, were required for risk in individuals of European ancestry. Eleven additional replicated effects (5 × 10(-8) < p(meta-Euro) < 9.99 × 10(-5)) were observed with CFHR1, CADM2, LOC730109/IL12A, LPP, LOC63920, SLU7, ADAMTSL1, C10orf64, OR8D4, FAM19A2, and STXBP6. The results of this study increase the number of confirmed SLE risk loci and identify others warranting further investigation.

Evaluation of TRAF6 in a large multiancestral lupus cohort.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease with significant immune system aberrations resulting from complex heritable genetics as well as environmental factors. We undertook to study the role of TRAF6 as a candidate gene for SLE, since it plays a major role in several signaling pathways that are important for immunity and organ development. METHODS: Fifteen single-nucleotide polymorphisms (SNPs) across TRAF6 were evaluated in 7,490 SLE patients and 6,780 control subjects from different ancestries. Population-based case-control association analyses and meta-analyses were performed. P values, false discovery rate q values, and odds ratios (ORs) with 95% confidence intervals (95% CIs) were calculated. RESULTS: Evidence of associations was detected in multiple SNPs. The best overall P values were obtained for SNPs rs5030437 and rs4755453 (P = 7.85 × 10(-5) and P = 4.73 × 10(-5) , respectively) without significant heterogeneity among populations (P = 0.67 and P = 0.50, respectively, in Q statistic). In addition, SNP rs540386, which was previously reported to be associated with rheumatoid arthritis (RA), was found to be in linkage disequilibrium with these 2 SNPs (r(2) = 0.95) and demonstrated evidence of association with SLE in the same direction (meta-analysis P = 9.15 × 10(-4) , OR 0.89 [95% CI 0.83-0.95]). The presence of thrombocytopenia improved the overall results in different populations (meta-analysis P = 1.99 × 10(-6) , OR 0.57 [95% CI 0.45-0.72], for rs5030470). Finally, evidence of family-based association in 34 African American pedigrees with the presence of thrombocytopenia was detected in 1 available SNP (rs5030437) with a Z score magnitude of 2.28 (P = 0.02) under a dominant model. CONCLUSION: Our data indicate the presence of association of TRAF6 with SLE, consistent with the previous report of association with RA. These data provide further support for the involvement of TRAF6 in the pathogenesis of autoimmunity.

Evidence for gene-gene epistatic interactions among susceptibility loci for systemic lupus erythematosus.

OBJECTIVE: Several confirmed genetic susceptibility loci for lupus have been described. To date, no clear evidence for genetic epistasis in lupus has been established. The aim of this study was to test for gene-gene interactions in a number of known lupus susceptibility loci. METHODS: Eighteen single-nucleotide polymorphisms tagging independent and confirmed lupus susceptibility loci were genotyped in a set of 4,248 patients with lupus and 3,818 normal healthy control subjects of European descent. Epistasis was tested by a 2-step approach using both parametric and nonparametric methods. The false discovery rate (FDR) method was used to correct for multiple testing. RESULTS: We detected and confirmed gene-gene interactions between the HLA region and CTLA4, IRF5, and ITGAM and between PDCD1 and IL21 in patients with lupus. The most significant interaction detected by parametric analysis was between rs3131379 in the HLA region and rs231775 in CTLA4 (interaction odds ratio 1.19, Z = 3.95, P = 7.8 × 10(-5) [FDR ≤0.05], P for multifactor dimensionality reduction = 5.9 × 10(-45)). Importantly, our data suggest that in patients with lupus, the presence of the HLA lupus risk alleles in rs1270942 and rs3131379 increases the odds of also carrying the lupus risk allele in IRF5 (rs2070197) by 17% and 16%, respectively (P = 0.0028 and P = 0.0047, respectively). CONCLUSION: We provide evidence for gene-gene epistasis in systemic lupus erythematosus. These findings support a role for genetic interaction contributing to the complexity of lupus heritability.

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(-06)) 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.

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(-8), OR = 1.18) and an intergenic SNP between CFHR1 and CFHR4 (rs16840639, P(meta) = 2.9×10(-7), 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(-7), OR = 1.47) conferred a higher risk of SLE than heterozygous deletion (P(meta) = 3.5×10(-4), 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.

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.

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⁻4), 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.