Genetic risk for schizophrenia is associated with increased proportion of indirect connections in brain networks revealed by a semi-metric analysis: evidence from population sample stratified for polygenic risk

Research studies based on tractography have revealed a prominent reduction of asymmetry in some key white-matter tracts in schizophrenia (SCZ). However, we know little about the influence of common genetic risk factors for SCZ on the efficiency of routing on structural brain networks (SBNs). Here, we use a novel recall-by-genotype approach, where we sample young adults from a population-based cohort (ALSPAC:N genotyped = 8,365) based on their burden of common SCZ risk alleles as defined by polygenic risk score (PRS). We compared 181 individuals at extremes of low (N = 91) or high (N = 90) SCZ-PRS under a robust diffusion MRI-based graph theoretical SBN framework. We applied a semi-metric analysis revealing higher SMR values for the high SCZ-PRS group compared with the low SCZ-PRS group in the left hemisphere. Furthermore, a hemispheric asymmetry index showed a higher leftward preponderance of indirect connections for the high SCZ-PRS group compared with the low SCZ-PRS group (PFDR < 0.05). These findings might indicate less efficient structural connectivity in the higher genetic risk group. This is the first study in a population-based sample that reveals differences in the efficiency of SBNs associated with common genetic risk variants for SCZ.

Identifying the Common Genetic Basis of Antidepressant Response

Background

Antidepressants are a first-line treatment for depression. However, only a third of individuals experience remission after the first treatment. Common genetic variation, in part, likely regulates antidepressant response, yet the success of previous genome-wide association studies has been limited by sample size. This study performs the largest genetic analysis of prospectively assessed antidepressant response in major depressive disorder to gain insight into the underlying biology and enable out-of-sample prediction.

Methods

Genome-wide analysis of remission (n remit = 1852, n nonremit = 3299) and percentage improvement (n = 5218) was performed. Single nucleotide polymorphism-based heritability was estimated using genome-wide complex trait analysis. Genetic covariance with eight mental health phenotypes was estimated using polygenic scores/AVENGEME. Out-of-sample prediction of antidepressant response polygenic scores was assessed. Gene-level association analysis was performed using MAGMA and transcriptome-wide association study. Tissue, pathway, and drug binding enrichment were estimated using MAGMA.

Results

Neither genome-wide association study identified genome-wide significant associations. Single nucleotide polymorphism-based heritability was significantly different from zero for remission (h 2 = 0.132, SE = 0.056) but not for percentage improvement (h 2 = -0.018, SE = 0.032). Better antidepressant response was negatively associated with genetic risk for schizophrenia and positively associated with genetic propensity for educational attainment. Leave-one-out validation of antidepressant response polygenic scores demonstrated significant evidence of out-of-sample prediction, though results varied in external cohorts. Gene-based analyses identified ETV4 and DHX8 as significantly associated with antidepressant response.

Conclusions

This study demonstrates that antidepressant response is influenced by common genetic variation, has a genetic overlap schizophrenia and educational attainment, and provides a useful resource for future research. Larger sample sizes are required to attain the potential of genetics for understanding and predicting antidepressant response.

Global Brain Flexibility During Working Memory Is Reduced in a High-Genetic-Risk Group for Schizophrenia

BACKGROUND

Altered functional brain connectivity has been proposed as an intermediate phenotype between genetic risk loci and clinical expression of schizophrenia. Genetic high-risk groups of healthy subjects are particularly suited for the investigation of this proposition because they can be tested in the absence of medication or other secondary effects of schizophrenia.

METHODS

Here, we applied dynamic functional connectivity analysis to functional magnetic resonance imaging data to reveal the reconfiguration of brain networks during a cognitive task. We recruited healthy carriers of common risk variants using the recall-by-genotype design. We assessed 197 individuals: 99 individuals (52 female, 47 male) with low polygenic risk scores (schizophrenia risk profile scores [SCZ-PRSs]) and 98 individuals (52 female, 46 male) with high SCZ-PRSs from both tails of the SCZ-PRS distribution from a genotyped population cohort, the Avon Longitudinal Study of Parents and Children (N = 8169). We compared groups both on conventional brain activation profiles, using the general linear model of the experiment, and on the neural flexibility index, which quantifies how frequent a brain region's community affiliation changes over experimental time.

RESULTS

Behavioral performance and standard brain activation profiles did not differ significantly between groups. High SCZ-PRS was associated with reduced flexibility index and network modularity across n-back levels. The whole-brain flexibility index and that of the frontoparietal working memory network was associated with n-back performance. We identified a dynamic network phenotype related to high SCZ-PRS.

CONCLUSIONS

Such neurophysiological markers can become important for the elucidation of biological mechanisms of schizophrenia and, particularly, the associated cognitive deficit.

The genetic architecture of the human cerebral cortex

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder.

The Relationship Between Polygenic Risk Scores and Cognition in Schizophrenia

BACKGROUND

Cognitive impairment is a clinically important feature of schizophrenia. Polygenic risk score (PRS) methods have demonstrated genetic overlap between schizophrenia, bipolar disorder (BD), major depressive disorder (MDD), educational attainment (EA), and IQ, but very few studies have examined associations between these PRS and cognitive phenotypes within schizophrenia cases.

METHODS

We combined genetic and cognitive data in 3034 schizophrenia cases from 11 samples using the general intelligence factor g as the primary measure of cognition. We used linear regression to examine the association between cognition and PRS for EA, IQ, schizophrenia, BD, and MDD. The results were then meta-analyzed across all samples. A genome-wide association studies (GWAS) of cognition was conducted in schizophrenia cases.

RESULTS

PRS for both population IQ (P = 4.39 × 10-28) and EA (P = 1.27 × 10-26) were positively correlated with cognition in those with schizophrenia. In contrast, there was no association between cognition in schizophrenia cases and PRS for schizophrenia (P = .39), BD (P = .51), or MDD (P = .49). No individual variant approached genome-wide significance in the GWAS.

CONCLUSIONS

Cognition in schizophrenia cases is more strongly associated with PRS that index cognitive traits in the general population than PRS for neuropsychiatric disorders. This suggests the mechanisms of cognitive variation within schizophrenia are at least partly independent from those that predispose to schizophrenia diagnosis itself. Our findings indicate that this cognitive variation arises at least in part due to genetic factors shared with cognitive performance in populations and is not solely due to illness or treatment-related factors, although our findings are consistent with important contributions from these factors.

Convergent Evidence That ZNF804A Is a Regulator of Pre-messenger RNA Processing and Gene Expression

Genome-wide association studies have linked common variation in ZNF804A with an increased risk of schizophrenia. However, little is known about the biology of ZNF804A and its role in schizophrenia. Here, we investigate the function of ZNF804A using a variety of complementary molecular techniques. We show that ZNF804A is a nuclear protein that interacts with neuronal RNA splicing factors and RNA-binding proteins including RBFOX1, which is also associated with schizophrenia, CELF3/4, components of the ubiquitin-proteasome system and the ZNF804A paralog, GPATCH8. GPATCH8 also interacts with splicing factors and is localized to nuclear speckles indicative of a role in pre-messenger RNA (mRNA) processing. Sequence analysis showed that GPATCH8 contains ultraconserved, alternatively spliced poison exons that are also regulated by RBFOX proteins. ZNF804A knockdown in SH-SY5Y cells resulted in robust changes in gene expression and pre-mRNA splicing converging on pathways associated with nervous system development, synaptic contact, and cell adhesion. We observed enrichment (P = 1.66 × 10-9) for differentially spliced genes in ZNF804A-depleted cells among genes that contain RBFOX-dependent alternatively spliced exons. Differentially spliced genes in ZNF804A-depleted cells were also enriched for genes harboring de novo loss of function mutations in autism spectrum disorder (P = 6.25 × 10-7, enrichment 2.16) and common variant alleles associated with schizophrenia (P = .014), bipolar disorder and schizophrenia (P = .003), and autism spectrum disorder (P = .005). These data suggest that ZNF804A and its paralogs may interact with neuronal-splicing factors and RNA-binding proteins to regulate the expression of a subset of synaptic and neurodevelopmental genes.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Oscillatory hyperactivity and hyperconnectivity in young APOE-ɛ4 carriers and hypoconnectivity in Alzheimer's disease

We studied resting-state oscillatory connectivity using magnetoencephalography in healthy young humans (N = 183) genotyped for APOE-ɛ4, the greatest genetic risk for Alzheimer's disease (AD). Connectivity across frequencies, but most prevalent in alpha/beta, was increased in APOE-ɛ4 in a set of mostly right-hemisphere connections, including lateral parietal and precuneus regions of the Default Mode Network. Similar regions also demonstrated hyperactivity, but only in gamma (40-160 Hz). In a separate study of AD patients, hypoconnectivity was seen in an extended bilateral network that partially overlapped with the hyperconnected regions seen in young APOE-ɛ4 carriers. Using machine-learning, AD patients could be distinguished from elderly controls with reasonable sensitivity and specificity, while young APOE-e4 carriers could also be distinguished from their controls with above chance performance. These results support theories of initial hyperconnectivity driving eventual profound disconnection in AD and suggest that this is present decades before the onset of AD symptomology.

Structural and Functional Neuroimaging of Polygenic Risk for Schizophrenia: A Recall-by-Genotype-Based Approach

Risk profile scores (RPS) derived from genome-wide association studies (GWAS) explain a considerable amount of susceptibility for schizophrenia (SCZ). However, little is known about how common genetic risk factors for SCZ influence the structure and function of the human brain, largely due to the constraints of imaging sample sizes. In the current study, we use a novel recall-by-genotype (RbG) methodological approach, where we sample young adults from a population cohort (Avon Longitudinal Study of Parents and Children: N genotyped = 8365) based on their SCZ-RPS. We compared 197 healthy individuals at extremes of low (N = 99) or high (N = 98) SCZ-RPS with behavioral tests, and structural and functional magnetic resonance imaging (fMRI). We first provide methodological details that will inform the design of future RbG studies for common SCZ genetic risk. We further provide an between group analysis of the RbG individuals (low vs high SCZ-RPS) who underwent structural neuroimaging data (T1-weighted scans) and fMRI data during a reversal learning task. While we found little evidence for morphometric differences between the low and high SCZ-RPS groups, we observed an impact of SCZ-RPS on blood oxygen level-dependent (BOLD) signal during reward processing in the ventral striatum (PFWE-VS-CORRECTED = .037), a previously investigated broader reward-related network (PFWE-ROIS-CORRECTED = .008), and across the whole brain (PFWE-WHOLE-BRAIN-CORRECTED = .013). We also describe the study strategy and discuss specific challenges of RbG for SCZ risk (such as SCZ-RPS related homoscedasticity). This study will help to elucidate the behavioral and imaging phenotypes that are associated with SCZ genetic risk.

Loss of offspring Peg3 reduces neonatal ultrasonic vocalizations and increases maternal anxiety in wild-type mothers

Depression and anxiety are the most common mental health conditions during pregnancy and can impair the normal development of mother-infant interactions. These adversities are associated with low birth weight and increased risk of behavioural disorders in children. We recently reported reduced expression of the imprinted gene PATERNALLY EXPRESSED GENE 3 (PEG3) in placenta of human infants born to depressed mothers. Expression of Peg3 in the brain has previously been linked maternal behaviour in rodents, at least in some studies, with mutant dams neglecting their pups. However, in our human study decreased expression was in the placenta derived from the fetus. Here, we examined maternal behaviour in response to reduced expression of Peg3 in the feto-placental unit. Prenatally we found novelty reactivity was altered in wild-type females carrying litters with a null mutation in Peg3. This behavioural alteration was short-lived and there were no significant differences the transcriptomes of either the maternal hypothalamus or hippocampus at E16.5. In contrast, while maternal gross maternal care was intact postnatally, the exposed dams were significantly slower to retrieve their pups and displayed a marked increase in anxiety. We also observed a significant reduction in the isolation-induced ultrasonic vocalizations (USVs) emitted by mutant pups separated from their mothers. USVs are a form of communication known to elicit maternal care suggesting Peg3 mutant pups drive the deficit in maternal behaviour. These data support the hypothesis that reduced placental PEG3 in human pregnancies occurs as a consequence of prenatal depression but leaves scope for feto-placental Peg3 dosage, during gestation, influencing aspects of maternal behaviour.

Fractional anisotropy of the uncinate fasciculus and cingulum in bipolar disorder type I, type II, unaffected siblings and healthy controls

BACKGROUND

Fractional anisotropy in the uncinate fasciculus and the cingulum may be biomarkers for bipolar disorder and may even be distinctly affected in different subtypes of bipolar disorder, an area in need of further research.AimsThis study aims to establish if fractional anisotropy in the uncinate fasciculus and cingulum shows differences between healthy controls, patients with bipolar disorder type I (BD-I) and type II (BD-II), and their unaffected siblings.

METHOD

Fractional anisotropy measures from the uncinate fasciculus, cingulum body and parahippocampal cingulum were compared with tractography methods in 40 healthy controls, 32 patients with BD-I, 34 patients with BD-II, 17 siblings of patients with BD-I and 14 siblings of patients with BD-II.

RESULTS

The main effects were found in both the right and left uncinate fasciculus, with patients with BD-I showing significantly lower fractional anisotropy than both patients with BD-II and healthy controls. Participants with BD-II did not differ from healthy controls. Siblings showed similar effects in the left uncinate fasciculus. In a subsequent complementary analysis, we investigated the association between fractional anisotropy in the uncinate fasciculus and polygenic risk for bipolar disorder and psychosis in a large cohort (n = 570) of healthy participants. However, we found no significant association.

CONCLUSIONS

Fractional anisotropy in the uncinate fasciculus differs significantly between patients with BD-I and patients with BD-II and healthy controls. This supports the hypothesis of differences in the physiological sub-tract between bipolar disorder subtypes. Similar results were found in unaffected siblings, suggesting the potential for this biomarker to represent an endophenotype for BD-I. However, fractional anisotropy in the uncinate fasciculus seems unrelated to polygenic risk for bipolar disorder or psychosis.Declaration of interestNone.

Effect of cytochrome CYP2C19 metabolizing activity on antidepressant response and side effects: Meta-analysis of data from genome-wide association studies

Cytochrome (CYP) P450 enzymes have a primary role in antidepressant metabolism and variants in these polymorphic genes are targets for pharmacogenetic investigation. This is the first meta-analysis to investigate how CYP2C19 polymorphisms predict citalopram/escitalopram efficacy and side effects. CYP2C19 metabolic phenotypes comprise poor metabolizers (PM), intermediate and intermediate+ metabolizers (IM; IM+), extensive and extensive+ metabolizers (EM [wild type]; EM+) and ultra-rapid metabolizers (UM) defined by the two most common CYP2C19 functional polymorphisms (rs4244285 and rs12248560) in Caucasians. These polymorphisms were genotyped or imputed from genome-wide data in four samples treated with citalopram or escitalopram (GENDEP, STAR*D, GenPod, PGRN-AMPS). Treatment efficacy was assessed by standardized percentage symptom improvement and by remission. Side effect data were available at weeks 2-4, 6 and 9 in three samples. A fixed-effects meta-analysis was performed using EM as the reference group. Analysis of 2558 patients for efficacy and 2037 patients for side effects showed that PMs had higher symptom improvement (SMD = 0.43, CI = 0.19-0.66) and higher remission rates (OR = 1.55, CI = 1.23-1.96) compared to EMs. At weeks 2-4, PMs showed higher risk of gastro-intestinal (OR = 1.26, CI = 1.08-1.47), neurological (OR = 1.28, CI = 1.07-1.53) and sexual side effects (OR = 1.52, CI = 1.23-1.87; week 6 values were similar). No difference was seen at week 9 or in total side effect burden. PMs did not have higher risk of dropout at week 4 compared to EMs. Antidepressant dose was not different among CYP2C19 groups. CYP2C19 polymorphisms may provide helpful information for guiding citalopram/escitalopram treatment, despite PMs being relatively rare among Caucasians (∼2%).

The Psychiatric Risk Gene Transcription Factor 4 (TCF4) Regulates Neurodevelopmental Pathways Associated With Schizophrenia, Autism, and Intellectual Disability

BACKGROUND

Common genetic variants in and around the gene encoding transcription factor 4 (TCF4) are associated with an increased risk of schizophrenia. Conversely, rare damaging TCF4 mutations cause Pitt-Hopkins syndrome and have also been found in individuals with intellectual disability (ID) and autism spectrum disorder (ASD).

METHODS

Chromatin immunoprecipitation and next generation sequencing were used to identify the genomic targets of TCF4. These data were integrated with expression, epigenetic and disease gene sets using a range of computational tools.

RESULTS

We identify 10604 TCF4 binding sites in the genome that were assigned to 5437 genes. De novo motif enrichment found that most TCF4 binding sites contained at least one E-box (5'-CAtcTG). Approximately 77% of TCF4 binding sites overlapped with the H3K27ac histone modification for active enhancers. Enrichment analysis on the set of TCF4 targets identified numerous, highly significant functional clusters for pathways including nervous system development, ion transport and signal transduction, and co-expression modules for genes associated with synaptic function and brain development. Importantly, we found that genes harboring de novo mutations in schizophrenia (P = 5.3 × 10-7), ASD (P = 2.5 × 10-4), and ID (P = 7.6 × 10-3) were also enriched among TCF4 targets. TCF4 binding sites were also found at other schizophrenia risk loci including the nicotinic acetylcholine receptor cluster, CHRNA5/CHRNA3/CHRNB4 and SETD1A.

CONCLUSIONS

These data demonstrate that TCF4 binding sites are found in a large number of neuronal genes that include many genetic risk factors for common neurodevelopmental disorders.

Examining cognition across the bipolar/schizophrenia diagnostic spectrum

BACKGROUND

Cognitive impairments are well-established features of schizophrenia, but there is ongoing debate about the nature and degree of cognitive impairment in patients with schizoaffective disorder and bipolar disorder. We hypothesized that there is a spectrum of increasing impairment from bipolar disorder to schizoaffective disorder bipolar type, to schizoaffective disorder depressive type and schizophrenia.

METHODS

We compared performance on the Measurement and Treatment Research to Improve Cognition in Schizophrenia (MATRICS) Consensus Cognitive Battery between participants with schizophrenia (n = 558), schizoaffective disorder depressive type (n = 112), schizoaffective disorder type (n = 76), bipolar disorder (n = 78) and healthy participants (n = 103) using analysis of covariance with post hoc comparisons. We conducted an ordinal logistic regression to examine whether cognitive impairments followed the hypothesized spectrum from bipolar disorder (least severe) to schizophrenia (most severe). In addition to categorical diagnoses, we addressed the influence of symptom domains, examining the association between cognition and mania, depression and psychosis.

RESULTS

Cognitive impairments increased in severity from bipolar disorder to schizoaffective disorder bipolar type, to schizophrenia and schizoaffective disorder depressive type. Participants with schizophrenia and schizoaffective disorder depressive type showed equivalent performance (d = 0.07, p = 0.90). The results of the ordinal logistic regression were consistent with a spectrum of deficits from bipolar disorder to schizoaffective disorder bipolar type, to schizophrenia/schizoaffective disorder depressive type (odds ratio = 1.98, p < 0.001). In analyses of the associations between symptom dimensions and cognition, higher scores on the psychosis dimension were associated with poorer performance (B = 0.015, standard error = 0.002, p < 0.001).

LIMITATIONS

There were fewer participants with schizoaffective disorder and bipolar disorder than schizophrenia. Despite this, our analyses were robust to differences in group sizes, and we were able to detect differences between groups.

CONCLUSION

Cognitive impairments represent a symptom dimension that cuts across traditional diagnostic boundaries.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the version of this article originally published, one of the two authors with the name Wei Zhao was omitted from the author list and the affiliations for both authors were assigned to the single Wei Zhao in the author list. In addition, the ORCID for Wei Zhao (Department of Biostatistics and Epidemiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA) was incorrectly assigned to author Wei Zhou. The errors have been corrected in the HTML and PDF versions of the article.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the published version of this paper, the name of author Emanuele Di Angelantonio was misspelled. This error has now been corrected in the HTML and PDF versions of the article.

Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression

Major depressive disorder (MDD) is a common illness accompanied by considerable morbidity, mortality, costs, and heightened risk of suicide. We conducted a genome-wide association meta-analysis based in 135,458 cases and 344,901 controls and identified 44 independent and significant loci. The genetic findings were associated with clinical features of major depression and implicated brain regions exhibiting anatomical differences in cases. Targets of antidepressant medications and genes involved in gene splicing were enriched for smaller association signal. We found important relationships of genetic risk for major depression with educational attainment, body mass, and schizophrenia: lower educational attainment and higher body mass were putatively causal, whereas major depression and schizophrenia reflected a partly shared biological etiology. All humans carry lesser or greater numbers of genetic risk factors for major depression. These findings help refine the basis of major depression and imply that a continuous measure of risk underlies the clinical phenotype.

Examining cognition across the bipolar/schizophrenia diagnostic spectrum

BACKGROUND

Cognitive impairments are well-established features of schizophrenia, but there is ongoing debate about the nature and degree of cognitive impairment in patients with schizoaffective disorder and bipolar disorder. We hypothesized that there is a spectrum of increasing impairment from bipolar disorder to schizoaffective disorder bipolar type, to schizoaffective disorder depressive type and schizophrenia.

METHODS

We compared performance on the Measurement and Treatment Research to Improve Cognition in Schizophrenia (MATRICS) Consensus Cognitive Battery between participants with schizophrenia (n = 558), schizoaffective disorder depressive type (n = 112), schizoaffective disorder type (n = 76), bipolar disorder (n = 78) and healthy participants (n = 103) using analysis of covariance with post hoc comparisons. We conducted an ordinal logistic regression to examine whether cognitive impairments followed the hypothesized spectrum from bipolar disorder (least severe) to schizophrenia (most severe). In addition to categorical diagnoses, we addressed the influence of symptom domains, examining the association between cognition and mania, depression and psychosis.

RESULTS

Cognitive impairments increased in severity from bipolar disorder to schizoaffective disorder bipolar type, to schizophrenia and schizoaffective disorder depressive type. Participants with schizophrenia and schizoaffective disorder depressive type showed equivalent performance (d = 0.07, p = 0.90). The results of the ordinal logistic regression were consistent with a spectrum of deficits from bipolar disorder to schizoaffective disorder bipolar type, to schizophrenia/schizoaffective disorder depressive type (odds ratio = 1.98, p < 0.001). In analyses of the associations between symptom dimensions and cognition, higher scores on the psychosis dimension were associated with poorer performance (B = 0.015, standard error = 0.002, p < 0.001).

LIMITATIONS

There were fewer participants with schizoaffective disorder and bipolar disorder than schizophrenia. Despite this, our analyses were robust to differences in group sizes, and we were able to detect differences between groups.

CONCLUSION

Cognitive impairments represent a symptom dimension that cuts across traditional diagnostic boundaries.

Genetic risk for Alzheimer's disease is concentrated in specific macrophage and microglial transcriptional networks

Background

Genome-wide association studies of Alzheimer’s disease (AD) have identified a number of significant risk loci, the majority of which lie in non-coding regions of the genome. The lack of causal alleles and considerable polygenicity remains a significant barrier to translation into mechanistic understanding. This includes identifying causal variants and the cell/tissue types in which they operate. A fuller understanding of the cell types and transcriptional networks involved in AD genetic risk mechanisms will provide important insights into pathogenesis.

Methods

We assessed the significance of the overlap between genome-wide significant AD risk variants and sites of open chromatin from data sets representing diverse tissue types. We then focussed on macrophages and microglia to investigate the role of open chromatin sites containing motifs for specific transcription factors. Partitioned heritability using LDscore regression was used to investigate the contribution of specific macrophage and microglia transcription factor motif-containing open chromatin sites to the heritability of AD.

Results

AD risk single nucleotide polymorphisms (SNPs) are preferentially located at sites of open chromatin in immune cells, particularly monocytes (z score = 4.43; corrected P = 5.88 × 10^{− 3}). Similar enrichments are observed for macrophages (z score = 4.10; corrected P < 2.40 × 10^{− 3}) and microglia (z score = 4.34, corrected P = 0.011). In both macrophages and microglia, AD risk variants are enriched at a subset of open chromatin sites that contain DNA binding motifs for specific transcription factors, e.g. SPI1 and MEF2. Genetic variation at many of these motif-containing sites also mediate a substantial proportion of AD heritability, with SPI1-containing sites capturing the majority of the common variant SNP-chip heritability (microglia enrichment = 16.28, corrected enrichment P = 0.0044).

Conclusions

AD risk alleles plausibly operate in immune cells, including microglia, and are concentrated in specific transcriptional networks. Combined with primary genetic association results, the SPI1 and MEF2 transcriptional networks appear central to AD risk mechanisms. Investigation of transcription factors targeting AD risk SNP associated regulatory elements could provide powerful insights into the molecular processes affected by AD polygenic risk. More broadly, our findings support a model of polygenic disease risk that arises from variants located in specific transcriptional networks.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0523-8) contains supplementary material, which is available to authorized users.

Formalising recall by genotype as an efficient approach to detailed phenotyping and causal inference

Detailed phenotyping is required to deepen our understanding of the biological mechanisms behind genetic associations. In addition, the impact of potentially modifiable risk factors on disease requires analytical frameworks that allow causal inference. Here, we discuss the characteristics of Recall-by-Genotype (RbG) as a study design aimed at addressing both these needs. We describe two broad scenarios for the application of RbG: studies using single variants and those using multiple variants. We consider the efficacy and practicality of the RbG approach, provide a catalogue of UK-based resources for such studies and present an online RbG study planner.

Enrichment of schizophrenia heritability in both neuronal and glia cell regulatory elements

Genome-wide association studies have identified over 100 robust risk loci for schizophrenia with thousands of variants mediating genetic heritability, the majority of which reside in non-coding regions. Analytical approaches have shown this heritability is strongly enriched at variants within regulatory elements identified from human post-mortem brain tissue. However, bulk post-mortem brain tissue has a heterogeneous cell composition, making biological interpretations difficult. We sought to refine the cell types mediating schizophrenia heritability by separating neuronal and glial signals using data from: (1) NeuN-sorted post-mortem brain and (2) cell culture systems. Schizophrenia heritability was partitioned using linkage disequilbrium (LD) score regression. Variants within genomic regions marked by H3K4me3 (marker of active promoters) from NeuN-positive (neuronal) and NeuN-negative (non-neuronal) cells explained a significant amount of schizophrenia heritability (P = 1.38 × 10^-10 and P = 7.97 × 10^-10). However, variants located in H3K4me3 sites specific to NeuN-positive (neuronal) cells were enriched (P = 3.13 × 10^-4), while those specific to NeuN-negative (non-neuronal) cells were not (P = 0.470). Data from cell culture systems mimicked this pattern of association. We show the previously observed enrichment of heritability from variants at brain H3K4me3 sites is mediated by both neuronal and non-neuronal brain cell types. However, only neuronal cell populations showed a unique contribution driven by cell-type specific regulatory elements. Cell culture systems recapitulate disease relevant gene-regulatory landscapes, validating them as a tool for future investigation of genetic mechanisms underlying schizophrenia. Identifying the cell types in which risk variants operate will greatly increase our understanding of schizophrenia pathobiology and aid in the development of novel model systems and therapies.

Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

An Analysis of Two Genome-wide Association Meta-analyses Identifies a New Locus for Broad Depression Phenotype

BACKGROUND

The genetics of depression has been explored in genome-wide association studies that focused on either major depressive disorder or depressive symptoms with mostly negative findings. A broad depression phenotype including both phenotypes has not been tested previously using a genome-wide association approach. We aimed to identify genetic polymorphisms significantly associated with a broad phenotype from depressive symptoms to major depressive disorder.

METHODS

We analyzed two prior studies of 70,017 participants of European ancestry from general and clinical populations in the discovery stage. We performed a replication meta-analysis of 28,328 participants. Single nucleotide polymorphism (SNP)-based heritability and genetic correlations were calculated using linkage disequilibrium score regression. Discovery and replication analyses were performed using a p-value-based meta-analysis. Lifetime major depressive disorder and depressive symptom scores were used as the outcome measures.

RESULTS

The SNP-based heritability of major depressive disorder was 0.21 (SE = 0.02), the SNP-based heritability of depressive symptoms was 0.04 (SE = 0.01), and their genetic correlation was 1.001 (SE = 0.2). We found one genome-wide significant locus related to the broad depression phenotype (rs9825823, chromosome 3: 61,082,153, p = 8.2 × 10^-9) located in an intron of the FHIT gene. We replicated this SNP in independent samples (p = .02) and the overall meta-analysis of the discovery and replication cohorts (1.0 × 10^-9).

CONCLUSIONS

This large study identified a new locus for depression. Our results support a continuum between depressive symptoms and major depressive disorder. A phenotypically more inclusive approach may help to achieve the large sample sizes needed to detect susceptibility loci for depression.

The genome-wide expression effects of escitalopram and its relationship to neurogenesis, hippocampal volume, and antidepressant response

Antidepressant-induced hippocampal neurogenesis (AHN) is hypothesized to contribute to increases in hippocampal volume among major depressive disorder patients after long-term treatment. Furthermore, rodent studies suggest AHN may be the cellular mechanism mediating the therapeutic benefits of antidepressants. Here, we perform the first investigation of genome-wide expression changes associated with AHN in human cells. We identify gene expression networks significantly activated during AHN, and we perform gene set analyses to probe the molecular relationship between AHN, hippocampal volume, and antidepressant response. The latter were achieved using genome-wide association summary data collected from 30,717 individuals as part of the ENIGMA Consortium (genetic predictors of hippocampal volume dataset), and data collected from 1,222 major depressed patients as part of the NEWMEDS Project (genetic predictors of response to antidepressants dataset). Our results showed that the selective serotonin reuptake inhibitor, escitalopram evoked AHN in human cells; dose-dependently increasing the differentiation of cells into neuroblasts, as well as increasing gliogenesis. Activated genome-wide expression networks relate to axon and microtubule formation, and ribosomal biogenesis. Gene set analysis revealed that gene expression changes associated with AHN were nominally enriched for genes predictive of hippocampal volume, but not for genes predictive of therapeutic response.

Association between C-reactive protein (CRP) with depression symptom severity and specific depressive symptoms in major depression

INTRODUCTION

Population-based studies have associated inflammation, particularly higher C-reactive protein (CRP), with depressive severity, but clinical trials in major depressive disorder were rather non-specific without examining the role of gender. We aimed to investigate the association between CRP and overall depression severity including specific depressive symptoms and to examine potential gender differences.

METHODS

We included 231 individuals with major depressive disorder from the Genome-Based Therapeutics Drugs for Depression (GENDEP) study. At baseline, we assessed high-sensitivity CRP levels and psychopathology with the Montgomery Aasberg Depression Rating Scale (MADRS). We performed linear regression analyses to investigate the association between baseline CRP levels with overall MADRS severity and specific symptoms at baseline and adjusted for age, gender, anti-inflammatory and psychotropic drug treatment, body mass index, smoking, inflammatory diseases, and recruitment center.

RESULTS

Higher CRP levels were significantly associated with greater overall MADRS symptom severity (p=0.02), which was significant among women (p=0.02) but not among men (p=0.68). Among women, higher CRP was associated with increased severity on observed mood, cognitive symptoms, interest-activity, and suicidality, but we found no significant associations among men. Interaction analyses showed no significant gender differences on the overall MADRS score or specific symptoms.

DISCUSSION

Our results support the sickness syndrome theory suggesting that chronic low-grade inflammation may be associated with a subtype of depression. The potential gender differences in psychopathology may be explained by biological and/or psychosocial factors, e.g. differential modulation of immune responses by sex hormones. Clinical studies should investigate symptom-specific and/or gender-specific treatment guided by peripheral inflammatory markers.

Rare and low-frequency coding variants alter human adult height

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.

Multimodal Brain Imaging Reveals Structural Differences in Alzheimer's Disease Polygenic Risk Carriers: A Study in Healthy Young Adults

BACKGROUND

Recent genome-wide association studies have identified genetic loci that jointly make a considerable contribution to risk of developing Alzheimer's disease (AD). Because neuropathological features of AD can be present several decades before disease onset, we investigated whether effects of polygenic risk are detectable by neuroimaging in young adults. We hypothesized that higher polygenic risk scores (PRSs) for AD would be associated with reduced volume of the hippocampus and other limbic and paralimbic areas. We further hypothesized that AD PRSs would affect the microstructure of fiber tracts connecting the hippocampus with other brain areas.

METHODS

We analyzed the association between AD PRSs and brain imaging parameters using T1-weighted structural (n = 272) and diffusion-weighted scans (n = 197).

RESULTS

We found a significant association between AD PRSs and left hippocampal volume, with higher risk associated with lower left hippocampal volume (p = .001). This effect remained when the APOE gene was excluded (p = .031), suggesting that the relationship between hippocampal volume and AD is the result of multiple genetic factors and not exclusively variability in the APOE gene. The diffusion tensor imaging analysis revealed that fractional anisotropy of the right cingulum was inversely correlated with AD PRSs (p = .009). We thus show that polygenic effects of AD risk variants on brain structure can already be detected in young adults.

CONCLUSIONS

This finding paves the way for further investigation of the effects of AD risk variants and may become useful for efforts to combine genotypic and phenotypic data for risk prediction and to enrich future prevention trials of AD.

Polygenic Risk of Psychosis and Ventral Striatal Activation During Reward Processing in Healthy Adolescents

IMPORTANCE

Psychotic disorders are characterized by attenuated activity in the brain's valuation system in key reward processing areas, such as the ventral striatum (VS), as measured with functional magnetic resonance imaging.

OBJECTIVE

To examine whether common risk variants for psychosis are associated with individual variation in the VS.

DESIGN, SETTING, AND PARTICIPANTS

A cross-sectional study of a large cohort of adolescents from the IMAGEN study (a European multicenter study of reinforcement sensitivity in adolescents) was performed from March 1, 2008, through December 31, 2011. Data analysis was conducted from October 1, 2015, to January 9, 2016. Polygenic risk profile scores (RPSs) for psychosis were generated for 1841 healthy adolescents. Sample size and characteristics varied across regression analyses, depending on mutual information available (N = 1524-1836).

MAIN OUTCOMES AND MEASURES

Reward-related brain function was assessed with blood oxygen level dependency (BOLD) in the VS using the monetary incentive delay (MID) task, distinguishing reward anticipation and receipt. Behavioral impulsivity, IQ, MID task performance, and VS BOLD were regressed against psychosis RPS at 4 progressive P thresholds (P < .01, P < .05, P < .10, and P < .50 for RPS models 1-4, respectively).

RESULTS

In a sample of 1841 healthy adolescents (mean age, 14.5 years; 906 boys and 935 girls), we replicated an association between increasing psychosis RPS and reduced IQ (matrix reasoning: corrected P = .003 for RPS model 2, 0.4% variance explained), supporting the validity of the psychosis RPS models. We also found a nominally significant association between increased psychosis RPS and reduced MID task performance (uncorrected P = .03 for RPS model 4, 0.2% variance explained). Our main finding was a positive association between psychosis RPS and VS BOLD during reward anticipation at all 4 psychosis RPS models and for 2 P thresholds for reward receipt (RPS models 1 and 3), correcting for the familywise error rate (0.8%-1.9% variance explained).

CONCLUSIONS AND RELEVANCE

These findings support an association between psychosis RPS and VS BOLD in adolescents. Genetic risk for psychosis may shape an individual's response to rewarding stimuli.

Evidence of Common Genetic Overlap Between Schizophrenia and Cognition

Cognitive impairment is a core feature of schizophrenia but there is limited understanding of the genetic relationship between cognition in the general population and schizophrenia. We examine how common variants associated with schizophreniaen massecontribute to childhood cognitive ability in a population-based sample, and the extent to which common genetic variants associated with childhood cognition explain variation in schizophrenia. Schizophrenia polygenic risk scores were derived from the Psychiatric Genomics Consortium (n= 69 516) and tested for association with IQ, attention, processing speed, working memory, problem solving, and social cognition in over 5000 children aged 8 from the Avon Longitudinal Study of Parents and Children birth cohort. Polygenic scores for these cognitive domains were tested for association with schizophrenia in a large UK schizophrenia sample (n= 11 853). Bivariate genome-wide complex trait analysis (GCTA) estimated the amount of shared genetic factors between schizophrenia and cognitive domains. Schizophrenia polygenic risk score was associated with lower performance IQ (P= .001) and lower full IQ (P= .013). Polygenic score for performance IQ was associated with increased risk for schizophrenia (P= 3.56E-04). Bivariate GCTA revealed moderate genetic correlation between schizophrenia and both performance IQ (rG= -.379,P= 6.62E-05) and full IQ (rG= -.202,P= 5.00E-03), with approximately 14% of the genetic component of schizophrenia shared with that for performance IQ. Our results support the presence of shared common genetic factors between schizophrenia and childhood cognitive ability. We observe a genetic relationship between schizophrenia and performance IQ but not verbal IQ or other cognitive variables, which may have implications for studies utilizing cognitive endophenotypes for psychosis.

CACNA1C risk variant is associated with increased amygdala volume

Genome-wide association studies suggest that genetic variation within L-type calcium channel subunits confer risk to psychosis. The single nucleotide polymorphism at rs1006737 in CACNA1C has been associated with both schizophrenia and bipolar disorder and with several intermediate phenotypes that may serve as neurobiological antecedents, linking psychosis to genetic aetiology. Amongst others, it has been implicated in alterations in amygdala structure and function. In the present study, we show that the risk allele (A) is associated with increased amygdala volume in healthy individuals (n = 258). This observation reinforces a hypothesis that genetic variation may confer risk to psychosis via alterations in limbic structures. Further study of CACNA1C using intermediate phenotypes for psychosis will determine the mechanisms by which variation in this gene confers risk.

Phenotypic Manifestation of Genetic Risk for Schizophrenia During Adolescence in the General Population

IMPORTANCE

Schizophrenia is a highly heritable, polygenic condition characterized by a relatively diverse phenotype and frequent comorbid conditions, such as anxiety and depression. At present, limited evidence explains how genetic risk for schizophrenia is manifest in the general population.

OBJECTIVE

To investigate the extent to which genetic risk for schizophrenia is associated with different phenotypes during adolescence in a population-based birth cohort.

DESIGN, SETTING, AND PARTICIPANTS

This cohort study used data from the Avon Longitudinal Study of Parents and Children (ALSPAC). Of 14,062 children in the birth cohort, genetic data were available for 9912 adolescents. Data were collected periodically from September 6, 1990, and collection is ongoing. Data were analyzed from March 4 to August 13, 2015.

EXPOSURES

Polygenic risk scores (PRSs) for schizophrenia generated for individuals in the ALSPAC cohort using results of the second Psychiatric Genomics Consortium Schizophrenia genome-wide association study as a training set.

MAIN OUTCOMES AND MEASURES

Logistic regression was used to assess associations between the schizophrenia PRS and (1) psychotic experiences (Psychosis-Like Symptom Interview at 12 and 18 years of age), (2) negative symptoms (Community Assessment of Psychic Experiences at 16.5 years of age), (3) depressive disorder (Development and Well-Being Assessment at 15.5 years of age), and (4) anxiety disorder (Development and Well-Being Assessment at 15.5 years of age) in adolescence.

RESULTS

Of the 8230 ALSPAC participants whose genetic data passed quality control checks (51.2% male, 48.8% female), 3676 to 5444 participated in assessments from 12 to 18 years of age. The PRSs created using single-nucleotide polymorphisms with a training-set P ≤ .05 threshold were associated with negative symptoms (odds ratio [OR] per SD increase in PRS, 1.21; 95% CI, 1.08-1.36; R(2) = 0.007) and anxiety disorder (OR per SD increase in PRS, 1.17; 95% CI, 1.06- 1.29; R(2) = 0.005). No evidence was found of an association between schizophrenia PRS and psychotic experiences (OR per SD increase in PRS, 1.08; 95% CI, 0.98-1.19; R(2) = 0.001) or depressive disorder (OR per SD increase in PRS, 1.02; 95% CI, 0.91-1.13; R(2) = 0.00005). Results were mostly consistent across different training-set P value thresholds and using different cutoffs and measures of the psychopathological outcomes.

CONCLUSIONS AND RELEVANCE

This study demonstrates polygenic overlaps between common genetic polymorphisms associated with schizophrenia and negative symptoms and anxiety disorder but not with psychotic experiences or depression. Because the genetic risk for schizophrenia appears to be manifest as anxiety and negative symptoms during adolescence, a greater focus on these phenotypes rather than on psychotic experiences might be required for prediction of transition in at-risk samples.

Phenotypic Association Analyses With Copy Number Variation in Recurrent Depressive Disorder

BACKGROUND

Defining the molecular genomic basis of the likelihood of developing depressive disorder is a considerable challenge. We previously associated rare, exonic deletion copy number variants (CNV) with recurrent depressive disorder (RDD). Sex chromosome abnormalities also have been observed to co-occur with RDD.

METHODS

In this reanalysis of our RDD dataset (N = 3106 cases; 459 screened control samples and 2699 population control samples), we further investigated the role of larger CNVs and chromosomal abnormalities in RDD and performed association analyses with clinical data derived from this dataset.

RESULTS

We found an enrichment of Turner's syndrome among cases of depression compared with the frequency observed in a large population sample (N = 34,910) of live-born infants collected in Denmark (two-sided p = .023, odds ratio = 7.76 [95% confidence interval = 1.79-33.6]), a case of diploid/triploid mosaicism, and several cases of uniparental isodisomy. In contrast to our previous analysis, large deletion CNVs were no more frequent in cases than control samples, although deletion CNVs in cases contained more genes than control samples (two-sided p = .0002).

CONCLUSIONS

After statistical correction for multiple comparisons, our data do not support a substantial role for CNVs in RDD, although (as has been observed in similar samples) occasional cases may harbor large variants with etiological significance. Genetic pleiotropy and sample heterogeneity suggest that very large sample sizes are required to study conclusively the role of genetic variation in mood disorders.

Transcriptomics and the mechanisms of antidepressant efficacy

The mechanisms by which antidepressants have their effects are not clear and the reasons for variability in treatment outcomes are also unknown. However, there is evidence from candidate gene research that indicates gene expression changes may be involved in antidepressant action. In this study, we examined antidepressant-induced alterations in gene expression on a transcriptome-wide scale, exploring associations with treatment response. Blood samples were taken from a subset of depressed patients from the GENDEP study (n=136) before and after eight weeks of treatment with either escitalopram or nortriptyline. Transcriptomic data were obtained from these samples using Illumina HumanHT-12 v4 Expression BeadChip microarrays. When analysing individual genes, we observed that changes in the expression of two genes (MMP28 and KXD1) were associated with better response to nortriptyline. Considering connectivity between genes, we identified modules of genes that were highly coexpressed. In the whole sample, changes in one of the ten identified coexpression modules showed significant correlation with treatment response (cor=0.27, p=0.0029). Using transcriptomic approaches, we have identified gene expression correlates of the therapeutic effects of antidepressants, highlighting possible molecular pathways involved in efficacious antidepressant treatment.

Association between genetic risk scoring for schizophrenia and bipolar disorder with regional subcortical volumes

Previous research has shown coincident abnormal regional brain volume in patients with schizophrenia (SCZ) and bipolar disorder (BD) compared with controls. Whether these abnormalities are genetically driven or explained by secondary effects of the disorder or environmental factors is unknown. We aimed to investigate the association between genetic risk scoring (GRS) for SCZ and BD with volume of brain areas previously shown to be different between these clinical groups and healthy controls. We obtained subcortical brain volume measures and GRS for SCZ and BD from a sample of 274 healthy volunteers (71.4% females, mean age 24.7 (s.d. 6.9)). Volume of the globus pallidus was associated with the shared GRS between SCZ and BD, and also with the independent GRS for each of these disorders. Volume of the amygdala was associated with the non-shared GRS between SCZ and BD, and with the independent GRS for BD. Our results for volume of the globus pallidus support the idea of SCZ and BD sharing a common underlying neurobiological abnormality associated with a common genetic risk for both these disorders. Results for volume of the amygdala, though, would suggest the existence of a distinct mechanism only associated with genetic risk for BD. Finally, the lack of association between genetic risk and volume of most subcortical structures suggests that the volumetric differences reported in patient-control comparisons may not be genetically driven, but a consequence of the disorder or co-occurring environmental factors.

Associations between polygenic risk for schizophrenia and brain function during probabilistic learning in healthy individuals

A substantial proportion of schizophrenia liability can be explained by additive genetic factors. Risk profile scores (RPS) directly index risk using a summated total of common risk variants weighted by their effect. Previous studies suggest that schizophrenia RPS predict alterations to neural networks that support working memory and verbal fluency. In this study, we apply schizophrenia RPS to fMRI data to elucidate the effects of polygenic risk on functional brain networks during a probabilistic‐learning neuroimaging paradigm. The neural networks recruited during this paradigm have previously been shown to be altered to unmedicated schizophrenia patients and relatives of schizophrenia patients, which may reflect genetic susceptibility. We created schizophrenia RPS using summary data from the Psychiatric Genetic Consortium (Schizophrenia Working Group) for 83 healthy individuals and explore associations between schizophrenia RPS and blood oxygen level dependency (BOLD) during periods of choice behavior (switch–stay) and reflection upon choice outcome (reward–punishment). We show that schizophrenia RPS is associated with alterations in the frontal pole (P_{WHOLE‐BRAIN‐CORRECTED} = 0.048) and the ventral striatum (P_{ROI‐CORRECTED} = 0.036), during choice behavior, but not choice outcome. We suggest that the common risk variants that increase susceptibility to schizophrenia can be associated with alterations in the neural circuitry that support the processing of changing reward contingencies. Hum Brain Mapp 37:491–500, 2016. © 2015 Wiley Periodicals, Inc.

Schizophrenia risk variants modulate white matter volume across the psychosis spectrum: evidence from two independent cohorts

Polygenic risk scores, based on risk variants identified in genome-wide-association-studies (GWAS), explain a considerable portion of the heritability for schizophrenia (SZ) and bipolar disorder (BD). However, little is known about the combined effects of these variants, although polygenic neuroimaging has developed into a powerful tool of translational neuroscience. In this study, we used genome wide significant SZ risk variants to test the predictive capacity of the polygenic model and explored potential associations with white matter volume, a key candidate in imaging phenotype for psychotic disorders. By calculating the combined additive schizophrenia risk of seven SNPs (significant hits from a recent schizophrenia GWAS study), we show that increased additive genetic risk for SZ was associated with reduced white matter volume in a group of participants (n = 94) consisting of healthy individuals, SZ first-degree relatives, SZ patients and BD patients. This effect was also seen in a second independent sample of healthy individuals (n = 89). We suggest that a moderate portion of variance (~4%) of white matter volume can be explained by the seven hits from the recent schizophrenia GWAS. These results provide evidence for associations between cumulative genetic risk for schizophrenia and intermediate neuroimaging phenotypes in models of psychosis. Our work contributes to a growing body of literature suggesting that polygenic risk may help to explain white matter alterations associated with familial risk for psychosis.

Schizophrenia genetics: emerging themes for a complex disorder

After two decades of frustration, genetic studies of schizophrenia have entered an era of spectacular success. Advances in genotyping technologies and high throughput sequencing, increasing analytic rigour and collaborative efforts on a global scale have generated a profusion of new findings. The broad conclusions from these studies are threefold: (1) schizophrenia is a highly polygenic disorder with a complex array of contributing risk loci across the allelic frequency spectrum; (2) many psychiatric illnesses share risk genes and alleles, specifically, schizophrenia has substantial overlaps with bipolar disorder, intellectual disability, major depressive disorder and autism spectrum disorders; and (3) some convergent biological themes are emerging from studies of schizophrenia and related disorders. In this commentary, we focus on the very recent findings that have emerged in the past 12 months, and in particular, the areas of convergence that are beginning to emerge from multiple study designs.

Alzheimer's disease risk variant in CLU is associated with neural inefficiency in healthy individuals

INTRODUCTION

Genome-wide association studies identify rs11136000 in the CLU gene, which codes for Apolipoprotein J/Clusterin, as a significant risk variant for Alzheimer's disease (AD). However, the mechanisms by which this variant confers susceptibility remain relatively unknown.

METHODS

Eighty-five healthy Caucasian participants underwent functional magnetic resonance imaging during a working memory (WM) task and were genotyped for CLU rs11136000/APOE loci.

RESULTS

Here we show that young individuals with the CLU rs11136000 risk variant (C) have higher activation levels in memory-related prefrontal and limbic areas during a WM task. We also found subtle reductions in gray matter in the right hippocampal formation in carriers of the risk variant.

DISCUSSION

We suggest that this pattern of multimodal imaging results may reflect incipient structural differences and inefficient functional activation. This study supports accumulating evidence suggesting that genetic risk for AD affects the neural networks associated with memory in healthy individuals.

An inflammatory biomarker as a differential predictor of outcome of depression treatment with escitalopram and nortriptyline

OBJECTIVE

Major depressive disorder has been linked with inflammatory processes, but it is unclear whether individual differences in levels of inflammatory biomarkers could help match patients to treatments that are most likely to be beneficial. The authors tested the hypothesis that C-reactive protein (CRP), a commonly available marker of systemic inflammation, predicts differential response to escitalopram (a serotonin reuptake inhibitor) and nortriptyline (a norepinephrine reuptake inhibitor).

METHOD

The hypothesis was tested in the Genome-Based Therapeutic Drugs for Depression (GENDEP) study, a multicenter open-label randomized clinical trial. CRP was measured with a high-sensitivity method in serum samples from 241 adult men and women with major depressive disorder randomly allocated to 12-week treatment with escitalopram (N=115) or nortriptyline (N=126). The primary outcome measure was the score on the Montgomery-Åsberg Depression Rating Scale (MADRS), administered weekly.

RESULTS

CRP level at baseline differentially predicted treatment outcome with the two antidepressants (CRP-drug interaction: β=3.27, 95% CI=1.65, 4.89). For patients with low levels of CRP (<1 mg/L), improvement on the MADRS score was 3 points higher with escitalopram than with nortriptyline. For patients with higher CRP levels, improvement on the MADRS score was 3 points higher with nortriptyline than with escitalopram. CRP and its interaction with medication explained more than 10% of individual-level variance in treatment outcome.

CONCLUSIONS

An easily accessible peripheral blood biomarker may contribute to improvement in outcomes of major depressive disorder by personalizing treatment choice.

Genetic relationships between suicide attempts, suicidal ideation and major psychiatric disorders: a genome-wide association and polygenic scoring study

Epidemiological studies have recognized a genetic diathesis for suicidal behavior, which is independent of other psychiatric disorders. Genome-wide association studies (GWAS) on suicide attempt (SA) and ideation have failed to identify specific genetic variants. Here, we conduct further GWAS and for the first time, use polygenic score analysis in cohorts of patients with mood disorders, to test for common genetic variants for mood disorders and suicide phenotypes. Genome-wide studies for SA were conducted in the RADIANT and GSK-Munich recurrent depression samples and London Bipolar Affective Disorder Case-Control Study (BACCs) then meta-analysis was performed. A GWAS on suicidal ideation during antidepressant treatment had previously been conducted in the Genome Based Therapeutic Drugs for Depression (GENDEP) study. We derived polygenic scores from each sample and tested their ability to predict SA in the mood disorder cohorts or ideation status in the GENDEP study. Polygenic scores for major depressive disorder, bipolar disorder and schizophrenia from the Psychiatric Genomics Consortium were used to investigate pleiotropy between psychiatric disorders and suicide phenotypes. No significant evidence for association was detected at any SNP in GWAS or meta-analysis. Polygenic scores for major depressive disorder significantly predicted suicidal ideation in the GENDEP pharmacogenetics study and also predicted SA in a combined validation dataset. Polygenic scores for SA showed no predictive ability for suicidal ideation. Polygenic score analysis suggests pleiotropy between psychiatric disorders and suicidal ideation whereas the tendency to act on such thoughts may have a partially independent genetic diathesis.

Functional polymorphism in the brain-derived neurotrophic factor gene interacts with stressful life events but not childhood maltreatment in the etiology of depression

BACKGROUND

We test the hypothesis that the functional Val66Met polymorphism of BDNF interacts with recent life events to produce onset of new depressive episodes. We also explore the possibility that the Met allele of this polymorphism interacts with childhood maltreatment to increase the risk of chronic depression.

METHODS

In a risk-enriched combined sample of unrelated women, childhood maltreatment and current life events were measured with the Childhood Experience of Care and Abuse, and Life Events and Difficulties Schedule interviews. Chronic episodes of depression (12 months or longer) during adulthood and onset of a major depressive episode during a 12-month follow-up were established with the Schedules for Clinical Assessment in Neuropsychiatry interview.

RESULTS

Met alleles of BDNF moderated the relationship between recent life events and adult onsets of depression in a significant gene-environment interaction (interaction risk difference 0.216, 95% CI 0.090-0.342; P =.0008). BDNF did not significantly influence the effect of childhood maltreatment on chronic depression in the present sample.

CONCLUSIONS

The Met allele of BDNF increases the risk of a new depressive episode following a severe life event. The BDNF and the serotonin transporter gene length polymorphism (5-HTTLPR) and BDNF may contribute to depression through distinct mechanisms involving interactions with childhood and adulthood adversity respectively, which may, in combination, be responsible for a substantial proportion of depression burden in the general population.

Putative transcriptomic biomarkers in the inflammatory cytokine pathway differentiate major depressive disorder patients from control subjects and bipolar disorder patients

Mood disorders consist of two etiologically related, but distinctly treated illnesses, major depressive disorder (MDD) and bipolar disorder (BPD). These disorders share similarities in their clinical presentation, and thus show high rates of misdiagnosis. Recent research has revealed significant transcriptional differences within the inflammatory cytokine pathway between MDD patients and controls, and between BPD patients and controls, suggesting this pathway may possess important biomarker properties. This exploratory study attempts to identify disorder-specific transcriptional biomarkers within the inflammatory cytokine pathway, which can distinguish between control subjects, MDD patients and BPD patients. This is achieved using RNA extracted from subject blood and applying synthesized complementary DNA to quantitative PCR arrays containing primers for 87 inflammation-related genes. Initially, we use ANOVA to test for transcriptional differences in a 'discovery cohort' (total n = 90) and then we use t-tests to assess the reliability of any identified transcriptional differences in a 'validation cohort' (total n = 35). The two most robust and reliable biomarkers identified across both the discovery and validation cohort were Chemokine (C-C motif) ligand 24 (CCL24) which was consistently transcribed higher amongst MDD patients relative to controls and BPD patients, and C-C chemokine receptor type 6 (CCR6) which was consistently more lowly transcribed amongst MDD patients relative to controls. Results detailed here provide preliminary evidence that transcriptional measures within inflammation-related genes might be useful in aiding clinical diagnostic decision-making processes. Future research should aim to replicate findings detailed in this exploratory study in a larger medication-free sample and examine whether identified biomarkers could be used prospectively to aid clinical diagnosis.

The inflammatory cytokines: molecular biomarkers for major depressive disorder?

Cytokines are pleotropic cell signaling proteins that, in addition to their role as inflammatory mediators, also affect neurotransmitter systems, brain functionality and mood. Here we explore the potential utility of cytokine biomarkers for major depressive disorder. Specifically, we explore how genetic, transcriptomic and proteomic information relating to the cytokines might act as biomarkers, aiding clinical diagnosis and treatment selection processes. We advise future studies to investigate whether cytokine biomarkers might differentiate major depressive disorder patients from other patient groups with overlapping clinical characteristics. Furthermore, we invite future pharmacogenetic studies to investigate whether early antidepressant-induced changes to cytokine mRNA or protein levels precede behavioral changes and act as longer-term predictors of clinical antidepressant response.

Interaction between stress and the BDNF Val66Met polymorphism in depression: a systematic review and meta-analysis

BACKGROUND

Major depression is a disabling psychiatric illness with complex origins. Life stress (childhood adversity and recent stressful events) is a robust risk factor for depression. The relationship between life stress and Val66Met polymorphism in the brain-derived neurotrophic factor (BDNF) gene has received much attention. The aim of the present work was to review and conduct a meta-analysis on the results from published studies examining this interaction.

METHODS

A literature search was conducted using PsychINFO and PubMed databases until 22 November 2013. A total of 22 studies with a pooled total of 14,233 participants met the inclusion criteria, the results of which were combined and a meta-analysis performed using the Liptak-Stouffer z-score method.

RESULTS

The results suggest that the Met allele of BDNF Val66Met significantly moderates the relationship between life stress and depression (P = 0.03). When the studies were stratified by type of environmental stressor, the evidence was stronger for an interaction with stressful life events (P = 0.01) and weaker for interaction of BDNF Val66Met with childhood adversity (P = 0.051).

CONCLUSIONS

The interaction between BDNF and life stress in depression is stronger for stressful life events rather than childhood adversity. Methodological limitations of existing studies include poor measurement of life stress.

Genetic susceptibility for bipolar disorder and response to antidepressants in major depressive disorder

The high heterogeneity of response to antidepressant treatment in major depressive disorder (MDD) makes individual treatment outcomes currently unpredictable. It has been suggested that resistance to antidepressant treatment might be due to undiagnosed bipolar disorder or bipolar spectrum features. Here, we investigate the relationship between genetic susceptibility for bipolar disorder and response to treatment with antidepressants in MDD. Polygenic scores indexing risk for bipolar disorder were derived from the Psychiatric Genomics Consortium Bipolar Disorder whole genome association study. Linear regressions tested the effect of polygenic risk scores for bipolar disorder on proportional reduction in depression severity in two large samples of individuals with MDD, treated with antidepressants, NEWMEDS (n=1,791) and STAR*D (n=1,107). There was no significant association between polygenic scores for bipolar disorder and response to treatment with antidepressants. Our data indicate that molecular measure of genetic susceptibility to bipolar disorder does not aid in understanding non-response to antidepressants.

Whole-exome sequencing identifies a polymorphism in the BMP5 gene associated with SSRI treatment response in major depression

Although antidepressants are widely used in the pharmacotherapy of major depressive disorder (MDD), their efficacy is still insufficient as approximately one-third of the patients do not fully recover even after several treatment trials. Inter-individual genetic differences are thought to contribute to the variability in antidepressant response; however, current findings from pharmacogenetic studies are uncertain or not clearly replicated. Here we report the first application of full exome sequencing for the analysis of pharmacogenomics on antidepressant treatment. After 12 weeks of treatment with the selective serotonin re-uptake inhibitor escitalopram, we selected five clear responders and five clear non-responders for exome sequencing. By comparing the allele counts of previously known single nucleotide polymorphisms and novel polymorphisms we selected 38 markers for further genotyping in two independent patient samples treated with escitalopram (n=116 and n=394). The A allele, carried by approximately 30% of the patients with MDD, of rs41271330 in the bone morphogenetic protein (BMP5) gene showed strong association with worse treatment response in both sample sets (p=0.001), indicating that this is an promising pharmacogenetic marker for prediction of antidepressant therapeutic outcome.

ATP-binding cassette sub-family F member 1 (ABCF1) is identified as a putative therapeutic target of escitalopram in the inflammatory cytokine pathway

The inflammatory cytokine pathway may be a potential therapeutic target for major depressive disorder (MDD). Previous reports suggest that antidepressants have anti-inflammatory properties and can cause a reduction in proinflammatory cytokines. Recent evidence suggests this might be mediated at the level of the transcriptome. The current study investigated the transcription of 86 genes in the inflammatory cytokine pathway both at baseline and after eight weeks of escitalopram treatment in MDD patients who were either clinical responders (n=25) or non-responders (n=21), using a subset of samples in the Genome-Based Therapeutic Drugs for Depression project (GENDEP). Changes in expression between baseline and eight weeks of treatment were assessed using two-tailed t-tests. To establish if any significant expression changes related to clinical response, the magnitude of the relative expression change between baseline and eight weeks of treatment was established and binary logistic regressions were used to compare differences between responders and non-responders. ATP-binding cassette sub-family F member 1 (ABCF1), a translational regulator of the inflammatory cytokine pathway showed a significant increase in expression after escitalopram treatment which was significantly greater in responders compared to non-responders, suggesting that ABCF1 may play a role in mediating antidepressant response.

Biomarkers predicting treatment outcome in depression: what is clinically significant?

AIM

To extend to biomarker studies the consensus clinical significance criterion of a three-point difference in Hamilton Rating Scale for Depression.

MATERIALS & METHODS

We simulated datasets modeled on large clinical trials.

RESULTS

In a typical clinical trial comparing active treatment and placebo, a difference of three Hamilton Rating Scale for Depression (HRSD) points at the end of treatment corresponds to 6.3% of variance in outcome explained. To achieve a similar explanatory power, genotypes with minor allele frequencies of 5, 10, 20, 30 and 50% need to attain a per allele difference of 4.7, 3.6, 2.8, 2.4 and 2.2 HRSD points, respectively. A normally distributed continuous biomarker will need an effect size of 1.5 HRSD points per standard deviation. A number needed to assess of three suggests that with this effect size, a biomarker will significantly improve the prediction of outcome in one out of every three patients assessed.

CONCLUSION

This report provides guidance on assessing clinical significance of biomarkers predictive of outcome in depression treatment.