- 1Henan Mental Hospital, The Second Affiliated Hospital of Xinxiang Medical University, Xinxiang, China
- 2Henan Key Lab of Biological Psychiatry, International Joint Research Laboratory for Psychiatry and Neuroscience of Henan, Xinxiang Medical University, Xinxiang, China
- 3Key Laboratory of Animal Models and Human Disease Mechanisms of the Chinese Academy of Sciences and Yunnan Province, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China
- 4Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, China
- 5KIZ-CUHK Joint Laboratory of Bioresources and Molecular Research in Common Diseases, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China
Epilepsy (affects about 70 million people worldwide) is one of the most prevalent brain disorders and imposes a huge economic burden on society. Epilepsy has a strong genetic component. In this study, we perform the largest genome-wide meta-analysis of epilepsy (N = 8,00,869 subjects) by integrating four large-scale genome-wide association studies (GWASs) of epilepsy. We identified three genome-wide significant (GWS) (p < 5 × 10–8) risk loci for epilepsy. The risk loci on 7q21.11 [lead single nucleotide polymorphism (SNP) rs11978015, p = 9.26 × 10–9] and 8p23.1 (lead SNP rs28634186, p = 4.39 × 10–8) are newly identified in the present study. Of note, rs11978015 resides in upstream of GRM3, which encodes glutamate metabotropic receptor 3. GRM3 has pivotal roles in neurotransmission and is involved in most aspects of normal brain function. In addition, we also identified three genes (TTC21B, RP11-375N15.2, and TNKS) whose cis-regulated expression level are associated with epilepsy, indicating that risk variants may confer epilepsy risk through regulating the expression of these genes. Our study not only provides new insights into genetic architecture of epilepsy but also prioritizes potential molecular targets (including GRM3 and TTC21B) for development of new drugs and therapeutics for epilepsy.
Introduction
Epilepsy is a common neurological disease characterized with recurrent unprovoked seizures. As one of the most prevalent brain disorders, epilepsy imposes a huge economic burden on society and affects about 70 million people worldwide (Thijs et al., 2019). Accumulating evidence indicate that epilepsy has a strong genetic component (Kjeldsen et al., 2003; Speed et al., 2014; Koeleman, 2018). Twin studies showed that genetic factors (total heritability) account for about 80% of the liability to epilepsy (Famula et al., 1997; Kjeldsen et al., 2003). Recent genome-wide association study (GWAS) estimated that the single nucleotide polymorphism (SNP) heritability (the proportion of variance in liability that can be attributed to SNPs, h2SNP) of genetic generalized epilepsy and focal epilepsy is about 32.1 and 9.2%, respectively (The International League Against Epilepsy Consortium on Complex Epilepsies, 2018), further indicating the pivotal role of common genetic variation in epilepsy. A recent study by Speed et al. (2014) also estimated that common variants collectively explain about 26% of phenotypic variation for all epilepsy. Despite the high heritability of epilepsy, only limited risk variants and loci have been identified by large-scale genetic studies to date (International League Against Epilepsy Consortium on Complex Epilepsies, 2014; The International League Against Epilepsy Consortium on Complex Epilepsies, 2018).
To further identify risk variants and to uncover the missing heritability of epilepsy, in this study, we report the largest genome-wide meta-analysis of epilepsy (N = 8,00,869 subjects) by integrating four large-scale GWASs of epilepsy. The first GWAS (14,534 cases and 24,218 controls) was from a recent study by the International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium) (The International League Against Epilepsy Consortium on Complex Epilepsies, 2018). The second GWAS was from the UK Biobank (Sudlow et al., 2015), genome-wide summary statistics of epilepsy in UK Biobank [5,087 epilepsy cases (Phecode:X345) and 3,95,209 controls] generated by the scalable and accurate implementation of generalized mixed model (SAIGE) (Zhou et al., 2018) were used in this study. The third GWAS was from a recent study by Ishigaki et al. (2020) (2,143 epilepsy cases and 2,10,310 controls). The fourth GWAS dataset was from the FINNGEN1 [4,588 epilepsy cases (phenocode: G6_EPLEPSY) and 1,44,780 controls].
Materials and Methods
Epilepsy GWAS From the ILAE Consortium
International League Against Epilepsy Consortium (The International League Against Epilepsy Consortium on Complex Epilepsies, 2018) conducted a large-scale trans-ethnic meta-analysis (15,212 epilepsy cases and 29,677 controls) by combining genome-wide associations conducted in Caucasians, Asians, and Africans. Epilepsy cases were classified into three broad categories and seven subtypes. In this study, we used the genome-wide associations from all epilepsy cases (including focal epilepsy, genetic generalized epilepsy, and unclassified epilepsy). In addition, considering that the number of cases in Asian (n = 531) and African (n = 147) GWASs were quite small, only associations from Europeans (14,534 cases and 24,218 controls) were included in our study. Detailed information about diagnosis and classification of epilepsy cases, genotyping, imputation, and quality controls have been described in the original study (The International League Against Epilepsy Consortium on Complex Epilepsies, 2018). The linear mixed model implemented in BOLT-LMM (Loh et al., 2015) was used to test the associations between genetic variants and epilepsy. As the effect size reported by BOLT-LMM is a Beta coefficient, we transformed Beta values into odds ratio (OR) using the method developed by Lloyd-Jones et al. (2018).2 The parameters used for transformation are as follows: k (the prevalence of epilepsy in the GWAS sample) = 0.38 (14,534/38,752). Allele 1, allele 2, allele frequency, Beta, SE, and p-values were adopted from the original study (The International League Against Epilepsy Consortium on Complex Epilepsies, 2018). The SE of ln(OR) was calculated using the following formula: SE = ln(OR)/Z (where Z-value was calculated using qnorm function implemented in R software (v3.4.1).
Epilepsy GWAS From UK Biobank
UK Biobank is an open access resource that aims to identify the associations between common genetic variation and multiple complex diseases and traits (Sudlow et al., 2015; Bycroft et al., 2018). Detailed information about UK Biobank have been described previously (Bycroft et al., 2018). In this study, we used the summary statistics of epilepsy GWAS from the PheWeb (Gagliano Taliun et al., 2020).3 Briefly, SAIGE (Zhou et al., 2018) was used to test the associations between genetic variants and phenotypes included in UK Biobank. The phenotype code for epilepsy in PheWeb is X345 (corresponding to ICD9 345, or ICD10 code G40). A total of 5,087 epilepsy and 3,95,209 controls were included in GWAS. Unlike the ILAE Consortium, epilepsy cases in UK Biobank were not divided into focal and generalized epilepsy. We transformed the Beta effect size into OR using the following formula: OR = exp(Beta). In addition, SNPs with minor allele frequency (MAF) of less than 0.01 were excluded. More detailed information about the UK Biobank and PheWeb can be found in previous papers (Bycroft et al., 2018; Zhou et al., 2018; Gagliano Taliun et al., 2020).
Epilepsy GWAS of Japanese Population
The genome-wide associations of Japanese population were from the study of Ishigaki et al. (2020). Briefly, Ishigaki et al. (2020) conducted a GWAS of 42 diseases in a large-scale Japanese population. Associations between genetic variants and diseases were assessed using SAIGE (Zhou et al., 2018). The GWAS of epilepsy included in the study of Ishigaki et al. (2020) contained 2,143 cases and 2,10,310 controls. Summary statistics were downloaded from the Japanese ENcyclopedia of GEnetic associations by Riken (JENGER).4 Detailed information about the subjects, diagnosis, genotyping, imputation, quality control, and statistical analysis have been described in the study of Ishigaki et al. (2020). We transformed the Beta effect size into OR using the following formula: OR = exp(Beta). Besides, SNPs with MAF of less than 0.01 were excluded for meta-analysis.
Epilepsy GWAS From FINNGEN
FINNGEN is a large-scale and open access resource which aims to improve human health through genetic research. Its ultimate goal is to identify new therapeutic targets and diagnostics for treating human diseases. It was launched in 2017, and it will collect genome information and digital healthcare data of about 5,00,000 Finnish people. By integrating genetic resources and digital healthcare data from multiple organizations, including Finnish universities, biobanks, hospitals, and so on, this project expects to achieve breakthroughs in disease diagnosis, prevention, and treatment. GWASs in FINNGEN were also performed using SAIGE (Zhou et al., 2018) and summary statistics from data freeze four results [4,588 epilepsy cases (phenocode: G6_EPLEPSY) and 1,44,780 controls] of FINNGEN were used in this study. Detailed information about the subjects, diagnosis, genotyping, imputation, quality control, and statistical analysis have been described in the homepage of FINNGEN: https://www.finngen.fi/en. We transformed the Beta effect size into OR using the following formula: OR = exp(Beta). Besides, SNPs with MAF of less than 0.01 were excluded for meta-analysis. Of note, the ILAE GWAS performed analysis on three subphenotypes of epilepsy (focal, GGE, and all epilepsy) (International League Against Epilepsy Consortium on Complex Epilepsies, 2014). In addition, cases in FinnGen included focal epilepsy, generalized epilepsy, and epilepsy (broader sense). However, the cases in Japanese and UK Biobank GWASs were not divided into subphenotypes (Ishigaki et al., 2020). Therefore, we did not divide the epilepsy into subphenotypes (i.e., we used all epilepsy cases) in this study.
Meta-Analysis
For the other three datasets (including UK Biobank, Japanese GWAS, and FINNGEN; effect size was calculated using SAIGE, which uses logistic repression to perform association test), the effect size of each SNP was converted into OR using the following formula: OR = expBeta. SNPs with MAF of less than 0.01 were excluded. Meta-analysis was performed using PLINK (V1.90) (Purcell et al., 2007), with the use of fixed-effect model. The inverse variance-based analysis to be implemented in PLINK was used for meta-analysis.
Tissue- and Cell-Type Enrichment Analysis
We used MAGMA (de Leeuw et al., 2015) [implemented in FUMA (Watanabe et al., 2017)] to perform tissue- and cell-type enrichment analysis. MAGMA is a powerful tool for gene and gene-set analysis, and it uses GWAS summary statistics as input. MAGMA first derives a gene-level p-value by using a multiple linear principal component regression model. The gene-level p-values were then used for further gene-set analysis. To test if the genetic associations from GWASs are enriched in specific tissues, MAGMA utilizes gene expression data from the GTEx (53 human tissues) for tissue enrichment analysis.
We also conducted single-cell enrichment analysis by MAGMA followed by the methods of a recent published paper about Parkinson’s disease (Bryois et al., 2020). The single-cell RNA-seq data was from mouse central nervous system which include 1,60,769 single cells in total (Zeisel et al., 2018); top 10% genes that ranked by gene expression specificity of each cell type was remained for the MAGMA gene-set enrichment analysis (de Leeuw et al., 2015). For further detailed information about single-cell data generation, data processing, and gene expression specificity calculation, please refer to the original paper (Zeisel et al., 2018; Bryois et al., 2020). Detailed information about MAGMA, FUMA, and tissue- and cell-type enrichment analyses can be found in the original papers (de Leeuw et al., 2015; Watanabe et al., 2017; Zeisel et al., 2018; Bryois et al., 2020) and the FUMA website.5
Transcriptome-Wide Association Study
Transcriptome-wide association study (TWAS) aims to identify genes whose genetically regulated expression level are associated with complex human diseases or traits (Gusev et al., 2016). TWAS firstly uses an external expression reference panel (which contained gene expression data and genome-wide SNPs) to establish SNP-expression weights (i.e., SNP-gene expression correlations). These SNP-expression weights are then used to predict the expression level of genes in individuals included in GWAS. Finally, statistical inferences are made to test if the expression level of a gene is associated with diseases or traits. Considering that epilepsy is a brain disorder, in this study, we utilized the gene expression and genotype data from the PsychENCODE (Gandal et al., 2018) as the reference panel to construct the SNP-gene expression weights. In brief, PsychENCODE integrated gene expression (from human brain tissues, most of tissues are the prefrontal cortex) and genotype data of over 2,000 human subjects. Gene expression was quantified with RNA sequencing, and genotypes were determined using SNP arrays. Among the subjects included in PsychENCODE, expression data, and genotypes of 1,321 indivuduals (only adult individuals with matching gene expression and genotypes can be used for expression quantitative trait loci (eQTL) analysis) were used for eQTL analysis and construction of SNP-gene expression weights. We used FUSION (Gusev et al., 2016) pipeline to prepare the SNP-gene expression weights. We integrated the constructed SNP-gene expression weights and GWAS summary statistic from the meta-analysis to conduct TWAS, with the use of default parameters and settings. Results of TWAS were corrected by the Bonferroni correction approach. More detailed information about TWAS and FUSION can be found in the original paper (Lloyd-Jones et al., 2018) and FUSION website.6
Interaction Analysis Between the Risk Genes and Drugs
To explore if the identified risk genes may be targeted as potential therapeutic targets, we explored the interaction between the identified epilepsy risk genes and drugs using the drug-gene interactions database (DGIdb7) (Freshour et al., 2020). Briefly, DGIdb collected interactions between genes and drugs from several well-characterized databases, including DrugBank, Drug Target Commons, PharmGKB, Chembl, etc. In addition, pathways and molecular functions of genes were also considered. These combined information were then used to predict the interactions between genes and drugs.
Results
We firstly examined the genome-wide significant (GWS) loci in each GWAS. For epilepsy GWAS from the ILAE Consortium, two loci reached GWS level (of note, this dataset contains three subtypes and we only included GWAS of all epilepsy) (The International League Against Epilepsy Consortium on Complex Epilepsies, 2018). For epilepsy GWAS from UK Biobank, Japanese population, and FINNGEN, no loci showed GWS associations with epilepsy (Supplementary Figures 1–3). We then conducted a genome-wide meta-analysis by combining the genome-wide association results of the four studies. The QQ plot is shown in Supplementary Figure 4. The lambdaGC (λ1,000, genomic control (GC) inflation lambda scaled for 1,000 cases and 1,000 controls) of our meta-analysis is 1.0027, indicating that the association signals were mainly driven by polynenicity rather than population structure. In this largest trans-ethnic meta-analysis (26,352 cases and 7,74,517 controls), we identified three GWS (p < 5 × 10–8) risk loci for epilepsy (Figure 1 and Table 1). The risk locus on 2q24.3 has been reported in a previous GWAS (The International League Against Epilepsy Consortium on Complex Epilepsies, 2018). Of note, the most significant variant (lead SNP rs11890028) on 2q24.3 resides in the intron 7 of SCN1A, a well-characterized risk gene for epilepsy (Claes et al., 2001; Parihar and Ganesh, 2013; International League Against Epilepsy Consortium on Complex Epilepsies, 2014; Figure 2B). However, the risk loci on 7q21.11 (lead SNP rs11978015, p = 9.26 × 10–9) and 8p23.1 (lead SNP rs28634186, p = 4.39 × 10–8) are newly identified in the present study (Figures 2B,C). Of note, rs11978015 resides in upstream of GRM3 (Figure 2B), which encodes glutamate metabotropic receptor 3. GRM3 has pivotal roles in neurotransmission and is involved in most aspects of normal brain function (Blacker et al., 2017; Chaki, 2017; Jin et al., 2018; Neale and Olszewski, 2019; Joffe et al., 2021; Kellner et al., 2021). Another GWS risk variant rs28634186 is located in intergenic region and the nearest gene for rs28634186 is TNKS (Figure 2C). Genetic correlation between epilepsy and other diseases [using LD Hub (Zheng et al., 2017)]8 showed significant correlations with amyotrophic lateral sclerosis (ALS), schizophrenia, and bipolar disorder (Supplementary Table 1).
Figure 1. Meta-analysis of four large-scale genome-wide association studies (GWASs) identified two novel risk loci for epilepsy. Manhattan plot of the meta-analysis of epilepsy (26,352 cases and 7,74,517 controls). The risk loci on 7q21.11 (lead SNP rs11978015, p = 9.26 × 10– 9) and 8p23.1 (lead SNP rs28634186, p = 4.39 × 10– 8) are newly identified in the present study.
Figure 2. Locuszoom plots of the three genome-wide significant (GWS) loci. (A) The lead SNP rs11890028 (P = 7.76 × 10– 13) of the 2q24.3 risk locus resides in intronic region of SCN1A. (B) The lead SNP rs11978015 (P = 9.26 × 10– 9) of the 7q21.11 risk locus resides upstream of GRM3. (C) The nearest gene for the lead SNP rs28634186 (P = 4.39 × 10– 8) on the 8p23.1 is TNKS. The read dotted line represents the genome-wide significance level (5.0 × 10– 8).
To identify the tissues and cell types that risk genes may exert their biological effects on epilepsy, we further performed tissue- and cell-type-specific enrichment analysis using MAGMA (see text footnote 5). As expected, the GWAS associations were significantly enriched in brain tissues (Figure 3A), with the highest enrichment in the cerebellar hemisphere and frontal cortex (Figure 3A). Cell-type-specific enrichment analysis showed significant enrichment of GWAS associations in telencephalon projecting inhibitory and excitatory neurons (Figure 3B).
Figure 3. Tissue- and cell-type enrichments of epilepsy GWAS associations. (A) Tissues that showed significant enrichment (corrected p < 0.05) are shown in red. (B) Cell types that showed significant enrichment (corrected p < 0.05) are shown in red.
To identify risk genes whose genetically regulated expression change are associated with epilepsy, we further conducted a TWAS by integrating genome-wide summary statistics of epilepsy (from meta-analysis) and SNP-gene expression weights from the PsychENCODE (N = 1,371).9 Three genes (TTC21B, RP11-375N15.2, and TNKS) showed transcriptome-wide significant (TWS) associations (Bonferroni corrected p < 0.05) with epilepsy (Figure 4), indicating that risk variants may confer epilepsy risk through regulating the expression of these genes.
Figure 4. Transcriptome-wide association study (TWAS) results of epilepsy. Three genes (TTC21B, RP11-375N15.2, and TNKS) showed transcriptome-wide significant (TWS) associations with epilepsy, indicating that the genetically regulated expression of these genes are associated with epilepsy. Transcriptome-wide significance (corrected by Bonferroni adjustment) is marked by a red dotted line.
Finally, to explore if the risk identified may be targeted for epilepsy treatment, we further examined the interactions between the identified risk genes and drugs using DGIdb (Freshour et al., 2020). Our analysis showed that SCN1A interacts with many drugs (Supplementary Table 2), suggesting this gene may be targeted as a potential therapeutic target. In addition, GRM3 and TNKS also show interactions with several drugs (Supplementary Tables 3, 4). These data suggest that these three genes may be targeted for epilepsy treatment.
Discussion
In summary, we performed the largest meta-analysis of epilepsy GWAS in this study and identified two new risk loci for epilepsy. We also showed that genetic associations of epilepsy are enriched in brain tissues and telencephalon projecting inhibitory and excitatory neurons. Of note, we identified three risk genes (TTC21B, RP11-375N15.2, and TNKS) whose expression perturbation may have a role in epilepsy. Interestingly, TTC21B resides in the 2q24.3 locus, a region that contains the well-characterized epilepsy risk gene SCN1A (Claes et al., 2001, 2009; Meng et al., 2013; Parihar and Ganesh, 2013; International League Against Epilepsy Consortium on Complex Epilepsies, 2014; Haigh et al., 2021). Although the genome-wide association signal (Figure 2A) and previous studies (Claes et al., 2001, 2009; Meng et al., 2013; Parihar and Ganesh, 2013; International League Against Epilepsy Consortium on Complex Epilepsies, 2014; Haigh et al., 2021) have clearly showed that SCN1A represents the most possible causal gene for this risk locus, our TWAS suggested that TTC21B may also have a potential role in epilepsy (Figure 4). In fact, previous studies also have revealed the potential role of TTC21B in epilepsy (Mirza et al., 2017; The International League Against Epilepsy Consortium on Complex Epilepsies, 2018). More work is needed to elucidate the role of TTC21B in epilepsy.
We noticed a second peak which is not in linkage disequilibrium (LD) with the main signal at 2q24.3 (located at the right of the top hit) (Figure 2A), suggesting two independent genetic signals at this locus. We thus performed a conditional analysis using genome-wide complex trait analysis (GCTA) (Yang et al., 2011). Conditional analysis on the top hit SNP rs11890028 suggested that the second peak might not be an independent signal. The p-value of rs11896706 (located at the right of the top hit, Figure 2A) conditioned on rs11890028 (the top hit) is 5.01 × 10–5. In addition, we also checked the LD pattern of genetic variants spanning this genomic region [using genotype data of Europeans from the 1,000 Genomes project (The 1000 Genomes Project Consortium, 2015)]. We found that the lead SNP (rs11890028) showed weak LD (r2 = 0.13) with SNP rs11896706 in this region, suggesting that the second signal peak in this locus was likely due to LD between the lead SNP rs11890028 and nearby variants. More work is needed to further investigate if two independent GWAS signals at this locus.
Intriguingly, our findings suggested that GRM3 may be a potential risk gene for epilepsy. GRM3 encodes glutamate metabotropic receptor 3 (mGluR3), a member of the family of G protein-coupled receptors. As a major receptor of glutamate, postsynaptic GRM3 is crucial for mGluR3-dependent long-term depression (LTD) (Joffe et al., 2021) and cognitive function (Jin et al., 2018; Neale and Olszewski, 2019). The binding of glutamate to mGluR3 results in activation of G protein-coupled receptor, which in turn regulates gene transcription, release of neurotransmitter, neuron activity, and synaptic transmission (Blacker et al., 2017; Chaki, 2017; Kellner et al., 2021). In fact, a recent study has proposed that modulation of astrocyte glutamate uptake (and/or mGluR activation) may represent a potential therapeutic approach for epilepsy treatment (Peterson and Binder, 2020). These lines of evidence suggest that GRM3 may have a role in epilepsy.
Our GWAS meta-analysis suggested that 8p23.1 is a risk locus for epilepsy. To further explore the potential risk gene in this locus, we searched literatures and performed additional analysis. Based on the results of GWAS and TWAS, we speculated that TNKS (also known as PARP5A) may be the potential risk gene at this locus. First, most of the significant variants identified in our GWAS meta-analysis are located in genomic region (including gene body and upstream of TNKS) containing TNKS (Figure 2C). Second, our TWAS results suggested that TNKS is a gene whose genetically regulated expression may have a role in epilepsy (Figure 4). Third, expression analysis using single-cell RNA-seq data showed that TNKS is highly expressed in different cell types of the brain (Supplementary Figures 5, 6). Fourth, the biological function also supports TNKS as a potential risk gene. TNKS encodes tankyrases [members of the poly(ADP-ribose) polymerase (PARP) family] that are involved in the regulation of Wnt/beta-catenin signaling (Bao et al., 2012; Kulak et al., 2015; Yang et al., 2019) and PTEN (Li et al., 2015). Considering the important role of Wnt/beta-catenin signaling (Bengoa-Vergniory and Kypta, 2015; Arnes and Casas Tinto, 2017; Noelanders and Vleminckx, 2017; Bem et al., 2019) and PTEN (Endersby and Baker, 2008; Zhou and Parada, 2012; Spina Nagy et al., 2021) in the brain, it is possible that TNKS confers risk of epilepsy by regulating Wnt/beta-catenin signaling pathway. In fact, recent studies also have showed that TNKS modulates TDP-43, a protein with a central role in ALS and frontotemporal degeneration (FTD) (McGurk et al., 2020; Tanji et al., 2021), further suggesting the pivotal role of TNKS in the human brain. Finally, Wnt/β-catenin signaling was proposed to be a potential target for epilepsy therapy (Huang et al., 2015; Hodges and Lugo, 2018). Collectively, these evidence suggest that TNKS may represent the potential risk gene at 8p23.1. However, further genetic studies and functional characterization are needed to validate if TNKS is a risk gene for epilepsy.
It should be noted that biobank data as a resource for epilepsy GWAS has weakness and limitations. Compared with traditional GWAS (GWASs usually use well-phenotyped cohort), the subjects included in biobank data are usually based on electronic records or questionnaire, which may affect the accurate phenotyping of the included individuals. In addition, the number of cases included in biobank data is usually much smaller than controls (case-control imbalance), which influences the statistic power of biobank data. Finally, population relatedness or structure also needs to be carefully considered and controlled in biobank data.
Conclusion
In conclusion, our study not only provides new insights into genetic architecture of epilepsy but also prioritizes potential molecular targets (including GRM3 and TNKS) for development of new drugs and therapeutics for epilepsy.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s. Custom codes used (including Perl, R, and Unix) were used for data processing. These custom codes can be made available from the corresponding author upon request.
Ethics Statement
The studies involving human participants were reviewed and approved by the ILAE Consortium, UK Biobank, JENGER, and FINNGEN. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
X-JL conceived, designed, and supervised the whole study and wrote and revised the manuscript. MS, X-JL, YY, LL, and WL collected the GWAS summary statistics, performed the data processing and transformation, and conducted the GWAS meta-analysis. JL performed the tissue and cell type-specific enrichments analyses, TWAS, and generated the Manhattan and locus zoom plots. LL provided critical comments for the manuscript improvement. All authors read this manuscript carefully, provided critical comments, and approved the manuscript.
Funding
This study was equally supported by the Innovative Research Team of Science and Technology Department of Yunnan Province (2019HC004) and the Distinguished Young Scientists grant of the Yunnan Province (202001AV070006) to X-JL. This study was also supported by the National Nature Science Foundation of China (31970561 to X-JL), the Western Light Innovative Research Team of Chinses Academy of Sciences, and the major science and technology projects of Henan Province (201300310200 to WL and LL).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We acknowledge the participants and investigators of FinnGen study (https://www.finngen.fi/en), UK Biobank study (https://www.ukbiobank.ac.uk/, ftp://share.sph.umich.edu/UKBB_SAIGE_HRC/), JENGER study (http://jenger.riken.jp/en/), and the International League Against Epilepsy Consortium on Complex Epilepsies. We also appreciate Wei Zhou (Massachusetts General Hospital, Harvard Medical School), Seunggeun Lee (Department of Computational Medicine and Bioinformatics, University of Michigan), and Luke R. Lloyd-Jones (Institute for Molecular Bioscience, University of Queensland) for their help in OR transformation.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2021.722592/full#supplementary-material
Footnotes
- ^ https://www.finngen.fi/en
- ^ https://github.com/lukelloydjones/ORShiny
- ^ http://pheweb.sph.umich.edu/SAIGE-UKB/
- ^ http://jenger.riken.jp/en/
- ^ https://fuma.ctglab.nl/
- ^ http://gusevlab.org/projects/fusion/
- ^ https://www.dgidb.org/
- ^ http://www.nealelab.is/tools-and-software
- ^ http://www.psychencode.org/
References
Arnes, M., and Casas Tinto, S. (2017). Aberrant Wnt signaling: a special focus in CNS diseases. J. Neurogenet. 31, 216–222. doi: 10.1080/01677063.2017.1338696
Bao, R., Christova, T., Song, S., Angers, S., Yan, X., and Attisano, L. (2012). Inhibition of tankyrases induces Axin stabilization and blocks Wnt signalling in breast cancer cells. PLoS One 7:e48670. doi: 10.1371/journal.pone.0048670
Bem, J., Brożko, N., Chakraborty, C., Lipiec, M. A., Koziński, K., Nagalski, A., et al. (2019). Wnt/beta-catenin signaling in brain development and mental disorders: keeping TCF7L2 in mind. FEBS Lett. 593, 1654–1674. doi: 10.1002/1873-3468.13502
Bengoa-Vergniory, N., and Kypta, R. M. (2015). Canonical and noncanonical Wnt signaling in neural stem/progenitor cells. Cell. Mol. Life Sci. 72, 4157–4172. doi: 10.1007/s00018-015-2028-6
Blacker, C. J., Lewis, C. P., Frye, M. A., and Veldic, M. (2017). Metabotropic glutamate receptors as emerging research targets in bipolar disorder. Psychiatry Res. 257, 327–337. doi: 10.1016/j.psychres.2017.07.059
Bryois, J., Skene, N. G., Hansen, T. F., Kogelman, L. J. A., Watson, H. J., Liu, Z., et al. (2020). Genetic identification of cell types underlying brain complex traits yields insights into the etiology of Parkinson’s disease. Nat. Genet. 52, 482–493. doi: 10.1038/s41588-020-0610-9
Bycroft, C., Freeman, C., Petkova, D., Band, G., Elliott, L. T., Sharp, K., et al. (2018). The UK Biobank resource with deep phenotyping and genomic data. Nature 562, 203–209. doi: 10.1038/s41586-018-0579-z
Chaki, S. (2017). mGlu2/3 receptor antagonists as novel antidepressants. Trends Pharmacol. Sci. 38, 569–580. doi: 10.1016/j.tips.2017.03.008
Claes, L., Del-Favero, J., Ceulemans, B., Lagae, L., Van Broeckhoven, C., and De Jonghe, P. (2001). De novo mutations in the sodium-channel gene SCN1A cause severe myoclonic epilepsy of infancy. Am. J. Hum. Genet. 68, 1327–1332. doi: 10.1086/320609
Claes, L. R., Deprez, L., Suls, A., Baets, J., Smets, K., Van Dyck, T., et al. (2009). The SCN1A variant database: a novel research and diagnostic tool. Hum. Mutat. 30, E904–E920.
de Leeuw, C. A., Mooij, J. M., Heskes, T., and Posthuma, D. (2015). MAGMA: generalized gene-set analysis of GWAS data. PLoS Comput. Biol. 11:e1004219. doi: 10.1371/journal.pcbi.1004219
Endersby, R., and Baker, S. J. (2008). PTEN signaling in brain: neuropathology and tumorigenesis. Oncogene 27, 5416–5430. doi: 10.1038/onc.2008.239
Famula, T. R., Oberbauer, A. M., and Brown, K. N. (1997). Heritability of epileptic seizures in the Belgian tervueren. J. Small Anim. Pract. 38, 349–352. doi: 10.1111/j.1748-5827.1997.tb03483.x
Freshour, S. L., Kiwala, S., Cotto, K. C., Coffman, A. C., McMichael, J. F., Song, J. J., et al. (2020). Integration of the Drug-Gene Interaction Database (DGIdb 4.0) with open crowdsource efforts. Nucleic Acids Res. 49, D1144–D1151.
Gagliano Taliun, S. A., VandeHaar, P., Boughton, A. P., Welch, R. P., Taliun, D., Schmidt, E. M., et al. (2020). Exploring and visualizing large-scale genetic associations by using PheWeb. Nat. Genet. 52, 550–552. doi: 10.1038/s41588-020-0622-5
Gandal, M. J., Zhang, P., Hadjimichael, E., Walker, R. L., Chen, C., Liu, S., et al. (2018). Transcriptome-wide isoform-level dysregulation in ASD, schizophrenia, and bipolar disorder. Science 362:eaat8127. doi: 10.1126/science.aat8127
Gusev, A., Ko, A., Shi, H., Bhatia, G., Chung, W., Penninx, B. W., et al. (2016). Integrative approaches for large-scale transcriptome-wide association studies. Nat. Genet. 48, 245–252.
Haigh, J. L., Adhikari, A., Copping, N. A., Stradleigh, T., Wade, A. A., Catta Preta, R., et al. (2021). Deletion of a non-canonical regulatory sequence causes loss of Scn1a expression and epileptic phenotypes in mice. Genome Med. 13:69. doi: 10.1186/s13073-021-00884-0
Hodges, S. L., and Lugo, J. N. (2018). Wnt/beta-catenin signaling as a potential target for novel epilepsy therapies. Epilepsy Res. 146, 9–16. doi: 10.1016/j.eplepsyres.2018.07.002
Huang, C., Fu, X. H., Zhou, D., and Li, J. M. (2015). The Role of Wnt/beta-catenin signaling pathway in disrupted hippocampal neurogenesis of temporal lobe epilepsy: a potential therapeutic target? Neurochem. Res. 40, 1319–1332. doi: 10.1007/s11064-015-1614-1
International League Against Epilepsy Consortium on Complex Epilepsies (2014). Genetic determinants of common epilepsies: a meta-analysis of genome-wide association studies. Lancet Neurol. 13, 893–903. doi: 10.1016/s1474-4422(14)70171-1
Ishigaki, K., Akiyama, M., Kanai, M., Takahashi, A., Kawakami, E., Sugishita, H., et al. (2020). Large-scale genome-wide association study in a Japanese population identifies novel susceptibility loci across different diseases. Nat. Genet. 52, 669–679.
Jin, L. E., Wang, M., Galvin, V. C., Lightbourne, T. C., Conn, P. J., Arnsten, A. F. T., et al. (2018). mGluR2 versus mGluR3 metabotropic glutamate receptors in primate dorsolateral prefrontal cortex: postsynaptic mGluR3 strengthen working memory networks. Cereb. Cortex 28, 974–987. doi: 10.1093/cercor/bhx005
Joffe, M. E., Santiago, C. I., Vermudez, S. A. D., Fisher, N. M., Dogra, S., Niswender, C. M., et al. (2021). Frontal cortex genetic ablation of metabotropic glutamate receptor subtype 3 (mGlu3) impairs postsynaptic plasticity and modulates affective behaviors. Neuropsychopharmacology doi: 10.1038/s41386-021-01041-2 [online ahead of print].
Kellner, V., Kersbergen, C. J., Li, S., Babola, T. A., Saher, G., and Bergles, D. E. (2021). Dual metabotropic glutamate receptor signaling enables coordination of astrocyte and neuron activity in developing sensory domains. Neuron doi: 10.1016/j.neuron.2021.06.010 [online ahead of print].
Kjeldsen, M. J., Corey, L. A., Christensen, K., and Friis, M. L. (2003). Epileptic seizures and syndromes in twins: the importance of genetic factors. Epilepsy Res. 55, 137–146. doi: 10.1016/s0920-1211(03)00117-7
Koeleman, B. P. C. (2018). What do genetic studies tell us about the heritable basis of common epilepsy? Polygenic or complex epilepsy? Neurosci. Lett. 667, 10–16. doi: 10.1016/j.neulet.2017.03.042
Kulak, O., Chen, H., Holohan, B., Wu, X., He, H., Borek, D., et al. (2015). Disruption of Wnt/beta-catenin signaling and telomeric shortening are inextricable consequences of tankyrase inhibition in human cells. Mol. Cell. Biol. 35, 2425–2435. doi: 10.1128/mcb.00392-15
Li, N., Zhang, Y., Han, X., Liang, K., Wang, J., Feng, L., et al. (2015). Poly-ADP ribosylation of PTEN by tankyrases promotes PTEN degradation and tumor growth. Genes Dev. 29, 157–170. doi: 10.1101/gad.251785.114
Lloyd-Jones, L. R., Robinson, M. R., Yang, J., and Visscher, P. M. (2018). Transformation of summary statistics from linear mixed model association on all-or-none traits to odds ratio. Genetics 208, 1397–1408. doi: 10.1534/genetics.117.300360
Loh, P. R., Tucker, G., Bulik-Sullivan, B. K., Vilhjálmsson, B. J., Finucane, H. K., Salem, R. M., et al. (2015). Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nat. Genet. 47, 284–290. doi: 10.1038/ng.3190
McGurk, L., Rifai, O. M., and Bonini, N. M. (2020). TDP-43, a protein central to amyotrophic lateral sclerosis, is destabilized by tankyrase-1 and -2. J. Cell Sci. 133:jcs245811. doi: 10.1242/jcs.245811
Meng, H., Xu, H. Q., Yu, L., Lin, G. W., He, N., Su, T., et al. (2013). The SCN1A mutation database: updating information and analysis of the relationships among genotype, functional alteration, and phenotype. Hum. Mutat. 36, 573–580. doi: 10.1002/humu.22782
Mirza, N., Appleton, R., Burn, S., du Plessis, D., Duncan, R., Farah, J. O., et al. (2017). Genetic regulation of gene expression in the epileptic human hippocampus. Hum. Mol. Genet. 26, 1759–1769. doi: 10.1093/hmg/ddx061
Neale, J. H., and Olszewski, R. (2019). A role for N-acetylaspartylglutamate (NAAG) and mGluR3 in cognition. Neurobiol. Learn. Mem. 158, 9–13. doi: 10.1016/j.nlm.2019.01.006
Noelanders, R., and Vleminckx, K. (2017). How Wnt signaling builds the brain: bridging development and disease. Neuroscientist 23, 314–329. doi: 10.1177/1073858416667270
Parihar, R., and Ganesh, S. (2013). The SCN1A gene variants and epileptic encephalopathies. J. Hum. Genet. 58, 573–580. doi: 10.1038/jhg.2013.77
Peterson, A. R., and Binder, D. K. (2020). Astrocyte glutamate uptake and signaling as novel targets for antiepileptogenic therapy. Front. Neurol. 11:1006.
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795
Speed, D., O’Brien, T. J., Palotie, A., Shkura, K., Marson, A. G., Balding, D. J., et al. (2014). Describing the genetic architecture of epilepsy through heritability analysis. Brain 137, 2680–2689. doi: 10.1093/brain/awu206
Spina Nagy, G., Kawamoto, E. M., and Bridi, J. C. (2021). The role of PTEN signaling in synaptic function: implications in autism spectrum disorder. Neurosci. Lett. 759:136015. doi: 10.1016/j.neulet.2021.136015
Sudlow, C., Gallacher, J., Allen, N., Beral, V., Burton, P., Danesh, J., et al. (2015). UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 12:e1001779. doi: 10.1371/journal.pmed.1001779
Tanji, K., Mori, F., Shirai, F., Fukami, T., Seimiya, H., Utsumi, J., et al. (2021). Novel tankyrase inhibitors suppress TDP-43 aggregate formation. Biochem. Biophys. Res. Commun. 537, 85–92. doi: 10.1016/j.bbrc.2020.12.037
The 1000 Genomes Project Consortium (2015). A global reference for human genetic variation. Nature 526, 68–74.
The International League Against Epilepsy Consortium on Complex Epilepsies (2018). Genome-wide mega-analysis identifies 16 loci and highlights diverse biological mechanisms in the common epilepsies. Nat. Commun. 9:5269.
Thijs, R. D., Surges, R., O’Brien, T. J., and Sander, J. W. (2019). Epilepsy in adults. Lancet 393, 689–701.
Watanabe, K., Taskesen, E., van Bochoven, A., and Posthuma, D. (2017). Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 8:1826.
Yang, H. Y., Shen, J. X., Wang, Y., Liu, Y., Shen, D. Y., and Quan, S. (2019). Tankyrase promotes aerobic glycolysis and proliferation of ovarian cancer through activation of Wnt/beta-catenin signaling. Biomed. Res. Int. 2019:2686340.
Yang, J., Lee, S. H., Goddard, M. E., and Visscher, P. M. (2011). GCTA: a tool for genome-wide complex trait analysis. Am. J. Hum. Genet. 88, 76–82. doi: 10.1016/j.ajhg.2010.11.011
Zeisel, A., Hochgerner, H., Lönnerberg, P., Johnsson, A., Memic, F., van der Zwan, J., et al. (2018). Molecular architecture of the mouse nervous system. Cell 174, 999–1014.e1022.
Zheng, J., Erzurumluoglu, A. M., Elsworth, B. L., Kemp, J. P., Howe, L., Haycock, P. C., et al. (2017). LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis. Bioinformatics 33, 272–279. doi: 10.1093/bioinformatics/btw613
Zhou, J., and Parada, L. F. (2012). PTEN signaling in autism spectrum disorders. Curr. Opin. Neurobiol. 22, 873–879. doi: 10.1016/j.conb.2012.05.004
Keywords: epilepsy, GWAS, meta-analysis, TWAS, GRM3
Citation: Song M, Liu J, Yang Y, Lv L, Li W and Luo X-J (2021) Genome-Wide Meta-Analysis Identifies Two Novel Risk Loci for Epilepsy. Front. Neurosci. 15:722592. doi: 10.3389/fnins.2021.722592
Received: 09 June 2021; Accepted: 19 July 2021;
Published: 12 August 2021.
Edited by:
Rossen Donev, MicroPharm Ltd., United KingdomReviewed by:
Ciaran Campbell, Royal College of Surgeons in Ireland, IrelandJeffrey Dennis Calhoun, Northwestern University, United States
Copyright © 2021 Song, Liu, Yang, Lv, Li and Luo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wenqiang Li, bHdxNzgxNjAzQDE2My5jb20=; Xiong-Jian Luo, bHVveGlvbmdqaWFuQG1haWwua2l6LmFjLmNu
†These authors have contributed equally to this work