Tugce Karaderi

Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility

Ankylosing spondylitis is a common form of inflammatory arthritis predominantly affecting the spine and pelvis that occurs in approximately 5 out of 1,000 adults of European descent. Here we report the identification of three variants in the RUNX3, LTBR-TNFRSF1A and IL12B regions convincingly associated with ankylosing spondylitis (P < 5 × 10−8 in the combined discovery and replication datasets) and a further four loci at PTGER4, TBKBP1, ANTXR2 and CARD9 that show strong association across all our datasets (P < 5 × 10−6 overall, with support in each of the three datasets studied). We also show that polymorphisms of ERAP1, which encodes an endoplasmic reticulum aminopeptidase involved in peptide trimming before HLA class I presentation, only affect ankylosing spondylitis risk in HLA-B27–positive individuals. These findings provide strong evidence that HLA-B27 operates in ankylosing spondylitis through a mechanism involving aberrant processing of antigenic peptides.

Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci

Ankylosing spondylitis is a common, highly heritable inflammatory arthritis affecting primarily the spine and pelvis. In addition to HLA-B*27 alleles, 12 loci have previously been identified that are associated with ankylosing spondylitis in populations of European ancestry, and 2 associated loci have been identified in Asians. In this study, we used the Illumina Immunochip microarray to perform a case-control association study involving 10,619 individuals with ankylosing spondylitis (cases) and 15,145 controls. We identified 13 new risk loci and 12 additional ankylosing spondylitis–associated haplotypes at 11 loci. Two ankylosing spondylitis–associated regions have now been identified encoding four aminopeptidases that are involved in peptide processing before major histocompatibility complex (MHC) class I presentation. Protective variants at two of these loci are associated both with reduced aminopeptidase function and with MHC class I cell surface expression.

Association between the interleukin 23 receptor and ankylosing spondylitis is confirmed by a new UK case–control study and meta-analysis of published series

OBJECTIVES It has been shown previously that IL-23R variants are associated with AS. We conducted an extended analysis in the UK population and a meta-analysis with the previously published studies, in order to refine these IL-23R associations with AS. METHODS The UK case-control study included 730 new cases and 1331 healthy controls. In the extended study, the 730 cases were combined with 1088 published cases. Allelic associations were analysed using contingency tables. In the meta-analysis, 3482 cases and 3150 controls from four different published studies and the new UK cases were combined. DerSimonian-Laird test was used to calculate random effects pooled odds ratios (ORs). RESULTS In the UK case-control study with new cases, four of the eight SNPs showed significant associations, whereas in the extended UK study, seven of the eight IL-23R SNPs showed significant associations (P < 0.05) with AS, maximal with rs11209032 (P < 10(-5), OR 1.3), when cases with IBD and/or psoriasis were excluded. The meta-analysis showed significant associations with all eight SNPs; the strongest associations were again seen not only with rs11209032 (P = 4.06 x 10(-9), OR approximately 1.2) but also with rs11209026 (P < 10(-10), OR approximately 0.6). CONCLUSIONS IL-23R polymorphisms are clearly associated with AS, but the primary causal association(s) is(are) still not established. These polymorphisms could contribute either increased or decreased susceptibility to AS; functional studies will be required for their full evaluation. Additionally, observed stronger associations with SNPs rs11209026 and rs11465804 upon exclusion of IBD and/or psoriasis cases may represent an independent association with AS.

Elucidating the chromosome 9 association with AS; CARD9 is a candidate gene

Ankylosing spondylitis (AS) is polygenic with contributions from the immunologically relevant genes HLA-B*27, ERAP1 and IL23R. A recent genome-wide association screen (GWAS) identified associations (P∼0.005) with the non-synonymous single-nucleotide polymorphisms (nsSNPs), rs4077515 and rs3812571, in caspase recruitment domain-containing protein 9 (CARD9) and small nuclear RNA-activating complex polypeptide 4 (SNAPC4) on chromosome 9q that had previously been linked to AS. We replicated these associations in a study of 730 AS patients compared with 2879 historic disease controls (rs4077515 P=0.0004, odds ratio (OR)=1.2, 95% confidence interval (CI)=1.1–1.4; rs3812571 P=0.0003, OR=1.2, 95% CI=1.1–1.4). Meta-analysis revealed strong associations of both SNPs with AS, rs4077515 P=0.000005, OR=1.2, 95% CI=1.1–1.3 and rs3812571 P=0.000006, OR=1.2, 95% CI=1.1–1.3. We then typed 1604 AS cases and 1020 controls for 13 tagging SNPs; 6 showed at least nominal association, 5 of which were in CARD9. We imputed genotypes for 13 additional SNPs but none was more strongly associated with AS than the tagging SNPs. Finally, interrogation of an mRNA expression database revealed that the SNPs most strongly associated with AS (or in strong linkage disequilibrium) were those most associated with CARD9 expression. CARD9 is a plausible candidate for AS given its central role in the innate immune response.

Insights into the Genetic Susceptibility to Type 2 Diabetes from Genome-Wide Association Studies of Obesity-Related Traits

Obesity and type 2 diabetes (T2D) are common and complex metabolic diseases, which are caused by an interchange between environmental and genetic factors. Recently, a number of large-scale genome-wide association studies (GWAS) have improved our knowledge of the genetic architecture and biological mechanisms of these diseases. Currently, more than ~250 genetic loci have been found for monogenic, syndromic, or common forms of T2D and/or obesity-related traits. In this review, we discuss the implications of these GWAS for obesity and T2D, and investigate the overlap of loci for obesity-related traits and T2D, highlighting potential mechanisms that affect T2D susceptibility.

The chromosome 16q region associated with ankylosing spondylitis includes the candidate gene tumour necrosis factor receptor type 1-associated death domain (TRADD)

Objective To replicate and refine the reported association of ankylosing spondylitis (AS) with two non-synonymous single nucleotide polymorphisms (nsSNPs) on chromosome 16q22.1. Methods Firstly, 730 independent UK patients with AS were genotyped for rs9939768 and rs6979 and allele frequencies were compared with 2879 previously typed historic disease controls. Secondly, the two data sets were combined in meta-analyses. Finally, 5 tagging SNPs, located between rs9939768 and rs6979, were analysed in 1604 cases and 1020 controls. Results The association of rs6979 with AS was replicated, p=0.03, OR=1.14 (95% CI 1.01 to 1.28), and a trend for association with rs9939768 detected, p=0.06, OR=1.25 (95% CI 0.99 to 1.57). Meta-analyses revealed association of both SNPs with AS, p=0.0008, OR=1.31 (95% CI 1.12 to 1.54) and p=0.0009, OR=1.15 (95% CI 1.06 to 1.23) for rs9939768 and rs6979, respectively. New associations with rs9033 and rs868213 (p=0.00002, OR=1.23 (95% CI 1.12 to 1.36) and p=0.00002 OR=1.45 (95% CI 1.22 to 1.72), respectively, were identified. Conclusions The region on chromosome 16 that has been replicated in the present work is interesting as the highly plausible candidate gene, tumour necrosis factor receptor type 1 (TNFR1)-associated death domain (TRADD), is located between rs9033 and rs868213. It will require additional work to identify the primary genetic association(s) with AS.

The histone demethylase JARID1A is associated with susceptibility to ankylosing spondylitis.

Associations with disease identified by genome-wide association studies (GWAS) must be replicated and refined to validate causative variants. In the Wellcome Trust Case Control Consortium (WTCCC) GWAS using 14 500 non-synonymous single nucleotide polymorphisms (nsSNP), rs11062385 (a nsSNP in JARID1A) showed nominal association with ankylosing spondylitis (AS) (P=0.0006, odds ratio (OR)=1.26, 95% confidence interval (95% CI)=1.1–1.4). To replicate and refine the association of JARID1A, rs11062385 was genotyped in 730 further cases and compared with allele frequencies in non-AS disease cohorts typed by WTCCC. We replicated the initial association (P=0.04, OR=1.16, 95% CI=1.01–1.34) and identified a strengthened association with AS in a meta-analysis of this new study combined with the original WTCCC study (P=0.0001, OR=1.21, 95% CI=1.10–1.33). We also genotyped nine further intronic tagging SNPs in JARID1A in 1604 AS cases and 1020 new control samples, but none was associated with AS. JARID1A or a locus in strong linkage disequilibrium with it is a positional candidate for susceptibility to AS.

Ankylosing spondylitis is associated with the anthrax toxin receptor 2 gene (ANTXR2)

Objectives ANTXR2 variants have been associated with ankylosing spondylitis (AS) in two previous genome-wide association studies (GWAS) (p∼9×10−8). However, a genome-wide significant association (p<5×10−8) was not observed. We conducted a more comprehensive analysis of ANTXR2 in an independent UK sample to confirm and refine this association. Methods A replication study was carried out with 2978 cases and 8365 controls. Then, these were combined with non-overlapping samples from the two previous GWAS in a meta-analysis. Human leukocyte antigen (HLA)-B27 stratification was also performed to test for ANTXR2-HLA-B27 interaction. Results Out of nine single nucleotide polymorphisms (SNP) in the study, five SNPs were nominally associated (p<0.05) with AS in the replication dataset. In the meta-analysis, eight SNPs showed evidence of association, the strongest being with rs12504282 (OR=0.88, p=6.7×10−9). Seven of these SNPs showed evidence for association in the HLA-B27-positive subgroup, but none was associated with HLA-B27-negative AS. However, no statistically significant interaction was detected between HLA-B27 and ANTXR2 variants. Conclusions ANTXR2 variants are clearly associated with AS. The top SNPs from two previous GWAS (rs4333130 and rs4389526) and this study (rs12504282) are in strong linkage disequilibrium (r2≥0.76). All are located near a putative regulatory region. Further studies are required to clarify the role played by these ANTXR2 variants in AS.

A common functional variant of endoplasmic reticulum aminopeptidase 2 (ERAP2) that reduces major histocompatibility complex class I expression is not associated with ankylosing spondylitis

Sir, The strong genetic association between AS and HLA-B27 has defied explanation for nearly 40 years. However, the additional discovery of a strong association between AS and ERAP1 (endoplasmic reticulum aminopeptidase 1), a gene that almost certainly operates in the trimming of peptides for optimal binding to MHC class I molecules, has rekindled hopes of rapid advances in this field [1]. It has been suggested that another aminopeptidase, endoplasmic reticulum aminopeptidase 2 (ERAP2), may act in concert with ERAP1, trimming residues inefficiently removed by ERAP1 [2]. The association described recently between ERAP2 and Crohn’s disease, which shares many clinical and genetic overlaps with AS [3], suggests that ERAP2 is worthy of further study in AS. An experiment of nature allows us to do this relatively simply. ERAP2 has evolved under balancing selection, similar to the MHC, and includes a high-frequency variant (∼50%) that influences antigen presentation [4]. A single nucleotide polymorphism (SNP), rs2248374 (A to G), located within the 5′ canonical splice site of exon 10, results in an alternatively spliced ERAP2 mRNA that is degraded by nonsense-mediated decay (NMD). Homozygosity for the minor G allele (carried by ∼25% of the population) results in failure to express ERAP2 protein; in turn this genotype is also associated with reduced surface MHC Class I expression on human B cells [4]. Such a dramatic phenotypic effect provides a great opportunity for studying the potential role of ERAP2 in AS, analogous to using a knockout model. We therefore specifically tested for differences in frequency of rs2248374 in 470 sporadic AS cases and 420 healthy, ethnically matched blood donors to determine whether this functionally important variant is implicated in AS. All patients were Caucasians, of UK origin and fulfilled the modified New York criteria for AS. Eighty-four per cent of the AS patients were HLA-B27 positive. All patients gave informed consent and ethical approval was obtained (Multicentre Research Ethics Committee 98/5/23). Genotyping involved restriction fragment length polymorphism (RFLP) analysis of PCR products performed under standard conditions. Primers (forward 5′-GCATCCATGGCTAATGTGCR and reverse 5′-GTTGTGGGAAAGCCGAACTA) amplified a 370-bp product that, in the presence of the G allele, was digested to 214 and 156 bp products by HphI. Genotype and allele frequencies in AS cases and controls were compared using the Cochrane–Armitage test of trend and the chi-squared test, respectively. This study had 80% power under a log-additive model to detect an odds ratio (OR) of 1.3 with a population risk of AS of 0.04% and a risk allele frequency of 46% (frequency in our controls) at a significance level of 0.05. There were no significant differences between allele or genotype frequencies in AS cases and controls at rs2248374 (Table 1). Additionally, due to the potential functional significance of homozygosity at the minor G allele (i.e. no functional ERAP2 protein), we also compared the frequency of GG homozygotes in AS cases and controls, but these proved almost identical again (Table 1). Table 1 Genotypes, minor allele frequency (MAF), OR, 95% CI and P-value for the ERAP2 SNP rs2248374 The association between ERAP1 and AS is the second only to HLA-B27, but formal testing for additional association with the neighbouring gene, ERAP2, had not been previously undertaken. A study of familial AS that typed one marker within ERAP2 reported that an ERAP1/ERAP2 haplotype was over-represented in affected family members [5]. However, such studies are frequently confounded by linkage disequilibrium effects. If genetic variation at ERAP2 were to play a role in susceptibility to AS, then a dramatic loss-of-function variant, such as rs2248374, would appear to be a prime candidate for the ‘causal variant’ at this locus. However, we have now excluded any significant effect arising from this important functional ERAP2 variant. Nonetheless, the possibility of other associations between ERAP2 and AS has not been excluded by this study. This could be done only by systematic mapping of the gene with tagging SNPs in a very large sample after controlling for the significant linkage disequilibrium between ERAP1 and ERAP2 would this be achievable. Our study was insufficiently powered to do this.

Sexual dimorphisms in genetic loci linked to body fat distribution.

Obesity is a chronic condition associated with increased morbidity and mortality and is a risk factor for a number of other diseases including type 2 diabetes and cardiovascular disease. Obesity confers an enormous, costly burden on both individuals and public health more broadly. Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes. Body fat distribution is distinct from overall obesity in measurement, but studies of body fat distribution can yield insights into the risk factors for and causes of overall obesity. Sexual dimorphism in body fat distribution is present throughout life. Though sexual dimorphism is subtle in early stages of life, it is attenuated in puberty and during menopause. This phenomenon could be, at least in part, due to the influence of sex hormones on the trait. Findings from recent large genome-wide association studies (GWAS) for various measures of body fat distribution (including waist-to-hip ratio, hip or waist circumference, trunk fat percentage and the ratio of android and gynoid fat percentage) emphasize the strong sexual dimorphism in the genetic regulation of fat distribution traits. Importantly, sexual dimorphism is not observed for overall obesity (as assessed by body mass index or total fat percentage). Notably, the genetic loci associated with body fat distribution, which show sexual dimorphism, are located near genes that are expressed in adipose tissues and/or adipose cells. Considering the epidemiological and genetic evidence, sexual dimorphism is a prominent feature of body fat distribution. Research that specifically focuses on sexual dimorphism in fat distribution can provide novel insights into human physiology and into the development of obesity and its comorbidities, as well as yield biological clues that will aid in the improvement of disease prevention and treatment.