An evolutionary timeline of the oxytocin signaling pathway

Oxytocin is a neuropeptide associated with both psychological and somatic processes like parturition and social bonding. Although oxytocin homologs have been identified in many species, the evolutionary timeline of the entire oxytocin signaling gene pathway has yet to be described. Using protein sequence similarity searches, microsynteny, and phylostratigraphy, we assigned the genes supporting the oxytocin pathway to different phylostrata based on when we found they likely arose in evolution. We show that the majority (64%) of genes in the pathway are 'modern'. Most of the modern genes evolved around the emergence of vertebrates or jawed vertebrates (540 - 530 million years ago, 'mya'), including OXTR, OXT and CD38. Of those, 45% were under positive selection at some point during vertebrate evolution. We also found that 18% of the genes in the oxytocin pathway are 'ancient', meaning their emergence dates back to cellular organisms and opisthokonta (3500-1100 mya). The remaining genes (18%) that evolved after ancient and before modern genes were classified as 'medium-aged'. Functional analyses revealed that, in humans, medium-aged oxytocin pathway genes are highly expressed in contractile organs, while modern genes in the oxytocin pathway are primarily expressed in the brain and muscle tissue.

Polygenic risk for schizophrenia and bipolar disorder in relation to cardiovascular biomarkers

Individuals with schizophrenia and bipolar disorder are at an increased risk of cardiovascular disease (CVD), and a range of biomarkers related to CVD risk have been found to be abnormal in these patients. Common genetic factors are a putative underlying mechanism, alongside lifestyle factors and antipsychotic medication. However, the extent to which the altered CVD biomarkers are related to genetic factors involved in schizophrenia and bipolar disorder is unknown. In a sample including 699 patients with schizophrenia, 391 with bipolar disorder, and 822 healthy controls, we evaluated 8 CVD risk biomarkers, including BMI, and fasting plasma levels of CVD biomarkers from a subsample. Polygenic risk scores (PGRS) were obtained from genome-wide associations studies (GWAS) of schizophrenia and bipolar disorder from the Psychiatric Genomics Consortium. The CVD biomarkers were used as outcome variables in linear regression models including schizophrenia and bipolar disorder PGRS as predictors, age, sex, diagnostic category, batch and 10 principal components as covariates, controlling for multiple testing by Bonferroni correction for the number of independent tests. Bipolar disorder PGRS was significantly (p = 0.03) negatively associated with BMI after multiple testing correction, and schizophrenia PGRS was nominally negatively associated with BMI. There were no other significant associations between bipolar or schizophrenia PGRS, and other investigated CVD biomarkers. Despite a range of abnormal CVD risk biomarkers in psychotic disorders, we only found a significant negative association between bipolar disorder PGRS and BMI. This has previously been shown for schizophrenia PGRS and BMI, and warrants further exploration.

Boosting Schizophrenia Genetics by Utilizing Genetic Overlap With Brain Morphology

BACKGROUND

Schizophrenia is a complex polygenic disorder with subtle, distributed abnormalities in brain morphology. There are indications of shared genetic architecture between schizophrenia and brain measures despite low genetic correlations. Through the use of analytical methods that allow for mixed directions of effects, this overlap may be leveraged to improve our understanding of underlying mechanisms of schizophrenia and enrich polygenic risk prediction outcome.

METHODS

We ran a multivariate genome-wide analysis of 175 brain morphology measures using data from 33,735 participants of the UK Biobank and analyzed the results in a conditional false discovery rate together with schizophrenia genome-wide association study summary statistics of the Psychiatric Genomics Consortium (PGC) Wave 3. We subsequently created a pleiotropy-enriched polygenic score based on the loci identified through the conditional false discovery rate approach and used this to predict schizophrenia in a nonoverlapping sample of 743 individuals with schizophrenia and 1074 healthy controls.

RESULTS

We found that 20% of the loci and 50% of the genes significantly associated with schizophrenia were also associated with brain morphology. The conditional false discovery rate analysis identified 428 loci, including 267 novel loci, significantly associated with brain-linked schizophrenia risk, with functional annotation indicating high relevance for brain tissue. The pleiotropy-enriched polygenic score explained more variance in liability than conventional polygenic scores across several scenarios.

CONCLUSIONS

Our results indicate strong genetic overlap between schizophrenia and brain morphology with mixed directions of effect. The results also illustrate the potential of exploiting polygenetic overlap between brain morphology and mental disorders to boost discovery of brain tissue-specific genetic variants and its use in polygenic risk frameworks.

Neural correlates of user learning during long-term BCI training for the Cybathlon competition

Background

Brain-computer interfaces (BCIs) are systems capable of translating human brain patterns, measured through electroencephalography (EEG), into commands for an external device. Despite the great advances in machine learning solutions to enhance the performance of BCI decoders, the translational impact of this technology remains elusive. The reliability of BCIs is often unsatisfactory for end-users, limiting their application outside a laboratory environment.

Methods

We present the analysis on the data acquired from an end-user during the preparation for two Cybathlon competitions, where our pilot won the gold medal twice in a row. These data are of particular interest given the mutual learning approach adopted during the longitudinal training phase (8 months), the long training break in between the two events (1 year) and the demanding evaluation scenario. A multifaceted perspective on long-term user learning is proposed: we enriched the information gathered through conventional metrics (e.g., accuracy, application performances) by investigating novel neural correlates of learning in different neural domains.

Results

First, we showed that by focusing the training on user learning, the pilot was capable of significantly improving his performance over time even with infrequent decoder re-calibrations. Second, we revealed that the analysis of the within-class modifications of the pilot’s neural patterns in the Riemannian domain is more effective in tracking the acquisition and the stabilization of BCI skills, especially after the 1-year break. These results further confirmed the key role of mutual learning in the acquisition of BCI skills, and particularly highlighted the importance of user learning as a key to enhance BCI reliability.

Conclusion

We firmly believe that our work may open new perspectives and fuel discussions in the BCI field to shift the focus of future research: not only to the machine learning of the decoder, but also in investigating novel training procedures to boost the user learning and the stability of the BCI skills in the long-term. To this end, the analyses and the metrics proposed could be used to monitor the user learning during training and provide a marker guiding the decoder re-calibration to maximize the mutual adaptation of the user to the BCI system.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12984-022-01047-x.

Dose-dependent transcriptional effects of lithium and adverse effect burden in a psychiatric cohort

Lithium is the first-line treatment for bipolar disorder (BD), but there is a large variation in response rate and adverse effects. Although the molecular effects of lithium have been studied extensively, the specific mechanisms of action remain unclear. In particular, the molecular changes underlying lithium adverse effects are little known. Multiple linear regression analyses of lithium serum concentrations and global gene expression levels in whole blood were carried out using a large case-control sample (n = 1450). Self-reported adverse effects of lithium were assessed with the "Udvalg for Kliniske Undersøgelser" (UKU) adverse effect rating scale, and regression analysis was used to identify significant associations between lithium-related genes and six of the most common adverse effects. Serum concentrations of lithium were significantly associated with the expression levels of 52 genes (FDR < 0.01), largely replicating previous results. We found 32 up-regulated genes and 20 down-regulated genes in lithium users compared to non-users. The down-regulated gene set was enriched for several processes related to the translational machinery. Two adverse effects were significantly associated (p < 0.01) with three or more lithium-associated genes: tremor (FAM13A-AS1, FAR2, ITGAX, RWDD1, and STARD10) and xerostomia (ANKRD13A, FAR2, RPS8, and RWDD1). The adverse effect association with the largest effect was between CAMK1D expression and nausea/vomiting. These results suggest putative transcriptional mechanisms that may predict lithium adverse effects, and could thus have a large potential for informing clinical practice.

Phenotype‐informed polygenic risk scores are associated with worse outcome in individuals at risk of Alzheimer's disease

Introduction

Patients with predementia Alzheimer's disease (AD) and at‐risk subjects are targets for promising disease‐modifying treatments, and improved polygenic risk scores (PRSs) could improve early‐stage case selection.

Methods

Phenotype‐informed PRSs were developed by selecting AD‐associated variants conditional on relevant inflammatory or cardiovascular traits. The primary outcome was longitudinal changes in measures of AD pathology, namely development of pathological amyloid deposition, medial temporal lobe atrophy, and cognitive decline in a prospective cohort study including 394 adults without AD dementia.

Results

High‐risk groups defined by phenotype‐informed AD PRSs had significantly steeper volume decline in medial temporal cortices, and the high‐risk group defined by the cardiovascular‐informed AD PRS had significantly increased hazard ratio of pathological amyloid deposition, compared to low‐risk groups.

Discussion

AD PRSs informed by inflammatory disorders or cardiovascular risk factors and diseases are associated with development of AD pathology markers and may improve identification of subjects at risk for progression of AD.

Oxytocin-pathway polygenic scores for severe mental disorders and metabolic phenotypes in the UK Biobank

Oxytocin is a neuromodulator and hormone that is typically associated with social cognition and behavior. In light of its purported effects on social cognition and behavior, research has investigated its potential as a treatment for psychiatric illnesses characterized by social dysfunction, such as schizophrenia and bipolar disorder. While the results of these trials have been mixed, more recent evidence suggests that the oxytocin system is also linked with cardiometabolic conditions for which individuals with severe mental disorders are at a higher risk for developing. To investigate whether the oxytocin system has a pleiotropic effect on the etiology of severe mental illness and cardiometabolic conditions, we explored oxytocin's role in the shared genetic liability of schizophrenia, bipolar disorder, type-2 diabetes, and several phenotypes linked with cardiovascular disease and type 2 diabetes risk using a polygenic pathway-specific approach. Analysis of a large sample with about 480,000 individuals (UK Biobank) revealed statistically significant associations across the range of phenotypes analyzed. By comparing these effects to those of polygenic scores calculated from 100 random gene sets, we also demonstrated the specificity of many of these significant results. Altogether, our results suggest that the shared effect of oxytocin-system dysfunction could help partially explain the co-occurrence of social and cardiometabolic dysfunction in severe mental illnesses.

Identification of genetic overlap and novel risk loci for attention-deficit/hyperactivity disorder and bipolar disorder

Differential diagnosis between childhood onset attention-deficit/hyperactivity disorder (ADHD) and bipolar disorder (BD) remains a challenge, mainly due to overlapping symptoms and high rates of comorbidity. Despite this, genetic correlation reported for these disorders is low and non-significant. Here we aimed to better characterize the genetic architecture of these disorders utilizing recent large genome-wide association studies (GWAS). We analyzed independent GWAS summary statistics for ADHD (19,099 cases and 34,194 controls) and BD (20,352 cases and 31,358 controls) applying the conditional/conjunctional false discovery rate (condFDR/conjFDR) statistical framework that increases the power to detect novel phenotype-specific and shared loci by leveraging the combined power of two GWAS. We observed cross-trait polygenic enrichment for ADHD conditioned on associations with BD, and vice versa. Leveraging this enrichment, we identified 19 novel ADHD risk loci and 40 novel BD risk loci at condFDR <0.05. Further, we identified five loci jointly associated with ADHD and BD (conjFDR < 0.05). Interestingly, these five loci show concordant directions of effect for ADHD and BD. These results highlight a shared underlying genetic risk for ADHD and BD which may help to explain the high comorbidity rates and difficulties in differentiating between ADHD and BD in the clinic. Improving our understanding of the underlying genetic architecture of these disorders may aid in the development of novel stratification tools to help reduce these diagnostic difficulties.

Genetic control of variability in subcortical and intracranial volumes

Sensitivity to external demands is essential for adaptation to dynamic environments, but comes at the cost of increased risk of adverse outcomes when facing poor environmental conditions. Here, we apply a novel methodology to perform genome-wide association analysis of mean and variance in ten key brain features (accumbens, amygdala, caudate, hippocampus, pallidum, putamen, thalamus, intracranial volume, cortical surface area, and cortical thickness), integrating genetic and neuroanatomical data from a large lifespan sample (n = 25,575 individuals; 8-89 years, mean age 51.9 years). We identify genetic loci associated with phenotypic variability in thalamus volume and cortical thickness. The variance-controlling loci involved genes with a documented role in brain and mental health and were not associated with the mean anatomical volumes. This proof-of-principle of the hypothesis of a genetic regulation of brain volume variability contributes to establishing the genetic basis of phenotypic variance (i.e., heritability), allows identifying different degrees of brain robustness across individuals, and opens new research avenues in the search for mechanisms controlling brain and mental health.

Dissecting the shared genetic basis of migraine and mental disorders using novel statistical tools

Migraine is three times more prevalent in people with bipolar disorder or depression. The relationship between schizophrenia and migraine is less certain although glutamatergic and serotonergic neurotransmission are implicated in both. A shared genetic basis to migraine and mental disorders has been suggested but previous studies have reported weak or non-significant genetic correlations and five shared risk loci. Using the largest samples to date and novel statistical tools, we aimed to determine the extent to which migraine’s polygenic architecture overlaps with bipolar disorder, depression and schizophrenia beyond genetic correlation, and to identify shared genetic loci.

Summary statistics from genome-wide association studies were acquired from large-scale consortia for migraine (n cases = 59 674; n controls = 316 078), bipolar disorder (n cases = 20 352; n controls = 31 358), depression (n cases = 170 756; n controls = 328 443) and schizophrenia (n cases = 40 675, n controls = 64 643). We applied the bivariate causal mixture model to estimate the number of disorder-influencing variants shared between migraine and each mental disorder, and the conditional/conjunctional false discovery rate method to identify shared loci. Loci were functionally characterized to provide biological insights.

Univariate MiXeR analysis revealed that migraine was substantially less polygenic (2.8 K disorder-influencing variants) compared to mental disorders (8100–12 300 disorder-influencing variants). Bivariate analysis estimated that 800 (SD = 300), 2100 (SD = 100) and 2300 (SD = 300) variants were shared between bipolar disorder, depression and schizophrenia, respectively. There was also extensive overlap with intelligence (1800, SD = 300) and educational attainment (2100, SD = 300) but not height (1000, SD = 100). We next identified 14 loci jointly associated with migraine and depression and 36 loci jointly associated with migraine and schizophrenia, with evidence of consistent genetic effects in independent samples. No loci were associated with migraine and bipolar disorder. Functional annotation mapped 37 and 298 genes to migraine and each of depression and schizophrenia, respectively, including several novel putative migraine genes such as L3MBTL2, CACNB2 and SLC9B1. Gene-set analysis identified several putative gene sets enriched with mapped genes including transmembrane transport in migraine and schizophrenia.

Most migraine-influencing variants were predicted to influence depression and schizophrenia, although a minority of mental disorder-influencing variants were shared with migraine due to the difference in polygenicity. Similar overlap with other brain-related phenotypes suggests this represents a pool of ‘pleiotropic’ variants that influence vulnerability to diverse brain-related disorders and traits. We also identified specific loci shared between migraine and each of depression and schizophrenia, implicating shared molecular mechanisms and highlighting candidate migraine genes for experimental validation.

Phenotype-specific differences in polygenicity and effect size distribution across functional annotation categories revealed by AI-MiXeR

Motivation

Determining the relative contributions of functional genetic categories is fundamental to understanding the genetic etiology of complex human traits and diseases. Here, we present Annotation Informed-MiXeR, a likelihood-based method for estimating the number of variants influencing a phenotype and their effect sizes across different functional annotation categories of the genome using summary statistics from genome-wide association studies.

Results

Extensive simulations demonstrate that the model is valid for a broad range of genetic architectures. The model suggests that complex human phenotypes substantially differ in the number of causal variants, their localization in the genome and their effect sizes. Specifically, the exons of protein-coding genes harbor more than 90% of variants influencing type 2 diabetes and inflammatory bowel disease, making them good candidates for whole-exome studies. In contrast, <10% of the causal variants for schizophrenia, bipolar disorder and attention-deficit/hyperactivity disorder are located in protein-coding exons, indicating a more substantial role of regulatory mechanisms in the pathogenesis of these disorders.

Availability and implementation

The software is available at: https://github.com/precimed/mixer.

Supplementary information

Supplementary data are available at Bioinformatics online.

Optimising classification of proximal arm strength impairment in wheelchair rugby: A proof of concept study

This study examined the relationship between proximal arm strength and mobility performance in wheelchair rugby (WR) athletes and examined whether a valid structure for classifying proximal arm strength impairment could be determined. Fifty-seven trained WR athletes with strength impaired arms and no trunk function performed six upper body isometric strength tests and three 10 m sprints in their rugby wheelchair. All strength measures correlated with 2 m and 10 m sprint times (r ≥ -0.43; p ≤ 0.0005) and were entered into k-means cluster analyses with 4-clusters (to mirror the current International Wheelchair Rugby Federation [IWRF] system) and 3-clusters. The 3-cluster structure provided a more valid structure than both the 4-cluster and existing IWRF system, as evidenced by clearer differences in strength (Effect sizes [ES] ≥ 1.0) and performance (ES ≥ 1.1) between adjacent clusters and stronger mean silhouette coefficient (0.64). Subsequently, the 3-cluster structure for classifying proximal arm strength impairment would result in less overlap between athletes from adjacent classes and reduce the likelihood of athletes being disadvantaged due to their impairment. This study demonstrated that the current battery of isometric strength tests and cluster analyses could facilitate the evidence-based development of classifying proximal arm strength impairment in WR.

In vivo hippocampal subfield volumes in bipolar disorder—A mega‐analysis from The Enhancing Neuro Imaging Genetics through Meta‐Analysis Bipolar Disorder Working Group

The hippocampus consists of anatomically and functionally distinct subfields that may be differentially involved in the pathophysiology of bipolar disorder (BD). Here we, the Enhancing NeuroImaging Genetics through Meta‐Analysis Bipolar Disorder workinggroup, study hippocampal subfield volumetry in BD. T1‐weighted magnetic resonance imaging scans from 4,698 individuals (BD = 1,472, healthy controls [HC] = 3,226) from 23 sites worldwide were processed with FreeSurfer. We used linear mixed‐effects models and mega‐analysis to investigate differences in hippocampal subfield volumes between BD and HC, followed by analyses of clinical characteristics and medication use. BD showed significantly smaller volumes of the whole hippocampus (Cohen's d = −0.20), cornu ammonis (CA)1 (d = −0.18), CA2/3 (d = −0.11), CA4 (d = −0.19), molecular layer (d = −0.21), granule cell layer of dentate gyrus (d = −0.21), hippocampal tail (d = −0.10), subiculum (d = −0.15), presubiculum (d = −0.18), and hippocampal amygdala transition area (d = −0.17) compared to HC. Lithium users did not show volume differences compared to HC, while non‐users did. Antipsychotics or antiepileptic use was associated with smaller volumes. In this largest study of hippocampal subfields in BD to date, we show widespread reductions in nine of 12 subfields studied. The associations were modulated by medication use and specifically the lack of differences between lithium users and HC supports a possible protective role of lithium in BD.

Identification of Genetic Loci Shared Between Attention-Deficit/Hyperactivity Disorder, Intelligence, and Educational Attainment

BACKGROUND

Attention-deficit/hyperactivity disorder (ADHD) is a neurodevelopmental disorder that is consistently associated with lower levels of educational attainment. A recent large genome-wide association study identified common gene variants associated with ADHD, but most of the genetic architecture remains unknown.

METHODS

We analyzed independent genome-wide association study summary statistics for ADHD (19,099 cases and 34,194 controls), educational attainment (N = 842,499), and general intelligence (N = 269,867) using a conditional/conjunctional false discovery rate (FDR) statistical framework that increases power of discovery by conditioning the FDR on overlapping associations. The genetic variants identified were characterized in terms of function, expression, and biological processes.

RESULTS

We identified 58 linkage disequilibrium-independent ADHD-associated loci (conditional FDR < 0.01), of which 30 were shared between ADHD and educational attainment or general intelligence (conjunctional FDR < 0.01) and 46 were novel risk loci for ADHD.

CONCLUSIONS

These results expand on previous genetic and epidemiological studies and support the hypothesis of a shared genetic basis between these phenotypes. Although the clinical utility of the identified loci remains to be determined, they can be used as resources to guide future studies aiming to disentangle the complex etiologies of ADHD, educational attainment, and general intelligence.

Shared Genetic Loci Between Body Mass Index and Major Psychiatric Disorders: A Genome-wide Association Study

Importance

People with major psychiatric disorders (MPDs) have a 10- to 20-year shorter life span than the rest of the population, and this difference is mainly due to comorbid cardiovascular diseases. Genome-wide association studies have identified common variants involved in schizophrenia (SCZ), bipolar disorder (BIP), and major depression (MD) and body mass index (BMI), a key cardiometabolic risk factor. However, genetic variants jointly influencing MPD and BMI remain largely unknown.

Objective

To assess the extent of the overlap between the genetic architectures of MPDs and BMI and identify genetic loci shared between them.

Design, Setting, and Participants

Using a conditional false discovery rate statistical framework, independent genome-wide association study data on individuals with SCZ (n = 82 315), BIP (n = 51 710), MD (n = 480 359), and BMI (n = 795 640) were analyzed. The UK Biobank cohort (n = 29 740) was excluded from the MD data set to avoid sample overlap. Data were collected from August 2017 to May 2018, and analysis began July 2018.

Main Outcomes and Measures

The primary outcomes were a list of genetic loci shared between BMI and MPDs and their functional pathways.

Results

Genome-wide association study data from 1 380 284 participants were analyzed, and the genetic correlation between BMI and MPDs varied (SCZ: r for genetic = -0.11, P = 2.1 × 10-10; BIP: r for genetic = -0.06, P = .0103; MD: r for genetic = 0.12, P = 6.7 × 10-10). Overall, 63, 17, and 32 loci shared between BMI and SCZ, BIP, and MD, respectively, were analyzed at conjunctional false discovery rate less than 0.01. Of the shared loci, 34% (73 of 213) in SCZ, 52% (36 of 69) in BIP, and 57% (56 of 99) in MD had risk alleles associated with higher BMI (conjunctional false discovery rate <0.05), while the rest had opposite directions of associations. Functional analyses indicated that the overlapping loci are involved in several pathways including neurodevelopment, neurotransmitter signaling, and intracellular processes, and the loci with concordant and opposite association directions pointed mostly to different pathways.

Conclusions and Relevance

In this genome-wide association study, extensive polygenic overlap between BMI and SCZ, BIP, and MD were found, and 111 shared genetic loci were identified, implicating novel functional mechanisms. There was mixture of association directions in SCZ and BMI, albeit with a preponderance of discordant ones.

Associations between schizophrenia polygenic risk and apathy in schizophrenia spectrum disorders and healthy controls

OBJECTIVE

Apathy is a central predictor of a poor functional outcome in schizophrenia. Schizophrenia polygenic risk scores (PRSs) are used to detect genetic associations to key clinical phenotypes in schizophrenia. We explored the associations between schizophrenia PRS and apathy levels in schizophrenia spectrum disorders (n = 281) and matched healthy controls (n = 298), and further how schizophrenia PRS contributed in predicting apathy when added to premorbid and clinical factors in the patient sample.

METHOD

Schizophrenia PRSs were computed for each participant. Apathy was assessed with the Apathy Evaluation Scale. Bivariate correlation analyses were used to investigate associations between schizophrenia PRS and apathy, and between apathy and premorbid and clinical factors. Multiple hierarchical regression analyses were employed to evaluate the contributions of clinical variables and schizophrenia PRS to apathy levels.

RESULTS

We found no significant associations between schizophrenia PRS and apathy in patients and healthy controls. Several premorbid and clinical characteristics significantly predicted apathy in patients, but schizophrenia PRS did not.

CONCLUSION

Since the PRSs are based on common genetic variants, our results do not preclude associations to other types of genetic factors. The results could also indicate that environmentally based biological or psychological factors contribute to apathy levels in schizophrenia.

Kinematic bidimensional analysis of the propulsion technique in wheelchair rugby athletes

Wheelchair rugby is a sport ideated for individuals with cervical spinal cord injury (CSCI) which is extremely important for maintaining their neuromuscular abilities and improving their social and psychological wellbeing. However, due to the frequent changes in direction and speed it considerably stresses the players’ upper limbs. 13 athletes have undergone two sports-related tests on an inertial drum bench and several kinematic parameters have been registered. Most athletes use a semi-circular pattern which is considered protective for the upper limb. With increasing speed, range of motion (ROM) increases. Release angles increment and contact angles reduce, displacing the push angle forward to increase speed. Instead, the more anterior late push angle used to increase velocity is a factor which further loads the shoulder joint. However, other factors affecting propulsion technique, such as posture and wheelchair set up should be studied to further reduce loading on the upper limb.

Childhood maltreatment and polygenic risk in bipolar disorders

BACKGROUND

Childhood maltreatment is a well-known risk factor for developing a more severe and complex form of bipolar disorders (BD). However, knowledge is scarce about the interactions between childhood maltreatment and underlying genetic vulnerability on the clinical expression of BD.

METHOD

We assigned a BD-polygenic risk score (BD-PRS), calculated from the Psychiatric Genomics Consortium, to each individual in a sample of 402 cases with BD. The lifetime clinical expression of BD was characterized using structured interviews and patients completed the Childhood Trauma Questionnaire (CTQ) to assess the severity of childhood maltreatment.

RESULTS

Cases who reported more severe childhood maltreatment had a lower BD-PRS (rho = -0.12, P = .01), especially when considering emotional abuse (rho = -0.16, P = .001). An interaction between BD-PRS and childhood maltreatment was observed for the risk of rapid cycling (P = .01). No further interactions between BD-PRS and childhood maltreatment were observed for other clinical characteristics (age at onset, suicide attempts, number of mood episodes, mixed features, substance use disorders and psychotic symptoms).

CONCLUSION

Our study is the first to show that less genetic risk may be needed to develop a more unstable form of BD when exposed to childhood maltreatment. Our study supports childhood trauma as an independent risk factor for BD.

Author Correction: Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer's disease risk

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

Brain scans from 21,297 individuals reveal the genetic architecture of hippocampal subfield volumes

The hippocampus is a heterogeneous structure, comprising histologically distinguishable subfields. These subfields are differentially involved in memory consolidation, spatial navigation and pattern separation, complex functions often impaired in individuals with brain disorders characterized by reduced hippocampal volume, including Alzheimer's disease (AD) and schizophrenia. Given the structural and functional heterogeneity of the hippocampal formation, we sought to characterize the subfields' genetic architecture. T1-weighted brain scans (n = 21,297, 16 cohorts) were processed with the hippocampal subfields algorithm in FreeSurfer v6.0. We ran a genome-wide association analysis on each subfield, co-varying for whole hippocampal volume. We further calculated the single-nucleotide polymorphism (SNP)-based heritability of 12 subfields, as well as their genetic correlation with each other, with other structural brain features and with AD and schizophrenia. All outcome measures were corrected for age, sex and intracranial volume. We found 15 unique genome-wide significant loci across six subfields, of which eight had not been previously linked to the hippocampus. Top SNPs were mapped to genes associated with neuronal differentiation, locomotor behaviour, schizophrenia and AD. The volumes of all the subfields were estimated to be heritable (h2 from 0.14 to 0.27, all p < 1 × 10^-16) and clustered together based on their genetic correlations compared with other structural brain features. There was also evidence of genetic overlap of subicular subfield volumes with schizophrenia. We conclude that hippocampal subfields have partly distinct genetic determinants associated with specific biological processes and traits. Taking into account this specificity may increase our understanding of hippocampal neurobiology and associated pathologies.

Examining the association between genetic liability for schizophrenia and psychotic symptoms in Alzheimer's disease

Psychosis (delusions or hallucinations) in Alzheimer's disease (AD + P) occurs in up to 50% of individuals and is associated with significantly worse clinical outcomes. Atypical antipsychotics, first developed for schizophrenia, are commonly used in AD + P, suggesting shared mechanisms. Despite this implication, little empirical research has been conducted to examine whether there are mechanistic similarities between AD + P and schizophrenia. In this study, we tested whether polygenic risk score (PRS) for schizophrenia was associated with AD + P. Schizophrenia PRS was calculated using Psychiatric Genomics Consortium data at ten GWAS p value thresholds (PT) in 3111 AD cases from 11 cohort studies characterized for psychosis using validated, standardized tools. Association between PRS and AD + P status was tested by logistic regression in each cohort individually and the results meta-analyzed. The schizophrenia PRS was associated with AD + P at an optimum PT of 0.01. The strongest association was for delusions where a one standard deviation increase in PRS was associated with a 1.18-fold increased risk (95% CI: 1.06-1.3; p = 0.001). These new findings point towards psychosis in AD-and particularly delusions-sharing some genetic liability with schizophrenia and support a transdiagnostic view of psychotic symptoms across the lifespan.

Author Correction: GBA and APOE ε4 associate with sporadic dementia with Lewy bodies in European genome wide association study

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

Common brain disorders are associated with heritable patterns of apparent aging of the brain

Common risk factors for psychiatric and other brain disorders are likely to converge on biological pathways influencing the development and maintenance of brain structure and function across life. Using structural MRI data from 45,615 individuals aged 3-96 years, we demonstrate distinct patterns of apparent brain aging in several brain disorders and reveal genetic pleiotropy between apparent brain aging in healthy individuals and common brain disorders.

Autism spectrum disorder polygenic scores are associated with every day executive function in children admitted for clinical assessment

Autism spectrum disorder (ASD) and other neurodevelopmental disorders (NDs) are behaviorally defined disorders with overlapping clinical features that are often associated with higher‐order cognitive dysfunction, particularly executive dysfunction. Our aim was to determine if the polygenic score (PGS) for ASD is associated with parent‐reported executive dysfunction in everyday life using the Behavior Rating Inventory of Executive Function (BRIEF). Furthermore, we investigated if PGS for general intelligence (INT) and attention deficit/hyperactivity disorder (ADHD) also correlate with BRIEF. We included 176 children, adolescents and young adults aged 5–22 years with full‐scale intelligence quotient (IQ) above 70. All were admitted for clinical assessment of ASD symptoms and 68% obtained an ASD diagnosis. We found a significant difference between low and high ASD PGS groups in the BRIEF behavior regulation index (BRI) (P = 0.015, Cohen's d = 0.69). A linear regression model accounting for age, sex, full‐scale IQ, Social Responsiveness Scale (SRS) total score, ASD, ADHD and INT PGS groups as well as genetic principal components, significantly predicted the BRI score; F(11,130) = 8.142, P < 0.001, R² = 0.41 (unadjusted). Only SRS total (P < 0.001), ASD PGS 0.1 group (P = 0.018), and sex (P = 0.022) made a significant contribution to the model. This suggests that the common ASD risk gene variants have a stronger association to behavioral regulation aspects of executive dysfunction than ADHD risk or INT variants in a clinical sample with ASD symptoms. Autism Res 2020, 13: 207–220. © 2019 International Society for Autism Research, Wiley Periodicals, Inc.

Lay Summary

People with autism spectrum disorder (ASD) often have difficulties with higher‐order cognitive processes that regulate thoughts and actions during goal‐directed behavior, also known as executive function (EF). We studied the association between genetics related to ASD and EF and found a relation between high polygenic score (PGS) for ASD and difficulties with behavior regulation aspects of EF in children and adolescents under assessment for ASD. Furthermore, high PGS for general intelligence was related to social problems.

Brain Heterogeneity in Schizophrenia and Its Association With Polygenic Risk

Importance

Between-individual variability in brain structure is determined by gene-environment interactions, possibly reflecting differential sensitivity to environmental and genetic perturbations. Magnetic resonance imaging (MRI) studies have revealed thinner cortices and smaller subcortical volumes in patients with schizophrenia. However, group-level comparisons may mask considerable within-group heterogeneity, which has largely remained unnoticed in the literature.

Objectives

To compare brain structural variability between individuals with schizophrenia and healthy controls and to test whether respective variability reflects the polygenic risk score (PRS) for schizophrenia in an independent sample of healthy controls.

Design, Setting, and Participants

This case-control and polygenic risk analysis compared MRI-derived cortical thickness and subcortical volumes between healthy controls and patients with schizophrenia across 16 cohorts and tested for associations between PRS and MRI features in a control cohort from the UK Biobank. Data were collected from October 27, 2004, through April 12, 2018, and analyzed from December 3, 2017, through August 1, 2018.

Main Outcomes and Measures

Mean and dispersion parameters were estimated using double generalized linear models. Vertex-wise analysis was used to assess cortical thickness, and regions-of-interest analyses were used to assess total cortical volume, total surface area, and white matter, subcortical, and hippocampal subfield volumes. Follow-up analyses included within-sample analysis, test of robustness of the PRS threshold, population covariates, outlier removal, and control for image quality.

Results

A comparison of 1151 patients with schizophrenia (mean [SD] age, 33.8 [10.6] years; 68.6% male [n = 790] and 31.4% female [n = 361]) with 2010 healthy controls (mean [SD] age, 32.6 [10.4] years; 56.0% male [n = 1126] and 44.0% female [n = 884]) revealed higher heterogeneity in schizophrenia for cortical thickness and area (t = 3.34), cortical (t = 3.24) and ventricle (t range, 3.15-5.78) volumes, and hippocampal subfields (t range, 2.32-3.55). In the UK Biobank sample of 12 490 participants (mean [SD] age, 55.9 [7.5] years; 48.2% male [n = 6025] and 51.8% female [n = 6465]), higher PRS was associated with thinner frontal and temporal cortices and smaller left CA2/3 (t = -3.00) but was not significantly associated with dispersion.

Conclusions and Relevance

This study suggests that schizophrenia is associated with substantial brain structural heterogeneity beyond the mean differences. These findings may reflect higher sensitivity to environmental and genetic perturbations in patients, supporting the heterogeneous nature of schizophrenia. A higher PRS was associated with thinner frontotemporal cortices and smaller hippocampal subfield volume, but not heterogeneity. This finding suggests that brain variability in schizophrenia results from interactions between environmental and genetic factors that are not captured by the PRS. Factors contributing to heterogeneity in frontotemporal cortices and hippocampus are key to furthering our understanding of how genetic and environmental factors shape brain biology in schizophrenia.

Analysis of differentially methylated regions in great apes and extinct hominids provides support for the evolutionary hypothesis of schizophrenia

INTRODUCTION

The persistence of schizophrenia in human populations separated by geography and time led to the evolutionary hypothesis that proposes schizophrenia as a by-product of the higher cognitive abilities of modern humans. To explore this hypothesis, we used here an evolutionary epigenetics approach building on differentially methylated regions (DMRs) of the genome.

METHODS

We implemented a polygenic enrichment testing pipeline using the summary statistics of genome-wide association studies (GWAS) of schizophrenia and 12 other phenotypes. We investigated the enrichment of association of these traits across genomic regions with variable methylation between modern humans and great apes (orangutans, chimpanzees and gorillas; great ape DMRs) and between modern humans and recently extinct hominids (Neanderthals and Denisovans; hominid DMRs).

RESULTS

Regions that are hypo-methylated in humans compared to great apes show enrichment of association with schizophrenia only if the major histocompatibility complex (MHC) region is included. With the MHC region removed from the analysis, only a modest enrichment for SNPs of low effect persists. The INRICH pipeline confirms this finding after rigorous permutation and bootstrapping procedures.

CONCLUSION

The analyses of regions with differential methylation changes in humans and great apes do not provide compelling evidence of enrichment of association with schizophrenia, in contrast to our previous findings on more recent methylation differences between modern humans, Neanderthals and Denisovans. Our results further support the evolutionary hypothesis of schizophrenia and indicate that the origin of some of the genetic susceptibility factors of schizophrenia may lie in recent human evolution.

Identification of common genetic risk variants for autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable and heterogeneous group of neurodevelopmental phenotypes diagnosed in more than 1% of children. Common genetic variants contribute substantially to ASD susceptibility, but to date no individual variants have been robustly associated with ASD. With a marked sample-size increase from a unique Danish population resource, we report a genome-wide association meta-analysis of 18,381 individuals with ASD and 27,969 controls that identified five genome-wide-significant loci. Leveraging GWAS results from three phenotypes with significantly overlapping genetic architectures (schizophrenia, major depression, and educational attainment), we identified seven additional loci shared with other traits at equally strict significance levels. Dissecting the polygenic architecture, we found both quantitative and qualitative polygenic heterogeneity across ASD subtypes. These results highlight biological insights, particularly relating to neuronal function and corticogenesis, and establish that GWAS performed at scale will be much more productive in the near term in ASD.

Alterations in Schizophrenia-Associated Genes Can Lead to Increased Power in Delta Oscillations

Genome-wide association studies have implicated many ion channels in schizophrenia pathophysiology. Although the functions of these channels are relatively well characterized by single-cell studies, the contributions of common variation in these channels to neurophysiological biomarkers and symptoms of schizophrenia remain elusive. Here, using computational modeling, we show that a common biomarker of schizophrenia, namely, an increase in delta-oscillation power, may be a direct consequence of altered expression or kinetics of voltage-gated ion channels or calcium transporters. Our model of a circuit of layer V pyramidal cells highlights multiple types of schizophrenia-related variants that contribute to altered dynamics in the delta-frequency band. Moreover, our model predicts that the same membrane mechanisms that increase the layer V pyramidal cell network gain and response to delta-frequency oscillations may also cause a deficit in a single-cell correlate of the prepulse inhibition, which is a behavioral biomarker highly associated with schizophrenia.

Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer's disease risk

Alzheimer's disease (AD) is highly heritable and recent studies have identified over 20 disease-associated genomic loci. Yet these only explain a small proportion of the genetic variance, indicating that undiscovered loci remain. Here, we performed a large genome-wide association study of clinically diagnosed AD and AD-by-proxy (71,880 cases, 383,378 controls). AD-by-proxy, based on parental diagnoses, showed strong genetic correlation with AD (r_g = 0.81). Meta-analysis identified 29 risk loci, implicating 215 potential causative genes. Associated genes are strongly expressed in immune-related tissues and cell types (spleen, liver, and microglia). Gene-set analyses indicate biological mechanisms involved in lipid-related processes and degradation of amyloid precursor proteins. We show strong genetic correlations with multiple health-related outcomes, and Mendelian randomization results suggest a protective effect of cognitive ability on AD risk. These results are a step forward in identifying the genetic factors that contribute to AD risk and add novel insights into the neurobiology of AD.

Biophysical Psychiatry-How Computational Neuroscience Can Help to Understand the Complex Mechanisms of Mental Disorders

The brain is the most complex of human organs, and the pathophysiology underlying abnormal brain function in psychiatric disorders is largely unknown. Despite the rapid development of diagnostic tools and treatments in most areas of medicine, our understanding of mental disorders and their treatment has made limited progress during the last decades. While recent advances in genetics and neuroscience have a large potential, the complexity and multidimensionality of the brain processes hinder the discovery of disease mechanisms that would link genetic findings to clinical symptoms and behavior. This applies also to schizophrenia, for which genome-wide association studies have identified a large number of genetic risk loci, spanning hundreds of genes with diverse functionalities. Importantly, the multitude of the associated variants and their prevalence in the healthy population limit the potential of a reductionist functional genetics approach as a stand-alone solution to discover the disease pathology. In this review, we outline the key concepts of a "biophysical psychiatry," an approach that employs large-scale mechanistic, biophysics-founded computational modelling to increase transdisciplinary understanding of the pathophysiology and strive toward robust predictions. We discuss recent scientific advances that allow a synthesis of previously disparate fields of psychiatry, neurophysiology, functional genomics, and computational modelling to tackle open questions regarding the pathophysiology of heritable mental disorders. We argue that the complexity of the increasing amount of genetic data exceeds the capabilities of classical experimental assays and requires computational approaches. Biophysical psychiatry, based on modelling diseased brain networks using existing and future knowledge of basic genetic, biochemical, and functional properties on a single neuron to a microcircuit level, may allow a leap forward in deriving interpretable biomarkers and move the field toward novel treatment options.

Meta-analysis of Alzheimer's disease on 9,751 samples from Norway and IGAP study identifies four risk loci

A large fraction of genetic risk factors for Alzheimer's Disease (AD) is still not identified, limiting the understanding of AD pathology and study of therapeutic targets. We conducted a genome-wide association study (GWAS) of AD cases and controls of European descent from the multi-center DemGene network across Norway and two independent European cohorts. In a two-stage process, we first performed a meta-analysis using GWAS results from 2,893 AD cases and 6,858 cognitively normal controls from Norway and 25,580 cases and 48,466 controls from the International Genomics of Alzheimer's Project (IGAP), denoted the discovery sample. Second, we selected the top hits (p < 1 × 10^-6) from the discovery analysis for replication in an Icelandic cohort consisting of 5,341 cases and 110,008 controls. We identified a novel genomic region with genome-wide significant association with AD on chromosome 4 (combined analysis OR = 1.07, p = 2.48 x 10^-8). This finding implicated HS3ST1, a gene expressed throughout the brain particularly in the cerebellar cortex. In addition, we identified IGHV1-68 in the discovery sample, previously not associated with AD. We also associated USP6NL/ECHDC3 and BZRAP1-AS1 to AD, confirming findings from a follow-up transethnic study. These new gene loci provide further evidence for AD as a polygenic disorder, and suggest new mechanistic pathways that warrant further investigation.

Cross-tissue eQTL enrichment of associations in schizophrenia

The genome-wide association study of the Psychiatric Genomics Consortium identified over one hundred schizophrenia susceptibility loci. The number of non-coding variants discovered suggests that gene regulation could mediate the effect of these variants on disease. Expression quantitative trait loci (eQTLs) contribute to variation in levels of mRNA. Given the co-occurrence of schizophrenia and several traits not involving the central nervous system (CNS), we investigated the enrichment of schizophrenia associations among eQTLs for four non-CNS tissues: adipose tissue, epidermal tissue, lymphoblastoid cells and blood. Significant enrichment was seen in eQTLs of all tissues: adipose (β = 0.18, p = 8.8 × 10-06), epidermal (β = 0.12, p = 3.1 × 10-04), lymphoblastoid (β = 0.19, p = 6.2 × 10-08) and blood (β = 0.19, p = 6.4 × 10-06). For comparison, we looked for enrichment of association with traits of known relevance to one or more of these tissues (body mass index, height, rheumatoid arthritis, systolic blood pressure and type-II diabetes) and found that schizophrenia enrichment was of similar scale to that observed when studying diseases in the context of a more likely causal tissue. To further investigate tissue specificity, we looked for differential enrichment of eQTLs with relevant Roadmap affiliation (enhancers and promoters) and varying distance from the transcription start site. Neither factor significantly contributed to the enrichment, suggesting that this is equally distributed in tissue-specific and cross-tissue regulatory elements. Our analyses suggest that functional correlates of schizophrenia risk are prevalent in non-CNS tissues. This could be because of pleiotropy or the effectiveness of variants affecting expression in different contexts. This suggests the utility of large, single-tissue eQTL experiments to increase eQTL discovery power in the study of schizophrenia, in addition to smaller, multiple-tissue approaches. Our results conform to the notion that schizophrenia is a systemic disorder involving many tissues.

Vitamin D levels, brain volume, and genetic architecture in patients with psychosis

BACKGROUND

Lower vitamin D levels are found in people with schizophrenia and depressive disorders, and also associated with neuroimaging abnormalities such as reduced brain volume in both animals and humans. Reduced whole brain and increased ventricular volume are also systematically reported in schizophrenia. Even though vitamin D deficiency has been proposed as a risk mechanism for schizophrenia there exist no studies to date of the association between vitamin D levels and brain volume in this population. Therefore, we investigated the relationship between vitamin D levels and brain phenotypes in psychotic disorders, and assessed possible interactions with genetic variants in vitamin D receptor (VDR) and other genetic variants that play a role in vitamin D levels in the body.

METHODS

Our sample consisted of 83 psychosis patients and 101 healthy controls. We measured vitamin D levels as serum 25-hydroxyvitamin D. All participants were genotyped and neuroimaging conducted by structural magnetic resonance imaging.

RESULTS

Vitamin D levels were significantly positively associated with peripheral grey matter volume in patients (β 860.6; 95% confidence interval (CI) 333.4-1466, p < .003). A significant interaction effect of BSML marker (rs1544410) was observed to mediate the association between patient status and both white matter volume (β 23603.3; 95% CI 2732.8-48708.6, p < .05) and whole brain volume (β 46670.6, 95% CI 8817.8-93888.3, p < .04). Vitamin D did not predict ventricular volume, which rather was associated with patient status (β 4423.3, 95% CI 1583.2-7267.8p < .002) and CYP24A1 marker (rs6013897) (β 2491.5, 95% CI 269.7-4978.5, p < .04).

CONCLUSIONS

This is the first study of the association between vitamin D levels and brain volume in patients with psychotic disorders that takes into account possible interaction with genetic polymorphisms. The present findings warrant replication in independent samples.

Enrichment of genetic markers of recent human evolution in educational and cognitive traits

Higher cognitive functions are regarded as one of the main distinctive traits of humans. Evidence for the cognitive evolution of human beings is mainly based on fossil records of an expanding cranium and an increasing complexity of material culture artefacts. However, the molecular genetic factors involved in the evolution are still relatively unexplored. Here, we investigated whether genomic regions that underwent positive selection in humans after divergence from Neanderthals are enriched for genetic association with phenotypes related to cognitive functions. We used genome wide association data from a study of college completion (N = 111,114), one of educational attainment (N = 293,623) and two different studies of general cognitive ability (N = 269,867 and 53,949). We found nominally significant polygenic enrichment of associations with college completion (p = 0.025), educational attainment (p = 0.043) and general cognitive ability (p = 0.015 and 0.025, respectively), suggesting that variants influencing these phenotypes are more prevalent in evolutionarily salient regions. The enrichment remained significant after controlling for other known genetic enrichment factors, and for affiliation to genes highly expressed in the brain. These findings support the notion that phenotypes related to higher order cognitive skills typical of humans have a recent genetic component that originated after the separation of the human and Neanderthal lineages.

Effects of autozygosity and schizophrenia polygenic risk on cognitive and brain developmental trajectories

Cognitive and brain development are determined by dynamic interactions between genes and environment across the lifespan. Aside from marker-by-marker analyses of polymorphisms, biologically meaningful features of the whole genome (derived from the combined effect of individual markers) have been postulated to inform on human phenotypes including cognitive traits and their underlying biological substrate. Here, estimates of inbreeding and genetic susceptibility for schizophrenia calculated from genome-wide data-runs of homozygosity (ROH) and schizophrenia polygenic risk score (PGRS)-are analyzed in relation to cognitive abilities (n = 4183) and brain structure (n = 516) in a general-population sample of European-ancestry participants aged 8-22, from the Philadelphia Neurodevelopmental Cohort. The findings suggest that a higher ROH burden and higher schizophrenia PGRS are associated with higher intelligence. Cognition-ROH and cognition-PGRS associations obtained in this cohort may, respectively, evidence that assortative mating influences intelligence, and that individuals with high schizophrenia genetic risk who do not transition to disease status are cognitively resilient. Neuroanatomical data showed that the effects of schizophrenia PGRS on cognition could be modulated by brain structure, although larger imaging datasets are needed to accurately disentangle the underlying neural mechanisms linking IQ with both inbreeding and the genetic burden for schizophrenia.

Stability of the Brain Functional Connectome Fingerprint in Individuals With Schizophrenia

This study combines data from 3 existing functional magnetic resonance imaging data sets of people with schizophrenia spectrum disorders and healthy controls to assess the association of schizophrenia spectrum disorders with connectome stability.

Genetic Overlap Between Schizophrenia and Volumes of Hippocampus, Putamen, and Intracranial Volume Indicates Shared Molecular Genetic Mechanisms

Schizophrenia (SCZ) is associated with differences in subcortical brain volumes and intracranial volume (ICV). However, little is known about the underlying etiology of these brain alterations. Here, we explored whether brain structure volumes and SCZ share genetic risk factors. Using conditional false discovery rate (FDR) analysis, we integrated genome-wide association study (GWAS) data on SCZ (n = 82315) and GWAS data on 7 subcortical brain volumes and ICV (n = 11840). By conditioning the FDR on overlapping associations, this statistical approach increases power to discover genetic loci. To assess the credibility of our approach, we studied the identified loci in larger GWAS samples on ICV (n = 26577) and hippocampal volume (n = 26814). We observed polygenic overlap between SCZ and volumes of hippocampus, putamen, and ICV. Based on conjunctional FDR < 0.05, we identified 2 loci shared between SCZ and ICV implicating genes FOXO3 (rs10457180) and ITIH4 (rs4687658), 2 loci shared between SCZ and hippocampal volume implicating SLC4A10 (rs4664442) and SPATS2L (rs1653290), and 2 loci shared between SCZ and volume of putamen implicating DCC (rs4632195) and DLG2 (rs11233632). The loci shared between SCZ and hippocampal volume or ICV had not reached significance in the primary GWAS on brain phenotypes. Proving our point of increased power, 2 loci did reach genome-wide significance with ICV (rs10457180) and hippocampal volume (rs4664442) in the larger GWAS. Three of the 6 identified loci are novel for SCZ. Altogether, the findings provide new insights into the relationship between SCZ and brain structure volumes, suggesting that their genetic architectures are not independent.

A molecule-based genetic association approach implicates a range of voltage-gated calcium channels associated with schizophrenia

Traditional genome-wide association studies (GWAS) have successfully detected genetic variants associated with schizophrenia. However, only a small fraction of heritability can be explained. Gene-set/pathway-based methods can overcome limitations arising from single nucleotide polymorphism (SNP)-based analysis, but most of them place constraints on size which may exclude highly specific and functional sets, like macromolecules. Voltage-gated calcium (Ca_v ) channels, belonging to macromolecules, are composed of several subunits whose encoding genes are located far away or even on different chromosomes. We combined information about such molecules with GWAS data to investigate how functional channels associated with schizophrenia. We defined a biologically meaningful SNP-set based on channel structure and performed an association study by using a validated method: SNP-set (sequence) kernel association test. We identified eight subtypes of Ca_v channels significantly associated with schizophrenia from a subsample of published data (N = 56,605), including the L-type channels (Ca_v 1.1, Ca_v 1.2, Ca_v 1.3), P-/Q-type Ca_v 2.1, N-type Ca_v 2.2, R-type Ca_v 2.3, T-type Ca_v 3.1, and Ca_v 3.3. Only genes from Ca_v 1.2 and Ca_v 3.3 have been implicated by the largest GWAS (N = 82,315). Each subtype of Ca_v channels showed relatively high chip heritability, proportional to the size of its constituent gene regions. The results suggest that abnormalities of Ca_v channels may play an important role in the pathophysiology of schizophrenia and these channels may represent appropriate drug targets for therapeutics. Analyzing subunit-encoding genes of a macromolecule in aggregate is a complementary way to identify more genetic variants of polygenic diseases. This study offers the potential of power for discovery the biological mechanisms of schizophrenia.

Recently evolved human-specific methylated regions are enriched in schizophrenia signals

Background

One explanation for the persistence of schizophrenia despite the reduced fertility of patients is that it is a by-product of recent human evolution. This hypothesis is supported by evidence suggesting that recently-evolved genomic regions in humans are involved in the genetic risk for schizophrenia. Using summary statistics from genome-wide association studies (GWAS) of schizophrenia and 11 other phenotypes, we tested for enrichment of association with GWAS traits in regions that have undergone methylation changes in the human lineage compared to Neanderthals and Denisovans, i.e. human-specific differentially methylated regions (DMRs). We used analytical tools that evaluate polygenic enrichment of a subset of genomic variants against all variants.

Results

Schizophrenia was the only trait in which DMR SNPs showed clear enrichment of association that passed the genome-wide significance threshold. The enrichment was not observed for Neanderthal or Denisovan DMRs. The enrichment seen in human DMRs is comparable to that for genomic regions tagged by Neanderthal Selective Sweep markers, and stronger than that for Human Accelerated Regions. The enrichment survives multiple testing performed through permutation (n = 10,000) and bootstrapping (n = 5000) in INRICH (p < 0.01). Some enrichment of association with height was observed at the gene level.

Conclusions

Regions where DNA methylation modifications have changed during recent human evolution show enrichment of association with schizophrenia and possibly with height. Our study further supports the hypothesis that genetic variants conferring risk of schizophrenia co-occur in genomic regions that have changed as the human species evolved. Since methylation is an epigenetic mark, potentially mediated by environmental changes, our results also suggest that interaction with the environment might have contributed to that association.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1177-2) contains supplementary material, which is available to authorized users.

Association of Heritable Cognitive Ability and Psychopathology With White Matter Properties in Children and Adolescents

IMPORTANCE

Many mental disorders emerge during adolescence, which may reflect a cost of the potential for brain plasticity offered during this period. Brain dysconnectivity has been proposed as a common factor across diagnostic categories.

OBJECTIVE

To investigate the hypothesis that brain dysconnectivity is a transdiagnostic phenotype in adolescence with increased susceptibility and symptoms of psychiatric disease.

DESIGN, SETTING, AND PARTICIPANTS

We investigated clinical symptoms as well as cognitive function in 6487 individuals aged 8 to 21 years from November 1, 2009, to November 30, 2011, in the Philadelphia Neurodevelopmental Cohort and analyzed diffusion magnetic resonance imaging brain scans for 748 of the participants.

MAIN OUTCOMES AND MEASURES

Independent component analysis was used to derive dimensional psychopathology scores, and genome-wide complex trait analysis was used to estimate its heritability. Multimodal fusion simultaneously modeled contributions of the diffusion magnetic resonance imaging metrics fractional anisotropy, mean diffusivity, radial diffusivity, L1 (the principal diffusion tensor imaging eigen value), mode of anisotropy, as well as dominant and secondary fiber orientations, and structural connectivity density, and their association with general psychopathology and cognition.

RESULTS

Machine learning with 10-fold cross-validation and permutation testing in 729 individuals (aged 8 to 22 years; mean [SD] age, 15.1 [3.3] years; 343 females [46%]) revealed significant association with general psychopathology levels (r = 0.24, P < .001) and cognition (r = 0.39, P < .001). A brain white matter pattern reflecting frontotemporal connectivity and crossing fibers in the uncinate fasciculus was the most associated feature for both traits. Univariate analysis across a range of clinical domains and cognitive test scores confirmed its transdiagnostic importance. Both the general psychopathology (16%; SE, 0.095; P = .05) and cognitive (18%; SE, 0.09; P = .01) factor were heritable and showed a negative genetic correlation.

CONCLUSION AND RELEVANCE

Dimensional and heritable general cognitive and psychopathology factors are associated with specific patterns of white matter properties, suggesting that dysconnectivity is a transdiagnostic brain-based phenotype in individuals with increased susceptibility and symptoms of psychiatric disorders.

Novel Loci Associated With Attention-Deficit/Hyperactivity Disorder Are Revealed by Leveraging Polygenic Overlap With Educational Attainment

OBJECTIVE

Attention-deficit/hyperactivity disorder (ADHD) is a common and highly heritable psychiatric condition. By exploiting the reported relationship between ADHD and educational attainment (EA), we aimed to improve discovery of ADHD-associated genetic variants and to investigate genetic overlap between these phenotypes.

METHOD

A conditional/conjunctional false discovery rate (condFDR/conjFDR) method was applied to genome-wide association study (GWAS) data on ADHD (2,064 trios, 896 cases, and 2,455 controls) and EA (n=328,917) to identify ADHD-associated loci and loci overlapping between ADHD and EA. Identified single nucleotide polymorphisms (SNPs) were tested for association in an independent population-based study of ADHD symptoms (n=17,666). Genetic correlation between ADHD and EA was estimated using LD score regression and Pearson correlation.

RESULTS

At levels of condFDR<0.01 and conjFDR<0.05, we identified 5 ADHD-associated loci, 3 of these being shared between ADHD and EA. None of these loci had been identified in the primary ADHD GWAS, demonstrating the increased power provided by the condFDR/conjFDR analysis. Leading SNPs for 4 of 5 identified regions are in introns of protein coding genes (KDM4A, MEF2C, PINK1, RUNX1T1), whereas the remaining one is an intergenic SNP on chromosome 2 at 2p24. Consistent direction of effects in the independent study of ADHD symptoms was shown for 4 of 5 identified loci. A polygenic overlap between ADHD and EA was supported by significant genetic correlation (rg=-0.403, p=7.90×10-8) and >10-fold mutual enrichment of SNPs associated with both traits.

CONCLUSION

We identified 5 novel loci associated with ADHD and provided evidence for a shared genetic basis between ADHD and EA. These findings could aid understanding of the genetic risk architecture of ADHD and its relation to EA.

Identification of shared genetic variants between schizophrenia and lung cancer

Epidemiology studies suggest associations between schizophrenia and cancer. However, the underlying genetic mechanisms are not well understood, and difficult to identify from epidemiological data. We investigated if there is a shared genetic architecture between schizophrenia and cancer, with the aim to identify specific overlapping genetic loci. First, we performed genome-wide enrichment analysis and second, we analyzed specific loci jointly associated with schizophrenia and cancer by the conjunction false discovery rate. We analyzed the largest genome-wide association studies of schizophrenia and lung, breast, prostate, ovary, and colon-rectum cancer including more than 220,000 subjects, and included genetic association with smoking behavior. Polygenic enrichment of associations with lung cancer was observed in schizophrenia, and weak enrichment for the remaining cancer sites. After excluding the major histocompatibility complex region, we identified three independent loci jointly associated with schizophrenia and lung cancer. The strongest association included nicotinic acetylcholine receptors and is an established pleiotropic locus shared between lung cancer and smoking. The two other loci were independent of genetic association with smoking. Functional analysis identified downstream pleiotropic effects on epigenetics and gene-expression in lung and brain tissue. These findings suggest that genetic factors may explain partly the observed epidemiological association of lung cancer and schizophrenia.

Psychotic patients who used cannabis frequently before illness onset have higher genetic predisposition to schizophrenia than those who did not

BACKGROUND

Schizophrenia (SZ) and bipolar disorder (BD) are heritable, polygenic disorders with shared clinical and genetic components, suggesting a psychosis continuum. Cannabis use is a well-documented environmental risk factor in psychotic disorders. In the current study, we investigated the relationship between SZ genetic load and cannabis use before illness onset in SZ and BD spectrums. Since frequent early cannabis use (age <18 years) is believed to increase the risk of developing psychosis more than later use, follow-up analyses were conducted comparing early use to later use and no use.

METHODS

We assigned a SZ-polygenic risk score (PGRS) to each individual in our independent sample (N = 381 SZ spectrum cases, 220 BD spectrum cases and 415 healthy controls), calculated from the results of the Psychiatric Genomics Consortium (PGC) SZ case-control study (N = 81 535). SZ-PGRS in patients who used cannabis weekly to daily in the period before first illness episode was compared with that of those who never or infrequently used cannabis.

RESULTS

Patients with weekly to daily cannabis use before illness onset had the highest SZ-PGRS (p = 0.02, Cohen's d = 0.33). The largest difference was found between patients with daily or weekly cannabis use before illness onset <18 years of age and patients with no or infrequent use of cannabis (p = 0.003, Cohen's d = 0.42).

CONCLUSIONS

Our study supports an association between high SZ-PGRS and frequent cannabis use before illness onset in psychosis continuum disorders.

Identification of Genetic Loci Jointly Influencing Schizophrenia Risk and the Cognitive Traits of Verbal-Numerical Reasoning, Reaction Time, and General Cognitive Function

IMPORTANCE

Schizophrenia is associated with widespread cognitive impairments. Although cognitive deficits are one of the factors most strongly associated with functional outcome in schizophrenia, current treatment strategies largely fail to ameliorate these impairments. To develop more efficient treatment strategies in patients with schizophrenia, a better understanding of the pathogenesis of these cognitive deficits is needed. Accumulating evidence indicates that genetic risk of schizophrenia may contribute to cognitive dysfunction.

OBJECTIVE

To identify genomic regions jointly influencing schizophrenia and the cognitive domains of reaction time and verbal-numerical reasoning, as well as general cognitive function, a phenotype that captures the shared variation in performance across cognitive domains.

DESIGN, SETTING, AND PARTICIPANTS

Combining data from genome-wide association studies from multiple phenotypes using conditional false discovery rate analysis provides increased power to discover genetic variants and could elucidate shared molecular genetic mechanisms. Data from the following genome-wide association studies, published from July 24, 2014, to January 17, 2017, were combined: schizophrenia in the Psychiatric Genomics Consortium cohort (n = 79 757 [cases, 34 486; controls, 45 271]); verbal-numerical reasoning (n = 36 035) and reaction time (n = 111 483) in the UK Biobank cohort; and general cognitive function in CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) (n = 53 949) and COGENT (Cognitive Genomics Consortium) (n = 27 888).

MAIN OUTCOMES AND MEASURES

Genetic loci identified by conditional false discovery rate analysis. Brain messenger RNA expression and brain expression quantitative trait locus functionality were determined.

RESULTS

Among the participants in the genome-wide association studies, 21 loci jointly influencing schizophrenia and cognitive traits were identified: 2 loci shared between schizophrenia and verbal-numerical reasoning, 6 loci shared between schizophrenia and reaction time, and 14 loci shared between schizophrenia and general cognitive function. One locus was shared between schizophrenia and 2 cognitive traits and represented the strongest shared signal detected (nearest gene TCF20; chromosome 22q13.2), and was shared between schizophrenia (z score, 5.01; P = 5.53 × 10-7), general cognitive function (z score, -4.43; P = 9.42 × 10-6), and verbal-numerical reasoning (z score, -5.43; P = 5.64 × 10-8). For 18 loci, schizophrenia risk alleles were associated with poorer cognitive performance. The implicated genes are expressed in the developmental and adult human brain. Replicable expression quantitative trait locus functionality was identified for 4 loci in the adult human brain.

CONCLUSIONS AND RELEVANCE

The discovered loci improve the understanding of the common genetic basis underlying schizophrenia and cognitive function, suggesting novel molecular genetic mechanisms.

Shared genetic risk between migraine and coronary artery disease: A genome-wide analysis of common variants

Migraine is a recurrent pain condition traditionally viewed as a neurovascular disorder, but little is known of its vascular basis. In epidemiological studies migraine is associated with an increased risk of cardiovascular disease, including coronary artery disease (CAD), suggesting shared pathogenic mechanisms. This study aimed to determine the genetic overlap between migraine and CAD, and to identify shared genetic risk loci, utilizing a conditional false discovery rate approach and data from two large-scale genome-wide association studies (GWAS) of CAD (C4D, 15,420 cases, 15,062 controls; CARDIoGRAM, 22,233 cases, 64,762 controls) and one of migraine (22,120 cases, 91,284 controls). We found significant enrichment of genetic variants associated with CAD as a function of their association with migraine, which was replicated across two independent CAD GWAS studies. One shared risk locus in the PHACTR1 gene (conjunctional false discovery rate for index SNP rs9349379 < 3.90 x 10⁻⁵), which was also identified in previous studies, explained much of the enrichment. Two further loci (in KCNK5 and AS3MT) showed evidence for shared risk (conjunctional false discovery rate < 0.05). The index SNPs at two of the three loci had opposite effect directions in migraine and CAD. Our results confirm previous reports that migraine and CAD share genetic risk loci in excess of what would be expected by chance, and highlight one shared risk locus in PHACTR1. Understanding the biological mechanisms underpinning this shared risk is likely to improve our understanding of both disorders.

Genome-wide Pleiotropy Between Parkinson Disease and Autoimmune Diseases

Importance

Recent genome-wide association studies (GWAS) and pathway analyses supported long-standing observations of an association between immune-mediated diseases and Parkinson disease (PD). The post-GWAS era provides an opportunity for cross-phenotype analyses between different complex phenotypes.

Objectives

To test the hypothesis that there are common genetic risk variants conveying risk of both PD and autoimmune diseases (ie, pleiotropy) and to identify new shared genetic variants and their pathways by applying a novel statistical framework in a genome-wide approach.

Design, Setting, and Participants

Using the conjunction false discovery rate method, this study analyzed GWAS data from a selection of archetypal autoimmune diseases among 138 511 individuals of European ancestry and systemically investigated pleiotropy between PD and type 1 diabetes, Crohn disease, ulcerative colitis, rheumatoid arthritis, celiac disease, psoriasis, and multiple sclerosis. NeuroX data (6927 PD cases and 6108 controls) were used for replication. The study investigated the biological correlation between the top loci through protein-protein interaction and changes in the gene expression and methylation levels. The dates of the analysis were June 10, 2015, to March 4, 2017.

Main Outcomes and Measures

The primary outcome was a list of novel loci and their pathways involved in PD and autoimmune diseases.

Results

Genome-wide conjunctional analysis identified 17 novel loci at false discovery rate less than 0.05 with overlap between PD and autoimmune diseases, including known PD loci adjacent to GAK, HLA-DRB5, LRRK2, and MAPT for rheumatoid arthritis, ulcerative colitis and Crohn disease. Replication confirmed the involvement of HLA, LRRK2, MAPT, TRIM10, and SETD1A in PD. Among the novel genes discovered, WNT3, KANSL1, CRHR1, BOLA2, and GUCY1A3 are within a protein-protein interaction network with known PD genes. A subset of novel loci was significantly associated with changes in methylation or expression levels of adjacent genes.

Conclusions and Relevance

The study findings provide novel mechanistic insights into PD and autoimmune diseases and identify a common genetic pathway between these phenotypes. The results may have implications for future therapeutic trials involving anti-inflammatory agents.

Distinct multivariate brain morphological patterns and their added predictive value with cognitive and polygenic risk scores in mental disorders

The brain underpinnings of schizophrenia and bipolar disorders are multidimensional, reflecting complex pathological processes and causal pathways, requiring multivariate techniques to disentangle. Furthermore, little is known about the complementary clinical value of brain structural phenotypes when combined with data on cognitive performance and genetic risk. Using data-driven fusion of cortical thickness, surface area, and gray matter density maps (GMD), we found six biologically meaningful patterns showing strong group effects, including four statistically independent multimodal patterns reflecting co-occurring alterations in thickness and GMD in patients, over and above two other independent patterns of widespread thickness and area reduction. Case-control classification using cognitive scores alone revealed high accuracy, and adding imaging features or polygenic risk scores increased performance, suggesting their complementary predictive value with cognitive scores being the most sensitive features. Multivariate pattern analyses reveal distinct patterns of brain morphology in mental disorders, provide insights on the relative importance between brain structure, cognitive and polygenetic risk score in classification of patients, and demonstrate the importance of multivariate approaches in studying the pathophysiological substrate of these complex disorders.

•

Linked ICA showed six independent multivariate morphology patterns sensitive to SZ.

•

Machine learning used to compare brain structure, cognitive and genetic scores.

•

Cognition showed highest prediction of SZ, boosted by brain structure or genetics.

Identification of Gene Loci That Overlap Between Schizophrenia and Educational Attainment

There is evidence for genetic overlap between cognitive abilities and schizophrenia (SCZ), and genome-wide association studies (GWAS) demonstrate that both SCZ and general cognitive abilities have a strong polygenic component with many single-nucleotide polymorphisms (SNPs) each with a small effect. Here we investigated the shared genetic architecture between SCZ and educational attainment, which is regarded as a "proxy phenotype" for cognitive abilities, but may also reflect other traits. We applied a conditional false discovery rate (condFDR) method to GWAS of SCZ (n = 82 315), college completion ("College," n = 95 427), and years of education ("EduYears," n = 101 069). Variants associated with College or EduYears showed enrichment of association with SCZ, demonstrating polygenic overlap. This was confirmed by an increased replication rate in SCZ. By applying a condFDR threshold <0.01, we identified 18 genomic loci associated with SCZ after conditioning on College and 15 loci associated with SCZ after conditioning on EduYears. Ten of these loci overlapped. Using conjunctional FDR, we identified 10 loci shared between SCZ and College, and 29 loci shared between SCZ and EduYears. The majority of these loci had effects in opposite directions. Our results provide evidence for polygenic overlap between SCZ and educational attainment, and identify novel pleiotropic loci. Other studies have reported genetic overlap between SCZ and cognition, or SCZ and educational attainment, with negative correlation. Importantly, our methods enable identification of bi-directional effects, which highlight the complex relationship between SCZ and educational attainment, and support polygenic mechanisms underlying both cognitive dysfunction and creativity in SCZ.

Bromodomain protein 4 discriminates tissue-specific super-enhancers containing disease-specific susceptibility loci in prostate and breast cancer

BACKGROUND

Epigenetic information can be used to identify clinically relevant genomic variants single nucleotide polymorphisms (SNPs) of functional importance in cancer development. Super-enhancers are cell-specific DNA elements, acting to determine tissue or cell identity and driving tumor progression. Although previous approaches have been tried to explain risk associated with SNPs in regulatory DNA elements, so far epigenetic readers such as bromodomain containing protein 4 (BRD4) and super-enhancers have not been used to annotate SNPs. In prostate cancer (PC), androgen receptor (AR) binding sites to chromatin have been used to inform functional annotations of SNPs.

RESULTS

Here we establish criteria for enhancer mapping which are applicable to other diseases and traits to achieve the optimal tissue-specific enrichment of PC risk SNPs. We used stratified Q-Q plots and Fisher test to assess the differential enrichment of SNPs mapping to specific categories of enhancers. We find that BRD4 is the key discriminant of tissue-specific enhancers, showing that it is more powerful than AR binding information to capture PC specific risk loci, and can be used with similar effect in breast cancer (BC) and applied to other diseases such as schizophrenia.

CONCLUSIONS

This is the first study to evaluate the enrichment of epigenetic readers in genome-wide associations studies for SNPs within enhancers, and provides a powerful tool for enriching and prioritizing PC and BC genetic risk loci. Our study represents a proof of principle applicable to other diseases and traits that can be used to redefine molecular mechanisms of human phenotypic variation.