Skip to main content

Two-stage case-control association study of dopamine-related genes and migraine

Abstract

Background

We previously reported risk haplotypes for two genes related with serotonin and dopamine metabolism: MAOA in migraine without aura and DDC in migraine with aura. Herein we investigate the contribution to migraine susceptibility of eight additional genes involved in dopamine neurotransmission.

Methods

We performed a two-stage case-control association study of 50 tag single nucleotide polymorphisms (SNPs), selected according to genetic coverage parameters. The first analysis consisted of 263 patients and 274 controls and the replication study was composed by 259 cases and 287 controls. All cases were diagnosed according to ICHD-II criteria, were Spanish Caucasian, and were sex-matched with control subjects.

Results

Single-marker analysis of the first population identified nominal associations of five genes with migraine. After applying a false discovery rate correction of 10%, the differences remained significant only for DRD2 (rs2283265) and TH (rs2070762). Multiple-marker analysis identified a five-marker T-C-G-C-G (rs12363125-rs2283265-rs2242592-rs1554929-rs2234689) risk haplotype in DRD2 and a two-marker A-C (rs6356-rs2070762) risk haplotype in TH that remained significant after correction by permutations. These results, however, were not replicated in the second independent cohort.

Conclusion

The present study does not support the involvement of the DRD1, DRD2, DRD3, DRD5, DBH, COMT, SLC6A3 and TH genes in the genetic predisposition to migraine in the Spanish population.

Peer Review reports

Background

Migraine is a highly prevalent neurological disorder involving multiple susceptibility genes and environmental factors [1, 2]. The current clinical classification follows the International Criteria for Headache Disorders (ICHD-II), with the two main categories of migraine without aura (MO) and migraine with aura (MA) [3]. The pathophysiology of migraine is not entirely understood, but a role for dopamine (DA) was already suggested thirty years ago [4]. The DA hypothesis relies on the observed signs of central DA hypersensitivity in migraine patients and the known capacity of DA receptors to regulate nociception, vascular tone and autonomic responses [5]. Studies in animal models revealed that DA receptors are present in the trigeminovascular pathway and showed that DA can act as an inhibitor of nociceptive trigeminovascular transmission in the rat brain [6]. Along this line, DA antagonists have proved useful in aborting migraine headache or associated symptoms [5]. However, DA antagonists are not always selective and may act through DA receptor-independent mechanisms [7]. Also, a review of pharmacological and therapeutic studies in migraine could not provide convincing evidence of a direct role of DA in migraine pathogenesis [8].

Several association studies in different populations have focused on genes encoding proteins of the dopaminergic neurotransmission system, including DA receptors, the DA transporter, and enzymes involved in the synthesis and catabolism of DA. These studies provided conflicting results [7], although a recent, most comprehensive analysis of 10 dopamine-related genes in MA suggested that DBH and SLC6A3, at least, might be involved in migraine pathogenesis [9].

In a previous study that evaluated the contribution of 19 serotonin-related genes to migraine susceptibility in our cohort of Spanish migraineurs, we reported risk haplotypes in MAOA for migraine without aura and in DDC for migraine with aura [10], both genes being key players in the serotonin and dopamine metabolic pathways. In order to further elucidate the involvement of the dopaminergic system in migraine liability, nine dopamine-related genes were selected for a two-stage case-control association study in the Spanish population.

Methods

Subjects

Our initial sample (population 1) was recruited between 2002 and 2006 in Spain and consisted of 271 migraineurs (mean age 37 +/- 16 years) and 285 unrelated migraine-free controls (mean age 55 +/- 18) matched for ethnicity (Caucasian Spanish) and sex frequency (76% women). The follow-up replication study (population 2) consisted of 272 patients and 302 healthy controls recruited subsequently between 2006 and 2007, also in Spain. All patients were diagnosed by clinical neurologists in the team (M.J.S., B.N., S.B. or A.M.) as having MO (55.9% in population 1 and 61.4% in population 2) or MA based on the International Criteria for Headache Disorders 2nd edition (ICHD-II) [3] after administration of a structured questionnaire and direct interview and examination. Patients were recruited from three centers (Hospital Universitari Vall d'Hebron, HUVH, Barcelona; Hospital Sant Joan de Déu, Manresa; Fundación Pública Galega de Medicina Xenómica, FPGMX, Santiago de Compostela). Patients with hemiplegic migraine, a MA variant usually showing monogenic inheritance, were excluded. The control samples consisted of Caucasian Spanish unrelated adult subjects (blood donors, individuals that underwent surgery unrelated to migraine or unaffected partners of migraine patients) that were matched for sex with patients and recruited in the same geographic areas (HUVH and FPGMX). Migraine and positive family history of migraine or any type of severe or recurrent headache in first - degree relatives were excluded in all control subjects through personal interview. Genomic DNA was isolated from peripheral leukocytes or saliva. Research was approved by the local Ethics Committees and all the adult participants, and the children or their parents gave their informed consent prior to participate in the study, according to the Helsinki Declaration.

SNP selection, SNPlex design, genotyping and quality control

We selected SNPs in nine candidate genes involved in dopaminergic neurotransmission. These genes encode five DA receptors (DRD1, DRD2, DRD3, DRD4 and DRD5), three enzymes involved in DA degradation (COMT and DBH) or synthesis (TH), and the DA transporter (SLC6A3 or DAT1; [see Additional File 1]).

In order to ensure a full genetic coverage of these genes and to minimize redundancy, we used the LD-select software [11] and the HapMap database (http://www.hapmap.org; release 20)[12] to evaluate the LD pattern of the region spanning each candidate gene plus three to five kb flanking regions. TagSNPs were selected at an r2 threshold of 0.85 and minimal allele frequency (MAF)>0.15 for genes with less than 15 tagSNPs and MAF>0.25 for those genes with 15 or more tagSNPs (COMT, DBH and SLC6A3). A total of 69 tagSNPs (26 in multi-loci bins and 43 singletons) were chosen [see Additional File 1]. Of these, five did not pass through the SNPlex design pipeline. After genotyping, MAF were determined in our control population 1 [see Additional File 2]. To ensure that no population stratification was present in the sample, 45 anonymous unlinked SNPs located at least 100 kb distant from known genes were also analyzed [13] by means of STRUCTURE [14], FSTAT [15] and the method by Pritchard and Rosenberg [16], as previously described [17].

Genotyping was performed at the Barcelona node of the National Genotyping Center http://www.cegen.org using the SNPlex technology [18]. Two CEPH DNA samples were included in the different genotyping assays, and a concordance rate of 100% with HapMap data was obtained. The allelic variants of the SNPs under study were named on the coding strand of each gene.

Statistical analyses

The minimal statistical power, calculated post hoc in population 1 using the Genetic Power Calculator software [19], assuming a disease prevalence of 0.12, an odds ratio (OR) of 1.7, a significance level (α) of 0.05 and a MAF of 0.123, the lowest in control population, was 85% and decreased to 74% for MO and 68% for MA. Given these estimates, we decided to begin our study by performing a joint analysis of the MO and MA groups, and to proceed to separate analyses of clinical subgroups only if a positive association was obtained in the whole sample.

Individuals with <40% successful genotypes were excluded from the analysis; SNPs with >10% missing genotypes were considered as failed; SNPs at r2 > 0.85 from any other studied SNP or showing deviation from Hardy-Weinberg equilibrium (HWE; threshold set at 0.01) as calculated in our control population 1 were also excluded. The SNPs analyzed in the follow-up population were also in HWE.

Single-marker analysis

The analysis of HWE as well as case-control comparisons of both allele and genotype frequencies under a codominant model were performed with the SNPassoc R library [20] initially in population 1, adjusting by sex. When a nominal association was identified (P < 0.05), dominant and recessive models were also analyzed. The significance threshold under the Bonferroni correction for multiple testing was set at P < 5e-04 upon consideration of 50 SNPs analyzed, genotype and allele comparisons and a single clinical group. Under a False Discovery Rate (FDR) of 10% the threshold was set at P < 0.0035, using the qvalue R library [21].

Multiple-marker analysis

Risk haplotypes were assessed in the whole migraine group (MO + MA) with the UNPHASED software [22], only for genes showing association in the single-marker analysis after the Bonferroni or FDR corrections. The best up to five-marker haplotype was selected as previously described [17]. Significance was estimated by a 10,000 permutation procedure with UNPHASED [22]. The specific assignment of haplotypes to individuals was performed independently in cases and controls with the PHASE 2.0 software [23]. The comparisons of risk haplotype carriers vs non carriers were performed using the SNPassoc R library [20] adjusting by sex. Subsequently, risk haplotypes originally identified in the all-migraine group and SNPs trespassing the FDR threshold were tested in the MO and MA subgroups.

Follow-up replication study

Risk haplotypes identified in population 1 were tested in the replication population 2. For SNPs showing nominal association with migraine in population 1 a comparison of genotype and allele frequencies was undertaken in population 2.

Results

Initially, 64 SNPs from nine candidate genes encoding proteins related with DA neurotransmission were genotyped [see Additional File 1]. Fourteen SNPs were excluded from statistical analysis after data depuration in the first population: ten had genotype call rates <90%, one was monomorphic and three were in strong LD with other SNPs (r2 > 0.85, [see Additional File 1]). The 50 remaining SNPs had MAF>0.12 and were in HWE in control population 1 (P > 0.01; [see Additional File 2]). After the exclusion of individuals with low genotyping rate, population 1 consisted of 263 patients and 274 controls, and population 2 was composed of 259 patients and 287 controls. No evidence of population stratification was found in any of the two populations studied by applying the STRUCTURE software [see Additional File 3], the Fst coefficient (theta = 0, 99%CI = 0.000-0.001 for population 1 and theta = 0, 99%CI = 0.000-0.002 for population 2) and the method by Pritchard and Rosenberg (P = 0.57 for population 1 and P = 0.05 for population 2).

In the single-marker analysis, genotype and allele frequencies were compared between patients and controls in population 1 [see Additional File 2]. Six SNPs within five genes (DRD1, DRD2, DRD3, DBH and TH) displayed P-values < 0.05 (table 1). Two of them, rs2283265 in DRD2 and rs2070762 in TH, remained significant after applying a FDR of 10% (Table 1A) and were further considered for the multiple-marker analysis. No SNP withstood the restrictive Bonferroni correction for multiple testing. For these two SNPs we sought to detect a specific association with either one of the clinical subtypes, MO or MA. We found that in population 1, DRD2 rs2283265 was associated wit both MO and MA, while TH rs2070762 was only associated with MO (Table 1B).

Table 1 Results of dopamine-related genes association studies

The analysis of all possible allelic combinations within the DRD2 gene revealed a five-marker haplotype (rs12363125-rs2283265-rs2242592-rs1554929-rs2234689) associated with migraine (best adjusted P-value = 0.00889; Table 2), with an over-representation of the T-C-G-C-G allelic combination in cases (OR = 1.85, 95%CI = 1.13-3.04, P = 0.0139) and the C-A-A-C-C haplotype in controls (OR = 1.88, 95%CI = 1.25-2.82, P = 0.00199; Table 3). The T-C-G-C-G haplotype carriers displayed an OR of 1.74 (95%CI = 1.06-2.88, P = 0.0277). In order to investigate if the association was specific of MO or MA, we compared the risk haplotype carrier frequencies between controls and each migraine subgroup separately. Although the frequencies were different (10.6% controls, 17.7% MO, 16.4% MA), they only reached borderline significance in MO (P = 0.042), while no significant differences were found in MA (P = 0.12). We performed a second-stage replication study in an independent Spanish case-control cohort to test these positive findings. The frequency of the risk haplotype T-C-G-C-G (rs12363125-rs2283265-rs2242592-rs1554929-rs2234689) carriers in population 2 was compared between 259 patients and 287 controls. Control carrier frequencies were very similar to those obtained in population 1 (10.4% control population 2 and 10.8% control population 1), whereas case carriers were more frequent in population 1 (17.1%) than in population 2 (14.3%). Thereby, the differences between cases and controls in population 2 did not reach significance (P = 0.22).

Table 2 Haplotype analysis of DRD2 and TH SNPs in 263 migraine patients and 274 unrelated non-migraine controls using the UNPHASED software.
Table 3 Haplotype distributions of DRD2 in populations 1 and 2 using the UNPHASED software.

Multiple-marker analysis of the two SNPs (rs6356 and rs2070762) in TH showed a different overall distribution between cases and controls (best adjusted P-value = 0.015, Table 2), due to an over-representation the A-C allelic combination in cases (P = 0.037, OR = 1.34, 95%CI = 0.94-1.90), and G-T in controls (P = 0.00573, OR = 1.47, 95%IC = 1.11-1.94; Table 3). However, individual haplotype assignation did not identify differences in the frequency of risk haplotype carriers between cases and controls. Moreover, the analysis of rs6356-rs2070762 haplotype distributions in cases and controls of population 2 found no evidence of association with migraine (table 3).

Finally, we aimed to determine whether those variants nominally associated with the disease phenotype in population 1 after the single-marker analysis, could be replicated in population 2, especially rs2283265 in the DRD2 gene and rs2070762 in the TH gene, which maintained significance after FDR correction in the initial analysis. The comparison of genotype and allele frequencies between cases - either the whole group of migraineurs, or MO or MA subgroups- and controls did not reveal significant differences for rs2283265 (codominant genotypes P = 0.62 and alleles P = 0.36), rs2070762 (codominant genotypes P = 0.44 and alleles P = 0.22) nor for any other SNP (table 1).

Discussion

We performed a two-stage case-control association study of eight dopamine-related genes in the Spanish population. In order to capture the common haplotype variation of these genes in the European population, we selected haplotype-tagging SNPs which covered each gene and its flanking regions. In population 1, a five-marker risk haplotype in the DRD2 gene and a single variant in the TH gene were found to be associated with migraine, and both remained significant after applying correction for multiple comparisons. In the initial single-marker analysis, pointing at five genes including the two above, no SNP withstood the Bonferroni correction. However, it is well known that this correction is often over-conservative as it assumes independence of all the tests performed, whereas many SNPs within the genes studied, although not in strong LD, are not independent. When markers found associated in population 1 were analyzed in the follow-up population, the results could not be replicated. As special attention was paid to rule out the existence of stratification and both populations were comparable in terms of size, gender distribution, ethnicity (Caucasians), geographical origin (Spain) and diagnostic criteria, failure to replicate the results suggests that the associations identified in population 1 may be spurious and that the genes analyzed here would not be involved in migraine susceptibility. However, these findings should be taken with caution, as the genetic coverage of some of the studied genes is not optimal for several reasons: First, SNPs with low frequencies, which would require very large sample sizes to produce significant results, were not selected. Second, some SNPs within the studied genes, for which no LD data were available in the HapMap database, were not included. And third, SNPlex design constraints and low genotyping call rates of some specific SNPs forced additional exclusions that left the DRD4 gene out of the study. Of note, the same Spanish cohort analyzed in the present work was previously scrutinized by us to detect association of MA or MO with genes related with serotonin neurotransmission [10]. Among the three genes that displayed significant association, two belong to the dopamine metabolic pathway: MAOA, found to be associated with MO, and DDC, which was associated with MA. However, these findings still await replication.

A number of association studies have focused on dopamine-related genes. The first susceptibility polymorphism identified in this system was the NcoI variant in the DRD2 gene (rs6275), with an over-representation of the C allele in MA [24]. Subsequent studies failed to replicate this association [2528] or that with other DRD2 polymorphisms [27, 29]. It is worth mentioning, however, that a (TG)n repeat variant in DRD2 was found associated with yawning and nausea in a small subgroup of migraine patients [30]. We analyzed DRD2 rs2242592, in strong LD with rs6275, that belonged to a risk haplotype identified in population 1 but not confirmed in population 2. Subsequent studies found association between migraine phenotypes and polymorphisms in DRD4 [31, 32] and DBH [9, 25, 33, 34], although negative associations have also been described [9, 30, 31, 33]. No associations have been identified in any of the polymorphisms analyzed in genes DRD1, DRD3, DRD5 or COMT [9, 27, 30, 3538]. The genetic marker set selected in the present analysis is in many respects not comparable with the polymorphisms analyzed in the previous studies. However, our results agree with previous negative findings in DRD1, DRD3 and one polymorphism in DBH.

A recent study carried in two German populations [9], analyzed the contribution of the nine dopamine-related genes we have examined, plus DDC, to MA susceptibility. In that study, MA was associated with three SNPs, SLC6A3 rs40184, DRD2 rs7131056 and DBH rs2097629. Overall, they analyzed 43 SNPs belonging to the nine genes studied in our work; of them, 23 SNPs are coincident or in strong LD (r2 > 0.85) with the ones we analyzed. The remaining 20 markers were not included in our analysis because their MAFs values were under the selected cut off (n = 13), lacked genotyping in the HapMap sample (n = 3), had SNPlex design constraints (n = 2) or failed in the genotyping step (n = 2). Conversely, our study included 27 SNPs with MAF>0.15 that were not analyzed by Todt et al. Three out of five nominal associations identified in our population 1, not replicated in population 2, also showed P-values > 0.05 in the first German population, thus reinforcing the likelihood of a spurious association in our population 1. Our study did not reveal association with rs7131056 in DRD2 or rs40184 in SLC6A3 at variance with the German study, while rs2097629 in DBH was not included in our study because of design constraints. In addition to differences in the respective SNP sets, our samples were composed of both MO and MA patients, and therefore a comparison of our results with those of Todt et al. is not altogether straightforward. Also, our analytical design set that the two population samples could only be grouped for analysis in case nominal associations were found in both populations 1 and 2, while in the German study their two samples were analyzed as a single group for all SNPs within the three genes that showed nominal association in only one population. This strategy produced significant associations despite lack of replication in their follow-up sample. Future studies combining both marker sets might help to reconcile these apparently discordant findings.

Much evidence points to dopamine hypersensitivity in migraineurs, particularly those displaying the premonitory symptoms of yawning or nausea. In our study, such specific symptoms could not be analyzed, since they were not available in the whole sample. To our knowledge, no well-powered association study has addressed the relationship between endophenotypes based on dopaminergic symptoms and genetic susceptibility using a pathway-based approach. Alternatively, latent class analysis of migraine symptoms, as used to enhance clinical homogeneity in genetic linkage analysis [39, 40], might define migraine phenotypes, not necessarily related to ICHD-II migraine subtype diagnoses, and thus uncover specific genetic susceptibility factors.

Conclusion

In summary, our results do not support the involvement of a set of dopamine-related genes in the genetic vulnerability to migraine in the Spanish population, albeit a previous association study in the same cohort identified DDC and MAOA as potential susceptibility genes [10]. Further studies in larger samples or family-based sets may help to clarify the contribution of dopamine-related genes to migraine genetic background.

References

  1. Wessman M, Terwindt GM, Kaunisto MA, Palotie A, Ophoff RA: Migraine: a complex genetic disorder. Lancet Neurol. 2007, 6: 521-32. 10.1016/S1474-4422(07)70126-6.

    Article  CAS  PubMed  Google Scholar 

  2. Pietrobon D: Familial hemiplegic migraine. Neurotherapeutics. 2007, 4: 274-84. 10.1016/j.nurt.2007.01.008.

    Article  CAS  PubMed  Google Scholar 

  3. International Headache Society: The International Classification of Headache Disorders. Cephalalgia. 2004, 24 (Suppl 1): 9-160. 2

    Google Scholar 

  4. Sicuteri F: Dopamine, the second putative protagonist in headache. Headache. 1977, 17: 129-31. 10.1111/j.1526-4610.1977.hed1703129.x.

    Article  CAS  PubMed  Google Scholar 

  5. Chen SC: Epilepsy and migraine: The dopamine hypotheses. Med Hypotheses. 2006, 66: 466-72. 10.1016/j.mehy.2005.09.045.

    Article  CAS  PubMed  Google Scholar 

  6. Bergerot A, Storer RJ, Goadsby PJ: Dopamine inhibits trigeminovascular transmission in the rat. Ann Neurol. 2007, 61: 251-62. 10.1002/ana.21077.

    Article  CAS  PubMed  Google Scholar 

  7. Akerman S, Goadsby PJ: Dopamine and migraine: biology and clinical implications. Cephalalgia. 2007, 27: 1308-14. 10.1111/j.1468-2982.2007.01478.x.

    Article  CAS  PubMed  Google Scholar 

  8. Mascia A, Afra J, Schoenen J: Dopamine and migraine: a review of pharmacological, biochemical, neurophysiological, and therapeutic data. Cephalalgia. 1998, 18: 174-82. 10.1046/j.1468-2982.1998.1804174.x.

    Article  CAS  PubMed  Google Scholar 

  9. Todt U, Netzer C, Toliat M, Heinze A, Goebel I, Nurnberg P, Gobel H, Freudenberg J, Kubisch J: New genetic evidence for involvement of the dopamine system in migraine with aura. Hum Genet. 2009, 125: 265-79. 10.1007/s00439-009-0623-z.

    Article  CAS  PubMed  Google Scholar 

  10. Corominas R, Sobrido MJ, Ribases M, Cuenca-Leon E, Blanco-Arias P, Narberhaus B, Roig M, Leira R, Lopez-Gonzalez J, Macaya A, Cormand B: Association study of the serotoninergic system in migraine in the Spanish population. Am J Med Genet B Neuropsychiatr Genet. 2009,

    Google Scholar 

  11. Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, Nickerson DA: Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet. 2004, 74: 106-20. 10.1086/381000.

    Article  CAS  PubMed  Google Scholar 

  12. Thorisson GA, Smith AV, Krishnan L, Stein LD: The International HapMap Project Web site. Genome Res. 2005, 15: 1592-3. 10.1101/gr.4413105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Sanchez JJ, Phillips C, Borsting C, Balogh K, Bogus M, Fondevila M, Harrison CD, Musgrave-Brown E, Salas A, Syndercombe-Court D, Schneider PM, Carracedo A, Morling N: A multiplex assay with 52 single nucleotide polymorphisms for human identification. Electrophoresis. 2006, 27: 1713-24. 10.1002/elps.200500671.

    Article  CAS  PubMed  Google Scholar 

  14. Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003, 164: 1567-87.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Goudet J: Fstat version 1.2: a computer program to calculate Fstatistics. J Hered. 1995, 86: 485-486.

    Google Scholar 

  16. Pritchard JK, Rosenberg NA: Use of unlinked genetic markers to detect population stratification in association studies. Am J Hum Genet. 1999, 65: 220-8. 10.1086/302449.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Ribases M, Hervas A, Ramos-Quiroga JA, Bosch R, Bielsa A, Gastaminza X, Fernandez-Anguiano M, Nogueira M, Gomez-Barros N, Valero S, Gratacos M, Estivill X, Casas M, Cormand B, Bayes M: Association Study of 10 Genes Encoding Neurotrophic Factors and Their Receptors in Adult and Child Attention-Deficit/Hyperactivity Disorder. Biol Psychiatry. 2008, 63: 935-45. 10.1016/j.biopsych.2007.11.004.

    Article  CAS  PubMed  Google Scholar 

  18. Tobler AR, Short S, Andersen MR, Paner TM, Briggs JC, Lambert SM, Wu PP, Wang Y, Spoonde AY, Koehler RT, Peyret N, Chen C, Broomer AJ, Ridzon DA, Zhou H, Hoo BS, Hayashibara KC, Leong LN, Ma CN, Rosenblum BB, Day JP, Ziegle JS, De La Vega FM, Rhodes MD, Hennessy KM, Wenz HM: The SNPlex genotyping system: a flexible and scalable platform for SNP genotyping. J Biomol Tech. 2005, 16: 398-406.

    PubMed  PubMed Central  Google Scholar 

  19. Purcell S, Cherny SS, Sham PC: Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits. Bioinformatics. 2003, 19: 149-50. 10.1093/bioinformatics/19.1.149.

    Article  CAS  PubMed  Google Scholar 

  20. Gonzalez JR, Armengol L, Sole X, Guino E, Mercader JM, Estivill X, Moreno V: SNPassoc: an R package to perform whole genome association studies. Bioinformatics. 2007, 23: 644-5. 10.1093/bioinformatics/btm025.

    PubMed  Google Scholar 

  21. Storey J: A direct approach to false discovery rates. Journal of the Royal Statistical Society, Series B. 2002, 64: 479-498. 10.1111/1467-9868.00346.

    Article  Google Scholar 

  22. Dudbridge F: Pedigree disequilibrium tests for multilocus haplotypes. Genet Epidemiol. 2003, 25: 115-21. 10.1002/gepi.10252.

    Article  PubMed  Google Scholar 

  23. Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68: 978-89. 10.1086/319501.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Peroutka SJ, Wilhoit T, Jones K: Clinical susceptibility to migraine with aura is modified by dopamine D2 receptor (DRD2) NcoI alleles. Neurology. 1997, 49: 201-6.

    Article  CAS  PubMed  Google Scholar 

  25. Lea RA, Dohy A, Jordan K, Quinlan S, Brimage PJ, Griffiths LR: Evidence for allelic association of the dopamine beta-hydroxylase gene (DBH) with susceptibility to typical migraine. Neurogenetics. 2000, 3: 35-40.

    CAS  PubMed  Google Scholar 

  26. McCarthy LC, Hosford DA, Riley JH, Bird MI, White NJ, Hewett DR, Peroutka SJ, Griffiths LR, Boyd PR, Lea RA, Bhatti SM, Hosking LK, Hood CM, Jones KW, Handley AR, Rallan R, Lewis KF, Yeo AJ, Williams PM, Priest RC, Khan P, Donnelly C, Lumsden SM, O'Sullivan J, See CG, Smart DH, Shaw-Hawkins S, Patel J, Langrish TC, Feniuk W, Knowles RG, Thomas M, Libri V, Montgomery DS, Manasco PK, Xu CF, Dykes C, Humphrey PP, Roses AD, Purvis IJ: Single-nucleotide polymorphism alleles in the insulin receptor gene are associated with typical migraine. Genomics. 2001, 78: 135-49. 10.1006/geno.2001.6647.

    Article  CAS  PubMed  Google Scholar 

  27. Stochino ME, Asuni C, Congiu D, Del Zompo M, Severino G: Association study between the phenotype migraine without aura-panic disorder and dopaminergic receptor genes. Pharmacol Res. 2003, 48: 531-4. 10.1016/S1043-6618(03)00213-5.

    Article  CAS  PubMed  Google Scholar 

  28. Rebaudengo N, Rainero I, Parziale A, Rosina F, Pavanelli E, Rubino E, Mazza E, Ostacoli L, Furlan PM: Lack of interaction between a polymorphism in the dopamine D2 receptor gene and the clinical features of migraine. Cephalalgia. 2004, 24: 503-7. 10.1111/j.1468-2982.2004.00689.x.

    Article  CAS  PubMed  Google Scholar 

  29. Maude S, Curtin J, Breen G, Collier D, Russell G, Shaw D, Clair DS: The -141C Ins/Del polymorphism of the dopamine D2 receptor gene is not associated with either migraine or Parkinson's disease. Psychiatr Genet. 2001, 11: 49-52. 10.1097/00041444-200103000-00010.

    Article  CAS  PubMed  Google Scholar 

  30. Del Zompo M, Cherchi A, Palmas MA, Ponti M, Bocchetta A, Gessa GL, Piccardi MP: Association between dopamine receptor genes and migraine without aura in a Sardinian sample. Neurology. 1998, 51: 781-6.

    Article  CAS  PubMed  Google Scholar 

  31. Mochi M, Cevoli S, Cortelli P, Pierangeli G, Soriani S, Scapoli C, Montagna P: A genetic association study of migraine with dopamine receptor 4, dopamine transporter and dopamine-beta-hydroxylase genes. Neurol Sci. 2003, 23: 301-5. 10.1007/s100720300005.

    Article  CAS  PubMed  Google Scholar 

  32. de Sousa SC, Karwautz A, Wober C, Wagner G, Breen G, Zesch HE, Konrad A, Zormann A, Wanner C, Kienbacher C, Collier DA, Wober-Bingol C: A dopamine D4 receptor exon 3 VNTR allele protecting against migraine without aura. Ann Neurol. 2007, 61: 574-8. 10.1002/ana.21140.

    Article  PubMed  Google Scholar 

  33. Fernandez F, Lea RA, Colson NJ, Bellis C, Quinlan S, Griffiths LR: Association between a 19 bp deletion polymorphism at the dopamine beta-hydroxylase (DBH) locus and migraine with aura. J Neurol Sci. 2006, 251: 118-23. 10.1016/j.jns.2006.09.013.

    Article  CAS  PubMed  Google Scholar 

  34. Fernandez F, Colson N, Quinlan S, Macmillan J, Lea RA, Griffiths LR: Association between migraine and a functional polymorphism at the dopamine beta-hydroxylase locus. Neurogenetics. 2009, 10: 199-208. 10.1007/s10048-009-0176-2.

    Article  CAS  PubMed  Google Scholar 

  35. Shepherd AG, Lea RA, Hutchins C, Jordan KL, Brimage PJ, Griffiths LR: Dopamine receptor genes and migraine with and without aura: an association study. Headache. 2002, 42: 346-51. 10.1046/j.1526-4610.2002.02105.x.

    Article  PubMed  Google Scholar 

  36. Hagen K, Pettersen E, Stovner LJ, Skorpen F, Zwart JA: The association between headache and Val158Met polymorphism in the catechol-O-methyltransferase gene: the HUNT Study. J Headache Pain. 2006, 7: 70-4. 10.1007/s10194-006-0281-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Karwautz A, Campos de Sousa S, Konrad A, Zesch HE, Wagner G, Zormann A, Wanner C, Breen G, Ray M, Kienbacher C, Natriashvili S, Collier DA, Wober C, Wober-Bingol C: Family-based association analysis of functional VNTR polymorphisms in the dopamine transporter gene in migraine with and without aura. J Neural Transm. 2008, 115: 91-5. 10.1007/s00702-007-0799-0.

    Article  CAS  PubMed  Google Scholar 

  38. McCallum LK, Fernandez F, Quinlan S, Macartney DP, Lea RA, Griffiths Lr: Association study of a functional variant in intron 8 of the dopamine transporter gene and migraine susceptibility. Eur J Neurol. 2007, 14: 706-7. 10.1111/j.1468-1331.2007.01800.x.

    Article  CAS  PubMed  Google Scholar 

  39. Nyholt DR, Gillespie NG, Heath AC, Merikangas KG, Duffy DL, Martin NG: Latent class and genetic analysis does not support migraine with aura and migraine without aura as separate entities. Genet Epidemiol. 2004, 26: 231-44. 10.1002/gepi.10311.

    Article  PubMed  Google Scholar 

  40. Anttila V, Nyholt DR, Kallela M, Artto V, Vepsalainen S, Jakkula E, Wennerstrom A, Tikka-Kleemola P, Kaunisto MA, Hamalainen E, Widen E, Terwilliger J, Merikangas K, Montgomery GW, Martin NG, Daly M, Kaprio J, Peltonen L, Farkkila M, Wessman M, Palotie A: Consistently replicating locus linked to migraine on 10q22-q23. Am J Hum Genet. 2008, 82: 1051-63. 10.1016/j.ajhg.2008.03.003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references

Acknowledgements

We are grateful to patients and controls for their participation in the study, to M. Dolors Castellar, Anna Daví, Bernat Narberhaus, Pilar Duocastella and Alba Corrons for their help in the recruitment of patients and control subjects, and to Miquel Casas for his support. This study was funded by the Spanish Ministry of Education and Science (SAF2006-13893-C02-01 and SAF2005-07978), Fundació La Marató de TV3 (061330) and Agència de Gestió d'Ajuts Universitaris i de Recerca-AGAUR (2005SGR00848). RC was funded by the Institut de Recerca Vall d'Hebron, MR is a recipient of a Juan de la Cierva contract, EC-L is funded by Fundació La Marató de TV3, MCT is supported by the Fundación Pedro Barrié de la Maza and MJS has a contract from the Fondo de Investigaciones Sanitarias (FIS), Instituto de Salud Carlos III. SNP genotyping services were provided by the Spanish "Centro Nacional de Genotipado" (CEGEN; http://www.cegen.org).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alfons Macaya.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

RC carried out the genotyping, analyzed the data and drafted the manuscript. MC, EC-L and MJS participated in genotyping. BC and AM designed the study and had the primary responsibility for writing the manuscript. MR helped in study design and supervised all the statistical analysis. JP, SB, MJS and AM were responsible for selecting and evaluating the patients in the respective centers. EC-L, MJS and RC participated in control recruitment. All authors have read and approved the final manuscript.

Electronic supplementary material

12881_2009_509_MOESM1_ESM.DOC

Additional file 1: Supplementary table S1. Description of the SNPs initially selected for the SNPlex analysis within 9 dopaminergic candidate genes for migraine. (DOC 208 KB)

12881_2009_509_MOESM2_ESM.DOC

Additional file 2: Supplementary table S2. Hardy-Weinberg equilibrium, minimal allele frequency (MAF) and nominal P-values observed when genotype and allele frequencies of 50 SNPs within 8 candidate genes were considered in 263 migraine cases and 274 unrelated migraine-free controls. (DOC 102 KB)

12881_2009_509_MOESM3_ESM.DOC

Additional file 3: Supplementary table S3. Assessment of population stratification using 45 unlinked anonymous SNPs and the Structure v2.1 software. (DOC 35 KB)

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Corominas, R., Ribases, M., Camiña, M. et al. Two-stage case-control association study of dopamine-related genes and migraine. BMC Med Genet 10, 95 (2009). https://doi.org/10.1186/1471-2350-10-95

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2350-10-95

Keywords