- Research article
- Open Access
- Open Peer Review
Common variants in the regulative regions of GRIA1 and GRIA3 receptor genes are associated with migraine susceptibility
© Formicola et al; licensee BioMed Central Ltd. 2010
- Received: 21 January 2010
- Accepted: 25 June 2010
- Published: 25 June 2010
Glutamate is the principal excitatory neurotransmitter in the central nervous system which acts by the activation of either ionotropic (AMPA, NMDA and kainate receptors) or G-protein coupled metabotropic receptors. Glutamate is widely accepted to play a major role in the path physiology of migraine as implicated by data from animal and human studies. Genes involved in synthesis, metabolism and regulation of both glutamate and its receptors could be, therefore, considered as potential candidates for causing/predisposing to migraine when mutated.
The association of polymorphic variants of GRIA1-GRIA4 genes which encode for the four subunits (GluR1-GluR4) of the alpha-amino-3- hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) receptor for glutamate was tested in migraineurs with and without aura (MA and MO) and healthy controls.
Two variants in the regulative regions of GRIA1 (rs2195450) and GRIA3 (rs3761555) genes resulted strongly associated with MA (P = 0.00002 and P = 0.0001, respectively), but not associated with MO, suggesting their role in cortical spreading depression. Whereas the rs548294 variant in GRIA1 gene showed association primarily with MO phenotype, supporting the hypothesis that MA and MO phenotypes could be genetically related. These variants modify binding sites for transcription factors altering the expression of GRIA1 and GRIA3 genes in different conditions.
This study represents the first genetic evidence of a link between glutamate receptors and migraine.
- Risk Allele
- Migraine Patient
- Migraine With Aura
- Cortical Spreading Depression
Migraine is a multifactorial disorder in which genetic factors play a relevant role in both predisposing and determining underlying mechanisms. Approaches to identified genes for monogenic subtypes migraine (e.g. Familial Hemiplegic Migraine - FHM) has been successful (CACNA1A, ATP1A2, SCN1A genes causing FMH 1, 2 and 3 respectively) [1–4]. Conversely, finding genes for the most frequent types of migraine (with and without aura) and defining their pathogenic role has proven much more difficult . However, synaptic hyper-excitability has been invocated as one major pathogenic mechanism in both common and monogenic forms of migraine, and data from animal (experimental cortical spreading depression; c-fos protein expression at the trigeminal nucleus caudalis, plasma protein extravasations at level of dura-mater circulation, electrophysiological studies) [6, 7] and human (assessment of glutamate concentration in plasma, platelets and cerebrospinal fluid in migraneurs vs. control, glutamate and migraine symptoms; glutamate and sensitization)  studies suggest that glutamate is one of the principal factors involved.
Glutamate is the major excitatory neurotransmitter in the central nervous system which exerts this action mainly through its interaction with its ionotropic or metabotropic receptors.
In this study we investigated the association of GRIA1, GRIA2, GRIA3 and GRIA4 genes that encode for the four subunits (GluR1-GluR4) of the alpha-amino-3-hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) ionotropic receptor and their variants to migraine with and without aura (MA and MO).
This is the first study that shows a significant association between genetic variants of glutamate receptor and MA.
Population Collection & Phenotype
Two hundred fifty outpatients consecutively observed from January 2006 to December 2008 at the Service of Diagnosis and Therapy of Headache of the Second University, Naples and two hundred sixty healthy volunteers consecutively enrolled during the same time interval were included in the study. Cases and controls were matched for sex, came from the same small geographical area and were Caucasian. All gave informed consent to the research procedures including molecular genetic analyses. The study protocol included history, direct clinical and neurological examination, MRI with both arterial and venous angiography sequences; venous blood collection for routine blood analysis as well as DNA extraction from lymphocytes according to standard salting out procedure . Migraineurs were diagnosed as having either migraine with (MA) or migraine without (MO) aura, according to the 2004 International Headache Society  diagnostic criteria. The study population was comprised of 244 migraineurs (M/F 56/188; age 14-64 years; MA/MO 135/109); 7 patients who resulted affected with Familial Hemiplegic Migraine were excluded from the study.
For the statistical analysis three experimental groups were considered: all migraineurs together (MA+MO= M) and MA and MO separately. Healthy control individuals (204 females and 56 males aged 21-64 years) of Caucasian origin matching for age and sex the migraineurs served as control for all three experimental groups.
Polymorphisms selection and Genotyping
SNP selection was performed inspecting the HapMap database PhaseIII/Rel 2, Feb09 on NCBI B36 assembly, dbSNP b126 http://www.hapmap.org/index.html.en. Default parameters in Haploview were used to create haplotype blocks (95% confidence intervals) as described by Gabriel et al (2002) .
TagSNPs were selected at an r 2 threshold of 0.80 from all SNPs with minor allele frequency (MAF) > 0.15 allowing us to capture more than 80% of the haplotypic diversity of these genes.
For the GRIA2 gene, that expands a genomic region of 143 kb on chromosome 4q32.1, two main blocks of LD were identified. We selected four SNPs (two for each block): rs13138842, rs11946396, rs28495775 and rs11728471 that were genotyped by restriction digestion using the enzymes MluI, AccII, HuncII, and HphI respectively (figure 1).
The GRIA3 gene spans a genomic region of about 300 kb consisting of 17 exons, and is located in Xq24-q28 region, a critical region for migraine . In a previous unpublished study we completely sequenced the coding and regulative regions of the GRIA3 gene in ten Italian migraine patients. We detected only a synonymous sequence variation in coding region C/T N400N (rs502434), and a number of SNPs located in intronic and promoter region. Taking in account SNPs distribution and LD-plot for GRIA3 gene, in a first step we selected 8 SNPs located in different LD-blocks: -1952T/C SNP (rs3761555) in the promoter region, the G/A SNP rs12559450 in intron 2, the G/A SNP rs4825840 in intron 2, the C/A SNP rs4825849 in intron 3, the T/C SNP rs5911598 in intron 5, the C/T SNP N400N (rs502434) in exon 9, the C/T SNP rs2354037 in intron 12 and the T/C SNP rs5910006 in intron 14 (figure 1). In a preliminary analysis, using a smaller size cohort of migraine patients and controls (100 cases vs 100 controls) only one SNPs rs3761555 showed a trend of allelic differences between cases and controls, than was analysed in more detail. Genotyping was performed by restriction digestion using the enzyme BglII.
Also for the GRIA4 gene, which spans over 370 kb on chromosome 11q22.3, only two main blocks of LD were identified. We selected four SNPs (two for each block): rs11226805, rs590928, rs649098 and rs589104 that were genotyped by restriction digestion using the enzymes DdeI, Eco57I, MboI, and NlaIII respectively (figure 1).
All PCRs were performed as described in our previous work .
Briefly, PCR reactions (10 μl final volume) containing 2 mmol/L MgCl2, 0.5 mol/L of each primer, 200 mol/L dNTPs, 1 unit of Taq polymerase (Fermentas) and approximately 20 ng of genomic DNA were undertaken for genotyping purposes.
primer sets, melting temperature (MT), PCRs size and genotyping conditions are reported for associated SNPs in GRIA1 and GRIA3 genes.
ALLELE 1 (size)
ALLELE 2 (size)
To detect association between the tested markers and migraine, we performed chi-square (χ2) analysis to test for significant differences in allele and genotype frequencies in case versus control results. χ2 provides the likelihood of a deviation in the distribution of the same attributes in different classes (e.g. allelic frequencies in controls versus affected subjects). If the probability (p-value) of an equal distribution between the two groups is below a determined significance level α (in percent), the statistical output will show enough significance to assume LD and therefore association. The p value of this test is asymptotically equal to the p value obtained in Armitage's trend test (ATT). ATT takes into account genotypes rather than alleles, avoiding a possible bias caused by doubling the sample size. It assumes an additive (or co-dominant) disease model where all disease alleles have equal and independent contribution to disease risk. The ATT is valid and powerful under a broad range of disease mechanisms. Second, allele and genotype frequencies, dominant and recessive genetic models as well as odds ratios were calculated to characterise the distribution of a distinct genotypes in different phenotypic subgroups of the population (using FINETTI, http://ihg.gsf.de/cgi-bin/hw/hwa1.pl). Hardy-Weinberg Equilibrium (HWE) was also calculated using FINETTI. The analysis of the statistical power was performed using the Genetic Power Calculator software http://pngu.mgh.harvard.edu/~purcell/gpc/cc2.html assuming a significance level (α) of 0.05, a lifetime risk of 2, and the MAF of each SNP as calculated in our control group. In the population-based association analysis, the nominal significance threshold was set at P < 0.05 and lowered to P < 0.0017 after the multiple comparison correction of Bonferroni considering 22 SNPs, the three clinical groups analyzed, genotypes-alleles and females-males. Confidence intervals at 95% are provided for all odds-ratio (OR) values. Haplotype frequencies were calculated using the Estimation Haplotype (EH) program (Jurg Ott, Rockefeller University). A P level of less than 0.05 was considered statistically significant.
Electrophoretic mobility shift assays
To determine whether the c.-1952T>C SNP in GRIA3 promoter region alters nuclear protein-DNA interactions, EMSAs was performed using radio-labelled oligonucleotide probes containing either -1952T or -1952C. Nuclear extracts were prepared from HEK 293 and SH5Y5Y cell lines as described by Granelli-Piperno et al. . Binding reactions mixture contained 10 ng 32P-labeled oligonucleotide probe (105 cpm), 15 μg nuclear extract, 2 μg poly (dI-dC), 100 mM KCl, 10 mM MgCl2, 20 mM HEPES (pH 7.9), 0.2 mM EDTA, 0.5 mM dithiothreitol, 20% glycerol and protease inhibitor, in a total vol of 20 μl. To perform super shift experiments, 2 μg of the anti-HSF1 antibody (Celbio) was added to the binding reaction mixture. In competition experiments, 100-fold molar excess of cold oligonucleotides was added to the binding reaction mixture. Reaction mixtures were incubated for 30 min at room temperature, resolved on non-denaturing 5% polyacrylamide gel, dried, and exposed to autoradiography. The same assay was performed for the c.-2012C>T SNP rs548294 in GRIA1 promoter region and for c.+561G>A SNP rs2195450 in intron 1.
Primers used for GRIA3 SNP were:
Primers used for GRIA1 SNP rs548294 were:
Primers used for GRIA1 SNP rs2195450 were:
To test the functional activity of the T>C polymorphism in the promoter region of the GRIA3 gene, we first identified the putative start site at -823 to the ATG, using the Promoter scan software vs2. A fragment of 1854 bp of the GRIA3 promoter encompassing base pairs -1411 to +325 (referred to the predicted start site) and including the T>C polymorphism was amplified from genomic DNA of migraine patients with genotype CC and controls with genotype TT. The primers used for generating this fragment were sense: 5'-GCGCGCCTCGAGTTTGAAGATGAGAGAACTGG-3' and antisense: 5'-GCGCGCAAGCTTCAAAAGAAAGAGAACGAAAG-3' which contain 5'-XhoI and 3'-HindIII restriction sites (highlighted bold sequences), respectively.
PCR products were double digested using XhoI and HindIII and then legated into a linearized pGL3-Basic vector (Promega, Madison USA) containing a Firefly Luciferase reporter gene. The inserts of the constructs were verified by sequencing before transient transfection.
Transfections were carried out using HEK293T cells. Cells were co-transfected with 1 μg of reporter plasmid and 10 ng of pRL-TK control vector (Promega, Madison USA) containing Renilla luciferase, which was used to correct for transfection efficiency. Transfection cells were harvested after 24 h and cell lysates were assayed sequentially for Firefly and Renilla luciferase activity using a dual-reporter assay system (Promega, Madison USA). Transfections were performed in quadruplicate for each plasmid construct. Changes in gene transcription were calculated relative to the corrected luciferase activity of the empty pGL3-Basic vector (control).
This research was reviewed and approved by the University of Naples Human Research Ethics Committee and all subjects participating in the study gave informed consent.
GRIA1 and GRIA3genes are associated with migraine
We selected SNPs for each gene (GRIA1-GRIA4) considering their distribution in LD blocks (figure 1). All variants were genotyped in a well-characterized panel of 244 cases and 260 controls (see methods). In controls, all SNPs analyzed are in Hardy-Weinberg equilibrium (HWE) (with the exception of rs2195450) and showed allele frequencies comparable to those observed in the Caucasian population and reported in the HapMap or SNP browser databases.
No significant association with the disease for any of the four SNPs in GRIA2 and GRIA4 genes tested was observed.
Allele and genotype frequencies of GRIA1 and GRIA3 polymorphisms in Italian patients and controls and association results
11 vs 12
11 vs 22
11 vs 12+22
ARMITAGE'S TEST -P VALUE
(OR 1,7/CI 1,2-2,1)
(OR 1,7/CI 1,1-2,5)
(OR 2,4/CI 1,5-4,1)
(OR 1,9/CI 1,3-2,7)
(OR 1,6/CI 1,1-2,1)
(OR 2,3/CI 1,3-4,1)
(OR 1,8/CI 1,3-2,5)
(OR 2,6/CI 1,6-4,5)
(OR 2,7/CI 1,3-5,3)
(OR 2,7/CI 1,6-4,3)
(OR 1,6/CI 1,2-2,2)
(OR 1,6/CI 1,0-2,5)
(OR 2,4/CI 1,3-4,4)
(OR 1,8/CI 1,2-2,7)
(OR 1,6/CI 1,2-2,0)
(OR 2,1/CI 1,3-3,3)
(OR 1,4/CI 1-2)
(OR 1,9/CI 1,4-2,5)
(OR 2,7/CI 1,6-4,6)
(OR 1,6/CI 1,0-2,4)
(OR 1,3/CI 1-1,8)
(OR 1,9/CI 1,1-3,2)
(OR 2,3/CI 1,3-4)
(OR 3,3/CI 1,3-8,6)
(OR 2,6/CI 1,2-5,7)
(OR 1,7/CI 1,2-2,3)
(OR 1,8/CI 1,2-2,8)
(OR 2,3/CI 1,1-4,6)
(OR 1,9/CI 1,3-2,9)
(OR 2/CI 1,4-2,8)
(OR 2,3/CI 1,3-3,7)
(OR 3/CI 1,4-6,4)
(OR 2,4/CI 1,5-3,8)
(OR 2,46/CI 1,07-5,6)
The analysis of the +561 G>A SNP (rs2195450) also showed significant association in both allelic and genotype distributions (p = 0.0002 and p = 0.0007, respectively) (p = 0.005 and p = 0.02 after Bonferroni corrections). Stratified analyses of migraine subtypes showed a specific association with MA subtype group for both allelic and genotypic frequencies (p = 0.00002 and p = 0.00008) (p = 0.0005 and p = 0.002 after Bonferroni corrections) (Table 2). However, the SNP rs2195450, was the only SNP found in Hardy-Weinberg disequilibrium (DHW) both in controls and in patients (p < 0.05).
Frequency of common GRIA1 haplotypes in migraine patients compared with healthy controls.
CASES N = 452
CONTROLS N = 504
95% CONFIDENCE INTERVAL
Combined genotypic frequencies of GRIA1 (CT) and GRIA3 (TC)
C ONTROLS N = 260
CASES N = 244
TWO RISK ALLELES
ONE RISK ALLELE (GRIA1)
ONE RISK ALLELE (GRIA3)
NO RISK ALLELE
-1952 T>C variant in GRIA3regulative region have functional significance
Functional analysis of the associated GRIA1variants
Bioinformatics predictions did not show any clear alteration of transcription factor binding sites both for SNP -2012C>T and +561G>A. To determine whether these variations could alter nuclear protein-DNA interactions we performed EMSAs assays using radiolabelled oligonucleotide probes containing either -2012C or -2012T and +561G or +561A (see materials and methods) (figure 2B). No binding affinity was revealed for both +561G and +561A alleles. In contrast, two distinct complexes were revealed on incubation of the probes -2012T and -2012C with nuclear extract from unstimulated HEK293T cells and in particular, the -2012T probe showed stronger DNA protein-binding activity compared with the -2012C probe (figure 2B). The binding specificity and the differences in binding affinity between the T and C allele was also demonstrated in competition assay using unlabeled -2012T or -2012C oligonucleotide. The results revealed that unlabeled -2012T oligonucleotide but not -2012C with 100- fold molar excess fully blocked the binding of the radiolabeled -2012T probe or -2012C probe with nuclear protein(s) (figure 2B).
Glutamate is a major excitatory neurotransmitter in the central nervous system, thought to be widely involved in migraine mechanisms. Indeed, migraine brain nociceptive centres, including the trigeminal ganglion, trigeminal nucleus caudalis, and thalamus, contain glutamate-positive neurons, and glutamate activates the trigeminal nucleus caudalis. Glutamate is implicated in cortical spreading depression (CSD), trigemino-vascular activation, and central sensitization . These observations argue for a strong link between migraine and the glutamatergic system. A link that is important to further characterization and understanding of migraine mechanisms in order to deliver effective therapies.
In this study we have investigated the role of glutamate receptors, ionotropic, AMPA (GRIA1, GRIA2 GRIA3 and GRIA4 genes) and their variants in MA and MO. We found two variants in GRIA1 gene (rs548294 and rs2195450) and one variant in GRIA3 gene (rs3761555), all located in the regulative region of the genes, to be associated with the migraine phenotype in an Italian population. Interestingly, both the GRIA1 and GRIA3 variants were specifically associated with MA suggesting a role of these genes in CSD, which is retained to be the path physiological mechanism-underlying aura. However, to date we cannot exclude that these variant can also be associated to MO phenotype due to limited size of MO subgroup. In fact, the rs548294 SNP in GRIA1 gene was primarily associated with MO rather than MA. Although MO and MA are separate categories according to the ICHD-II criteria these two subtypes of migraine share some path physiological mechanisms . Several studies support that the same SNP or different SNPs in the same gene could be associated with both MO and MA phenotypes .
We analyzed the distribution of the haplotypes formed by the associated SNPs in GRIA1 gene in cases and controls. Haplotype-specific testing showed that the putative risk haplotype T_A was more frequent in cases than in controls (13% and 7%, respectively, P = 0.001) and was associated with increased risk of migraine (OR = 2,10; 95% CI (1.30- 3,12)) compared with the risk of the two alleles taken separately. Furthermore, our analysis of the combined effects of the two SNPs in GRIA1 (rs548294) and GRIA3 (rs3761555) genes not supported a synergy of the two SNPs in the genetic susceptibility to migraine. Indeed, in migraine patients, the combined two risk allele genotypes exhibited an OR that was more or less the sum of the ORs of each risk allele taken separately, suggesting that, when both risk alleles are present, the risk to develop migraine is merely the sum of the effect carried by each single risk allele.
Although our data need to be confirmed in more numerous series as well as in populations of different origin, the identification of a statistically significant genetic link between glutamate receptor genes and MA needs to be highlighted. A growing number of scientific report based on preclinical and clinical data argue strongly for a role of glutamatergic receptor activation in migraine and in particular in CSD . Despite the remarkable advancing in the pharmacotherapy, there is still a tremendous need for the research of more effective treatments. Glutamate receptors antagonists have been recently proposed as major target addressed for the treatment of migraine. Glutamate receptors indeed represent a promising target for a valuable, non-vasoactive oriented approach to the treatment of migraine .
The SNP rs2195450, which shows the strongest association, is in Hardy-Weinberg disequilibrium (DHW) both in controls and in patients. DHW may be due to several factors including genotyping errors, inappropriate population stratification and selection, inbreeding, presence of a causative allele or simply the chance. We can exclude with high confidence that DHW as well as patient-control differences could be due to genotyping errors. For all SNPs, about 20 subjects were randomly selected and sequenced; for rs2195450 the sample was increased to 100 subjects. The original classification was confirmed in all samples. Moreover, for the autosome genes, the presence of DHW in different direction in controls and in patients, as observed for rs2195450, is compatible with the segregation of a causative variant according to a dominant model . However, we cannot exclude that DHW is due to other confounding factors, and then the association of this SNP must be considered with prudence.
Further remarkable in our cohort is that GRIA1 and GRIA3 variants are associated to migraine in female but not in male subjects. To date we cannot assert that AMPA receptor genes and genetic variations within them may act differently in the two sexes, in light of the small sample size of case males that do not provide enough statistical power to detect association.
However, the underlying mechanism by which these variants confer susceptibility to migraine remains unclear. Bioinformatics predictions showed that the C variant of GRIA3 gene could affect putative binding sites for transcription factors altering the consensus sequence. We strengthened this hypothesis by transfecting the HEK293T cells with a promoter region of GRIA3 gene. It should be highlighted that, under stress condition the luciferase activity generated by the C construct (associated allele) was 1.3-fold greater than that generated by the T construct. Moreover, EMSA assays suggest that the -1952T>C polymorphism can affect the binding affinity of the GRIA3 promoter with transcription factor(s) and the variant T allele has a stronger binding strength compared with the C allele. Considering that the T-allele displays lower promoter transcription activity, we infer that the nuclear factor(s) bound to the region covering the -1952T>C polymorphism may function as negative regulator(s) for GRIA3 transcription. Therefore, the -1952T>C polymorphism may affect the GRIA3 promoter binding affinity with nuclear proteins and in turn the GRIA3 expression, which consequently modulates the individual's susceptibility to migraine. Future genetic and functional studies are necessary to further elucidate the role of GRIA3 in migraine susceptibility. Preliminary results using an Australian case-control panel indicate that the same -1952T>C polymorphism showed a trend of association with migraine (p = 0.04) (L. Griffiths personal communication).
More complex is the case of GRIA1 gene where only the -2012C>T variant seems to affect the binding affinity of the GRIA1 promoter with transcription factor(s) and the variant T allele has a stronger binding strength compared with the C allele. Considering that the T-allele is also the risk allele we can infer that nuclear factor/s bound to the region covering the -2012C>T polymorphism may function as positive regulator(s) for GRIA1 transcription.
These results indicate an involvement of GRIA1 and GRIA3 genes in the susceptibility to migraine with aura and encourage replication in other data sets. Moreover, functional investigations to explain the role of glutamate receptors variants in MA might be useful.
Migraine Research Foundation to T. Esposito supported this research. The Short-term Mobility Program from the Italian National Research Council also supported F. Gianfrancesco and T. Esposito. The authors would like to express their most sincere gratitude to Dr. R.T. Ngomba for critical reading of the manuscript, Dr. Serena Mirra for support in luciferase assay and Dr. Maria Mancini for support in EMSA assay.
- Kors EE, van den Maagdenberg AM, Plomp JJ, Frants RR, Ferrari MD: Calcium channel mutations and migraine. Curr Opin Neurol. 2002, 15 (3): 311-316. 10.1097/00019052-200206000-00014.View ArticlePubMedGoogle Scholar
- De Fusco M, Marconi R, Silvestri L, et al: Haploinsufficiency of ATP1A2 encoding the Na+/K+ pump alpha2 subunit associated with familial hemiplegic migraine type 2. Nat Genet. 2003, 33 (2): 192-196. 10.1038/ng1081.View ArticlePubMedGoogle Scholar
- Dichgans M, Freilinger T, Eckstein G, et al: Mutation in the neuronal voltage-gated sodium channel SCN1A in familial hemiplegic migraine. Lancet. 2005, 366 (9483): 371-374. 10.1016/S0140-6736(05)66786-4.View ArticlePubMedGoogle Scholar
- de Vries B, Frants RR, Ferrari MD, van den Maagdenberg AMJM: Molecular genetics of migraine. Hum Genet. 2009, DOI 10.1007/s00439-009-0684-zGoogle Scholar
- Colson NJ, Fernandez F, Lea RA, Griffiths LR: The search for migraine genes: an overview of current knowledge. Cell Mol Life Sci. 2007, 64: 331-344. 10.1007/s00018-006-5592-y.View ArticlePubMedGoogle Scholar
- Bergerot A, Holland PR, Akerman S, et al: Animal models of migraine: looking at the component parts of a complex disorder. Eur J Neurosci. 2006, 24 (6): 1517-34. 10.1111/j.1460-9568.2006.05036.x.View ArticlePubMedGoogle Scholar
- Tottene A, Conti R, Fabbro A, et al: Enhanced excitatory transmission at cortical synapses as the basis for facilitated spreading depression in Ca(v)2.1 knockin migraine mice. Neuron. 2009, 61: 762-773. 10.1016/j.neuron.2009.01.027.View ArticlePubMedGoogle Scholar
- Ramadan NM: The link between glutamate and migraine. CNS Spectr. 2003, 8: 446-449.View ArticlePubMedGoogle Scholar
- Miller SA, Dykes DD, Polesky HF: A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988, 16: 1215-10.1093/nar/16.3.1215.View ArticlePubMedPubMed CentralGoogle Scholar
- Headache Classification Committee of the International Headache Society: The International Classification of Headache Disorders, 2nd Edition. Cephalalgia. 2004, 24 (Supplement 1): 1-52.Google Scholar
- Gabriel SB, Schaffner SF, Nguyen H, et al: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-9. 10.1126/science.1069424.View ArticlePubMedGoogle Scholar
- Nyholt DR, Curtain RP, Griffiths LR: Familial typical migraine: significant linkage and localization of a gene to Xq24-28. Hum Genet. 2000, 107: 18-23. 10.1007/s004390050004.PubMedGoogle Scholar
- Fernandez F, Esposito T, Lea RA, Colson NJ, Ciccodicola A, Gianfrancesco F, Griffiths LR: Investigation of GABA A Receptors Genes and Migraine Susceptibility. BMC Medical Genetics. 2008, 16: 9-109.Google Scholar
- Purcell S, Cherny SS, Sham PC: Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits. Bioinformatics. 2003, 19: 149-150. 10.1093/bioinformatics/19.1.149.View ArticlePubMedGoogle Scholar
- Granelli-Piperno A, Nolan P, Inaba K, Steinman RM: The effect of immunosuppressive agents on the induction of nuclear factors that bind to sites on the interleukin 2 promoter. J Exp Med. 1990, 172: 1869-1872. 10.1084/jem.172.6.1869.View ArticlePubMedGoogle Scholar
- Heinemeyer T, Wingender E, Reuter I, et al: Databases on transcriptional regulation: TRANSFAC, TRRD and COMPEL. Nucleic Acids Res. 1998, 26: 362-367. 10.1093/nar/26.1.362.View ArticlePubMedPubMed CentralGoogle Scholar
- Vikelis M, Mitsikostas DD: The role of glutamate and its receptors in migraine. CNS Neurol Disord Drug Targets. 2007, 6: 251-7. 10.2174/187152707781387279.View ArticlePubMedGoogle Scholar
- Nyholt DR, Gillespie NG, Heath AC, Merikangas KR, 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-244. 10.1002/gepi.10311.View ArticlePubMedGoogle Scholar
- Dehbandi S, Speckmann EJ, Pape HC, Gorji A: Cortical spreading depression modulates synaptic transmission of the rat lateral amygdala. Eur J Neurosci. 2008, 27 (8): 2057-65. 10.1111/j.1460-9568.2008.06188.x.View ArticlePubMedGoogle Scholar
- Andreou AP, Goadsby PJ: Therapeutic potential of novel glutamate receptor antagonists in migraine. Expert Opin Investig Drugs. 2009.Google Scholar
- Wittke-Thompson JK, Pluzhnikov A, Cox NJ: Rational inferences about departures from Hardy-Weinberg equilibrium. Am J Hum Genet. 2005, 76: 967-986. 10.1086/430507.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/11/103/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.