Epilepsy, hippocampal sclerosis and febrile seizures linked by common genetic variation around SCN1A

Epilepsy comprises several syndromes, amongst the most common being mesial temporal lobe epilepsy with hippocampal sclerosis. Seizures in mesial temporal lobe epilepsy with hippocampal sclerosis are typically drug-resistant, and mesial temporal lobe epilepsy with hippocampal sclerosis is frequently associated with important co-morbidities, mandating the search for better understanding and treatment. The cause of mesial temporal lobe epilepsy with hippocampal sclerosis is unknown, but there is an association with childhood febrile seizures. Several rarer epilepsies featuring febrile seizures are caused by mutations in SCN1A, which encodes a brain-expressed sodium channel subunit targeted by many anti-epileptic drugs. We undertook a genome-wide association study in 1018 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 7552 control subjects, with validation in an independent sample set comprising 959 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 3591 control subjects. To dissect out variants related to a history of febrile seizures, we tested cases with mesial temporal lobe epilepsy with hippocampal sclerosis with (overall n = 757) and without (overall n = 803) a history of febrile seizures. Meta-analysis revealed a genome-wide significant association for mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures at the sodium channel gene cluster on chromosome 2q24.3 [rs7587026, within an intron of the SCN1A gene, P = 3.36 × 10−9, odds ratio (A) = 1.42, 95% confidence interval: 1.26–1.59]. In a cohort of 172 individuals with febrile seizures, who did not develop epilepsy during prospective follow-up to age 13 years, and 6456 controls, no association was found for rs7587026 and febrile seizures. These findings suggest SCN1A involvement in a common epilepsy syndrome, give new direction to biological understanding of mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures, and open avenues for investigation of prognostic factors and possible prevention of epilepsy in some children with febrile seizures.

Epilepsy comprises several syndromes, amongst the most common being mesial temporal lobe epilepsy with hippocampal sclerosis. Seizures in mesial temporal lobe epilepsy with hippocampal sclerosis are typically drug-resistant, and mesial temporal lobe epilepsy with hippocampal sclerosis is frequently associated with important co-morbidities, mandating the search for better understanding and treatment. The cause of mesial temporal lobe epilepsy with hippocampal sclerosis is unknown, but there is an association with childhood febrile seizures. Several rarer epilepsies featuring febrile seizures are caused by mutations in SCN1A, which encodes a brain-expressed sodium channel subunit targeted by many anti-epileptic drugs. We undertook a genome-wide association study in 1018 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 7552 control subjects, with validation in an independent sample set comprising 959 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 3591 control subjects. To dissect out variants related to a history of febrile seizures, we tested cases with mesial temporal lobe epilepsy with hippocampal sclerosis with (overall n = 757) and without (overall n = 803) a history of febrile seizures. Meta-analysis revealed a genome-wide significant association for mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures at the sodium channel gene cluster on chromosome 2q24.3 [rs7587026, within an intron of the SCN1A gene, P = 3.36 Â 10 À 9 , odds ratio (A) = 1.42, 95% confidence interval: 1.26-1.59]. In a cohort of 172 individuals with febrile seizures, who did not develop epilepsy during prospective follow-up to age 13 years, and 6456 controls, no association was found for rs7587026 and febrile seizures. These findings suggest SCN1A involvement in a common epilepsy syndrome, give new direction to biological understanding of mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures, and open avenues for investigation of prognostic factors and possible prevention of epilepsy in some children with febrile seizures.

Introduction
Mesial temporal lobe epilepsy with hippocampal sclerosis (MTLEHS) is typically a serious epilepsy syndrome and the most common drug-resistant epilepsy (Berg et al., 2010). It is associated with burdensome co-morbidities, such as memory and psychiatric disorders. MTLEHS is the epilepsy most considered for therapeutic neurosurgery. Although surgery is a proven therapy, only $50% of patients have sustained postoperative seizure freedom (de Tisi et al., 2011), and surgery can have important adverse consequences. Better treatment options, or even prevention, of MTLEHS are therefore needed, but rational therapy for MTLEHS remains elusive because its causes are obscure (O'Dell et al., 2012).
MTLEHS is associated with a history of febrile seizures in childhood (Pittau et al., 2009;O'Dell et al., 2012). About 3% of children have febrile seizures; why only some go on to develop epilepsy, including MTLEHS, is unknown. There are a number of rare, genetically-determined, epilepsy syndromes in which febrile seizures are a prominent feature, such as Dravet syndrome and 'genetic epilepsy with febrile seizures plus' (GEFS + ) (Oliva et al., 2012). MTLEHS has rarely been described in families with GEFS + (Abou-Khalil et al., 2001) or familial febrile seizures (Mantegazza et al., 2005) associated with SCN1A mutations. In familial mesial temporal lobe epilepsy, some family members may have hippocampal sclerosis (Labate et al., 2011). A cluster of families with mesial temporal lobe epilepsy with hippocampal changes has been described in Brazil (Andrade-Valença et al., 2008). Together, this evidence implies genetic susceptibility to MTLEHS, although its heritability is unknown.
We hypothesized that MTLEHS, or MTLEHS with febrile seizures, as common epilepsy syndromes, might be associated with common genetic variation, and tested this 'common diseasecommon variant' hypothesis in a genetic association study.

Materials and methods
All aspects of the study were approved by the relevant institutional review board. All participants gave written informed consent.

Subjects
Patients were recruited during clinical appointments. MTLEHS was defined as in Wieser (2004). The diagnosis was made and/or reviewed by a consultant epileptologist who was part of this study, with access to history and investigation results. Patients with bilateral hippocampal sclerosis or dual pathology were excluded. One thousand and eighteen patients were included in the discovery stage and 959 patients in the replication. The number of patients by country is shown in Table 1, with further details in Supplementary Table 1. A history of presence or absence of febrile seizures was accepted only if contemporary medical records or a parental account was available; otherwise it was considered unknown, and not eligible for analysis. Population-based controls (n = 7552) were included in the discovery stage, and 3591 in the replication (Table 1 and Supplementary Table 1).
We also studied 542 individuals who had had febrile seizures but by the last follow-up had not had unprovoked seizures. These came from three groups: a German group; an Austrian group and the ALSPAC (Avon Longitudinal Study of Parents and Children) cohort, the latter followed to age 13 years (Supplementary material); MTLEHS after febrile seizures almost always develops by the age of 15 (Neligan et al., 2012). These cases were compared with 7387 control subjects from three relevant populations (Table 1). For the German and Austrian samples, the same controls as in the MTLEHS study were used.
To minimize population stratification, only individuals of white European ancestry were included. In the discovery stage, a combination of self-identified ancestry and EIGENSTRAT principal component methods was used to determine European ancestry. In the replication and febrile seizures analyses, only self-reported white individuals of European ancestry were included. More detailed ancestry data were available from all sources except Austria, allowing exclusion of individuals self-reported as coming from countries other than those where they were recruited.

Genotyping and quality control
In the discovery stage, all but the Austrian samples and Belgian controls comprised a subset of a previously described data set (Kasperaviciū te et al., 2010), genotyped on Illumina genome-wide genotyping chips, mostly on Illumina Human610-Quadv1/Human1-2M-DuoCustom. One hundred and fifty-seven Austrian patients and 332 controls were genotyped on Illumina HumanCNV370duo, and 285 Belgian controls were genotyped on Illumina HumanHap300 genotyping chips. Gender and relatedness checks were performed on all samples. The cluster plots of the top-associated single nucleotide polymorphisms were inspected manually. Details are given in Kasperaviciū te et al. (2010) and in the online Supplementary material. For replication analysis, several methods were used for genotyping.

Statistical analysis
In the discovery stage, genome-wide association analysis was performed using PLINK. Only single nucleotide polymorphisms present on both Illumina Human610-Quadv1 and Human1-2M-DuoCustom were analysed. In the discovery stage, we performed logistic regression using an additive model, including all significant EIGENSTRAT axes (assessed using the Tracy-Widom statistic with P 5 0.05) as covariates. Only single nucleotide polymorphisms with minor allele frequency of 5 1% were analysed. Since the replication samples did not have genome-wide data available to calculate EIGENSTRAT axes, we performed stratified analysis using the Cochran-Mantel-Haenszel test for 2 Â 2 Â 8 stratified case-control subsamples deriving from eight different recruitment countries and self-identified ancestry, using R. The Woolf test was used to assess effect heterogeneity. Meta-analysis of discovery and replication studies was performed using the inverse variance-weighted fixed-effects model as implemented in the GWAMA software (Mä gi and Morris, 2010). We considered an association to be genome-wide significant at P 5 5 Â 10 À 8 .
To fine map the association signal in the discovery stage, we imputed single nucleotide polymorphisms in the 10 Mb region surrounding rs7587026. Imputation was performed using MINIMAC (Howie et al., 2012), and 1000 Genomes Project data (1000Genomes Project Consortium et al., 2010 as the reference data set. Subsequent association analysis was performed using MACH2DAT (Li et al., 2010) using significant EIGENSTRAT axes as covariates.
Power calculations were performed using Genetic Power Calculator (Purcell et al., 2003).

Expression analysis
We tested association between genotypes of the two top single nucleotide polymorphisms rs7587026 and rs11692675 and SCN1A exons and gene expression in the middle temporal cortex (Brodmann areas 20 and 21) from 78 patients with MTLEHS who had undergone surgical resection, compared with 78 neurologically normal individuals from the MRC Sudden Death Brain and Tissue Bank. We specifically chose not to study the hippocampus to avoid confounding due to tissue changes such as cell loss and gliosis. All samples were randomly hybridized to Affymetrix Human Exon 1.0 ST arrays. Differential expression of SCN1A transcripts incorporating the 'neonatal' or 'adult' exon 5 form (5N or 5A exon, respectively), and expression of non-coding exons 1a and 1b (GenBank accession numbers DQ993522 and DQ993523, respectively) (Martin et al., 2007) in the 5' region of SCN1A, were tested by quantitative RT-PCR as they are not covered by the array. Details are provided in the Supplementary material.
Further, we tested whether the associated single nucleotide polymorphisms have an effect on expression or splicing of any genes in the genome in post-mortem tissue of nine brain regions from 134 control individuals (Supplementary material).

Genome-wide association analyses
We performed a two-stage study. For discovery, we first investigated genome-wide association between all MTLEHS and 531 164 single nucleotide polymorphisms in 1018 MTLEHS cases and 7552 controls from seven populations of European descent (Table 1 and  Supplementary Table 1). Using logistic regression analysis and correcting for population stratification, suggestive association emerged for three single nucleotide polymorphisms in a region of strong linkage disequilibrium on chromosome 2q24.3 encompassing SCN1A and other sodium channel genes ( Supplementary  Fig. 1). The most strongly associated single nucleotide polymorphism, rs11692675, is within intron 3 of the SCN1A full-length SCN1A encodes brain-expressed voltage-gated sodium channel type I, alpha subunit. It bears the largest number of known epilepsy-related mutations, some associated with febrile seizures (Oliva et al., 2012). The common SCN1A single nucleotide polymorphism rs3812718, affecting splicing (Heinzen et al., 2007), has also been associated with febrile seizures (Schlachter et al., 2009), though replication has failed (Petrovski et al., 2009). Retrospective studies show association between MTLE and febrile seizures (Pittau et al., 2009;O'Dell et al., 2012). Whether febrile seizures cause MTLEHS (Koyama et al., 2012) or whether pre-existing hippocampal abnormalities predispose to febrile seizures (Cendes, 2004), which may then also be injurious, is unknown. Clinical differences between patients with and without a history of febrile seizures suggest MTLEHS is heterogeneous (Thom et al., 2010). This evidence motivated our pre-analysis collection of febrile seizure data, and our previous study of febrile seizures (Petrovski et al., 2009). We performed analysis of patients in the discovery cohort with a known history of presence of childhood febrile seizures (MTLEHS + FS, n = 341) ( Table 2 and Fig. 1). The strongest association was for rs7587026, P = 2.64 Â 10 À 8 [OR(A) = 1.59, 95% CI: 1.35-1.87] and rs580041, P = 8.91 Â 10 À7 [OR(A) = 1.56, 95% CI: 1.33-1.84], whereas the signal for rs11692675 was slightly weaker, P = 1.25 Â 10 À 6 [OR(G) = 1.49, 95% CI: 1.27-1.75]. No association was seen in patients with MTLEHS without febrile seizures (MTLEHS À FS), despite similar sample size.
To refine the association signal, we performed regional imputation in the discovery data set of a 10 Mb region surrounding rs7587026 using the 1000 Genomes reference panel. Two single nucleotide polymorphisms had slightly lower P-values in the MTLEHS + FS analysis than the original single nucleotide polymorphisms [rs16851603 (P = 2.23 Â 10 À8 ) and rs3919196 (P = 2.26 Â 10 À8 )], but neither were significantly stronger than the original associations, and these signals reflected regional linkage disequilibrium structure ( Supplementary Fig. 2). No known functional variants in SCN1A, nor in other genes in the region, were in high linkage disequilibrium with rs7587026. The association signal is localized within one linkage disequilibrium block that also spans the promoter and 5' UTR region of SCN1A ( Supplementary Fig. 2).

Replication and combined analyses
We selected the two top single nucleotide polymorphisms, rs7587026 and rs11692675, for replication in an independent sample of 959 patients with MTLEHS, of whom 416 had MTLEHS + FS, and 3591 population-matched controls of European descent from eight populations (Table 1 and  Supplementary Table 1). We did not study rs580041 because of its perfect linkage disequilibrium with rs7587026 in white Europeans (r 2 = 1). We detected an association between rs7587026 and MTLEHS + FS, P = 5.88 Â 10 À3 [OR(A) = 1.26, 95% CI: 1.07-1.48; Table 2]; this value remains significant at a revised alpha threshold of 6.3 Â 10 À 3 after Bonferroni correction for multiple comparisons in the replication cohort. No association was detected for MTLEHSÀ FS. = 1.022. (C) Regional association results for the chromosome 2q24.3 locus. The left y-axis represents À log 10 (P-values) for association with MTLEHS, the right y-axis represents the recombination rate, and the x-axis represents base-pair positions along the chromosome (human genome Build 37). The top single nucleotide polymorphism, rs7587026, is shown in purple, the rest of the single nucleotide polymorphisms are coloured according to their linkage disequilibrium r 2 value with rs7587026.

Febrile seizure analysis
To explore whether the observed association with rs7587026 is specific for MTLEHS + FS, or is specific for febrile seizures in general, we examined a total of 7387 controls and three data sets of patients (totalling 542) who had had febrile seizures but had not developed epilepsy by the time of the latest follow-up. It has been shown that almost all children who go on to develop any type of epilepsy after febrile seizures have done so by the age of 15 years (Neligan et al., 2012). Therefore, to address the specificity of the association for MTLEHS + FS rather than febrile seizures alone, the ideal febrile seizures cohort would have been followed to age 15 at least. Of the three data sets available to us, only one met this criterion closely. The ALSPAC prospective cohort, which has the most comprehensive phenotypic data of the three data sets, followed children to age 13: there was no association of rs7587026 with febrile seizures in 171 individuals who did not go on to develop epilepsy (Table 3) in comparison with 6443 controls from the same cohort. The two other cohorts of children with febrile seizures, from Austria [samples partially overlapping with those reported in Schlachter et al. (2009)] and Germany, were ascertained at a young age (to six years of age only) and had no follow-up to establish whether the children had febrile seizures only, or febrile seizures in the context of subsequent epilepsy including MTLEHS, and are therefore not best suited to address the question, but were examined as few cohorts overall are available. Bearing this key caveat in mind, in these two data sets there was an observed association of febrile seizures with rs7587026 ( Table 3). The previously reported association in the Austrian population with an SCN1A functional splice-site single nucleotide polymorphism, rs3812718, was also seen in our Austrian sample [which is not unexpected as there is partial overlap of cases with those in the original report (Schlachter et al., 2009)] and was present in the German sample. The observed association of rs7587026 with febrile seizures disappeared in both Austrian and German data sets when analysis was conditioned on rs3812718 (P 4 0.19; Table 3). Moreover, although the association of febrile seizures with rs3812718 may be thought to be of interest for pure febrile seizures alone, we note there is no association of rs3812718 with febrile seizures in the best characterized cohort, from ALSPAC (Table 3), nor in a published sample (Petrovski et al., 2009).
Thus, although other single nucleotide polymorphisms in or near SCN1A may predispose to pure febrile seizures, the signal we observed in MTLEHS + FS is very unlikely to be due to the history of febrile seizures alone. Moreover, no significant association was detected in a group of patients with other partial epilepsies with a history of febrile seizures [data set from Kasperaviciute et al. The sample for this analysis was smaller (177 patients; 7552 controls), but had 81% power to detect association of OR 5 1.42 (as seen in MTLEHS + FS group combined analysis) under 0.05 significance level. Collectively, we found no evidence that the MTLEHS + FS association was due to febrile seizures, or that it holds for all partial epilepsies with febrile seizures.

SCN1A expression in the human brain
The observed association could act by modulating SCN1A gene expression. The associated region harbours several alternative untranslated SCN1A exons (Martin et al., 2007;Nakayama et al., 2010). We did not detect association between rs7587026 and any protein-coding exon except one (see below) or total SCN1A expression, or with expression of untranslated 5' exons 1a and 1b (Martin et al., 2007) (data not shown) in 78 patients and 78 control subjects.
The presence or absence of transcripts incorporating the 'neonatal' SCN1A exon 5 ('5N') was significantly different according to genotype of the two top single nucleotide polymorphisms (rs11692675 and rs7587026, P-values 1.08 Â 10 À9 and 1.17 Â 10 À6 , respectively; Supplementary material). For rs11692675 and rs7587026, respectively, none and 1% of the Figure 2 Forest plot for association of rs7587026 with MTLEHS + FS. The confidence interval for each study population is given by a horizontal line, and the point estimate is given by a square whose area is inversely proportional to the standard error of the estimate. The combined odds ratio is drawn as a diamond with horizontal limits at the confidence limits and width inversely proportional to its standard error. The study populations are ordered in descending order by the number of MTLEHS + FS cases.
individuals with the GG and AA genotype showed SCN1A transcripts in the neonatal form, compared with 83% and 81% with the genotype AA or CC. This alternative splicing event is influenced by rs3812718 (Heinzen et al., 2007). The association of alternative splicing with rs922224 (r 2 = 1 with rs3812718) was stronger, P = 2.33 Â 10 À31 . The level of expression of SCN1A exon 5N was also significantly different according to genotype (P = 1.62 Â 10 À11 for rs11692675, 2.70 Â 10 À 6 for rs7587026, 7.40 Â 10 À34 for rs922224). In conditional analyses including all three single nucleotide polymorphisms, only rs922224 remained significant (P = 1.08 Â 10 À25 ). Finally, expression quantitative trait loci analyses for subsets of patients according to a known history of presence (n = 46) or absence (n = 27) of febrile seizures in childhood for rs11692675 or rs7587026 showed significant differences in the level of expression of 5N exon according to genotype in both MTLEHS + FS and MTLEHS À FS. Including both rs11692675 or rs11692675 and rs922224 in the regression models, only rs922224 remained significant in both MTLEHS + FS and MTLEHS À FS groups (Supplementary material).
We cannot exclude the possibility that rs7587026 (or another single nucleotide polymorphism in the high linkage disequilibrium region) may act as an additional splicing controller to rs3812718, but our data are consistent with rs7587026 having no solo effect on 5N splicing. We also did not detect any correlation using a significance level of P 5 5 Â 10 À5 between rs7587026 and expression/splicing of any other genes across the genome (Supplementary material).

Discussion
We show that common variation in and near SCN1A may increase susceptibility to MTLEHS + FS. Our previously published larger genome-wide association study for a broader range of focal epilepsies did not identify any single-single nucleotide polymorphism association (Kasperaviciū te et al., 2010), but the findings here demonstrate that associated variants may exist for more narrowly-defined syndromes. Because the biology of most of the epilepsies is poorly understood, there are few a priori data upon which to base selection of the range of phenotypes to include in studies of possible genetic causation. Our findings suggest that focussing on clinically recognized syndromes or constellations (Berg et al., 2010) may prove fruitful by reducing heterogeneity before genomic analyses.
Our association seems to be specific for MTLEHS + FS, with no association for MTLEHS À FS, febrile seizures alone or non-MTLEHS partial epilepsies with febrile seizures. Our findings suggest that there is genetic susceptibility to MTLEHS, and that it, or hippocampal sclerosis, may not necessarily be only acquired. The results support the concept of heterogeneity in MTLEHS, beyond that already documented clinico-pathologically (Tassi et al., 2009;Thom et al., 2010;Blü mcke et al., 2012). However, further work will be needed to confirm the specificity of our findings, as we did not formally establish a significant difference in odds ratios between MTLEHS + FS and MTLEHS À FS. It would also be interesting to explore, in a suitably-powered study, whether there is any association with MTLE without hippocampal sclerosis.
The notably weaker association in the replication stage could be due to several factors, the most important of which is the 'winner's curse' (Ioannidis et al., 2009); there may be a large number of weak but real associations in the data, some of which achieve genome-wide significance in a particular study through random stochastic chance, but will not do so in another study. The association in our discovery cohort was replicated in the second independent sample, but it is nevertheless important that other studies are undertaken to further replicate our findings. Other limitations of our study are the lack of genome-wide data in the replication sample, preventing direct population stratification assessment, though self-identification closely corresponds to genetically-determined ancestry (Lao et al., 2008;Wang et al., 2010), a phenomenon we confirmed in the discovery stage, and the small size of some of the replication groups, reducing replication power, and magnifying effects of undetected population admixture. As for many genome-wide association studies, we could not narrow the association to a single gene or functional variant. There are other genes designated 'SCNxA' in the vicinity: SCN3A, SCN2A, SCN9A and SCN7A (this last does not show any sodium channel activity in exogenous expression systems) (Meisler et al., 2010). Among these genes, SCN2A has the most published evidence to support its role in the epilepsies. We cannot exclude the possibility that the association is driven by deleterious variants in these or other nearby genes. SCN1A, however, emerges as the most plausible candidate, due both to its proximity to the associated region and its role in other epilepsies with febrile seizures. Notably, our association is with a syndrome involving hippocampal damage, whereas typically no hippocampal damage is observed in patients with Dravet syndrome caused by deleterious changes affecting SCN1A (Catarino et al., 2011), suggesting that SCN1A might influence epileptogenesis through various mechanisms.
The location of the associated variants within SCN1A and overlapping its promoter regions (Long et al., 2008), was suggestive of possible roles in SCN1A expression modulation. In fact, we did not detect a definitive effect on expression of SCN1A or its exons in temporal neocortex. However, this analysis may have been confounded by many factors: effects may be brain-region or cellpopulation specific, as in SCN1A-related Dravet syndrome, where consequences are only found in interneurons (Ogiwara et al., 2007); our whole-tissue expression analysis would not detect such subtle signals. Moreover, noting the febrile seizures association, the effects may be temporally or spatially restricted, acting only in childhood or/and in the stress of febrile seizures (Koyama et al., 2012). Further studies will be needed to explore possible functional effects.
The detected association could act in different ways, predisposing to MTLEHS + FS as a distinct syndrome, or to the specific development of MTLEHS in the context of remote febrile seizures. If the association does indeed relate to SCN1A and function of the encoded protein, new lines of investigation may prove possible in the context of the existing deep knowledge of SCN1A, experimental models of MTLE and in vitro study of mechanisms of hippocampal dysfunction in epilepsy, as well as intriguing reports of the role of SCN1A in many epilepsies, such as the suggestion that mutations in SCN1A in Dravet Syndrome may protect against hippocampal sclerosis (Auvin et al., 2008;Catarino et al., 2011). Stratifying by febrile seizures type could also prove illuminating, as prolonged, lateralized or repeated febrile seizures within a short interval may have different effects to 'uncomplicated' febrile seizures. Our retrospective febrile seizures data were insufficiently resolved to permit such analysis. This is an important avenue for further investigation, because no predictors exist for the development of epilepsy in the 3% of all the children who have febrile seizures, and because established MTLEHS can have devastating consequences. Eventual reliable prediction of significant risk of MTLEHS after febrile seizures could lead to novel preventative measures in at-risk individuals: here, we note that SCN1A encodes an important anti-epileptic drug target and that it is possible to pharmacologically prevent the development of epilepsy after febrile seizures in an animal model (Koyama et al., 2012). Our findings suggest that further work on SCN1A variation may contribute to understanding the risk of developing MTLEHS after febrile seizures. 064154-01 to R.J.B. and H.H.); the collection of the Irish patient cohort was supported by the Irish Higher Education Authority Programme for Research in Third Level Institutions (PRTLI3) through a Science Foundation Ireland Research Frontiers Programme award (08/RFP/GEN1538) and a Medical Research Charities Group of Ireland/Health Research Board award (2009/ 001) from Brainwave-the Irish Epilepsy Association. GlaxoSmithKline funded the recruitment and phenotypic data collection of the GenEpA Consortium samples used in this study and contributed to the genotyping costs associated with their study. The collection of the Belgian patients was supported by the Fonds National de la Recherche Scientifique, grant n. FC 63574/ 3.4.620.06 F, and the Fonds Erasme, Université Libre de Bruxelles. The collection of the Belgian control cohort was supported by the Walloon Region and the French-Speaking Community of Belgium, the Belgian Science Policy and the University of Liè ge. J.N. is supported by the Swiss National Science Foundation-Fellowships for prospective researchers and the SICPA Foundation, Prilly, Switzerland. Computing facilities used at King's College London were supported by the National This study makes use of data generated by the Wellcome Trust Case-Control Consortium (a full list of the investigators who contributed to the generation of the data is available from www. wtccc.org.uk; funding for the project was provided by the Wellcome Trust under award 076113 and 085475) and the Irish Amyotrophic Lateral Sclerosis Study (funding support was provided by Muscular Dystrophy Association, USA, Irish Institute of Clinical Neurosciences Travel Award, and National Institutes of Health, USA; additional genotyping was provided by the National Institute of Neurological Disorders and Stroke (NINDS); the dataset used for the analyses described in this manuscript were obtained from the NINDS Database found at http://www.ncbi. nlm.nih.gov/gap through dbGaP accession number phs000127.v1.p1).