Genome-Wide Association Studies of Asthma in Population-Based Cohorts Confirm Known and Suggested Loci and Identify an Additional Association near HLA

Rationale: Asthma has substantial morbidity and mortality and a strong genetic component, but identification of genetic risk factors is limited by availability of suitable studies. Objectives: To test if population-based cohorts with self-reported physician-diagnosed asthma and genome-wide association (GWA) data could be used to validate known associations with asthma and identify novel associations. Methods: The APCAT (Analysis in Population-based Cohorts of Asthma Traits) consortium consists of 1,716 individuals with asthma and 16,888 healthy controls from six European-descent population-based cohorts. We examined associations in APCAT of thirteen variants previously reported as genome-wide significant (P<5x10-8) and three variants reported as suggestive (P<5×10-7). We also searched for novel associations in APCAT (Stage 1) and followed-up the most promising variants in 4,035 asthmatics and 11,251 healthy controls (Stage 2). Finally, we conducted the first genome-wide screen for interactions with smoking or hay fever. Main Results: We observed association in the same direction for all thirteen previously reported variants and nominally replicated ten of them. One variant that was previously suggestive, rs11071559 in RORA, now reaches genome-wide significance when combined with our data (P = 2.4×10-9). We also identified two genome-wide significant associations: rs13408661 near IL1RL1/IL18R1 (PStage1+Stage2 = 1.1x10-9), which is correlated with a variant recently shown to be associated with asthma (rs3771180), and rs9268516 in the HLA region (PStage1+Stage2 = 1.1x10-8), which appears to be independent of previously reported associations in this locus. Finally, we found no strong evidence for gene-environment interactions with smoking or hay fever status. Conclusions: Population-based cohorts with simple asthma phenotypes represent a valuable and largely untapped resource for genetic studies of asthma. © 2012 Ramasamy et al.

Genomewide screens in ankylosing spondylitis

Efforts to identify genes other than HLA-B27 in AS have been driven by the strength of the evidence from genetic epidemiology studies indicating that HLA-B27, although a major gene in AS, is clearly not the only significant gene operating. This is the case for both genetic determinants of disease-susceptibility and phenotypic characteristics such as disease severity and associated disease features. In this chapter the genetic epidemiology of AS and the gene-mapping studies performed to date will be reviewed and the future direction of research in this field discussed.

Whole-genome screening in ankylosing spondylitis: Evidence of non-MHC genetic-susceptibility loci

Ankylosing spondylitis (AS) is a common inflammatory arthritis predominantly affecting the axial skeleton. Susceptibility to the disease is thought to be oligogenic. To identify the genes involved, we have performed a genomewide scan in 185 families containing 255 affected sibling pairs. Two-point and multipoint nonparametric linkage analysis was performed. Regions were identified showing "suggestive" or stronger linkage with the disease on chromosomes 1p, 2q, 6p, 9q, 10q, 16q, and 19q. The MHC locus was identified as encoding the greatest component of susceptibility, with an overall LOD score of 15.6. The strongest non-MHC linkage lies on chromosome 16q (overall LOD score 4.7). These results strongly support the presence of non-MHC genetic-susceptibility factors in AS and point to their likely locations.

The interleukin 1 gene cluster contains a major susceptibility locus for ankylosing spondylitis

Ankylosing spondylitis (AS) is a common and highly heritable inflammatory arthropathy. Although the gene HLA-B27 is almost essential for the inheritance of the condition, it alone is not sufficient to explain the pattern of familial recurrence of the disease. We have previously demonstrated suggestive linkage of AS to chromosome 2q13, a region containing the interleukin 1 (IL-1) family gene cluster, which includes several strong candidates for involvement in the disease. In the current study, we describe strong association and transmission of IL-1 family gene cluster single-nucleotide polymorphisms and haplotypes with AS.

Five years of GWAS discovery

The past five years have seen many scientific and biological discoveries made through the experimental design of genome-wide association studies (GWASs). These studies were aimed at detecting variants at genomic loci that are associated with complex traits in the population and, in particular, at detecting associations between common single-nucleotide polymorphisms (SNPs) and common diseases such as heart disease, diabetes, auto-immune diseases, and psychiatric disorders. We start by giving a number of quotes from scientists and journalists about perceived problems with GWASs. We will then briefly give the history of GWASs and focus on the discoveries made through this experimental design, what those discoveries tell us and do not tell us about the genetics and biology of complex traits, and what immediate utility has come out of these studies. Rather than giving an exhaustive review of all reported findings for all diseases and other complex traits, we focus on the results for auto-immune diseases and metabolic diseases. We return to the perceived failure or disappointment about GWASs in the concluding section. © 2012 The American Society of Human Genetics.

Matters Arising: Zhu et al, "A novel gene variation of TNFα associated with ankylosing spondylitis: A reconfirmed study"

In a recently published study1 involving 79 ankylosing spondylitis (AS) patients and 132 unrelated healthy blood donors, Zhu et al report association of a tumour necrosis factor α (TNFα) single nucleotide polymorphism (SNP), rs1799724, with AS.1 I have many concerns with this paper and its conclusions...

Recurrence risk modelling of the genetic susceptibility to ankylosing spondylitis

Objectives - It has long been suspected that susceptibility to ankylosing spondylitis (AS) is influenced by genes lying distant to the major histocompatibility complex. This study compares genetic models of AS to assess the most likely mode of inheritance, using recurrence risk ratios in relatives of affected subjects. Methods - Recurrence risk ratios in different degrees of relatives were determined using published data from studies specifically designed to address the question. The methods of Risch were used to determine the expected recurrence risk ratios in different degrees of relatives, assuming equal first degree relative recurrence risk between models. Goodness of fit was determined by χ2 comparison of the expected number of affected subjects with the observed number, given equal numbers of each type of relative studied. Results - The recurrence risks in different degrees of relatives were: monozygotic (MZ) twins 63% (17/27), first degree relatives 8.2% (441/5390), second degree relatives 1.0% (8/834), and third degree relatives 0.7% (7/997). Parent-child recurrence risk (7.9%, 37/466) was not significantly different from the sibling recurrence risk (8.2%, 404/4924), excluding a significant dominance genetic component to susceptibility. Poor fitting models included single gene, genetic heterogeneity, additive, two locus multiplicative, and one locus and residual polygenes (χ2 > 32 (two degrees of freedom), p < 10-6 for all models). The best fitting model studied was a five locus model with multiplicative interaction between loci (χ2 = 1.4 (two degrees of freedom), p = 0.5). Oligogenic multiplicative models were the best fitting over a range of population prevalences and first degree recurrence risk rates. Conclusions - This study suggests that of the genetic models tested, the most likely model operating in AS is an oligogenic model with predominantly multiplicative interaction between loci.

Association of variants in MMEL1 and CTLA4 with rheumatoid arthritis in the Han Chinese population

Background: The genome-wide association study era has made great progress in identifying susceptibility genes and genetic loci for rheumatoid arthritis (RA) in populations of White European ancestry. However, few studies have tried to dissect disease aetiopathogenesis in other ethnic populations. Objective: To investigate these associations in the Han Chinese population. Methods: Haplotypes from the HapMap database Chinese population were used to select tag-single-nucleotide polymorphisms (SNPs) (r2 =0.8) across 19 distinct RA genomic regions. A two phase case-control association study was performed, with 169 SNPs genotyped in phase I (n=571 cases, n=880 controls), and 64 SNPs achieving p<0.2 in the first phase being genotyped in phase II (n=464 cases, n=822 controls). Association statistics were calculated using permutation tests both unadjusted and adjusted for the number of markers studied. Results: Robust association was detected for MMEL1 and CTLA4 , and modest association was identified for another six loci: PADI4 , STAT4 , PRDM1 , CDK6 , TRAF1-C5 and KIF5A-PIP4K2C. All three markers genotyped in MMEL1 demonstrated association, with peak signal for rs3890745 (p=2.6×10 -5unadjusted, p=0.003 adjusted, OR=0.79). For CTLA4 , significance was detected for three of five variants showing association, with peak association for marker rs12992492 (p=4.3×10-5 unadjusted, p=0.0021 adjusted, OR=0.77). Lack of association of common variants in PTPN22 with RA in Han Chinese was confirmed. Conclusion: This study identifies MMEL1 and CTLA4 as RA susceptibility genes, provides suggestive evidence of association for a further six loci in the Han Chinese population and confirms lack of PTPN22 association in Asian populations. It also confirms the value of multiethnic population studies to help dissect disease aetiopathogenesis.

Association of STAT3 and TNFRSF1A with ankylosing spondylitis in Han Chinese

Objectives: Recent association studies by the Australo-Anglo-American Spondyloarthritis Consortium (TASC) in Caucasian European populations from Australia, North America and the UK have identified a number of genes as being associated with ankylosing spondylitis (AS). A candidate gene study in a Han Chinese population was performed based on these findings to identify associated genes in this population. Methods: A case-control study was performed in a Han Chinese population of patients with AS (n=775) and controls (n=1587) from Shanghai and Nanjing. All patients met the modified New York criteria for AS. The cases and controls were genotyped for 115 single nucleotide polymorphisms (SNPs) tagging IL23R, ERAP1, STAT3, JAK2, TNFRSF1A and TRADD, as well as other confirmation SNPs from the TASC study, using the Sequenom iPlex and the ABI OpenArray platforms. Statistical analysis of genotyped SNPs was performed using the Cochran - Armitage test for trend and meta-analysis was performed using METAL. SNPs in AS-associated genes in this study were then imputed using MaCH, and association with AS tested by logistic regression. Results: SNPs in TNFRSF1A (rs4149577, p=8.2×10-4), STAT3 (rs2293152, p=0.0015; rs1053005, p=0.017) and ERAP1 (rs27038, p=0.0091; rs27037, p=0.0092) were significantly associated with AS in Han Chinese. Association was also observed between AS and the intergenic region 2p15 (rs10865331, p=0.023). The lack of association between AS and IL23R in Han Chinese was confirmed (all SNPs p>0.1). Conclusions: The study results demonstrate for the first time that genetic polymorphisms in STAT3, TNFRSF1A and 2p15 are associated with AS in Han Chinese, suggesting common pathogenic mechanisms for the disease in Chinese and Caucasian European populations. Furthermore, previous findings demonstrating that ERAP1, but not IL23R, is associated with AS in Chinese patients were confirmed.

Matters Arising: HLA-DR-DQ haplotypes and genotypes in Finnish patients with rheumatoid arthritis

We read with great interest the paper by Laivoranta-Nyman et al.1 They studied Finnish patients with rheumatoid arthritis (RA) and controls, and suggested that there exist susceptibility, neutral and protective HLA-DR-DQ haplotypes that do not match the predictions of current hypotheses for the mechanism of association of the HLA class II region and RA...

Finnish HLA studies confirm the increased risk conferred by HLA-B27 homozygosity in ankylosing spondylitis

Objective: To determine the influence of HLA-B27 homozygosity and HLA-DRB1 alleles in the susceptibility to, and severity of, ankylosing spondylitis in a Finnish population. Methods: 673 individuals from 261 families with ankylosing spondylitis were genotyped for HLA-DRB1 alleles and HLA-B27 heterozygosity/ homozygosity. The frequencies of HLA-B27 homozygotes in probands from these families were compared with the expected number of HLA-B27 homozygotes in controls under Hardy-Weinberg equilibrium (HWE). The effect of HLA-DRB1 alleles was assessed using a logistic regression procedure conditioned on HLA-B27 and case-control analysis. Results: HLA-B27 was detected in 93% of cases of ankylosing spondylitis. An overrepresentation of HLA-B27 homozygotes was noted in ankylosing spondylitis (11%) compared with the expected number of HLA-B27 homozygotes under HWE (4%) (odds ratio (OR) = 3.3 (95% confidence interval, 1.6 to 6.8), p = 0.002). HLA-B27 homozygosity was marginally associated with reduced BASDAI (HLA-B27 homozygotes, 4.5 (1.6); HLA-B27 heterozygotes, 5.4 (1.8) (mean (SD)), p = 0.05). Acute anterior uveitis (AAU) was present in significantly more HLA-B27 positive cases (50%) than HLA-B27 negative cases (16%) (OR = 5.4 (1.7 to 17), p<0.004). HLA-B27 positive cases had a lower average age of symptom onset (26.7 (8.0) years) compared with HLA-B27 negative cases (35.7 (11.2) years) (p<0.0001). Conclusions: HLA-627 homozygosity is associated with a moderately increased risk of ankylosing spondylitis compared with HLA-β27 heterozygosity. HLA-B27 positive cases had an earlier age of onset of ankylosing spondylitis than HLA-B27 negative cases and were more likely to develop AAU. HLA-DRB1 alleles may influence the age of symptom onset of ankylosing spondylitis.

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.

The F158V polymorphism in FcγRIIIA shows disparate associations with rheumatoid arthritis in two genetically distinct populations

Objectives: To investigate the association of the FcγRIIIA gene with rheumatoid orthritis (RA) in two genetically distinct groups: a white group from the United Kingdom and a northern Indian group. Methods: The distributions of the two alleles of the FcγRIIIA F158V polymorphism were determined in 398 white patients from the United Kingdom and 63 Indian patients with RA and compared with those from 289 United Kingdom and 93 Indian healthy controls, respectively. Results: Among the Indian patients, the frequency of the rare 158V allele and the proportion of 158VV homozygotes were reduced (relative risk (RR)=0.3, 95% confidence interval (95% CI) 0.1 to 1.1, p<0.06), reaching statistical significance for carrying the 158VV phenotype relative to 158FV or FF (RR=0.2, 95% CI 0.05-0.9, p<0.02). Conversely, no significant deviation in allelic frequencies was noted between the patients and controls from the United Kingdom. Conclusions: The 158VV phenotype showed a weak protective effect against developing RA in the Indian group. However, this sample was small (resulting in a low power for statistical analysis) and no independent confirmation was found in the larger white United Kingdom group. Thus the FcγRIIIA locus is unlikely to be of major importance in causing RA.

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 nonsynonymous 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.

Prospective meta-analysis of interleukin 1 gene complex polymorphisms confirms associations with ankylosing spondylitis

Objectives: The aim of the current study was to determine the contribution of interleukin (IL) 1 gene cluster polymorphisms previously implicated in susceptibility for ankylosing spondylitis (AS) to AS susceptibility in different populations worldwide. Methods: Nine polymorphisms in the IL1 gene cluster members IL1A (rs2856836, rs17561 and rs1894399), IL1B (rs16944), IL1F10 (rs3811058) and IL1RN (rs419598, the IL1RA VNTR, rs315952 and rs315951) were genotyped in 2675 AS cases and 2592 healthy controls recruited in 12 different centres in 10 countries. Association of variants with AS was tested by Mantel-Haenszel random effects analysis. Results: Strong association was observed with three single nucleotide polymorphisms (SNPs) in the IL1A gene (rs2856836, rs17561, rs1894399, p = 0.0036, 0.000019 and 0.0003, respectively). There was no evidence of significant heterogeneity of effects between centres, and no evidence of non-combinability of findings. The population attributable risk fraction of these variants in Caucasians is estimated at 4-6%. Conclusions: This study confirms that IL1A is associated with susceptibility to AS. Association of the other IL1 gene complex members could not be excluded in specific populations. Prospective meta-analysis is a useful tool in confirmation studies of genes associated with complex genetic disorders such as AS, providing sufficiently large sample sizes to produce robust findings often not achieved in smaller individual cohorts.