A bioinformatics workflow for detecting signatures of selection in genomic data

The detection of “signatures of selection” is now possible on a genome-wide scale in many plant and animal species, and can be performed in a population-specific manner due to the wealth of per-population genome-wide genotype data that is available. With genomic regions that exhibit evidence of having been under selection shown to also be enriched for genes associated with biologically important traits, detection of evidence of selective pressure is emerging as an additional approach for identifying novel gene-trait associations. While high-density genotype data is now relatively easy to obtain, for many researchers it is not immediately obvious how to go about identifying signatures of selection in these data sets. Here we describe a basic workflow, constructed from open source tools, for detecting and examining evidence of selection in genomic data. Code to install and implement the pipeline components, and instructions to run a basic analysis using the workflow described here, can be downloaded from our public GitHub repository: http://www.github.com/smilefreak/selectionTools

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

What predicts regression from pre-diabetes to normal glucose regulation following a primary care nurse-delivered dietary intervention? A study protocol for a prospective cohort study

Introduction Pre-diabetes is a high-risk state for the development of type 2 diabetes mellitus (T2DM) and cardiovascular disease. Regression to normoglycaemia, even if transient, significantly reduces the risk of developing T2DM. The primary aim of this mixed-methods study is to determine if there are clinically relevant differences among those with pre-diabetes and excess weight who regress to normoglycaemia, those who have persistent pre-diabetes and those who progress to T2DM following participation in a 6-month primary care nurse-delivered pre-diabetes dietary intervention. Incidence of T2DM at 2 years will be examined. Methods and analysis Four hundred participants with pre-diabetes (New Zealand definition glycated haemoglobin 41–49 mmol/mol) and a body mass index \u3e25 kg/m2 will be recruited through eight primary care practices in Hawke’s Bay, New Zealand. Trained primary care nurses will deliver a 6-month structured dietary intervention, followed by quarterly reviews for 18 months post-intervention. Clinical data, data on lifestyle factors and health-related quality of life (HR-QoL) and blood samples will be collected at baseline, 6 months, 12 months and 24 months. Sixty participants purposefully selected will complete a semi-structured interview following the 6-month intervention. Poisson regression with robust standard errors and clustered by practice will be used to identify predictors of regression or progression at 6 months, and risk factors for developing T2DM at 2 years. Qualitative data will be analysed thematically. Changes in HR-QoL will be described and potential cost savings will be estimated from a funder’s perspective at 2 years. Ethics and dissemination This study was approved by the Northern A Health and Disability Ethics Committee, New Zealand (Ethics Reference: 17/NTA/24). Study results will be presented to participants, published in peer-reviewed journals and presented at relevant conferences. Trial registration number ACTRN12617000591358; Pre-results

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

Only one independent genetic association with rheumatoid arthritis within the KIAA1109-TENR-IL2-IL21 locus in Caucasian sample sets: confirmation of association of rs6822844 with rheumatoid arthritis at a genome-wide level of significance

INTRODUCTION: The single nucleotide polymorphism (SNP) rs6822844 within the KIAA1109-TENR-IL2-IL21 gene cluster has been associated with rheumatoid arthritis (RA). Other variants within this cluster, including rs17388568 that is not in linkage disequilibrium (LD) with rs6822844, and rs907715 that is in moderate LD with rs6822844 and rs17388568, have been associated with a number of autoimmune phenotypes, including type 1 diabetes (T1D). Here we aimed to: one, confirm at a genome-wide level of significance association of rs6822844 with RA and, two, evaluate whether or not there were effects independent of rs6822844 on RA at the KIAA1109-TENR-IL2-IL21 locus. METHODS: A total of 842 Australasian RA patients and 1,115 controls of European Caucasian ancestry were genotyped for rs6822844, rs17388568 and rs907715. Meta-analysis of these data with published and publicly-available data was conducted using STATA. RESULTS: No statistically significant evidence for association was observed in the Australasian sample set for rs6822844 (odds ratio (OR)=0.95 (0.80 to 1.12), P=0.54), or rs17388568 (OR=1.03 (0.90 to 1.19), P=0.65) or rs907715 (OR=0.98 (0.86 to 1.12), P=0.69). When combined in a meta-analysis using data from a total of 9,772 cases and 10,909 controls there was a genome-wide level of significance supporting association of rs6822844 with RA (OR=0.86 (0.82 to 0.91), P=8.8x10(-8), P=2.1x10(-8) including North American Rheumatoid Arthritis Consortium data). Meta-analysis of rs17388568, using a total of 6,585 cases and 7,528 controls, revealed no significant association with RA (OR=1.03, (0.98 to 1.09); P=0.22) and meta-analysis of rs907715 using a total of 2,689 cases and 4,045 controls revealed a trend towards association (OR=0.93 (0.87 to 1.00), P=0.07). However, this trend was not independent of the association at rs6822844. CONCLUSIONS: The KIAA1109-TENR-IL2-IL21 gene cluster, that encodes an interleukin (IL-21) that plays an important role in Th17 cell biology, is the 20th locus for which there is a genome-wide (P ≤ 5x10(-8)) level of support for association with RA. As for most other autoimmune diseases, with the notable exception of T1D, rs6822844 is the dominant association in the locus. The KIAA1109-TENR-IL2-IL21 locus also confers susceptibility to other autoimmune phenotypes with a heterogeneous pattern of association

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