Computational dissection of human episodic memory reveals mental process-specific genetic profiles.

Episodic memory performance is the result of distinct mental processes, such as learning, memory maintenance, and emotional modulation of memory strength. Such processes can be effectively dissociated using computational models. Here we performed gene set enrichment analyses of model parameters estimated from the episodic memory performance of 1,765 healthy young adults. We report robust and replicated associations of the amine compound SLC (solute-carrier) transporters gene set with the learning rate, of the collagen formation and transmembrane receptor protein tyrosine kinase activity gene sets with the modulation of memory strength by negative emotional arousal, and of the L1 cell adhesion molecule (L1CAM) interactions gene set with the repetition-based memory improvement. Furthermore, in a large functional MRI sample of 795 subjects we found that the association between L1CAM interactions and memory maintenance revealed large clusters of differences in brain activity in frontal cortical areas. Our findings provide converging evidence that distinct genetic profiles underlie specific mental processes of human episodic memory. They also provide empirical support to previous theoretical and neurobiological studies linking specific neuromodulators to the learning rate and linking neural cell adhesion molecules to memory maintenance. Furthermore, our study suggests additional memory-related genetic pathways, which may contribute to a better understanding of the neurobiology of human memory.

Human genome-guided identification of memory-modulating drugs.

In the last decade there has been an exponential increase in knowledge about the genetic basis of complex human traits, including neuropsychiatric disorders. It is not clear, however, to what extent this knowledge can be used as a starting point for drug identification, one of the central hopes of the human genome project. The aim of the present study was to identify memory-modulating compounds through the use of human genetic information. We performed a multinational collaborative study, which included assessment of aversive memory--a trait central to posttraumatic stress disorder--and a gene-set analysis in healthy individuals. We identified 20 potential drug target genes in two genomewide-corrected gene sets: the neuroactive ligand-receptor interaction and the long-term depression gene set. In a subsequent double-blind, placebo-controlled study in healthy volunteers, we aimed at providing a proof of concept for the genome-guided identification of memory modulating compounds. Pharmacological intervention at the neuroactive ligand-receptor interaction gene set led to significant reduction of aversive memory. The findings demonstrate that genome information, along with appropriate data mining methodology, can be used as a starting point for the identification of memory-modulating compounds.

Epigenetic modification of the glucocorticoid receptor gene is linked to traumatic memory and post-traumatic stress disorder risk in genocide survivors.

Recent evidence suggests that altered expression and epigenetic modification of the glucocorticoid receptor gene (NR3C1) are related to the risk of post-traumatic stress disorder (PTSD). The underlying mechanisms, however, remain unknown. Because glucocorticoid receptor signaling is known to regulate emotional memory processes, particularly in men, epigenetic modifications of NR3C1 might affect the strength of traumatic memories. Here, we found that increased DNA methylation at the NGFI-A (nerve growth factor-induced protein A) binding site of the NR3C1 promoter was associated with less intrusive memory of the traumatic event and reduced PTSD risk in male, but not female survivors of the Rwandan genocide. NR3C1 methylation was not significantly related to hyperarousal or avoidance symptoms. We further investigated the relationship between NR3C1 methylation and memory functions in a neuroimaging study in healthy subjects. Increased NR3C1 methylation-which was associated with lower NR3C1 expression-was related to reduced picture recognition in male, but not female subjects. Furthermore, we found methylation-dependent differences in recognition memory-related brain activity in men. Together, these findings indicate that an epigenetic modification of the glucocorticoid receptor gene promoter is linked to interindividual and gender-specific differences in memory functions and PTSD risk.

PKCα is genetically linked to memory capacity in healthy subjects and to risk for posttraumatic stress disorder in genocide survivors.

Strong memory of a traumatic event is thought to contribute to the development and symptoms of posttraumatic stress disorder (PTSD). Therefore, a genetic predisposition to build strong memories could lead to increased risk for PTSD after a traumatic event. Here we show that genetic variability of the gene encoding PKCα (PRKCA) was associated with memory capacity--including aversive memory--in nontraumatized subjects of European descent. This finding was replicated in an independent sample of nontraumatized subjects, who additionally underwent functional magnetic resonance imaging (fMRI). fMRI analysis revealed PRKCA genotype-dependent brain activation differences during successful encoding of aversive information. Further, the identified genetic variant was also related to traumatic memory and to the risk for PTSD in heavily traumatized survivors of the Rwandan genocide. Our results indicate a role for PKCα in memory and suggest a genetic link between memory and the risk for PTSD.

Association of KIBRA with episodic and working memory: a meta-analysis.

WWC1 was first implicated in human cognition through a genome wide association study in 2006 that reported an association of the intronic single nucleotide polymorphism (SNP) rs17070145 with episodic memory performance. WWC1 encodes the protein KIBRA, which is almost ubiquitously expressed. Together with its binding partners, KIBRA is assumed to play a role in synaptic plasticity. T-allele carriers of SNP rs17070145 have been reported to outperform individuals that are homozygous for the C-allele in episodic memory tasks. Here we report two random effects meta-analyses testing the association of rs17070145 with episodic and working memory. All currently available population-based association studies that investigated effects of rs17070145 on episodic or working memory were included in the analyses. Where performance measures for multiple domain-specific tasks were available for a given study population, averaged effect size estimates were calculated. The performed meta-analyses relied on 17 samples that were tested for episodic memory performance (N = 8,909) and 9 samples that had performed working memory tasks (N = 4,696). We report a significant association of rs17070145 with both episodic (r = 0.068, P = 0.001) and working memory (r = 0.035, P = 0.018). In summary, our findings indicate that SNP rs17070145 located within KIBRA explains 0.5% of the variance for episodic memory tasks and 0.1% of the variance for working memory tasks in samples of primarily Caucasian background.

Rare variants in TMEM132D in a case-control sample for panic disorder.

Genome-wide association studies have identified common variants associated with common diseases. Most variants, however, explain only a small proportion of the estimated heritability, suggesting that rare variants might contribute to a larger extent to common diseases than assumed to date. Here, we use next-generation sequencing to test whether such variants contribute to the risk for anxiety disorders by re-sequencing 40 kb including all exons of the TMEM132D locus which we have previously shown to be associated with panic disorder and anxiety severity measures. DNA from 300 patients suffering from anxiety disorders, mostly panic disorder (84.7%), and 300 healthy controls was screened for the presence of genetic variants using next-generation re-sequencing in a pooled approach. Results were verified by individual re-genotyping. We identified 371 variants of which 247 had not been reported before, including 15 novel non-synonymous variants. The majority, 76% of these variants had a minor allele frequency less than 5%. While we did not identify additional common variants in TMEM132D associated with panic disorders, we observed an overrepresentation of presumably functional coding variants in healthy controls as compared to cases as well as a higher rate of private coding variants in cases, with one non-synonymous coding variant present in four patients but not in any of the matched controls nor in over 5,500 individuals of different ethnic origins from publicly available re-sequencing datasets. Our data suggest that not only common but also putatively functional and/or rare variants within TMEM132D might contribute to the risk to develop anxiety disorders.

Polymorphisms in the GAD2 gene-region are associated with susceptibility for unipolar depression and with a risk factor for anxiety disorders.

Glutamate decarboxylase (GAD) is the rate limiting enzyme for conversion of glutamic acid to gamma-aminobutyric acid (GABA). The GAD 65 kDa isoform is encoded by the gene GAD2 and is mainly expressed in synaptic terminals. It serves as an apoenzyme, which shows enhanced availability in situations of stress, responding to short-term demands for GABA. We analyzed 18 single nucleotide polymorphisms (SNPs) in the GAD2-gene region for associations with psychiatric diagnosis and behavioral inhibition (BI) derived from the personality traits neuroticism and extraversion as defined by the Eysenck Personality Questionaire (EPQ). A total of 268 patients with anxiety disorder (AD), 541 with unipolar depression (MD), and 541 healthy controls were included. We observe associations for five tag-SNPs with BI in the AD- and control samples as well as two additional case-control associations in the MD-sample. The associated SNPs lie within a 16KB linkage disequilibrium-block, including putative 5' GAD2-promoter-elements as well as the 3' end of the gene MYO3A. Using open access mRNA-expression data, we could show that BI-associated SNPs appear to be associated with differences in MYO3A- but not GAD2 lymphoblastoid-mRNA expression levels. These results support earlier studies that suggest associations of polymorphisms within the GAD2 locus with anxiety and affective disorders. However, data from expression studies imply that these polymorphisms could tag functional effects on the neighboring gene MYO3A, which is also expressed in the brain, including the cingulate cortex and the amygdala.

Polymorphisms in the angiotensin-converting enzyme gene region predict coping styles in healthy adults and depressed patients.

Dispositional coping styles are important moderators of the stress reaction and are altered in stress-related disorders like cardiovascular diseases and affective disorders. Heritability studies suggest a considerable genetic contribution to the interindividual variability in coping styles. Since the angiotensin-converting enzyme (ACE) gene has been described to be associated with the vulnerability for stress-related disorders and with altered stress hormone regulation, we investigated the ACE gene as potential candidate gene for coping styles. Five hundred forty one mentally healthy subjects and 194 patients suffering from depression participating in the Munich Antidepressant Response Signature (MARS) project were examined. Coping styles were assessed with a self-report questionnaire (German Stress Coping Questionnaire SVF78) measuring the individual coping style pattern in response to stressful situations. We genotyped 15 single nucleotide polymorphisms (SNPs) and the insertion/Deletion (I/D)-polymorphism in the ACE gene region and investigated their associations with coping styles. In healthy subjects, the highest association was observed between rs8066276, an intronic SNP of the ACE gene, and the coping factor Distraction. A further intronic SNP rs4305, not in linkage disequilibrium with rs8066276, showed an association with Devaluation/Defense. All associated copying styles can be categorized as potentially stress reducing factors (positive coping). Both SNPs were also found to be associated with positive coping styles in the patient sample; rs8066276 was associated with Devaluation/Defense, and rs4305 showed associations with Control. These results suggest that the ACE gene is involved in the development of coping strategies.

Identification of Two Distinct Working Memory-Related Brain Networks in Healthy Young Adults.

Working memory (WM) is an important cognitive domain for everyday life functioning and is often disturbed in neuropsychiatric disorders. Functional magnetic resonance imaging (fMRI) studies in humans show that distributed brain areas typically described as fronto-parietal regions are implicated in WM tasks. Based on data from a large sample of healthy young adults (N = 1369), we applied independent component analysis (ICA) to the WM-fMRI signal and identified two distinct networks that were relevant for differences in individual WM task performance. A parietally-centered network was particularly relevant for individual differences in task measures related to WM performance ("WM dependent") and a frontally-centered network was relevant for differences in attention-dependent task performance. Importantly, frontal areas that are typically considered as key regions for WM were either involved in both WM-dependent and attention-dependent performance, or in attention-dependent performance only. The networks identified here are provided as publicly available datasets. These networks can be applied in future studies to derive a low-dimensional representation of the overall WM brain activation.

Exhaustive search for epistatic effects on the human methylome.

Studies assessing the existence and magnitude of epistatic effects on complex human traits provide inconclusive results. The study of such effects is complicated by considerable increase in computational burden, model complexity, and model uncertainty, which in concert decrease model stability. An additional source introducing significant uncertainty with regard to the detection of robust epistasis is the biological distance between the genetic variation and the trait under study. Here we studied CpG methylation, a genetically complex molecular trait that is particularly close to genomic variation, and performed an exhaustive search for two-locus epistatic effects on the CpG-methylation signal in two cohorts of healthy young subjects. We detected robust epistatic effects for a small number of CpGs (N = 404). Our results indicate that epistatic effects explain only a minor part of variation in DNA-CpG methylation. Interestingly, these CpGs were more likely to be associated with gene-expression of nearby genes, as also shown by their overrepresentation in DNase I hypersensitivity sites and underrepresentation in CpG islands. Finally, gene ontology analysis showed a significant enrichment of these CpGs in pathways related to HPV-infection and cancer.

The NCAM1 gene set is linked to depressive symptoms and their brain structural correlates in healthy individuals.

Depressive symptoms exist on a continuum, the far end of which is found in depressive disorders. Utilizing the continuous spectrum of depressive symptoms may therefore contribute to the understanding of the biological underpinnings of depression. Gene set enrichment analysis (GSEA) is an important tool for the identification of gene groups linked to complex traits, and was applied in the present study on genome-wide association study (GWAS) data of depression scores and their brain-level structural correlates in healthy young individuals. On symptom level (i.e. depression scores), robust enrichment was identified for two gene sets: NCAM1 Interactions and Collagen Formation. Depression scores were also associated with decreased fractional anisotropy (FA) - a brain white matter property - within the forceps minor and the left superior temporal longitudinal fasciculus. Within each of these tracts, mean FA value of depression score-associated voxels was used as a phenotype in a subsequent GSEA. The NCAM1 Interactions gene set was significantly enriched in these tracts. By linking the NCAM1 Interactions gene set to depression scores and their structural brain correlates in healthy participants, the current study contributes to the understanding of the molecular underpinnings of depressive symptomatology.

Common epigenetic variation in a European population of mentally healthy young adults.

DNA methylation represents an important link between structural genetic variation and complex phenotypes. The study of genome-wide CpG methylation and its relation to traits relevant to psychiatry has become increasingly important. Here, we analyzed quality metrics of 394,043 CpG sites in two samples of 568 and 319 mentally healthy young adults. For 25% of all CpGs we observed medium to large common epigenetic variation. These CpGs were overrepresented in open sea and shore regions, as well as in intergenic regions. They also showed a strong enrichment of significant hits in association analyses. Furthermore, a significant proportion of common DNA methylation is at least partially genetically driven and thus may be observed similarly across tissues. These findings could be of particular relevance for studies of complex neuropsychiatric traits, which often rely on proxy tissues.

Genetic Analysis of Association Between Calcium Signaling and Hippocampal Activation, Memory Performance in the Young and Old, and Risk for Sporadic Alzheimer Disease.

IMPORTANCE: Human episodic memory performance is linked to the function of specific brain regions, including the hippocampus; declines as a result of increasing age; and is markedly disturbed in Alzheimer disease (AD), an age-associated neurodegenerative disorder that primarily affects the hippocampus. Exploring the molecular underpinnings of human episodic memory is key to the understanding of hippocampus-dependent cognitive physiology and pathophysiology. OBJECTIVE: To determine whether biologically defined groups of genes are enriched in episodic memory performance across age, memory encoding-related brain activity, and AD. DESIGN, SETTING, AND PARTICIPANTS: In this multicenter collaborative study, which began in August 2008 and is ongoing, gene set enrichment analysis was done by using primary and meta-analysis data from 57 968 participants. The Swiss cohorts consisted of 3043 healthy young adults assessed for episodic memory performance. In a subgroup (n = 1119) of one of these cohorts, functional magnetic resonance imaging was used to identify gene set-dependent differences in brain activity related to episodic memory. The German Study on Aging, Cognition, and Dementia in Primary Care Patients cohort consisted of 763 elderly participants without dementia who were assessed for episodic memory performance. The International Genomics of Alzheimer's Project case-control sample consisted of 54 162 participants (17 008 patients with sporadic AD and 37 154 control participants). Analyses were conducted between January 2014 and June 2015. Gene set enrichment analysis in all samples was done using genome-wide single-nucleotide polymorphism data. MAIN OUTCOMES AND MEASURES: Episodic memory performance in the Swiss cohort and German Study on Aging, Cognition, and Dementia in Primary Care Patients cohort was quantified by picture and verbal delayed free recall tasks. In the functional magnetic resonance imaging experiment, activation of the hippocampus during encoding of pictures served as the phenotype of interest. In the International Genomics of Alzheimer's Project sample, diagnosis of sporadic AD served as the phenotype of interest. RESULTS: In the discovery sample, we detected significant enrichment for genes constituting the calcium signaling pathway, especially those related to the elevation of cytosolic calcium (P = 2 × 10-4). This enrichment was replicated in 2 additional samples of healthy young individuals (P = .02 and .04, respectively) and a sample of healthy elderly participants (P = .004). Hippocampal activation (P = 4 × 10-4) and the risk for sporadic AD (P = .01) were also significantly enriched for genes related to the elevation of cytosolic calcium. CONCLUSIONS AND RELEVANCE: By detecting consistent significant enrichment in independent cohorts of young and elderly participants, this study identified that calcium signaling plays a central role in hippocampus-dependent human memory processes in cognitive health and disease, contributing to the understanding and potential treatment of hippocampus-dependent cognitive pathology.

BAIAP2 is related to emotional modulation of human memory strength.

Memory performance is the result of many distinct mental processes, such as memory encoding, forgetting, and modulation of memory strength by emotional arousal. These processes, which are subserved by partly distinct molecular profiles, are not always amenable to direct observation. Therefore, computational models can be used to make inferences about specific mental processes and to study their genetic underpinnings. Here we combined a computational model-based analysis of memory-related processes with high density genetic information derived from a genome-wide study in healthy young adults. After identifying the best-fitting model for a verbal memory task and estimating the best-fitting individual cognitive parameters, we found a common variant in the gene encoding the brain-specific angiogenesis inhibitor 1-associated protein 2 (BAIAP2) that was related to the model parameter reflecting modulation of verbal memory strength by negative valence. We also observed an association between the same genetic variant and a similar emotional modulation phenotype in a different population performing a picture memory task. Furthermore, using functional neuroimaging we found robust genotype-dependent differences in activity of the parahippocampal cortex that were specifically related to successful memory encoding of negative versus neutral information. Finally, we analyzed cortical gene expression data of 193 deceased subjects and detected significant BAIAP2 genotype-dependent differences in BAIAP2 mRNA levels. Our findings suggest that model-based dissociation of specific cognitive parameters can improve the understanding of genetic underpinnings of human learning and memory.

Converging genetic and functional brain imaging evidence links neuronal excitability to working memory, psychiatric disease, and brain activity.

Working memory, the capacity of actively maintaining task-relevant information during a cognitive task, is a heritable trait. Working memory deficits are characteristic for many psychiatric disorders. We performed genome-wide gene set enrichment analyses in multiple independent data sets of young and aged cognitively healthy subjects (n = 2,824) and in a large schizophrenia case-control sample (n = 32,143). The voltage-gated cation channel activity gene set, consisting of genes related to neuronal excitability, was robustly linked to performance in working memory-related tasks across ages and to schizophrenia. Functional brain imaging in 707 healthy participants linked this gene set also to working memory-related activity in the parietal cortex and the cerebellum. Gene set analyses may help to dissect the molecular underpinnings of cognitive dimensions, brain activity, and psychopathology.

The BclI polymorphism of the glucocorticoid receptor gene is associated with emotional memory performance in healthy individuals.

Glucocorticoids, stress hormones released from the adrenal cortex, are important players in the regulation of emotional memory. Specifically, in animals and in humans, glucocorticoids enhance memory consolidation of emotionally arousing experiences, but impair memory retrieval. These glucocorticoid actions are partly mediated by glucocorticoid receptors in the hippocampus, amygdala and prefrontal cortex, key brain regions for emotional memory. In a recent study in patients who underwent cardiac surgery, the BclI polymorphism of the glucocorticoid receptor gene (NR3C1) was associated with traumatic memories and posttraumatic stress disorder symptoms after intensive care therapy. Based on this finding, we investigated if the BclI polymorphism is also associated with emotional memory in healthy young subjects (N=841). We used a picture-learning task consisting of learning and recalling neutral and emotional photographs on two consecutive days. The BclI variant was associated with short-delay recall of emotional pictures on both days, with GG carriers showing increased emotional memory performance as compared to GC and CC carriers. We did not detect a genotype-dependent difference in recall performance for neutral pictures. These findings suggest that the Bcll polymorphism contributes to inter-individual differences in emotional memory also in healthy humans.

Possible associations of NTRK2 polymorphisms with antidepressant treatment outcome: findings from an extended tag SNP approach.

BACKGROUND: Data from clinical studies and results from animal models suggest an involvement of the neurotrophin system in the pathology of depression and antidepressant treatment response. Genetic variations within the genes coding for the brain-derived neurotrophic factor (BDNF) and its key receptor Trkb (NTRK2) may therefore influence the response to antidepressant treatment. METHODS: We performed a single and multi-marker association study with antidepressant treatment outcome in 398 depressed Caucasian inpatients participating in the Munich Antidepressant Response Signature (MARS) project. Two Caucasian replication samples (N = 249 and N = 247) were investigated, resulting in a total number of 894 patients. 18 tagging SNPs in the BDNF gene region and 64 tagging SNPs in the NTRK2 gene region were genotyped in the discovery sample; 16 nominally associated SNPs were tested in two replication samples. RESULTS: In the discovery analysis, 7 BDNF SNPs and 9 NTRK2 SNPs were nominally associated with treatment response. Three NTRK2 SNPs (rs10868223, rs1659412 and rs11140778) also showed associations in at least one replication sample and in the combined sample with the same direction of effects (Pcorr  = .018, Pcorr  = .015 and Pcorr  = .004, respectively). We observed an across-gene BDNF-NTRK2 SNP interaction for rs4923468 and rs1387926. No robust interaction of associated SNPs was found in an analysis of BDNF serum protein levels as a predictor for treatment outcome in a subset of 93 patients. CONCLUSIONS/LIMITATIONS: Although not all associations in the discovery analysis could be unambiguously replicated, the findings of the present study identified single nucleotide variations in the BDNF and NTRK2 genes that might be involved in antidepressant treatment outcome and that have not been previously reported in this context. These new variants need further validation in future association studies.

Polymorphisms within the metabotropic glutamate receptor 1 gene are associated with depression phenotypes.

OBJECTIVES: Glutamate has been implicated in the pathophysiology and treatment of mood disorders possibly by affecting the regulation of the hypothalamus-pituitary-adrenocortical (HPA) axis. Growing evidence suggests an important role of the metabotropic glutamate receptor 1 (mGlu1) in depression-related phenotypes. To test whether these findings can also be supported by human genetics data, we explored polymorphisms within the metabotropic glutamate receptor 1 gene (GRM1) for their association with unipolar depression (UPD) as well as with biological phenotypes of this disorder. METHODS: We first tested the association of 43 tag-SNPs covering the GRM1 locus with UPD in 350 patients and 370 matched controls. We then investigated the effects of the associated SNPs on hippocampal glutamate levels estimated using ¹H-MR-spectroscopy (¹H-MRS) and on endocrine measures from the combined dexamethasone-suppression/CRH stimulation (dex/CRH) test. RESULTS: Within the GRM1 locus, 22 SNPs showed nominally significant association with UPD, of which 6 withstood corrections for multiple testing (rs2268666 with best allelic p=7.0×10⁻5). Supportive evidence for an association with UPD was gained from a second independent sample with 904 patients and 1012 controls. Furthermore, patients homozygous for the non-risk genotypes showed reduced hippocampal glutamate levels as measured by ¹H-MRS, a more pronounced normalization of HPA-axis hyperactivity as well as a better antidepressant treatment outcome. CONCLUSIONS: These results suggest that the combination of genetic and biological markers may allow to subgroup patients into etiopathogenetically more relevant subcategories which could guide clinicians in their antidepressant treatment choices.

Further evidence for executive dysfunction in subjects with RLS from a non-clinical sample.

OBJECTIVE AND BACKGROUND: Previous studies exploring cognitive functioning in RLS have either relied on medication free subjects sampled within a clinical context or on subjects with RLS symptoms identified within population samples. However, in contrast to clinical samples, population studies so far have not excluded the use of antidepressants, hypnotics, or RLS relevant medication, and study subjects were exclusively older in age. We therefore report on cognitive functioning in predominantly middle-aged individuals with RLS symptoms sampled from the general population and free of mental disorders and of hypnotic, psychopharmacological, or RLS relevant medication. METHODS: Participants with RLS symptoms and individually matched controls were identified within the MARS control study, a non-clinical control group study of 550 participants between 18 and 75 years. Cognitive functioning was assessed with the Trail Making Test A and B and a computerized German version of the Wisconsin Card Sorting Test (WCST). Performance was compared between 41 participants with RLS and 133 controls, and between a subgroup of 10 participants with frequent RLS symptoms (≥ 2/week) and 36 matched controls. RESULTS: There was no difference in cognitive functioning for the complete group of participants with RLS and controls. However, participants with frequent RLS symptoms showed impaired performance in the WCST. CONCLUSION: The results of this study add to the evidence that executive functioning is impaired in individuals with frequent RLS.

Statistical epistasis and functional brain imaging support a role of voltage-gated potassium channels in human memory.

Despite the current progress in high-throughput, dense genome scans, a major portion of complex traits' heritability still remains unexplained, a phenomenon commonly termed "missing heritability." The negligence of analytical approaches accounting for gene-gene interaction effects, such as statistical epistasis, is probably central to this phenomenon. Here we performed a comprehensive two-way SNP interaction analysis of human episodic memory, which is a heritable complex trait, and focused on 120 genes known to show differential, memory-related expression patterns in rat hippocampus. Functional magnetic resonance imaging was also used to capture genotype-dependent differences in memory-related brain activity. A significant, episodic memory-related interaction between two markers located in potassium channel genes (KCNB2 and KCNH5) was observed (P(nominal combined)=0.000001). The epistatic interaction was robust, as it was significant in a screening (P(nominal)=0.0000012) and in a replication sample (P(nominal)=0.01). Finally, we found genotype-dependent activity differences in the parahippocampal gyrus (P(nominal)=0.001) supporting the behavioral genetics finding. Our results demonstrate the importance of analytical approaches that go beyond single marker statistics of complex traits.