Genome-wide association studies

Genome-wide association studies (GWAS) are a powerful hypothesis-free tool for the dissection of susceptibility to common heritable human diseases, including osteoporosis. To date, more than 2000 loci for common human diseases have been identified by GWAS. Success using the GWAS model depends on genetic risk being determined by shared stretches of DNA carried with different frequencies in cases and controls, inherited from ancient ancestors, termed the “common disease–common variant” hypothesis. Not all disease risk is caused by common variants, however, and thus GWAS will not detect all variants involved. Successful GWAS performance requires careful quality control, especially as the effect sizes under study are modest, and there are multiple potential sources of error. Conservative interpretation, use of stringent significance thresholds, and replication in independent cohorts are required to ensure results are robust. Despite these challenging parameters, much has been learnt from GWAS and, as the approach matures and is modified to identify a wider range of variants, significantly more will be learnt about the etiopathogenesis of common diseases such as osteoporosis.

Genome-wide association study to identify genetic determinants of severe asthma

Background: The genetic basis for developing asthma has been extensively studied. However, association studies to date have mostly focused on mild to moderate disease and genetic risk factors for severe asthma remain unclear. Objective: To identify common genetic variants affecting susceptibility to severe asthma. Methods: A genome-wide association study was undertaken in 933 European ancestry individuals with severe asthma based on Global Initiative for Asthma (GINA) criteria 3 or above and 3346 clean controls. After standard quality control measures, the association of 480 889 genotyped single nucleotide polymorphisms (SNPs) was tested. To improve the resolution of the association signals identified, non-genotyped SNPs were imputed in these regions using a dense reference panel of SNP genotypes from the 1000 Genomes Project. Then replication of SNPs of interest was undertaken in a further 231 cases and 1345 controls and a meta-analysis was performed to combine the results across studies. Results: An association was confirmed in subjects with severe asthma of loci previously identified for association with mild to moderate asthma. The strongest evidence was seen for the ORMDL3/GSDMB locus on chromosome 17q12-21 (rs4794820, p=1.03×10 (-8)following meta-analysis) meeting genome-wide significance. Strong evidence was also found for the IL1RL1/IL18R1 locus on 2q12 (rs9807989, p=5.59×10 (-8) following meta-analysis) just below this threshold. No novel loci for susceptibility to severe asthma met strict criteria for genome-wide significance. Conclusions: The largest genome-wide association study of severe asthma to date was carried out and strong evidence found for the association of two previously identified asthma susceptibility loci in patients with severe disease. A number of novel regions with suggestive evidence were also identified warranting further study.

Dissection of the genetics of Parkinson's disease identifies an additional association 5' of SNCA and multiple associated haplotypes at 17q21

We performed a genome-wide association study (GWAS) in 1705 Parkinson's disease (PD) UK patients and 5175 UK controls, the largest sample size so far for a PD GWAS. Replication was attempted in an additional cohort of 1039 French PD cases and 1984 controls for the 27 regions showing the strongest evidence of association (P < 10 4). We replicated published associations in the 4q22/SNCA and 17q21/MAPT chromosome regions (P < 10 10) and found evidence for an additional independent association in 4q22/SNCA.A detailed analysis of the haplotype structure at 17q21 showed that there are three separate risk groups within this region. We found weak but consistent evidence of association for common variants located in three previously published associated regions (4p15/BST1, 4p16/GAK and 1q32/PARK16). We found no support for the previously reported SNP association in 12q12/LRRK2. We also found an association of the two SNPs in 4q22/SNCA with the age of onset of the disease. © The Author 2010. Published by Oxford University Press.

Genome-wide association study identifies five new schizophrenia loci

We examined the role of common genetic variation in schizophrenia in a genome-wide association study of substantial size: a stage 1 discovery sample of 21,856 individuals of European ancestry and a stage 2 replication sample of 29,839 independent subjects. The combined stage 1 and 2 analysis yielded genome-wide significant associations with schizophrenia for seven loci, five of which are new (1p21.3, 2q32.3, 8p23.2, 8q21.3 and 10q24.32-q24.33) and two of which have been previously implicated (6p21.32-p22.1 and 18q21.2). The strongest new finding (P = 1.6 × 10 -11) was with rs1625579 within an intron of a putative primary transcript for MIR137 (microRNA 137), a known regulator of neuronal development. Four other schizophrenia loci achieving genome-wide significance contain predicted targets of MIR137, suggesting MIR137-mediated dysregulation as a previously unknown etiologic mechanism in schizophrenia. In a joint analysis with a bipolar disorder sample (16,374 affected individuals and 14,044 controls), three loci reached genome-wide significance: CACNA1C (rs4765905, P = 7.0 × 10 -9), ANK3 (rs10994359, P = 2.5 × 10 -8) and the ITIH3-ITIH4 region (rs2239547, P = 7.8 × 10 -9).

Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1

We report a genome-wide association study for open-angle glaucoma (OAG) blindness using a discovery cohort of 590 individuals with severe visual field loss (cases) and 3,956 controls. We identified associated loci at TMCO1 (rs4656461[G] odds ratio (OR) = 1.68, P = 6.1 × 10-10) and CDKN2B-AS1 (rs4977756[A] OR = 1.50, P = 4.7 × 10-9). We replicated these associations in an independent cohort of cases with advanced OAG (rs4656461 P = 0.010; rs4977756 P = 0.042) and two additional cohorts of less severe OAG (rs4656461 combined discovery and replication P = 6.00 × 10-14, OR = 1.51, 95% CI 1.35-1.68; rs4977756 combined P = 1.35 × 10-14, OR = 1.39, 95% CI 1.28-1.51). We show retinal expression of genes at both loci in human ocular tissues. We also show that CDKN2A and CDKN2B are upregulated in the retina of a rat model of glaucoma. © 2011 Nature America, Inc. All rights reserved.

Genome-wide association study of ankylosing spondylitis identifies non-MHC susceptibility loci

To identify susceptibility loci for ankylosing spondylitis, we undertook a genome-wide association study in 2,053 unrelated ankylosing spondylitis cases among people of European descent and 5,140 ethnically matched controls, with replication in an independent cohort of 898 ankylosing spondylitis cases and 1,518 controls. Cases were genotyped with Illumina HumHap370 genotyping chips. In addition to strong association with the major histocompatibility complex (MHC; P 10 800), we found association with SNPs in two gene deserts at 2p15 (rs10865331; combined P = 1.9 × 10 19) and 21q22 (rs2242944; P = 8.3 × 10 20), as well as in the genes ANTXR2 (rs4333130; P = 9.3 × 10 8) and IL1R2 (rs2310173; P = 4.8 × 10 7). We also replicated previously reported associations at IL23R (rs11209026; P = 9.1 × 10 14) and ERAP1 (rs27434; P = 5.3 × 10 12). This study reports four genetic loci associated with ankylosing spondylitis risk and identifies a major role for the interleukin (IL)-23 and IL-1 cytokine pathways in disease susceptibility. © 2010 Nature America, Inc. All rights reserved.

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-B27, ERAP1 and IL23R. A recent genome-wide association screen (GWAS) identified associations (P0.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 P0.0004, odds ratio (OR)1.2, 95% confidence interval (CI)1.1-1.4; rs3812571 P0.0003, OR1.2, 95% CI1.1-1.4). Meta-analysis revealed strong associations of both SNPs with AS, rs4077515 P0.000005, OR1.2, 95% CI1.1-1.3 and rs3812571 P0.000006, OR1.2, 95% CI1.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.

Predominant association of HLA-B*2704 with ankylosing spondylitis in Chinese Han patients

The HLA-B27 subtypes have a varied racial and ethnic prevalence throughout the world. However, the association of B27-subtypes with ankylosing spondylitis (AS) in the mainland China is unknown. To determine the association of B27-subtypes with AS in the Mainland Chinese Han population, a total of unrelated 153 patients with AS were enrolled in a large case-control association study, and 1545 unrelated, healthy, ethnically matched blood donors were included as controls. The genotyping of B27 and its subtypes was performed using the polymerase chain reaction with sequence specific primers (PCR-SSP). A total of 130 (84.97%) AS patients and 61 (3.95%) healthy controls were B27 positive. Three B27-subtypes, B*2704, B*2705 and B*2710, were further identified, of which both B*2704 and B*2705 were strongly AS associated. B*2710 was only detected in one AS patient and two other healthy controls. Considering only B27-positive cases and controls, a statistically different frequency of B27-subtypes was observed, with an over-representation of B*2704 (P = 0.018). B*2704 was clearly more strongly associated than B*2705 with AS [odds ratio (OR) = 2.4, P = 0.011]. Furthermore, a combined analysis including three previous studies of B27-subtype distributions in Chinese AS cases confirmed the stronger association of B*2704 with AS than B*2705 (OR = 2.5, P = 0.00094).

Genome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls

Copy number variants (CNVs) account for a major proportion of human genetic polymorphism and have been predicted to have an important role in genetic susceptibility to common disease. To address this we undertook a large, direct genome-wide study of association between CNVs and eight common human diseases. Using a purpose-designed array we typed 19,000 individuals into distinct copy-number classes at 3,432 polymorphic CNVs, including an estimated 50% of all common CNVs larger than 500 base pairs. We identified several biological artefacts that lead to false-positive associations, including systematic CNV differences between DNAs derived from blood and cell lines. Association testing and follow-up replication analyses confirmed three loci where CNVs were associated with diseaseIRGM for Crohns disease, HLA for Crohns disease, rheumatoid arthritis and type 1 diabetes, and TSPAN8 for type 2 diabetesalthough in each case the locus had previously been identified in single nucleotide polymorphism (SNP)-based studies, reflecting our observation that most common CNVs that are well-typed on our array are well tagged by SNPs and so have been indirectly explored through SNP studies. We conclude that common CNVs that can be typed on existing platforms are unlikely to contribute greatly to the genetic basis of common human diseases. © 2010 Macmillan Publishers Limited. All rights reserved.

Genome-wide association study identifies eight loci associated with blood pressure

Elevated blood pressure is a common, heritable cause of cardiovascular disease worldwide. To date, identification of common genetic variants influencing blood pressure has proven challenging. We tested 2.5 million genotyped and imputed SNPs for association with systolic and diastolic blood pressure in 34,433 subjects of European ancestry from the Global BPgen consortium and followed up findings with direct genotyping (N 71,225 European ancestry, N 12,889 Indian Asian ancestry) and in silico comparison (CHARGE consortium, N = 29,136). We identified association between systolic or diastolic blood pressure and common variants in eight regions near the CYP17A1 (P = 7 × 10 24), CYP1A2 (P = 1 × 10 23), FGF5 (P = 1 × 10 21), SH2B3 (P = 3 × 10 18), MTHFR (P = 2 × 10 13), c10orf107 (P = 1 × 10 9), ZNF652 (P = 5 × 10 9) and PLCD3 (P = 1 × 10 8) genes. All variants associated with continuous blood pressure were associated with dichotomous hypertension. These associations between common variants and blood pressure and hypertension offer mechanistic insights into the regulation of blood pressure and may point to novel targets for interventions to prevent cardiovascular disease.

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 × 10-9, OR ∼1.2) but also with rs11209026 (P < 10-10, OR ∼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. © The Author 2009. Published by Oxford University Press on behalf of the British Society for Rheumatology. All rights reserved.

Investigating the genetic association between ERAP1 and ankylosing spondylitis

A strong association between ERAP1 and ankylosing spondylitis (AS) was recently identified by the Wellcome Trust Case Control Consortium and the Australo-Anglo-American Spondylitis Consortium (WTCCC-TASC) study. ERAP1 is highly polymorphic with strong linkage disequilibrium evident across the gene. We therefore conducted a series of experiments to try to identify the primary genetic association(s) with ERAP1. We replicated the original associations in an independent set of 730 patients and 1021 controls, resequenced ERAP1 to define the full extent of coding polymorphisms and tested all variants in additional association studies. The genetic association with ERAP1 was independently confirmed; the strongest association was with rs30187 in the replication set (P = 3.4 × 103). When the data were combined with the original WTCCC-TASC study the strongest association was with rs27044 (P = 1.1 × 10-9). We identified 33 sequence polymorphisms in ERAP1, including three novel and eight known non-synonymous polymorphisms. We report several new associations between AS and polymorphisms distributed across ERAP1 from the extended case-control study, the most significant of which was with rs27434 (P = 4.7 × 10-7). Regression analysis failed to identify a primary association clearly; we therefore used data from HapMap to impute genotypes for an additional 205 non-coding SNPs located within and adjacent to ERAP1. A number of highly significant associations (P < 5 × 10-9) were identified in regulatory sequences which are good candidates for causing susceptibility to AS, possibly by regulating ERAP1 expression. © 2009 The Author(s).

Genome-wide association study of ulcerative colitis identifies three new susceptibility loci, including the HNF4A region

Ulcerative colitis is a common form of inflammatory bowel disease with a complex etiology. As part of the Wellcome Trust Case Control Consortium 2, we performed a genome-wide association scan for ulcerative colitis in 2,361 cases and 5,417 controls. Loci showing evidence of association at P 1 × 10 5 were followed up by genotyping in an independent set of 2,321 cases and 4,818 controls. We find genome-wide significant evidence of association at three new loci, each containing at least one biologically relevant candidate gene, on chromosomes 20q13 (HNF4A; P = 3.2 × 10 17), 16q22 (CDH1 and CDH3; P = 2.8 × 10 8) and 7q31 (LAMB1; P = 3.0 × 10 8). Of note, CDH1 has recently been associated with susceptibility to colorectal cancer, an established complication of longstanding ulcerative colitis. The new associations suggest that changes in the integrity of the intestinal epithelial barrier may contribute to the pathogenesis of ulcerative colitis. © 2009 Nature America, Inc. All rights reserved.