Common variants in KCNK5 and FHL5 genes contributed to the susceptibility of migraine without aura in Han Chinese population

A recent genome-wide meta study suggested that rs67338227 in the FHL5 gene and rs10456100 in the KCNK5 gene are associated with migraine from 27 population-based cohorts excluding Chinese population. Given that migraine without aura (MO) is the most common subtype of migraine, our aim was to systematically investigate the relationship of common variants in FHL5 and KCNK5 genes with the susceptibility to MO and provide clues as to the nature of the mechanisms involved in the etiology of migraine. A total of 3306 subjects including 1042 patients with MO and 2264 controls were recruited for the discovery stage, and 2530 individuals including 842 patients with MO and 1688 controls for the replication stage. Twenty-two tag SNPs (7 from FHL5 and 15 from KCNK5) were selected for genotyping. Genetic associations were analyzed at both single-marker and haplotype levels. Potential functional consequences of the significant SNPs were analyzed using gene expression data obtained from the GTEx database. Two SNPs, rs10456100 (KCNK5, P = 9.01 × 10–9) and rs7775721 (FHL5, P = 6.86 × 10–13), were determined to be significantly associated with MO in the discovery sample and were then replicated in another sample. In the combined sample set, the T allele of both SNPs was significantly associated with the increased risk of MO. Significant eQTL signals were identified for both SNP rs10456100 and rs7775721. Our findings suggest that the T allele carriers of SNP rs10456100 and rs7775721 are at increased risk of migraine.


Scientific Reports
| (2021) 11:6807 | https://doi.org/10.1038/s41598-021-86374-0 www.nature.com/scientificreports/ pain-sensing. Then, the ion channel KCNK5 gene was identified as a migraine susceptibility locus. Furthermore, Olsen' research suggested that the vessels dilate during MO and contract during MA 12 . Hence, vascular etiologies play significant roles in common migraine, and the genes implicated in vascular function could also be associated with migraine susceptibility, including the four-and-a-half LIM domains protein 5 (FHL5) gene. For its multiple polyadenylation sites, FHL5 gene may stimulates vascular smooth muscle proliferation and migration via cAMP activation regulation to result in migraine. A recent GWA meta study suggested that rs67338227 in the FHL5 gene and rs10456100 in the KCNK5 gene are associated with migraine from 27 population-based cohorts excluding Chinese population 10 . Then, Lin et al. found that the A allele of rs13208321 in FHL5 gene was associated with reduced risk of migraine in a Chinese population 13 , suggesting that the FHL5 gene is a potential risk gene for migraine. An et al. also performed a multilocus analysis to reveal the association between 17 candidate genes (including FHL5 gene) and migraine susceptibility in Chinese population, but rs13208321 in FHL5 gene showed no significant association signal for the susceptibility to migraine 14 . In addition, the relationship between KNCK5 gene polymorphisms and the risk of migraine remains unknown. Given that MO is the most common subtype of migraine, we conducted this case-control study to evaluate the associations of FHL5 and KCNK5 genes with the risk of MO in a Han Chinese population. Our aim was to systematically investigate the relationship of common variants in FHL5 and KCNK5 genes with the susceptibility to MO and provide clues as to the nature of the mechanisms involved in the etiology of migraine.

Methods
Study subjects. We implemented a two-stage study design. In the discovery stage, a total of 22 SNPs were genotyped in 3306 study subjects. Significant SNPs identified in discovery stage were then genotyped in the independent samples of the replication stage comprising of 2530 study subjects. Subjects including 1042 patients with MO and 2264 healthy age-matched controls were recruited in the headache specialized clinic at the Department of Neurology of Xijing Hospital of Air Force Medical University and Xianyang Hospital of Yan'an University (Xi'an City). Subjects including 842 patients with MO and 1688 healthy age-matched controls were enrolled from the headache specialized clinic at the Department of Neurology of the Ninth People's Hospital of Chongqing (Chongqing City). Two-stage subjects included in the study were unrelated Han Chinese individuals from Xi'an City and Chongqing City. All patients with MO were diagnosed by headache specialists based on the strict neurological examination results according to the criteria of the International Classification of Headache disorders (ICHD-3rd edition). Subjects with severe organic disease, tumors, neurological disorders, histories of depression or other comorbid psychiatric disorders involving migraines and nonmigrainous headaches were excluded. In addition, patients with secondary MO post-head injury were also excluded from the study. All healthy controls with no personal chronic headache were recruited from attendees at the physical examination centers in the hospitals mentioned above. Detailed demographic and medical history questionnaires were collected from all subjects, including age, gender, family history, unilateral pain, pulsating pain, nausea, vomiting, photophobia, phonophobia, and aggravation by physical activity. Written informed consent forms were obtained from all subjects. This study was carried out according to the ethical guidelines of the Declaration of Helsinki (version 2002) and was approved by the Medical Ethics Committee of Xi'an Jiaotong University Health Science Center.
SNP selection and genotyping. SNPs located within gene regions of FHL5 and KCNK5 were extracted for genotyping. Based on 1000 genome Chinese Han Beijing (CHB) data, we selected 21 tag SNPs (from 162 SNPs with MAF > 0.05) for genotyping (6 from FHL5 and 15 from KCNK5) using the criteria of r 2 ≥ 0.5. In addition, SNP rs67338227 was also included for purposes of replication for a previous meta-study 10 . All 22 SNPs were located at intronic regions of the two genes (Supplemental Table S1). The peripheral blood samples of the study subjects were utilized for genomic DNA extractions based on the Tiangen DNA extraction kit (Tiangen Biotech Co. Ltd, Beijing, China). The Sequenom MassARRAY platform with iPLEX GOLD chemistry (Sequenom, San Diego, CA, USA) was used for SNP genotyping. To minimize artificial factors, laboratory personnel were blinded to labels of case and control for each sample during the experiments 15,16 . Data entry was also reviewed independently by two authors. Finally, 5% of the random sample was replicated for SNP genotyping, and the concordance rate was 100%.

Statistical power analyses.
We have utilized GAS power calculator (version 4.2.5, http:// csg. sph. umich. edu/ abeca sis/ gas_ power_ calcu lator/) to conduct the statistical power analyses. The parameters utilized for the power analyses were summarized in Supplemental Table S2. The results indicated that our study setting could achieve 82.8% statistical power in detecting a genetic association signal with the genotypic relative risk of 1.2 (Supplemental Figure S1).

Statistical and bioinformatic analyses. Hardy-Weinberg equilibrium (HWE) was tested in controls.
Single marker-based genetic associations were analyzed by performing χ 2 tests. Linkage disequilibrium (LD) blocks were constructed and haplotype-based association analyses were conducted for each LD block. In addition, distributions of clinical variables for patients with migraine in genotype groups of the significant SNPs were also examined. Bonferroni corrections were applied to address multiple comparisons. The threshold of P values chosen for discovery stage was 0.05/22≈0.002. Plink 17 and Haploview 18 were utilized for genetic association analyses. R project was used for general statistical analyses 19 .
To examine the potential functional consequences of the significant SNPs, we extracted gene expression data from the GTEx database to investigate the tissue-specific expression quantitative traits loci (eQTL) and

Results
Demographic and clinical characteristics of the study subjects. A total of 3306 study subjects including 1042 patients with MO and 2264 controls were recruited for discovery stage. A total of 2530 individuals including 842 patients with MO and 1688 controls were enrolled for the replication stage. Age, gender and family history of migraine were compared between controls and patients with migraine in the discovery and replication samples, and no significant differences were identified (Table 1).
Genetic association between MO and SNPs. All 22 genotyped SNPs were in HWE in controls (Supplemental Table S1). Two SNPs were identified to be significantly associated with migraine in the discovery sample and were then replicated in the replication sample (Table 2). In the combined sample set, SNP rs10456100 (KCNK5) was significantly associated with migraine (χ 2 = 37.05, P Value = 9.01 × 10 -9 , Table 2). The T allele of this SNP was significantly associated with the increased risk of migraine (OR TT Table 2). As shown in Table 2, SNP rs7775721 (FHL5) was also significantly associated with migraine (χ 2 = 56.02, P value = 6.86 × 10 -13 ) in the combined sample set, and its T allele was significantly associated with the increased risk of migraine (OR TT Table S3. LD plots for both genes were constructed (Fig. 1) and only one LD block was identified. It is a two-SNP LD block (rs3892126-rs2815118) located within gene region of KCNK5. No significant signals were identified in haplotype-based association analyses (Supplemental Table S4). For both significant SNPs, there is no significant difference for clinical variables in different genotype groups (Supplemental Table S5). To investigate the potential systematic inflation of association signals caused by population stratification, we generated Q-Q plots based on the results of allelic association analyses and the result is shown in Supplemental Figure S2. No systematic inflation for the association signals could be identified.  www.nature.com/scientificreports/ Significant eQTL and sQTL signals captured on rs10456100 and rs7775721. Significant eQTL signals were identified for both SNP rs10456100 and rs7775721 (Supplemental Table S6 and S7). The threshold of P values was 0.05/46 = 0.001. The T allele of rs10456100 was significantly associated with decreased gene expression level of KCNK5 in subcutaneous adipose tissue ( Fig. 2A, P = 2.70 × 10 -10 ). On the other hand, the T allele of rs7775721 was significantly associated with decreased gene expression of FHL5 in the putamen of the brain (Fig. 2B, P = 2.70 × 10 -10 ). In addition, we identified widespread significant sQTL signals for rs7775721 on FHL5 ( Table 3). The T allele of rs7775721 was related to the decreased level of alternative splicing of the pre-mRNA for gene FHL5 from multiple types of human tissues.

Discussion
In the present study, we identified two significant genetic association signals, SNP rs10456100 in KCNK5 and SNP rs7775721 in FHL5, of MO in a Chinese Han population. For SNP rs10456100, although a recent large-scale GWA meta-study reported this association signal, a Chinese population was not included in that study 10 . In this sense, our study could be considered a successful replication of this previous study in a Chinese population. The effect direction of SNP rs10456100 was the same for both studies, even though the effect size was considerably greater in the present study. Furthermore, SNP rs67338227 was also reported to be significant in the GWA meta-study. However, in the present study, we did not identify this signal. In addition, another study based on European population indicated that SNP rs9394578 in KCNK5 was a significant association signal for risk of MO 21 . However, this association signal was also not replicated in our Chinese samples. These inconsistencies could be explained, at least in part, by differences in genetic background of the study populations. For example, the minor allele frequency (the A allele) of SNP rs9394578 was 0.23 in a European population (1000 genome data) compared to 0.08 in our Chinese Han data. For SNP rs7775721, to the best of our knowledge, no previous study has reported this genetic association signal, although multiple genetic association studies have been conducted for genetic polymorphisms of FHL5 www.nature.com/scientificreports/ and risk of MO. SNP rs2983896 was identified to be significantly associated with risk of MO in a study based on European population in 2016 22 ; however, this genetic association signal was not replicated in the present study. Similar to the association signals in KCNK5, we believe that this inconsistency could also be explained by the genetic heterogeneity of different populations. In addition, at least two previous studies have investigated the associations between genetic polymorphisms in FHL5 and risk of MO based on Chinese populations 13,14 . SNP rs13208321 was found to contribute to the risk of MO in both studies. In the present study, this SNP was not included because it was located out of the gene region of FHL5.  www.nature.com/scientificreports/ It is beyond the scope of this study to investigate the pathological significance of the association signals. However, we may still conduct a general discussion based on previous studies and the present one. In general, association signals between an SNP and a disorder could follow one of the two following situations. The first is that this SNP might have true effects on the pathological mechanisms of the disorder. Another is that this SNP is a surrogate of some underlying ungenotyped DNA variant that has true effects on the onset and development of the disorder. One thing interesting to note is that, for both loci, multiple different SNPs have been reported to be associated with MO. Therefore, it is probable that both SNP rs10456100 in KCNK5 and SNP rs7775721 in FHL5 might only be surrogates without true effects. With the widespread application of high-throughput technology, many susceptible genes were identified to be contributed to the risk of complex diseases [23][24][25] , but it is difficult to draw solid conclusions only from SNP results [26][27][28] . Thus, more research on the genetic architecture of the two loci are still needed in the future to further investigate the pathological significance of the two association signals.
The KCNK5 gene is located at 6p21.2 with 5 exons. The protein encoded by KCNK5 gene is a member of the super family K of potassium channel proteins. Multiple lines of evidence have shown that other members of the super family K are also indicated to be linked with migraine susceptibility 29 . K + plays an important role in a variety of metabolic activities, and it has to get into the cytoplasm via potassium channels first before its functional implementation. Therefore, KCNK5 gene possibly regulates the K + -related membrane potential and excitability to induce the occurrence and development of migraine 30 . Moreover, high extracellular K + can cause vascular smooth muscle contraction, vasospasm, ischemia, and cortex neuron activity inhibition, which in turn leads to migraine 31,32 . Interestingly, in the present study, we found that the T allele of SNP rs10456100 was associated with increased risk of migraine and decreased level of KCNK5 gene expression. Based on these findings, we hypothesize that this allele might increase extracellular K + levels by reducing the transcription of the KCNK5 gene, in turn increasing the risk of migraine for individuals. More studies in future are still needed to investigate the underlying mechanisms of KCNK5 gene and pathology of migraine.
The FHL5 gene is located at 6q16.1 with 8 exons. The protein encoded by FHL5 gene is coordinately expressed with activator of cAMP-responsive element modulator (CREM) and confers a powerful transcriptional activation function. In the present study, we showed that the T allele of SNP rs7775721 was associated with increased risk of MO and decreased gene expression level of FHL5. In addition, bioinformatics evidence also indicated that the T allele of rs7775721 was related to decreased levels of alternative splicing of the pre-mRNA for gene FHL5 from multiple types of human tissues. Nevertheless, it has not been determined how to connect the functional consequences of rs7775721 and the pathology of migraine. Additional studies should be conducted for FHL5 gene to unravel the underlying pathological mechanisms.
This study suffers from several limitations. As a candidate gene-based study, it is difficult for us to adjust for population stratifications. Although the Q-Q plot showed no significant inflation of these association signals, we still cannot rule out the possibility that the significant findings might be spurious. In addition, we investigated the potential functional consequences of SNP rs10456100 and SNP rs7775721 through eQTL and sQTL data obtained from public available database. However, we need to be careful to interpret these data because samples collected by the GTEx database are not patients with MO. Gene expression patterns could be quite different in tissues of patients with migraine compared to those of healthy individuals.
In conclusion, our findings indicated that the T allele carriers of SNP rs10456100 and rs7775721 would increase the risk of MO in a Han Chinese population, providing evidence indicating that the KCNK5 and FHL5 genes contribute to the susceptibility of migraine. These findings help to elucidate the pathogenesis of migraine.