VarWatch-A stand-alone software tool for variant matching.

Massively parallel DNA sequencing of clinical samples holds great promise for the gene-based diagnosis of human inherited diseases because it allows rapid detection of putatively causative mutations at genome-wide level. Without additional evidence complementing their initial bioinformatics evaluation, however, the clinical relevance of such candidate genetic variants often remains unclear. In consequence, dedicated 'matching' services have been established in recent years that aim at the discovery of other, comparable case reports to facilitate individual diagnoses. However, legal concerns have been raised about the global sharing of genetic data, particularly in Europe where the recently enacted General Data Protection Regulation EU-2016/679 classifies genetic data as highly sensitive. Hence, unrestricted sharing of genetic data from clinical cases on platforms outside the national jurisdiction increasingly may be perceived as problematic. To allow collaborative data producers, particularly large consortia of diagnostic laboratories, to acknowledge these concerns while still practicing efficient case matching internally, novel tools are required. To this end, we developed VarWatch, an easy-to-deploy and highly scalable case matching software that provides users with comprehensive programmatic tools and a user-friendly interface to fulfil said purpose.

Novel multiple sclerosis susceptibility loci implicated in epigenetic regulation.

We conducted a genome-wide association study (GWAS) on multiple sclerosis (MS) susceptibility in German cohorts with 4888 cases and 10,395 controls. In addition to associations within the major histocompatibility complex (MHC) region, 15 non-MHC loci reached genome-wide significance. Four of these loci are novel MS susceptibility loci. They map to the genes L3MBTL3, MAZ, ERG, and SHMT1. The lead variant at SHMT1 was replicated in an independent Sardinian cohort. Products of the genes L3MBTL3, MAZ, and ERG play important roles in immune cell regulation. SHMT1 encodes a serine hydroxymethyltransferase catalyzing the transfer of a carbon unit to the folate cycle. This reaction is required for regulation of methylation homeostasis, which is important for establishment and maintenance of epigenetic signatures. Our GWAS approach in a defined population with limited genetic substructure detected associations not found in larger, more heterogeneous cohorts, thus providing new clues regarding MS pathogenesis.

Successful Replication of GWAS Hits for Multiple Sclerosis in 10,000 Germans Using the Exome Array.

Genome-wide association studies (GWAS) successfully identified various chromosomal regions to be associated with multiple sclerosis (MS). The primary aim of this study was to replicate reported associations from GWAS using an exome array in a large German study. German MS cases (n = 4,476) and German controls (n = 5,714) were genotyped using the Illumina HumanExome v1-Chip. Genotype calling was performed with the Illumina Genome Studio(TM) Genotyping Module, followed by zCall. Single-nucleotide polymorphisms (SNPs) in seven regions outside the human leukocyte antigen (HLA) region showed genome-wide significant associations with MS (P values < 5 × 10(-8) ). These associations have been reported previously. In addition, SNPs in three previously reported regions outside the HLA region yielded P values < 10(-5) . The effect of nine SNPs in the HLA region remained (P < 10(-5) ) after adjustment for other significant SNPs in the HLA region. All of these findings have been reported before or are driven by known risk loci. In summary, findings from previous GWAS for MS could be successfully replicated. We conclude that the regions identified in previous GWAS are also associated in the German population. This reassures the need for detailed investigations of the functional mechanisms underlying the replicated associations.

Connecting Anxiety and Genomic Copy Number Variation: A Genome-Wide Analysis in CD-1 Mice.

Genomic copy number variants (CNVs) have been implicated in multiple psychiatric disorders, but not much is known about their influence on anxiety disorders specifically. Using next-generation sequencing (NGS) and two additional array-based genotyping approaches, we detected CNVs in a mouse model consisting of two inbred mouse lines showing high (HAB) and low (LAB) anxiety-related behavior, respectively. An influence of CNVs on gene expression in the central (CeA) and basolateral (BLA) amygdala, paraventricular nucleus (PVN), and cingulate cortex (Cg) was shown by a two-proportion Z-test (p = 1.6 x 10-31), with a positive correlation in the CeA (p = 0.0062), PVN (p = 0.0046) and Cg (p = 0.0114), indicating a contribution of CNVs to the genetic predisposition to trait anxiety in the specific context of HAB/LAB mice. In order to confirm anxiety-relevant CNVs and corresponding genes in a second mouse model, we further examined CD-1 outbred mice. We revealed the distribution of CNVs by genotyping 64 CD 1 individuals using a high-density genotyping array (Jackson Laboratory). 78 genes within those CNVs were identified to show nominally significant association (48 genes), or a statistical trend in their association (30 genes) with the time animals spent on the open arms of the elevated plus-maze (EPM). Fifteen of them were considered promising candidate genes of anxiety-related behavior as we could show a significant overlap (permutation test, p = 0.0051) with genes within HAB/LAB CNVs. Thus, here we provide what is to our knowledge the first extensive catalogue of CNVs in CD-1 mice and potential corresponding candidate genes linked to anxiety-related behavior in mice.

A Candidate Gene Association Study Identifies DAPL1 as a Female-Specific Susceptibility Locus for Age-Related Macular Degeneration (AMD).

Age-related macular degeneration (AMD) is the leading cause of blindness among white caucasians over the age of 50 years with a prevalence rate expected to increase markedly with an anticipated increase in the life span of the world population. To further expand our knowledge of the genetic architecture of the disease, we pursued a candidate gene approach assessing 25 genes and a total of 109 variants. Of these, synonymous single nucleotide polymorphism (SNP) rs17810398 located in death-associated protein-like 1 (DAPL1) was found to be associated with AMD in a joint analysis of 3,229 cases and 2,835 controls from five studies [combined PADJ = 1.15 × 10(-6), OR 1.332 (1.187-1.496)]. This association was characterized by a highly significant sex difference (Pdiff = 0.0032) in that it was clearly confined to females with genome-wide significance [PADJ = 2.62 × 10(-8), OR 1.541 (1.324-1.796); males: PADJ = 0.382, OR 1.084 (0.905-1.298)]. By targeted resequencing of risk and non-risk associated haplotypes in the DAPL1 locus, we identified additional potentially functional risk variants, namely a common 897-bp deletion and a SNP predicted to affect a putative binding site of an exonic splicing enhancer. We show that the risk haplotype correlates with a reduced retinal transcript level of two, less frequent, non-canonical DAPL1 isoforms. DAPL1 plays a role in epithelial differentiation and may be involved in apoptotic processes thereby suggesting a possible novel pathway in AMD pathogenesis.

Genome-wide meta-analysis of myopia and hyperopia provides evidence for replication of 11 loci.

Refractive error (RE) is a complex, multifactorial disorder characterized by a mismatch between the optical power of the eye and its axial length that causes object images to be focused off the retina. The two major subtypes of RE are myopia (nearsightedness) and hyperopia (farsightedness), which represent opposite ends of the distribution of the quantitative measure of spherical refraction. We performed a fixed effects meta-analysis of genome-wide association results of myopia and hyperopia from 9 studies of European-derived populations: AREDS, KORA, FES, OGP-Talana, MESA, RSI, RSII, RSIII and ERF. One genome-wide significant region was observed for myopia, corresponding to a previously identified myopia locus on 8q12 (p = 1.25×10(-8)), which has been reported by Kiefer et al. as significantly associated with myopia age at onset and Verhoeven et al. as significantly associated to mean spherical-equivalent (MSE) refractive error. We observed two genome-wide significant associations with hyperopia. These regions overlapped with loci on 15q14 (minimum p value = 9.11×10(-11)) and 8q12 (minimum p value 1.82×10(-11)) previously reported for MSE and myopia age at onset. We also used an intermarker linkage- disequilibrium-based method for calculating the effective number of tests in targeted regional replication analyses. We analyzed myopia (which represents the closest phenotype in our data to the one used by Kiefer et al.) and showed replication of 10 additional loci associated with myopia previously reported by Kiefer et al. This is the first replication of these loci using myopia as the trait under analysis. "Replication-level" association was also seen between hyperopia and 12 of Kiefer et al.'s published loci. For the loci that show evidence of association to both myopia and hyperopia, the estimated effect of the risk alleles were in opposite directions for the two traits. This suggests that these loci are important contributors to variation of refractive error across the distribution.

Genetic relationship between five psychiatric disorders estimated from genome-wide SNPs.

Most psychiatric disorders are moderately to highly heritable. The degree to which genetic variation is unique to individual disorders or shared across disorders is unclear. To examine shared genetic etiology, we use genome-wide genotype data from the Psychiatric Genomics Consortium (PGC) for cases and controls in schizophrenia, bipolar disorder, major depressive disorder, autism spectrum disorders (ASD) and attention-deficit/hyperactivity disorder (ADHD). We apply univariate and bivariate methods for the estimation of genetic variation within and covariation between disorders. SNPs explained 17-29% of the variance in liability. The genetic correlation calculated using common SNPs was high between schizophrenia and bipolar disorder (0.68 ± 0.04 s.e.), moderate between schizophrenia and major depressive disorder (0.43 ± 0.06 s.e.), bipolar disorder and major depressive disorder (0.47 ± 0.06 s.e.), and ADHD and major depressive disorder (0.32 ± 0.07 s.e.), low between schizophrenia and ASD (0.16 ± 0.06 s.e.) and non-significant for other pairs of disorders as well as between psychiatric disorders and the negative control of Crohn's disease. This empirical evidence of shared genetic etiology for psychiatric disorders can inform nosology and encourages the investigation of common pathophysiologies for related disorders.

Meta-analysis of genome-wide association studies in five cohorts reveals common variants in RBFOX1, a regulator of tissue-specific splicing, associated with refractive error.

Visual refractive errors (REs) are complex genetic traits with a largely unknown etiology. To date, genome-wide association studies (GWASs) of moderate size have identified several novel risk markers for RE, measured here as mean spherical equivalent (MSE). We performed a GWAS using a total of 7280 samples from five cohorts: the Age-Related Eye Disease Study (AREDS); the KORA study ('Cooperative Health Research in the Region of Augsburg'); the Framingham Eye Study (FES); the Ogliastra Genetic Park-Talana (OGP-Talana) Study and the Multiethnic Study of Atherosclerosis (MESA). Genotyping was performed on Illumina and Affymetrix platforms with additional markers imputed to the HapMap II reference panel. We identified a new genome-wide significant locus on chromosome 16 (rs10500355, P = 3.9 × 10(-9)) in a combined discovery and replication set (26 953 samples). This single nucleotide polymorphism (SNP) is located within the RBFOX1 gene which is a neuron-specific splicing factor regulating a wide range of alternative splicing events implicated in neuronal development and maturation, including transcription factors, other splicing factors and synaptic proteins.

A mega-analysis of genome-wide association studies for major depressive disorder.

Prior genome-wide association studies (GWAS) of major depressive disorder (MDD) have met with limited success. We sought to increase statistical power to detect disease loci by conducting a GWAS mega-analysis for MDD. In the MDD discovery phase, we analyzed more than 1.2 million autosomal and X chromosome single-nucleotide polymorphisms (SNPs) in 18 759 independent and unrelated subjects of recent European ancestry (9240 MDD cases and 9519 controls). In the MDD replication phase, we evaluated 554 SNPs in independent samples (6783 MDD cases and 50 695 controls). We also conducted a cross-disorder meta-analysis using 819 autosomal SNPs with P<0.0001 for either MDD or the Psychiatric GWAS Consortium bipolar disorder (BIP) mega-analysis (9238 MDD cases/8039 controls and 6998 BIP cases/7775 controls). No SNPs achieved genome-wide significance in the MDD discovery phase, the MDD replication phase or in pre-planned secondary analyses (by sex, recurrent MDD, recurrent early-onset MDD, age of onset, pre-pubertal onset MDD or typical-like MDD from a latent class analyses of the MDD criteria). In the MDD-bipolar cross-disorder analysis, 15 SNPs exceeded genome-wide significance (P<5 × 10(-8)), and all were in a 248 kb interval of high LD on 3p21.1 (chr3:52 425 083-53 822 102, minimum P=5.9 × 10(-9) at rs2535629). Although this is the largest genome-wide analysis of MDD yet conducted, its high prevalence means that the sample is still underpowered to detect genetic effects typical for complex traits. Therefore, we were unable to identify robust and replicable findings. We discuss what this means for genetic research for MDD. The 3p21.1 MDD-BIP finding should be interpreted with caution as the most significant SNP did not replicate in MDD samples, and genotyping in independent samples will be needed to resolve its status.

Estimating the age of the most common Italian GRN mutation: walking back to Canossa times.

Mutations in the progranulin gene (GRN) were first implicated in frontotemporal lobar degeneration in 2006. The GRN p.Leu271LeufsX10 mutation is one of the most common GRN mutations worldwide. To gain further insight into the origin of this mutation in Italy, we performed a haplotype sharing analysis (32 families, residents of Lombardy) and refined the GRN p.Leu271LeufsX10 mutation dating. We showed that almost all families (30/32) can be traced to a single founder. We further estimated the age of this mutation using different methods and population growth rates both for Italy and Lombardy. Using DMLE, we dated the origin of this mutation to the Middle Ages, at the turn of the first millennium (phased families only, Italy: 39 and Lombardy: 32 generations ago; all families Italy: 45 and Lombardy 38 generations ago). Mutation dating was slightly postdated using Estiage (phased families only: 15 generations ago; all families: 20 generation ago). From a translational perspective, targeting mutation carriers offers a unique model to test disease-modifying drugs in clinical trials.

Resistance to antidepressant treatment is associated with polymorphisms in the leptin gene, decreased leptin mRNA expression, and decreased leptin serum levels.

Leptin, a peptide hormone from adipose tissue and key player in weight regulation, has been suggested to be involved in sleep and cognition and to exert antidepressant-like effects, presumably via its action on the HPA-axis and hippocampal function. This led us to investigate whether genetic variants in the leptin gene, the level of leptin mRNA-expression and leptin serum concentrations are associated with response to antidepressant treatment. Our sample consisted of inpatients from the Munich Antidepressant Response Signature (MARS) project with weekly Hamilton Depression ratings, divided into two subsamples. In the exploratory sample (n=251) 17 single nucleotide polymorphisms (SNPs) covering the leptin gene region were genotyped. We found significant associations of several SNPs with impaired antidepressant treatment outcome and impaired cognitive performance after correction for multiple testing. The SNP (rs10487506) showing the highest association with treatment response (p=3.9×10(-5)) was analyzed in the replication sample (n=358) and the association could be verified (p=0.021) with response to tricyclic antidepressants. In an additional meta-analysis combining results from the MARS study with data from the Genome-based Therapeutic Drugs for Depression (GENDEP) and the Sequenced Treatment Alternatives to Relieve Depression (STAR(⁎)D) studies, nominal associations of several polymorphisms in the upstream vicinity of rs10487506 with treatment outcome were detected (p=0.001). In addition, we determined leptin mRNA expression in lymphocytes and leptin serum levels in subsamples of the MARS study. Unfavorable treatment outcome was accompanied with decreased leptin mRNA and leptin serum levels. Our results suggest an involvement of leptin in antidepressant action and cognitive function in depression with genetic polymorphisms in the leptin gene, decreased leptin gene expression and leptin deficiency in serum being risk factors for resistance to antidepressant therapy in depressed patients.

The endocrine stress response is linked to one specific locus on chromosome 3 in a mouse model based on extremes in trait anxiety.

BACKGROUND: The hypothalamic-pituitary-adrenal (HPA) axis is essential to control physiological stress responses in mammals. Its dysfunction is related to several mental disorders, including anxiety and depression. The aim of this study was to identify genetic loci underlying the endocrine regulation of the HPA axis. METHOD: High (HAB) and low (LAB) anxiety-related behaviour mice were established by selective inbreeding of outbred CD-1 mice to model extremes in trait anxiety. Additionally, HAB vs. LAB mice exhibit comorbid characteristics including a differential corticosterone response upon stress exposure. We crossbred HAB and LAB lines to create F1 and F2 offspring. To identify the contribution of the endocrine phenotypes to the total phenotypic variance, we examined multiple behavioural paradigms together with corticosterone secretion-based phenotypes in F2 mice by principal component analysis. Further, to pinpoint the genomic loci of the quantitative trait of the HPA axis stress response, we conducted genome-wide multipoint oligogenic linkage analyses based on Bayesian Markov chain Monte Carlo approach as well as parametric linkage in three-generation pedigrees, followed by a two-dimensional scan for epistasis and association analysis in freely segregating F2 mice using 267 single-nucleotide polymorphisms (SNPs), which were identified to consistently differ between HAB and LAB mice as genetic markers. RESULTS: HPA axis reactivity measurements and behavioural phenotypes were represented by independent principal components and demonstrated no correlation. Based on this finding, we identified one single quantitative trait locus (QTL) on chromosome 3 showing a very strong evidence for linkage (2ln (L-score) > 10, LOD > 23) and significant association (lowest Bonferroni adjusted p < 10-28) to the neuroendocrine stress response. The location of the linkage peak was estimated at 42.3 cM (95% confidence interval: 41.3 - 43.3 cM) and was shown to be in epistasis (p-adjusted < 0.004) with the locus at 35.3 cM on the same chromosome. The QTL harbours genes involved in steroid synthesis and cardiovascular effects. CONCLUSION: The very prominent effect on stress-induced corticosterone secretion of the genomic locus on chromosome 3 and its involvement in epistasis highlights the critical role of this specific locus in the regulation of the HPA axis.

Apoptotic cleavage of DNA in human lymphocyte chromatin shows high sequence specificity.

Apoptotic digestion of human lymphocyte chromatin results in the appearance of large amounts of nucleosome size DNA fragments. Sequencing of these fragments and analysis of the distribution of bases around the apoptotic nucleases' cutting sites revealed a rather strong consensus sequence, not observed earlier. The consensus TAAAgTAcTTTA is characterized by complementary symmetry, resembling prokaryotic restriction sites. This consensus also possesses three TA dinucleotide steps, separated by five bases (corresponding to a half-period of the DNA double helix), suggesting strong bending of the DNA at the cut sites which is perhaps required for cutting.

Genotyping of fanconi anemia patients by whole exome sequencing: advantages and challenges.

Fanconi anemia (FA) is a rare genomic instability syndrome. Disease-causing are biallelic mutations in any one of at least 15 genes encoding members of the FA/BRCA pathway of DNA-interstrand crosslink repair. Patients are diagnosed based upon phenotypical manifestations and the diagnosis of FA is confirmed by the hypersensitivity of cells to DNA interstrand crosslinking agents. Customary molecular diagnostics has become increasingly cumbersome, time-consuming and expensive the more FA genes have been identified. We performed Whole Exome Sequencing (WES) in four FA patients in order to investigate the potential of this method for FA genotyping. In search of an optimal WES methodology we explored different enrichment and sequencing techniques. In each case we were able to identify the pathogenic mutations so that WES provided both, complementation group assignment and mutation detection in a single approach. The mutations included homozygous and heterozygous single base pair substitutions and a two-base-pair duplication in FANCJ, -D1, or -D2. Different WES strategies had no critical influence on the individual outcome. However, database errors and in particular pseudogenes impose obstacles that may prevent correct data perception and interpretation, and thus cause pitfalls. With these difficulties in mind, our results show that WES is a valuable tool for the molecular diagnosis of FA and a sufficiently safe technique, capable of engaging increasingly in competition with classical genetic approaches.

Nucleosome positioning patterns derived from human apoptotic nucleosomes.

This communication reports on the nucleosome positioning patterns (bendability matrices) for the human genome, derived from over 8_million nucleosome DNA sequences obtained from apoptotically digested lymphocytes. This digestion procedure is used here for the first time for the purpose of extraction and sequencing of the nucleosome DNA fragments. The dominant motifs suggested by the matrices of DNA bendability calculated for light and heavy isochores are significantly different. Both, however, are in full agreement with the linear description YRRRRRYYYYYR, and with earlier derivations by N-gram extensions. Thus, the choice of the nucleosome positioning patterns crucially depends on the G + C composition of the analyzed sequences.

Human nucleosomes: special role of CG dinucleotides and Alu-nucleosomes.

BACKGROUND: The periodical occurrence of dinucleotides with a period of 10.4 bases now is undeniably a hallmark of nucleosome positioning. Whereas many eukaryotic genomes contain visible and even strong signals for periodic distribution of dinucleotides, the human genome is rather featureless in this respect. The exact sequence features in the human genome that govern the nucleosome positioning remain largely unknown. RESULTS: When analyzing the human genome sequence with the positional autocorrelation method, we found that only the dinucleotide CG shows the 10.4 base periodicity, which is indicative of the presence of nucleosomes. There is a high occurrence of CG dinucleotides that are either 31 (10.4 × 3) or 62 (10.4 × 6) base pairs apart from one another - a sequence bias known to be characteristic of Alu-sequences. In a similar analysis with repetitive sequences removed, peaks of repeating CG motifs can be seen at positions 10, 21 and 31, the nearest integers of multiples of 10.4. CONCLUSIONS: Although the CG dinucleotides are dominant, other elements of the standard nucleosome positioning pattern are present in the human genome as well.The positional autocorrelation analysis of the human genome demonstrates that the CG dinucleotide is, indeed, one visible element of the human nucleosome positioning pattern, which appears both in Alu sequences and in sequences without repeats. The dominant role that CG dinucleotides play in organizing human chromatin is to indicate the involvement of human nucleosomes in tuning the regulation of gene expression and chromatin structure, which is very likely due to cytosine-methylation/-demethylation in CG dinucleotides contained in the human nucleosomes. This is further confirmed by the positions of CG-periodical nucleosomes on Alu sequences. Alu repeats appear as monomers, dimers and trimers, harboring two to six nucleosomes in a run. Considering the exceptional role CG dinucleotides play in the nucleosome positioning, we hypothesize that Alu-nucleosomes, especially, those that form tightly positioned runs, could serve as "anchors" in organizing the chromatin in human cells.

The neuronal transporter gene SLC6A15 confers risk to major depression.

Major depression (MD) is one of the most prevalent psychiatric disorders and a leading cause of loss in work productivity. A combination of genetic and environmental risk factors probably contributes to MD. We present data from a genome-wide association study revealing a neuron-specific neutral amino acid transporter (SLC6A15) as a susceptibility gene for MD. Risk allele carrier status in humans and chronic stress in mice were associated with a downregulation of the expression of this gene in the hippocampus, a brain region implicated in the pathophysiology of MD. The same polymorphisms also showed associations with alterations in hippocampal volume and neuronal integrity. Thus, decreased SLC6A15 expression, due to genetic or environmental factors, might alter neuronal circuits related to the susceptibility for MD. Our convergent data from human genetics, expression studies, brain imaging, and animal models suggest a pathophysiological mechanism for MD that may be accessible to drug targeting.

Nucleosome DNA sequence structure of isochores.

BACKGROUND: Significant differences in G+C content between different isochore types suggest that the nucleosome positioning patterns in DNA of the isochores should be different as well. RESULTS: Extraction of the patterns from the isochore DNA sequences by Shannon N-gram extension reveals that while the general motif YRRRRRYYYYYR is characteristic for all isochore types, the dominant positioning patterns of the isochores vary between TAAAAATTTTTA and CGGGGGCCCCCG due to the large differences in G+C composition. This is observed in human, mouse and chicken isochores, demonstrating that the variations of the positioning patterns are largely G+C dependent rather than species-specific. The species-specificity of nucleosome positioning patterns is revealed by dinucleotide periodicity analyses in isochore sequences. While human sequences are showing CG periodicity, chicken isochores display AG (CT) periodicity. Mouse isochores show very weak CG periodicity only. CONCLUSIONS: Nucleosome positioning pattern as revealed by Shannon N-gram extension is strongly dependent on G+C content and different in different isochores. Species-specificity of the pattern is subtle. It is reflected in the choice of preferentially periodical dinucleotides.

SLX4, a coordinator of structure-specific endonucleases, is mutated in a new Fanconi anemia subtype.

DNA interstrand crosslink repair requires several classes of proteins, including structure-specific endonucleases and Fanconi anemia proteins. SLX4, which coordinates three separate endonucleases, was recently recognized as an important regulator of DNA repair. Here we report the first human individuals found to have biallelic mutations in SLX4. These individuals, who were previously diagnosed as having Fanconi anemia, add SLX4 as an essential component to the FA-BRCA genome maintenance pathway.

Genome-wide association-, replication-, and neuroimaging study implicates HOMER1 in the etiology of major depression.

BACKGROUND: Genome-wide association studies are a powerful tool for unravelling the genetic background of complex disorders such as major depression. METHODS: We conducted a genome-wide association study of 604 patients with major depression and 1364 population based control subjects. The top hundred findings were followed up in a replication sample of 409 patients and 541 control subjects. RESULTS: Two SNPs showed nominally significant association in both the genome-wide association study and the replication samples: 1) rs9943849 (p(combined) = 3.24E-6) located upstream of the carboxypeptidase M (CPM) gene and 2) rs7713917 (p(combined) = 1.48E-6), located in a putative regulatory region of HOMER1. Further evidence for HOMER1 was obtained through gene-wide analysis while conditioning on the genotypes of rs7713917 (p(combined) = 4.12E-3). Homer1 knockout mice display behavioral traits that are paradigmatic of depression, and transcriptional variants of Homer1 result in the dysregulation of cortical-limbic circuitry. This is consistent with the findings of our subsequent human imaging genetics study, which revealed that variation in single nucleotide polymorphism rs7713917 had a significant influence on prefrontal activity during executive cognition and anticipation of reward. CONCLUSION: Our findings, combined with evidence from preclinical and animal studies, suggest that HOMER1 plays a role in the etiology of major depression and that the genetic variation affects depression via the dysregulation of cognitive and motivational processes.