Investigating the Genetics of Hippocampal Volume in Older Adults without Dementia

Hippocampal atrophy is observed with ageing and age-related neurodegenerative disease. Identification of the genetic correlates of hippocampal volume (HV) and atrophy may assist in elucidating the mechanisms of ageing and age-related neurodegeneration. Using two community cohorts of older Caucasians we estimated the heritability of HV and examined associations of HV with previously identified single nucleotide polymorphisms (SNPs). In addition we undertook genome-association studies (GWAS) examining HV and HV atrophy. Participants were community-dwelling non-demented older adults from the longitudinal Sydney Memory and Ageing Study (Sydney MAS) (N = 498) and the Older Australian Twins Study (OATS) (N = 351) aged 65 and over. HV was measured using T1-weighted magnetic resonance images. Heritability of HV was estimated in OATS. Genome-wide genotyping was imputed using the 1K Genomes reference set. Associations with HV-candidate and Alzheimer’s disease (AD)-related SNPs were investigated. A GWAS examining HV (in both cohorts) and an exploratory GWAS of HV atrophy over two years (in Sydney MAS only) were also undertaken. HV heritability was estimated at 62–65%. The previously identified GWAS HV SNP (rs6581612) and the candidate BDNF SNP (rs6265) were nominally significant (p = 0.047 and p = 0.041 respectively). No AD-related SNPs, including the APOE ε4 polymorphism, were significant. No significant results were observed for either of the GWAS undertaken. Despite our estimate of a high heritability of HV, our results are consistent with a highly polygenic model suggesting that SNPs identified from prior studies, including GWAS meta-analyses, can be difficult to replicate in smaller samples of older adults.     

Alzheimer’s Disease: Analyzing the Missing Heritability

Alzheimer’s disease (AD) is a complex disorder influenced by environmental and genetic factors. Recent work has identified 11 AD markers in 10 loci. We used Genome-wide Complex Trait Analysis to analyze >2 million SNPs for 10,922 individuals from the Alzheimer’s Disease Genetics Consortium to assess the phenotypic variance explained first by known late-onset AD loci, and then by all SNPs in the Alzheimer’s Disease Genetics Consortium dataset. In all, 33% of total phenotypic variance is explained by all common SNPs. APOE alone explained 6% and other known markers 2%, meaning more than 25% of phenotypic variance remains unexplained by known markers, but is tagged by common SNPs included on genotyping arrays or imputed with HapMap genotypes. Novel AD markers that explain large amounts of phenotypic variance are likely to be rare and unidentifiable using genome-wide association studies. Based on our findings and the current direction of human genetics research, we suggest specific study designs for future studies to identify the remaining heritability of Alzheimer’s disease.     

A Multi-Marker Genetic Association Test Based on the Rasch Model Applied to Alzheimer’s Disease

Results from Genome-Wide Association Studies (GWAS) have shown that the genetic basis of complex traits often include many genetic variants with small to moderate effects whose identification remains a challenging problem. In this context multi-marker analysis at the gene and pathway level can complement traditional point-wise approaches that treat the genetic markers individually. In this paper we propose a novel statistical approach for multi-marker analysis based on the Rasch model. The method summarizes the categorical genotypes of SNPs by a generalized logistic function into a genetic score that can be used for association analysis. Through different sets of simulations, the false-positive rate and power of the proposed approach are compared to a set of existing methods, and shows good performances. The application of the Rasch model on Alzheimer’s Disease (AD) ADNI GWAS dataset also allows a coherent interpretation of the results. Our analysis supports the idea that APOE is a major susceptibility gene for AD. In the top genes selected by proposed method, several could be functionally linked to AD. In particular, a pathway analysis of these genes also highlights the metabolism of cholesterol, that is known to play a key role in AD pathogenesis. Interestingly, many of these top genes can be integrated in a hypothetic signalling network.     

The Triangulation WIthin a STudy (TWIST) framework for causal inference within pharmacogenetic research

In this paper we review the methodological underpinnings of the general pharmacogenetic approach for uncovering genetically-driven treatment effect heterogeneity. This typically utilises only individuals who are treated and relies on fairly strong baseline assumptions to estimate what we term the ‘genetically moderated treatment effect’ (GMTE). When these assumptions are seriously violated, we show that a robust but less efficient estimate of the GMTE that incorporates information on the population of untreated individuals can instead be used. In cases of partial violation, we clarify when Mendelian randomization and a modified confounder adjustment method can also yield consistent estimates for the GMTE. A decision framework is then described to decide when a particular estimation strategy is most appropriate and how specific estimators can be combined to further improve efficiency. Triangulation of evidence from different data sources, each with their inherent biases and limitations, is becoming a well established principle for strengthening causal analysis. We call our framework ‘Triangulation WIthin a STudy’ (TWIST)’ in order to emphasise that an analysis in this spirit is also possible within a single data set, using causal estimates that are approximately uncorrelated, but reliant on different sets of assumptions. We illustrate these approaches by re-analysing primary-care-linked UK Biobank data relating to CYP2C19 genetic variants, Clopidogrel use and stroke risk, and data relating to APOE genetic variants, statin use and Coronary Artery Disease.     

A genome-wide association study identifies 5 loci associated with frozen shoulder and implicates diabetes as a causal risk factor

Frozen shoulder is a painful condition that often requires surgery and affects up to 5% of individuals aged 40–60 years. Little is known about the causes of the condition, but diabetes is a strong risk factor. To begin to understand the biological mechanisms involved, we aimed to identify genetic variants associated with frozen shoulder and to use Mendelian randomization to test the causal role of diabetes. We performed a genome-wide association study (GWAS) of frozen shoulder in the UK Biobank using data from 10,104 cases identified from inpatient, surgical and primary care codes. We used data from FinnGen for replication and meta-analysis. We used one-sample and two-sample Mendelian randomization approaches to test for a causal association of diabetes with frozen shoulder. We identified five genome-wide significant loci. The most significant locus (lead SNP rs28971325; OR = 1.20, [95% CI: 1.16–1.24], p = 5x10^-29) contained WNT7B. This variant was also associated with Dupuytren’s disease (OR = 2.31 [2.24, 2.39], p<1x10^-300) as were a further two of the frozen shoulder associated variants. The Mendelian randomization results provided evidence that type 1 diabetes is a causal risk factor for frozen shoulder (OR = 1.03 [1.02–1.05], p = 3x10^-6). There was no evidence that obesity was causally associated with frozen shoulder, suggesting that diabetes influences risk of the condition through glycemic rather than mechanical effects. We have identified genetic loci associated with frozen shoulder. There is a large overlap with Dupuytren’s disease associated loci. Diabetes is a likely causal risk factor. Our results provide evidence of biological mechanisms involved in this common painful condition.     

Genetic Dissection of the Drosophila melanogaster Female Head Transcriptome Reveals Widespread Allelic Heterogeneity

Modern genetic mapping is plagued by the “missing heritability” problem, which refers to the discordance between the estimated heritabilities of quantitative traits and the variance accounted for by mapped causative variants. One major potential explanation for the missing heritability is allelic heterogeneity, in which there are multiple causative variants at each causative gene with only a fraction having been identified. The majority of genome-wide association studies (GWAS) implicitly assume that a single SNP can explain all the variance for a causative locus. However, if allelic heterogeneity is prevalent, a substantial amount of genetic variance will remain unexplained. In this paper, we take a haplotype-based mapping approach and quantify the number of alleles segregating at each locus using a large set of 7922 eQTL contributing to regulatory variation in the Drosophila melanogaster female head. Not only does this study provide a comprehensive eQTL map for a major community genetic resource, the Drosophila Synthetic Population Resource, but it also provides a direct test of the allelic heterogeneity hypothesis. We find that 95% of cis-eQTLs and 78% of trans-eQTLs are due to multiple alleles, demonstrating that allelic heterogeneity is widespread in Drosophila eQTL. Allelic heterogeneity likely contributes significantly to the missing heritability problem common in GWAS studies.     

A One-Degree-of-Freedom Test for Supra-Multiplicativity of SNP Effects

Deviation from multiplicativity of genetic risk factors is biologically plausible and might explain why Genome-wide association studies (GWAS) so far could unravel only a portion of disease heritability. Still, evidence for SNP-SNP epistasis has rarely been reported, suggesting that 2-SNP models are overly simplistic. In this context, it was recently proposed that the genetic architecture of complex diseases could follow limiting pathway models. These models are defined by a critical risk allele load and imply multiple high-dimensional interactions. Here, we present a computationally efficient one-degree-of-freedom “supra-multiplicativity-test” (SMT) for SNP sets of size 2 to 500 that is designed to detect risk alleles whose joint effect is fortified when they occur together in the same individual. Via a simulation study we show that the SMT is powerful in the presence of threshold models, even when only about 30–45% of the model SNPs are available. In addition, we demonstrate that the SMT outperforms standard interaction analysis under recessive models involving just a few SNPs. We apply our test to 10 consensus Alzheimer’s disease (AD) susceptibility SNPs that were previously identified by GWAS and obtain evidence for supra-multiplicativity () that is not attributable to either two-way or three-way interaction.     

Association between APOC1 Polymorphism and Alzheimer’s Disease: A Case-Control Study and Meta-Analysis

Background

Previous association studies examining the relationship between the APOC1 polymorphism and susceptibility to Alzheimer’s disease (AD) have shown conflicting results, and it is not clear if an APOC1 variant acts as a genetic risk factor in AD etiology across multiple populations.

Methods

To confirm the risk association between APOC1 and AD, we designed a case-control study and also performed a meta-analysis of previously published studies.

Results

Seventy-nine patients with AD and one hundred fifty-six unrelated controls were included in case-control study. No association was found between the variation of APOC1 and AD in stage 1 of our study. However, our meta-analysis pooled a total of 2092 AD patients and 2685 controls. The APOC1 rs11568822 polymorphism was associated with increased AD risk in Caucasians, Asians and Caribbean Hispanics, but not in African Americans. APOE ε4 carriers harboring the APOC1 insertion allele, were more prevalent in AD patients than controls (χ² = 119.46, OR = 2.79, 95% CI = 2.31–3.36, P<0.01).

Conclusions

The APOC1 insertion allele, in combination with APOE ε4, likely serves as a potential risk factor for developing AD.     

Genetic Loci Associated with Alzheimer’s Disease and Cerebrospinal Fluid Biomarkers in a Finnish Case-Control Cohort

Objectives

To understand the relation between risk genes for Alzheimer’s disease (AD) and their influence on biomarkers for AD, we examined the association of AD in the Finnish cohort with single nucleotide polymorphisms (SNPs) from top AlzGene loci, genome-wide association studies (GWAS), and candidate gene studies; and tested the correlation between these SNPs and AD markers Aβ_1–42, total tau (t-tau), and phosphorylated tau (p-tau) in cerebrospinal fluid (CSF).

Methods

We tested 25 SNPs for genetic association with clinical AD in our cohort comprised of 890 AD patients and 701-age matched healthy controls using logistic regression. For the correlational study with biomarkers, we tested 36 SNPs in a subset of 222 AD patients with available CSF using mixed models. Statistical analyses were adjusted for age, gender and APOE status. False discovery rate for multiple testing was applied. All participants were from academic hospital and research institutions in Finland.

Results

APOE-ε4, CLU rs11136000, and MS4A4A rs2304933 correlated with significantly decreased Aβ_1–42 (corrected p<0.05). At an uncorrected p<0.05, PPP3R1 rs1868402 and MAPT rs2435211 were related with increased t-tau; while SORL1 rs73595277 and MAPT rs16940758, with increased p-tau. Only TOMM40 rs2075650 showed association with clinical AD after adjusting for APOE-ε4 (p = 0.007), but not after multiple test correction (p>0.05).

Conclusions

We provide evidence that APOE-ε4, CLU and MS4A4A, which have been identified in GWAS to be associated with AD, also significantly reduced CSF Aβ_1–42 in AD. None of the other AlzGene and GWAS loci showed significant effects on CSF tau. The effects of other SNPs on CSF biomarkers and clinical AD diagnosis did not reach statistical significance. Our findings suggest that APOE-ε4, CLU and MS4A4A influence both AD risk and CSF Aβ_1–42.     

The Nature of Genetic Variation for Complex Traits Revealed by GWAS and Regional Heritability Mapping Analyses

We use computer simulations to investigate the amount of genetic variation for complex traits that can be revealed by single-SNP genome-wide association studies (GWAS) or regional heritability mapping (RHM) analyses based on full genome sequence data or SNP chips. We model a large population subject to mutation, recombination, selection, and drift, assuming a pleiotropic model of mutations sampled from a bivariate distribution of effects of mutations on a quantitative trait and fitness. The pleiotropic model investigated, in contrast to previous models, implies that common mutations of large effect are responsible for most of the genetic variation for quantitative traits, except when the trait is fitness itself. We show that GWAS applied to the full sequence increases the number of QTL detected by as much as 50% compared to the number found with SNP chips but only modestly increases the amount of additive genetic variance explained. Even with full sequence data, the total amount of additive variance explained is generally below 50%. Using RHM on the full sequence data, a slightly larger number of QTL are detected than by GWAS if the same probability threshold is assumed, but these QTL explain a slightly smaller amount of genetic variance. Our results also suggest that most of the missing heritability is due to the inability to detect variants of moderate effect (∼0.03–0.3 phenotypic SDs) segregating at substantial frequencies. Very rare variants, which are more difficult to detect by GWAS, are expected to contribute little genetic variation, so their eventual detection is less relevant for resolving the missing heritability problem.     

Alzheimer’s Disease Risk Polymorphisms Regulate Gene Expression in the ZCWPW1 and the CELF1 Loci

Late onset Alzheimer’s disease (LOAD) is a genetically complex and clinically heterogeneous disease. Recent large-scale genome wide association studies (GWAS) have identified more than twenty loci that modify risk for AD. Despite the identification of these loci, little progress has been made in identifying the functional variants that explain the association with AD risk. Thus, we sought to determine whether the novel LOAD GWAS single nucleotide polymorphisms (SNPs) alter expression of LOAD GWAS genes and whether expression of these genes is altered in AD brains. The majority of LOAD GWAS SNPs occur in gene dense regions under large linkage disequilibrium (LD) blocks, making it unclear which gene(s) are modified by the SNP. Thus, we tested for brain expression quantitative trait loci (eQTLs) between LOAD GWAS SNPs and SNPs in high LD with the LOAD GWAS SNPs in all of the genes within the GWAS loci. We found a significant eQTL between rs1476679 and PILRB and GATS, which occurs within the ZCWPW1 locus. PILRB and GATS expression levels, within the ZCWPW1 locus, were also associated with AD status. Rs7120548 was associated with MTCH2 expression, which occurs within the CELF1 locus. Additionally, expression of several genes within the CELF1 locus, including MTCH2, were highly correlated with one another and were associated with AD status. We further demonstrate that PILRB, as well as other genes within the GWAS loci, are most highly expressed in microglia. These findings together with the function of PILRB as a DAP12 receptor supports the critical role of microglia and neuroinflammation in AD risk.     

Quantifying the underestimation of relative risks from genome-wide association studies

Genome-wide association studies (GWAS) have identified hundreds of associated loci across many common diseases. Most risk variants identified by GWAS will merely be tags for as-yet-unknown causal variants. It is therefore possible that identification of the causal variant, by fine mapping, will identify alleles with larger effects on genetic risk than those currently estimated from GWAS replication studies. We show that under plausible assumptions, whilst the majority of the per-allele relative risks (RR) estimated from GWAS data will be close to the true risk at the causal variant, some could be considerable underestimates. For example, for an estimated RR in the range 1.2-1.3, there is approximately a 38% chance that it exceeds 1.4 and a 10% chance that it is over 2. We show how these probabilities can vary depending on the true effects associated with low-frequency variants and on the minor allele frequency (MAF) of the most associated SNP. We investigate the consequences of the underestimation of effect sizes for predictions of an individual's disease risk and interpret our results for the design of fine mapping experiments. Although these effects mean that the amount of heritability explained by known GWAS loci is expected to be larger than current projections, this increase is likely to explain a relatively small amount of the so-called "missing" heritability.     

Fine mapping of five loci associated with low-density lipoprotein cholesterol detects variants that double the explained heritability

Complex trait genome-wide association studies (GWAS) provide an efficient strategy for evaluating large numbers of common variants in large numbers of individuals and for identifying trait-associated variants. Nevertheless, GWAS often leave much of the trait heritability unexplained. We hypothesized that some of this unexplained heritability might be due to common and rare variants that reside in GWAS identified loci but lack appropriate proxies in modern genotyping arrays. To assess this hypothesis, we re-examined 7 genes (APOE, APOC1, APOC2, SORT1, LDLR, APOB, and PCSK9) in 5 loci associated with low-density lipoprotein cholesterol (LDL-C) in multiple GWAS. For each gene, we first catalogued genetic variation by re-sequencing 256 Sardinian individuals with extreme LDL-C values. Next, we genotyped variants identified by us and by the 1000 Genomes Project (totaling 3,277 SNPs) in 5,524 volunteers. We found that in one locus (PCSK9) the GWAS signal could be explained by a previously described low-frequency variant and that in three loci (PCSK9, APOE, and LDLR) there were additional variants independently associated with LDL-C, including a novel and rare LDLR variant that seems specific to Sardinians. Overall, this more detailed assessment of SNP variation in these loci increased estimates of the heritability of LDL-C accounted for by these genes from 3.1% to 6.5%. All association signals and the heritability estimates were successfully confirmed in a sample of ～10,000 Finnish and Norwegian individuals. Our results thus suggest that focusing on variants accessible via GWAS can lead to clear underestimates of the trait heritability explained by a set of loci. Further, our results suggest that, as prelude to large-scale sequencing efforts, targeted re-sequencing efforts paired with large-scale genotyping will increase estimates of complex trait heritability explained by known loci.     

Partitioning the Heritability of Tourette Syndrome and Obsessive Compulsive Disorder Reveals Differences in Genetic Architecture

The direct estimation of heritability from genome-wide common variant data as implemented in the program Genome-wide Complex Trait Analysis (GCTA) has provided a means to quantify heritability attributable to all interrogated variants. We have quantified the variance in liability to disease explained by all SNPs for two phenotypically-related neurobehavioral disorders, obsessive-compulsive disorder (OCD) and Tourette Syndrome (TS), using GCTA. Our analysis yielded a heritability point estimate of 0.58 (se = 0.09, p = 5.64e-12) for TS, and 0.37 (se = 0.07, p = 1.5e-07) for OCD. In addition, we conducted multiple genomic partitioning analyses to identify genomic elements that concentrate this heritability. We examined genomic architectures of TS and OCD by chromosome, MAF bin, and functional annotations. In addition, we assessed heritability for early onset and adult onset OCD. Among other notable results, we found that SNPs with a minor allele frequency of less than 5% accounted for 21% of the TS heritability and 0% of the OCD heritability. Additionally, we identified a significant contribution to TS and OCD heritability by variants significantly associated with gene expression in two regions of the brain (parietal cortex and cerebellum) for which we had available expression quantitative trait loci (eQTLs). Finally we analyzed the genetic correlation between TS and OCD, revealing a genetic correlation of 0.41 (se = 0.15, p = 0.002). These results are very close to previous heritability estimates for TS and OCD based on twin and family studies, suggesting that very little, if any, heritability is truly missing (i.e., unassayed) from TS and OCD GWAS studies of common variation. The results also indicate that there is some genetic overlap between these two phenotypically-related neuropsychiatric disorders, but suggest that the two disorders have distinct genetic architectures.     

A genome‐wide association study of the frailty index highlights brain pathways in ageing

Frailty is a common geriatric syndrome and strongly associated with disability, mortality and hospitalization. Frailty is commonly measured using the frailty index (FI), based on the accumulation of a number of health deficits during the life course. The mechanisms underlying FI are multifactorial and not well understood, but a genetic basis has been suggested with heritability estimates between 30 and 45%. Understanding the genetic determinants and biological mechanisms underpinning FI may help to delay or even prevent frailty. We performed a genome‐wide association study (GWAS) meta‐analysis of a frailty index in European descent UK Biobank participants (n = 164,610, 60–70 years) and Swedish TwinGene participants (n = 10,616, 41–87 years). FI calculation was based on 49 or 44 self‐reported items on symptoms, disabilities and diagnosed diseases for UK Biobank and TwinGene, respectively. 14 loci were associated with the FI (p < 5*10⁻⁸). Many FI‐associated loci have established associations with traits such as body mass index, cardiovascular disease, smoking, HLA proteins, depression and neuroticism; however, one appears to be novel. The estimated single nucleotide polymorphism (SNP) heritability of the FI was 11% (0.11, SE 0.005). In enrichment analysis, genes expressed in the frontal cortex and hippocampus were significantly downregulated (adjusted p < 0.05). We also used Mendelian randomization to identify modifiable traits and exposures that may affect frailty risk, with a higher educational attainment genetic risk score being associated with a lower degree of frailty. Risk of frailty is influenced by many genetic factors, including well‐known disease risk factors and mental health, with particular emphasis on pathways in the brain.     

Effect of APOE and CD33 on Cognitive Decline

ObjectiveAn Alzheimer’s disease (AD) diagnosis is preceded by a long period of cognitive decline. We previously demonstrated increased risk of decline among individuals possessing one or more APOE ε4 alleles together with a family history of AD. The objective of this study is to investigate the possibility that such an increased risk might be due to AD risk genes with small effects in combination with APOE.MethodsParticipants in the Health and Retirement Study (HRS) over the age of 65, who contributed DNA, and had two or more evaluations with an abbreviated version of the modified Telephone Interview for Cognitive Status (TICS-m) were eligible for the study (n = 7451). A genetic risk score (g-score) was derived using AD risk genes’ meta-analyses data, assigning risk according to the number of risk alleles and summed over all the risk genes. Trajectories of cognitive function were modeled in four groups of Caucasian participants with and without one or more APOE ε4 alleles and either a high or low g-score: APOE ε4-/low g-score; APOE ε4-/high g-score; APOE ε4+/low g-score; and APOE ε4+/high g-score. Post hoc analyses evaluated interactions between individual genes and APOE.ResultsIndividuals in the APOE ε4+/high g-score group exhibited the greatest cognitive decline over time (p<.0001). This risk appeared to be greater than the sum of the effects of either high g-score or APOE ε4 alone. When gene interactions were individually tested with APOE, a statistically significant interaction with CD33 was discovered (p = 0.04) although the interaction was no longer significant when adjusted for multiple comparisons.ConclusionsIndividuals with multiple AD risk genes in addition to having one or more APOE ε4 alleles are at greater risk of cognitive decline than individuals with either APOE ε4 or a high genetic risk score. Among those with one or more APOE ε4 alleles, having one or more copies of the CD33 C (risk) allele may further increase the risk of cognitive decline.     

Promise and pitfalls of the Immunochip

Genomewide association studies (GWAS) have proven a powerful hypothesis-free method to identify common disease-associated variants. Even quite large GWAS, however, have only at best identified moderate proportions of the genetic variants contributing to disease heritability. To provide cost-effective genotyping of common and rare variants to map the remaining heritability and to fine-map established loci, the Immunochip Consortium has developed a 200,000 SNP chip that has been produced in very large numbers for a fraction of the cost of GWAS chips. This chip provides a powerful tool for immunogenetics gene mapping.     

Genetics of Callous-Unemotional Behavior in Children

Callous-unemotional behavior (CU) is currently under consideration as a subtyping index for conduct disorder diagnosis. Twin studies routinely estimate the heritability of CU as greater than 50%. It is now possible to estimate genetic influence using DNA alone from samples of unrelated individuals, not relying on the assumptions of the twin method. Here we use this new DNA method (implemented in a software package called Genome-wide Complex Trait Analysis, GCTA) for the first time to estimate genetic influence on CU. We also report the first genome-wide association (GWA) study of CU as a quantitative trait. We compare these DNA results to those from twin analyses using the same measure and the same community sample of 2,930 children rated by their teachers at ages 7, 9 and 12. GCTA estimates of heritability were near zero, even though twin analysis of CU in this sample confirmed the high heritability of CU reported in the literature, and even though GCTA estimates of heritability were substantial for cognitive and anthropological traits in this sample. No significant associations were found in GWA analysis, which, like GCTA, only detects additive effects of common DNA variants. The phrase ‘missing heritability’ was coined to refer to the gap between variance associated with DNA variants identified in GWA studies versus twin study heritability. However, GCTA heritability, not twin study heritability, is the ceiling for GWA studies because both GCTA and GWA are limited to the overall additive effects of common DNA variants, whereas twin studies are not. This GCTA ceiling is very low for CU in our study, despite its high twin study heritability estimate. The gap between GCTA and twin study heritabilities will make it challenging to identify genes responsible for the heritability of CU.     

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.