Search

BACKGROUND: Geographical variations in mood and psychotic disorders have been found in upper-income countries. We looked for geographic variation in these disorders in Colombia, a middle-income country. We analyzed electronic health records from the Clínica San Juan de Dios Manizales (CSJDM), which provides comprehensive mental healthcare for the one million inhabitants of Caldas. METHODS: We constructed a friction surface map of Caldas and used it to calculate the travel-time to the CSJDM for 16,295 patients who had received an initial diagnosis of mood or psychotic disorder. Using a zero-inflated negative binomial regression model, we determined the relationship between travel-time and incidence, stratified by disease severity. We employed spatial scan statistics to look for patient clusters. RESULTS: We show that travel-times (for driving) to the CSJDM are less than 1 h for ~50% of the population and more than 4 h for ~10%. We find a distance-decay relationship for outpatients, but not for inpatients: for every hour increase in travel-time, the number of expected outpatient cases decreases by 20% (RR = 0.80, 95% confidence interval [0.71, 0.89], p = 5.67E-05). We find nine clusters/hotspots of inpatients. CONCLUSIONS: Our results reveal inequities in access to healthcare: many individuals requiring only outpatient treatment may live too far from the CSJDM to access healthcare. Targeting of resources to comprehensively identify severely ill individuals living in the observed hotspots could further address treatment inequities and enable investigations to determine factors generating these hotspots.

Polygenic scores (PGSs) have limited portability across different groupings of individuals (for example, by genetic ancestries and/or social determinants of health), preventing their equitable use1-3. PGS portability has typically been assessed using a single aggregate population-level statistic (for example, R2)4, ignoring inter-individual variation within the population. Here, using a large and diverse Los Angeles biobank5 (ATLAS, n = 36,778) along with the UK Biobank6 (UKBB, n = 487,409), we show that PGS accuracy decreases individual-to-individual along the continuum of genetic ancestries7 in all considered populations, even within traditionally labelled 'homogeneous' genetic ancestries. The decreasing trend is well captured by a continuous measure of genetic distance (GD) from the PGS training data: Pearson correlation of -0.95 between GD and PGS accuracy averaged across 84 traits. When applying PGS models trained on individuals labelled as white British in the UKBB to individuals with European ancestries in ATLAS, individuals in the furthest GD decile have 14% lower accuracy relative to the closest decile; notably, the closest GD decile of individuals with Hispanic Latino American ancestries show similar PGS performance to the furthest GD decile of individuals with European ancestries. GD is significantly correlated with PGS estimates themselves for 82 of 84 traits, further emphasizing the importance of incorporating the continuum of genetic ancestries in PGS interpretation. Our results highlight the need to move away from discrete genetic ancestry clusters towards the continuum of genetic ancestries when considering PGSs.

Mapping genetic variants that regulate gene expression (eQTL mapping) in large-scale RNA sequencing (RNA-seq) studies is often employed to understand functional consequences of regulatory variants. However, the high cost of RNA-seq limits sample size, sequencing depth, and, therefore, discovery power in eQTL studies. In this work, we demonstrate that, given a fixed budget, eQTL discovery power can be increased by lowering the sequencing depth per sample and increasing the number of individuals sequenced in the assay. We perform RNA-seq of whole-blood tissue across 1,490 individuals at low coverage (5.9 million reads/sample) and show that the effective power is higher than that of an RNA-seq study of 570 individuals at moderate coverage (13.9 million reads/sample). Next, we leverage synthetic datasets derived from real RNA-seq data (50 million reads/sample) to explore the interplay of coverage and number individuals in eQTL studies, and show that a 10-fold reduction in coverage leads to only a 2.5-fold reduction in statistical power to identify eQTLs. Our work suggests that lowering coverage while increasing the number of individuals in RNA-seq is an effective approach to increase discovery power in eQTL studies.

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.

BackgroundStudying phenotypic and genetic characteristics of age at onset (AAO) and polarity at onset (PAO) in bipolar disorder can provide new insights into disease pathology and facilitate the development of screening tools.AimsTo examine the genetic architecture of AAO and PAO and their association with bipolar disorder disease characteristics.MethodGenome-wide association studies (GWASs) and polygenic score (PGS) analyses of AAO (n = 12 977) and PAO (n = 6773) were conducted in patients with bipolar disorder from 34 cohorts and a replication sample (n = 2237). The association of onset with disease characteristics was investigated in two of these cohorts.ResultsEarlier AAO was associated with a higher probability of psychotic symptoms, suicidality, lower educational attainment, not living together and fewer episodes. Depressive onset correlated with suicidality and manic onset correlated with delusions and manic episodes. Systematic differences in AAO between cohorts and continents of origin were observed. This was also reflected in single-nucleotide variant-based heritability estimates, with higher heritabilities for stricter onset definitions. Increased PGS for autism spectrum disorder (β = -0.34 years, s.e. = 0.08), major depression (β = -0.34 years, s.e. = 0.08), schizophrenia (β = -0.39 years, s.e. = 0.08), and educational attainment (β = -0.31 years, s.e. = 0.08) were associated with an earlier AAO. The AAO GWAS identified one significant locus, but this finding did not replicate. Neither GWAS nor PGS analyses yielded significant associations with PAO.ConclusionsAAO and PAO are associated with indicators of bipolar disorder severity. Individuals with an earlier onset show an increased polygenic liability for a broad spectrum of psychiatric traits. Systematic differences in AAO across cohorts, continents and phenotype definitions introduce significant heterogeneity, affecting analyses.

The most significant common variant association for schizophrenia (SCZ) reflects increased expression of the complement component 4A (C4A). Yet, it remains unclear how C4A interacts with other SCZ risk genes or whether the complement system more broadly is implicated in SCZ pathogenesis. Here, we integrate several existing, large-scale genetic and transcriptomic datasets to interrogate the functional role of the complement system and C4A in the human brain. Unexpectedly, we find no significant genetic enrichment among known complement system genes for SCZ. Conversely, brain co-expression network analyses using C4A as a seed gene reveal that genes downregulated when C4A expression increases exhibit strong and specific genetic enrichment for SCZ risk. This convergent genomic signal reflects synaptic processes, is sexually dimorphic and most prominent in frontal cortical brain regions, and is accentuated by smoking. Overall, these results indicate that synaptic pathways-rather than the complement system-are the driving force conferring SCZ risk.

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.

Psychotic symptoms are not only an important feature of severe neuropsychiatric disorders, but are also common in the general population, especially in youth. The genetic etiology of psychosis symptoms in youth remains poorly understood. To characterize genetic risk for psychosis spectrum symptoms (PS), we leverage a community-based multiethnic sample of children and adolescents aged 8-22 years, the Philadelphia Neurodevelopmental Cohort (n = 7225, 20% PS). Using an elastic net regression model, we aim to classify PS status using polygenic scores (PGS) based on a range of heritable psychiatric and brain-related traits in a multi-PGS model. We also perform univariate PGS associations and evaluate age-specific effects. The multi-PGS analyses do not improve prediction of PS status over univariate models, but reveal that the attention deficit hyperactivity disorder (ADHD) PGS is robustly and uniquely associated with PS (OR 1.12 (1.05, 1.18) P = 0.0003). This association is driven by subjects of European ancestry (OR = 1.23 (1.14, 1.34), P = 4.15 × 10-7) but is not observed in African American subjects (P = 0.65). We find a significant interaction of ADHD PGS with age (P = 0.01), with a stronger association in younger children. The association is independent of phenotypic overlap between ADHD and PS, not indirectly driven by substance use or childhood trauma, and appears to be specific to PS rather than reflecting general psychopathology in youth. In an independent sample, we replicate an increased ADHD PGS in 328 youth at clinical high risk for psychosis, compared to 216 unaffected controls (OR 1.06, CI(1.01, 1.11), P = 0.02). Our findings suggest that PS in youth may reflect a different genetic etiology than psychotic symptoms in adulthood, one more akin to ADHD, and shed light on how genetic risk can be investigated across early disease trajectories.

Large-scale studies spanning diverse project sites, populations, languages, and measurements are increasingly important to relate psychological to biological variables. National and international consortia already are collecting and executing mega-analyses on aggregated data from individuals, with different measures on each person. In this research, we show that Asparouhov and Muthén's alignment method can be adapted to align data from disparate item sets and response formats. We argue that with these adaptations, the alignment method is well suited for combining data across multiple sites even when they use different measurement instruments. The approach is illustrated using data from the Whole Genome Sequencing in Psychiatric Disorders consortium and a real-data-based simulation is used to verify accurate parameter recovery. Factor alignment appears to increase precision of measurement and validity of scores with respect to external criteria. The resulting parameter estimates may further inform development of more effective and efficient methods to assess the same constructs in prospectively designed studies.

Many common illnesses, for reasons that have not been identified, differentially affect men and women. For instance, the autoimmune diseases systemic lupus erythematosus (SLE) and Sjögren's syndrome affect nine times more women than men1, whereas schizophrenia affects men with greater frequency and severity relative to women2. All three illnesses have their strongest common genetic associations in the major histocompatibility complex (MHC) locus, an association that in SLE and Sjögren's syndrome has long been thought to arise from alleles of the human leukocyte antigen (HLA) genes at that locus3-6. Here we show that variation of the complement component 4 (C4) genes C4A and C4B, which are also at the MHC locus and have been linked to increased risk for schizophrenia7, generates 7-fold variation in risk for SLE and 16-fold variation in risk for Sjögren's syndrome among individuals with common C4 genotypes, with C4A protecting more strongly than C4B in both illnesses. The same alleles that increase risk for schizophrenia greatly reduce risk for SLE and Sjögren's syndrome. In all three illnesses, C4 alleles act more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for SLE, 31-fold variation in risk for Sjögren's syndrome, and 1.7-fold variation in schizophrenia risk among men (versus 6-fold, 15-fold and 1.26-fold variation in risk among women, respectively). At a protein level, both C4 and its effector C3 were present at higher levels in cerebrospinal fluid and plasma8,9 in men than in women among adults aged between 20 and 50 years, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help to explain the more potent effects of C4 alleles in men, women's greater risk of SLE and Sjögren's syndrome and men's greater vulnerability to schizophrenia. These results implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

BackgroundSevere mental illness diagnoses have overlapping symptomatology and shared genetic risk, motivating cross-diagnostic investigations of disease-relevant quantitative measures. We analysed relationships between neurocognitive performance, symptom domains, and diagnoses in a large sample of people with severe mental illness not ascertained for a specific diagnosis (cases), and people without mental illness (controls) from a single, homogeneous population.MethodsIn this case-control study, cases with severe mental illness were ascertained through electronic medical records at Clínica San Juan de Dios de Manizales (Manizales, Caldas, Colombia) and the Hospital Universitario San Vicente Fundación (Medellín, Antioquía, Colombia). Participants were assessed for speed and accuracy using the Penn Computerized Neurocognitive Battery (CNB). Cases had structured interview-based diagnoses of schizophrenia, bipolar 1, bipolar 2, or major depressive disorder. Linear mixed models, using CNB tests as repeated measures, modelled neurocognition as a function of diagnosis, sex, and all interactions. Follow-up analyses in cases included symptom factor scores obtained from exploratory factor analysis of symptom data as main effects.FindingsBetween Oct 1, 2017, and Nov 1, 2019, 2406 participants (1689 cases [schizophrenia n=160; bipolar 1 disorder n=519; bipolar 2 disorder n=204; and major depressive disorder n=806] and 717 controls; mean age 39 years (SD 14); and 1533 female) were assessed. Participants with bipolar 1 disorder and schizophrenia had similar impairments in accuracy and speed across cognitive domains. Participants with bipolar 2 disorder and major depressive disorder performed similarly to controls, with subtle deficits in executive and social cognition. A three-factor model (psychosis, mania, and depression) best represented symptom data. Controlling for diagnosis, premorbid IQ, and disease severity, high lifetime psychosis scores were associated with reduced accuracy and speed across cognitive domains, whereas high depression scores were associated with increased social cognition accuracy.InterpretationCross-diagnostic investigations showed that neurocognitive function in severe mental illness is characterised by two distinct profiles (bipolar 1 disorder and schizophrenia, and bipolar 2 disorder and major depressive disorder), and is associated with specific symptom domains. These results suggest the utility of this design for elucidating severe mental illness causes and trajectories.FundingUS National Institute of Mental Health.

Carriers of large recurrent copy number variants (CNVs) have a higher risk of developing neurodevelopmental disorders. The 16p11.2 distal CNV predisposes carriers to e.g., autism spectrum disorder and schizophrenia. We compared subcortical brain volumes of 12 16p11.2 distal deletion and 12 duplication carriers to 6882 non-carriers from the large-scale brain Magnetic Resonance Imaging collaboration, ENIGMA-CNV. After stringent CNV calling procedures, and standardized FreeSurfer image analysis, we found negative dose-response associations with copy number on intracranial volume and on regional caudate, pallidum and putamen volumes (β = -0.71 to -1.37; P

Prior to and following the publication of this article the authors noted that the complete list of authors was not included in the main article and was only present in Supplementary Table 1. The author list in the original article has now been updated to include all authors, and Supplementary Table 1 has been removed. All other supplementary files have now been updated accordingly. Furthermore, in Table 1 of this Article, the replication cohort for the row Close relative in data set, n (%) was incorrect. All values have now been corrected to 0(0%). The publishers would like to apologise for this error and the inconvenience it may have caused.

BackgroundCognitive 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.MethodsWe 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.ResultsPRS 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.ConclusionsCognition 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.

Psychosis spectrum disorders are conceptualized as neurodevelopmental disorders accompanied by disruption of large-scale functional brain networks. Dynamic functional dysconnectivity has been described in patients with schizophrenia and in help-seeking individuals at clinical high risk for psychosis. Less is known, about developmental aspects of dynamic functional network connectivity (dFNC) associated with psychotic symptoms (PS) in the general population. Here, we investigate resting state functional magnetic resonance imaging data using established dFNC methods in the Philadelphia Neurodevelopmental Cohort (ages 8-22 years), including 129 participants experiencing PS and 452 participants without PS (non-PS). Functional networks were identified using group spatial independent component analysis. A sliding window approach and k-means clustering were applied to covariance matrices of all functional networks to identify recurring whole-brain connectivity states. PS-associated dysconnectivity of default mode, salience, and executive networks occurred only in a few states, whereas dysconnectivity in the sensorimotor and visual systems in PS youth was more pervasive, observed across multiple states. This study provides new evidence that disruptions of dFNC are present even at the less severe end of the psychosis continuum in youth, complementing previous work on help-seeking and clinically diagnosed cohorts that represent the more severe end of this spectrum.

Current evidence from case/control studies indicates that genetic risk for psychiatric disorders derives primarily from numerous common variants, each with a small phenotypic impact. The literature describing apparent segregation of bipolar disorder (BP) in numerous multigenerational pedigrees suggests that, in such families, large-effect inherited variants might play a greater role. To identify roles of rare and common variants on BP, we conducted genetic analyses in 26 Colombia and Costa Rica pedigrees ascertained for bipolar disorder 1 (BP1), the most severe and heritable form of BP. In these pedigrees, we performed microarray SNP genotyping of 838 individuals and high-coverage whole-genome sequencing of 449 individuals. We compared polygenic risk scores (PRS), estimated using the latest BP1 genome-wide association study (GWAS) summary statistics, between BP1 individuals and related controls. We also evaluated whether BP1 individuals had a higher burden of rare deleterious single-nucleotide variants (SNVs) and rare copy number variants (CNVs) in a set of genes related to BP1. We found that compared with unaffected relatives, BP1 individuals had higher PRS estimated from BP1 GWAS statistics (P = 0.001 ~ 0.007) and displayed modest increase in burdens of rare deleterious SNVs (P = 0.047) and rare CNVs (P = 0.002 ~ 0.033) in genes related to BP1. We did not observe rare variants segregating in the pedigrees. These results suggest that small-to-moderate effect rare and common variants are more likely to contribute to BP1 risk in these extended pedigrees than a few large-effect rare variants.

ObjectiveThe 2-year risk of psychosis in persons who meet research criteria for a high-risk syndrome is about 15%-25%; improvements in risk prediction accuracy would benefit the development and implementation of preventive interventions. The authors sought to assess polygenic risk score (PRS) prediction of subsequent psychosis in persons at high risk and to determine the impact of adding the PRS to a previously validated psychosis risk calculator.MethodsPersons meeting research criteria for psychosis high risk (N=764) and unaffected individuals (N=279) were followed for up to 2 years. The PRS was based on the latest schizophrenia and bipolar genome-wide association studies. Variables in the psychosis risk calculator included stressful life events, trauma, disordered thought content, verbal learning, information processing speed, and family history of psychosis.ResultsFor Europeans, the PRS varied significantly by group and was higher in the psychosis converter group compared with both the nonconverter and unaffected groups, but was similar for the nonconverter group compared with the unaffected group. For non-Europeans, the PRS varied significantly by group; the difference between the converters and nonconverters was not significant, but the PRS was significantly higher in converters than in unaffected individuals, and it did not differ between nonconverters and unaffected individuals. The R2liability (R2 adjusted for the rate of disease risk in the population being studied, here assuming a 2-year psychosis risk between 10% and 30%) for Europeans varied between 9.2% and 12.3% and for non-Europeans between 3.5% and 4.8%. The amount of risk prediction information contributed by the addition of the PRS to the risk calculator was less than severity of disordered thoughts and similar to or greater than for other variables. For Europeans, the PRS was correlated with risk calculator variables of information processing speed and verbal memory.ConclusionsThe PRS discriminates psychosis converters from nonconverters and modestly improves individualized psychosis risk prediction when added to a psychosis risk calculator. The schizophrenia PRS shows promise in enhancing risk prediction in persons at high risk for psychosis, although its potential utility is limited by poor performance in persons of non-European ancestry.

BackgroundRapid, preoperative identification of patients with the highest risk for medical complications is necessary to ensure that limited infrastructure and human resources are directed towards those most likely to benefit. Existing risk scores either lack specificity at the patient level or utilise the American Society of Anesthesiologists (ASA) physical status classification, which requires a clinician to review the chart.MethodsWe report on the use of machine learning algorithms, specifically random forests, to create a fully automated score that predicts postoperative in-hospital mortality based solely on structured data available at the time of surgery. Electronic health record data from 53 097 surgical patients (2.01% mortality rate) who underwent general anaesthesia between April 1, 2013 and December 10, 2018 in a large US academic medical centre were used to extract 58 preoperative features.ResultsUsing a random forest classifier we found that automatically obtained preoperative features (area under the curve [AUC] of 0.932, 95% confidence interval [CI] 0.910-0.951) outperforms Preoperative Score to Predict Postoperative Mortality (POSPOM) scores (AUC of 0.660, 95% CI 0.598-0.722), Charlson comorbidity scores (AUC of 0.742, 95% CI 0.658-0.812), and ASA physical status (AUC of 0.866, 95% CI 0.829-0.897). Including the ASA physical status with the preoperative features achieves an AUC of 0.936 (95% CI 0.917-0.955).ConclusionsThis automated score outperforms the ASA physical status score, the Charlson comorbidity score, and the POSPOM score for predicting in-hospital mortality. Additionally, we integrate this score with a previously published postoperative score to demonstrate the extent to which patient risk changes during the perioperative period.

BackgroundCommon psychiatric disorders are characterized by complex disease architectures with many small genetic effects that contribute and complicate biological understanding of their etiology. There is therefore a pressing need for in vitro experimental systems that allow for interrogation of polygenic psychiatric disease risk to study the underlying biological mechanisms.MethodsWe have developed an analytical framework that integrates genome-wide disease risk from genome-wide association studies with longitudinal in vitro gene expression profiles of human neuronal differentiation.ResultsWe demonstrate that the cumulative impact of risk loci of specific psychiatric disorders is significantly associated with genes that are differentially expressed and upregulated during differentiation. We find the strongest evidence for schizophrenia, a finding that we replicate in an independent dataset. A longitudinal gene cluster involved in synaptic function primarily drives the association with schizophrenia risk.ConclusionsThese findings reveal that in vitro human neuronal differentiation can be used to translate the polygenic architecture of schizophrenia to biologically relevant pathways that can be modeled in an experimental system. Overall, this work emphasizes the use of longitudinal in vitro transcriptomic signatures as a cellular readout and the application to the genetics of complex traits.

The active form of vitamin B6, pyridoxal phosphate (PLP), is essential for human metabolism. The brain is dependent on vitamin B6 for its neurotransmitter balance. To obtain insight into the genetic determinants of vitamin B6 homeostasis, we conducted a genome-wide association study (GWAS) of the B6 vitamers pyridoxal (PL), PLP and the degradation product of vitamin B6, pyridoxic acid (PA). We collected a unique sample set of cerebrospinal fluid (CSF) and plasma from the same healthy human subjects of Dutch ancestry (n = 493) and included concentrations and ratios in and between these body fluids in our analysis. Based on a multivariate joint analysis of all B6 vitamers and their ratios, we identified a genome-wide significant association at a locus on chromosome 1 containing the ALPL (alkaline phosphatase) gene (minimal p = 7.89 × 10-10, rs1106357, minor allele frequency (MAF) = 0.46), previously associated with vitamin B6 levels in blood. Subjects homozygous for the minor allele showed a 1.4-times-higher ratio between PLP and PL in plasma, and even a 1.6-times-higher ratio between PLP and PL in CSF than subjects homozygous for the major allele. In addition, we observed a suggestive association with the CSF:plasma ratio of PLP on chromosome 15 (minimal p = 7.93 × 10-7, and MAF = 0.06 for rs28789220). Even though this finding is not reaching genome-wide significance, it highlights the potential of our experimental setup for studying transport and metabolism across the blood⁻CSF barrier. This GWAS of B6 vitamers identifies alkaline phosphatase as a key regulator in human vitamin B6 metabolism in CSF as well as plasma. Furthermore, our results demonstrate the potential of genetic studies of metabolites in plasma and CSF to elucidate biological aspects underlying metabolite generation, transport and degradation.

The active form of vitamin B6, pyridoxal phosphate (PLP), is essential for human metabolism. The brain is dependent on vitamin B6 for its neurotransmitter balance. To obtain insight into the genetic determinants of vitamin B6 homeostasis, we conducted a genome-wide association study (GWAS) of the B6 vitamers pyridoxal (PL), PLP and the degradation product of vitamin B6, pyridoxic acid (PA). We collected a unique sample set of cerebrospinal fluid (CSF) and plasma from the same healthy human subjects of Dutch ancestry (n = 493) and included concentrations and ratios in and between these body fluids in our analysis. Based on a multivariate joint analysis of all B6 vitamers and their ratios, we identified a genome-wide significant association at a locus on chromosome 1 containing the ALPL (alkaline phosphatase) gene (minimal p = 7.89 × 10-10, rs1106357, minor allele frequency (MAF) = 0.46), previously associated with vitamin B6 levels in blood. Subjects homozygous for the minor allele showed a 1.4-times-higher ratio between PLP and PL in plasma, and even a 1.6-times-higher ratio between PLP and PL in CSF than subjects homozygous for the major allele. In addition, we observed a suggestive association with the CSF:plasma ratio of PLP on chromosome 15 (minimal p = 7.93 × 10-7, and MAF = 0.06 for rs28789220). Even though this finding is not reaching genome-wide significance, it highlights the potential of our experimental setup for studying transport and metabolism across the blood⁻CSF barrier. This GWAS of B6 vitamers identifies alkaline phosphatase as a key regulator in human vitamin B6 metabolism in CSF as well as plasma. Furthermore, our results demonstrate the potential of genetic studies of metabolites in plasma and CSF to elucidate biological aspects underlying metabolite generation, transport and degradation.

The co-occurrence of a copy number variant (CNV) and a functional variant on the other allele may be a relevant genetic mechanism in schizophrenia. We hypothesized that the cumulative burden of such double hits-in particular those composed of a deletion and a coding single-nucleotide variation (SNV)-is increased in patients with schizophrenia. We combined CNV data with coding variants data in 795 patients with schizophrenia and 474 controls. To limit false CNV-detection, only CNVs called by two algorithms were included. CNV-affected genes were subsequently examined for coding SNVs, which we termed "CNV-SNVs." Correcting for total queried sequence, we assessed the CNV-SNV-burden and the combined predicted deleterious effect. We estimated P-values by permutation of the phenotype. We detected 105 CNV-SNVs; 67 in duplicated and 38 in deleted genic sequence. Although the difference in CNV-SNVs rates was not significant, the combined deleteriousness inferred by CNV-SNVs in deleted sequence was almost 4-fold higher in cases compared with controls (nominal P = 0.009). This effect may be driven by a higher number of CNV-SNVs and/or by a higher degree of predicted deleteriousness of CNV-SNVs. No such effect was observed for duplications. We provide early evidence that deletions co-occurring with a functional variant may be relevant, albeit of modest impact, for the genetic etiology of schizophrenia. Large-scale consortium studies are required to validate our findings. Sequence-based analyses would provide the best resolution for detection of CNVs as well as coding variants genome-wide.

Schizophrenia and bipolar disorder are two distinct diagnoses that share symptomology. Understanding the genetic factors contributing to the shared and disorder-specific symptoms will be crucial for improving diagnosis and treatment. In genetic data consisting of 53,555 cases (20,129 bipolar disorder [BD], 33,426 schizophrenia [SCZ]) and 54,065 controls, we identified 114 genome-wide significant loci implicating synaptic and neuronal pathways shared between disorders. Comparing SCZ to BD (23,585 SCZ, 15,270 BD) identified four genomic regions including one with disorder-independent causal variants and potassium ion response genes as contributing to differences in biology between the disorders. Polygenic risk score (PRS) analyses identified several significant correlations within case-only phenotypes including SCZ PRS with psychotic features and age of onset in BD. For the first time, we discover specific loci that distinguish between BD and SCZ and identify polygenic components underlying multiple symptom dimensions. These results point to the utility of genetics to inform symptomology and potential treatment.

The role of the human microbiome in health and disease is increasingly appreciated. We studied the composition of microbial communities present in blood across 192 individuals, including healthy controls and patients with three disorders affecting the brain: schizophrenia, amyotrophic lateral sclerosis, and bipolar disorder. By using high-quality unmapped RNA sequencing reads as candidate microbial reads, we performed profiling of microbial transcripts detected in whole blood. We were able to detect a wide range of bacterial and archaeal phyla in blood. Interestingly, we observed an increased microbial diversity in schizophrenia patients compared to the three other groups. We replicated this finding in an independent schizophrenia case-control cohort. This increased diversity is inversely correlated with estimated cell abundance of a subpopulation of CD8+ memory T cells in healthy controls, supporting a link between microbial products found in blood, immunity and schizophrenia.

Genetic isolated populations have features that may facilitate genetic analyses and can be leveraged to improve power of mapping genes to complex traits. Our aim was to test the extent to which a population with a former history of geographic isolation and religious endogamy, and currently with one of the highest fertility rates in The Netherlands, shows signs of genetic isolation. For this purpose, genome-wide genotype data was collected of 72 unrelated individuals from this population as well as in a sample of 104 random control subjects from The Netherlands. Additional reference data from different populations and population isolates was available through HapMap and the Human Genome Diversity Project. We performed a number of analyses to compare the genetic structure between these populations: we calculated the pairwise genetic distance between populations, examined the extent of identical-by-descent (IBD) sharing and estimated the effective population size. Genetic analysis of this population showed consistent patterns of a population isolate at all levels tested. We confirmed that this population is most closely related to the Dutch control subjects, and detected high levels of IBD sharing and runs of homozygosity at equal or even higher levels than observed in previously described population isolates. The effective population size of this population was estimated to be several orders of magnitude smaller than that of the Dutch control sample. We conclude that the geographic isolation of this population combined with rapid population growth has resulted in a genetic isolate with great potential value for future genetic studies.

The hippocampal formation is a brain structure integrally involved in episodic memory, spatial navigation, cognition and stress responsiveness. Structural abnormalities in hippocampal volume and shape are found in several common neuropsychiatric disorders. To identify the genetic underpinnings of hippocampal structure here we perform a genome-wide association study (GWAS) of 33,536 individuals and discover six independent loci significantly associated with hippocampal volume, four of them novel. Of the novel loci, three lie within genes (ASTN2, DPP4 and MAST4) and one is found 200 kb upstream of SHH. A hippocampal subfield analysis shows that a locus within the MSRB3 gene shows evidence of a localized effect along the dentate gyrus, subiculum, CA1 and fissure. Further, we show that genetic variants associated with decreased hippocampal volume are also associated with increased risk for Alzheimer's disease (rg=-0.155). Our findings suggest novel biological pathways through which human genetic variation influences hippocampal volume and risk for neuropsychiatric illness.

Intracranial volume reflects the maximally attained brain size during development, and remains stable with loss of tissue in late life. It is highly heritable, but the underlying genes remain largely undetermined. In a genome-wide association study of 32,438 adults, we discovered five previously unknown loci for intracranial volume and confirmed two known signals. Four of the loci were also associated with adult human stature, but these remained associated with intracranial volume after adjusting for height. We found a high genetic correlation with child head circumference (ρgenetic = 0.748), which indicates a similar genetic background and allowed us to identify four additional loci through meta-analysis (Ncombined = 37,345). Variants for intracranial volume were also related to childhood and adult cognitive function, and Parkinson's disease, and were enriched near genes involved in growth pathways, including PI3K-AKT signaling. These findings identify the biological underpinnings of intracranial volume and their link to physiological and pathological traits.

Schizophrenia is a common complex disorder with polygenic inheritance. Here we show that by using an approach that compares the individual loads of rare variants in 1,042 schizophrenia cases and 961 controls, schizophrenia cases carry an increased burden of deleterious mutations. At a genome-wide level, our results implicate non-synonymous, splice site as well as stop-altering single-nucleotide variations occurring at minor allele frequency of ≥ 0.01% in the population. In an independent replication sample of 5,585 schizophrenia cases and 8,103 controls of European ancestry we confirm an enrichment in cases of the alleles identified in our study. In addition, the genes implicated by the increased burden of rare coding variants highlight the involvement of neurodevelopment in the aetiology of schizophrenia.

The highly complex structure of the human brain is strongly shaped by genetic influences. Subcortical brain regions form circuits with cortical areas to coordinate movement, learning, memory and motivation, and altered circuits can lead to abnormal behaviour and disease. To investigate how common genetic variants affect the structure of these brain regions, here we conduct genome-wide association studies of the volumes of seven subcortical regions and the intracranial volume derived from magnetic resonance images of 30,717 individuals from 50 cohorts. We identify five novel genetic variants influencing the volumes of the putamen and caudate nucleus. We also find stronger evidence for three loci with previously established influences on hippocampal volume and intracranial volume. These variants show specific volumetric effects on brain structures rather than global effects across structures. The strongest effects were found for the putamen, where a novel intergenic locus with replicable influence on volume (rs945270; P = 1.08 × 10(-33); 0.52% variance explained) showed evidence of altering the expression of the KTN1 gene in both brain and blood tissue. Variants influencing putamen volume clustered near developmental genes that regulate apoptosis, axon guidance and vesicle transport. Identification of these genetic variants provides insight into the causes of variability in human brain development, and may help to determine mechanisms of neuropsychiatric dysfunction.