Association of the HLA locus and TNF with type I autoimmune hepatitis susceptibility in New Zealand Caucasians

PURPOSE: The precise etiology of autoimmune hepatitis (AIH) remains unknown, although a number of genetic loci have been implicated in the susceptibility of type 1 AIH. The purpose of this study was to test for association of these loci with type 1 AIH in New Zealand Caucasians. METHODS: 77 AIH patients and 485 healthy controls were genotyped for the SNPs rs2187668 (HLA-DRB*03:01), rs660895 (HLA-DRB*04:01), rs3749971 (HLA-A1-B8-DR3), rs231775 (CLTLA4), rs1800629 (TNF), and rs1800682 (FAS) using predesigned TaqMan SNP genotyping assays. Chi square analysis was used to test for association of allele and genotype with overall AIH, and with severe fibrosis and ALT levels at 6 months. RESULTS: Significant risk of AIH was conferred by the minor alleles of rs2187668 (OR = 2.45, 95% CI 1.65-3.61, p < 0.0001), rs3749971 (OR = 1.89, 95% CI 1.21-2.94, p = 0.004) and rs1800629 (OR = 2.06, 95% CI 1.41-3.01, p = 0.0001). Multivariate analysis showed that rs2187668 was independently associated with type 1 AIH susceptibility (OR = 2.40, 95% CI 1.46-3.93, p = 0.001). The C allele of FAS SNP rs1800682 was associated with increased risk of severe fibrosis at diagnosis (OR = 2.03, 95% CI 1.05-3.93, p = 0.035) and with incomplete normalization of ALT levels at 6 months post-diagnosis (OR = 3.94, 95% CI 1.62-9.54, p = 0.0015). CONCLUSIONS: This is the first population-based study to investigate genetic risk loci for type 1 AIH in New Zealand Caucasians. We report significant independent association of HLA-DRB1*03:01 with overall susceptibility to type 1 AIH, as well as FAS with a more aggressive disease phenotype

The PTPN22 Locus and Rheumatoid Arthritis: No Evidence for an Effect on Risk Independent of Arg620Trp

The Trp(620) allotype of PTPN22 confers susceptibility to rheumatoid arthritis (RA) and certain other classical autoimmune diseases. There has been a report of other variants within the PTPN22 locus that alter risk of RA; protective haplotype '5', haplotype group '6-10' and susceptibility haplotype '4', suggesting the possibility of other PTPN22 variants involved in the pathogenesis of RA independent of R620W (rs2476601). Our aim was to further investigate this possibility.A total of 4,460 RA cases and 4,481 controls, all European, were analysed. Single nucleotide polymorphisms rs3789607, rs12144309, rs3811021 and rs12566340 were genotyped over New Zealand (NZ) and UK samples. Publically-available Wellcome Trust Case Control Consortium (WTCCC) genotype data were used.The protective effect of haplotype 5 was confirmed (rs3789607; (OR = 0.91, P = 0.016), and a second protective effect (possibly of haplotype 6) was observed (rs12144309; OR = 0.90, P = 0.021). The previously reported susceptibility effect of haplotype 4 was not replicated; instead a protective effect was observed (rs3811021; OR = 0.85, P = 1.4×10(-5)). Haplotypes defined by rs3789607, rs12144309 and rs3811021 coalesced with the major allele of rs12566340 within the adjacent BFK (B-cell lymphoma 2 (BCL2) family kin) gene. We, therefore, tested rs12566340 for association with RA conditional on rs2476601; there was no evidence for an independent effect at rs12566340 (P = 0.76). Similarly, there was no evidence for an independent effect at rs12566340 in type 1 diabetes (P = 0.85).We have no evidence for a common variant additional to rs2476601 within the PTPN22 locus that influences the risk of RA. Arg620Trp is almost certainly the single common causal variant

The ITGAV rs3738919 variant and susceptibility to rheumatoid arthritis in four Caucasian sample sets

INTRODUCTION: Angiogenesis is an important process in the development of destructive synovial pannus in rheumatoid arthritis (RA). The ITGAV +gene encodes a cell cycle-associated antigen, integrin alphanubeta 3, which plays a role in RA angiogenesis. Previously, two independent studies identified an association between the major allele of the ITGAV single-nucleotide polymorphism (SNP) rs3738919 and RA. We therefore tested this association in an independent study using New Zealand (NZ) and Oxford (UK) RA case control samples. METHODS: We compared genotype frequencies in 740 NZ Caucasian RA patients and 553 controls genotyped for rs3738919, using a polymerase chain reaction-restriction fragment length polymorphism assay. A TaqMan genotyping SNP assay was used to type 713 Caucasian RA patients and 515 control samples from Oxford for the rs3738919 variant. Association of rs3738919 with RA was tested in these two sample sets using the chi-square goodness-of-fit test. The Mantel-Haenszel test was used to perform a meta-analysis, combining the genetic results from four independent Caucasian case control cohorts, consisting of 3,527 cases and 4,126 controls. Haplotype analysis was also performed using SNPs rs3911238, rs10174098 and rs3738919 in the Wellcome Trust Case Control Consortium, NZ and Oxford case control samples. RESULTS: We found no evidence for association between ITGAV and RA in either the NZ or Oxford sample set (odds ratio [OR] = 0.88, P(allelic) = 0.11 and OR = 1.18, P(allelic) = 0.07, respectively). Inclusion of these data in a meta-analysis (random effects) of four independent cohorts (3,527 cases and 4,126 controls) weakens support for the hypothesis that rs3738919 plays a role in the development of RA (OR(combined) = 0.92, 95% confidence interval 0.80 to 1.07; P = 0.29). No consistent haplotype associations were evident. CONCLUSIONS: Association of ITGAV SNP rs7378919 with RA was not replicated in NZ or Oxford case control sample sets. Meta-analysis of these and previously published data lends limited support for a role for the ITGAV in RA in Caucasians of European ancestry

The Oxytocin Receptor Gene (OXTR) Variant rs53576 Is Not Related to Emotional Traits or States in Young Adults

Background: To understand the genetic underpinnings of emotion, researchers have studied genetic variants in the oxytocin system, a hormone and neurotransmitter important to socio-emotional functioning. The oxytocin receptor gene (OXTR) variant rs53576 has been associated with emotional traits such as positive affect and related constructs such as optimism and self-esteem. Individuals carrying the A allele (AG and AA genotypes) of rs53576 have been found to score lower in these traits when compared to GG homozygotes, although not always. Given recent mixed evidence regarding this polymorphism, replication of these associations is critical.Methods: Using a cross-sectional design, the present study tested the association between rs53576 and a wide variety of emotional traits and states in a sample of 611 young adults ages 18 – 25 of various ethnicities (European, Asian, Māori/Pacific Islander, other). Participants completed standard trait measures of positive and negative affect, depressive symptoms, life engagement, psychological well-being, optimism, and self-esteem. They also completed state measures of positive and negative affect and life engagement for 13-days using Internet daily diaries.Results: Controlling for ethnicity and gender, variation at the OXTR variant rs53576 obtained from blood samples was not related to any of the emotional traits or states. This null finding occurred despite measuring emotions in “near to real time” using daily diaries and having sufficient power to detect a medium effect size difference between homozygous genotype groups.Conclusion: These findings suggest that variation at the rs53576 locus may not be as involved in emotional differences as initial studies suggested

Abundant local interactions in the 4p16.1 region suggest functional mechanisms underlying SLC2A9 associations with human serum uric acid

Human serum uric acid concentration (SUA) is a complex trait. A recent meta-analysis of multiple genome-wide association studies (GWAS) identified 28 loci associated with SUA jointly explaining only 7.7% of the SUA variance, with 3.4% explained by two major loci (SLC2A9 and ABCG2). Here we examined whether gene–gene interactions had any roles in regulating SUA using two large GWAS cohorts included in the meta-analysis [the Atherosclerosis Risk in Communities study cohort (ARIC) and the Framingham Heart Study cohort (FHS)]. We found abundant genome-wide significant local interactions in ARIC in the 4p16.1 region located mostly in an intergenic area near SLC2A9 that were not driven by linkage disequilibrium and were replicated in FHS. Taking the forward selection approach, we constructed a model of five SNPs with marginal effects and three epistatic SNP pairs in ARIC—three marginal SNPs were located within SLC2A9 and the remaining SNPs were all located in the nearby intergenic area. The full model explained 1.5% more SUA variance than that explained by the lead SNP alone, but only 0.3% was contributed by the marginal and epistatic effects of the SNPs in the intergenic area. Functional analysis revealed strong evidence that the epistatically interacting SNPs in the intergenic area were unusually enriched at enhancers active in ENCODE hepatic (HepG2, P = 4.7E−05) and precursor red blood (K562, P = 5.0E−06) cells, putatively regulating transcription of WDR1 and SLC2A9. These results suggest that exploring epistatic interactions is valuable in uncovering the complex functional mechanisms underlying the 4p16.1 region

A functional genomics approach reveals suggestive quantitative trait loci associated with combined TLR4 and BCP crystal-induced inflammation and osteoarthritis

Objective: Basic calcium phosphate (BCP) crystals can activate the NLRP3 inflammasome and are potentially involved in the pathogenesis of osteoarthritis (OA). In order to elucidate relevant inflammatory mechanisms in OA, we used a functional genomics approach to assess genetic variation influencing BCP crystal-induced cytokine production. Method: Peripheral blood mononuclear cells (PBMCs) were isolated from healthy volunteers who were previously genotyped and stimulated with BCP crystals and/or lipopolysaccharide (LPS) after which cytokines release was assessed. Cytokine quantitative trait locus (cQTL) mapping was performed. For in vitro validation of the cQTL located in anoctamin 3 (ANO3), PBMCs were incubated with Tamoxifen and Benzbromarone prior to stimulation. Additionally, we performed co-localisation analysis of our top cQTLs with the most recent OA meta-analysis of genome-wide association studies (GWAS). Results: We observed that BCP crystals and LPS synergistically induce IL-1β in human PBMCs. cQTL analysis revealed several suggestive loci influencing cytokine release upon stimulation, among which are quantitative trait locus annotated to ANO3 and GLIS3. As functional validation, anoctamin inhibitors reduced IL-1β release in PBMCs after stimulation. Co-localisation analysis showed that the GLIS3 locus was shared between LPS/BCP crystal-induced IL-1β and genetic association with Knee OA. Conclusions: We identified and functionally validated a new locus, ANO3, associated with LPS/BCP crystal-induced inflammation in PBMCs. Moreover, the cQTL in the GLIS3 locus co-localises with the previously found locus associated with Knee OA, suggesting that this Knee OA locus might be explained through an inflammatory mechanism. These results form a basis for further exploration of inflammatory mechanisms in OA.</p

A genome-wide association study of rheumatoid arthritis without antibodies against citrullinated peptides

Introduction. Rheumatoid arthritis (RA) patients can be classified based on presence or absence of anticitrullinated peptide antibodies (ACPA) in their serum. This heterogeneity among patients may reflect important biological differences underlying the disease process. To date, the majority of genetic studies have focused on the ACPA-positive group. Therefore, our goal was to analyse the genetic risk factors that contribute to ACPA-negative RA. Methods. We performed a large-scale genome-wide association study (GWAS) in three Caucasian European cohorts comprising 1148 ACPA-negative RA patients and 6008 controls. All patients were screened using the Illumina Human Cyto-12 chip, and controls were genotyped using different genome-wide platforms. Population-independent analyses were carried out by means of logistic regression. Meta-analysis with previously published data was performed as follow-up for selected signals (reaching a total of 1922 ACPA-negative RA patients and 7087 controls). Imputation of classical HLA alleles, aminoacid residues and single nucleotide polymorphisms was undertaken. Results. The combined analysis of the studied cohorts resulted in identification of a peak of association in the HLA-region and several suggestive non-HLA associations. Meta-analysis with previous reports confirmed the association of the HLA region with this subset and an observed association in the CLYBL locus remained suggestive. The imputation and deep interrogation of the HLA region led to identification of a two aminoacid model (HLA-B at position 9 and HLA-DRB1 at position 11) that accounted for the observed genome-wide associations in this region. Conclusions. Our study shed light on the influence of the HLA region in ACPA-negative RA and identified a suggestive risk locus for this condition

Multiplicative interaction of functional inflammasome genetic variants in determining the risk of gout

Introduction: The acute gout flare results from a localised self-limiting innate immune response to monosodium urate (MSU) crystals deposited in joints in hyperuricaemic individuals. Activation of the caspase recruitment domain-containing protein 8 (CARD8) NOD-like receptor pyrin-containing 3 (NLRP3) inflammasome by MSU crystals and production of mature interleukin-1 beta (IL-1 beta) is central to acute gouty arthritis. However very little is known about genetic control of the innate immune response involved in acute gouty arthritis. Therefore our aim was to test functional single nucleotide polymorphism (SNP) variants in the toll-like receptor (TLR)-inflammasome-IL-1 beta axis for association with gout. Methods: 1,494 gout cases of European and 863 gout cases of New Zealand (NZ) Polynesian (Maori and Pacific Island) ancestry were included. Gout was diagnosed by the 1977 ARA gout classification criteria. There were 1,030 Polynesian controls and 10,942 European controls including from the publicly-available Atherosclerosis Risk in Communities (ARIC) and Framingham Heart (FHS) studies. The ten SNPs were either genotyped by Sequenom MassArray or by Affymetrix SNP array or imputed in the ARIC and FHS datasets. Allelic association was done by logistic regression adjusting by age and sex with European and Polynesian data combined by meta-analysis. Sample sets were pooled for multiplicative interaction analysis, which was also adjusted by sample set. Results: Eleven SNPs were tested in the TLR2, CD14, IL1B, CARD8, NLRP3, MYD88, P2RX7, DAPK1 and TNXIP genes. Nominally significant (P <0.05) associations with gout were detected at CARD8 rs2043211 (OR = 1.12, P = 0.007), IL1B rs1143623 (OR = 1.10, P = 0.020) and CD14 rs2569190 (OR = 1.08; P = 0.036). There was significant multiplicative interaction between CARD8 and IL1B (P = 0.005), with the IL1B risk genotype amplifying the risk effect of CARD8. Conclusion: There is evidence for association of gout with functional variants in CARD8, IL1B and CD14. The gout-associated allele of IL1B increases expression of IL-1 beta-the multiplicative interaction with CARD8 would be consistent with a synergy of greater inflammasome activity (resulting from reduced CARD8) combined with higher levels of pre-IL-1 beta expression leading to increased production of mature IL-1 beta in gout