A powerful and efficient two-stage method for detecting gene-to-gene interactions in GWAS.

For over a decade functional gene-to-gene interaction (epistasis) has been suspected to be a determinant in the "missing heritability" of complex traits. However, searching for epistasis on the genome-wide scale has been challenging due to the prohibitively large number of tests which result in a serious loss of statistical power as well as computational challenges. In this article, we propose a two-stage method applicable to existing case-control data sets, which aims to lessen both of these problems by pre-assessing whether a candidate pair of genetic loci is involved in epistasis before it is actually tested for interaction with respect to a complex phenotype. The pre-assessment is based on a two-locus genotype independence test performed in the sample of cases. Only the pairs of loci that exhibit non-equilibrium frequencies are analyzed via a logistic regression score test, thereby reducing the multiple testing burden. Since only the computationally simple independence tests are performed for all pairs of loci while the more demanding score tests are restricted to the most promising pairs, genome-wide association study (GWAS) for epistasis becomes feasible. By design our method provides strong control of the type I error. Its favourable power properties especially under the practically relevant misspecification of the interaction model are illustrated. Ready-to-use software is available. Using the method we analyzed Parkinson's disease in four cohorts and identified possible interactions within several SNP pairs in multiple cohorts.

Use of support vector machines for disease risk prediction in genome-wide association studies: concerns and opportunities.

The success of genome-wide association studies (GWAS) in deciphering the genetic architecture of complex diseases has fueled the expectations whether the individual risk can also be quantified based on the genetic architecture. So far, disease risk prediction based on top-validated single-nucleotide polymorphisms (SNPs) showed little predictive value. Here, we applied a support vector machine (SVM) to Parkinson disease (PD) and type 1 diabetes (T1D), to show that apart from magnitude of effect size of risk variants, heritability of the disease also plays an important role in disease risk prediction. Furthermore, we performed a simulation study to show the role of uncommon (frequency 1-5%) as well as rare variants (frequency <1%) in disease etiology of complex diseases. Using a cross-validation model, we were able to achieve predictions with an area under the receiver operating characteristic curve (AUC) of ~0.88 for T1D, highlighting the strong heritable component (∼90%). This is in contrast to PD, where we were unable to achieve a satisfactory prediction (AUC ~0.56; heritability ~38%). Our simulations showed that simultaneous inclusion of uncommon and rare variants in GWAS would eventually lead to feasible disease risk prediction for complex diseases such as PD. The used software is available at http://www.ra.cs.uni-tuebingen.de/software/MACLEAPS/.

Extended kindred with recessive late-onset Alzheimer disease maps to locus 8p22-p21.2: a genome-wide linkage analysis.

Late-onset Alzheimer disease (LOAD) is a complex genetic disorder. Although genes involved in early-onset forms were discovered more than a decade ago, LOAD research has only been able to point out small effect loci, with the exception of APOE. We mapped the gene predisposing to LOAD in an extended inbred family coming from a genetically isolated region (24 sampled individuals, 12 of whom are affected), completing a genome-wide screen with an Affymetrix10 K single nucleotide polymorphism microarray. Genotyping results were evaluated under model-dependent (dominant and recessive) and model-free analysis. We obtained a maximum nonparametric linkage score of 3.24 (P=0.00006) on chromosome 8p22-p21.2. The same genomic position also yielded the highest multipoint heterogeneity LOD (HLOD) under a recessive model (HLOD=3.04). When we compared the results of the model-dependent analysis, a higher score was obtained in the recessive model (3.04) than in the dominant model (1.0). This is a new locus identified in LOAD, in chromosome 8p22-p21.2 and encompassing several candidate genes, among them CLU and PPP3CC that were excluded by sequencing. The finding of a recessive model of inheritance, consistent with the assumption of inbreeding as a morbidity factor in this population, supports the notion of a role of recessive genes in LOAD.

Probing the roles of LRR RLK genes in Arabidopsis thaliana roots using a custom T-DNA insertion set.

Leucine-rich repeat receptor-like protein kinases (LRR RLKs) represent the largest group of Arabidopsis RLKs with approximately 235 members. A minority of these LRR RLKs have been assigned to diverse roles in development, pathogen resistance and hormone perception. Using a reverse genetics approach, a collection of homozygous T-DNA insertion lines for 69 root expressed LRR RLK genes was screened for root developmental defects and altered response after exposure to environmental, hormonal/chemical and abiotic stress. The obtained data demonstrate that LRR RLKs play a role in a wide variety of signal transduction pathways related to hormone and abiotic stress responses. The described collection of T-DNA insertion mutants provides a valuable tool for future research into the function of LRR RLK genes.

Transcriptional plasticity of a soil arthropod across different ecological conditions.

Ecological functional genomics, dealing with the responses of organisms to their natural environment is confronted with a complex pattern of variation and a large number of confounding environmental factors. For gene expression studies to provide meaningful information on conditions deviating from normal, a baseline or normal operating range (NOR) response needs to be established which indicates how an organism's transcriptome reacts to naturally varying ecological factors. Here we determine the transcriptional plasticity of a soil arthropod, Folsomia candida, exposed to various natural environments, as part of a first attempt in establishing such a NOR. Animals were exposed to 26 different field soils after which gene expression levels were measured. The main factor found to regulate gene expression was soil-type (sand or clay). Cell homeostasis and DNA replication were affected in collembolans exposed to sandy soil, indicating general stress. Multivariate analysis identified soil fertility as the main factor influencing gene expression. Regarding land-use, only forest soils showed an expression pattern deviating from the others. No significant effect of land-use, agricultural practice or soil type on fitness was observed, but arsenic concentration was negatively correlated with reproductive output. In conclusion, transcriptional responses remained within a limited range across the different land-uses but were significantly affected by soil-type. This may be caused by the contrasting soil physicochemical properties to which F. candida strongly responds. The broad range of conditions over which this soil-living detritivore is able to survive and reproduce, indicates a strategy of high plasticity, which comes with extensive gene expression regulation.

Do consanguineous parents of a child affected by an autosomal recessive disease have more DNA identical-by-descent than similarly-related parents with healthy offspring? Design of a case-control study.

BACKGROUND: The offspring of consanguineous relations have an increased risk of congenital/genetic disorders and early mortality. Consanguineous couples and their offspring account for approximately 10% of the global population. The increased risk for congenital/genetic disorders is most marked for autosomal recessive disorders and depends on the degree of relatedness of the parents. For children of first cousins the increased risk is 2-4%. For individual couples, however, the extra risk can vary from zero to 25% or higher, with only a minority of these couples having an increased risk of at least 25%. It is currently not possible to differentiate between high-and low-risk couples. The quantity of DNA identical-by-descent between couples with the same degree of relatedness shows a remarkable variation. Here we hypothesize that consanguineous partners with children affected by an autosomal recessive disease have more DNA identical-by-descent than similarly-related partners who have only healthy children. The aim of the study is thus to establish whether the amount of DNA identical-by-descent in consanguineous parents of children with an autosomal recessive disease is indeed different from its proportion in consanguineous parents who have healthy children only. METHODS/DESIGN: This project is designed as a case-control study. Cases are defined as consanguineous couples with one or more children with an autosomal recessive disorder and controls as consanguineous couples with at least three healthy children and no affected child. We aim to include 100 case couples and 100 control couples. Control couples are matched by restricting the search to the same family, clan or ethnic origin as the case couple. Genome-wide SNP arrays will be used to test our hypothesis. DISCUSSION: This study contains a new approach to risk assessment in consanguineous couples. There is no previous study on the amount of DNA identical-by-descent in consanguineous parents of affected children compared to the consanguineous parents of healthy children. If our hypothesis proves to be correct, further studies are needed to obtain different risk figure estimates for the different proportions of DNA identical-by-descent. With more precise information about their risk status, empowerment of couples can be improved when making reproductive decisions.

The role of the COMT Val(158)Met polymorphism in the phenotypic expression of obsessive-compulsive disorder.

Obsessive-Compulsive Disorder (OCD) is characterized by the presence of obsessions and compulsions, and shows considerable phenotypic variability. Family and twin studies have indicated a genetic component in the etiology of OCD, and the catechol-O-methyl transferase (COMT) gene is an important candidate gene for OCD. This study investigates the influence of the functional COMT Val158Met polymorphism on the phenotypic expression of OCD, using an item-level factor-analytic approach in a large sample. The COMT Val158Met variant was genotyped in 373 patients and 462 controls. It was tested whether there was an association between the COMT Val158Met polymorphism and OCD or dimensional phenotypes such as YBOCS severity score, age of onset of obsessive-compulsive symptoms and six symptom dimensions recently found in a large item-level factor-analytic study [Katerberg et al., submitted]. We further investigated possible sex-specific associations between the COMT Val158Met polymorphism and OCD or dimensional phenotypes. There was a trend for an association of the COMT 158Met allele with OCD in males, and an interaction between the COMT Val158Met genotype and sex on the somatic and sensory phenomena symptom dimension, with females showing lower scores. In conclusion, a dimensional approach seems fruitful in detecting genes of importance for OCD.

Adaptive differences in gene expression associated with heavy metal tolerance in the soil arthropod Orchesella cincta.

Field-selected tolerance to heavy metals has been reported for Orchesella cincta (Arthropoda: Collembola) populations occurring at metal-contaminated mining sites. This tolerance correlated with heritable increase in metal excretion efficiency, less pronounced cadmium (Cd)-induced growth reduction and overexpression of the metallothionein gene. We applied transcriptomics to determine differential gene expression caused by this abiotic stress in reference and Cd-tolerant populations. Many cDNAs responded to Cd exposure in the reference population. Significantly fewer clones were Cd responsive in tolerant animals. Analysis of variance revealed transcripts that interact between Cd exposure and population. Hierarchical cluster analysis of these clones identified two major groups. The first one contained cDNAs that were up-regulated by Cd in the reference culture but non-responsive or down-regulated in tolerant animals. This cluster was also characterized by elevated constitutive expression in the tolerant population. Gene ontology analysis revealed that these cDNAs were involved in structural integrity of the cuticle, anti-microbial defence, calcium channel-blocking, sulphur assimilation and chromatin remodelling. The second group consisted of cDNAs down-regulated in reference animals but not responding or slightly up-regulated in tolerant animals. Their functions involved carbohydrate metabolic processes, Ca(2+)-dependent stress signalling, redox state, proteolysis and digestion. The reference population showed a strong signature of stress-induced genome-wide perturbation of gene expression, whereas the tolerant animals maintained normal gene expression upon Cd exposure. We confirmed the micro-evolutionary processes occurring in soil arthropod populations and suggest a major contribution of gene regulation to the evolution of a stress-adapted phenotype.

A Functional polymorphism under positive evolutionary selection in ADRB2 is associated with human intelligence with opposite effects in the young and the elderly.

Comparative genomics offers a novel approach to unravel the genetic basis of complex traits. We performed a two stage analysis where genes ascertained for enhanced protein evolution in primates are subsequently searched for the presence of non-synonymous coding SNPs in the current human population at amino acid sites that differ between humans and chimpanzee. Positively selected genes among primates are generally presumed to determine phenotypic differences between humans and chimpanzee, such as the enhanced cognitive ability of our species. Amino acid substitutions segregating in humans at positively selected amino acid sites are expected to affect phenotypic differences among humans. Therefore we conducted an association study in two family based cohorts and one population based cohort between cognitive ability and the most likely candidate gene among the five that harbored more than one such polymorphism. The derived, human-specific allele of the beta-2 adrenergic receptor Arg16Gly polymorphism was found to be the increaser allele for performance IQ in the young, family based cohort but the decreaser allele for two different measures of cognition in the large Scottish cohort of unrelated individuals. The polymorphism is known to affect signaling activity and modulation of beta-2 adrenergic signaling has been shown to adjust memory consolidation, a trait related to cognition. The opposite effect of the polymorphism on cognition in the two age classes observed in the different cohorts resembles the effect of ADRB2 on hypertension, which also has been reported to be age dependent. This result illustrates the relevance of comparative genomics to detect genes that are involved in human behavior.

Glucocorticoid receptor gene polymorphisms and childhood adversity are associated with depression: New evidence for a gene-environment interaction.

The hypothalamic-pituitary-adrenal (HPA)-axis regulates the response to stressful events and is expected to be involved in the pathogenesis of depression. The glucocorticoid receptor (GR) regulates the activity of the HPA-axis. Both GR gene polymorphisms and childhood adversity are known to be associated with increased risk for depression. In the Longitudinal Aging Study Amsterdam, a large population based sample of older men and women, 906 subjects were genotyped. An association study was performed to determine the relationship between GR gene polymorphisms, childhood adversity, HPA-axis markers and depressive symptoms. A gene-environment interaction between the GR polymorphisms 22/23EK and 9beta and childhood adversity resulted in an increased risk of clinically relevant depressive symptoms. Without childhood adversity no increased risk was present. The 22/23EK variant was also associated with a lower Free Cortisol Index in the presence of childhood adversity. Persons that are heterozygous for the BclI variant, in contrast with wild-type and BclI-homozygotes, had lower serum levels of cortisol binding globulin and had no increased risk of recurrent depressive symptoms in the presence of childhood adversity. We found a gene-environment (G x E) interaction between common variants of the GR gene and childhood adversity, demonstrating a vulnerable phenotype for developing clinically relevant depressive symptoms at old age. This G x E interaction also influenced HPA-axis markers providing support for the involvement of the HPA-axis in both stress regulation and the pathogenesis of depression.

The role of the brain-derived neurotrophic factor (BDNF) val66met variant in the phenotypic expression of obsessive-compulsive disorder (OCD).

Evidence suggests that the Val66Met variant of the brain-derived neurotrophic factor (BDNF) gene may play a role in the etiology of Obsessive-Compulsive Disorder (OCD). In this study, the role of the BDNF Val66Met variant in the etiology and the phenotypic expression of OCD is investigated. Associations between the BDNF Val66Met variant and OCD, obsessive-compulsive symptom dimensions, Yale-Brown Obsessive Compulsive Scale (YBOCS) severity scores, age of onset and family history of obsessive-compulsive symptoms were assessed. The BDNF Val66Met variant was genotyped in 419 patients with sub-/clinical OCD and 650 controls. No differences in allele or genotype frequency were observed between cases and controls. In females with OCD, the Met66Met genotype was associated with later age of onset and a trend for a negative family history, whereas the Val66Val genotype was associated with a trend for lower YBOCS severity scores. Item-level factor analysis revealed six factors: 1) Contamination/cleaning; 2) Aggressive obsessions/checking; 3) Symmetry obsessions, counting, ordering and repeating; 4) Sexual/religious obsessions; 5) Hoarding and 6) Somatic obsessions/checking. A trend was found for a positive association between Factor 4 (Sexual/religious obsessions) and the BDNF Val66Val genotype. The results suggest that BDNF function may be implicated in the mediation of OCD. We found that for the BDNF Met66Met genotype may be associated with a milder phenotype in females and a possible role for the BDNF Val66Val genotype and the BDNF Val66 allele in the sexual/religious obsessions.

Empirical assessment of the validity of the 'fundamental theorem of the HapMap' in the light of 'cryptic' tagging of multiple susceptibility loci.

Underestimation of the sample size needed to detect genetic association may occur as a result of deviations from the 'fundamental theorem of the HapMap'. A biologically plausible mechanism that might cause this deviation is 'cryptic' tagging of multiple susceptibility loci by the same neutral marker. For complex disorders, the existence of multiple susceptibility loci on the same chromosome is probably the rule rather than the exception. Our results show that conditional on the known haplotype structure of the genome the probability that a tagging SNP that is in linkage disequilibrium (LD) with a susceptibility gene is also in LD with another susceptibility gene is not negligible. Consequently, we were able to estimate the extent and the prevalence of the bias in the necessary sample size to find association induced by 'cryptic' tagging. In general, the underestimation of the necessary sample size is modest: 5% of all association studies will underestimate the sample size by 5-30%. On the basis of our results, a safe bet is to use a sample that is 10% larger than otherwise deemed necessary.

Testing replication of a 5-SNP set for general cognitive ability in six population samples.

A 5-single nucleotide polymorphism (SNP) set has been associated with general cognitive ability in 5000 7-year-old children from the Twins Early Development Study (TEDS). Four of these SNPs were identified through a 10 K microarray analysis and one was identified through a targeted analysis of brain-expressed genes. The present study tested this association with general cognitive ability in six population samples of varying size and age from Australia, the UK (Scotland and England) and the Netherlands. Results from the largest sample (N=1310) approached significance (P=0.06) in the direction of the original finding, but results from the other samples (N=205-758) were mixed. A meta-analysis of the results--allowing for effect size heterogeneity between samples--yielded a non-significant correlation (r=-0.01, P=0.57), indicating that this SNP set was not associated with general cognitive ability in the populations studied.

Experimental evolution in Drosophila melanogaster: interaction of temperature and food quality selection regimes.

In Drosophila, both the phenotypic and evolutionary effect of temperature on adult size involves alterations to larval resource processing and affects other life-history traits, that is, development time but most notably, larval survival. Therefore, thermal evolution of adult body size might not be independent of simultaneous adaptation of larval traits to resource availability. Using experimental evolution lines adapted to high and low temperatures at different levels of food, we show that selection pressures interact in shaping larval resource processing. Evolution on poor food invariably leads to lower resource acquisition suggesting a cost to feeding behavior. However, following low temperature selection, lower resource acquisition led to a higher adult body size, probably by more efficient allocation to growth. In contrast, following high temperature selection, low resource acquisition benefited larval survival, possibly by reducing feeding-associated costs. We show that evolved differences to larval resource processing provide a possible proximate mechanism to variation in a suite of correlated life-history traits during adaptation to different climates. The implication for natural populations is that in nature, thermal evolution drives populations to opposite ends of an adult size versus larval survival trade-off by altering resource processing, if resource availability is limited.

Antagonistic pleiotropy for life-history traits at the gene expression level.

Life-history trade-offs prevent different components of fitness from being maximized simultaneously. Although the existence of trade-offs has been clearly demonstrated, the 'classical' mechanism of adaptive resource allocation that should underlie them has recently received criticism. In this study, we explore the molecular mechanisms of life-history trade-offs by applying a quantitative genomic approach. Analysis of global gene expression in Drosophila melanogaster revealed 34 genes whose expression coincided with the genetic trade-off between larval survival and adult size. The joint expression of these candidate 'trade-off' genes explained 86.3% of the trade-off. Fourteen of these genes have known functions which suggest that the larval survival-adult size trade-off could be the result of resource allocation at the organismal level, but at the level of cellular metabolism the trade-off would reduce to a shift between energy metabolism versus protein biosynthesis, regulated by the RAS signalling pathway.

Latitudinal clines in Drosophila melanogaster: body size, allozyme frequencies, inversion frequencies, and the insulin-signalling pathway.

Many latitudinal clines exist in Drosophila melanogaster: in adult body size, in allele frequency at allozyme loci, and in frequencies of common cosmopolitan inversions. The question is raised whether these latitudinal clines are causally related. This review aims to connect data from two very different fields of study, evolutionary biology and cell biology, in explaining such natural genetic variation in D. melanogaster body size and development time. It is argued that adult body size clines, inversion frequency clines, and clines in allele frequency at loci involved in glycolysis and glycogen storage are part of the same adaptive strategy. Selection pressure is expected to differ at opposite ends of the clines. At high latitudes, selection on D. melanogaster would favour high larval growth rate at low temperatures, and resource storage in adults to survive winter. At low latitudes selection would favour lower larval critical size to survive crowding, and increased male activity leading to high male reproductive success. Studies of the insulin-signalling pathway in D. melanogaster point to the involvement of this pathway in metabolism and adult body size. The genes involved in the insulin-signalling pathway are associated with common cosmopolitan inversions that show latitudinal clines. Each chromosome region connected with a large common cosmopolitan inversion possesses a gene of the insulin transmembrane complex, a gene of the intermediate pathway and a gene of the TOR branch. The hypothesis is presented that temperate D. melanogaster populations have a higher frequency of a 'thrifty' genotype corresponding to high insulin level or high signal level, while tropical populations possess a more 'spendthrift' genotype corresponding to low insulin or low signal level.

Accurate prediction of a minimal region around a genetic association signal that contains the causal variant.

In recent years, genome-wide association studies have been very successful in identifying loci for complex traits. However, typically these findings involve noncoding and/or intergenic SNPs without a clear functional effect that do not directly point to a gene. Hence, the challenge is to identify the causal variant responsible for the association signal. Typically, the first step is to identify all genetic variation in the locus region, usually by resequencing a large number of case chromosomes. Among all variants, the causal one needs to be identified in further functional studies. Because the experimental follow up can be very laborious, restricting the number of variants to be scrutinized can yield a great advantage. An objective method for choosing the size of the region to be followed up would be highly valuable. Here, we propose a simple method to call the minimal region around a significant association peak that is very likely to contain the causal variant. We model linkage disequilibrium (LD) in cases from the observed single SNP association signals, and predict the location of the causal variant by quantifying how well this relationship fits the data. Simulations showed that our approach identifies genomic regions of on average ∼50 kb with up to 90% probability to contain the causal variant. We apply our method to two genome-wide association data sets and localize both the functional variant REP1 in the α-synuclein gene that conveys susceptibility to Parkinson's disease and the APOE gene responsible for the association signal in the Alzheimer's disease data set.

Resequencing three candidate genes for major depressive disorder in a Dutch cohort.

Major depressive disorder (MDD) is a psychiatric disorder, characterized by periods of low mood of more than two weeks, loss of interest in normally enjoyable activities and behavioral changes. MDD is a complex disorder and does not have a single genetic cause. In 2009 a genome wide association study (GWAS) was performed on the Dutch GAIN-MDD cohort. Many of the top signals of this GWAS mapped to a region spanning the gene PCLO, and the non-synonymous coding single nucleotide polymorphism (SNP) rs2522833 in the PCLO gene became genome wide significant after post-hoc analysis. We performed resequencing of PCLO, GRM7, and SLC6A4 in 50 control samples from the GAIN-MDD cohort, to detect new genomic variants. Subsequently, we genotyped these variants in the entire GAIN-MDD cohort and performed association analysis to investigate if rs2522833 is the causal variant or simply in linkage disequilibrium with a more associated variant. GRM7 and SLC6A4 are both candidate genes for MDD from literature. We aimed to gather more evidence that rs2522833 is indeed the causal variant in the GAIN-MDD cohort or to find a previously undetected common variant in either PCLO, GRM7, or SLC6A4 with a higher association in this cohort. After next generation sequencing and association analysis we excluded the possibility of an undetected common variant to be more associated. For neither PCLO nor GRM7 we found a more associated variant. For SLC6A4, we found a new SNP that showed a lower P-value (P = 0.07) than in the GAIN-MDD GWAS (P = 0.09). However, no evidence for genome-wide significance was found. Although we did not take into account rare variants, we conclude that our results provide further support for the hypothesis that the non-synonymous coding SNP rs2522833 in the PCLO gene is indeed likely to be the causal variant in the GAIN-MDD cohort.

Cervical dystonia and genetic common variation in the dopamine pathway.

Cervical dystonia, a late onset focal dystonia, has a complex genetic background. Multiple lines of evidence point to a role for aberrant dopamine levels in dystonia. We assessed whether common variation within genes that regulate brain dopamine levels and in key genes of the dopamine metabolic pathway, modulate the risk for cervical dystonia. DNA was collected from 363 Dutch CD patients and a cohort of Dutch control individuals. Haplotype-tagging single nucleotide polymorphisms (SNPs) complemented with selected variants of functional importance in COMT, DAT, TH, MAO-A and -B, DDC and DBH were investigated. We tested the 143 markers in single-SNP, haplotype and epistasis analyses. We did not find an association with any of the selected 143 SNPs in these key dopamine genes. Our data shows that common variations in key genes of the dopamine pathway do not contribute to dystonia risk in the Dutch population. Possibly, risk alleles in this pathway may be rarer than detectable in this study, or might be located in downstream dopamine signaling pathway. Alternatively, found dopamine level changes are secondary to the dystonia disease processes.

Modulatory effects of the piccolo genotype on emotional memory in health and depression.

Major depressive disorder (MDD) has been associated with biased memory formation for mood-congruent information, which may be related to altered monoamine levels. The piccolo (PCLO) gene, involved in monoaminergic neurotransmission, has previously been linked to depression in a genome-wide association study. Here, we investigated the role of the PCLO risk allele on functional magnetic resonance imaging (MRI) correlates of emotional memory in a sample of 89 MDD patients (64 PCLO risk allele carriers) and 29 healthy controls (18 PCLO risk allele carriers). During negative word encoding, risk allele carriers showed significant lower activity relative to non-risk allele carriers in the insula, and trend-wise in the anterior cingulate cortex and inferior frontal gyrus. Moreover, depressed risk allele carriers showed significant lower activity relative to non-risk allele carriers in the striatum, an effect which was absent in healthy controls. Finally, amygdalar response during processing new positive words vs. known words was blunted in healthy PCLO+ carriers and in MDD patients irrespective of genotype, which may indicate that signalling of salient novel information does not occur to the same extent in PCLO+ carriers and MDD patients. The PCLO risk allele may increase vulnerability for MDD by modulating local brain function with regard to responsiveness to salient stimuli (i.e. insula) and processing novel negative information. Also, depression-specific effects of PCLO on dorsal striatal activation during negative word encoding and the absence of amygdalar salience signalling for novel positive information further suggest a role of PCLO in symptom maintenance in MDD.