SATB2 organizes the 3D genome architecture of cognition in cortical neurons

The DNA-binding protein SATB2 is genetically linked to human intelligence. We studied its influence on the three-dimensional (3D) epigenome by mapping chromatin interactions and accessibility in control versus SATB2-deficient cortical neurons. We find that SATB2 affects the chromatin looping between enhancers and promoters of neuronal-activity-regulated genes, thus influencing their expression. It also alters A/B compartments, topologically associating domains, and frequently interacting regions. Genes linked to SATB2-dependent 3D genome changes are implicated in highly specialized neuronal functions and contribute to cognitive ability and risk for neuropsychiatric and neurodevelopmental disorders. Non-coding DNA regions with a SATB2-dependent structure are enriched for common variants associated with educational attainment, intelligence, and schizophrenia. Our data establish SATB2 as a cell-type-specific 3D genome modulator, which operates both independently and in cooperation with CCCTC-binding factor (CTCF) to set up the chromatin landscape of pyramidal neurons for cognitive processes.

Corpus callosum microstructural organization mediates the effects of physical neglect on social cognition in schizophrenia

Exposure to early life adversity is associated with both increased risk of developing schizophrenia and poorer performance on measures of social cognitive functioning. In this study, we examined whether interleukin-6 (IL-6) and Corpus Callosum (CC) microstructure mediated the association between childhood physical neglect and social cognition. Fifty-eight patients with a diagnosis of schizophrenia were included. The CANTAB emotion recognition task (unbiased hit rate) was used to assess social cognition. We found that the microstructural organization of the CC significantly mediated the association between physical neglect and emotion recognition. Furthermore, in a sequential mediation analysis that also considered the role of inflammatory response, the association between physical neglect, and lower emotion recognition performance was sequentially mediated by higher IL-6 and lower fractional anisotropy of the CC. This mediating effect of IL-6 was only present when simultaneously considering the effects of CC microstructural organization and remained significant while controlling for the effects of sex, BMI and medication dosage (but not age). Overall, the findings suggest that the association between physical neglect and poorer emotion recognition in schizophrenia occurs, at least in part, via its association with white matter microstructure.

Neuropsychiatric disorders, chronotype and sleep: A narrative review of GWAS findings and the application of Mendelian randomization to investigate causal relationships

Genome-wide association studies (GWAS) have been important for characterizing the genetic component and enhancing our understanding of the biological aetiology of both neuropsychiatric disorders and sleep-related phenotypes such as chronotype, which is our preference for morning or evening time. Mendelian randomization (MR) is a post-GWAS analysis that is used to infer causal relationships between potential risk factors and outcomes. MR uses genetic variants as instrumental variants for exposures to study the effect on outcomes. This review details the main results from GWAS of neuropsychiatric disorders and sleep-related phenotypes, and the application of MR to investigate their bidirectional relationship. The main results from MR studies of neuropsychiatric disorders and sleep-related phenotypes are summarized. These MR studies have identified 37 causal relationships between neuropsychiatric disorders and sleep-related phenotypes. MR studies identified evidence of a causal role for five neuropsychiatric disorders and symptoms (attention deficit hyperactivity disorder, bipolar disorder, depressive symptoms, major depressive disorder and schizophrenia) on sleep-related phenotypes and evidence of a causal role for five sleep-related phenotypes (daytime napping, insomnia, morning person, long sleep duration and sleep duration) on risk for neuropsychiatric disorders. These MR results show a bidirectional relationship between neuropsychiatric disorders and sleep-related phenotypes and identify potential risk factors for follow-up studies.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers

BACKGROUND

Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure.

METHODS

Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference.

RESULTS

For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness.

CONCLUSIONS

We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Genetic and inflammatory effects on childhood trauma and cognitive functioning in patients with schizophrenia and healthy participants

Recent studies have reported a negative association between exposure to childhood trauma, including physical neglect, and cognitive functioning in patients with schizophrenia. Childhood trauma has been found to influence immune functioning, which may contribute to the risk of schizophrenia and cognitive symptoms of the disorder. In this study, we aimed to test the hypothesis that physical neglect is associated with cognitive ability, and that this association is mediated by a combined latent measure of inflammatory response, and moderated by higher genetic risk for schizophrenia. The study included 279 Irish participants, comprising 102 patients and 177 healthy participants. Structural equation modelling was used to perform mediation and moderation analyses. Inflammatory response was measured via basal plasma levels of IL-6, TNF-α, and CRP, and cognitive performance was assessed across three domains: full-scale IQ, logical memory, and the emotion recognition task. Genetic variation for schizophrenia was estimated using a genome-wide polygenic score based on genome-wide association study summary statistics. The results showed that inflammatory response mediated the association between physical neglect and all measures of cognitive functioning, and explained considerably more variance than any of the inflammatory markers alone. Furthermore, genetic risk for schizophrenia was observed to moderate the direct pathway between physical neglect and measures of non-social cognitive functioning in both patient and healthy participants. However, genetic risk did not moderate the mediated pathway associated with inflammatory response. Therefore, we conclude that the mediating role of inflammatory response and the moderating role of higher genetic risk may independently influence the association between adverse early life experiences and cognitive function in patients and healthy participants.

Social isolation-induced transcriptomic changes in mouse hippocampus impact the synapse and show convergence with human genetic risk for neurodevelopmental phenotypes

Early life stress (ELS) can impact brain development and is a risk factor for neurodevelopmental disorders such as schizophrenia. Post-weaning social isolation (SI) is used to model ELS in animals, using isolation stress to disrupt a normal developmental trajectory. We aimed to investigate how SI affects the expression of genes in mouse hippocampus and to investigate how these changes related to the genetic basis of neurodevelopmental phenotypes. BL/6J mice were exposed to post-weaning SI (PD21-25) or treated as group-housed controls (n = 7-8 per group). RNA sequencing was performed on tissue samples from the hippocampus of adult male and female mice. Four hundred and 1,215 differentially-expressed genes (DEGs) at a false discovery rate of < 0.05 were detected between SI and control samples for males and females respectively. DEGS for both males and females were significantly overrepresented in gene ontologies related to synaptic structure and function, especially the post-synapse. DEGs were enriched for common variant (SNP) heritability in humans that contributes to risk of neuropsychiatric disorders (schizophrenia, bipolar disorder) and to cognitive function. DEGs were also enriched for genes harbouring rare de novo variants that contribute to autism spectrum disorder and other developmental disorders. Finally, cell type analysis revealed populations of hippocampal astrocytes that were enriched for DEGs, indicating effects in these cell types as well as neurons. Overall, these data suggest a convergence between genes dysregulated by the SI stressor in the mouse and genes associated with neurodevelopmental disorders and cognitive phenotypes in humans.

Bilateral anterior corona radiata microstructure organisation relates to impaired social cognition in schizophrenia

OBJECTIVE

The Corona Radiata (CR) is a large white matter tract in the brain comprising of the anterior CR (aCR), superior CR (sCR), and posterior CR (pCR), which have associations with cognition, self-regulation, and, in schizophrenia, positive symptom severity. This study tested the hypothesis that the microstructural organisation of the aCR, as measured by Fractional Anisotropy (FA) using Diffusion Tensor Imaging (DTI), would relate to poorer social cognitive outcomes and higher positive symptom severity for people with schizophrenia, when compared to healthy participants. We further hypothesised that increased positive symptoms would relate to poorer social cognitive outcomes.

METHODS

Data were derived from n = 178 healthy participants (41 % females; 36.11 ± 12.36 years) and 58 people with schizophrenia (30 % females; 42.4 ± 11.1 years). The Positive and Negative Symptom Severity Scale measured clinical symptom severity. Social Cognition was measured using the Reading the Mind in the Eyes Test (RMET) Total Score, as well as the Positive, Neutral, and Negative stimuli valence. The ENIGMA-DTI protocol tract-based spatial statistics (TBSS) was used.

RESULTS

There was a significant difference in FA for the CR, in individuals with schizophrenia compared to healthy participants. On stratification, both the aCR and pCR were significantly different between groups, with patients showing reduced white matter tract microstructural organisation. Significant negative correlations were observed between positive symptomatology and reduced microstructural organisation of the aCR. Performance for RMET negative valence items was significantly correlated bilaterally with the aCR, but not the sCR or pCR, and no relationship to positive symptoms was observed.

CONCLUSIONS

These data highlight specific and significant microstructural white-matter differences for people with schizophrenia, which relates to positive clinical symptomology and poorer performance on social cognition stimuli. While reduced FA is associated with higher positive symptomatology in schizophrenia, this study shows the specific associated with anterior frontal white matter tracts and reduced social cognitive performance. The aCR may have a specific role to play in frontal-disconnection syndromes, psychosis, and social cognitive profile within schizophrenia, though further research requires more sensitive, specific, and detailed consideration of social cognition outcomes.

Genes positively regulated by Mef2c in cortical neurons are enriched for common genetic variation associated with IQ and educational attainment

The myocyte enhancer factor 2 C (MEF2C) gene encodes a transcription factor important for neurogenesis and synapse development and contains common variants associated with intelligence (IQ) and educational attainment (EA). Here, we took gene expression data from the mouse cortex of a Mef2c mouse model with a heterozygous DNA binding-deficient mutation of Mef2c (Mef2c-het) and combined these data with MEF2C ChIP-seq data from cortical neurons and single-cell data from the mouse brain. This enabled us to create a set of genes that were differentially regulated in Mef2c-het mice, represented direct target genes of MEF2C and had elevated in expression in cortical neurons. We found this gene-set to be enriched for genes containing common genetic variation associated with IQ and EA. Genes within this gene-set that were down-regulated, i.e. have reduced expression in Mef2c-het mice versus controls, were specifically significantly enriched for both EA and IQ associated genes. These down-regulated genes were enriched for functionality in the adenylyl cyclase signalling system, which is known to positively regulate synaptic transmission and has been linked to learning and memory. Within the adenylyl cyclase signalling system, three genes regulated by MEF2C, CRHR1, RGS6, and GABRG3, are associated at genome-wide significant levels with IQ and/or EA. Our results indicate that genetic variation in MEF2C and its direct target genes within cortical neurons contribute to variance in cognition within the general population, and the molecular mechanisms involved include the adenylyl cyclase signalling system's role in synaptic function.

The impact of early adversity and education on genetic and brain morphological predictors of cognitive ability

Cognitive ability is a strong predictor of occupational achievement, quality of life and physical health. While variation in cognition is strongly heritable and has been robustly associated with early environment and brain morphology, little is known about how these factors combine and interact to explain this variation in cognition. To address this, we modelled the relationship between common genetic variation, grey matter volume, early life adversity and education and cognitive ability in a UK Biobank sample of N = 5237 individuals using structural equation modelling. We tested the hypotheses that total grey matter volume would mediate the association between genetic variation and cognitive ability, and that early life adversity and educational attainment would moderate this relationship. Common genetic variation, grey matter volume and early life adversity were each significant predictors in the model, explaining ~15% of variation in cognitive ability. Contrary to our hypothesis, grey matter volume did not mediate the relation between genetic variation and cognition performance. Neither did early life adversity or educational attainment moderate this relation, although educational attainment was observed to moderate the relationship between grey matter volume and cognitive performance. We interpret these findings in terms of the modest explanatory value of currently estimated polygenic scores accounting for variation in cognitive performance (~5%), making potential mediating and moderating variables difficult to confirm.

Early life stress, low-grade systemic inflammation and weaker suppression of the default mode network (DMN) during face processing in Schizophrenia

Childhood trauma (CT) is associated with lower cognitive and social cognitive function in schizophrenia. Recent evidence suggests that the relationship between CT and cognition is mediated by both low-grade systemic inflammation and reduced connectivity of the default mode network (DMN) during resting state. This study sought to test whether the same pattern of associations was observed for DMN connectivity during task based activity. Fifty-three individuals with schizophrenia (SZ) or schizoaffective disorder (SZA) and one hundred and seventy six healthy participants were recruited from the Immune Response and Social Cognition (iRELATE) project. A panel of pro-inflammatory markers that included IL-6, IL-8, IL-10, tumour necrosis factor-alpha (TNFa), and C-reactive protein (CRP), were measured in plasma using ELISA. DMN connectivity was measured during an fMRI social cognitive face processing task. Patients showed evidence of low grade systemic inflammation and significantly increased connectivity between the left lateral parietal (LLP) cortex-cerebellum and LLP-left angular gyrus compared to healthy participants. Across the entire sample, IL-6 predicted increased connectivity between LLP-cerebellum, LLP-precuneus, and mPFC-bilateral-precentral-gyri and left postcentral gyrus. In turn, and again in the entire sample, IL-6 (but no other inflammatory marker) mediated the relationship between childhood physical neglect and LLP-cerebellum. Physical neglect scores also significantly predicted the positive association between IL-6 and LLP-precuneus connectivity. This is to our knowledge the first study that provides evidence that higher plasma IL-6 mediates the association between higher childhood neglect and increased DMN connectivity during task based activity. Consistent with our hypothesis, exposure to trauma is associated with weaker suppression of the DMN during a face processing task, and this association was mediated via increased inflammatory response. The findings may represent part of the biological mechanism by which CT and cognitive performance are related.

Independent Associated SNPs at SORCS3 and Its Protein Interactors for Multiple Brain-Related Disorders and Traits

Sortilin-related vacuolar protein sorting 10 (VPS10) domain containing receptor 3 (SORCS3) is a neuron-specific transmembrane protein involved in the trafficking of proteins between intracellular vesicles and the plasma membrane. Genetic variation at SORCS3 is associated with multiple neuropsychiatric disorders and behavioural phenotypes. Here, we undertake a systematic search of published genome-wide association studies to identify and catalogue associations between SORCS3 and brain-related disorders and traits. We also generate a SORCS3 gene-set based on protein-protein interactions and investigate the contribution of this gene-set to the heritability of these phenotypes and its overlap with synaptic biology. Analysis of association signals at SORSC3 showed individual SNPs to be associated with multiple neuropsychiatric and neurodevelopmental brain-related disorders and traits that have an impact on the experience of feeling, emotion or mood or cognitive function, while multiple LD-independent SNPs were associated with the same phenotypes. Across these SNPs, alleles associated with the more favourable outcomes for each phenotype (e.g., decreased risk of neuropsychiatric illness) were associated with increased expression of the SORCS3 gene. The SORCS3 gene-set was enriched for heritability contributing to schizophrenia (SCZ), bipolar disorder (BPD), intelligence (IQ) and education attainment (EA). Eleven genes from the SORCS3 gene-set were associated with more than one of these phenotypes at the genome-wide level, with RBFOX1 associated with SCZ, IQ and EA. Functional annotation revealed that the SORCS3 gene-set is enriched for multiple ontologies related to the structure and function of synapses. Overall, we find many independent association signals at SORCS3 with brain-related disorders and traits, with the effect possibly mediated by reduced gene expression, resulting in a negative impact on synaptic function.

Whole blood transcriptome analysis in dairy calves experimentally challenged with bovine herpesvirus 1 (BoHV-1) and comparison to a bovine respiratory syncytial virus (BRSV) challenge

Bovine herpesvirus 1 (BoHV-1), is associated with several clinical syndromes in cattle, among which bovine respiratory disease (BRD) is of particular significance. Despite the importance of the disease, there is a lack of information on the molecular response to infection via experimental challenge with BoHV-1. The objective of this study was to investigate the whole-blood transcriptome of dairy calves experimentally challenged with BoHV-1. A secondary objective was to compare the gene expression results between two separate BRD pathogens using data from a similar challenge study with BRSV. Holstein-Friesian calves (mean age (SD) = 149.2 (23.8) days; mean weight (SD) = 174.6 (21.3) kg) were either administered BoHV-1 inoculate (1 × 10⁷/mL × 8.5 mL) (n = 12) or were mock challenged with sterile phosphate buffered saline (n = 6). Clinical signs were recorded daily from day (d) -1 to d 6 (post-challenge), and whole blood was collected in Tempus RNA tubes on d six post-challenge for RNA-sequencing. There were 488 differentially expressed (DE) genes (p < 0.05, False Discovery rate (FDR) < 0.10, fold change ≥2) between the two treatments. Enriched KEGG pathways (p < 0.05, FDR <0.05); included Influenza A, Cytokine-cytokine receptor interaction and NOD-like receptor signalling. Significant gene ontology terms (p < 0.05, FDR <0.05) included defence response to virus and inflammatory response. Genes that are highly DE in key pathways are potential therapeutic targets for the treatment of BoHV-1 infection. A comparison to data from a similar study with BRSV identified both similarities and differences in the immune response to differing BRD pathogens.

Genetic variants associated with longitudinal changes in brain structure across the lifespan

Human brain structure changes throughout the lifespan. Altered brain growth or rates of decline are implicated in a vast range of psychiatric, developmental and neurodegenerative diseases. In this study, we identified common genetic variants that affect rates of brain growth or atrophy in what is, to our knowledge, the first genome-wide association meta-analysis of changes in brain morphology across the lifespan. Longitudinal magnetic resonance imaging data from 15,640 individuals were used to compute rates of change for 15 brain structures. The most robustly identified genes GPR139, DACH1 and APOE are associated with metabolic processes. We demonstrate global genetic overlap with depression, schizophrenia, cognitive functioning, insomnia, height, body mass index and smoking. Gene set findings implicate both early brain development and neurodegenerative processes in the rates of brain changes. Identifying variants involved in structural brain changes may help to determine biological pathways underlying optimal and dysfunctional brain development and aging.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Rare coding variants in ten genes confer substantial risk for schizophrenia

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

Thirteen Independent Genetic Loci Associated with Preserved Processing Speed in a Study of Cognitive Resilience in 330,097 Individuals in the UK Biobank

Cognitive resilience is the ability to withstand the negative effects of stress on cognitive functioning and is important for maintaining quality of life while aging. The UK Biobank does not have measurements of the same cognitive phenotype at distal time points. Therefore, we used education years (EY) as a proxy phenotype for past cognitive performance and current cognitive performance was based on processing speed. This represented an average time span of 40 years between past and current cognitive performance in 330,097 individuals. A confounding factor was that EY is highly polygenic and masked the genetics of resilience. To overcome this, we employed Genomics Structural Equation Modelling (GenomicSEM) to perform a genome-wide association study (GWAS)-by-subtraction using two GWAS, one GWAS of EY and resilience and a second GWAS of EY but not resilience, to generate a GWAS of Resilience. Using independent discovery and replication samples, we found 13 independent genetic loci for Resilience. Functional analyses showed enrichment in several brain regions and specific cell types. Gene-set analyses implicated the biological process "neuron differentiation", the cellular component "synaptic part" and the "WNT signalosome". Mendelian randomisation analysis showed a causative effect of white matter volume on cognitive resilience. These results may contribute to the neurobiological understanding of resilience.

Current psychosocial stress, childhood trauma and cognition in patients with schizophrenia and healthy participants

BACKGROUND

Cognitive difficulties are experienced frequently in schizophrenia (SZ) and are strongly predictive of functional outcome. Although severity of cognitive difficulties has been robustly associated with early life adversity, whether and how they are affected by current stress is unknown. The present study investigated whether acute stress reactivity as measured by heart rate and mood changes predict cognitive performance in patients with schizophrenia and healthy individuals, and whether this is moderated by diagnosis and previous childhood trauma exposure.

METHODS

One hundred and four patients with schizophrenia and 207 healthy participants were administered a battery of tasks assessing cognitive performance after psychosocial stress induction (Trier Social Stress Test; TSST). Mood states (Profile of Mood States; POMS) and heart rate were assessed at baseline, immediately before, and after the TSST.

RESULTS

Both healthy participants and patients showed increases in POMS Tension and Total Mood Disturbance scores between Time Point 2 (pre-TSST) and Time Point 3 (post-TSST). These changes were not associated with variation in cognition. Although childhood trauma exposure was associated with higher stress reactivity and poorer cognitive function in all participants, childhood trauma did not moderate the association between stress reactivity and cognition. Neither was diagnosis a moderator of this relationship.

DISCUSSION

These findings suggest that while chronic stress exposure explains significant variation in cognition, acute stress reactivity (measured by changes in Tension and Total Mood Disturbance) did not. In the context of broader developmental processes, we conclude that stressful events that occur earlier in development, and with greater chronicity, are likely to be more strongly associated with cognitive variation than acute transient stressors experienced in adulthood.

Microglial-expressed genetic risk variants, cognitive function and brain volume in patients with schizophrenia and healthy controls

Changes in immune function are associated with variance in cognitive functioning in schizophrenia. Given that microglia are the primary innate immune cells in the brain, we examined whether schizophrenia risk-associated microglial genes (measured via polygenic score analysis) explained variation in cognition in patients with schizophrenia and controls (n = 1,238) and tested whether grey matter mediated this association. We further sought to replicate these associations in an independent sample of UK Biobank participants (n = 134,827). We then compared the strength of these microglial associations to that of neuronal and astroglial (i.e., other brain-expressed genes) polygenic scores, and used MAGMA to test for enrichment of these gene-sets with schizophrenia risk. Increased microglial schizophrenia polygenic risk was associated with significantly lower performance across several measures of cognitive functioning in both samples; associations which were then found to be mediated via total grey matter volume in the UK Biobank. Unlike neuronal genes which did show evidence of enrichment, the microglial gene-set was not significantly enriched for schizophrenia, suggesting that the relevance of microglia may be for neurodevelopmental processes related more generally to cognition. Further, the microglial polygenic score was associated with performance on a range of cognitive measures in a manner comparable to the neuronal schizophrenia polygenic score, with fewer cognitive associations observed for the astroglial score. In conclusion, our study supports the growing evidence of the importance of immune processes to understanding cognition and brain structure in both patients and in the healthy population.

Meta-Analysis of Brain Gene Expression Data from Mouse Model Studies of Maternal Immune Activation Using Poly(I:C)

Maternal immune activation (MIA) is a known risk factor for schizophrenia (SCZ) and autism spectrum disorder (ASD) and is often modelled in animal studies in order to study the effect of prenatal infection on brain function including behaviour and gene expression. Although the effect of MIA on gene expression are highly heterogeneous, combining data from multiple gene expression studies in a robust method may shed light on the true underlying biological effects caused by MIA and this could inform studies of SCZ and ASD. This study combined four RNA-seq and microarray datasets in an overlap analysis and ranked meta-analysis in order to investigate genes, pathways and cell types dysregulated in the MIA mouse models. Genes linked to SCZ and ASD and crucial in neurodevelopmental processes including neural tube folding, regulation of cellular stress and neuronal/glial cell differentiation were among the most consistently dysregulated in these ranked analyses. Gene ontologies including K+ ion channel function, neuron and glial cell differentiation, synaptic structure, axonal outgrowth, cilia function and lipid metabolism were also strongly implicated. Single-cell analysis identified excitatory and inhibitory cell types in the cortex, hippocampus and striatum that may be affected by MIA and are also enriched for genes associated with SCZ, ASD and cognitive phenotypes. This points to the cellular location of molecular mechanisms that may be consistent between the MIA model and neurodevelopmental disease, improving our understanding of its utility to study prenatal infection as an environmental stressor.

Childhood trauma, brain structure and emotion recognition in patients with schizophrenia and healthy participants

Childhood trauma, and in particular physical neglect, has been repeatedly associated with lower performance on measures of social cognition (e.g. emotion recognition tasks) in both psychiatric and non-clinical populations. The neural mechanisms underpinning this association have remained unclear. Here, we investigated whether volumetric changes in three stress-sensitive regions—the amygdala, hippocampus and anterior cingulate cortex (ACC)—mediate the association between childhood trauma and emotion recognition in a healthy participant sample (N = 112) and a clinical sample of patients with schizophrenia (N = 46). Direct effects of childhood trauma, specifically physical neglect, on Emotion Recognition Task were observed in the whole sample. In healthy participants, reduced total and left ACC volumes were observed to fully mediate the association between both physical neglect and total childhood trauma score, and emotion recognition. No mediating effects of the hippocampus and amygdala volumes were observed for either group. These results suggest that reduced ACC volume may represent part of the mechanism by which early life adversity results in poorer social cognitive function. Confirmation of the causal basis of this association would highlight the importance of resilience-building interventions to mitigate the detrimental effects of childhood trauma on brain structure and function.

Altered gene regulation as a candidate mechanism by which ciliopathy gene SDCCAG8 contributes to schizophrenia and cognitive function

Mutations in genes that encode centrosomal/ciliary proteins cause severe cognitive deficits, while common single-nucleotide polymorphisms in these genes are associated with schizophrenia (SZ) and cognition in genome-wide association studies. The role of these genes in neuropsychiatric disorders is unknown. The ciliopathy gene SDCCAG8 is associated with SZ and educational attainment (EA). Genome editing of SDCCAG8 caused defects in primary ciliogenesis and cilium-dependent cell signalling. Transcriptomic analysis of SDCCAG8-deficient cells identified differentially expressed genes that are enriched in neurodevelopmental processes such as generation of neurons and synapse organization. These processes are enriched for genes associated with SZ, human intelligence (IQ) and EA. Phenotypic analysis of SDCCAG8-deficent neuronal cells revealed impaired migration and neuronal differentiation. These data implicate ciliary signalling in the aetiology of SZ and cognitive dysfunction. We found that centrosomal/ciliary genes are enriched for association with IQ, suggesting altered gene regulation as a general model for neurodevelopmental impacts of centrosomal/ciliary genes.

Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

Changes in Default-Mode Network Associated With Childhood Trauma in Schizophrenia

BACKGROUND

There is considerable evidence of dysconnectivity within the default-mode network (DMN) in schizophrenia, as measured during resting-state functional MRI (rs-fMRI). History of childhood trauma (CT) is observed at a higher frequency in schizophrenia than in the general population, but its relationship to DMN functional connectivity has yet to be investigated.

METHODS

CT history and rs-fMRI data were collected in 65 individuals with schizophrenia and 132 healthy controls. Seed-based functional connectivity between each of 4 a priori defined seeds of the DMN (medial prefrontal cortex, right and left lateral parietal lobes, and the posterior cingulate cortex) and all other voxels of the brain were compared across groups. Effects of CT on functional connectivity were examined using multiple regression analyses. Where significant associations were observed, regression analyses were further used to determine whether variance in behavioral measures of Theory of Mind (ToM), previously associated with DMN recruitment, were explained by these associations.

RESULTS

Seed-based analyses revealed evidence of widespread reductions in functional connectivity in patients vs controls, including between the left/right parietal lobe (LP) and multiple other regions, including the parietal operculum bilaterally. Across all subjects, increased CT scores were associated with reduced prefrontal-parietal connectivity and, in patients, with increased prefrontal-cerebellar connectivity also. These CT-associated differences in DMN connectivity also predicted variation in behavioral measures of ToM.

CONCLUSIONS

These findings suggest that CT history is associated with variation in DMN connectivity during rs-fMRI in patients with schizophrenia and healthy participants, which may partly mediate associations observed between early life adversity and cognitive performance.

1q21.1 distal copy number variants are associated with cerebral and cognitive alterations in humans

Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMA-CNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers-the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function.

DNA methylation meta-analysis reveals cellular alterations in psychosis and markers of treatment-resistant schizophrenia

We performed a systematic analysis of blood DNA methylation profiles from 4483 participants from seven independent cohorts identifying differentially methylated positions (DMPs) associated with psychosis, schizophrenia, and treatment-resistant schizophrenia. Psychosis cases were characterized by significant differences in measures of blood cell proportions and elevated smoking exposure derived from the DNA methylation data, with the largest differences seen in treatment-resistant schizophrenia patients. We implemented a stringent pipeline to meta-analyze epigenome-wide association study (EWAS) results across datasets, identifying 95 DMPs associated with psychosis and 1048 DMPs associated with schizophrenia, with evidence of colocalization to regions nominated by genetic association studies of disease. Many schizophrenia-associated DNA methylation differences were only present in patients with treatment-resistant schizophrenia, potentially reflecting exposure to the atypical antipsychotic clozapine. Our results highlight how DNA methylation data can be leveraged to identify physiological (e.g., differential cell counts) and environmental (e.g., smoking) factors associated with psychosis and molecular biomarkers of treatment-resistant schizophrenia.

SATB2-LEMD2 interaction links nuclear shape plasticity to regulation of cognition-related genes

SATB2 is a schizophrenia risk gene and is genetically associated with human intelligence. How it affects cognition at molecular level is currently unknown. Here, we show that interactions between SATB2, a chromosomal scaffolding protein, and the inner nuclear membrane protein LEMD2 orchestrate the response of pyramidal neurons to neuronal activation. Exposure to novel environment in vivo causes changes in nuclear shape of CA1 hippocampal neurons via a SATB2-dependent mechanism. The activity-driven plasticity of the nuclear envelope requires not only SATB2, but also its protein interactor LEMD2 and the ESCRT-III/VPS4 membrane-remodeling complex. Furthermore, LEMD2 depletion in cortical neurons, similar to SATB2 ablation, affects neuronal activity-dependent regulation of multiple rapid and delayed primary response genes. In human genetic data, LEMD2-regulated genes are enriched for de novo mutations reported in intellectual disability and schizophrenia and are, like SATB2-regulated genes, enriched for common variants associated with schizophrenia and cognitive function. Hence, interactions between SATB2 and the inner nuclear membrane protein LEMD2 influence gene expression programs in pyramidal neurons that are linked to cognitive ability and psychiatric disorder etiology.

Effects of complement gene-set polygenic risk score on brain volume and cortical measures in patients with psychotic disorders and healthy controls

Multiple genome-wide association studies of schizophrenia have reported associations between genetic variants within the MHC region and disease risk, an association that has been partially accounted for by alleles of the complement component 4 (C4) gene. Following on previous findings of association between both C4 and other complement-related variants and memory function, we tested the hypothesis that polygenic scores calculated based on identified schizophrenia risk alleles within the "complement" system would be broadly associated with memory function and associated brain structure. We tested this using a polygenic risk score (PRS) calculated for complement genes, but excluding C4 variants. Higher complement-based PRS scores were observed to be associated with lower memory scores for the sample as a whole (N = 620, F change = 8.25; p = .004). A significant association between higher PRS and lower hippocampal volume was also observed (N = 216, R² change = 0.016, p = .015). However, after correcting for further testing of association with the more general indices of cortical thickness, surface area or total brain volume, none of which were associated with complement, the association with hippocampal volume became non-significant. A post-hoc analysis of hippocampal subfields suggested an association between complement PRS and several hippocampal subfields, findings that appeared to be particularly driven by the patient sample. In conclusion, our study yielded suggestive evidence of association between complement-based schizophrenia PRS and variation in memory function and hippocampal volume.

Genes influenced by MEF2C contribute to neurodevelopmental disease via gene expression changes that affect multiple types of cortical excitatory neurons

Myocyte enhancer factor 2 C (MEF2C) is an important transcription factor during neurodevelopment. Mutation or deletion of MEF2C causes intellectual disability (ID), and common variants within MEF2C are associated with cognitive function and schizophrenia risk. We investigated if genes influenced by MEF2C during neurodevelopment are enriched for genes associated with neurodevelopmental phenotypes and if this can be leveraged to identify biological mechanisms and individual brain cell types affected. We used a set of 1055 genes that were differentially expressed in the adult mouse brain following early embryonic deletion of Mef2c in excitatory cortical neurons. Using genome-wide association studies data, we found these differentially expressed genes (DEGs) to be enriched for genes associated with schizophrenia, intelligence and educational attainment but not autism spectrum disorder (ASD). For this gene set, genes that overlap with target genes of the Fragile X mental retardation protein (FMRP) are a major driver of these enrichments. Using trios data, we found these DEGs to be enriched for genes containing de novo mutations reported in ASD and ID, but not schizophrenia. Using single-cell RNA sequencing data, we identified that a number of different excitatory glutamatergic neurons in the cortex were enriched for these DEGs including deep layer pyramidal cells and cells in the retrosplenial cortex, entorhinal cortex and subiculum, and these cell types are also enriched for FMRP target genes. The involvement of MEF2C and FMRP in synapse elimination suggests that disruption of this process in these cell types during neurodevelopment contributes to cognitive function and risk of neurodevelopmental disorders.

Genes regulated by BCL11B during T-cell development are enriched for de novo mutations found in schizophrenia patients

While abnormal neurodevelopment contributes to schizophrenia (SCZ) risk, there is also evidence to support a role for immune dysfunction in SCZ. BCL11B, associated with SCZ in genome-wide association study (GWAS), is a transcription factor that regulates the differentiation and development of cells in the central nervous and immune systems. Here, we use functional genomics data from studies of BCL11B to investigate the contribution of neuronal and immune processes to SCZ pathophysiology. We identified the gene targets of BCL11B in brain striatal cells (n = 223 genes), double negative 4 (DN4) developing T cells (n = 114 genes) and double positive (DP) developing T cells (n = 518 genes) using an integrated analysis of RNA-seq and ChIP-seq data. No gene-set was enriched for genes containing common variants associated with SCZ but the DP gene-set was enriched for genes containing missense de novo mutations (DNMs; p = .001) using data from 3,447 SCZ trios. Post hoc analysis revealed the enrichment to be stronger for DP genes negatively regulated by BCL11B. Biological processes enriched for genes negatively regulated by BCL11B in DP gene-set included immune system development and cytokine signaling. These analyses, leveraging a GWAS-identified SCZ risk gene and data on gene expression and transcription factor binding, indicate that DNMs in immune pathways contribute to SCZ risk.

Association of Copy Number Variation of the 15q11.2 BP1-BP2 Region With Cortical and Subcortical Morphology and Cognition

Importance

Recurrent microdeletions and duplications in the genomic region 15q11.2 between breakpoints 1 (BP1) and 2 (BP2) are associated with neurodevelopmental disorders. These structural variants are present in 0.5% to 1.0% of the population, making 15q11.2 BP1-BP2 the site of the most prevalent known pathogenic copy number variation (CNV). It is unknown to what extent this CNV influences brain structure and affects cognitive abilities.

Objective

To determine the association of the 15q11.2 BP1-BP2 deletion and duplication CNVs with cortical and subcortical brain morphology and cognitive task performance.

Design, Setting, and Participants

In this genetic association study, T1-weighted brain magnetic resonance imaging were combined with genetic data from the ENIGMA-CNV consortium and the UK Biobank, with a replication cohort from Iceland. In total, 203 deletion carriers, 45 247 noncarriers, and 306 duplication carriers were included. Data were collected from August 2015 to April 2019, and data were analyzed from September 2018 to September 2019.

Main Outcomes and Measures

The associations of the CNV with global and regional measures of surface area and cortical thickness as well as subcortical volumes were investigated, correcting for age, age2, sex, scanner, and intracranial volume. Additionally, measures of cognitive ability were analyzed in the full UK Biobank cohort.

Results

Of 45 756 included individuals, the mean (SD) age was 55.8 (18.3) years, and 23 754 (51.9%) were female. Compared with noncarriers, deletion carriers had a lower surface area (Cohen d = -0.41; SE, 0.08; P = 4.9 × 10-8), thicker cortex (Cohen d = 0.36; SE, 0.07; P = 1.3 × 10-7), and a smaller nucleus accumbens (Cohen d = -0.27; SE, 0.07; P = 7.3 × 10-5). There was also a significant negative dose response on cortical thickness (β = -0.24; SE, 0.05; P = 6.8 × 10-7). Regional cortical analyses showed a localization of the effects to the frontal, cingulate, and parietal lobes. Further, cognitive ability was lower for deletion carriers compared with noncarriers on 5 of 7 tasks.

Conclusions and Relevance

These findings, from the largest CNV neuroimaging study to date, provide evidence that 15q11.2 BP1-BP2 structural variation is associated with brain morphology and cognition, with deletion carriers being particularly affected. The pattern of results fits with known molecular functions of genes in the 15q11.2 BP1-BP2 region and suggests involvement of these genes in neuronal plasticity. These neurobiological effects likely contribute to the association of this CNV with neurodevelopmental disorders.

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.

Correction to: Diagnostic yield of a custom-designed multi-gene cancer panel in Irish patients with breast cancer

The originally published version of this article contained typesetting errors in Figs. 2 and 3 legends. The correct figure legends are presented here. The original article has been corrected.

Diagnostic yield of a custom-designed multi-gene cancer panel in Irish patients with breast cancer

BACKGROUND

Breast cancer is genetically heterogeneous, and parellel multi-gene sequencing is the most cost- and time-efficient manner to investigate breast cancer predisposition. Numerous multi-gene panels (MGPs) are commercially available, but many include genes with weak/unproven associaton with breast cancer, or with predisposition to cancer of other types. This study investigates the utility of a custom-designed multi-gene panel in an Irish cohort with breast cancer.

METHODS

A custom panel comprising 83 genes offered by 19 clinical "breast cancer predisposition" MGPs was designed and applied to germline DNA from 91 patients with breast cancer and 77 unaffected ethnicially matched controls. Variants were identified and classified using a custom pipeline.

RESULTS

Nineteen loss-of-function (LOF) and 334 missense variants were identified. After removing common and/or benign variants, 15 LOF and 30 missense variants were analysed. Variants in known breast cancer susceptibility genes were identified, including in BRCA1 and ATM in cases, and in NF1 and CHEK2 in controls. Most variants identified were in genes associated with predisposition to cancers other than breast cancer (BRIP1, RAD50, MUTYH, and mismatch repair genes), or in genes with unknown or unproven association with cancer.

CONCLUSION

Using multi-gene panels enables rapid, cost-effective identification of individuals with high-risk cancer predisposition syndromes. However, this approach also leads to an increased amount of uncertain results. Clinical management of individuals with particular genetic variants in the absence of a matching phenotype/family history is challenging. Further population and functional evidence is required to fully elucidate the clinical relevance of variants in genes of uncertain significance.

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.

The effect of breed and diet type on the global transcriptome of hepatic tissue in beef cattle divergent for feed efficiency

Background

Feed efficiency is an important economic and environmental trait in beef production, which can be measured in terms of residual feed intake (RFI). Cattle selected for low-RFI (feed efficient) have similar production levels but decreased feed intake, while also emitting less methane. RFI is difficult and expensive to measure and is not widely adopted in beef production systems. However, development of DNA-based biomarkers for RFI may facilitate its adoption in genomic-assisted breeding programmes. Cattle have been shown to re-rank in terms of RFI across diets and age, while also RFI varies by breed. Therefore, we used RNA-Seq technology to investigate the hepatic transcriptome of RFI-divergent Charolais (CH) and Holstein-Friesian (HF) steers across three dietary phases to identify genes and biological pathways associated with RFI regardless of diet or breed.

Results

Residual feed intake was measured during a high-concentrate phase, a zero-grazed grass phase and a final high-concentrate phase. In total, 322 and 33 differentially expressed genes (DEGs) were identified across all diets for CH and HF steers, respectively. Three genes, GADD45G, HP and MID1IP1, were differentially expressed in CH when both the high-concentrate zero-grazed grass diet were offered. Two canonical pathways were enriched across all diets for CH steers. These canonical pathways were related to immune function.

Conclusions

The absence of common differentially expressed genes across all dietary phases and breeds in this study supports previous reports of the re-ranking of animals in terms of RFI when offered differing diets over their lifetime. However, we have identified biological processes such as the immune response and lipid metabolism as potentially associated with RFI divergence emphasising the previously reported roles of these biological processes with respect to RFI.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5906-8) contains supplementary material, which is available to authorized users.

Genome-wide association study identifies 30 loci associated with bipolar disorder

Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study (GWAS) including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P < 1 × 10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (P < 5 × 10-8) in the discovery GWAS were not genome-wide significant in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis, 30 loci were genome-wide significant, including 20 newly identified loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene sets, including regulation of insulin secretion and endocannabinoid signaling. Bipolar I disorder is strongly genetically correlated with schizophrenia, driven by psychosis, whereas bipolar II disorder is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential biological mechanisms for bipolar disorder.

Targeted Sequencing of 10,198 Samples Confirms Abnormalities in Neuronal Activity and Implicates Voltage-Gated Sodium Channels in Schizophrenia Pathogenesis

BACKGROUND

Sequencing studies have pointed to the involvement in schizophrenia of rare coding variants in neuronally expressed genes, including activity-regulated cytoskeleton-associated protein (ARC) and N-methyl-D-aspartate receptor (NMDAR) complexes; however, larger samples are required to reveal novel genes and specific biological mechanisms.

METHODS

We sequenced 187 genes, selected for prior evidence of association with schizophrenia, in a new dataset of 5207 cases and 4991 controls. Included among these genes were members of ARC and NMDAR postsynaptic protein complexes, as well as voltage-gated sodium and calcium channels. We performed a rare variant meta-analysis with published sequencing data for a total of 11,319 cases, 15,854 controls, and 1136 trios.

RESULTS

While no individual gene was significantly associated with schizophrenia after genome-wide correction for multiple testing, we strengthen the evidence that rare exonic variants in the ARC (p = 4.0 × 10-4) and NMDAR (p = 1.7 × 10-5) synaptic complexes are risk factors for schizophrenia. In addition, we found that loss-of-function variants and missense variants at paralog-conserved sites were enriched in voltage-gated sodium channels, particularly the alpha subunits (p = 8.6 × 10-4).

CONCLUSIONS

In one of the largest sequencing studies of schizophrenia to date, we provide novel evidence that multiple voltage-gated sodium channels are involved in schizophrenia pathogenesis and confirm the involvement of ARC and NMDAR postsynaptic complexes.

Population-based identity-by-descent mapping combined with exome sequencing to detect rare risk variants for schizophrenia

Genome-wide association studies (GWASs) are highly effective at identifying common risk variants for schizophrenia. Rare risk variants are also important contributors to schizophrenia etiology but, with the exception of large copy number variants, are difficult to detect with GWAS. Exome and genome sequencing, which have accelerated the study of rare variants, are expensive so alternative methods are needed to aid detection of rare variants. Here we re-analyze an Irish schizophrenia GWAS dataset (n = 3,473) by performing identity-by-descent (IBD) mapping followed by exome sequencing of individuals identified as sharing risk haplotypes to search for rare risk variants in coding regions. We identified 45 rare haplotypes (>1 cM) that were significantly more common in cases than controls. By exome sequencing 105 haplotype carriers, we investigated these haplotypes for functional coding variants that could be tested for association in independent GWAS samples. We identified one rare missense variant in PCNT but did not find statistical support for an association with schizophrenia in a replication analysis. However, IBD mapping can prioritize both individual samples and genomic regions for follow-up analysis but genome rather than exome sequencing may be more effective at detecting risk variants on rare haplotypes.

Genes encoding SATB2-interacting proteins in adult cerebral cortex contribute to human cognitive ability

During CNS development, the nuclear protein SATB2 is expressed in superficial cortical layers and determines projection neuron identity. In the adult CNS, SATB2 is expressed in pyramidal neurons of all cortical layers and is a regulator of synaptic plasticity and long-term memory. Common variation in SATB2 locus confers risk of schizophrenia, whereas rare, de novo structural and single nucleotide variants cause severe intellectual disability and absent or limited speech. To characterize differences in SATB2 molecular function in developing vs adult neocortex, we isolated SATB2 protein interactomes at the two ontogenetic stages and identified multiple novel SATB2 interactors. SATB2 interactomes are highly enriched for proteins that stabilize de novo chromatin loops. The comparison between the neonatal and adult SATB2 protein complexes indicates a developmental shift in SATB2 molecular function, from transcriptional repression towards organization of chromosomal superstructure. Accordingly, gene sets regulated by SATB2 in the neocortex of neonatal and adult mice show limited overlap. Genes encoding SATB2 protein interactors were grouped for gene set analysis of human GWAS data. Common variants associated with human cognitive ability are enriched within the genes encoding adult but not neonatal SATB2 interactors. Our data support a shift in the function of SATB2 in cortex over lifetime and indicate that regulation of spatial chromatin architecture by the SATB2 interactome contributes to cognitive function in the general population.

Cortical Brain Abnormalities in 4474 Individuals With Schizophrenia and 5098 Control Subjects via the Enhancing Neuro Imaging Genetics Through Meta Analysis (ENIGMA) Consortium

BACKGROUND

The profile of cortical neuroanatomical abnormalities in schizophrenia is not fully understood, despite hundreds of published structural brain imaging studies. This study presents the first meta-analysis of cortical thickness and surface area abnormalities in schizophrenia conducted by the ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) Schizophrenia Working Group.

METHODS

The study included data from 4474 individuals with schizophrenia (mean age, 32.3 years; range, 11-78 years; 66% male) and 5098 healthy volunteers (mean age, 32.8 years; range, 10-87 years; 53% male) assessed with standardized methods at 39 centers worldwide.

RESULTS

Compared with healthy volunteers, individuals with schizophrenia have widespread thinner cortex (left/right hemisphere: Cohen's d = -0.530/-0.516) and smaller surface area (left/right hemisphere: Cohen's d = -0.251/-0.254), with the largest effect sizes for both in frontal and temporal lobe regions. Regional group differences in cortical thickness remained significant when statistically controlling for global cortical thickness, suggesting regional specificity. In contrast, effects for cortical surface area appear global. Case-control, negative, cortical thickness effect sizes were two to three times larger in individuals receiving antipsychotic medication relative to unmedicated individuals. Negative correlations between age and bilateral temporal pole thickness were stronger in individuals with schizophrenia than in healthy volunteers. Regional cortical thickness showed significant negative correlations with normalized medication dose, symptom severity, and duration of illness and positive correlations with age at onset.

CONCLUSIONS

The findings indicate that the ENIGMA meta-analysis approach can achieve robust findings in clinical neuroscience studies; also, medication effects should be taken into account in future genetic association studies of cortical thickness in schizophrenia.

GWAS and eQTL analysis identifies a SNP associated with both residual feed intake and GFRA2 expression in beef cattle

Residual feed intake (RFI), a measure of feed efficiency, is an important economic and environmental trait in beef production. Selection of low RFI (feed efficient) cattle could maintain levels of production, while decreasing feed costs and methane emissions. However, RFI is a difficult and expensive trait to measure. Identification of single nucleotide polymorphisms (SNPs) associated with RFI may enable rapid, cost effective genomic selection of feed efficient cattle. Genome-wide association studies (GWAS) were conducted in multiple breeds followed by meta-analysis to identify genetic variants associated with RFI and component traits (average daily gain (ADG) and feed intake (FI)) in Irish beef cattle (n = 1492). Expression quantitative trait loci (eQTL) analysis was conducted to identify functional effects of GWAS-identified variants. Twenty-four SNPs were associated (P < 5 × 10^-5) with RFI, ADG or FI. The variant rs43555985 exhibited strongest association for RFI (P = 8.28E-06). An eQTL was identified between this variant and GFRA2 (P = 0.0038) where the allele negatively correlated with RFI was associated with increased GFRA2 expression in liver. GFRA2 influences basal metabolic rates, suggesting a mechanism by which genetic variation may contribute to RFI. This study identified SNPs that may be useful both for genomic selection of RFI and for understanding the biology of feed efficiency.

Genes regulated by SATB2 during neurodevelopment contribute to schizophrenia and educational attainment

SATB2 is associated with schizophrenia and is an important transcription factor regulating neocortical organization and circuitry. Rare mutations in SATB2 cause a syndrome that includes developmental delay, and mouse studies identify an important role for SATB2 in learning and memory. Interacting partners BCL11B and GATAD2A are also schizophrenia risk genes indicating that other genes interacting with or are regulated by SATB2 are making a contribution to schizophrenia and cognition. We used data from Satb2 mouse models to generate three gene-sets that contain genes either functionally related to SATB2 or targeted by SATB2 at different stages of development. Each was tested for enrichment using the largest available genome-wide association studies (GWAS) datasets for schizophrenia and educational attainment (EA) and enrichment analysis was also performed for schizophrenia and other neurodevelopmental disorders using data from rare variant sequencing studies. These SATB2 gene-sets were enriched for genes containing common variants associated with schizophrenia and EA, and were enriched for genes containing rare variants reported in studies of schizophrenia, autism and intellectual disability. In the developing cortex, genes targeted by SATB2 based on ChIP-seq data, and functionally affected when SATB2 is not expressed based on differential expression analysis using RNA-seq data, show strong enrichment for genes associated with EA. For genes expressed in the hippocampus or at the synapse, those targeted by SATB2 are more strongly enriched for genes associated EA than gene-sets not targeted by SATB2. This study demonstrates that single gene findings from GWAS can provide important insights to pathobiological processes. In this case we find evidence that genes influenced by SATB2 and involved in synaptic transmission, axon guidance and formation of the corpus callosum are contributing to schizophrenia and cognition.

Schizophrenia is a complex disorder caused by many genes. Using new gene discoveries to understand pathobiology is a foundation for development of new treatments. Current drugs for schizophrenia are only partially effective and do not treat cognitive deficits, which are key factors for explaining disability, leading to unemployment, homelessness and social isolation. Genome-wide association studies (GWAS) of schizophrenia have been effective at identifying individual SNPs and genes that contribute to risk but have struggled to immediately uncover the bigger picture of the underlying biology of the disorder. Here we take an individual gene identified in a schizophrenia GWAS called SATB2, which on its own is a very important regulator of brain development. We use functional genomics data from mouse studies to identify sets of others genes that are influenced by SATB2 during development. We show that these gene sets are enriched for common variants associated with schizophrenia and educational attainment (used as a proxy for cognition), and for rare variants that increase risk of various neurodevelopmental disorders. This study provides evidence that the molecular mechanisms that underpin schizophrenia and cognitive function include disruption of biological processes influenced by SATB2 as the brain is being organized and wired during development.

Genetically predicted complement component 4A expression: effects on memory function and middle temporal lobe activation

BACKGROUND

The longstanding association between the major histocompatibility complex (MHC) locus and schizophrenia (SZ) risk has recently been accounted for, partially, by structural variation at the complement component 4 (C4) gene. This structural variation generates varying levels of C4 RNA expression, and genetic information from the MHC region can now be used to predict C4 RNA expression in the brain. Increased predicted C4A RNA expression is associated with the risk of SZ, and C4 is reported to influence synaptic pruning in animal models.

METHODS

Based on our previous studies associating MHC SZ risk variants with poorer memory performance, we tested whether increased predicted C4A RNA expression was associated with reduced memory function in a large (n = 1238) dataset of psychosis cases and healthy participants, and with altered task-dependent cortical activation in a subset of these samples.

RESULTS

We observed that increased predicted C4A RNA expression predicted poorer performance on measures of memory recall (p = 0.016, corrected). Furthermore, in healthy participants, we found that increased predicted C4A RNA expression was associated with a pattern of reduced cortical activity in middle temporal cortex during a measure of visual processing (p < 0.05, corrected).

CONCLUSIONS

These data suggest that the effects of C4 on cognition were observable at both a cortical and behavioural level, and may represent one mechanism by which illness risk is mediated. As such, deficits in learning and memory may represent a therapeutic target for new molecular developments aimed at altering C4's developmental role.

Estimation of Genetic Correlation via Linkage Disequilibrium Score Regression and Genomic Restricted Maximum Likelihood

Genetic correlation is a key population parameter that describes the shared genetic architecture of complex traits and diseases. It can be estimated by current state-of-art methods, i.e., linkage disequilibrium score regression (LDSC) and genomic restricted maximum likelihood (GREML). The massively reduced computing burden of LDSC compared to GREML makes it an attractive tool, although the accuracy (i.e., magnitude of standard errors) of LDSC estimates has not been thoroughly studied. In simulation, we show that the accuracy of GREML is generally higher than that of LDSC. When there is genetic heterogeneity between the actual sample and reference data from which LD scores are estimated, the accuracy of LDSC decreases further. In real data analyses estimating the genetic correlation between schizophrenia (SCZ) and body mass index, we show that GREML estimates based on ∼150,000 individuals give a higher accuracy than LDSC estimates based on ∼400,000 individuals (from combined meta-data). A GREML genomic partitioning analysis reveals that the genetic correlation between SCZ and height is significantly negative for regulatory regions, which whole genome or LDSC approach has less power to detect. We conclude that LDSC estimates should be carefully interpreted as there can be uncertainty about homogeneity among combined meta-datasets. We suggest that any interesting findings from massive LDSC analysis for a large number of complex traits should be followed up, where possible, with more detailed analyses with GREML methods, even if sample sizes are lesser.

The Genetics of Endophenotypes of Neurofunction to Understand Schizophrenia (GENUS) consortium: A collaborative cognitive and neuroimaging genetics project

BACKGROUND

Schizophrenia has a large genetic component, and the pathways from genes to illness manifestation are beginning to be identified. The Genetics of Endophenotypes of Neurofunction to Understand Schizophrenia (GENUS) Consortium aims to clarify the role of genetic variation in brain abnormalities underlying schizophrenia. This article describes the GENUS Consortium sample collection.

METHODS

We identified existing samples collected for schizophrenia studies consisting of patients, controls, and/or individuals at familial high-risk (FHR) for schizophrenia. Samples had single nucleotide polymorphism (SNP) array data or genomic DNA, clinical and demographic data, and neuropsychological and/or brain magnetic resonance imaging (MRI) data. Data were subjected to quality control procedures at a central site.

RESULTS

Sixteen research groups contributed data from 5199 psychosis patients, 4877 controls, and 725 FHR individuals. All participants have relevant demographic data and all patients have relevant clinical data. The sex ratio is 56.5% male and 43.5% female. Significant differences exist between diagnostic groups for premorbid and current IQ (both p<1×10^-10). Data from a diversity of neuropsychological tests are available for 92% of participants, and 30% have structural MRI scans (half also have diffusion-weighted MRI scans). SNP data are available for 76% of participants. The ancestry composition is 70% European, 20% East Asian, 7% African, and 3% other.

CONCLUSIONS

The Consortium is investigating the genetic contribution to brain phenotypes in a schizophrenia sample collection of >10,000 participants. The breadth of data across clinical, genetic, neuropsychological, and MRI modalities provides an important opportunity for elucidating the genetic basis of neural processes underlying schizophrenia.