Next Article in Journal
Signaling Pathways Governing Cardiomyocyte Differentiation
Previous Article in Journal
Genetic Analysis of Egg Production Traits in Luhua Chickens: Insights from a Multi-Trait Animal Model and a Genome-Wide Association Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Polymorphisms in EAAT1 Glutamate Transporter Gene SLC1A3 Associated with Reduced Migraine Risk

1
Genomics Research Centre, Centre for Genomics and Personalised Health, School of Biomedical Sciences, Queensland University of Technology (QUT), 60 Musk Ave., Kelvin Grove, QLD 4059, Australia
2
Stem Cell and Neurogenesis Group, Genomics Research Centre, Centre for Genomics and Personalised Health, School of Biomedical Sciences, Queensland University of Technology (QUT), Kelvin Grove, QLD 4059, Australia
3
ARC Training Centre for Cell and Tissue Engineering Technologies, Queensland University of Technology(QUT), Kelvin Grove, QLD 4059, Australia
4
Max Planck Queensland Centre for the Materials Science of Extracellular Matrices, Kelvin Grove, QLD 4059, Australia
5
Centre for Biomedical Technologies, Queensland University of Technology (QUT), 60 Musk Ave., Kelvin Grove, QLD 4059, Australia
*
Author to whom correspondence should be addressed.
Work is included as part of a doctoral thesis: Albury, C.L. (2018). Using whole exome sequencing and genetic association studies to investigate common and complex migraine. (122954). [PhD thesis, QUT]. QUT ePrints. https://eprints.qut.edu.au/122954/.
Genes 2024, 15(6), 797; https://doi.org/10.3390/genes15060797
Submission received: 5 May 2024 / Revised: 11 June 2024 / Accepted: 14 June 2024 / Published: 18 June 2024
(This article belongs to the Section Molecular Genetics and Genomics)

Abstract

:
Dysfunction in ion channels or processes involved in maintaining ionic homeostasis is thought to lower the threshold for cortical spreading depression (CSD), and plays a role in susceptibility to associated neurological disorders, including pathogenesis of a migraine. Rare pathogenic variants in specific ion channels have been implicated in monogenic migraine subtypes. In this study, we further examined the channelopathic nature of a migraine through the analysis of common genetic variants in three selected ion channel or transporter genes: SLC4A4, SLC1A3, and CHRNA4. Using the Agena MassARRAY platform, 28 single-nucleotide polymorphisms (SNPs) across the three candidate genes were genotyped in a case–control cohort comprised of 182 migraine cases and 179 matched controls. Initial results identified significant associations between migraine and rs3776578 (p = 0.04) and rs16903247 (p = 0.05) genotypes within the SLC1A3 gene, which encodes the EAAT1 glutamate transporter. These SNPs were subsequently genotyped in an independent cohort of 258 migraine cases and 290 controls using a high-resolution melt assay, and association testing supported the replication of initial findings—rs3776578 (p = 0.0041) and rs16903247 (p = 0.0127). The polymorphisms are in linkage disequilibrium and localise within a putative intronic enhancer region of SLC1A3. The minor alleles of both SNPs show a protective effect on migraine risk, which may be conferred via influencing the expression of SLC1A3.

1. Introduction

The feedback of neuronal excitability and suppression is essential for maintaining functional brain activities. Elements that affect neuronal management can either lower or elevate the threshold for neuronal firing and may play a role in the pathological cause of a number of neurological conditions including a migraine, epilepsy, stroke, and transient global amnesia [1,2,3,4,5]. An altered threshold may trigger concurrent, self-propagating waves of depolarisation accompanied by massive ion fluxes, which is termed cortical spreading depression (CSD). CSD action originates in the occipital lobe, travels through interconnected neurons and glia, and terminates within the prefrontal cortex [6]. In a normoxic brain, CSD action causes momentary neural suppression, providing an in-built protective mechanism based on negative feedback loops altering excitatory amino acid release, cerebral vascular function and blood flow, cellular metabolism, O2 consumption, and protein synthesis [7,8]. In approximately one-third of people with migraines, an aura phase, which may manifest visual or other sensory disturbances, precedes the migraine headache, and is thought to be correlated with CSD [9].
Genetic variants affecting protein function and gene transcription can disrupt the regulation of these processes, and a causative association between a migraine and dysfunction in ion channels and the processes involved in maintaining ionic homeostasis has been found in some rare migraine subtypes [10,11]. In particular, CACNA1A, ATP1A2, and SCN1A are known causative genes of familial hemiplegic migraine (FHM), a rare autosomal dominant form of a migraine, which features motor aura [12,13,14]. Pathogenic variants in these genes are believed to influence the initiation and propagation of CSD action through hyperexcitability of the neuronal network [15,16,17]: calcium channel variants (CACNA1A) have been found to cause a gain-of-function effect in FHM type 1 (FHM1) through enhanced glutamate release and facilitated CSD action [15]; variants in the sodium potassium ATPase gene (ATP1A2) conversely cause a loss-of-function effect in FHM type 2 (FHM2) due to extracellular K+ accumulation and CSD formation [15,16,17]; and sodium channel (SCN1A) variants have been linked to FHM type 3 (FHM3) due to increased firing and hyperexcitability of GABAergic neurons, thereby elevating extracellular K+ and prompting CSD initiation [18,19].
In addition to the FHM genes, other genes with roles in neuronal excitability have been associated with familial migraine-related disorders [20,21]. For example, pathogenic variants in PRRT2, which are commonly involved in movement disorders, have also been found in hemiplegic migraine families and cases [22]. Similarly, pathogenic variants in the SLC4A4 and SLC1A3 genes, encoding astrocytic transporters for bicarbonate and glutamate, respectively, have also been identified in patients with hemiplegic migraines as well as other disorders [23,24,25].
A migraine is a complex neurological disorder with both rare monogenic and more common polygenic forms [21]. With respect to the latter, over 180 susceptibility variants for a migraine have now been identified, largely via genome-wide association studies (GWASs) with increasingly larger and diverse population samples [26,27]. The genes that contribute to polygenic migraine aetiology include those with neural and vascular functions. While the pathogenic variants that cause the rare Mendelian migraine disorders are mainly in genes with ion channel function or involved in ion transport in the brain, recently, some crossover has been observed of the monogenic genes with genes implicated in common migraine susceptibility, e.g., CACNA1A and PRRT2 [27,28]. Given this precedence, it is possible that variants in other ion channel or transporter genes may play a role in susceptibility to common migraine forms. SLC4A4, SLC1A3, and CHRNA4 were selected as potential candidate genes for the investigation of polymorphisms in migraine susceptibility.
The SLC4A4 gene on chromosome (chr) 4q13.3 encodes for the universally expressed sodium/bicarbonate co-transporter, NBCe1. One of its functions is the control of neuronal hyperexcitability by regulating localised synaptic pH levels within connected astrocyte structures [23]. Interestingly, a study by Suzuki et al. examining the genetic basis of proximal renal tubular acidosis (pRTA) identified five rare homozygous variants in SLC4A4 in patients with pRTA who also had hemiplegic migraines [23]. Functional studies found that constructs carrying these variants in C6 cells revealed a significant correlation between the apparent lack of NBCe1 membrane expression and migraine occurrence [23]. These findings suggest that the near-complete loss of NBCe1 function in astrocytes may cause a migraine through the dysregulation of the synaptic pH, triggering the excessive release and stimulation via neurotransmitters—ultimately leading to neuronal hyperexcitability.
The SLC1A3 gene (chr 5p13.3) encodes the glutamate transporter Excitatory Amino Acid Transporter 1 (EAAT1). EAAT1 is a master ion channel with two major functions: (i) the regulation of neurotransmitter concentrations at the excitatory glutamatergic synapses within the central nervous system [29]; and (ii) a glutamate-activated anion channel, which directly influences the excitability of presynaptic terminals in certain neurons [30]. An investigation into a patient with a sporadic hemiplegic migraine, seizures, and ataxia identified a variant (P290R) in the SLC1A3 gene causing neuronal hyperexcitability [24]. Functional expression studies of the P290R variant demonstrated a significant decrease in the expression of the mutated EAAT1 with a markedly reduced capacity for glutamate uptake [24]. Additional in vitro electrophysiological studies illustrated that this variant caused a gain-of-function effect when noticeably increased anion currents in both the presence and absence of glutamate were observed, supporting findings of EAAT1 malfunction in epilepsy [24].
The ion channel gene CHRNA4 (chr 20q13.33) encodes for an α subunit of the pentameric neuronal nicotinic acetylcholine receptor (nAchR), which is responsible for the control of sodium, potassium, and calcium concentration gradients, membrane depolarisation, and norepinephrine, serotonin, and dopamine release [31,32]. nAchR dysfunction may lower the threshold for CSD initiation and cause central sensitisation through the repeated activation of second- and third-order neurons [33]. When central sensitisation occurs, it is believed to attribute symptoms of allodynia, phonophobia, osmophobia, nausea, and vomiting [34,35,36]. While CHRNA4 dysfunction has not been previously implicated in a migraine, other studies have attributed its dysfunction to diseases including cognitive dysfunction and attention disorders, epilepsy, depression, Attention Deficit Hyperactivity Disorder (ADHD), Tourette’s syndrome, and other neurodegenerative disorders such as Alzheimer’s disease [37,38]. Interestingly, a study investigating whether the second-generation anticholinesterase drugs may have prophylactic activity in pain disorders found that the acetylcholinesterase inhibitor Donepezil was effective in preventing a migraine [39].
Given the ability of SLC4A4, SLC1A3, and CHRNA4 to affect neurotransmission and initiate the progression of neuronal hyperexcitability, we investigated whether common polymorphisms in these three ion channel genes may be associated with a migraine. In this study, a panel of 28 target SNPs from SLC4A4, SLC1A3, and CHRNA4 were genotyped in migraine cases and controls and tested for association with migraine susceptibility. The minor alleles of two SNPs within a putative intronic enhancer region of SLC1A3 were found to show a protective effect on migraine risk.

2. Materials and Methods

2.1. Sample Populations

Migraine-affected and -unaffected individuals were recruited by the Genomics Research Centre (GRC) from the South East Queensland region through standard advertisement—radio, newspapers, and flyers. All collected samples were of Caucasian origin and diagnosed, by a neurologist, with migraine with aura disorder (MA), with migraine without aura disorder (MO), or as unaffected according to the International Headache Society (IHS) standards of Classification. All migraine samples were matched with a corresponding age-, sex-, and ethnically matched control. Blood collection followed informed participation consent. DNA was extracted and purified from blood using a standard salting-out procedure. From this database, a total of 182 cases and 179 controls were used to complete the Agena MassARRAY genotyping study. A second cohort consisting of 258 cases and 290 controls was used to replicate significant findings. This research was approved by the QUT ethics committee (Ethics Approval Number 1800000611).

2.2. Candidate Gene SNP Selection and Genotyping

A total of 28 SNPs across the three genes of interest (SLC4A4, SLC1A3, CHRNA4) were selected based on four principal considerations: (i) literature reports regarding genes/polymorphisms associated with neurological disease, (ii) tagSNPs within the three candidate genes on HapMap [35] with minor allele frequencies below 25%, (iii) functional SNPs predicted by SNPredict2, and (iv) genome-wide associated SNPs highlighted by previous studies. For example, SLC4A4 rs1041453 was associated in a GWAS with vitamin D-binding protein levels [40] and rs16846575 with vitamin D levels [41]. A multiplex assay was designed using the Assay Design Suite (Agena Bioscience TM). The final design encompassed 12 SNPs in SLC4A4 (rs2602070, rs2602072, rs4353873, rs4458426, rs10938142, rs1031452, rs1453450, rs12504851, rs16846575, rs4254735, rs9997127, rs1031453), 11 SNPs in SLC1A3 (rs7728680, rs7729389, rs3776565, rs3776567, rs1428968, rs3776578, rs3776581, rs16903247, rs10491374, rs2269272, rs1529461), and 5 SNPs in CHRNA4 (rs4809538, rs4522666, rs2273504, rs6010918, rs6122429). Amplification and extension primers were purchased from Integrated DNA Technologies (IDT, Singapore). Primer sequences are available on request. PCR amplification and MassARRAY® genotyping was performed according to the manufacturer’s instructions (Agena Bioscience, San Diego, CA, USA. Briefly, initial genotyping was performed using MassARRAY® matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry on a 96-well platform (Agena Bioscience, San Diego, CA, USA). The replication genotyping study was performed using High-Resolution Melt (HRM) protocols with the pre-designed MassARRAY primers according to the manufacturer’s instructions on the QuantStudio 7 Flex Real-Time PCR System (ThermoFisher Scientific, Waltham, MA, USA).

2.3. Statistical Analysis

SPSS Statistics 18 (SPSS Inc., Chicago, IL, USA) was used to determine baseline population characteristics. Mean ± SD metrics were determined for the parametric variable age; size and respective population percentages were calculated for the migraine and control cohorts in terms of gender and migraine subtypes. PLINK v1.7 [36] was used for the statistical analysis of all genotypic data. Two SNPs (rs1031453 and rs3776581) were removed from a further analysis due to strong deviations from Hardy–Weinberg Equilibrium [HWE > 0.05]. Pearson’s chi-square tests were used to determine significant associations between selected SNPs and a migraine. A logistic regression analysis was performed to minimise covariate (age/sex) influence. LDlink [37] was then used to calculate linkage disequilibrium (LD) associations between gene-specific SNPs.

3. Results

3.1. Characteristics of Study Cohorts

The demographics for the initial (cohort 1, N = 361) and replication (cohort 2, N = 548) case–control cohort can be viewed in Table 1. In cohort 1, mean age was 44 and females were predominant in both the migraine population (N = 144, 79.1%) and control population (N = 121, 67.6%). In cohort 2, both cases and controls involved a prominent number of females (~86%) and a mean age of 54. The distribution of gender across the two populations closely resembles the literature-defined 3 females/1 male migraine prevalence ratio [38].

3.2. Genotyping of 26 SNPs across SLC4A4, SLC1A3, and CHRNA4 in Migraine Samples versus Controls

The Agena MassARRAY genotyping plex was run with cohort 1 DNA samples. Genotype details and frequencies of examined SNPs in SLC4A4, SLC1A3, and CHRNA4 are included in Table 2. Proportions of wildtype homozygotes (AA), heterozygotes (AB), and mutant homozygotes (BB) were compared within migraine and control samples for each SNP.
Results from Pearson’s chi-square test of association for the SNP alleles tested are summarized in Table 3. Two SNPs, rs3776578 [X2 (1, N = 361) = 4.54, p = 0.03311] and rs16903247 [X2 (1, N = 361) = 3.891, p = 0.04855], located in gene SLC1A3 on chromosome 5, were found to show a significant association (α ≤ 0.05) with unaffected individuals when compared to those with a migraine. The odds ratios (ORs) for the two SNPs rs3776578 and rs16903247 are 0.4285 (95% CI = 0.192, 0.954) and 0.4523 (95% CI = 0.201, 1.014), respectively. An OR < 1 indicates that the minor allele ‘C’ reduces disease risk and incurs a protective function in both instances.
Both the rs3776578 and rs16903247 SNPs were found to be in Hardy–Weinberg Equilibrium (p > 0.05) and the minor allele frequency was significantly higher in healthy individuals (Pearson’s chi-square test of association ≤ 0.05). To avoid age effects and gender bias when assessing the sole protective effect of each polymorphism, a logistic regression analysis was performed using age and gender as covariates, summarized in Table 4. The rs3776578 (p-value = 0.0338) and rs16903247 (p = 0.049) SNPs were found to be significantly associated with healthy individuals (α ≤ 0.05) under the unadjusted model. However, when adjusted for effects of age and gender, rs16903247 was no longer significant (p = 0.056), while rs3776578 remained significantly associated (p = 0.043) with healthy individuals. Interestingly, a significant association was identified for both SNPs with sex and a migraine [rs16903247 p = 0.014, rs3776578 p = 0.016], with this correlation previously reported [42]. Smoking has also been reported to be causally associated with risk of developing a migraine [43]; however, this information was not available for participants to be able to add it as a covariate in this analysis.

3.3. Linkage Disequilibrium Analysis

The linkage disequilibrium (LD) analysis was conducted to determine LD between target SNPs rs3776578 and rs16903247 and between target SNPs and other rare SNPs within SLC1A3 (MAF < 0.05). There were strong LD associations—a D’ value of 1.0 and an r2 correlation value of 0.960 were identified between the rs3776578 and rs16903247 SNPs. The examination of the LD association analysis between rs3776578 and rs16903247 and other rare SNPs within SLC1A3 identified mild-to-moderate LD associations (D’ > 0.80, r2 > 0.40) with five additional SLC1A3 intronic SNPs (MAF > 0.05): rs3776563, rs3836803, rs79029484, rs3756468, and rs3776583.

3.4. Replication Analysis

Genotyping for significant SNPs rs3776578 and rs16903247 was undertaken in a larger, independent migraine case–control cohort using an HRM assay. For SNP rs3776578, 252 cases and 267 controls had the wildtype ‘TT’ and 6 cases and 23 controls had the heterozygous ‘TC’ genotype; no mutant ‘CC’ genotypes were detected in either cases or controls. For SNP rs16903247, 252 cases and 270 controls had the wildtype ‘TT’ and 6 cases and 20 controls had the heterozygous ‘TC’ genotype, with no mutant ‘CC’ genotypes found in either cases or controls. Pearson’s chi-square test on the replication cohort is summarised in Table 5. The SLC1A3 SNPs rs3776578 [X2 (1, N = 548) = 8.232, p = 0.0041] and rs16903247 [X2 (1, N = 548) = 6.197, p = 0.0127] were found to be significantly associated with migraine disorder at α ≤ 0.05. Furthermore, the odds ratio for both SNPs was <1, which supports the finding from the initial study cohort that the minor ‘C’ allele is protective with respect to migraine risk.
A haplotype analysis was then used to further confirm rs3776578 and rs16903247 migraine association. Genotyping data from cohort 1 and 2 were combined for this analysis (N = 909). As seen in Table 6, both haplotypes H1 (CC) and H2 (TT) formed by SNPs rs3776578 and rs16903247 [X2 (1, N = 909) = 8.866, p = 0.002] were significantly associated with migraine (α ≤ 0.01), supporting the results observed for the SNP associations.

4. Discussion

The examination of common SNPs within candidate genes SLC4A4, SLC1A3, and CHRNA4, and their relationship to migraine susceptibility, found significant associations for two SLC1A3 gene SNPs, rs3776578 and rs16903247. These two SNPs were found to be in high LD, and odds ratios for each suggested a protective effect of the minor allele variant with respect to migraine susceptibility. No statistically significant associations were detected for the other two candidate ion channel genes CHRNA4 and SLC4A4. Additional genotyping of SLC1A3 rs16903247 and rs3776578 variants in a replication case–control cohort further supported association between the minor alleles at the SNP and reduced migraine risk.
The SLC1A3 gene encodes a Na+/K+-dependent glutamate transporter, which plays an integral role in excitatory neurotransmission regulation and downstream antioxidant synthesis, and actively provides protection against localised excitotoxic damage [44]. Functionally, the EAAT1 receptor regulates neuronal glutamate-dependent excitation through secondary-active glutamate transport to recapture glutamate from the synaptic cleft, while also simultaneously facilitating ion regulation as a Na+/K+ ion channel [45]. Rare pathogenic genetic variations within SLC1A3 have been found to be causal for episodic ataxia, type 6 (OMIM #612656), and may be associated with other central nervous system (CNS) pathologies, including epilepsy and a hemiplegic migraine [24,25,46,47]. Our findings suggest that common SNPs in SLC1A3 may also play a role in susceptibility to the more common migraine type with a polygenic basis.
Although the functional effects of both implicated intronic SNP sites are unknown, in silico prediction methods may suggest a mechanism behind the association. Imputation of the two SNP sites into Ensembl’s Variant Effect predictor platform [48] revealed that both intronic SNPs rs3776578 and rs16903247 are predicted to have a ‘modifier impact’ within an enhancer regulatory region. This is further supported by the observed histone H3 lysine 27 acetylation marks and clusters of DNaseI hypersensitivity documented in the ENCODE database [49]. According to the University of California Santa Cruz (UCSC) genome browser CAVIAR track [50], which presents high-confidence sets of genetic variants that are expression quantitative trait loci (eQTLs), SNPs in the region encompassing rs3776578 and rs16903247 affect SLC1A3 expression in a range of tissues. Furthermore, the JASPAR CORE 2024 track shows that rs16903247 is located in predicted transcription factor binding sites for ZNF214 and HNF1A, with the SNP site directly affecting the consensus binding site for the latter. It is therefore possible that these SNPs affect migraine susceptibility via influencing either alternative mRNA splicing or expression levels of SLC1A3 transcripts.
The finding of strong LD between both SNPs supports not only the integrity of the genotyping results but also supports the premise that protective SNPs within SLC1A3 gene enhancer regions play a role in limiting migraine susceptibility. Correlations between polymorphisms in enhancer regions and effects on gene expression have been noted in numerous functional studies. For example, an interesting study of the human renin gene found that two SNPs in strong LD, located within an enhancer region, were associated with a 45% increase in enhancer activity [51].
A study conducted by de Vries et al. (2016) evaluated genes from candidate gene association studies in migraines using the International Headache Genetics Consortium (IHGC) genome-wide association (GWA) data set (5175 patients with a migraine and 13,972 controls) [52]. The study focused on 21 genes from published candidate gene association studies and 6 additional genes from other non-GWA studies in migraines, which included the SLC1A3 gene. Results from this study found no SNPs in or near the 27 targeted genes, including SNPs previously identified as migraine-associated, to be statistically significant after Bonferroni correction [52]. De Vries et al. suggests that the lack of the replication of genetic markers may be a result of a low power and small effect size in candidate gene association studies, or that the GWAS was unable to sufficiently reflect specific LD patterns [52], either of which may apply to our targeted SNPs of interest, rs3776578 and rs16903247.
Despite the lack of supporting association evidence of rs3776578 and rs16903247 in migraine GWASs, SNP correlations with biological processes can be drawn from a GWAS performed on the British 1958 birth cohort [15]. Genotyping data were subjected to a quantitative trait association analysis based on measurements collected during a follow-up study by Strachan et al. in 2007 [53]. The quantitative trait analysis identified 959,443 significant markers in the GWAS central database, including nominal associations between SNPs rs3776578 (p = 0.0004) and rs16903247 (p = 0.001) and blood plasma fibrinogen levels. Therefore, the same SLC1A3 SNPs we identified as associated with a migraine are also associated with blood fibrinogen levels. Interestingly, relationships between fibrinogen levels and migraine susceptibility have been reported previously, although not always consistent; nevertheless, a two-sample Mendelian randomization study found a causal relationship between genetically predicted decreased fibrinogen and a migraine, particularly a migraine with aura [54].
A limitation of this study was that only a limited number of variants and ion channel or ion transporter-related genes were considered. While SNPs in the CHRNA4 locus were investigated in the current study, CHRNA7 may be a better target as some recent findings suggest that α7nAChR channels may play a role in a migraine, modulating pain and inflammatory pathways via the attenuation of glial and astrocytic cell activation [55,56].

5. Conclusions

In summary, we have identified associations between SLC1A3 polymorphisms (rs3776578 and rs16903247) and migraine risk. When we consider the location of both polymorphisms and the neuroprotective function of the SLC1A3 gene, it is conceivable that an enhanced gene function due to genotypic variation may increase the threshold for CSD action, thereby affecting and reducing susceptibility to migraine attacks. Our results suggest that further studies examining additional SLC1A3 SNPs, or other neuroprotective ion channel genes, may provide insight into migraine aetiology.

Author Contributions

Conceptualization, C.L.A., H.G.S., L.M.H. and L.R.G.; Data curation, C.L.A. and A.W.Y.L.; Formal analysis, C.L.A., A.W.Y.L. and N.K.T.; Funding acquisition, L.R.G.; Investigation, C.L.A. and H.G.S.; Methodology, C.L.A., H.G.S. and R.A.L.; Project administration, L.R.G.; Resources, R.A.L. and L.R.G.; Supervision, L.M.H. and L.R.G.; Validation, A.W.Y.L. and N.K.T.; Writing—original draft, C.L.A.; Writing—review and editing, H.G.S., A.W.Y.L., N.K.T., R.A.L., L.M.H. and L.R.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by an Australian NHMRC Project Grant, APP1122387. Cassie Albury was a recipient of an Australian Postgraduate Award (APA).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the Human Research Ethics Committee of the Queensland University of Technology (approval number: 1800000611, 2 November 2018), with appropriate consent for the patient cohort in place.

Informed Consent Statement

All participants provided informed consent before taking part in this study.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

We would like to thank all patients and control subjects who participated in this study.

Conflicts of Interest

The authors declare that they have no competing interests. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Quinette, P.; Constans, J.M.; Hainselin, M.; Desgranges, B.; Eustache, F.; Viader, F. Hippocampal modifications in transient global amnesia. Rev. Neurol. 2015, 171, 282–288. [Google Scholar] [CrossRef] [PubMed]
  2. Cui, Y.; Kataoka, Y.; Watanabe, Y. Role of cortical spreading depression in the pathophysiology of migraine. Neurosci. Bull. 2014, 30, 812–822. [Google Scholar] [CrossRef] [PubMed]
  3. Kowa, H.; Takigawa, H.; Nakashima, K. Cortical spreading depression and pain: A missing link in the pathophysiology of migraine? Rinsho Shinkeigaku = Clin. Neurol. 2014, 54, 1006–1008. [Google Scholar] [CrossRef] [PubMed]
  4. Fabricius, M.; Fuhr, S.; Willumsen, L.; Dreier, J.P.; Bhatia, R.; Boutelle, M.G.; Hartings, J.A.; Bullock, R.; Strong, A.J.; Lauritzen, M. Association of seizures with cortical spreading depression and peri-infarct depolarisations in the acutely injured human brain. Clin. Neurophysiol. 2008, 119, 1973–1984. [Google Scholar] [CrossRef] [PubMed]
  5. Lauritzen, M.; Dreier, J.P.; Fabricius, M.; Hartings, J.A.; Graf, R.; Strong, A.J. Clinical relevance of cortical spreading depression in neurological disorders: Migraine, malignant stroke, subarachnoid and intracranial hemorrhage, and traumatic brain injury. J. Cereb. Blood Flow Metab. 2011, 31, 17–35. [Google Scholar] [CrossRef] [PubMed]
  6. Petrusic, I.; Zidverc-Trajkovic, J. Cortical spreading depression: Origins and paths as inferred from the sequence of events during migraine aura. Funct. Neurol. 2014, 29, 207–212. [Google Scholar] [PubMed]
  7. Shibata, M.; Suzuki, N. Exploring the role of microglia in cortical spreading depression in neurological disease. J. Cereb. Blood Flow Metab. 2017, 37, 1182–1191. [Google Scholar] [CrossRef]
  8. Barreto, G.E.; Capani, F.; Cabezas, R. Cortical spreading depression and mitochondrial dysfunction with aging: Lessons from ethanol abuse. Front. Aging Neurosci. 2014, 6, 117. [Google Scholar] [CrossRef]
  9. Dodick, D.W. A Phase-by-Phase Review of Migraine Pathophysiology. Headache 2018, 58 (Suppl. 1), 4–16. [Google Scholar] [CrossRef]
  10. Yan, J.; Dussor, G. Ion channels and migraine. Headache 2014, 54, 619–639. [Google Scholar] [CrossRef]
  11. Estevez, M. Invertebrate modeling of a migraine channelopathy. Headache 2006, 46 (Suppl. 1), S25–S31. [Google Scholar] [CrossRef]
  12. Ducros, A.; Denier, C.; Joutel, A.; Vahedi, K.; Michel, A.; Darcel, F.; Madigand, M.; Guerouaou, D.; Tison, F.; Julien, J.; et al. Recurrence of the T666M calcium channel CACNA1A gene mutation in familial hemiplegic migraine with progressive cerebellar ataxia. Am. J. Hum. Genet. 1999, 64, 89–98. [Google Scholar] [CrossRef] [PubMed]
  13. Gritz, S.M.; Radcliffe, R.A. Genetic effects of ATP1A2 in familial hemiplegic migraine type II and animal models. Human. Genom. 2013, 7, 8. [Google Scholar] [CrossRef] [PubMed]
  14. Dichgans, M.; Freilinger, T.; Eckstein, G.; Babini, E.; Lorenz-Depiereux, B.; Biskup, S.; Ferrari, M.D.; Herzog, J.; van den Maagdenberg, A.M.; Pusch, M.; et al. Mutation in the neuronal voltage-gated sodium channel SCN1A in familial hemiplegic migraine. Lancet 2005, 366, 371–377. [Google Scholar] [CrossRef] [PubMed]
  15. Power, C.; Elliott, J. Cohort profile: 1958 British birth cohort (National Child Development Study). Int. J. Epidemiol. 2006, 35, 34–41. [Google Scholar] [CrossRef]
  16. Silberstein, S.D.; Dodick, D.W. Migraine genetics: Part II. Headache 2013, 53, 1218–1229. [Google Scholar] [CrossRef] [PubMed]
  17. de Vries, B.; Frants, R.R.; Ferrari, M.D.; van den Maagdenberg, A.M. Molecular genetics of migraine. Hum. Genet. 2009, 126, 115–132. [Google Scholar] [CrossRef] [PubMed]
  18. Cestele, S.; Scalmani, P.; Rusconi, R.; Terragni, B.; Franceschetti, S.; Mantegazza, M. Self-limited hyperexcitability: Functional effect of a familial hemiplegic migraine mutation of the Nav1.1 (SCN1A) Na+ channel. J. Neurosci. 2008, 28, 7273–7283. [Google Scholar] [CrossRef] [PubMed]
  19. Auffenberg, E.; Hedrich, U.B.; Barbieri, R.; Miely, D.; Groschup, B.; Wuttke, T.V.; Vogel, N.; Luhrs, P.; Zanardi, I.; Bertelli, S.; et al. Hyperexcitable interneurons trigger cortical spreading depression in an Scn1a migraine model. J. Clin. Investig. 2021, 131, e142202. [Google Scholar] [CrossRef]
  20. de Boer, I.; Harder, A.V.E.; Ferrari, M.D.; van den Maagdenberg, A.; Terwindt, G.M. Genetics of migraine: Delineation of contemporary understanding of the genetic underpinning of migraine. Handb. Clin. Neurol. 2023, 198, 85–103. [Google Scholar] [CrossRef]
  21. Sutherland, H.G.; Jenkins, B.; Griffiths, L.R. Genetics of migraine: Complexity, implications, and potential clinical applications. Lancet Neurol. 2024, 23, 429–446. [Google Scholar] [CrossRef] [PubMed]
  22. Riant, F.; Roos, C.; Roubertie, A.; Barbance, C.; Hadjadj, J.; Auvin, S.; Baille, G.; Beltramone, M.; Boulanger, C.; Cahn, A.; et al. Hemiplegic Migraine Associated With PRRT2 Variations: A Clinical and Genetic Study. Neurology 2022, 98, e51–e61. [Google Scholar] [CrossRef] [PubMed]
  23. Suzuki, M.; Van Paesschen, W.; Stalmans, I.; Horita, S.; Yamada, H.; Bergmans, B.A.; Legius, E.; Riant, F.; De Jonghe, P.; Li, Y.; et al. Defective membrane expression of the Na+-HCO3 cotransporter NBCe1 is associated with familial migraine. Proc. Natl. Acad. Sci. USA 2010, 107, 15963–15968. [Google Scholar] [CrossRef] [PubMed]
  24. Jen, J.C.; Wan, J.; Palos, T.P.; Howard, B.D.; Baloh, R.W. Mutation in the glutamate transporter EAAT1 causes episodic ataxia, hemiplegia, and seizures. Neurology 2005, 65, 529–534. [Google Scholar] [CrossRef] [PubMed]
  25. Sutherland, H.G.; Maksemous, N.; Albury, C.L.; Ibrahim, O.; Smith, R.A.; Lea, R.A.; Haupt, L.M.; Jenkins, B.; Tsang, B.; Griffiths, L.R. Comprehensive Exonic Sequencing of Hemiplegic Migraine-Related Genes in a Cohort of Suspected Probands Identifies Known and Potential Pathogenic Variants. Cells 2020, 9, 2368. [Google Scholar] [CrossRef] [PubMed]
  26. Grangeon, L.; Lange, K.S.; Waliszewska-Prosol, M.; Onan, D.; Marschollek, K.; Wiels, W.; Mikulenka, P.; Farham, F.; Gollion, C.; Ducros, A.; et al. Genetics of migraine: Where are we now? J. Headache Pain. 2023, 24, 12. [Google Scholar] [CrossRef]
  27. Bjornsdottir, G.; Chalmer, M.A.; Stefansdottir, L.; Skuladottir, A.T.; Einarsson, G.; Andresdottir, M.; Beyter, D.; Ferkingstad, E.; Gretarsdottir, S.; Halldorsson, B.V.; et al. Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura. Nat. Genet. 2023, 55, 1843–1853. [Google Scholar] [CrossRef] [PubMed]
  28. Hautakangas, H.; Winsvold, B.S.; Ruotsalainen, S.E.; Bjornsdottir, G.; Harder, A.V.E.; Kogelman, L.J.A.; Thomas, L.F.; Noordam, R.; Benner, C.; Gormley, P.; et al. Genome-wide analysis of 102,084 migraine cases identifies 123 risk loci and subtype-specific risk alleles. Nat. Genet. 2022, 54, 152–160. [Google Scholar] [CrossRef] [PubMed]
  29. Greenamyre, J.T.; Porter, R.H. Anatomy and physiology of glutamate in the CNS. Neurology 1994, 44, S7. [Google Scholar]
  30. Koch, H.P.; Lane Brown, R.; Larsson, H.P. The Glutamate-Activated Anion Conductance in Excitatory Amino Acid Transporters Is Gated Independently by the Individual Subunits. J. Neurosci. 2007, 27, 2943–2947. [Google Scholar] [CrossRef]
  31. Corringer, P.-J.; Taly, A.; Lestage, P.; Changeux, J.-P.; Guedin, D. Nicotinic receptors: Allosteric transitions and therapeutic targets in the nervous system. Nat. Rev. Drug Discov. 2009, 8, 733–750. [Google Scholar] [CrossRef] [PubMed]
  32. Rossman, A.C. The physiology of the nicotinic acetylcholine receptor and its importance in the administration of anesthesia. AANA J. 2011, 79, 433–440. [Google Scholar]
  33. Chapman, C.R.; Tuckett, R.P.; Song, C.W. Pain and Stress in a Systems Perspective: Reciprocal Neural, Endocrine, and Immune Interactions. J. Pain 2008, 9, 122–145. [Google Scholar] [CrossRef]
  34. Fejes, A.; Pardutz, A.; Toldi, J.; Vecsei, L. Kynurenine Metabolites and Migraine: Experimental Studies and Therapeutic Perspectives. Curr. Neuropharmacol. 2011, 9, 376–387. [Google Scholar] [CrossRef]
  35. Woolf, C.J. Central sensitization: Implications for the diagnosis and treatment of pain. Pain 2011, 152, S2–S15. [Google Scholar] [CrossRef] [PubMed]
  36. Lovati, C.; Giani, L.; Castoldi, D.; Mariotti D’Alessandro, C.; DeAngeli, F.; Capiluppi, E.; D’Amico, D.; Mariani, C. Osmophobia in allodynic migraineurs: Cause or consequence of central sensitization? Neurol. Sci. 2015, 36 (Suppl. 1), 145–147. [Google Scholar] [CrossRef]
  37. Schaaf, C.P. Nicotinic acetylcholine receptors in human genetic disease. Genet. Med. 2014, 16, 649–656. [Google Scholar] [CrossRef]
  38. Becchetti, A.; Grandi, L.C.; Cerina, M.; Amadeo, A. Nicotinic acetylcholine receptors and epilepsy. Pharmacol. Res. 2023, 189, 106698. [Google Scholar] [CrossRef] [PubMed]
  39. Nicolodi, M.; Galeotti, N.; Ghelardini, C.; Bartolini, A.; Sicuteri, F. Central cholinergic challenging of migraine by testing second-generation anticholinesterase drugs. Headache 2002, 42, 596–602. [Google Scholar] [CrossRef]
  40. Albinana, C.; Zhu, Z.; Borbye-Lorenzen, N.; Boelt, S.G.; Cohen, A.S.; Skogstrand, K.; Wray, N.R.; Revez, J.A.; Prive, F.; Petersen, L.V.; et al. Genetic correlates of vitamin D-binding protein and 25-hydroxyvitamin D in neonatal dried blood spots. Nat. Commun. 2023, 14, 852. [Google Scholar] [CrossRef]
  41. Sinnott-Armstrong, N.; Tanigawa, Y.; Amar, D.; Mars, N.; Benner, C.; Aguirre, M.; Venkataraman, G.R.; Wainberg, M.; Ollila, H.M.; Kiiskinen, T.; et al. Genetics of 35 blood and urine biomarkers in the UK Biobank. Nat. Genet. 2021, 53, 185–194. [Google Scholar] [CrossRef] [PubMed]
  42. Burstein, R.; Noseda, R.; Borsook, D. Migraine: Multiple processes, complex pathophysiology. J. Neurosci. 2015, 35, 6619–6629. [Google Scholar] [CrossRef] [PubMed]
  43. Yuan, S.; Daghlas, I.; Larsson, S.C. Alcohol, coffee consumption, and smoking in relation to migraine: A bidirectional Mendelian randomization study. Pain 2022, 163, e342–e348. [Google Scholar] [CrossRef] [PubMed]
  44. Danbolt, N.C. Glutamate uptake. Prog. Neurobiol. 2001, 65, 1–105. [Google Scholar] [PubMed]
  45. Grewer, C.; Gameiro, A.; Rauen, T. SLC1 glutamate transporters. Pflügers Arch.-Eur. J. Physiol. 2014, 466, 3–24. [Google Scholar] [CrossRef] [PubMed]
  46. de Vries, B.; Mamsa, H.; Stam, A.H.; Wan, J.; Bakker, S.L.; Vanmolkot, K.R.; Haan, J.; Terwindt, G.M.; Boon, E.M.; Howard, B.D.; et al. Episodic ataxia associated with EAAT1 mutation C186S affecting glutamate reuptake. Arch. Neurol. 2009, 66, 97–101. [Google Scholar] [CrossRef] [PubMed]
  47. Beart, P.M.; O’Shea, R.D. Transporters for L-glutamate: An update on their molecular pharmacology and pathological involvement. Br. J. Pharmacol. 2007, 150, 5–17. [Google Scholar] [CrossRef] [PubMed]
  48. Variant Effect Predictor. Available online: http://www.ensembl.org/Tools/VEP (accessed on 20 June 2017).
  49. ENCODE: Encyclopedia of DNA Elements. Available online: https://www.encodeproject.org/ (accessed on 10 June 2024).
  50. Nassar, L.R.; Barber, G.P.; Benet-Pages, A.; Casper, J.; Clawson, H.; Diekhans, M.; Fischer, C.; Gonzalez, J.N.; Hinrichs, A.S.; Lee, B.T.; et al. The UCSC Genome Browser database: 2023 update. Nucleic Acids Res. 2023, 51, D1188–D1195. [Google Scholar] [CrossRef] [PubMed]
  51. Fuchs, S.; Philippe, J.; Germain, S.; Mathieu, F.; Jeunemaitre, X.; Corvol, P.; Pinet, F. Functionality of two new polymorphisms in the human renin gene enhancer region. J. Hypertens. 2002, 20, 2391–2398. [Google Scholar] [CrossRef]
  52. de Vries, B.; Anttila, V.; Freilinger, T.; Wessman, M.; Kaunisto, M.A.; Kallela, M.; Artto, V.; Vijfhuizen, L.S.; Gobel, H.; Dichgans, M.; et al. Systematic re-evaluation of genes from candidate gene association studies in migraine using a large genome-wide association data set. Cephalalgia 2016, 36, 604–614. [Google Scholar] [CrossRef]
  53. Strachan, D.P.; Rudnicka, A.R.; Power, C.; Shepherd, P.; Fuller, E.; Davis, A.; Gibb, I.; Kumari, M.; Rumley, A.; Macfarlane, G.J.; et al. Lifecourse influences on health among British adults: Effects of region of residence in childhood and adulthood. Int. J. Epidemiol. 2007, 36, 522–531. [Google Scholar] [CrossRef] [PubMed]
  54. Guo, Y.; Rist, P.M.; Sabater-Lleal, M.; de Vries, P.; Smith, N.; Ridker, P.M.; Kurth, T.; Chasman, D.I. Association Between Hemostatic Profile and Migraine: A Mendelian Randomization Analysis. Neurology 2021, 96, e2481–e2487. [Google Scholar] [CrossRef] [PubMed]
  55. Yasin, S.; Gorucu Yilmaz, S.; Geyik, S.; Oguzkan Balci, S. The holistic approach to the CHRNA7 gene, hsa-miR-3158-5p, and 15q13.3 hotspot CNVs in migraineurs. Mol. Pain 2023, 19, 17448069231152104. [Google Scholar] [CrossRef] [PubMed]
  56. Liu, Q.; Liu, C.; Jiang, L.; Li, M.; Long, T.; He, W.; Qin, G.; Chen, L.; Zhou, J. alpha7 Nicotinic acetylcholine receptor-mediated anti-inflammatory effect in a chronic migraine rat model via the attenuation of glial cell activation. J. Pain. Res. 2018, 11, 1129–1140. [Google Scholar] [CrossRef] [PubMed]
Table 1. Descriptive characteristics of cohorts 1 and 2. Cohort 1 was used in the initial MassARRAY genotyping study and cohort 2 was used in the replication genotyping study.
Table 1. Descriptive characteristics of cohorts 1 and 2. Cohort 1 was used in the initial MassARRAY genotyping study and cohort 2 was used in the replication genotyping study.
Cohort 1
Migraine
N = 182
Control
N = 179
Female N, (%)144 (79.1)121 (67.6)
Male N, (%)38 (20.9)58 (32.4)
Age (years), mean ± SD44 ± 1644 ± 14
Migraine with aura104 (57.1)-
Migraine without aura 78 (42.9)-
Cohort 2
Migraine
N = 258
Control
N = 290
Female N, (%)221 (85.6)250 (86.2)
Male N, (%)37 (14.4)40 (13.8)
Age (years), mean ± SD54 ± 1354 ± 14
Migraine with aura213 (82.5)-
Migraine without aura45 (17.5)-
N—Number, SD—Standard Deviation.
Table 2. Prevalence of 26 SNPs across SLC4A4, SLC1A3, and CHRNA4 in migraine versus controls.
Table 2. Prevalence of 26 SNPs across SLC4A4, SLC1A3, and CHRNA4 in migraine versus controls.
SNPChromosome Position *AlleleWildtype Homozygotes,
AA (%)
Heterozygotes, AB (%)Mutant Homozygotes,
BB (%)
ABMigraineControlMigraineControlsMigraineControls
SLC4A4Chr 4:
rs260207071,250,885CA180 (98.9)173 (96.6)2 (1.1)6 (3.4)0 (0.0)0 (0.0)
rs260207271,256,316GT179 (98.4)172 (96.1)3 (1.6)7 (3.9)0 (0.0)0 (0.0)
rs435387371,428,816AG123 (67.6)123 (69.5)54 (29.7)49 (27.7)5 (2.7)5 (2.8)
rs445842671,475,673TA116 (63.7)124 (69.3)59 (32.4)51 (28.5)7 (3.8)4 (2.2)
rs1093814271,496,462GA154 (84.6)147 (82.1)28 (15.4)29 (16.2)0 (0.0)3 (1.7)
rs103145271,519,326CT95 (52.2)90 (50.6)74 (40.7)78 (43.8)13 (7.1)10 (5.6)
rs145345071,533,996CA140 (76.9)133 (74.3)38 (20.9)41 (22.9)4 (2.2)5 (2.8)
rs1250485171,545,167AG112 (61.5)103 (57.9)60 (33.0)61 (34.3)10 (5.5)14 (7.9)
rs1684657571,555,992TC118 (64.8)116 (64.8)56 (30.8)60 (33.5)8 (4.4)3 (7.7)
rs425473571,564,768TC179 (98.9)175 (97.8)2 (1.1)4 (2.2)0 (0.0)0 (0.0)
rs999712771,570,781AG180 (98.9)175 (97.8)2 (1.1)4 (2.2)0 (0.0)0 (0.0)
SLC1A3Chr 5:
rs772868036,612,951AG88 (48.4)80 (44.7)78 (42.8)72 (40.2)16 (8.8)27 (15.1)
rs772938936,613,099AG117 (64.3)116 (64.8)58 (31.9)57 (31.8)7 (3.8)6 (3.4)
rs377656536,631,551TC123 (67.6)129 (72.1)54 (29.7)44 (24.6)5 (2.7)6 (3.3)
rs377656736,632,373AC124 (84.3)137 (91.3)22 (15.0)13 (8.7)1 (0.7)0 (0.0)
rs142896836,646,844CT133 (73.5)127 (71.3)45 (24.9)43 (24.2)3 (1.6)8 (4.5)
rs377657836,651,884TC173 (95.1)159 (88.8)9 (4.9)20 (11.2)0 (0.0)0 (0.0)
rs1690324736,653,353TC173 (95.1)160 (89.4)9 (4.9)19 (10.6)0 (0.0)0 (0.0)
rs1049137436,667,477CG48 (26.4)45 (25.3)93 (51.1)90 (50.6)41 (22.5)43 (24.2)
rs226927236,687,754CT112 (61.9)120 (67.0)65 (35.9)49 (27.4)4 (2.2)10 (5.6)
rs152946136,689,261GA104 (57.1)116 (64.8)72 (39.6)55 (30.7)6 (3.3)8 (4.5)
CHRNA4Chr 20:
rs480953863,338,824GA102 (56.4)100 (56.8)73 (40.3)69 (39.2)6 (3.3)7 (4.0)
rs452266663,343,128AG71 (39.0)64 (35.7)93 (51.1)90 (50.3)18 (9.9)25 (14.0)
rs227350463,356,709GA128 (70.3)120 (67.4)52 (28.6)52 (29.2)2 (1.1)6 (3.4)
rs601091863,358,149GA116 (95.1)103 (88.8)6 (4.9)13 (11.2)0 (0.0)0 (0.0)
rs612242963,361,854CT134 (73.6)133 (74.3)45 (24.7)44 (24.6)3 (1.6)2 (1.1)
SNP—Single-nucleotide polymorphism, SLC4A4—Solute carrier family 4 member 4, SLC1A3—Solute carrier family 1 member 3, CHRNA4—Cholinergic receptor nicotinic α 4 subunit, N—Number, A—Major allele, B—Minor allele; * genomic positions refer to GRCh38.p7/hg38 assembly.
Table 3. Association testing of SLC4A4, SLC1A3, and CHRNA4 SNPs with migraine using Pearson’s Chi-Square tests.
Table 3. Association testing of SLC4A4, SLC1A3, and CHRNA4 SNPs with migraine using Pearson’s Chi-Square tests.
SNPAlleleF_AF_UX2p-ValueORL95U95
AB
SLC4A4
rs2602070CA0.00540.01672.09000.14820.32410.0641.617
rs2602072GT0.00820.01951.69100.19350.41670.1061.624
rs4353873AG0.17580.16670.10600.74471.06700.7231.573
rs4458426TA0.20050.16481.54400.21411.27100.8701.857
rs10938142GA0.07690.09770.98440.32110.76900.4571.293
rs1031452CT0.27470.27530.00020.98670.99720.7181.383
rs1453450CA0.12640.14250.40150.52630.87080.5671.336
rs12504851AG0.21980.25000.91500.33880.84510.5981.193
rs16846575TC0.19780.18440.21100.64601.09100.7521.581
rs4254735TC0.00550.01110.69490.40450.49170.0892.701
rs9997127AG0.00540.01110.70620.40070.48900.0882.686
SLC1A3
rs7728680AG0.30220.35202.03100.15410.79740.5831.089
rs7729389AG0.19780.19270.02940.86371.03300.7141.492
rs3776565TC0.17580.15640.49020.48391.15000.7761.704
rs3776567AC0.08160.04333.72900.05341.96200.9793.932
rs1428968CT0.14090.16570.85410.35540.82550.5491.24
rs3776578TC0.02470.05584.54000.0331 *0.42850.1920.954
rs16903247TC0.02470.05303.89100.0485 *0.45230.2011.014
rs10491374CG0.48080.49440.13350.71480.94700.6841.279
rs2269272CT0.20170.19270.09040.76361.05800.7061.268
rs1529461GA0.23080.19831.12700.28841.21300.7321.527
CHRNA4
rs4809538GA0.23480.23580.00090.97520.99450.8491.732
rs4522666AG0.35440.39111.03800.30830.85480.7031.405
rs2273504GA0.15380.17980.87120.35060.82950.6321.156
rs6010918GA0.02450.05372.74400.09760.44410.5601.229
rs6122429CT0.14010.13410.05550.81381.05200.1661.188
rs6122429CT0.14010.13410.05550.81381.05200.6881.608
SNP—Single-Nucleotide Polymorphism, SLC4A4—Solute Carrier Family 4 Member 4, SLC1A3—Solute Carrier Family 1 Member 3, CHRNA4—Cholinergic Receptor Nicotinic α 4 Subunit, A—Major Allele, B—Minor Allele, F_A—Frequency of Minor Allele in Migraine, F_U—Frequency of Minor Allele in Controls, X2—Chi-Square Value, OR—Odds Ratio, L95—Lower Bound of 95% Confidence Interval, U95—Upper Bound of 95% Confidence Interval. * Meets significance threshold (α ≤ 0.05).
Table 4. Logistic regression analysis results, adjusting for age/sex and independently determining covariate association with disease.
Table 4. Logistic regression analysis results, adjusting for age/sex and independently determining covariate association with disease.
TESTORSEL95U95STATp-Value
SLC1A3 rs16903247
Logistic Regression (Unadjusted)0.43810.41930.19260.996−1.9690.0490 *
Adjusted by Sex and Age 0.4470.42290.19511.024−1.9040.0569
Sex Association with Disease 1.8120.24431.12202.9242.4330.01497 *
Age Association with Disease 1.0060.00780.99081.0220.7870.4310
SLC1A3 rs3776578
Logistic Regression (Unadjusted)0.41360.41610.18300.935−2.1220.03387 *
Adjusted by Sex and Age 0.42820.41980.18810.975−2.0200.04336 *
Sex Association with Disease 1.7970.24451.11302.9022.3970.01654 *
Age Association with Disease 1.0060.00780.99081.0220.7780.4360
SNP—Single-nucleotide polymorphism, SLC1A3—Solute carrier family 1 member 3, OR—Odds ratio, SE—Standard error, L95—Lower bound of 95% confidence interval, U95—Upper bound of 95% confidence interval, STAT—Coefficient t-statistic. * Meets significance threshold (α ≤ 0.05).
Table 5. Results from Pearson’s chi-square test of association performed on the replication cohort.
Table 5. Results from Pearson’s chi-square test of association performed on the replication cohort.
SNPAlleleF_AF_UX2p-ValueORL95U95
AB
SLC1A3
rs3776578TC0.01190.04048.2320.0041 *0.28660.11580.7907
rs16903247TC0.01160.03476.1970.0127 *0.32830.13080.8241
SNP—Single-Nucleotide Polymorphism, SLC1A3—Solute Carrier Family 1 Member 3, A—Major Allele, B—Minor Allele, F_A—Frequency of Minor Allele in Migraine, F_U—Frequency of Minor Allele in Controls, X2—Chi-Square Value, OR—Odds Ratio, L95—Lower Bound of 95% Confidence Interval, U95—Upper Bound of 95% Confidence Interval. * Meets significance threshold (α ≤ 0.05).
Table 6. Results from haplotype analysis performed on the combined genotyping data from cohort 1 and cohort 2.
Table 6. Results from haplotype analysis performed on the combined genotyping data from cohort 1 and cohort 2.
SNPsHaplotypeF_AF_UX2DFp-Value
rs16903247/rs3776578H1-CC0.0180.0428.86610.002 *
H2-TT0.9820.9588.86610.002 *
H3-TC0.0000.000---
H4-CT0.0000.000---
SNP—Single-Nucleotide Polymorphism, H1—Haplotype 1, H2—Haplotype 2, H3—Haplotype 3, H4—Haplotype 4, F_A—Frequency of Minor Allele in Migraine, F_U—Frequency of Minor Allele in Controls, X2—Chi-Square Value, DF—Degree of Freedom. * Meets significance threshold (α ≤ 0.01).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Albury, C.L.; Sutherland, H.G.; Lam, A.W.Y.; Tran, N.K.; Lea, R.A.; Haupt, L.M.; Griffiths, L.R. Identification of Polymorphisms in EAAT1 Glutamate Transporter Gene SLC1A3 Associated with Reduced Migraine Risk. Genes 2024, 15, 797. https://doi.org/10.3390/genes15060797

AMA Style

Albury CL, Sutherland HG, Lam AWY, Tran NK, Lea RA, Haupt LM, Griffiths LR. Identification of Polymorphisms in EAAT1 Glutamate Transporter Gene SLC1A3 Associated with Reduced Migraine Risk. Genes. 2024; 15(6):797. https://doi.org/10.3390/genes15060797

Chicago/Turabian Style

Albury, Cassie L., Heidi G. Sutherland, Alexis W. Y. Lam, Ngan K. Tran, Rod A. Lea, Larisa M. Haupt, and Lyn R. Griffiths. 2024. "Identification of Polymorphisms in EAAT1 Glutamate Transporter Gene SLC1A3 Associated with Reduced Migraine Risk" Genes 15, no. 6: 797. https://doi.org/10.3390/genes15060797

APA Style

Albury, C. L., Sutherland, H. G., Lam, A. W. Y., Tran, N. K., Lea, R. A., Haupt, L. M., & Griffiths, L. R. (2024). Identification of Polymorphisms in EAAT1 Glutamate Transporter Gene SLC1A3 Associated with Reduced Migraine Risk. Genes, 15(6), 797. https://doi.org/10.3390/genes15060797

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop