Next Article in Journal
Maintaining the Quality of Mechanical Thrombectomy after Acute Ischemic Stroke in COVID(-)19 Patients
Previous Article in Journal
Neuroprotective Effect of Angiopoietin2 Is Associated with Angiogenesis in Mouse Brain Following Ischemic Stroke
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Whole-Exome Sequencing Reveals Migraine-Associated Novel Functional Variants in Arab Ancestry Females: A Pilot Study

1
Department of Medical Laboratory Sciences, College of Applied Medical Sciences, Majmaah University, Majmaah 11952, Saudi Arabia
2
Health and Basic Sciences Research Center, Majmaah University, Majmaah 11952, Saudi Arabia
3
Department of Physiology, College of Medicine, Imam Abdulrahman Bin Faisal University, Dammam 31541, Saudi Arabia
4
Department of Neurology, College of Medicine, Imam Abdulrahman Bin Faisal University, Dammam 31952, Saudi Arabia
5
Department of Genetic Research, Institute for Research and Medical Consultations (IRMC), Imam Abdulrahman Bin Faisal University, Dammam 31441, Saudi Arabia
6
Department of Biomedical Sciences, Oregon State University, Corvallis, OR 97331, USA
7
Department of Epidemic Diseases Research, Institute for Research and Medical Consultations (IRMC), Imam Abdulrahman Bin Faisal University, Dammam 31441, Saudi Arabia
*
Author to whom correspondence should be addressed.
Brain Sci. 2022, 12(11), 1429; https://doi.org/10.3390/brainsci12111429
Submission received: 20 September 2022 / Revised: 13 October 2022 / Accepted: 20 October 2022 / Published: 24 October 2022
(This article belongs to the Special Issue The Effects of Neurotransmitters on Stem Cell Differentiation)

Abstract

:
Migraine, as the seventh most disabling neurological disease with 26.9% prevalence in Saudi females, lacks studies on identifying associated genes and pathways with migraines in the Arab population. This case control study aims to identify the migraine-associated novel genes and risk variants. More than 1900 Arab ancestry young female college students were screened: 103 fulfilled the ICHD-3 criteria for migraine and 20 cases confirmed in the neurology clinic were included for the study with age-matched healthy controls. DNA from blood samples were subjected to paired-end whole-exome sequencing. After quality control, 3365343 missense, frameshift, missense splice region variants and insertion–deletion (indels) polymorphisms were tested for association with migraine. Significant variants were validated using Sanger sequencing. A total of 17 (p-value 9.091 × 10−05) functional variants in 12 genes (RETNLB, SCAI, ADH4, ESPL1, CPT2, FLG, PPP4R1, SERPINB5, ZNF66, ETAA1, EXO1 and CPA6) were associated with higher migraine risk, including a stop-gained frameshift (-13-14*SX) variant in the gene RETNLB (rs5851607; p-value 3.446 × 10−06). Gene analysis revealed that half of the significant novel migraine risk genes were expressed in the temporal lobe (p-value 0.0058) of the cerebral cortex. This is the first study exploring the migraine risk of 17 functional variants in 12 genes among Saudi female migraineurs of Arab ancestry using whole-exome sequencing. Half of the significant genes were expressed in the temporal lobe, which expands migraine pathophysiology and early identification using biomarkers for research possibilities on personalised genetics.

1. Introduction

Migraine is a multifactorial chronic neurological disorder with a global prevalence between 14 to 21.7% [1]. It creates a burden on the economy as well as on the normal lives of individuals [2]. It is considered the seventh most disabling disease around the world as it causes 2.9% of years of life lost to disability (YLD) [3,4,5]. In Saudi Arabia, the prevalence of migraines in females is around 26.9% and is considered higher than the global average [6,7]. In addition, recent studies in this population revealed that high oestrogen levels could be involved in mediating the non-menstrual-related migraine among young Saudi females [8], and serum ApoE was reported as an excellent diagnostic marker for the same subjects with migraine in ictal or interictal phase [9]. However, genetic studies are deficient concerning migraines in this population [8].
Migraine is classified mainly into two forms: migraine with aura (MA) and migraine without aura (MO). It can also be classified into a chronic and episodic migraine. Hemiplegic migraine is another type of MA and is a severe and rare condition that affects one side of the body and causes temporary numbness [9]. It once was believed that migraine is a vascular disease; however, many recent types of research confirmed the involvement of multiple mechanisms related to the brain structures and processes [10]. Furthermore, accumulated evidence demonstrates a genetic background for the transmission of the disease. Many epidemiological studies showed the passage of migraine in families with a concordance rate of 1- to 2-fold in monozygotic twins in comparison to dizygotic twins [11]. Different approaches were used to determine the genetic factors and the DNA variants that might be responsible for migraine. Some studies focused on candidate-gene-associated studies (CGAS) to identify specific genetic markers associated with migraine but the negative result of these studies and the development of cost-effective techniques resulted in genome-wide association studies [12]. In these studies, hundreds of thousands of single-nucleotide polymorphisms (SNPs) were screened for their association with migraine. Studies from headache clinics of the Netherlands, Germany and Finland on a 10,747 people led to the identification of rs835740, a single, significant SNP [13]. A similar recent study identified 38 different SNPs and 7 loci, TRPM8, LRP1, FHL5, TSPAN2, ASTN2, near FGF6, and PHACTR1, but none of them were found related to migraine without aura [14]. However, there are no detailed studies on migraine-associated genes from the Arab population [8]. Focused studies on SNPs especially in females are needed, as they are the most affected patients of migraine. Identification of these migraine-associated SNPs needs collaboration from different areas of the world to collect sufficient information from populations with different ethnicities and draw the full picture of the genetic predisposition of migraine. Therefore, this study was designed to identify the associated gene variants and single-nucleotide polymorphisms with migraine in young Saudi females using whole-exome sequencing analysis.

2. Materials and Methods

The study protocol was approved by the Institutional Review Board of Imam Abdulrahman Bin Faisal University (IRB number: IRB-2021-01-250) and was conducted following the Declaration of Helsinki. More than 1900 Arab ancestry young female college students were screened, 103 fulfilled the ICHD-3 criteria for migraine and 20 cases confirmed in the neurology clinic were included for the study with age-matched controls. This is an exome-wide association study involving 40 participants: 20 controls (healthy subjects) and 20 cases (migraineurs) conducted in the female campus of the College of Medicine and the Institute for Research and Medical Consultations (IRMC) of Imam Abdulrahman bin Faisal University (IAU), Dammam, Saudi Arabia. Subjects were recruited during the year 2021 by convenience sampling and were all Saudi female college students with an age range of 18–30. For the cases to be included in the study, they had to be diagnosed with migraine by a neurologist and satisfy the criteria of the International Classification of Headache Disorders 3rd edition (ICHD). Controls were healthy female subjects of the same age group and with no complaints of headache.
All participants filled a written informed consent for their enrolment in the study, and then they were interviewed and asked to fill an electronic datasheet (Supplementary Material S1). The datasheet included the demographic characteristics (age, marital status, college level, height, weight, BMI). In addition, it included specific questions related to migraine such as frequency of attacks/month, the severity of the attack (using the visual scale 1–10), associated symptoms, presence of triggers (stress, lack of sleep, missed meal or fasting, physical activity, noise, smell, strong lights, fluctuation of weather or temperature, food, relation to the menstrual cycle), type of migraine (with aura, without aura), family history of migraine, history (other chronic diseases) and use of medications (for migraine, for other diseases).

2.1. Whole-Exome Sequencing and Statistical Data Analysis

Blood samples were obtained from the study participants in EDTA-vacutainers. Deoxyribonucleic acid was extracted using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany), DNA purity was checked using nanodrop and concentration was determined using a qubit fluorometer. DNA integrity was tested using agarose gel electrophoresis. Paired-end whole-exome sequencing was conducted for all the samples and subjected to quality screening. The qualified DNA sample was randomly fragmented (150 to 200 bp) and adapters were ligated to both ends and purified by the AMPure beads to excise about 200 bp fragments. Fragmented DNA was subjected to ligation-mediated polymerase chain reaction, followed by the SureSelect Library for enrichment. Captured ligation-mediated polymerase chain reaction products were estimated using Bioanalyzer and loaded on Hiseq platform for high-throughput sequencing. Sample depth means ≥ 7.5, sample variant call rate ≥ 0.5 and sample genotype quality mean ≥ 28 were considered for good quality sample. The conditions such as Phred score quality ≥ 30, raw read depth ≥ 10 and mapping quality ≥ 30 were applied for filtering the good quality variants. The Hail standard (python package) was used for the entire genome-wide association study pipeline analysis as well as for the quality check of the samples, variants and genotypes. Manhattan plot and QQ plot of the association of SNPs with migraine as statistical significance in terms of p-values on a genomic scale were constructed. Gene-level linkage disequilibrium analysis of the SNPs in the significant gene was conducted using Haploview 4.2.

2.2. Statistical Analysis

Data were analysed using a statistical package for the social sciences (SPSS) software version 21. Student’s t-test was used for identifying the significant difference among the demographic characteristics between migraineurs and the controls. The most significant top 10 genes were analysed using the gene functional classification tool DAVID to identify the significance in the site of expression under the category of GNF_U133A_QUARTILE (p-value < 0.05) [15]. The functional annotation of the top 50 highly associated genes (with the lowest p-value < 0.00023) was performed using the Uniprot database. Further, the GO and pathway enrichment were carried out using enrichR server, and pathway involvement of the genes was performed by KEGG search and colour pathway server. All the significant markers identified from the exome-wide association study and genes were selected for the expression profile analysis in the brain and related tissues using DAVID. Brain and related tissue-expressed genes were separated and analysed for KEGG pathway enrichment.

3. Results

3.1. Study Population

Study participants (Table 1) were drawn from Arab ancestries. The demographic characteristics of the 40 selected Saudi Arabian subjects including 20 migraine patients and controls (n = 20) are presented in Table 1. The clinical characteristics and the frequency of the precipitating factors of headache attack in the migraineurs of the study are presented in Table 2.

3.2. Single-Variant Analysis

After quality controls for the variants were obtained in the whole-exome sequencing, 3,365,343 variants were satisfied for further exome-wide association analyses. Our study highlights the added influence of considering the functional variants such as missense variants, frameshift variant and missense splice region variant in the analysis: p-values < 0.00001 were set from the exome-wide association analysis (Table 3) to prioritise the top migraine-associated variants in Arab ancestry. The entire list of migraine-associated (p-value < 0.00001) variants identified through exome sequencing are shown in Supplementary Table S1. Seventeen variants were found to be the most significant (9.091 × 10−5) functional variants distributed among twelve genes (RETNLB, SCAI, ADH4, ESPL1, CPT2, FLG, PPP4R1, SERPINB5, ZNF66, ETAA1, EXO1 and CPA6) (Table 3 and Figure 1). The stop-gained frameshift (-13-14*SX) variant in the gene RETNLB was found to be the most significant functional variant (rs5851607; p-value = 3.446 × 10−06).

3.3. Gene Analysis

The most significant top 12 genes were analysed using the gene functional classification tool DAVID to identify the significance in the site of expression. The analysis revealed that 6 out of 12 genes (Table 3) were significantly expressed in the temporal lobe (p-value = 0.00582) (Figure 2). The functional variants with a p-value between 9.091 × 10−05 to 0.05 on the genes identified as expressed in the temporal lobe revealed the significance of FLG gene with 37 functional variants (Table 4). The significant SNPs observed in the FLG gene, and their amino acid position, are presented in Figure 3. The gene-level linkage disequilibrium analysis (Haploview 4.2) of the 53 SNPs in the FLG gene reveals significant association for all these SNPs (Chi-Square = 8.556; p-value = 0.0034) (Figure 3). The most significant three markers haplotype rs3126075G, rs7532285T and rs7540123G (Chi-Square = 7.64; p-value = 0.0057) appear to associate significantly with migraine, while the opposite alleles rs3126075C, rs7532285C and rs7540123C (Chi-Square = 3.81; p-value = 0.0407) are the protective type of haplotype. The presence of variants at FLG was confirmed using Sanger sequencing with the designed primers (forward primer, FLGF: 5′ CCTCTACCAGGTGAGCACTCATGAACAGTCTG 3′, and reverse primer, FLGR: 5′ TCTCTGACTGCAGATGAAGCTTGTCCGTGCC 3′). Sanger sequencing revealed additional variants in the gene and the need to check their association with migraine (Figure 4).

3.4. Pathways Analysis

Gene ontology and pathway enrichment of the top 50 genes revealed that the migraine subjects are moderately significant for the organic hydroxy compound catabolic process (GO:1901616; p-value = 9.12893 × 10−05; adjusted p-value = 0.026108731) and quinone metabolic process (GO:1901661; p-value = 0.000272123; adjusted p-value = 0.038913602) (Table 5; Supplementary Table S2). All the significant markers having genes were checked for the expression profile in the brain-related tissues using DAVID. A total of 6305 genes were present for DAVID among 11958 genes presented for expression analysis using DAVID. A total of 1349 genes were separated based on the expression in the brain and related tissues were analysed for KEGG pathway enrichment; the results revealed that 34 genes were found to be significantly (term p-value = 8.071 × 10−12; adjusted p-value = 2.364 × 10−09) associated with systemic lupus erythematosus (Table 6; Supplementary Table S3; Figure S3). Furthermore, the pathways such as focal adhesion, ECM-receptor interaction, human papillomavirus infection, alcoholism, pathways in cancer, pi3k-akt signalling pathway and cholesterol metabolism are significantly associated with migraine in the Saudis (adjusted p-value ≤ 2.192 × 10−05) (Table 5).

4. Discussion

Genetics play a significant part in migraines, in addition to other factors [16]. However, migraine genetic predisposition does not follow a direct Mendelian pattern. The common form of migraine is most probably polygenic and involves multiple variants at several genetic loci that possibly interact with multiple environmental factors. Exome-wide association study is the most successful method to identify the genes involved in a disease. In this methodology, cohorts of migraine cases and controls are explored for any differences in allele frequencies of single-nucleotide polymorphisms (SNPs) to identify genetic risk factors. There is no single genetic variant that can explain migraine heterogeneity across populations. We performed the first exome-wide association study of migraine in Arab ancestry. Through exome sequencing, we identified an entire list of migraine-associated (p-value < 0.00001) variants and prioritised 17 as the most significant (9.091 × 10−05) functional variants distributed among 12 genes (RETNLB, SCAI, ADH4, ESPL1, CPT2, FLG, PPP4R1, SERPINB5, ZNF66, ETAA1, EXO1 and CPA6) in the Saudi females suffering from migraine. All of these were novel and have not been documented in earlier studies involving other populations such as Europeans [17] and Chinese [18].
The stop-gained frameshift variant in the gene RETNLB is the most significant functional variant, rs5851607; this gene encodes a bactericidal protein, Resistin-like molecule β (RELMβ), which is released from colonic cells to destroy Gram-negative bacteria. Migraine may be associated with diseases such as irritable bowel syndrome (IBS), inflammatory bowel syndrome, and celiac disease [19]. The ADH4 gene encodes the alcohol dehydrogenase enzyme and variations in this gene are associated with alcohol dependence [20]. The CPA6 gene encodes the Carboxypeptidase A6 enzyme, and its mutations can predispose to various types of epilepsy [21]. As shown in Figure 2 and Table 3, a total of 36 functional variants were found to be significant in the gene FLG. This gene encodes a protein called profilaggrin present in the epidermis of the skin. This protein is important for the skin’s barrier function. Functional variations in this gene can cause sensitisation or atopic dermatitis [22] and might be the underlying mechanism of migraine-associated allodynia [23]. The ETAA1 gene encodes the protein Ewing tumour-associated antigen 1. This protein functions as a DNA replication stress response protein [24]. The CPT2 gene encodes the carnitine palmitoyltransferase 2 enzyme, which is essential for fatty acid oxidation. The ESPL1 gene codes a protease separase/separin which causes separation of sister chromatids in mitosis [25]. The SERPINB5 gene encodes a protein Maspin, a tumour suppressor that binds directly to extracellular matrix components and inhibits tumour-induced angiogenesis, invasion and metastatic spread [26]. SCAI encodes a protein that suppresses cancer cell invasion [27]. The significant association of multiple genes to migraine might help to explain the wide spectrum of migraine phenotypes.
The morphological changes in the temporal lobe were reported to be associated with migraines [27]. Recently, a reduction in grey matter volume in the temporal lobe was observed in migraine patients [28]. Gene analysis revealed that 6 of the significant 12 novel migraine risk genes were expressed in the temporal lobe of the cerebral cortex. The present study adds a molecular insight into the observations on the temporal lobe [29] and migraine-associated genes [30,31] and opens new avenues for migraine research. The current study will help in power calculations in the future and will provide potential loci to look for in replication studies. This may facilitate a thorough understanding of migraine pathophysiology and its underlying molecular mechanism, and open avenues for more precise diagnosis and therapeutic strategies targeting migraine patients of Arab ancestry. The study may also help in the polygenic risk scoring of the patients. The number of study subjects is one of the notable limitations in the study.

5. Conclusions

Our study is the first one exploring migraine genetic variations in Arab ancestry. Seventeen significant functional variants including a stop gained in twelve genes are the migraine risk variants in Arab ancestry. Half of the significant novel migraine risk genes are expressed in the temporal lobe of the brain.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/brainsci12111429/s1, Figure S1: The p53 signaling pathway with the significant genes from the study; Figure S2: The Fatty acid degradation with the significant genes is the study; Figure S3: Systemic lupus erythematosus (hsa05322) pathway and significant genes (highlighted) identified from the study. Table S1: List of migraine associated (p < 0.00001) variants identified through exome sequencing; Table S2: Gene ontology pathway analysis of the top 50 genes; Table S3: KEGG pathway enrichment from the 1349 genes based on the expression in brain related tissues associated in the GWAS analysis.

Author Contributions

J.K., L.A.A., S.A. and J.F.B. conceived and designed the research and analysed the experiments. L.A.A., J.K., M.A. (Majed Aloqaily), S.A. and J.F.B. performed and analysed the experiments. S.A. and J.F.B. performed the whole-genome analysis. L.A.A., A.A.S., N.R., R.L. and M.A. (Majed Alabdali) performed clinical analysis. J.K., L.A.A. and J.F.B. wrote the paper with the contributions of A.A.S., N.R., R.L., M.A. (Majed Alabdali), A.A., M.A. (Majed Aloqaily), S.A., S.J., S.B., J.K., L.A.A., A.A.S., N.R., R.L., M.A. (Majed Alabdali), A.A., M.A. (Majed Aloqaily), S.A., S.J., S.B., and J.F.B. reviewed and approved the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The author extends their appreciation to the deputyship for Research and Innovation, Ministry of Education in Saudi Arabia for funded this research work through the project number (IFP-2020-46).

Institutional Review Board Statement

The study protocol was approved by the Institutional Review Board of Imam Abdulrahman Bin Faisal University (IRB number: IRB-2021-01-250) and was conducted following the Declaration of Helsinki.

Informed Consent Statement

Informed consentwas obtained from all subjects involved in the study.

Data Availability Statement

All data will be available on reasonable request from the corresponding author.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number (IFP-2020-46). The authors thank the Dean, Institute for Research and Medical Consultations (IRMC), Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia, for her continuous support and encouragement. The authors thank Ranilo, M. Tumbaga, and Horace T. Pacifico for their assistance. The preprint of the paper is available in the Research Square preprint platform “https://www.researchsquare.com/article/rs-885676/v1” (accessed on 16 September 2021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Woldeamanuel, Y.W.; Cowan, R.P. Migraine affects 1 in 10 people worldwide featuring recent rise: A systematic review and meta-analysis of community-based studies involving 6 million participants. J. Neurol. Sci. 2017, 372, 307–315. [Google Scholar] [CrossRef] [PubMed]
  2. Khan, J.; Al Asoom, L.I.; Al Sunni, A.; Rafique, N.; Latif, R.; Al Saif, S.; Almandil, N.B.; Almohazey, D.; AbdulAzeez, S.; Borgio, J.F. Genetics, pathophysiology, diagnosis, treatment, management, and prevention of migraine. Biomed. Pharmacother. 2021, 139, 111557. [Google Scholar] [CrossRef] [PubMed]
  3. Pietrobon, D. Familial hemiplegic migraine. Neurotherapeutics 2007, 4, 274–284. [Google Scholar] [CrossRef] [PubMed]
  4. Vos, T.; Flaxman, A.D.; Naghavi, M.; Lozano, R.; Michaud, C.; Ezzati, M.; Shibuya, K.; Salomon, J.A.; Abdalla, S.; Aboyans, V.; et al. Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990–2010: A systematic analysis for the Global Burden of Disease Study 2010. Lancet 2012, 380, 2163–2196. [Google Scholar] [CrossRef]
  5. Al Asoom, L.; Khan, J.; Al Sunni, A.; Rafique, N.; Latif, R.; Alabdali, M.; AbdulAzeez, S.; Borgio, J.F. A Pilot Mitochondrial Genome-Wide Association on Migraine among Saudi Arabians. Int. J. Gen. Med. 2022, 15, 6249–6258. [Google Scholar] [CrossRef]
  6. Charles, A. The pathophysiology of migraine: Implications for clinical management. Lancet Neurol. 2018, 17, 174–182. [Google Scholar] [CrossRef]
  7. Almalki, Z.A.; Alzhrani, M.A.G.; Altowairqi, A.T.; Aljawi, Y.A.; Fallatah, S.A.; Assaedi, L.M.; Aljawi, M.A.; Alqusair, S.A. Prevalence of Migraine Headache in Taif City, Saudi Arabia. J. Clin. Med. Res. 2018, 10, 125–133. [Google Scholar] [CrossRef] [Green Version]
  8. Borgio, J.F. Heterogeneity in biomarkers, mitogenome and genetic disorders of Arab population with special emphasis on large-scale whole-exome sequencing. Arch. Med. Sci. 2021. [Google Scholar] [CrossRef]
  9. Latif, R.; Rafique, N.; Al Asoom, L.; Alsunni, A.A.; Salem, A.M.; AlNoaimi, D.A.; AlHarbi, W.M.; Shaikh, A.Y.; AlRaddadi, S.A.; AlMuhaish, L.A.; et al. Diagnostic Accuracy of Serum Calcitonin Gene-Related Peptide and Apolipoprotein E in Migraine: A Preliminary Study. Int. J. Gen. Med. 2021, 14, 851–856. [Google Scholar] [CrossRef]
  10. Al Asoom, L.I.; Alajmi, M.S.; Alsudairi, R.R.; AlShamlan, A.A.; Almomaten, A.A.; Alqarni, A.A.; Alshammari, M.H.; Rafique, N.; Latif, R.; Alsunni, A.A.; et al. Association between sex hormones and migraine in young Saudi females. Saudi Med. J. 2021, 42, 793–797. [Google Scholar] [CrossRef]
  11. Al Jumah, M.; Al Khathaami, A.M.; Kojan, S.; Hussain, M.; Thomas, H.; Steiner, T.J. The prevalence of primary headache disorders in Saudi Arabia: A cross-sectional population-based study. J. Headache Pain 2020, 21, 11. [Google Scholar] [CrossRef] [PubMed]
  12. Dodick, D.W. A Phase-by-Phase Review of Migraine Pathophysiology. Headache J. Head Face Pain 2018, 58, 4–16. [Google Scholar] [CrossRef] [Green Version]
  13. Lipton, R.B.; Bigal, M.E. The epidemiology of migraine. Am. J. Med. Suppl. 2005, 118, 3–10. [Google Scholar] [CrossRef]
  14. Samuel Appadurai Alexander, R. Study of Migraine Headache Spectrum in patients treated in the Tertiary Care Hospital. Ph.D. Thesis, Stanley Medical College, Chennai, India, 2014. [Google Scholar]
  15. Huang, D.W.; Sherman, B.T.; Tan, Q.; Collins, J.R.; Alvord, W.G.; Roayaei, J.; Stephens, R.; Baseler, M.W.; Lane, H.C.; Lempicki, R.A. The DAVID Gene Functional Classification Tool: A novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol. 2007, 8, R183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Sutherland, H.G.; Griffiths, L.R. Genetics of Migraine: Insights into the Molecular Basis of Migraine Disorders. Headache J. Head Face Pain 2017, 57, 537–569. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Gormley, P.; Anttila, V.; Winsvold, B.S.; Palta, P.; Esko, T.; Pers, T.H.; Farh, K.-H.; Cuenca-Leon, E.; Muona, M.; Palotie, A.; et al. Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine. Nat. Genet. 2016, 48, 856–866. [Google Scholar] [CrossRef] [Green Version]
  18. Chang, X.; Pellegrino, R.; Garifallou, J.; March, M.; Snyder, J.; Mentch, F.; Li, J.; Hou, C.; Liu, Y.; Sleiman, P.M.A.; et al. Common variants at 5q33.1 predispose to migraine in African-American children. J. Med. Genet. 2018, 55, 831–836. [Google Scholar] [CrossRef]
  19. Jiang, Z.; Zhao, L.; Zhang, X.; Zhang, W.; Feng, Y.; Li, T. Common variants in KCNK5 and FHL5 genes contributed to the susceptibility of migraine without aura in Han Chinese population. Sci. Rep. 2021, 11, 6807. [Google Scholar] [CrossRef]
  20. Van Hemert, S.; Breedveld, A.C.; Rovers, J.M.; Vermeiden, J.P.; Witteman, B.J.; Smits, M.G.; de Roos, N.M. Migraine associated with gastrointestinal disorders: Review of the literature and clinical implications. Front. Neurol. 2014, 5, 241. [Google Scholar] [CrossRef] [Green Version]
  21. Luo, X.; Kranzler, H.R.; Zuo, L.; Lappalainen, J.; Yang, B.-Z.; Gelernter, J. ADH4 Gene Variation is Associated with Alcohol Dependence and Drug Dependence in European Americans: Results from HWD Tests and Case–Control Association Studies. Neuropsychopharmacology 2005, 31, 1085–1095. [Google Scholar] [CrossRef]
  22. Sapio, M.R.; Vessaz, M.; Thomas, P.; Genton, P.; Fricker, L.D.; Salzmann, A. Novel Carboxypeptidase A6 (CPA6) Mutations Identified in Patients with Juvenile Myoclonic and Generalized Epilepsy. PLoS ONE 2015, 10, e0123180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Weidinger, S.; Illig, T.; Baurecht, H.; Irvine, A.D.; Rodriguez, E.; Diaz-Lacava, A.; Klopp, N.; Wagenpfeil, S.; Zhao, Y.; Liao, H.; et al. Loss-of-function variations within the filaggrin gene predispose for atopic dermatitis with allergic sensitizations. J. Allergy Clin. Immunol. 2006, 118, 214–219. [Google Scholar] [CrossRef] [PubMed]
  24. Lipton, R.B.; Bigal, M.E.; Ashina, S.; Burstein, R.; Silberstein, S.; Reed, M.L.; Ma, D.S.; Stewart, W.F. American Migraine Prevalence Prevention Advisory Group Cutaneous Allodynia in the Migraine Population. Ann. Neurol. 2008, 63, 148–158. [Google Scholar] [CrossRef] [PubMed]
  25. Lee, Y.-C.; Zhou, Q.; Chen, J.; Yuan, J. RPA-Binding Protein ETAA1 Is an ATR Activator Involved in DNA Replication Stress Response. Curr. Biol. 2016, 26, 3257–3268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Wirth, K.G.; Wutz, G.; Kudo, N.; Desdouets, C.; Zetterberg, A.; Taghybeeglu, S.; Seznec, J.; Ducos, G.M.; Ricci, R.; Firnberg, N.; et al. Separase: A universal trigger for sister chromatid disjunction but not chromosome cycle progression. J. Cell Biol. 2006, 172, 847–860. [Google Scholar] [CrossRef]
  27. Sheng, S.; Carey, J.; Seftor, E.A.; Dias, L.; Hendrix, M.J.; Sager, R. Maspin acts at the cell membrane to inhibit invasion and motility of mammary and prostatic cancer cells. Proc. Natl. Acad. Sci. USA 1996, 93, 11669–11674. [Google Scholar] [CrossRef] [Green Version]
  28. Brandt, D.T.; Baarlink, C.; Kitzing, T.M.; Kremmer, E.; Ivaska, J.; Nollau, P.; Grosse, R. SCAI acts as a suppressor of cancer cell invasion through the transcriptional control of β1-integrin. Nat. Cell Biol. 2009, 11, 557–568. [Google Scholar] [CrossRef] [PubMed]
  29. Kursun, O.; Yemisci, M.; van den Maagdenberg, A.M.J.M.; Karatas, H. Migraine and neuroinflammation: The inflammasome perspective. J. Headache Pain 2021, 22, 55. [Google Scholar] [CrossRef] [PubMed]
  30. Schwedt, T.J.; Berisha, V.; Chong, C.D. Temporal Lobe Cortical Thickness Correlations Differentiate the Migraine Brain from the Healthy Brain. PLoS ONE 2015, 10, e0116687. [Google Scholar] [CrossRef]
  31. Al Asoom, L.I.; Khan, J.; Al Sunni, A.; Rafique, N.; Latif, R.; Alabdali, M.M.; Alhariri, A.; Aloqaily, M.; AbdulAzeez, S.; Borgio, J.F. Whole-Exome Sequencing Reveals Migraine Associated Novel Functional Variants. 2021. Available online: https://doi.org/10.21203/rs.3.rs-885676/v1 (accessed on 1 July 2022).
Figure 1. Manhattan plot (A) and QQ plot (B) of the association of SNPs with migraine as statistical significance in terms of p-values on a genomic scale. SNP numbers with blue colour indicate the highly associated (p < 3.44623 × 10−06) stop-gained frameshift variant. SNP numbers with red colour indicate the highly associated (p < 0.00001) missense variants. SNP numbers with ash colour indicate the highly associated (p < 0.00001) intronic and 5 prime variants. Chr15: rs8037668 [‘T’, ‘C’] (p-value 2.7135 × 10−07), LOC400464 intron variant, non-coding transcript variant; and Chr 2: rs772174 [‘A’, ‘G’], (p-value 5.2006 × 10−07), ITPRIPL1 protein-coding 5 prime UTR variant showed the highest p-value. (C) Consequences of coding variants. Full list of variants (p-value < 0.00001) is listed in Table S1.
Figure 1. Manhattan plot (A) and QQ plot (B) of the association of SNPs with migraine as statistical significance in terms of p-values on a genomic scale. SNP numbers with blue colour indicate the highly associated (p < 3.44623 × 10−06) stop-gained frameshift variant. SNP numbers with red colour indicate the highly associated (p < 0.00001) missense variants. SNP numbers with ash colour indicate the highly associated (p < 0.00001) intronic and 5 prime variants. Chr15: rs8037668 [‘T’, ‘C’] (p-value 2.7135 × 10−07), LOC400464 intron variant, non-coding transcript variant; and Chr 2: rs772174 [‘A’, ‘G’], (p-value 5.2006 × 10−07), ITPRIPL1 protein-coding 5 prime UTR variant showed the highest p-value. (C) Consequences of coding variants. Full list of variants (p-value < 0.00001) is listed in Table S1.
Brainsci 12 01429 g001
Figure 2. Significant functional variants of genes expressed in the temporal lobe (p-value = 0.00582) for migraine in Saudis. The most significant genes (top 12 as listed in Table 3) associated functionally were used as input to identify the expression nature of them using the gene functional classification tool DAVID with the ‘GNF U133A QUARTILE’. The most significant (p-value < 9.091 × 10−05) functional variants are presented in the middle path. The functional variants with a p-value between 9.091 × 10−05 to 0.05 on the genes identified as expressed in the temporal lobe are also presented. A total of 37 functional variants are found to be significant in the gene FLG (full list presented in Table 4).
Figure 2. Significant functional variants of genes expressed in the temporal lobe (p-value = 0.00582) for migraine in Saudis. The most significant genes (top 12 as listed in Table 3) associated functionally were used as input to identify the expression nature of them using the gene functional classification tool DAVID with the ‘GNF U133A QUARTILE’. The most significant (p-value < 9.091 × 10−05) functional variants are presented in the middle path. The functional variants with a p-value between 9.091 × 10−05 to 0.05 on the genes identified as expressed in the temporal lobe are also presented. A total of 37 functional variants are found to be significant in the gene FLG (full list presented in Table 4).
Brainsci 12 01429 g002
Figure 3. Top: variants observed in FLG gene. Pink lollipop indicates silent mutation; green lollipop indicates missense variations. Amino acid substitution in blue colour indicates the most significant variant. Box coloured green indicates the calcium-binding domain. Bottom: Haplomap of gene-level linkage disequilibrium analysis of the SNPs in the FLG gene. Pink line: the most significant three markers rs3126075G, rs7532285T and rs7540123G.
Figure 3. Top: variants observed in FLG gene. Pink lollipop indicates silent mutation; green lollipop indicates missense variations. Amino acid substitution in blue colour indicates the most significant variant. Box coloured green indicates the calcium-binding domain. Bottom: Haplomap of gene-level linkage disequilibrium analysis of the SNPs in the FLG gene. Pink line: the most significant three markers rs3126075G, rs7532285T and rs7540123G.
Brainsci 12 01429 g003
Figure 4. Electropherogram of FLG gene with (A) NM_002016.2:c.1676A>G p.H559R and (B) NM_002016.2:c.1741A>T p.T581S. The blue arrow indicates mutant. The green arrow indicates wild type.
Figure 4. Electropherogram of FLG gene with (A) NM_002016.2:c.1676A>G p.H559R and (B) NM_002016.2:c.1741A>T p.T581S. The blue arrow indicates mutant. The green arrow indicates wild type.
Brainsci 12 01429 g004
Table 1. Demographic characteristics of the migraineurs and the controls.
Table 1. Demographic characteristics of the migraineurs and the controls.
Control n = 20Migraineurs n = 20p-Value *
Age (years)21.86 ± 1.7522.10 ± 3.630.818
Bodyweight (kg)56.21 ± 14.0263.00 ± 12.580.155
BMI21.48 ± 5.3924.66 ± 5.180.095
* Student’s t-test.
Table 2. Clinical characteristics of the migraineurs of the study and the frequency of the precipitating factors.
Table 2. Clinical characteristics of the migraineurs of the study and the frequency of the precipitating factors.
VariableDescription
Type of migraineWith aura 11 (57.9%), without aura 9 (42.1%)
Family historyYes 11 (55%), no 9 (45%)
Use of medicationsNo medication 7 (35%), pain killer 10 (50%), prophylaxis 3 (15%)
Number of attacks/monthsRanges from 2–28 attacks/month with 2 attacks/month the most common
Duration of the attackRanges from 4 h to >72 h with the most common 12–24 h
The severity of the attack (Visual scale 1–10)Ranges from 7–10
Precipitating factors
YesNo
Sleep disturbances17 (85%)3 (15%)
Stress16 (80%)4 (20%)
Bright light18 (90%)2 (10%)
Excessive noise15 (75%)5 (25%)
Strong smells12 (60%)8 (40%)
Weather changes8 (40%)12 (60%)
Skipped meal10 (50%)10 (50%)
Physical exertion5 (25%)15 (75%)
Certain types of food5 (25%)15 (75%)
Coffee2 (10%)18 (90%)
Table 3. List of migraine-associated (p-value < 0.00001) exonic functional variants identified through exome sequencing.
Table 3. List of migraine-associated (p-value < 0.00001) exonic functional variants identified through exome sequencing.
Locus. ContigVariation IDGene SYMBOLMAFAllelesAmino AcidsCodonsp-Value
1chr3rs5851607 *RETNLB0.18[‘G’, ‘GGGGGATTA’]-13-14*SX-/TAATCCCC3.446 × 10−06
2chr9rs589292SCAI0.2[‘C’, ‘T’]A37TGct/Act1.169 × 10−05
3chr4rs1126671ADH40.26[‘T’, ‘C’]I309VAtt/Gtt1.575 × 10−05
4chr4rs1126673 $ADH40.26[‘C’, ‘T’]V393IGtc/Atc1.575 × 10−05
5chr12rs6580942ESPL10.3[‘C’, ‘A’]A25DgCc/gAc3.413 × 10−05
6chr1rs1799821CPT20.28[‘G’, ‘A’]V368IGtc/Atc3.585 × 10−05
7chr1rs3126075FLG0.23[‘G’, ‘C’]T3579RaCg/aGg4.266 × 10−05
8chr18rs329003PPP4R10.28[‘T’, ‘C’]I381VAta/Gta4.500 × 10−05
9chr18rs2289520SERPINB50.14[‘G’, ‘C’]V187LGtc/Ctc4.556 × 10−05
10chr12rs56358776ESPL10.26[‘G’, ‘A’]R1561QcGg/cAg5.923 × 10−05
11chr19rs10413187ZNF660.14[‘C’, ‘A’]Q66KCag/Aag7.328 × 10−05
12chr19rs432839ZNF660.22[‘G’, ‘T’]C173FtGc/tTc7.320 × 10−05
13chr19rs383038ZNF660.14[‘T’, ‘C’]F188LTtt/Ctt7.328 × 10−05
14chr19rs370551ZNF660.14[‘A’, ‘G’]T420AAct/Gct7.328 × 10−05
15chr2rs61740794ETAA10.4[‘G’, ‘A’]E673KGaa/Aaa8.508 × 10−05
16chr1rs735943EXO10.42[‘A’, ‘G’]H354RcAt/cGt8.978 × 10−05
17chr8rs17343819 $CPA60.22[‘T’, ‘C’]N249SaAt/aGt9.091 × 10−05
* Stop-gained frameshift variant; $ missense splice region variant. Significant (p-value < 0.00001) functional and other variants identified through exome sequencing are listed in Table S1. MAF: minor allele frequency.
Table 4. List of significant migraine-associated missense variants of FLG gene in chromosome 1.
Table 4. List of significant migraine-associated missense variants of FLG gene in chromosome 1.
SNPAllelesMAFAmino AcidsProtein PositionCodonsp-Value
1rs3126075[‘G’, ‘C’]0.23T/R3579aCg/aGg4.27 × 10−05
2rs7532285[‘T’, ‘C’]0.05Q/R3568cAg/cGg0.000148
3rs7540123[‘G’, ‘C’]0.05Q/E3568Cag/Gag0.000148
4rs2065955[‘C’, ‘G’]0.3G/A3436gGa/gCa0.000592
5rs3126079[‘G’, ‘T’]0.3H/Q1961caC/caA0.000592
6rs58001094[‘G’, ‘C’]0.3A/G1167gCa/gGa0.000592
7rs11582087[‘T’, ‘G’]0.04S/R2836Agt/Cgt0.00092
8rs71625202[‘C’, ‘G’]0.08S/T2366aGt/aCt0.00092
9rs139476473[‘C’, ‘T’]0.04D/N2339Gac/Aac0.00092
10rs2065957[‘A’, ‘C’]0.17V/G3179gTg/gGg0.001166
11rs12083389[‘C’, ‘G’]0.15E/D3593gaG/gaC0.002327
12rs3126072[‘C’, ‘T’]0.22G/R2545Gga/Aga0.003517
13rs3126074[‘G’, ‘C’]0.22H/Q2507caC/caG0.003517
14rs2011331[‘T’, ‘C’]0.22T/A454Aca/Gca0.003517
15rs66954353[‘T’, ‘G’]0.13K/Q2192Aaa/Caa0.003521
16rs2184953[‘A’, ‘G’]0.32Y/H2194Tat/Cat0.004922
17rs140376327[‘G’, ‘A’]0.04R/W2430Cgg/Tgg0.010229
18rs12135040[‘C’, ‘G’]0.04G/R1936Ggg/Cgg0.010229
19rs138721961[‘C’, ‘T’]0.04R/H402cGc/cAc0.010229
20rs78179835[‘C’, ‘G’]0.08E/D2297gaG/gaC0.014364
21rs113544881[‘A’, ‘T’]0.07L/H1943cTt/cAt0.016921
22rs74129452[‘T’, ‘G’]0.18Q/H2154caA/caC0.018862
23rs7512553[‘A’, ‘G’]0.18Y/H2119Tat/Cat0.018862
24rs7522925[‘G’, ‘A’]0.18A/V2108gCg/gTg0.018862
25rs7512857[‘A’, ‘C’]0.18S/A2020Tca/Gca0.018862
26rs12407807[‘C’, ‘T’]0.16R/H1684cGc/cAc0.023383
27rs75235053[‘C’, ‘G’]0.06S/T3662aGt/aCt0.026472
28rs199888588[‘A’, ‘G’]0.1W/R962Tgg/Cgg0.027008
29rs74129455[‘T’, ‘G’]0.1K/Q2064Aaa/Caa0.031164
30rs149817134[‘G’, ‘T’]0.04H/N1880Cac/Aac0.03179
31rs55650366[‘A’, ‘G’]0.14L/S2481tTg/tCg0.035772
32rs71625200[‘T’, ‘C’]0.14K/E2444Aag/Gag0.035772
33rs71625201[‘C’, ‘G’]0.14E/Q2398Gag/Cag0.035772
34rs11581433[‘T’, ‘C’]0.14R/G1376Aga/Gga0.035772
35rs74129461[‘C’, ‘T’]0.14E/K755Gaa/Aaa0.035772
36rs11584340[‘G’, ‘A’]0.14P/S478Cct/Tct0.035772
37rs41267154[‘C’, ‘A’]0.14G/V332gGc/gTc0.035772
MAF: minor allele frequency.
Table 5. KEGG pathway enrichment from the top 50 genes associated with the exome wide association study analysis.
Table 5. KEGG pathway enrichment from the top 50 genes associated with the exome wide association study analysis.
Termp-ValueNumber of Genes from Top 50 ListTotal Genes InvolvedGenes from the Top 50 List
p53 signalling pathway0.0007765713 *72TNFRSF10B; SERPINB5; ATR
Fatty acid degradation0.0054187662 $44ADH4; CPT2
Metabolism of xenobiotics by cytochrome P4500.014752518274ADH4; AKR1C1
Cell cycle0.0384970722124ESPL1; ATR
Mannose type O-glycan biosynthesis0.055975781123FKRP
Mismatch repair0.055975781123EXO1
Tyrosine metabolism0.086243722136ADH4
Human T-cell leukaemia virus 1 infection0.103954352219ESPL1; ATR
Vasopressin-regulated water reabsorption0.104395242144AQP4
Fanconi anaemia pathway0.126588324154ATR
Steroid hormone biosynthesis0.139644523160AKR1C1
Cytosolic DNA-sensing pathway0.146100703163CCL4L2
Retinol metabolism0.154635136167ADH4
Glycolysis/Gluconeogenesis0.156755649168ADH4
* The p53 signalling pathway with the significant genes is presented in Figure S1. $ The fatty acid degradation with the significant genes is presented in Figure S2.
Table 6. KEGG pathway enrichment from the 1349 genes based on the expression in brain-related tissues associated with the exome-wide association analysis.
Table 6. KEGG pathway enrichment from the 1349 genes based on the expression in brain-related tissues associated with the exome-wide association analysis.
Termp-ValueAdjusted p-ValueGenes from 1349 List *Total Genes
Systemic lupus erythematosus8.071 × 10−122.364 × 10−0934133
Focal adhesion2.665 × 10−113.905 × 10−0942199
ECM-receptor interaction4.693 × 10−104.583 × 10−082482
Human papillomavirus infection3.167 × 10−092.319 × 10−0753330
Alcoholism1.142 × 10−086.693 × 10−0735180
Pathways in cancer2.191 × 10−081.069 × 10−0671530
PI3K-Akt signalling pathway9.237 × 10−083.866 × 10−0652354
Cholesterol metabolism5.986 × 10−072.192 × 10−051550
* The list of significant genes is presented in Supplementary Table S3.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khan, J.; Al Asoom, L.; Al Sunni, A.; Rafique, N.; Latif, R.; Alabdali, M.; Alhariri, A.; Aloqaily, M.; AbdulAzeez, S.; Jahan, S.; et al. Whole-Exome Sequencing Reveals Migraine-Associated Novel Functional Variants in Arab Ancestry Females: A Pilot Study. Brain Sci. 2022, 12, 1429. https://doi.org/10.3390/brainsci12111429

AMA Style

Khan J, Al Asoom L, Al Sunni A, Rafique N, Latif R, Alabdali M, Alhariri A, Aloqaily M, AbdulAzeez S, Jahan S, et al. Whole-Exome Sequencing Reveals Migraine-Associated Novel Functional Variants in Arab Ancestry Females: A Pilot Study. Brain Sciences. 2022; 12(11):1429. https://doi.org/10.3390/brainsci12111429

Chicago/Turabian Style

Khan, Johra, Lubna Al Asoom, Ahmad Al Sunni, Nazish Rafique, Rabia Latif, Majed Alabdali, Azhar Alhariri, Majed Aloqaily, Sayed AbdulAzeez, Sadaf Jahan, and et al. 2022. "Whole-Exome Sequencing Reveals Migraine-Associated Novel Functional Variants in Arab Ancestry Females: A Pilot Study" Brain Sciences 12, no. 11: 1429. https://doi.org/10.3390/brainsci12111429

APA Style

Khan, J., Al Asoom, L., Al Sunni, A., Rafique, N., Latif, R., Alabdali, M., Alhariri, A., Aloqaily, M., AbdulAzeez, S., Jahan, S., Banawas, S., & Borgio, J. F. (2022). Whole-Exome Sequencing Reveals Migraine-Associated Novel Functional Variants in Arab Ancestry Females: A Pilot Study. Brain Sciences, 12(11), 1429. https://doi.org/10.3390/brainsci12111429

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