Identification of a regulatory pathway governing TRAF1 via an arthritis-associated non-coding variant.

TRAF1/C5 was among the first loci shown to confer risk for inflammatory arthritis in the absence of an associated coding variant, but its genetic mechanism remains undefined. Using Immunochip data from 3,939 patients with juvenile idiopathic arthritis (JIA) and 14,412 control individuals, we identified 132 plausible common non-coding variants, reduced serially by single-nucleotide polymorphism sequencing (SNP-seq), electrophoretic mobility shift, and luciferase studies to the single variant rs7034653 in the third intron of TRAF1. Genetically manipulated experimental cells and primary monocytes from genotyped donors establish that the risk G allele reduces binding of Fos-related antigen 2 (FRA2), encoded by FOSL2, resulting in reduced TRAF1 expression and enhanced tumor necrosis factor (TNF) production. Conditioning on this JIA variant eliminated attributable risk for rheumatoid arthritis, implicating a mechanism shared across the arthritis spectrum. These findings reveal that rs7034653, FRA2, and TRAF1 mediate a pathway through which a non-coding functional variant drives risk of inflammatory arthritis in children and adults.

Association of SLCO1B1 *14 Allele with Poor Response to Methotrexate in Juvenile Idiopathic Arthritis Patients.

OBJECTIVE: Variants in the SLCO1B1 gene, encoding a hepatic methotrexate (MTX) transporter, affect clearance of high-dose MTX. We tested whether in the *14 and *15 alleles of SLCO1B1 influenced the response to low-dose MTX in juvenile idiopathic arthritis (JIA) patients. METHODS: The study included 310 JIA patients genotyped for three single nucleotide polymorphisms (SNPs) in SLCO1B1 (rs4149056, rs2306283, and rs11045819). A patient's SLCO1B1 diplotype was determined by combining the SNPs into the *1a, *1b, *4, *5, *14, and *15 alleles. Number of active joints at follow-up (visit closest to 6 months of treatment and prior to starting a tumor necrosis factor inhibitor) was used as the dependent variable in a negative binomial regression model that included active joint count at baseline as a covariate. RESULTS: The SLCO1B1*14 allele was associated with less response to MTX (P = 0.024) and the *15 allele was not associated with response to MTX (P = 0.392). CONCLUSION: SLCO1B1 alleles may be associated with poor response to MTX in JIA patients. The *14 allele has been associated with fast clearance (low exposure) after high-dose MTX in patients with leukemia. Thus, the SLCO1B1 gene may be informative for precision dosing of MTX in JIA patients. Patients carrying the *14 allele may require a higher dose than noncarriers to achieve a similar response to MTX.

Trans-Ethnic Mapping of BANK1 Identifies Two Independent SLE-Risk Linkage Groups Enriched for Co-Transcriptional Splicing Marks.

BANK1 is a susceptibility gene for several systemic autoimmune diseases in several populations. Using the genome-wide association study (GWAS) data from Europeans (EUR) and African Americans (AA), we performed an extensive fine mapping of ankyrin repeats 1 (BANK1). To increase the SNP density, we used imputation followed by univariate and conditional analysis, combined with a haplotypic and expression quantitative trait locus (eQTL) analysis. The data from Europeans showed that the associated region was restricted to a minimal and dependent set of SNPs covering introns two and three, and exon two. In AA, the signal found in the Europeans was split into two independent effects. All of the major risk associated SNPs were eQTLs, and the risks were associated with an increased BANK1 gene expression. Functional annotation analysis revealed the enrichment of repressive B cell epigenomic marks (EZH2 and H3K27me3) and a strong enrichment of splice junctions. Furthermore, one eQTL located in intron two, rs13106926, was found within the binding site for RUNX3, a transcriptional activator. These results connect the local genome topography, chromatin structure, and the regulatory landscape of BANK1 with co-transcriptional splicing of exon two. Our data defines a minimal set of risk associated eQTLs predicted to be involved in the expression of BANK1 modulated through epigenetic regulation and splicing. These findings allow us to suggest that the increased expression of BANK1 will have an impact on B-cell mediated disease pathways.

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.

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.

Genetic variation in the peroxisome proliferator activated receptor-gamma gene is associated with histologically advanced NAFLD.

BACKGROUND: The peroxisome proliferator activated receptor-gamma (PPARG) is a nuclear receptor that regulates adipocyte differentiation, insulin sensitivity and lipid metabolism, thus, it represents a good candidate gene for non-alcoholic fatty liver disease (NAFLD). PURPOSE AND METHOD: We investigated the association of two PPARG variants (Pro12Ala and C1431T) with NAFLD and its histological features. DNA was extracted from 274 archived, formalin-fixed liver biopsy specimens from 212 patients with NAFLD and 62 controls with normal liver histology. RESULTS: Individual SNPs did not show significant association with NAFLD or its histological features. A haplotype comprised of both minor alleles (GT) was less enriched whereas a haplotype comprised of the two major alleles (CC) was more enriched in subjects with NAFLD compared to controls [9.3% vs. 28.1% for GT (P = 0.001, OR 0.26 (range 0.14-0.48) and 80.4% vs. 64.8% for CC (P = 0.037, OR 2.23 (range 1.30-3.81)]. Both haplotypes were significantly associated with steatosis and fibrosis. The GT haplotype was also associated with lobular inflammation. CONCLUSIONS: Genetic variation in PPARG is associated with NAFLD, and the GT haplotype is associated with inflammatory and fibrotic changes that denote histologically advanced NAFLD.

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.

Identification of candidate loci at 6p21 and 21q22 in a genome-wide association study of cardiac manifestations of neonatal lupus.

OBJECTIVE: Cardiac manifestations of neonatal lupus, comprising atrioventricular conduction defects and cardiomyopathy, occur in fetuses exposed to anti-Ro/SSA antibodies, and carry substantial mortality. There is strong evidence of a genetic contribution to the risk. This study was undertaken to evaluate single-nucleotide polymorphisms (SNPs) for associations with cardiac neonatal lupus. METHODS: Children of European ancestry with cardiac neonatal lupus (n = 116) were genotyped using the Illumina 370K SNP platform and merged with 3,351 controls. Odds ratios (ORs) and 95% confidence intervals (95% CIs) for association with cardiac neonatal lupus were determined. RESULTS: The 17 most significant associations with cardiac neonatal lupus were found in the HLA region. The region near the MICB gene showed the strongest variant (rs3099844; P(dom) = 4.52 × 10(-10) , OR 3.34 [95% CI 2.29-4.89]), followed by a missense variant within C6orf10 (rs7775397; P(dom) = 1.35 × 10(-9) , OR 3.30), which lies between NOTCH4 and BTNL2, and several SNPs near the tumor necrosis factor α gene, including rs2857595 (P(add) = 1.96 × 10(-9) , OR 2.37), rs2230365 (P(add) = 1.00 × 10(-3) , OR 0.46), and rs3128982 (P(add) = 6.40 × 10(-6) , OR 1.86). Outside the HLA region, an association was detected at 21q22, upstream of the transcription regulator ets-related isoform 1 (rs743446; P = 5.45 × 10(-6) , OR 2.40). HLA notwithstanding, no individual locus previously implicated in autoimmune diseases achieved genome-wide significance. CONCLUSION: These results suggest that variation near genes related to inflammatory and apoptotic responses may promote cardiac injury initiated by passively acquired autoantibodies.