Meta-analysis of epigenetic aging in schizophrenia reveals multifaceted relationships with age, sex, illness duration, and polygenic risk

BACKGROUND

The study of biological age acceleration may help identify at-risk individuals and reduce the rising global burden of age-related diseases. Using DNA methylation (DNAm) clocks, we investigated biological aging in schizophrenia (SCZ), a mental illness that is associated with an increased prevalence of age-related disabilities and morbidities. In a whole blood DNAm sample of 1090 SCZ cases and 1206 controls across four European cohorts, we performed a meta-analysis of differential aging using three DNAm clocks (i.e., Hannum, Horvath, and Levine). To dissect how DNAm aging contributes to SCZ, we integrated information on duration of illness and SCZ polygenic risk, as well as stratified our analyses by chronological age and biological sex.

RESULTS

We found that blood-based DNAm aging is significantly altered in SCZ independent from duration of the illness since onset. We observed sex-specific and nonlinear age effects that differed between clocks and point to possible distinct age windows of altered aging in SCZ. Most notably, intrinsic cellular age (Horvath clock) is decelerated in SCZ cases in young adulthood, while phenotypic age (Levine clock) is accelerated in later adulthood compared to controls. Accelerated phenotypic aging was most pronounced in women with SCZ carrying a high polygenic burden with an age acceleration of + 3.82 years (CI 2.02-5.61, P = 1.1E-03). Phenotypic aging and SCZ polygenic risk contributed additively to the illness and together explained up to 14.38% of the variance in disease status.

CONCLUSIONS

Our study contributes to the growing body of evidence of altered DNAm aging in SCZ and points to intrinsic age deceleration in younger adulthood and phenotypic age acceleration in later adulthood in SCZ. Since increased phenotypic age is associated with increased risk of all-cause mortality, our findings indicate that specific and identifiable patient groups are at increased mortality risk as measured by the Levine clock. Our study did not find that DNAm aging could be explained by the duration of illness of patients, but we did observe age- and sex-specific effects that warrant further investigation. Finally, our results show that combining genetic and epigenetic predictors can improve predictions of disease outcomes and may help with disease management in schizophrenia.

Accelerated brain aging as a biomarker for staging in bipolar disorder: an exploratory study

BACKGROUND

Two established staging models outline the longitudinal progression in bipolar disorder (BD) based on episode recurrence or inter-episodic functioning. However, underlying neurobiological mechanisms and corresponding biomarkers remain unexplored. This study aimed to investigate if global and (sub)cortical brain structures, along with brain-predicted age difference (brain-PAD) reflect illness progression as conceptualized in these staging models, potentially identifying brain-PAD as a biomarker for BD staging.

METHODS

In total, 199 subjects with bipolar-I-disorder and 226 control subjects from the Dutch Bipolar Cohort with a high-quality T1-weighted magnetic resonance imaging scan were analyzed. Global and (sub)cortical brain measures and brain-PAD (the difference between biological and chronological age) were estimated. Associations between individual brain measures and the stages of both staging models were explored.

RESULTS

A higher brain-PAD (higher biological age than chronological age) correlated with an increased likelihood of being in a higher stage of the inter-episodic functioning model, but not in the model based on number of mood episodes. However, after correcting for the confounding factors lithium-use and comorbid anxiety, the association lost significance. Global and (sub)cortical brain measures showed no significant association with the stages.

CONCLUSIONS

These results suggest that brain-PAD may be associated with illness progression as defined by impaired inter-episodic functioning. Nevertheless, the significance of this association changed after considering lithium-use and comorbid anxiety disorders. Further research is required to disentangle the intricate relationship between brain-PAD, illness stages, and lithium intake or anxiety disorders. This study provides a foundation for potentially using brain-PAD as a biomarker for illness progression.

Polygenic Risk Associations with Clinical Characteristics and Recurrence of Dupuytren Disease

BACKGROUND

Dupuytren disease (DD) is a common complex trait, with varying severity and incompletely understood cause. Genome-wide association studies (GWAS) have identified risk loci. In this article, we examine whether genetic risk profiles of DD in patients are associated with clinical variation and disease severity and with patient genetic risk profiles of genetically correlated traits, including body mass index (BMI), triglycerides, high-density lipoproteins, type 2 diabetes mellitus, and endophenotypes fasting glucose and glycated hemoglobin.

METHODS

The authors used a well-characterized cohort of 1461 DD patients with available phenotypic and genetic data. Phenotype data include age at onset, recurrence, and family history of disease. Polygenic risk scores (PRSs) of DD, BMI, triglycerides, high-density lipoprotein, type 2 diabetes, fasting glucose, and hemoglobin A1c using various significance thresholds were calculated with PRSice using the most recent GWAS summary statistics. Control data from LifeLines were used to determine P value cutoffs for PRS generation explaining most variance.

RESULTS

The PRS for DD was significantly associated with a positive family history for DD, age at onset, disease onset before the age of 50, and recurrence. We also found a significant negative correlation between the PRSs for DD and BMI.

CONCLUSIONS

Although GWAS studies of DD are designed to identify genetic risk factors distinguishing case/control status, we show that the genetic risk profile for DD also explains part of its clinical variation and disease severity. The PRS may therefore aid in accurate prognostication, choosing initial treatment and in personalized medicine in the future.

CLINICAL QUESTION/LEVEL OF EVIDENCE

Risk, III.

Fine-mapping genomic loci refines bipolar disorder risk genes

Bipolar disorder (BD) is a heritable mental illness with complex etiology. While the largest published genome-wide association study identified 64 BD risk loci, the causal SNPs and genes within these loci remain unknown. We applied a suite of statistical and functional fine-mapping methods to these loci, and prioritized 22 likely causal SNPs for BD. We mapped these SNPs to genes, and investigated their likely functional consequences by integrating variant annotations, brain cell-type epigenomic annotations, brain quantitative trait loci, and results from rare variant exome sequencing in BD. Convergent lines of evidence supported the roles of SCN2A, TRANK1, DCLK3, INSYN2B, SYNE1, THSD7A, CACNA1B, TUBBP5, PLCB3, PRDX5, KCNK4, AP001453.3, TRPT1, FKBP2, DNAJC4, RASGRP1, FURIN, FES, YWHAE, DPH1, GSDMB, MED24, THRA, EEF1A2, and KCNQ2 in BD. These represent promising candidates for functional experiments to understand biological mechanisms and therapeutic potential. Additionally, we demonstrated that fine-mapping effect sizes can improve performance and transferability of BD polygenic risk scores across ancestrally diverse populations, and present a high-throughput fine-mapping pipeline (https://github.com/mkoromina/SAFFARI).

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Childhood abuse v. neglect and risk for major psychiatric disorders

BACKGROUND

Childhood maltreatment (CM) is a strong risk factor for psychiatric disorders but serves in its current definitions as an umbrella for various fundamentally different childhood experiences. As first step toward a more refined analysis of the impact of CM, our objective is to revisit the relation of abuse and neglect, major subtypes of CM, with symptoms across disorders.

METHODS

Three longitudinal studies of major depressive disorder (MDD, N = 1240), bipolar disorder (BD, N = 1339), and schizophrenia (SCZ, N = 577), each including controls (N = 881), were analyzed. Multivariate regression models were used to examine the relation between exposure to abuse, neglect, or their combination to the odds for MDD, BD, SCZ, and symptoms across disorders. Bidirectional Mendelian randomization (MR) was used to probe causality, using genetic instruments of abuse and neglect derived from UK Biobank data (N = 143 473).

RESULTS

Abuse was the stronger risk factor for SCZ (OR 3.51, 95% CI 2.17-5.67) and neglect for BD (OR 2.69, 95% CI 2.09-3.46). Combined CM was related to increased risk exceeding additive effects of abuse and neglect for MDD (RERI = 1.4) and BD (RERI = 1.1). Across disorders, abuse was associated with hallucinations (OR 2.16, 95% CI 1.55-3.01) and suicide attempts (OR 2.16, 95% CI 1.55-3.01) whereas neglect was associated with agitation (OR 1.24, 95% CI 1.02-1.51) and reduced need for sleep (OR 1.64, 95% CI 1.08-2.48). MR analyses were consistent with a bidirectional causal effect of abuse with SCZ (IVWforward = 0.13, 95% CI 0.01-0.24).

CONCLUSIONS

Childhood abuse and neglect are associated with different risks to psychiatric symptoms and disorders. Unraveling the origin of these differences may advance understanding of disease etiology and ultimately facilitate development of improved personalized treatment strategies.

GWAS Meta-Analysis of Suicide Attempt: Identification of 12 Genome-Wide Significant Loci and Implication of Genetic Risks for Specific Health Factors

OBJECTIVE

Suicidal behavior is heritable and is a major cause of death worldwide. Two large-scale genome-wide association studies (GWASs) recently discovered and cross-validated genome-wide significant (GWS) loci for suicide attempt (SA). The present study leveraged the genetic cohorts from both studies to conduct the largest GWAS meta-analysis of SA to date. Multi-ancestry and admixture-specific meta-analyses were conducted within groups of significant African, East Asian, and European ancestry admixtures.

METHODS

This study comprised 22 cohorts, including 43,871 SA cases and 915,025 ancestry-matched controls. Analytical methods across multi-ancestry and individual ancestry admixtures included inverse variance-weighted fixed-effects meta-analyses, followed by gene, gene-set, tissue-set, and drug-target enrichment, as well as summary-data-based Mendelian randomization with brain expression quantitative trait loci data, phenome-wide genetic correlation, and genetic causal proportion analyses.

RESULTS

Multi-ancestry and European ancestry admixture GWAS meta-analyses identified 12 risk loci at p values <5×10-8. These loci were mostly intergenic and implicated DRD2, SLC6A9, FURIN, NLGN1, SOX5, PDE4B, and CACNG2. The multi-ancestry SNP-based heritability estimate of SA was 5.7% on the liability scale (SE=0.003, p=5.7×10-80). Significant brain tissue gene expression and drug set enrichment were observed. There was shared genetic variation of SA with attention deficit hyperactivity disorder, smoking, and risk tolerance after conditioning SA on both major depressive disorder and posttraumatic stress disorder. Genetic causal proportion analyses implicated shared genetic risk for specific health factors.

CONCLUSIONS

This multi-ancestry analysis of suicide attempt identified several loci contributing to risk and establishes significant shared genetic covariation with clinical phenotypes. These findings provide insight into genetic factors associated with suicide attempt across ancestry admixture populations, in veteran and civilian populations, and in attempt versus death.

Quantitative trait loci mapping of circulating metabolites in cerebrospinal fluid to uncover biological mechanisms involved in brain-related phenotypes

Genomic studies of molecular traits have provided mechanistic insights into complex disease, though these lag behind for brain-related traits due to the inaccessibility of brain tissue. We leveraged cerebrospinal fluid (CSF) to study neurobiological mechanisms in vivo , measuring 5,543 CSF metabolites, the largest panel in CSF to date, in 977 individuals of European ancestry. Individuals originated from two separate cohorts including cognitively healthy subjects (n=490) and a well-characterized memory clinic sample, the Amsterdam Dementia Cohort (ADC, n=487). We performed metabolite quantitative trait loci (mQTL) mapping on CSF metabolomics and found 126 significant mQTLs, representing 65 unique CSF metabolites across 51 independent loci. To better understand the role of CSF mQTLs in brain-related disorders, we performed a metabolome-wide association study (MWAS), identifying 40 associations between CSF metabolites and brain traits. Similarly, over 90% of significant mQTLs demonstrated colocalized associations with brain-specific gene expression, unveiling potential neurobiological pathways.

Large-Scale Whole Genome Sequence Analysis of >22,000 Subjects Provides no Evidence of FMR1 Premutation Allele Involvement in Autism Spectrum Disorder

Expansion of a CGG repeat in the Fragile X Messenger Ribonucleoprotein 1 (FMR1) gene on the X chromosome is the cause of Fragile X Syndrome (FXS). The repeat length of unaffected individuals varies between 5-40 repeats, whereas >200 repeats are observed in cases of FXS. The intermediate range between 55-200 repeats is considered the premutation range and is observed in roughly 1:300 females and 1:900 males in the general population. With the availability of large-scale whole genome sequence (WGS) data and the development of computational tools to detect repeat expansions, we systematically examined the role of FMR1 premutation alleles in autism spectrum disorder (ASD) susceptibility, assess the prevalence, and consider the allelic stability between parents and offspring. We analyzed the WGS data of 22,053 subjects, including 32 FXS positive controls, 1359 population controls, and 5467 ASD families. We observed no FMR1 full mutation range repeats among the ASD parent-offspring families but identified 180 family members with premutation range alleles, which represents a higher prevalence compared to the independent WGS control sample and previous reports in the literature. A sex-specific analysis between probands and unaffected siblings did not reveal a significant increase in the burden of premutation alleles in either males or females with ASD. PCR validation, however, suggests an overestimation of the frequency of FMR1 premutation range alleles through computational analysis of WGS data. Overall, we show the utility of large-scale repeat expansion screening in WGS data and conclude that there is no apparent evidence of FMR1 premutation alleles contributing to ASD susceptibility.

Mega-analysis of association between obesity and cortical morphology in bipolar disorders: ENIGMA study in 2832 participants

BACKGROUND

Obesity is highly prevalent and disabling, especially in individuals with severe mental illness including bipolar disorders (BD). The brain is a target organ for both obesity and BD. Yet, we do not understand how cortical brain alterations in BD and obesity interact.

METHODS

We obtained body mass index (BMI) and MRI-derived regional cortical thickness, surface area from 1231 BD and 1601 control individuals from 13 countries within the ENIGMA-BD Working Group. We jointly modeled the statistical effects of BD and BMI on brain structure using mixed effects and tested for interaction and mediation. We also investigated the impact of medications on the BMI-related associations.

RESULTS

BMI and BD additively impacted the structure of many of the same brain regions. Both BMI and BD were negatively associated with cortical thickness, but not surface area. In most regions the number of jointly used psychiatric medication classes remained associated with lower cortical thickness when controlling for BMI. In a single region, fusiform gyrus, about a third of the negative association between number of jointly used psychiatric medications and cortical thickness was mediated by association between the number of medications and higher BMI.

CONCLUSIONS

We confirmed consistent associations between higher BMI and lower cortical thickness, but not surface area, across the cerebral mantle, in regions which were also associated with BD. Higher BMI in people with BD indicated more pronounced brain alterations. BMI is important for understanding the neuroanatomical changes in BD and the effects of psychiatric medications on the brain.

Significant variation in the performance of DNA methylation predictors across data preprocessing and normalization strategies

Background

DNA methylation (DNAm)-based predictors hold great promise to serve as clinical tools for health interventions and disease management. While these algorithms often have high prediction accuracy, the consistency of their performance remains to be determined. We therefore conduct a systematic evaluation across 101 different DNAm data preprocessing and normalization strategies and assess how each analytical strategy affects the consistency of 41 DNAm-based predictors.

Results

Our analyses are conducted in a large EPIC DNAm array dataset from the Jackson Heart Study (N = 2053) that included 146 pairs of technical replicate samples. By estimating the average absolute agreement between replicate pairs, we show that 32 out of 41 predictors (78%) demonstrate excellent consistency when appropriate data processing and normalization steps are implemented. Across all pairs of predictors, we find a moderate correlation in performance across analytical strategies (mean rho = 0.40, SD = 0.27), highlighting significant heterogeneity in performance across algorithms. Successful or unsuccessful removal of technical variation furthermore significantly impacts downstream phenotypic association analysis, such as all-cause mortality risk associations.

Conclusions

We show that DNAm-based algorithms are sensitive to technical variation. The right choice of data processing strategy is important to achieve reproducible estimates and improve prediction accuracy in downstream phenotypic association analyses. For each of the 41 DNAm predictors, we report its degree of consistency and provide the best performing analytical strategy as a guideline for the research community. As DNAm-based predictors become more and more widely used, our work helps improve their performance and standardize their implementation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02793-w.

Common Genetic Variation and Age of Onset of Anorexia Nervosa

Background

Genetics and biology may influence the age of onset of anorexia nervosa (AN). The aims of this study were to determine whether common genetic variation contributes to age of onset of AN and to investigate the genetic associations between age of onset of AN and age at menarche.

Methods

A secondary analysis of the Psychiatric Genomics Consortium genome-wide association study (GWAS) of AN was performed, which included 9335 cases and 31,981 screened controls, all from European ancestries. We conducted GWASs of age of onset, early-onset AN (<13 years), and typical-onset AN, and genetic correlation, genetic risk score, and Mendelian randomization analyses.

Results

Two loci were genome-wide significant in the typical-onset AN GWAS. Heritability estimates (single nucleotide polymorphism-h 2) were 0.01-0.04 for age of onset, 0.16-0.25 for early-onset AN, and 0.17-0.25 for typical-onset AN. Early- and typical-onset AN showed distinct genetic correlation patterns with putative risk factors for AN. Specifically, early-onset AN was significantly genetically correlated with younger age at menarche, and typical-onset AN was significantly negatively genetically correlated with anthropometric traits. Genetic risk scores for age of onset and early-onset AN estimated from independent GWASs significantly predicted age of onset. Mendelian randomization analysis suggested a causal link between younger age at menarche and early-onset AN.

Conclusions

Our results provide evidence consistent with a common variant genetic basis for age of onset and implicate biological pathways regulating menarche and reproduction.

Diagnosis of bipolar disorders and body mass index predict clustering based on similarities in cortical thickness-ENIGMA study in 2436 individuals

AIMS

Rates of obesity have reached epidemic proportions, especially among people with psychiatric disorders. While the effects of obesity on the brain are of major interest in medicine, they remain markedly under-researched in psychiatry.

METHODS

We obtained body mass index (BMI) and magnetic resonance imaging-derived regional cortical thickness, surface area from 836 bipolar disorders (BD) and 1600 control individuals from 14 sites within the ENIGMA-BD Working Group. We identified regionally specific profiles of cortical thickness using K-means clustering and studied clinical characteristics associated with individual cortical profiles.

RESULTS

We detected two clusters based on similarities among participants in cortical thickness. The lower thickness cluster (46.8% of the sample) showed thinner cortex, especially in the frontal and temporal lobes and was associated with diagnosis of BD, higher BMI, and older age. BD individuals in the low thickness cluster were more likely to have the diagnosis of bipolar disorder I and less likely to be treated with lithium. In contrast, clustering based on similarities in the cortical surface area was unrelated to BD or BMI and only tracked age and sex.

CONCLUSIONS

We provide evidence that both BD and obesity are associated with similar alterations in cortical thickness, but not surface area. The fact that obesity increased the chance of having low cortical thickness could explain differences in cortical measures among people with BD. The thinner cortex in individuals with higher BMI, which was additive and similar to the BD-associated alterations, may suggest that treating obesity could lower the extent of cortical thinning in BD.

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

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

Mapping genomic loci implicates genes and synaptic biology in schizophrenia

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

Rare coding variants in ten genes confer substantial risk for schizophrenia

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

Powerful eQTL mapping through low coverage RNA sequencing

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.

Dissecting the Shared Genetic Architecture of Suicide Attempt, Psychiatric Disorders, and Known Risk Factors

BACKGROUND

Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders.

METHODS

We conducted a genome-wide association study (GWAS) of 29,782 suicide attempt (SA) cases and 519,961 controls in the International Suicide Genetics Consortium (ISGC). The GWAS of SA was conditioned on psychiatric disorders using GWAS summary statistics via multitrait-based conditional and joint analysis, to remove genetic effects on SA mediated by psychiatric disorders. We investigated the shared and divergent genetic architectures of SA, psychiatric disorders, and other known risk factors.

RESULTS

Two loci reached genome-wide significance for SA: the major histocompatibility complex and an intergenic locus on chromosome 7, the latter of which remained associated with SA after conditioning on psychiatric disorders and replicated in an independent cohort from the Million Veteran Program. This locus has been implicated in risk-taking behavior, smoking, and insomnia. SA showed strong genetic correlation with psychiatric disorders, particularly major depression, and also with smoking, pain, risk-taking behavior, sleep disturbances, lower educational attainment, reproductive traits, lower socioeconomic status, and poorer general health. After conditioning on psychiatric disorders, the genetic correlations between SA and psychiatric disorders decreased, whereas those with nonpsychiatric traits remained largely unchanged.

CONCLUSIONS

Our results identify a risk locus that contributes more strongly to SA than other phenotypes and suggest a shared underlying biology between SA and known risk factors that is not mediated by psychiatric disorders.

What we learn about bipolar disorder from large-scale neuroimaging: Findings and future directions from the ENIGMA Bipolar Disorder Working Group

MRI-derived brain measures offer a link between genes, the environment and behavior and have been widely studied in bipolar disorder (BD). However, many neuroimaging studies of BD have been underpowered, leading to varied results and uncertainty regarding effects. The Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Bipolar Disorder Working Group was formed in 2012 to empower discoveries, generate consensus findings and inform future hypothesis-driven studies of BD. Through this effort, over 150 researchers from 20 countries and 55 institutions pool data and resources to produce the largest neuroimaging studies of BD ever conducted. The ENIGMA Bipolar Disorder Working Group applies standardized processing and analysis techniques to empower large-scale meta- and mega-analyses of multimodal brain MRI and improve the replicability of studies relating brain variation to clinical and genetic data. Initial BD Working Group studies reveal widespread patterns of lower cortical thickness, subcortical volume and disrupted white matter integrity associated with BD. Findings also include mapping brain alterations of common medications like lithium, symptom patterns and clinical risk profiles and have provided further insights into the pathophysiological mechanisms of BD. Here we discuss key findings from the BD working group, its ongoing projects and future directions for large-scale, collaborative studies of mental illness.

Sex-Dependent Shared and Nonshared Genetic Architecture Across Mood and Psychotic Disorders

BACKGROUND

Sex differences in incidence and/or presentation of schizophrenia (SCZ), major depressive disorder (MDD), and bipolar disorder (BIP) are pervasive. Previous evidence for shared genetic risk and sex differences in brain abnormalities across disorders suggest possible shared sex-dependent genetic risk.

METHODS

We conducted the largest to date genome-wide genotype-by-sex (G×S) interaction of risk for these disorders using 85,735 cases (33,403 SCZ, 19,924 BIP, and 32,408 MDD) and 109,946 controls from the PGC (Psychiatric Genomics Consortium) and iPSYCH.

RESULTS

Across disorders, genome-wide significant single nucleotide polymorphism-by-sex interaction was detected for a locus encompassing NKAIN2 (rs117780815, p = 3.2 × 10-8), which interacts with sodium/potassium-transporting ATPase (adenosine triphosphatase) enzymes, implicating neuronal excitability. Three additional loci showed evidence (p < 1 × 10-6) for cross-disorder G×S interaction (rs7302529, p = 1.6 × 10-7; rs73033497, p = 8.8 × 10-7; rs7914279, p = 6.4 × 10-7), implicating various functions. Gene-based analyses identified G×S interaction across disorders (p = 8.97 × 10-7) with transcriptional inhibitor SLTM. Most significant in SCZ was a MOCOS gene locus (rs11665282, p = 1.5 × 10-7), implicating vascular endothelial cells. Secondary analysis of the PGC-SCZ dataset detected an interaction (rs13265509, p = 1.1 × 10-7) in a locus containing IDO2, a kynurenine pathway enzyme with immunoregulatory functions implicated in SCZ, BIP, and MDD. Pathway enrichment analysis detected significant G×S interaction of genes regulating vascular endothelial growth factor receptor signaling in MDD (false discovery rate-corrected p < .05).

CONCLUSIONS

In the largest genome-wide G×S analysis of mood and psychotic disorders to date, there was substantial genetic overlap between the sexes. However, significant sex-dependent effects were enriched for genes related to neuronal development and immune and vascular functions across and within SCZ, BIP, and MDD at the variant, gene, and pathway levels.

A summary-statistics-based approach to examine the role of serotonin transporter promoter tandem repeat polymorphism in psychiatric phenotypes

In genetic studies of psychiatric disorders in the pre-genome-wide association study (GWAS) era, one of the most commonly studied loci is the serotonin transporter (SLC6A4) promoter polymorphism, a 43-base-pair insertion/deletion polymorphism in the promoter region (5-HTTLPR). The genetic association signals between 5-HTTLPR and psychiatric phenotypes, however, have been inconsistent across many studies. Since the polymorphism cannot be tested via available SNP arrays, we had previously proposed an efficient machine learning algorithm to predict the genotypes of 5-HTTLPR based on the genotypes of eight nearby SNPs, which requires access to individual-level genotype and phenotype data. To utilize the advantage of publicly available GWAS summary statistics obtained from studies with very large sample sizes, we develop a GWAS summary-statistics-based approach for testing the variable number of tandem repeat (VNTR) associations with various phenotypes. We first cross-verify the accuracy of the summary-statistics-based approach for 61 phenotypes in the UK Biobank. Since we observed a strong similarity between the predicted individual-level 5-HTTLPR genotype-based approach and the summary-statistics-based approach, we applied our method to the available neurobehavioral GWAS summary statistics data obtained from large-scale GWAS. We found no genome-wide significant evidence for association between 5-HTTLPR and any of the neurobehavioral traits. We did observe, however, genome-wide significant evidence for association between this locus and human adult height, BMI, and total cholesterol. Our summary-statistics-based approach provides a systematic way to examine the role of VNTRs and related types of genetic polymorphisms in disease risk and trait susceptibility of phenotypes for which large-scale GWAS summary statistics data are available.

Characterisation of age and polarity at onset in bipolar disorder

BACKGROUND

Studying 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.

AIMS

To examine the genetic architecture of AAO and PAO and their association with bipolar disorder disease characteristics.

METHOD

Genome-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.

RESULTS

Earlier 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.

CONCLUSIONS

AAO 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.

Common and rare variant association analyses in amyotrophic lateral sclerosis identify 15 risk loci with distinct genetic architectures and neuron-specific biology

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with a lifetime risk of one in 350 people and an unmet need for disease-modifying therapies. We conducted a cross-ancestry genome-wide association study (GWAS) including 29,612 patients with ALS and 122,656 controls, which identified 15 risk loci. When combined with 8,953 individuals with whole-genome sequencing (6,538 patients, 2,415 controls) and a large cortex-derived expression quantitative trait locus (eQTL) dataset (MetaBrain), analyses revealed locus-specific genetic architectures in which we prioritized genes either through rare variants, short tandem repeats or regulatory effects. ALS-associated risk loci were shared with multiple traits within the neurodegenerative spectrum but with distinct enrichment patterns across brain regions and cell types. Of the environmental and lifestyle risk factors obtained from the literature, Mendelian randomization analyses indicated a causal role for high cholesterol levels. The combination of all ALS-associated signals reveals a role for perturbations in vesicle-mediated transport and autophagy and provides evidence for cell-autonomous disease initiation in glutamatergic neurons.

Bipolar episodes after reproductive events in women with bipolar I disorder, A study of 919 pregnancies

BACKGROUND

Women with bipolar I disorder are at high risk for severe episodes after childbirth, but there is no study that provides an overview on bipolar episode risk both during pregnancy and after childbirth, miscarriage and induced abortion. The aim of this study was to determine the episode risk during all pregnancy outcomes subdivided by first and subsequent pregnancies.

METHODS

Participants were 436 women with bipolar I disorder from the Dutch Bipolar Cohort, having 919 pregnancies of which 762 resulted in a live childbirth, 118 ended in a miscarriage and 39 ended in induced abortion. Women reported on the occurrence of manic or depressed episodes during the perinatal period. Information about medication use was obtained by questionnaires.

RESULTS

Episode risk was 5.2% during pregnancy, and 30.1% in the postpartum period, with a peak in the early postpartum period. Risk of an episode was highest after live birth (34.4%), and lower after miscarriage (15.2%) and induced abortion (27.8%). Women with an episode during pregnancy or postpartum were less likely to have a second child compared to women with an uneventful first pregnancy (cOR=0.34; 95%CI: 0.22-0.51; p<0.001); if they had a second child their risk of an episode was significantly elevated with a subsequent pregnancy (cOR=6.17; 95%CI: 3.64-10.45; p<0.001).

LIMITATIONS

Retrospective cross-sectional design with assessment (partial) through self-report in a homogeneous population.

CONCLUSIONS

Women with bipolar I disorder have a six times higher risk of an episode after delivery compared to during pregnancy, therefore preventive strategies are particularly important immediately after delivery.

Association between body mass index and subcortical brain volumes in bipolar disorders-ENIGMA study in 2735 individuals

Individuals with bipolar disorders (BD) frequently suffer from obesity, which is often associated with neurostructural alterations. Yet, the effects of obesity on brain structure in BD are under-researched. We obtained MRI-derived brain subcortical volumes and body mass index (BMI) from 1134 BD and 1601 control individuals from 17 independent research sites within the ENIGMA-BD Working Group. We jointly modeled the effects of BD and BMI on subcortical volumes using mixed-effects modeling and tested for mediation of group differences by obesity using nonparametric bootstrapping. All models controlled for age, sex, hemisphere, total intracranial volume, and data collection site. Relative to controls, individuals with BD had significantly higher BMI, larger lateral ventricular volume, and smaller volumes of amygdala, hippocampus, pallidum, caudate, and thalamus. BMI was positively associated with ventricular and amygdala and negatively with pallidal volumes. When analyzed jointly, both BD and BMI remained associated with volumes of lateral ventricles  and amygdala. Adjusting for BMI decreased the BD vs control differences in ventricular volume. Specifically, 18.41% of the association between BD and ventricular volume was mediated by BMI (Z = 2.73, p = 0.006). BMI was associated with similar regional brain volumes as BD, including lateral ventricles, amygdala, and pallidum. Higher BMI may in part account for larger ventricles, one of the most replicated findings in BD. Comorbidity with obesity could explain why neurostructural alterations are more pronounced in some individuals with BD. Future prospective brain imaging studies should investigate whether obesity could be a modifiable risk factor for neuroprogression.

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

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

Accelerated aging in the brain, epigenetic aging in blood, and polygenic risk for schizophrenia

Schizophrenia patients show signs of accelerated aging in cognitive and physiological domains. Both schizophrenia and accelerated aging, as measured by MRI brain images and epigenetic clocks, are correlated with increased mortality. However, the association between these aging measures have not yet been studied in schizophrenia patients. In schizophrenia patients and healthy subjects, accelerated aging was assessed in brain tissue using a longitudinal MRI (N = 715 scans; mean scan interval 3.4 year) and in blood using two epigenetic age clocks (N = 172). Differences ('gaps') between estimated ages and chronological ages were calculated, as well as the acceleration rate of brain aging. The correlations between these aging measures as well as with polygenic risk scores for schizophrenia (PRS; N = 394) were investigated. Brain aging and epigenetic aging were not significantly correlated. Polygenic risk for schizophrenia was significantly correlated with brain age gap, brain age acceleration rate, and negatively correlated with DNAmAge gap, but not with PhenoAge gap. However, after controlling for disease status and multiple comparisons correction, these effects were no longer significant. Our results imply that the (accelerated) aging observed in the brain and blood reflect distinct biological processes. Our findings will require replication in a larger cohort.

Genome-wide analyses of smoking behaviors in schizophrenia: Findings from the Psychiatric Genomics Consortium

While 17% of US adults use tobacco regularly, smoking rates among persons with schizophrenia are upwards of 60%. Research supports a shared etiological basis for smoking and schizophrenia, including findings from genome-wide association studies (GWAS). However, few studies have directly tested whether the same or distinct genetic variants also influence smoking behavior among schizophrenia cases. Using data from the Psychiatric Genomics Consortium (PGC) study of schizophrenia (35476 cases, 46839 controls), we estimated genetic correlations between these traits and tested whether polygenic risk scores (PRS) constructed from the results of smoking behaviors GWAS were associated with schizophrenia risk or smoking behaviors among schizophrenia cases. Results indicated significant genetic correlations of schizophrenia with smoking initiation (rg = 0.159; P = 5.05 × 10-10), cigarettes-smoked-per-day (rg = 0.094; P = 0.006), and age-of-onset of smoking (rg = 0.10; P = 0.009). Comparing smoking behaviors among schizophrenia cases to the general population, we observe positive genetic correlations for smoking initiation (rg = 0.624, P = 0.002) and cigarettes-smoked-per-day (rg = 0.689, P = 0.120). Similarly, TAG-based PRS for smoking initiation and cigarettes-smoked-per-day were significantly associated with smoking initiation (P = 3.49 × 10-5) and cigarettes-smoked-per-day (P = 0.007) among schizophrenia cases. We performed the first GWAS of smoking behavior among schizophrenia cases and identified a novel association with cigarettes-smoked-per-day upstream of the TMEM106B gene on chromosome 7p21.3 (rs148253479, P = 3.18 × 10-8, n = 3520). Results provide evidence of a partially shared genetic basis for schizophrenia and smoking behaviors. Additionally, genetic risk factors for smoking behaviors were largely shared across schizophrenia and non-schizophrenia populations. Future research should address mechanisms underlying these associations to aid both schizophrenia and smoking treatment and prevention efforts.

Gene Expression Imputation Across Multiple Tissue Types Provides Insight Into the Genetic Architecture of Frontotemporal Dementia and Its Clinical Subtypes

BACKGROUND

The etiology of frontotemporal dementia (FTD) is poorly understood. To identify genes with predicted expression levels associated with FTD, we integrated summary statistics with external reference gene expression data using a transcriptome-wide association study approach.

METHODS

FUSION software was used to leverage FTD summary statistics (all FTD: n = 2154 cases, n = 4308 controls; behavioral variant FTD: n = 1337 cases, n = 2754 controls; semantic dementia: n = 308 cases, n = 616 controls; progressive nonfluent aphasia: n = 269 cases, n = 538 controls; FTD with motor neuron disease: n = 200 cases, n = 400 controls) from the International FTD-Genomics Consortium with 53 expression quantitative loci tissue type panels (n = 12,205; 5 consortia). Significance was assessed using a 5% false discovery rate threshold.

RESULTS

We identified 73 significant gene-tissue associations for FTD, representing 44 unique genes in 34 tissue types. Most significant findings were derived from dorsolateral prefrontal cortex splicing data (n = 19 genes, 26%). The 17q21.31 inversion locus contained 23 significant associations, representing 6 unique genes. Other top hits included SEC22B (a gene involved in vesicle trafficking), TRGV5, and ZNF302. A single gene finding (RAB38) was observed for behavioral variant FTD. For other clinical subtypes, no significant associations were observed.

CONCLUSIONS

We identified novel candidate genes (e.g., SEC22B) and previously reported risk regions (e.g., 17q21.31) for FTD. Most significant associations were observed in dorsolateral prefrontal cortex splicing data despite the modest sample size of this reference panel. This suggests that our findings are specific to FTD and are likely to be biologically relevant highlights of genes at different FTD risk loci that are contributing to the disease pathology.

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

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

Genetic analysis of activity, brain and behavioral associations in extended families with heavy genetic loading for bipolar disorder

BACKGROUND

Disturbed sleep and activity are prominent features of bipolar disorder type I (BP-I). However, the relationship of sleep and activity characteristics to brain structure and behavior in euthymic BP-I patients and their non-BP-I relatives is unknown. Additionally, underlying genetic relationships between these traits have not been investigated.

METHODS

Relationships between sleep and activity phenotypes, assessed using actigraphy, with structural neuroimaging (brain) and cognitive and temperament (behavior) phenotypes were investigated in 558 euthymic individuals from multi-generational pedigrees including at least one member with BP-I. Genetic correlations between actigraphy-brain and actigraphy-behavior associations were assessed, and bivariate linkage analysis was conducted for trait pairs with evidence of shared genetic influences.

RESULTS

More physical activity and longer awake time were significantly associated with increased brain volumes and cortical thickness, better performance on neurocognitive measures of long-term memory and executive function, and less extreme scores on measures of temperament (impulsivity, cyclothymia). These associations did not differ between BP-I patients and their non-BP-I relatives. For nine activity-brain or activity-behavior pairs there was evidence for shared genetic influence (genetic correlations); of these pairs, a suggestive bivariate quantitative trait locus on chromosome 7 for wake duration and verbal working memory was identified.

CONCLUSIONS

Our findings indicate that increased physical activity and more adequate sleep are associated with increased brain size, better cognitive function and more stable temperament in BP-I patients and their non-BP-I relatives. Additionally, we found evidence for pleiotropy of several actigraphy-behavior and actigraphy-brain phenotypes, suggesting a shared genetic basis for these traits.

Genetic and clinical analyses of psychosis spectrum symptoms in a large multiethnic youth cohort reveal significant link with ADHD

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.

Genome-wide association study identifies 48 common genetic variants associated with handedness

Handedness has been extensively studied because of its relationship with language and the over-representation of left-handers in some neurodevelopmental disorders. Using data from the UK Biobank, 23andMe and the International Handedness Consortium, we conducted a genome-wide association meta-analysis of handedness (N = 1,766,671). We found 41 loci associated (P < 5 × 10-8) with left-handedness and 7 associated with ambidexterity. Tissue-enrichment analysis implicated the CNS in the aetiology of handedness. Pathways including regulation of microtubules and brain morphology were also highlighted. We found suggestive positive genetic correlations between left-handedness and neuropsychiatric traits, including schizophrenia and bipolar disorder. Furthermore, the genetic correlation between left-handedness and ambidexterity is low (rG = 0.26), which implies that these traits are largely influenced by different genetic mechanisms. Our findings suggest that handedness is highly polygenic and that the genetic variants that predispose to left-handedness may underlie part of the association with some psychiatric disorders.

Virtual Histology of Cortical Thickness and Shared Neurobiology in 6 Psychiatric Disorders

IMPORTANCE

Large-scale neuroimaging studies have revealed group differences in cortical thickness across many psychiatric disorders. The underlying neurobiology behind these differences is not well understood.

OBJECTIVE

To determine neurobiologic correlates of group differences in cortical thickness between cases and controls in 6 disorders: attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder (BD), major depressive disorder (MDD), obsessive-compulsive disorder (OCD), and schizophrenia.

DESIGN, SETTING, AND PARTICIPANTS

Profiles of group differences in cortical thickness between cases and controls were generated using T1-weighted magnetic resonance images. Similarity between interregional profiles of cell-specific gene expression and those in the group differences in cortical thickness were investigated in each disorder. Next, principal component analysis was used to reveal a shared profile of group difference in thickness across the disorders. Analysis for gene coexpression, clustering, and enrichment for genes associated with these disorders were conducted. Data analysis was conducted between June and December 2019. The analysis included 145 cohorts across 6 psychiatric disorders drawn from the ENIGMA consortium. The numbers of cases and controls in each of the 6 disorders were as follows: ADHD: 1814 and 1602; ASD: 1748 and 1770; BD: 1547 and 3405; MDD: 2658 and 3572; OCD: 2266 and 2007; and schizophrenia: 2688 and 3244.

MAIN OUTCOMES AND MEASURES

Interregional profiles of group difference in cortical thickness between cases and controls.

RESULTS

A total of 12 721 cases and 15 600 controls, ranging from ages 2 to 89 years, were included in this study. Interregional profiles of group differences in cortical thickness for each of the 6 psychiatric disorders were associated with profiles of gene expression specific to pyramidal (CA1) cells, astrocytes (except for BD), and microglia (except for OCD); collectively, gene-expression profiles of the 3 cell types explain between 25% and 54% of variance in interregional profiles of group differences in cortical thickness. Principal component analysis revealed a shared profile of difference in cortical thickness across the 6 disorders (48% variance explained); interregional profile of this principal component 1 was associated with that of the pyramidal-cell gene expression (explaining 56% of interregional variation). Coexpression analyses of these genes revealed 2 clusters: (1) a prenatal cluster enriched with genes involved in neurodevelopmental (axon guidance) processes and (2) a postnatal cluster enriched with genes involved in synaptic activity and plasticity-related processes. These clusters were enriched with genes associated with all 6 psychiatric disorders.

CONCLUSIONS AND RELEVANCE

In this study, shared neurobiologic processes were associated with differences in cortical thickness across multiple psychiatric disorders. These processes implicate a common role of prenatal development and postnatal functioning of the cerebral cortex in these disorders.

Whole blood transcriptome analysis in bipolar disorder reveals strong lithium effect

BACKGROUND

Bipolar disorder (BD) is a highly heritable mood disorder with complex genetic architecture and poorly understood etiology. Previous transcriptomic BD studies have had inconsistent findings due to issues such as small sample sizes and difficulty in adequately accounting for confounders like medication use.

METHODS

We performed a differential expression analysis in a well-characterized BD case-control sample (Nsubjects = 480) by RNA sequencing of whole blood. We further performed co-expression network analysis, functional enrichment, and cell type decomposition, and integrated differentially expressed genes with genetic risk.

RESULTS

While we observed widespread differential gene expression patterns between affected and unaffected individuals, these effects were largely linked to lithium treatment at the time of blood draw (FDR < 0.05, Ngenes = 976) rather than BD diagnosis itself (FDR < 0.05, Ngenes = 6). These lithium-associated genes were enriched for cell signaling and immune response functional annotations, among others, and were associated with neutrophil cell-type proportions, which were elevated in lithium users. Neither genes with altered expression in cases nor in lithium users were enriched for BD, schizophrenia, and depression genetic risk based on information from genome-wide association studies, nor was gene expression associated with polygenic risk scores for BD.

CONCLUSIONS

These findings suggest that BD is associated with minimal changes in whole blood gene expression independent of medication use but emphasize the importance of accounting for medication use and cell type heterogeneity in psychiatric transcriptomic studies. The results of this study add to mounting evidence of lithium's cell signaling and immune-related mechanisms.

Extensions of Multiple-Group Item Response Theory Alignment: Application to Psychiatric Phenotypes in an International Genomics Consortium

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.

Lithium Use during Pregnancy and the Risk of Miscarriage

Recent studies have provided new data on the teratogenicity of lithium. Less is known about the risk of miscarriage after lithium use during pregnancy. The aim of this study was to investigate the association between lithium use during pregnancy and miscarriage. Participants were women with bipolar I disorder and one or more pregnancies, of which information on medication use and pregnancy outcome was available (n = 443). The unadjusted odds ratios for miscarriage after lithium use during pregnancy was calculated. Multilevel logistic regression was used to calculate the odds ratio, adjusted for the age at conception and the clustering of pregnancies per woman. Miscarriages occurred in 20.8% of the lithium-exposed pregnancies (16/77), compared with 10.9% of the unexposed pregnancies (40/366) (OR = 2.14; 95% CI: 1.13-4.06). The adjusted odds ratio of miscarriage after lithium use during pregnancy was 2.94 (95% CI: 1.39-6.22). Lithium use during pregnancy may increase the risk of miscarriage.

Complement genes contribute sex-biased vulnerability in diverse disorders

Many common illnesses differentially affect men and women for unknown reasons. The autoimmune diseases lupus and Sjögren’s syndrome affect nine times more women than men¹, whereas schizophrenia affects men more frequently and severely². All three illnesses have their strongest common genetic associations in the Major Histocompatibility Complex (MHC) locus, an association that in lupus and Sjögren’s syndrome has long been thought to arise from alleles of the human leukocyte antigen (HLA) genes at that locus^3–6. Here we show that the complement component 4 (C4) genes, which are also in the MHC locus and were recently found to increase risk for schizophrenia⁷, generate 7-fold variation in risk for lupus (95% CI: 5.88–8.61; p < 10⁻¹¹⁷ in total) and 16-fold variation in risk for Sjögren’s syndrome (95% CI: 8.59–30.89; p < 10⁻²³ in total) 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 reduced risk for lupus and Sjögren’s syndrome. In all three illnesses, C4 alleles acted more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for lupus, 31-fold variation in risk for Sjögren’s syndrome, and 1.7-fold variation in schizophrenia risk among men (vs. 6-fold, 15-fold, and 1.26-fold among women respectively). At a protein level, both C4 and its effector C3 were present at greater levels in men than women in cerebrospinal fluid (p < 10⁻⁵ for both C4 and C3) and plasma^8,9 among adults ages 20–50, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help explain the larger effects of C4 alleles in men, women’s greater risk of SLE and Sjögren’s, and men’s greater vulnerability in schizophrenia. These results implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

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

Importance

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

Objective

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

Design, Setting, and Participants

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

Main Outcomes and Measures

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

Results

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

Conclusions and Relevance

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

The genetic architecture of the human cerebral cortex

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

Transcriptomic abnormalities in peripheral blood in bipolar disorder, and discrimination of the major psychoses

We performed a transcriptome-wide meta-analysis and gene co-expression network analysis to identify genes and gene networks dysregulated in the peripheral blood of bipolar disorder (BD) cases relative to unaffected comparison subjects, and determined the specificity of the transcriptomic signatures of BD and schizophrenia (SZ). Nineteen genes and 4 gene modules were significantly differentially expressed in BD cases. Thirteen gene modules were shown to be differentially expressed in a combined case-group of BD and SZ subjects called "major psychosis", including genes biologically linked to apoptosis, reactive oxygen, chromatin remodeling, and immune signaling. No modules were differentially expressed between BD and SZ cases. Machine-learning classifiers trained to separate diagnostic classes based solely on gene expression profiles could distinguish BD cases from unaffected comparison subjects with an area under the curve (AUC) of 0.724, as well as BD cases from SZ cases with AUC = 0.677 in withheld test samples. We introduced a novel and straightforward method called "polytranscript risk scoring" that could distinguish BD cases from unaffected subjects (AUC = 0.672) and SZ cases (AUC = 0.607) significantly better than expected by chance. Taken together, our results highlighted gene expression alterations common to BD and SZ that involve biological processes of inflammation, oxidative stress, apoptosis, and chromatin regulation, and highlight disorder-specific changes in gene expression that discriminate the major psychoses.

Contribution of common and rare variants to bipolar disorder susceptibility in extended pedigrees from population isolates

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.

Exploring the clinical utility of two staging models for bipolar disorder

OBJECTIVE

To assess the clinical utility of two staging models for bipolar disorder by examining distribution, correlation, and the relationship to external criteria. These are primarily defined by the recurrence of mood episodes (model A), or by intra-episodic functioning (model B).

METHODS

In the Dutch Bipolar Cohort, stages according to models A and B were assigned to all patients with bipolar-I-disorder (BD-I; N = 1396). The dispersion of subjects over the stages was assessed and the association between the two models calculated. For both models, change in several clinical markers were concordant with the stage was investigated.

RESULTS

Staging was possible in 87% of subjects for model A and 75% for model B. For model A, 1079 participants (93%) were assigned to stage 3c (recurrent episodes). Subdividing stage 3c with cut-offs at 5 and 10 episodes resulted in subgroups containing 242, 510, and 327 subjects. For model B, most participants were assigned to stage II (intra-episodic symptoms, N = 431 (41%)) or stage III (inability to work, N = 451 (43%)). A low association between models was found. For both models, the clinical markers "age at onset," "treatment resistance," and "episode acceleration" changed concordant with the stages.

CONCLUSION

The majority of patients with BD-I clustered in recurrent stage 3 of Model A. Model B showed a larger dispersion. The stepwise change in several clinical markers supports the construct validity of both models. Combining the two staging models and sub-differentiating the recurrent stage into categories with cut-offs at 5 and 10 lifetime episodes improves the clinical utility of staging for individual patients.

No association between anti-thyroidperoxidase antibodies and bipolar disorder: a study in the Dutch Bipolar Cohort and a meta-analysis

BACKGROUND

Thyroid autoimmunity has been associated with bipolar disorder (BD). However, results from previous studies on the seroprevalence of anti-thyroid peroxidase antibodies (TPO-abs) in BD are inconsistent.

OBJECTIVES

The aim of the present study is to investigate whether the seroprevalence and titer levels of TPO-abs are related to BD.

METHOD

TPO-abs were measured in plasma samples of 760 patients with bipolar disorder, 261 first-degree relatives and 363 controls by enzyme-linked immunosorbent assay (ELISA). To address methodological limitations of previous studies, we assessed clinical characteristics with several (self-reported) questionnaires to investigate whether TPO-abs positivity is related to particular clinical subgroups of BD patients. We performed an additional meta-analysis of seroprevalences of TPO-abs in BD patients including data from present and previous studies.

RESULTS

Seroprevalence or titer levels of TPO-abs did not significantly differ between patients with BD, their first-degree relatives, and controls. In BD patients, the prevalence of TPO-abs was unrelated to specific clinical factors, including lithium use. Our meta-analysis of twelve studies showed an overall odds ratio of 1.3 (CI 95 %: 0.7-2.3; p = 0.30), reaffirming the absence of an association of BD with TPO-abs.

CONCLUSIONS

In the largest study of TPO-abs in BD to date, our findings indicate that TPO-abs are not associated with (the risk for) bipolar disorder.

State-Dependent Functional Dysconnectivity in Youth With Psychosis Spectrum Symptoms

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.

The Relationship Between Polygenic Risk Scores and Cognition in Schizophrenia

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Double hits in schizophrenia

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.

The association between antibodies to neurotropic pathogens and bipolar disorder : A study in the Dutch Bipolar (DB) Cohort and meta-analysis

Exposure to neurotropic pathogens has been hypothesized to be a risk factor for the development of bipolar disorder (BD). However, evidence so far is inconsistent. We, therefore, analyzed the seroprevalence and titer levels of IgG antibodies against several herpesviruses and Toxoplasma gondii (T. gondii) in plasma of 760 patients with a bipolar disorder, 144 first-degree matched relatives and 132 controls of the Dutch Bipolar (DB) Cohort using ELISA. In addition, we performed a literature-based meta-analysis on the seroprevalence of IgG antibodies against these pathogens (n = 14). Our results in the DB Cohort and subsequent meta-analysis (n = 2364 BD patients, n = 5101 controls) show no association between exposure to herpesviruses and bipolar disorder (HSV-1 [adjusted OR 0.842, 95% CI 0.567-1.230], HSV-2 [adjusted OR 0.877, 95% CI 0.437-1.761], CMV [adjusted OR 0.884 95% CI 0.603-1.295], EBV [adjusted OR 0.968 95% CI 0.658-1.423]). In the DB Cohort, we did not find an association between bipolar disorder and T. gondii titer or seroprevalence either [adjusted OR 1.018, 95% CI 0.672-1.542]. The overall OR was not significant for T. gondii [OR: 1.4, 95% CI 0.95-1.90, p = 0.09), but subgroup analyses in age groups below 40 years showed a significantly increased seroprevalence of T. gondii IgGs in BD [OR: 1.8 (95% CI 1.10-2.89, p = 0.021]. Our meta-analysis indicates that T. gondii exposure may be a risk factor for BD in certain subpopulations.