Within-sibship genome-wide association analyses decrease bias in estimates of direct genetic effects.

Estimates from genome-wide association studies (GWAS) of unrelated individuals capture effects of inherited variation (direct effects), demography (population stratification, assortative mating) and relatives (indirect genetic effects). Family-based GWAS designs can control for demographic and indirect genetic effects, but large-scale family datasets have been lacking. We combined data from 178,086 siblings from 19 cohorts to generate population (between-family) and within-sibship (within-family) GWAS estimates for 25 phenotypes. Within-sibship GWAS estimates were smaller than population estimates for height, educational attainment, age at first birth, number of children, cognitive ability, depressive symptoms and smoking. Some differences were observed in downstream SNP heritability, genetic correlations and Mendelian randomization analyses. For example, the within-sibship genetic correlation between educational attainment and body mass index attenuated towards zero. In contrast, analyses of most molecular phenotypes (for example, low-density lipoprotein-cholesterol) were generally consistent. We also found within-sibship evidence of polygenic adaptation on taller height. Here, we illustrate the importance of family-based GWAS data for phenotypes influenced by demographic and indirect genetic effects.

Identification of pleiotropy at the gene level between psychiatric disorders and related traits.

Major mental disorders are highly prevalent and make a substantial contribution to the global disease burden. It is known that mental disorders share clinical characteristics, and genome-wide association studies (GWASs) have recently provided evidence for shared genetic factors as well. Genetic overlaps are usually identified at the single-marker level. Here, we aimed to identify genetic overlaps at the gene level between 7 mental disorders (schizophrenia, autism spectrum disorder, major depressive disorder, anorexia nervosa, ADHD, bipolar disorder and anxiety), 8 brain morphometric traits, 2 cognitive traits (educational attainment and general cognitive function) and 9 personality traits (subjective well-being, depressive symptoms, neuroticism, extraversion, openness to experience, agreeableness and conscientiousness, children's aggressive behaviour, loneliness) based on publicly available GWASs. We performed systematic conditional regression analyses to identify independent signals and select loci associated with more than one trait. We identified 48 genes containing independent markers associated with several traits (pleiotropy at the gene level). We also report 9 genes with different markers that show independent associations with single traits (allelic heterogeneity). This study demonstrates that mental disorders and related traits do show pleiotropy at the gene level as well as the single-marker level. The identification of these genes might be important for prioritizing further deep genotyping, functional studies, or drug targeting.

Author Correction: Study of 300,486 individuals identifies 148 independent genetic loci influencing general cognitive function.

Christina M. Lill, who contributed to analysis of data, was inadvertently omitted from the author list in the originally published version of this article. This has now been corrected in both the PDF and HTML versions of the article.

Meta-analysis of Alzheimer's disease on 9,751 samples from Norway and IGAP study identifies four risk loci.

A large fraction of genetic risk factors for Alzheimer's Disease (AD) is still not identified, limiting the understanding of AD pathology and study of therapeutic targets. We conducted a genome-wide association study (GWAS) of AD cases and controls of European descent from the multi-center DemGene network across Norway and two independent European cohorts. In a two-stage process, we first performed a meta-analysis using GWAS results from 2,893 AD cases and 6,858 cognitively normal controls from Norway and 25,580 cases and 48,466 controls from the International Genomics of Alzheimer's Project (IGAP), denoted the discovery sample. Second, we selected the top hits (p < 1 × 10-6) from the discovery analysis for replication in an Icelandic cohort consisting of 5,341 cases and 110,008 controls. We identified a novel genomic region with genome-wide significant association with AD on chromosome 4 (combined analysis OR = 1.07, p = 2.48 x 10-8). This finding implicated HS3ST1, a gene expressed throughout the brain particularly in the cerebellar cortex. In addition, we identified IGHV1-68 in the discovery sample, previously not associated with AD. We also associated USP6NL/ECHDC3 and BZRAP1-AS1 to AD, confirming findings from a follow-up transethnic study. These new gene loci provide further evidence for AD as a polygenic disorder, and suggest new mechanistic pathways that warrant further investigation.

Genetic variation in 117 myelination-related genes in schizophrenia: Replication of association to lipid biosynthesis genes.

Schizophrenia is a serious psychotic disorder with high heritability. Several common genetic variants, rare copy number variants and ultra-rare gene-disrupting mutations have been linked to disease susceptibility, but there is still a large gap between the estimated and explained heritability. Since several studies have indicated brain myelination abnormalities in schizophrenia, we aimed to examine whether variants in myelination-related genes could be associated with risk for schizophrenia. We established a set of 117 myelination genes by database searches and manual curation. We used a combination of GWAS (SCZ_N = 35,476; CTRL_N = 46,839), exome chip (SCZ_N = 269; CTRL_N = 336) and exome sequencing data (SCZ_N = 2,527; CTRL_N = 2,536) from schizophrenia cases and healthy controls to examine common and rare variants. We found that a subset of lipid-related genes was nominally associated with schizophrenia (p = 0.037), but this signal did not survive multiple testing correction (FWER = 0.16) and was mainly driven by the SREBF1 and SREBF2 genes that have already been linked to schizophrenia. Further analysis demonstrated that the lowest nominal p-values were p = 0.0018 for a single common variant (rs8539) and p = 0.012 for burden of rare variants (LRP1 gene), but none of them survived multiple testing correction. Our findings suggest that variation in myelination-related genes is not a major risk factor for schizophrenia.

Recently evolved human-specific methylated regions are enriched in schizophrenia signals.

BACKGROUND: One explanation for the persistence of schizophrenia despite the reduced fertility of patients is that it is a by-product of recent human evolution. This hypothesis is supported by evidence suggesting that recently-evolved genomic regions in humans are involved in the genetic risk for schizophrenia. Using summary statistics from genome-wide association studies (GWAS) of schizophrenia and 11 other phenotypes, we tested for enrichment of association with GWAS traits in regions that have undergone methylation changes in the human lineage compared to Neanderthals and Denisovans, i.e. human-specific differentially methylated regions (DMRs). We used analytical tools that evaluate polygenic enrichment of a subset of genomic variants against all variants. RESULTS: Schizophrenia was the only trait in which DMR SNPs showed clear enrichment of association that passed the genome-wide significance threshold. The enrichment was not observed for Neanderthal or Denisovan DMRs. The enrichment seen in human DMRs is comparable to that for genomic regions tagged by Neanderthal Selective Sweep markers, and stronger than that for Human Accelerated Regions. The enrichment survives multiple testing performed through permutation (n = 10,000) and bootstrapping (n = 5000) in INRICH (p < 0.01). Some enrichment of association with height was observed at the gene level. CONCLUSIONS: Regions where DNA methylation modifications have changed during recent human evolution show enrichment of association with schizophrenia and possibly with height. Our study further supports the hypothesis that genetic variants conferring risk of schizophrenia co-occur in genomic regions that have changed as the human species evolved. Since methylation is an epigenetic mark, potentially mediated by environmental changes, our results also suggest that interaction with the environment might have contributed to that association.

Study of 300,486 individuals identifies 148 independent genetic loci influencing general cognitive function.

General cognitive function is a prominent and relatively stable human trait that is associated with many important life outcomes. We combine cognitive and genetic data from the CHARGE and COGENT consortia, and UK Biobank (total N = 300,486; age 16-102) and find 148 genome-wide significant independent loci (P < 5 × 10-8) associated with general cognitive function. Within the novel genetic loci are variants associated with neurodegenerative and neurodevelopmental disorders, physical and psychiatric illnesses, and brain structure. Gene-based analyses find 709 genes associated with general cognitive function. Expression levels across the cortex are associated with general cognitive function. Using polygenic scores, up to 4.3% of variance in general cognitive function is predicted in independent samples. We detect significant genetic overlap between general cognitive function, reaction time, and many health variables including eyesight, hypertension, and longevity. In conclusion we identify novel genetic loci and pathways contributing to the heritability of general cognitive function.

Analysis of the joint effect of SNPs to identify independent loci and allelic heterogeneity in schizophrenia GWAS data.

We have tested published methods for capturing allelic heterogeneity and identifying loci of joint effects to uncover more of the "hidden heritability" of schizophrenia (SCZ). We used two tools, cojo-GCTA and multi-SNP, to analyze meta-statistics from the latest genome-wide association study (GWAS) on SCZ by the Psychiatric Genomics Consortium (PGC). Stepwise regression on markers with p values <10-7 in cojo-GCTA identified 96 independent signals. Eighty-five passed the genome-wide significance threshold. Cross-validation of cojo-GCTA by CLUMP was 76%, i.e., 26 of the loci identified by the PGC using CLUMP were found to be dependent on another locus by cojo-GCTA. The overlap between cojo-GCTA and multi-SNP was better (up to 92%). Three markers reached genome-wide significance (5 × 10-8) in a joint effect model. In addition, two loci showed possible allelic heterogeneity within 1-Mb genomic regions, while CLUMP analysis had identified 16 such regions. Cojo-GCTA identified fewer independent loci than CLUMP and seems to be more conservative, probably because it accounts for long-range LD and interaction effects between markers. These findings also explain why fewer loci with possible allelic heterogeneity remained significant after cojo-GCTA analysis. With multi-SNP, 86 markers were selected at the threshold 10-7. Multi-SNP identifies fewer independent signals, due to splitting of the data and use of smaller samples. We recommend that cojo-GCTA and multi-SNP are used for post-GWAS analysis of all traits to call independent loci. We conclude that only a few loci in SCZ show joint effects or allelic heterogeneity, but this could be due to lack of power for that data set.

Novel genetic loci associated with hippocampal volume.

The hippocampal formation is a brain structure integrally involved in episodic memory, spatial navigation, cognition and stress responsiveness. Structural abnormalities in hippocampal volume and shape are found in several common neuropsychiatric disorders. To identify the genetic underpinnings of hippocampal structure here we perform a genome-wide association study (GWAS) of 33,536 individuals and discover six independent loci significantly associated with hippocampal volume, four of them novel. Of the novel loci, three lie within genes (ASTN2, DPP4 and MAST4) and one is found 200 kb upstream of SHH. A hippocampal subfield analysis shows that a locus within the MSRB3 gene shows evidence of a localized effect along the dentate gyrus, subiculum, CA1 and fissure. Further, we show that genetic variants associated with decreased hippocampal volume are also associated with increased risk for Alzheimer's disease (rg=-0.155). Our findings suggest novel biological pathways through which human genetic variation influences hippocampal volume and risk for neuropsychiatric illness.

A genetic association study of CSMD1 and CSMD2 with cognitive function.

The complement cascade plays a role in synaptic pruning and synaptic plasticity, which seem to be involved in cognitive functions and psychiatric disorders. Genetic variants in the closely related CSMD1 and CSMD2 genes, which are implicated in complement regulation, are associated with schizophrenia. Since patients with schizophrenia often show cognitive impairments, we tested whether variants in CSMD1 and CSMD2 are also associated with cognitive functions per se. We took a discovery-replication approach, using well-characterized Scandinavian cohorts. A total of 1637 SNPs in CSMD1 and 206 SNPs in CSMD2 were tested for association with cognitive functions in the NCNG sample (Norwegian Cognitive NeuroGenetics; n=670). Replication testing of SNPs with p-value<0.001 (7 in CSMD1 and 3 in CSMD2) was carried out in the TOP sample (Thematically Organized Psychosis; n=1025) and the BETULA sample (Betula Longitudinal Study on aging, memory and dementia; n=1742). Finally, we conducted a meta-analysis of these SNPs using all three samples. The previously identified schizophrenia marker in CSMD1 (SNP rs10503253) was also included. The strongest association was observed between the CSMD1 SNP rs2740931 and performance in immediate episodic memory (p-value=5×10-6, minor allele A, MAF 0.48-0.49, negative direction of effect). This association reached the study-wide significance level (p⩽1.2×10-5). SNP rs10503253 was not significantly associated with cognitive functions in our samples. In conclusion, we studied n=3437 individuals and found evidence that a variant in CSMD1 is associated with cognitive function. Additional studies of larger samples with cognitive phenotypes will be needed to further clarify the role of CSMD1 in cognitive phenotypes in health and disease.

Novel genetic loci underlying human intracranial volume identified through genome-wide association.

Intracranial volume reflects the maximally attained brain size during development, and remains stable with loss of tissue in late life. It is highly heritable, but the underlying genes remain largely undetermined. In a genome-wide association study of 32,438 adults, we discovered five previously unknown loci for intracranial volume and confirmed two known signals. Four of the loci were also associated with adult human stature, but these remained associated with intracranial volume after adjusting for height. We found a high genetic correlation with child head circumference (ρgenetic = 0.748), which indicates a similar genetic background and allowed us to identify four additional loci through meta-analysis (Ncombined = 37,345). Variants for intracranial volume were also related to childhood and adult cognitive function, and Parkinson's disease, and were enriched near genes involved in growth pathways, including PI3K-AKT signaling. These findings identify the biological underpinnings of intracranial volume and their link to physiological and pathological traits.

Genetics of structural connectivity and information processing in the brain.

Understanding the genetic factors underlying brain structural connectivity is a major challenge in imaging genetics. Here, we present results from genome-wide association studies (GWASs) of whole-brain white matter (WM) fractional anisotropy (FA), an index of microstructural coherence measured using diffusion tensor imaging. Data from independent GWASs of 355 Swedish and 250 Norwegian healthy adults were integrated by meta-analysis to enhance power. Complementary GWASs on behavioral data reflecting processing speed, which is related to microstructural properties of WM pathways, were performed and integrated with WM FA results via multimodal analysis to identify shared genetic associations. One locus on chromosome 17 (rs145994492) showed genome-wide significant association with WM FA (meta P value = 1.87 × 10-08). Suggestive associations (Meta P value <1 × 10-06) were observed for 12 loci, including one containing ZFPM2 (lowest meta P value = 7.44 × 10-08). This locus was also implicated in multimodal analysis of WM FA and processing speed (lowest Fisher P value = 8.56 × 10-07). ZFPM2 is relevant in specification of corticothalamic neurons during brain development. Analysis of SNPs associated with processing speed revealed association with a locus that included SSPO (lowest meta P value = 4.37 × 10-08), which has been linked to commissural axon growth. An intergenic SNP (rs183854424) 14 kb downstream of CSMD1, which is implicated in schizophrenia, showed suggestive evidence of association in the WM FA meta-analysis (meta P value = 1.43 × 10-07) and the multimodal analysis (Fisher P value = 1 × 10-07). These findings provide novel data on the genetics of WM pathways and processing speed, and highlight a role of ZFPM2 and CSMD1 in information processing in the brain.

Common variants in the ARC gene are not associated with cognitive abilities.

INTRODUCTION: The Activity-Regulated Cytoskeleton-associated (ARC) gene encodes a protein that is critical for the consolidation of synaptic plasticity and long-term memory formation. Given ARC's key role in synaptic plasticity, we hypothesized that genetic variations in ARC may contribute to interindividual variability in human cognitive abilities or to attention-deficit hyperactivity disorder (ADHD) susceptibility, where cognitive impairment often accompanies the disorder. METHODS: We tested whether ARC variants are associated with six measures of cognitive functioning in 670 healthy subjects in the Norwegian Cognitive NeuroGenetics (NCNG) by extracting data from its Genome-Wide Association Study (GWAS). In addition, the Swedish Betula sample of 1800 healthy subjects who underwent similar cognitive testing was also tested for association with 19 tag SNPs. RESULTS: No ARC variants show association at the study-wide level, but several markers show a trend toward association with human cognitive functions. We also tested for association between ARC SNPs and ADHD in a Norwegian sample of cases and controls, but found no significant associations. CONCLUSION: This study suggests that common genetic variants located in ARC do not account for variance in human cognitive abilities, though small effects cannot be ruled out.

Common genetic variants influence human subcortical brain structures.

The highly complex structure of the human brain is strongly shaped by genetic influences. Subcortical brain regions form circuits with cortical areas to coordinate movement, learning, memory and motivation, and altered circuits can lead to abnormal behaviour and disease. To investigate how common genetic variants affect the structure of these brain regions, here we conduct genome-wide association studies of the volumes of seven subcortical regions and the intracranial volume derived from magnetic resonance images of 30,717 individuals from 50 cohorts. We identify five novel genetic variants influencing the volumes of the putamen and caudate nucleus. We also find stronger evidence for three loci with previously established influences on hippocampal volume and intracranial volume. These variants show specific volumetric effects on brain structures rather than global effects across structures. The strongest effects were found for the putamen, where a novel intergenic locus with replicable influence on volume (rs945270; P = 1.08 × 10(-33); 0.52% variance explained) showed evidence of altering the expression of the KTN1 gene in both brain and blood tissue. Variants influencing putamen volume clustered near developmental genes that regulate apoptosis, axon guidance and vesicle transport. Identification of these genetic variants provides insight into the causes of variability in human brain development, and may help to determine mechanisms of neuropsychiatric dysfunction.

Lack of association of the rs1344706 ZNF804A variant with cognitive functions and DTI indices of white matter microstructure in two independent healthy populations.

The rs1344706 single nucleotide polymorphism within intron 2 of the ZNF804A gene is strongly associated with schizophrenia and bipolar disorder. This variant has also been associated in some studies with a range of cognitive and neuroimaging phenotypes, but several studies have reported no effect on the same phenotypes in other samples. Here, we genotyped 670 healthy adult Norwegian subjects and 1753 healthy adult Swedish subjects for rs1344706, and tested for associations with cognitive phenotypes including general intellectual abilities, memory functions and cognitive inhibition. We also tested whether rs1344706 is associated with white matter microstructural properties using diffusion tensor imaging (DTI) data from 250 to 340 of the Norwegian and Swedish subjects, respectively. Whole-brain voxel-wise statistical modeling of the effect of the ZNF804A variant on two DTI indices, fractional anisotropy (FA) and radial diffusivity (RD), was performed using tract-based spatial statistics (TBSS), and commonly reported effect sizes were calculated within several large-scale white matter pathways based on neuroanatomical atlases. No significant associations were found between rs1344706 and the cognitive traits or white matter microstructure. We conclude that the rs1344706 SNP has no significant effect on these phenotypes in our two reasonably powered samples.

A genetic deconstruction of neurocognitive traits in schizophrenia and bipolar disorder.

BACKGROUND: Impairments in cognitive functions are common in patients suffering from psychiatric disorders, such as schizophrenia and bipolar disorder. Cognitive traits have been proposed as useful for understanding the biological and genetic mechanisms implicated in cognitive function in healthy individuals and in the dysfunction observed in psychiatric disorders. METHODS: Sets of genes associated with a range of cognitive functions often impaired in schizophrenia and bipolar disorder were generated from a genome-wide association study (GWAS) on a sample comprising 670 healthy Norwegian adults who were phenotyped for a broad battery of cognitive tests. These gene sets were then tested for enrichment of association in GWASs of schizophrenia and bipolar disorder. The GWAS data was derived from three independent single-centre schizophrenia samples, three independent single-centre bipolar disorder samples, and the multi-centre schizophrenia and bipolar disorder samples from the Psychiatric Genomics Consortium. RESULTS: The strongest enrichments were observed for visuospatial attention and verbal abilities sets in bipolar disorder. Delayed verbal memory was also enriched in one sample of bipolar disorder. For schizophrenia, the strongest evidence of enrichment was observed for the sets of genes associated with performance in a colour-word interference test and for sets associated with memory learning slope. CONCLUSIONS: Our results are consistent with the increasing evidence that cognitive functions share genetic factors with schizophrenia and bipolar disorder. Our data provides evidence that genetic studies using polygenic and pleiotropic models can be used to link specific cognitive functions with psychiatric disorders.

DCLK1 variants are associated across schizophrenia and attention deficit/hyperactivity disorder.

Doublecortin and calmodulin like kinase 1 (DCLK1) is implicated in synaptic plasticity and neurodevelopment. Genetic variants in DCLK1 are associated with cognitive traits, specifically verbal memory and general cognition. We investigated the role of DCLK1 variants in three psychiatric disorders that have neuro-cognitive dysfunctions: schizophrenia (SCZ), bipolar affective disorder (BP) and attention deficit/hyperactivity disorder (ADHD). We mined six genome wide association studies (GWASs) that were available publically or through collaboration; three for BP, two for SCZ and one for ADHD. We also genotyped the DCLK1 region in additional samples of cases with SCZ, BP or ADHD and controls that had not been whole-genome typed. In total, 9895 subjects were analysed, including 5308 normal controls and 4,587 patients (1,125 with SCZ, 2,496 with BP and 966 with ADHD). Several DCLK1 variants were associated with disease phenotypes in the different samples. The main effect was observed for rs7989807 in intron 3, which was strongly associated with SCZ alone and even more so when cases with SCZ and ADHD were combined (P-value = 4 × 10(-5) and 4 × 10(-6), respectively). Associations were also observed with additional markers in intron 3 (combination of SCZ, ADHD and BP), intron 19 (SCZ+BP) and the 3'UTR (SCZ+BP). Our results suggest that genetic variants in DCLK1 are associated with SCZ and, to a lesser extent, with ADHD and BP. Interestingly the association is strongest when SCZ and ADHD are considered together, suggesting common genetic susceptibility. Given that DCLK1 variants were previously found to be associated with cognitive traits, these results are consistent with the role of DCLK1 in neurodevelopment and synaptic plasticity.

Linkage-disequilibrium-based binning affects the interpretation of GWASs.

Genome-wide association studies (GWASs) are critically dependent on detailed knowledge of the pattern of linkage disequilibrium (LD) in the human genome. GWASs generate lists of variants, usually SNPs, ranked according to the significance of their association to a trait. Downstream analyses generally focus on the gene or genes that are physically closest to these SNPs and ignore their LD profile with other SNPs. We have developed a flexible R package (LDsnpR) that efficiently assigns SNPs to genes on the basis of both their physical position and their pairwise LD with other SNPs. We used the positional-binning and LD-based-binning approaches to investigate whether including these "LD-based" SNPs would affect the interpretation of three published GWASs on bipolar affective disorder (BP) and of the imputed versions of two of these GWASs. We show how including LD can be important for interpreting and comparing GWASs. In the published, unimputed GWASs, LD-based binning effectively "recovered" 6.1%-8.3% of Ensembl-defined genes. It altered the ranks of the genes and resulted in nonnegligible differences between the lists of the top 2,000 genes emerging from the two binning approaches. It also improved the overall gene-based concordance between independent BP studies. In the imputed datasets, although the increases in coverage (>0.4%) and rank changes were more modest, even greater concordance between the studies was observed, attesting to the potential of LD-based binning on imputed data as well. Thus, ignoring LD can result in the misinterpretation of the GWAS findings and have an impact on subsequent genetic and functional studies.

Gene-based analysis of regionally enriched cortical genes in GWAS data sets of cognitive traits and psychiatric disorders.

BACKGROUND: Despite its estimated high heritability, the genetic architecture leading to differences in cognitive performance remains poorly understood. Different cortical regions play important roles in normal cognitive functioning and impairment. Recently, we reported on sets of regionally enriched genes in three different cortical areas (frontomedial, temporal and occipital cortices) of the adult rat brain. It has been suggested that genes preferentially, or specifically, expressed in one region or organ reflect functional specialisation. Employing a gene-based approach to the analysis, we used the regionally enriched cortical genes to mine a genome-wide association study (GWAS) of the Norwegian Cognitive NeuroGenetics (NCNG) sample of healthy adults for association to nine psychometric tests measures. In addition, we explored GWAS data sets for the serious psychiatric disorders schizophrenia (SCZ) (n = 3 samples) and bipolar affective disorder (BP) (n = 3 samples), to which cognitive impairment is linked. PRINCIPAL FINDINGS: At the single gene level, the temporal cortex enriched gene RAR-related orphan receptor B (RORB) showed the strongest overall association, namely to a test of verbal intelligence (Vocabulary, P = 7.7E-04). We also applied gene set enrichment analysis (GSEA) to test the candidate genes, as gene sets, for enrichment of association signal in the NCNG GWAS and in GWASs of BP and of SCZ. We found that genes differentially expressed in the temporal cortex showed a significant enrichment of association signal in a test measure of non-verbal intelligence (Reasoning) in the NCNG sample. CONCLUSION: Our gene-based approach suggests that RORB could be involved in verbal intelligence differences, while the genes enriched in the temporal cortex might be important to intellectual functions as measured by a test of reasoning in the healthy population. These findings warrant further replication in independent samples on cognitive traits.

Genetic profile of patients with epilepsy on first-line antiepileptic drugs and potential directions for personalized treatment.

BACKGROUND: The first-line antiepileptic drugs, although affordable and effective in the control of seizures, are associated with adverse drug effects, and there is large interindividual variability in the appropriate dose at which patients respond favorably. This variability may partly be explained by functional consequences of genetic polymorphisms in the drug-metabolizing enzymes, such as the CYP450 family, microsomal epoxide hydrolase and UDP-glucuronosyltransferases, drug transporters, mainly ATP-binding cassette transporters, and drug targets, including sodium channels. The purpose of this study was to determine the allele and genotype frequencies of such genetic variants in patients with epilepsy from North India administered first-line antiepileptic drugs, such as phenobarbitone, phenytoin, carbamazepine and valproic acid, and compare them with worldwide epilepsy populations. MATERIALS & METHODS: SNP screening of 19 functional variants from 12 genes in 392 patients with epilepsy was carried out, and the patients were classified with respect to the metabolizing rate of their drug-metabolizing enzymes, efflux rate of drug transporters and sensitivity of drug targets. RESULTS: A total of 16 SNPs were found to be polymorphic, and the allelic frequencies for these SNPs were in conformance with Hardy-Weinberg equilibrium. Among all the polymorphisms studied, functional variants from genes encoding CYP2C19, EPHX1, ABCB1 and SCN1A were highly polymorphic in North Indian epilepsy patients, and might account for differential drug response to first-line antiepileptic drugs. CONCLUSION: Interethnic differences were elucidated for several polymorphisms that might be responsible for differential serum drug levels and optimal dose requirement for efficacious treatment.