Meta-analysis of genome-wide association for migraine in six population-based European cohorts

Migraine is a common neurological disorder with a genetically complex background. This paper describes a meta-analysis of genome-wide association (GWA) studies on migraine, performed by the Dutch-Icelandic migraine genetics (DICE) consortium, which brings together six population-based European migraine cohorts with a total sample size of 10,980 individuals (2446 cases and 8534 controls). A total of 32 SNPs showed marginal evidence for association at a P-value<10(-5). The best result was obtained for SNP rs9908234, which had a P-value of 8.00 × 10(-8). This top SNP is located in the nerve growth factor receptor (NGFR) gene. However, this SNP did not replicate in three cohorts from the Netherlands and Australia. Of the other 31 SNPs, 18 SNPs were tested in two replication cohorts, but none replicated. In addition, we explored previously identified candidate genes in the meta-analysis data set. This revealed a modest gene-based significant association between migraine and the metadherin (MTDH) gene, previously identified in the first clinic-based GWA study (GWAS) for migraine (Bonferroni-corrected gene-based P-value=0.026). This finding is consistent with the involvement of the glutamate pathway in migraine. Additional research is necessary to further confirm the involvement of glutamate.     

Genome-wide association study of coronary artery disease in the Japanese

A new understanding of the genetic basis of coronary artery disease (CAD) has recently emerged from genome-wide association (GWA) studies of common single-nucleotide polymorphisms (SNPs), thus far performed mostly in European-descent populations. To identify novel susceptibility gene variants for CAD and confirm those previously identified mostly in populations of European descent, a multistage GWA study was performed in the Japanese. In the discovery phase, we first genotyped 806 cases and 1337 controls with 451 382 SNP markers and subsequently assessed 34 selected SNPs with direct genotyping (541 additional cases) and in silico comparison (964 healthy controls). In the replication phase, involving 3052 cases and 6335 controls, 12 SNPs were tested; CAD association was replicated and/or verified for 4 (of 12) SNPs from 3 loci: near BRAP and ALDH2 on 12q24 (P=1.6 × 10⁻³⁴), HLA-DQB1 on 6p21 (P=4.7 × 10⁻⁷), and CDKN2A/B on 9p21 (P=6.1 × 10⁻¹⁶). On 12q24, we identified the strongest association signal with the strength of association substantially pronounced for a subgroup of myocardial infarction cases (P=1.4 × 10⁻⁴⁰). On 6p21, an HLA allele, DQB1^*0604, could show one of the most prominent association signals in an ∼8-Mb interval that encompasses the LTA gene, where an association with myocardial infarction had been reported in another Japanese study. CAD association was also identified at CDKN2A/B, as previously reported in different populations of European descent and Asians. Thus, three loci confirmed in the Japanese GWA study highlight the likely presence of risk alleles with two types of genetic effects – population specific and common – on susceptibility to CAD.     

Screening for replication of genome-wide SNP associations in sporadic ALS

We recently reported a joint analysis of genome-wide association (GWA) data on 958 sporadic amyotrophic lateral sclerosis (ALS) cases and 932 controls from Ireland and the publicly available data sets from the United States and the Netherlands. The strongest pooled association was rs10260404 in the dipeptidyl-peptidase 6 (DPP6) gene. Here, we sought confirmation of joint analysis signals in both an expanded Irish and a Polish ALS cohort. Among 287 522 autosomal single-nucleotide polymorphisms (SNPs), 27 were commonly associated on joint analysis of the Irish, US and Dutch GWAs. These 27 SNPs were genotyped in an expanded Irish cohort (312 patients with SALS; 259 controls) and an additional Polish cohort (218 patients; 356 controls). Eleven SNPs, including rs10260404, reached a final P-value below 0.05 in the Irish cohort. In the Polish cohort, only one SNP, rs6299711, showed nominal association with ALS. Pooling of data for 1267 patients with ALS and 1336 control subjects did not identify any association reaching Bonferroni significance (P<1.74 x 10(-7)). The present strategy did not reveal any consistently associated SNP across four populations. The result for DPP6 is surprising, as it has been replicated elsewhere. We discuss the possible interpretations and implications of these findings for future ALS GWA studies both within and between populations.     

Identifying genetic risk variants for coronary heart disease in familial hypercholesterolemia: an extreme genetics approach

Mutations in the low-density lipoprotein receptor (LDLR) gene cause familial hypercholesterolemia (FH), a disorder characterized by coronary heart disease (CHD) at young age. We aimed to apply an extreme sampling method to enhance the statistical power to identify novel genetic risk variants for CHD in individuals with FH. We selected cases and controls with an extreme contrast in CHD risk from 17 000 FH patients from the Netherlands, whose functional LDLR mutation was unequivocally established. The genome-wide association (GWA) study was performed on 249 very young FH cases with CHD and 217 old FH controls without CHD (above 65 years for males and 70 years of age for females) using the Illumina HumanHap550K chip. In the next stage, two independent samples (one from the Netherlands and one from Italy, Norway, Spain, and the United Kingdom) of FH patients were used as replication samples. In the initial GWA analysis, we identified 29 independent single nucleotide polymorphisms (SNPs) with suggestive associations with premature CHD (P<1 × 10⁻⁴). We examined the association of these SNPs with CHD risk in the replication samples. After Bonferroni correction, none of the SNPs either replicated or reached genome-wide significance after combining the discovery and replication samples. Therefore, we conclude that the genetics of CHD risk in FH is complex and even applying an ‘extreme genetics'' approach we did not identify new genetic risk variants. Most likely, this method is not as effective in leveraging effect size as anticipated, and may, therefore, not lead to significant gains in statistical power.     

Genome-wide association study identifies novel alleles associated with risk of cutaneous basal cell carcinoma and squamous cell carcinoma

We conducted a genome-wide association study on cutaneous basal cell carcinoma (BCC) among 2045 cases and 6013 controls of European ancestry, with follow-up replication in 1426 cases and 4845 controls. A non-synonymous SNP in the MC1R gene (rs1805007 encoding Arg151Cys substitution), a previously well-documented pigmentation gene, showed the strongest association with BCC risk in the discovery set (rs1805007[T]: OR (95% CI) for combined discovery set and replication set [1.55 (1.45-1.66); P= 4.3 × 10(-17)]. We identified that an SNP rs12210050 at 6p25 near the EXOC2 gene was associated with an increased risk of BCC [rs12210050[T]: combined OR (95% CI), 1.24 (1.17-1.31); P= 9.9 × 10(-10)]. In the locus on 13q32 near the UBAC2 gene encoding ubiquitin-associated domain-containing protein 2, we also identified a variant conferring susceptibility to BCC [rs7335046 [G]; combined OR (95% CI), 1.26 (1.18-1.34); P= 2.9 × 10(-8)]. We further evaluated the associations of these two novel SNPs (rs12210050 and rs7335046) with squamous cell carcinoma (SCC) risk as well as melanoma risk. We found that both variants, rs12210050[T] [OR (95% CI), 1.35 (1.16-1.57); P= 7.6 × 10(-5)] and rs7335046 [G] [OR (95% CI), 1.21 (1.02-1.44); P= 0.03], were associated with an increased risk of SCC. These two variants were not associated with melanoma risk. We conclude that 6p25 and 13q32 are novel loci conferring susceptibility to non-melanoma skin cancer.     

Variation within the CLEC16A gene shows consistent disease association with both multiple sclerosis and type 1 diabetes in Sardinia

Variation within intron 19 of the CLEC16A (KIAA0350) gene region was recently found to be unequivocally associated with type 1 diabetes (T1D) in genome-wide association (GWA) studies in Northern European populations. A variant in intron 22 that is nearly independent of the intron 19 variant showed suggestive evidence of association with multiple sclerosis (MS). Here, we genotyped the rs725613 polymorphism, representative of the earlier reported associations with T1D within CLEC16A, in 1037 T1D cases, 1498 MS cases and 1706 matched controls, all from the founder, autoimmunity-prone Sardinian population. In these Sardinian samples, allele A of rs725613 is positively associated not only with T1D (odds ratio=1.15, P one-tail=5.1 x 10(-3)) but also, and with a comparable effect size, with MS (odds ratio=1.21, P one-tail 6.7 x 10(-5)). Taken together these data provide evidence of joint disease association in T1D and MS within CLEC16A and underline a shared disease pathway.     

A genome-wide association study identifies a potential novel gene locus for keratoconus, one of the commonest causes for corneal transplantation in developed countries

Keratoconus is a condition in which the cornea progressively thins over time, and is a major cause for cornea transplantation. To identify keratoconus susceptibility regions, we performed a comprehensive genome-wide association study (GWAS) using a discovery and replication design. A discovery panel of 222 keratoconus Caucasian patients and 3324 Caucasian controls was genotyped using Illumina 370K beadchips. Further associated and fine-mapping single nucleotide polymorphisms (SNPs) (n= 4905) were genotyped in an independent replication case-control panel of 304 cases and 518 controls and a family panel of 307 subjects in 70 families. Logistic regression models implemented in PLINK were performed to test associations in case-control samples with and without principal component (PC) adjustments. Generalized estimation equation models accounting for familial correlations implemented in GWAF were used for association testing in families. No genome-wide associations were identified in the discovery GWAS panel. From the initial testing without adjustments for PCs, the top three SNPs located at 3p26 (rs6442925), 2q21.3 (rs4954218) and 19q13.3 (rs1428642) were identified with unadjusted P-values of 6.5 × 10(-8), 2.4 × 10(-7) and 3.1 × 10(-7), respectively. After adjustments for PCs, rs1428642 became the most significant through the genome with a P-value of 1.4 × 10(-6), while rs6442925 and rs4954218 were less significant (P= 1.9 × 10(-5) and 2.6 × 10(-4)). SNP rs4954218 was confirmed in two independent replication panels with P-values of 0.004 and 0.009, respectively. Meta-analysis revealed a highest association at rs4954218 with adjusted P= 1.6 × 10(-7) (unadjusted P= 1.2 × 10(-9)). These findings suggest SNP rs4954218, located near the RAB3GAP1 gene, previously reported to be associated with corneal malformation, is a potential susceptibility locus for keratoconus.     

Genetic analysis of autoimmune regulator haplotypes in alopecia areata

Alopecia areata is an immune-mediated disorder, occurring with the highest observed frequency in the rare recessive autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED) syndrome caused by mutations of the autoimmune regulator (AIRE) gene on chromosome 21q22.3. We have previously detected association between alopecia areata and a single nucleotide polymorphism (SNP) in the AIRE gene in patients without APECED, and we now report the findings of an extended examination of the association of alopecia areata with haplotype analysis including six SNPs in the AIRE gene: C-103T, C4144G, T5238C, G6528A, T7215C and T11787C. In Caucasian groups of 295 patients and 363 controls, we found strong association between the AIRE 7215C allele and AA [P = 3.8 x 10(-8), OR (95% CI): 2.69 (1.8-4.0)]. The previously reported association between AA and the AIRE 4144G allele was no longer significant on correction for multiple testing. The AIRE haplotypes CCTGCT and CGTGCC showed a highly significant association with AA [P = 6.05 x 10(-6), 9.47 (2.91-30.8) and P = 0.001, 3.51 (1.55-7.95), respectively]. To select the haplotypes most informative for analysis, we tagged the polymorphisms using SNPTag software. Employing AIRE C-103T, G6528A, T7215C and T11787C as tag SNPs, two haplotypes were associated with AA; AIRE CGCT and AIRE CGCC [P = 3.84 x 10(-7), 11.40 (3.53-36.9) and P = 3.94 x 10(-4), 2.13 (1.39-3.24) respectively]. The AIRE risk haplotypes identified in this study potentially account for a major component of the genetic risk of developing alopecia areata.     

A genome-wide association study identifies novel susceptibility genetic variation for thyrotoxic hypokalemic periodic paralysis

Several lines of evidence have pointed out that genetic components have roles in thyrotoxic hypokalemic periodic paralysis (TTPP). In this study, for the first time we performed genome-wide association study (GWAS) in male hyperthyroid subjects in order to identify genetic loci conferring susceptibility to TTPP. We genotyped 78 Thai male TTPP cases and 74 Thai male hyperthyroid patients without hypokalemia as controls with Illumina Human-Hap610 Genotyping BeadChip. Among the SNPs analyzed in the GWAS, rs312729 at chromosome 17q revealed the lowest P-value for association (P=2.09 × 10(-7)). After fine mapping for linkage disequilibrium blocks surrounding the landmark SNP, we found a significant association of rs623011; located at 75?kb downstream of KCNJ2 on chromosome 17q, reached the GWAS significance after Bonferroni's adjustment (P=3.23 × 10(-8), odds ratio (OR)=6.72; 95% confidence interval (CI)=3.11-14.5). The result was confirmed in an independent cohort of samples consisting of 28 TTPP patients and 48 controls using the same clinical criteria diagnosis (replication analysis P=3.44 × 10(-5), OR=5.13; 95% CI=1.87-14.1; combined-analysis P=3.71 × 10(-12), OR=5.47; 95% CI=3.04-9.83).     

Genome-wide association study confirms extant PD risk loci among the Dutch

In view of the population-specific heterogeneity in reported genetic risk factors for Parkinson's disease (PD), we conducted a genome-wide association study (GWAS) in a large sample of PD cases and controls from the Netherlands. After quality control (QC), a total of 514,799 SNPs genotyped in 772 PD cases and 2024 controls were included in our analyses. Direct replication of SNPs within SNCA and BST1 confirmed these two genes to be associated with PD in the Netherlands (SNCA, rs2736990: P = 1.63 × 10(-5), OR = 1.325 and BST1, rs12502586: P = 1.63 × 10(-3), OR = 1.337). Within SNCA, two independent signals in two different linkage disequilibrium (LD) blocks in the 3' and 5' ends of the gene were detected. Besides, post-hoc analysis confirmed GAK/DGKQ, HLA and MAPT as PD risk loci among the Dutch (GAK/DGKQ, rs2242235: P = 1.22 × 10(-4), OR = 1.51; HLA, rs4248166: P = 4.39 × 10(-5), OR = 1.36; and MAPT, rs3785880: P = 1.9 × 10(-3), OR = 1.19).     

“Replicated” genome wide association for dependence on illegal substances: Genomic regions identified by overlapping clusters of nominally positive SNPs

Declaring “replication” from results of genome wide association (GWA) studies is straightforward when major gene effects provide genome‐wide significance for association of the same allele of the same SNP in each of multiple independent samples. However, such unambiguous replication may be unlikely when phenotypes display polygenic genetic architecture, allelic heterogeneity, locus heterogeneity, and when different samples display linkage disequilibria with different fine structures. We seek chromosomal regions that are tagged by clustered SNPs that display nominally significant association in each of several independent samples. This approach provides one “nontemplate” approach to identifying overall replication of groups of GWA results in the face of difficult genetic architectures. We apply this strategy to 1 million (1M) SNP Affymetrix and Illumina GWA results for dependence on illegal substances. This approach provides high confidence in rejecting the null hypothesis that chance alone accounts for the extent to which clustered, nominally significant SNPs from samples of the same racial/ethnic background identify the same chromosomal regions. There is more modest confidence in: (a) identification of individual chromosomal regions and genes and (b) overlap between results from samples of different racial/ethnic backgrounds. The strong overlap identified among the samples with similar racial/ethnic backgrounds, together with prior work that identified overlapping results in samples of different racial/ethnic backgrounds, support contributions to individual differences in vulnerability to addictions that come from both relatively older allelic variants that are common in many current human populations and newer allelic variants that are common in fewer current human populations. © 2010 Wiley‐Liss, Inc.     

Sequencing of the ANKYRIN 3 gene (ANK3) encoding ankyrin G in bipolar disorder reveals a non-conservative amino acid change in a short isoform of ankyrin G

Significant association between polymorphisms at the ANK3 gene with bipolar disorder has previously been reported and confirmed in several samples. Here we report on association between ANK3 and bipolar disorder in a new sample of 593 patients and 642 controls (UCL2) as well as the results of sequencing of the exons and flanking regions of ANK3 from bipolar patients. Single nucleotide polymorphisms (SNPs) associated with bipolar disorder in our original GWA study (UCL1) were genotyped and tested for association in the new sample. Novel SNPs found by sequencing were genotyped in both samples to test for association with bipolar disorder. None of the SNPs previously associated with bipolar disorder were associated in the UCL2 sample. One of the four SNPs associated in the UCL1 sample, rs1938526, was still significantly associated with bipolar disorder when the UCL1 and UCL2 samples were combined (P?=?0.0095). The results demonstrate the impact of heterogeneity on replication of allelic associations even within well-defined ancestral populations. DNA sequencing revealed a novel low frequency (0.007) ANK3 SNP (ss469104599) which causes a non-conservative amino acid change at position 794 in the shorter isoforms of the ankyrin G protein. Protein-function analysis software predicted the amino acid change to be "probably damaging" and it could therefore be detrimental to the function of this isoform. Given that there was only a modest increase in the allele frequency of ss469104599 in cases compared to controls further association studies are needed in additional samples to establish a possible etiological role for this amino acid change.     

Test of association between 10 single nucleotide polymorphisms in the oxytocin receptor gene and conduct disorder

Animal and human studies have implicated oxytocin in affiliative and prosocial behaviors. We tested whether genetic variation in the oxytocin receptor (OXTR) gene is associated with conduct disorder (CD). Utilizing a family-based sample of adolescent probands recruited from an adolescent substance abuse treatment program, control probands and their families (total sample, n=1750), we conducted three tests of association with CD and 10 single nucleotide polymorphisms (SNPs) in the OXTR gene: (a) family-based comparison utilizing the entire sample; (b) within-Whites, case-control comparison of adolescent patients with CD and controls without CD; and (c) within-Whites case-control comparison of parents of patients and parents of controls. Family-based association tests failed to show significant results (no results P<0.05). While strictly correcting for the number of tests (α=0.002), adolescent patients with CD did not differ significantly from adolescent controls in genotype frequency for the OXTR SNPs tested; similarly, comparison of OXTR genotype frequencies for parents failed to differentiate patient and control family type, except a trend association for rs237889 (P=0.004). We concluded that in this sample, 10 SNPs in the OXTR gene were not significantly associated with CD.     

Population-specific coding variant underlies genome-wide association with adiponectin level

Adiponectin is a protein hormone that can affect major metabolic processes including glucose regulation and fat metabolism. Our previous genome-wide association (GWA) study of circulating plasma adiponectin levels in Filipino women from the Cebu Longitudinal Health and Nutrition Survey (CLHNS) detected a 100 kb two-SNP haplotype at KNG1-ADIPOQ associated with reduced adiponectin (frequency = 0.050, P = 1.8 × 10(-25)). Subsequent genotyping of CLHNS young adult offspring detected an uncommon variant [minor allele frequency (MAF) = 0.025] located ~800 kb from ADIPOQ that showed strong association with lower adiponectin levels (P = 2.7 × 10(-15), n = 1695) and tagged a subset of KNG1-ADIPOQ haplotype carriers with even lower adiponectin levels. Sequencing of the ADIPOQ-coding region detected variant R221S (MAF = 0.015, P = 2.9 × 10(-69)), which explained 17.1% of the variance in adiponectin levels and largely accounted for the initial GWA signal in Filipinos. R221S was not present in 12 514 Europeans with previously sequenced exons. To explore the mechanism of this substitution, we re-measured adiponectin level in 20 R221S offspring carriers and 20 non-carriers using two alternative antibodies and determined that the presence of R221S resulted in artificially low quantification of adiponectin level using the original immunoassay. These data provide an example of an uncommon variant responsible for a GWA signal and demonstrate that genetic associations with phenotypes measured by antibody-based quantification methods can be affected by uncommon coding SNPs residing in the antibody target region.     

Dense mapping of IL18 shows no association in SLE

Systemic lupus erythematosus (SLE) is an autoimmune disease which behaves as a complex genetic trait. At least 20 SLE risk susceptibility loci have been mapped using both candidate gene and genome-wide association strategies. The gene encoding the pro-inflammatory cytokine, IL18, has been reported as a candidate gene showing an association with SLE. This pleiotropic cytokine is expressed in a range of immune cells and has been shown to induce interferon-γ and tumour necrosis factor-α. Serum interleukin-18 has been reported to be elevated in patients with SLE. Here we aimed to densely map single nucleotide polymorphisms (SNPs) across IL18 to investigate the association across this locus. We genotyped 36 across IL18 by Illumina bead express in 372 UK SLE trios. We also genotyped these SNPs in a further 508 non-trio UK cases and were able to accurately impute a dense marker set across IL18 in WTCCC2 controls with a total of 258 SNPs. To improve the study's power, we also imputed a total of 158 SNPs across the IL18 locus using data from an SLE genome-wide association study and performed association testing. In total, we analysed 1818 cases and 10 770 controls in this study. Our large well-powered study (98% to detect odds ratio = 1.5, with respect to rs360719) showed that no individual SNP or haplotype was associated with SLE in any of the cohorts studied. We conclude that we were unable to replicate the SLE association with rs360719 located upstream of IL18. No evidence for association with any other common variant at IL18 with SLE was found.     

High-density association study and nomination of susceptibility genes for hypertension in the Japanese National Project

Essential hypertension is one of the most common, complex diseases, of which considerable efforts have been made to unravel the pathophysiological mechanisms. Over the last decade, multiple genome-wide linkage analyses have been conducted using 300-900 microsatellite markers but no single study has yielded definitive evidence for 'principal' hypertension susceptibility gene(s). Here, we performed a three-tiered, high-density association study of hypertension, which has been recently made possible. For tier 1, we genotyped 80 795 SNPs distributed throughout the genome in 188 male hypertensive subjects and two general population control groups (752 subjects per group). For tier 2 (752 hypertensive and 752 normotensive subjects), we genotyped a panel of 2676 SNPs selected (odds ratio >or= 1.4 and P 相似文献   

Haplotype association of IL-8 gene with Behcet's disease

Interleukin-8 (IL-8), a CXC chemokine that recruits and activates inflammatory cells, plays a critical role in the pathogenesis of Behcet's disease (BD). To investigate the association of the genetic polymorphism of IL-8 and BD, we genotyped IL-8 -845 T/C, -738 T/A, -353 A/T, -251 A/T, +293 G/T, +678 T/C and receptors CXCR-1 +2607 G/C and CXCR-2 +785 C/T polymorphisms in 119 Korean patients with BD and 119 age- and sex-matched healthy blood donors. Then, single nucleotide polymorphisms (SNPs) and haplotypes were analyzed between patients and controls. There were no SNPs associated with BD. However, the frequency of haplotype TAT inferred from SNPs, IL-8 -353 A/T, -251 A/T and +678 T/C, was significantly higher in patients with BD than controls (5.9 vs 0.0%, P = 0.0001), as was haplotype ATC (6.7 vs 0.0%, P < 0.0001). The haplotype difference was still valid in human leukocyte antigen-B51-negative subjects. In conclusion, we found a significant difference in the distribution of IL-8 gene haplotypes between patients with BD and healthy controls. These results suggest that the genetic polymorphisms of proinflammatory chemokine IL-8 can contribute to the pathogenesis of BD.