Genome-wide association study reveals genetic risk factors for trigeminal neuralgia
===================================================================================

* Andrew T. Hale
* Jing He
* Oluwatoyin Akinnusotu
* Rebecca L. Sale
* Janey Wang
* Lisa Bastarache
* Eric R. Gamazon

## Abstract

**Background** While many clinical risk factors of trigeminal neuralgia (TN) have been identified, the genetic basis of TN is largely unknown. Here, we perform the first genome-wide association study (GWAS) for TN using three independent DNA biobanks – BioVU, the UK Biobank, and Finngen.

**Objective** To elucidate the genetic basis of TN.

**Methods** Using GWAS summary statistics generated from BioVU, the UK Biobank, and Finngen, we performed fixed-effect meta-analysis across 490,912 individuals (1,188 TN cases and 489,724 controls) to identify genetic risk factors for TN. Genome-wide significance was defined as p < 5.0×10−8.

**Results** We identify an intergenic locus on chromosome 1p22.2 flanked by *ZNF326* and *SNORD3G* containing 5 SNPs (rs77449572, rs543311093, rs35117749, rs71666259, and rs116010656) reaching genome-wide significance (p < 5.0 x 10−8), where rs77449572 is the sentinel variant (p = 1.72 x 10−9). The SNP rs77449572 overlaps an enhancer element in cortex-derived neurospheres. In addition, rs71666259 and rs116010656 are located in enhancer elements in embryonic stem cells (HUES48), suggesting potential functional consequences of this locus. We also identify a second locus on chromosome 5q35.1 containing sentinel variant rs62376947 reaching genome-wide significance (p = 2.49 x 10−8).

**Conclusions** To our knowledge, we perform the first GWAS of TN. Future studies should be aimed at understanding the extent to which genetic variation stratifies response to neuropathic pain medication and whether genetic information may be used to identify patients who are likely to benefit (or not) from surgical intervention.

Keywords
*   GWAS
*   Genetics
*   Trigeminal Neuralgia
*   Neuropathic pain

## Introduction

Trigeminal neuralgia (TN) is a debilitating neuropathic pain disease characterized by unilateral paroxysmal, lancinating pain along the distribution of 1 or more branches of the trigeminal nerve. Pain can occur stochastically or be triggered by noxious stimuli such as changes in temperature, light touch, or other sensory stimuli.1 TN affects approximately 27 out of 100,000 individuals, although the true incidence may be significantly higher depending on the population.2 First-line treatment of TN is the sodium channel blocker carbamazepine, based on the hypothesis that neuronal hyperexcitability underlies the pathophysiology of TN.3, 4 However, ∼50% of patients do not experience significant symptom reduction with carbamazepine.5 Agents that modulate GABA signaling (i.e., gabapentin) have been used, especially in individuals with multiple sclerosis.6 When medication fails to result in clinical improvement, surgical treatments such as microvascular decompression (MVD)7, 8 or radiosurgery, among others, are performed.9 However, population-based anatomical studies have identified many asymptomatic individuals with vascular compression of the trigeminal nerve10 and ∼30% of individuals undergoing MVD experience recurrent symptoms,7 suggesting that the etiology of at least some portion of individuals with TN is not due to aberrant anatomy. Furthermore, the clinical course of TN patients without radiographic evidence of trigeminal nerve compression may be more severe (i.e., earlier age at clinical presentation).11 Overall, these data indicate that additional factors contribute to TN risk.

We hypothesized that, analogous to other neuropathic pain disorders,12 there is a substantial genetic component underlying TN risk. Genetic studies of TN have identified 1 single nucleotide polymorphism (SNP) in a serotonin transporter gene in a series of 244 TN patients, loss of function mutations in the GABAA receptor channel gamma-1 subunit (GABRG1), and significant variant burden in the α-1 subunit of Cav3.2 (CACNA1H) in a series of 290 TN patients using whole-exome sequencing.13, 14 However, no genome-wide association studies (GWAS) of TN have been performed, significantly impairing our understanding of the genetic basis of the disease. Thus, to expand on these studies, we perform the first GWAS and largest genetic study of TN to date across 1,188 TN cases from three independent biobanks: BioVU,15 FinnGen,16 and the UK Biobank.17 Here we provide evidence for genetic variation conferring significant risk to TN and suggest avenues for incorporating genetic information into clinical care.

## Methods

### Patient cohort description

We used BioVU,18 a deidentified health record linked to genetic data, to identify 561 patients with TN and 7,359 control individuals. TN cases in BioVU were first identified using ICD9 (350.1) and ICD10 (G50.0) billing codes before being manually reviewed and selected based on a neurologist or neurosurgeon’s diagnosis in the electronic health record. GWAS summary statistics from FinnGen (G6_TRINEU, 438 cases and 121,360 controls)16 and the UK Biobank (189 cases and 361,005 controls)17 were also included. Thus, the total cohort consisted of 490,912 individuals (1,188 TN cases and 489,724 controls). Institutional review board approval was provided by Vanderbilt University.

### Genome wide association study analysis

Whole-genome data in BioVU was obtained using the MEGA array and imputed using the 1000 genome reference panel (phase 3). GWAS quality control and imputation were performed as previously described.18, 19 Age, sex, and five genotype-based principal components (genomic ancestry) were used as covariates in logistic regression models to identify TN-associated SNPs. We used GWAS summary statistics in the UK Biobank17, which had been genotyped using Illumina HumanCoreExome array and imputed using the Haplotype Reference Consortium panel (< 50 cases and 15,582 controls; [http://pheweb.sph.umich.edu/pheno/343](http://pheweb.sph.umich.edu/pheno/343)). SNP association statistics had been calculated using saddlepoint approximation, adjusted for age, sex, and four principal components. Additional information can be found at [http://pheweb.sph.umich.edu/pheno/343](http://pheweb.sph.umich.edu/pheno/343). We also utilized GWAS summary statistics in FinnGen16, which had been genotyped using Illumina and Affymetrix chip arrays and imputed using the SISu v3 imputation reference panel of 3,775 whole genomes. SNP associations had been determined after controlling for age, sex, 10 principal components, and genotype batch as covariates. For additional information on the construction and data elements within the Finngen consortium, see [https://finngen.gitbook.io/documentation/](https://finngen.gitbook.io/documentation/).

We performed fixed-effect meta-analysis and appropriate genomic control using METAL.20 Genome-wide significance was defined as p < 5.0×10−8. Linkage disequilibrium (LD, r2) was determined using the LD matrix tool (NCI) leveraging the 1000 genome European ancestry reference dataset.21

## Results

We performed a genome-wide association study (GWAS) of TN utilizing whole-genome genetic information from 490,912 individuals (1,188 TN cases and 489,724 controls) across three independent DNA biobanks linked to phenotypic information: BioVU, FinnGen, and the UK Biobank. We identified an intergenic locus on chromosome 1p22.2 flanked by two genes – zinc finger protein 326 (*ZNF326)* and small nucleolar RNA C/D box 36 (*SNORD3G*). This locus contains 5 single nucleotide polymorphisms (SNPs; rs77449572, rs543311093, rs35117749, rs71666259, and rs116010656) reaching genome-wide significance (p < 5.0 x 10−8), where rs77449572 is the sentinel variant in the locus (p = 1.72 x 10−9) (Figure 1A). As highlighted by the heatmap of r2 statistic22 values, a number of SNPs in the 1p22.2 locus reaching genome-wide significance with TN are in moderate linkage disequilibrium (LD) (Figure 1B).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/10/2021.02.08.21251349/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/02/10/2021.02.08.21251349/F1)

Figure 1. 
Genome-wide association study of trigeminal neuralgia (TN). (A) Manhattan plot of single nucleotide polymorphisms (SNP) associated with TN (1,188 TG cases and 489,724 controls). (B) Heatmap of r2 values for linkage disequilibrium between SNPs within the locus on chromosome 1 reaching genome-wide significance. (C) Q-Q plot of the distribution of SNP associations with TN. (D) LocusZoom plot of region on chromosome 1p22.2 around the sentinel variant rs77449572 associated with TN. (E) LocusZoom plot of region on chromosome 5q35.1 around the sentinel variant rs62376947 associated with TN.

To infer functional consequences of these variants in biological samples relevant to TN pathophysiology, we utilized Haploreg,23 a resource for exploring functional annotations across non-coding regions of the genome in diverse array of cell and tissue types. This approach identified rs77449572 overlapping an enhancer element in cortex derived neurospheres and germinal matrix. In addition, rs71666259 and rs116010656 overlap enhancer elements in embryonic stem cells (HUES48), suggesting that these genetic variants may confer risk to TN by regulating gene expression. We also identified a second locus on chromosome 5q35.1 containing sentinel variant rs62376947in an intronic region of *SLIT3* reaching genome-wide significance (Figure 1A, p = 2.49 x 10−8).

Population stratification has been adequately accounted for in our association analysis, as confirmed by a Q-Q plot which demonstrates that there is no global inflation in the test statistic – i.e., no early separation of the expected from the observed (Figure 1C LocusZoom plots show the magnitude of the genome-wide significant associations for the sentinel variants and additional variants in the 1p22.2 (Figure 1D) and the 5q35.1 (Figure 1E) loci. Collectively, these data highlight genetic factors conferring risk to TN.

## Discussion

We perform the first GWAS, to our knowledge, and largest genetic study of TN to date, revealing two independent loci reaching genome-wide significance with TN. We exploit 3 of the largest electronic health records databases linked to genetic information data available and identify two independent genome-wide significant loci. Functional genomics analysis using the Haploreg resource23 suggests that variants within the 1p22.2 locus augment gene expression in tissues of neuronal origin. While further molecular studies are needed to determine the precise mechanism by which variation in these loci contribute to TN, this study identifies genetic loci for further experimental testing to understand how genetic variation contributes to the pathophysiology of TN.

*The locus on chromosome* 1p22.2 is located in an intergenic region between *ZNF326* and *SNORD3G. ZNF326* codes for a protein that is a core component of the DBIRD complex that coordinates transcript elongation and alternative splicing at the interface of mRNPs and RNA polymerase II.24 In addition, *ZNF326* is highly expressed in the brain and spinal cord (gtexportal.org) and has been shown to play a role in neuronal differentiation, neurodegeneration through regulation of hippocampal volume, and cancer.25-27 On the other hand, *SNORD3G* codes for a small nucleolar RNA of largely unknown function. Furthermore, *SLIT3* does not appear to be highly expressed in neuronal tissues (gtexportal.org), but has been shown to be a critical regulator of angiogenesis, monocyte migration, and axon regeneration during peripheral nerve damage.28-30 Thus, functional validation of these genes may enable identification of a mechanistic basis for TN pathophysiology and infer development of novel therapeutics.

A number of critical questions remain. First, what is the role of pharmacogenes in regulating response to neuropathic pain medications used to treat TN? Clinical management of TN is challenging, and elucidation of the genetic basis of the disease may inform strategies to understand the role of genetic variation in clinical management.31 Given that ∼50% of TN patients do not experience significant symptom reduction with carbamazepine5 and ∼30% of individuals with TN undergoing MVD have recurrent symptoms,7 it is likely that genetic variation may play a large role, consistent with other neuropathic pain disorders.32 Furthermore, it is possible that use of genetic information may assist in identifying appropriate treatment approaches. Specifically, while HLA-B∗1502 is associated with severe adverse effects related to carbamazepine use,33 variation in cytochrome p450 genes (specifically, CYP3A4) also regulates metabolism of the drug.34, 35 There are likely complex epistatic interactions between these genes, coupled with heretofore unknown rare and common genetic variants influencing carbamazepine (and other neuropathic pain medications) efficacy.

Next, what is the role of genetic variation in defining subtypes of TN? For instance, are patients that present earlier in life, have bilateral symptoms, with no obvious vascular compression of the trigeminal nerve, and/or patients that still experience symptoms post-MVD genetically distinct subtypes? Additional phenotypic refinement and subgroup analysis may reveal a role for genetic risk factors within these populations. For instance, do patients with TN who present earlier harbor large effect-size damaging variants? Further resolution of TN subtypes - stratified by clinical presentation, genetic risk factors, or both – may expand the applicability of genetic information in clinical management. In addition, understanding how genetic variation confers on pharmacologically actionable targets and pathways will improve clinical care. Finally, these data may identify candidates for surgical intervention sooner, therefore limiting the already severe morbidity associated with TN.

## Limitations

While we had access to a deidentified electronic medical record for the subjects in BioVU, we did not have similar phenotyping details for the UK Biobank and Finngen individuals which rely on ICD classifications and summary statistics.16, 17 In addition, we were not able to interrogate the contribution of rare (i.e. *de novo* mutations) which have been shown to play a role in TN.13, 14, 36 There is likely to be environmental contributions and phenotypic heterogeneity37 necessitating evolving approaches to understand the genetic basis of the disease.

## Conclusions

We present results from the first GWAS and largest genetic study of TN to date. These data identify genetic risk factors for TN which will inform our understanding of the genetic basis of the disease and potential avenues for treatment. Future studies should be aimed at understanding the extent to which genetic variation predicts (as part of a polygenic risk score) disease liability across the genome, stratifies response to neuropathic pain medication, and identifies individuals who are likely to benefit from surgical intervention.

## Data Availability

All data will be made available by the corresponding author, Dr. Andrew Hale (andrew.hale{at}vanderbilt.edu).

## Footnotes

*   **Conflict of Interest:** Dr. Gamazon receives an honorarium from the journal *Circulation Research* of the American Heart Association, as a member of the Editorial Board.

*   **Disclosure of Funding:** Dr. Hale receives salary support from the National Institutes of Health (NIH, F30HL143826 and R35HG010718). Processing and imputation of genetic data was supported by a grant awarded to Dr. Hale from the Vanderbilt Institute for Clinical and Translational Research (VR53934.1). L.B. is supported by R01LM010685. Dr. Gamazon is supported by the National Human Genome Research Institute of the National Institutes of Health under Award Numbers R35HG010718 and R01HG011138. Dr. Gamazon has also significantly benefitted from a Fellowship at Clare Hall, University of Cambridge (UK) and is grateful to the President and Fellows of the college for a stimulating intellectual home.

*   Received February 8, 2021.
*   Revision received February 8, 2021.
*   Accepted February 10, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Zakrzewska JM, Linskey ME. Trigeminal neuralgia. BMJ : British Medical Journal. 2014;348:g474.
    
    

2.  2.Hall GC, Carroll D, Parry D, McQuay HJ. Epidemiology and treatment of neuropathic pain: The UK primary care perspective. PAIN. 2006/05/01/ 2006;122(1):156–162.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pain.2006.01.030&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16545908&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000237549100024&link_type=ISI) 

3.  3.Burchiel KJ, Baumann TK. Pathophysiology of trigeminal neuralgia: new evidence from a trigeminal ganglion intraoperative microneurographic recording. Case report. J Neurosurg. Nov 2004;101(5):872–873.
    
    

4.  4.Devor M, Amir R, Rappaport ZH. Pathophysiology of Trigeminal Neuralgia: The Ignition Hypothesis. The Clinical Journal of Pain. 2002;18(1).
    
    

5.  5.Sindrup SH, Jensen TS. Pharmacotherapy of Trigeminal Neuralgia. The Clinical Journal of Pain. 2002;18(1).
    
    

6.  6.Khan OA. Gabapentin relieves trigeminal neuralgia in multiple sclerosis patients. Neurology. 1998;51(2):611–614.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.51.2.611&link_type=DOI) 

7.  7.Barker FG, Jannetta PJ, Bissonette DJ, Larkins MV, Jho HD. The Long-Term Outcome of Microvascular Decompression for Trigeminal Neuralgia. New England Journal of Medicine. 1996;334(17):1077–1084.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM199604253341701&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8598865&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996UG45700001&link_type=ISI) 

8.  8.Holste K, Chan AY, Rolston JD, Englot DJ. Pain Outcomes Following Microvascular Decompression for Drug-Resistant Trigeminal Neuralgia: A Systematic Review and Meta-Analysis. Neurosurgery. 2019;86(2):182–190.
    
    

9.  9.Bick SKB, Eskandar EN. Surgical Treatment of Trigeminal Neuralgia. Neurosurgery clinics of North America. Jul 2017;28(3):429–438.
    
    

10. 10.Jani RH, Hughes MA, Gold MS, Branstetter BF, Ligus ZE, Sekula RF, Jr.. Trigeminal Nerve Compression Without Trigeminal Neuralgia: Intraoperative vs Imaging Evidence. Neurosurgery. Jan 1 2019;84(1):60–65.
    
    

11. 11.Ko AL, Lee A, Raslan AM, Ozpinar A, McCartney S, Burchiel KJ. Trigeminal neuralgia without neurovascular compression presents earlier than trigeminal neuralgia with neurovascular compression. J Neurosurg. Dec 2015;123(6):1519–1527.
    
    

12. 12.Mogil JS. Pain genetics: past, present and future. Trends Genet. Jun 2012;28(6):258–266.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tig.2012.02.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22464640&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305094000002&link_type=ISI) 

13. 13.Cui W, Yu X, Zhang H. The serotonin transporter gene polymorphism is associated with the susceptibility and the pain severity in idiopathic trigeminal neuralgia patients. J Headache Pain. Jun 20 2014;15(1):42.
    
    

14. 14.Dong W, Jin SC, Allocco A, et al. Exome Sequencing Implicates Impaired GABA Signaling and Neuronal Ion Transport in Trigeminal Neuralgia. iScience. 2020/10/23/ 2020;23(10):101552.
    
    

15. 15.Roden DM, Pulley JM, Basford MA, et al. Development of a large-scale de-identified DNA biobank to enable personalized medicine. Clinical pharmacology and therapeutics. Sep 2008;84(3):362–369.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/clpt.2008.89&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18500243&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258582700015&link_type=ISI) 

16. 16.Locke AE, Steinberg KM, Chiang CWK, et al. Exome sequencing of Finnish isolates enhances rare-variant association power. Nature. 2019;572(7769):323–328.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-019-1457-z&link_type=DOI) 

17. 17.Bycroft C, Freeman C, Petkova D, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. Oct 2018;562(7726):203–209.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0579-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30305743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 

18. 18.Unlu G, Qi X, Gamazon ER, et al. Phenome-based approach identifies RIC1-linked Mendelian syndrome through zebrafish models, biobank associations and clinical studies. Nature medicine. 2020/01/01 2020;26(1):98–109.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-019-0705-y&link_type=DOI) 

19. 19.Hale AT, Zhou D, Bastarache L, et al. The genetic architecture of human infectious diseases and pathogen-induced cellular phenotypes. medRxiv. 2020:2020.2007.2019.20157404.
    
    

20. 20.Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. Sep 1 2010;26(17):2190–2191.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq340&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20616382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281738900017&link_type=ISI) 

21. 21.Alexander TA, Machiela MJ. LDpop: an interactive online tool to calculate and visualize geographic LD patterns. BMC bioinformatics. Jan 10 2020;21(1):14.
    
    

22. 22.VanLiere JM, Rosenberg NA. Mathematical properties of the r2 measure of linkage disequilibrium. Theor Popul Biol. Aug 2008;74(1):130–137.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tpb.2008.05.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18572214&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257912400014&link_type=ISI) 

23. 23.Ward LD, Kellis M. HaploReg v4: systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease. Nucleic Acids Res. Jan 4 2016;44(D1):D877–881.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkv1340&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26657631&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 

24. 24.Close P, East P, Dirac-Svejstrup AB, et al. DBIRD complex integrates alternative mRNA splicing with RNA polymerase II transcript elongation. Nature. Mar 25 2012;484(7394):386–389.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature10925&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22446626&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302946500035&link_type=ISI) 

25. 25.Lee JY, Nakane Y, Koshikawa N, Nakayama K, Hayashi M, Takenaga K. Characterization of a zinc finger protein ZAN75: nuclear localization signal, transcriptional activator activity, and expression during neuronal differentiation of P19 cells. DNA Cell Biol. Apr 2000;19(4):227–234.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/104454900314492&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10798446&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000086607300004&link_type=ISI) 

26. 26.Stein JL, Hua X, Morra JH, et al. Genome-wide analysis reveals novel genes influencing temporal lobe structure with relevance to neurodegeneration in Alzheimer’s disease. NeuroImage. Jun 2010;51(2):542–554.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2010.02.068&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20197096&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000277141200005&link_type=ISI) 

27. 27.Yu X, Wang M, Wu J, Han Q, Zhang X. ZNF326 promotes malignant phenotype of glioma by up-regulating HDAC7 expression and activating Wnt pathway. J Exp Clin Cancer Res. Jan 28 2019;38(1):40.
    
    

28. 28.Dun XP, Carr L, Woodley PK, et al. Macrophage-Derived Slit3 Controls Cell Migration and Axon Pathfinding in the Peripheral Nerve Bridge. Cell Rep. Feb 5 2019;26(6):1458-1472.e1454.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2018.12.081&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30726731&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 

29. 29.Geutskens SB, Hordijk PL, van Hennik PB. The chemorepellent Slit3 promotes monocyte migration. J Immunol. Dec 15 2010;185(12):7691–7698.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTE6IjE4NS8xMi83NjkxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDIvMTAvMjAyMS4wMi4wOC4yMTI1MTM0OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

30. 30.Zhang B, Dietrich UM, Geng JG, Bicknell R, Esko JD, Wang L. Repulsive axon guidance molecule Slit3 is a novel angiogenic factor. Blood. Nov 5 2009;114(19):4300–4309.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czoxMToiMTE0LzE5LzQzMDAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMi8xMC8yMDIxLjAyLjA4LjIxMjUxMzQ5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

31. 31.Yamamoto PA, Costa ACC, Lauretti GR,  Moraes NVd. Pharmacogenomics in chronic pain therapy: from disease to treatment and challenges for clinical practice. Pharmacogenomics. 2019;20(13):971–982.
    
    

32. 32.Zorina-Lichtenwalter K, Parisien M, Diatchenko L. Genetic studies of human neuropathic pain conditions: a review. Pain. Mar 2018;159(3):583–594.
    
    

33. 33.Ferrell PB, Jr.., McLeod HL. Carbamazepine, HLA-B∗1502 and risk of Stevens-Johnson syndrome and toxic epidermal necrolysis: US FDA recommendations. Pharmacogenomics. Oct 2008;9(10):1543–1546.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2217/14622416.9.10.1543&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18855540&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000261016200022&link_type=ISI) 

34. 34.Pelkonen O, Myllynen P, Taavitsainen P, et al. Carbamazepine: a ‘blind’ assessment of CVP-associated metabolism and interactions in human liver-derived in vitro systems. Xenobiotica. Jun 2001;31(6):321–343.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00498250110055479&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11513246&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000170215000002&link_type=ISI) 

35. 35.Yuki H, Honma T, Hata M, Hoshino T. Prediction of sites of metabolism in a substrate molecule, instanced by carbamazepine oxidation by CYP3A4. Bioorganic & medicinal chemistry. Jan 15 2012;20(2):775–783.
    
    

36. 36.Tanaka BS, Zhao P, Dib-Hajj FB, et al. A gain-of-function mutation in Nav1.6 in a case of trigeminal neuralgia. Mol Med. Sep 2016;22:338–348.
    
    

37. 37.Love S, Coakham HB. Trigeminal neuralgia: Pathology and pathogenesis. Brain. 2001;124(12):2347–2360.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/124.12.2347&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11701590&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F10%2F2021.02.08.21251349.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000172803400002&link_type=ISI)