Monozygotic twins discordant for schizophrenia differ in maturation and synaptic transmission

Schizophrenia affects approximately 1% of the world population. Genetics, epigenetics, and environmental factors are known to play a role in this psychiatric disorder. While there is a high concordance in monozygotic twins, about half of twin pairs are discordant for schizophrenia. To address the question of how and when concordance in monozygotic twins occur, we have obtained fibroblasts from two pairs of schizophrenia discordant twins (one sibling with schizophrenia while the second one is unaffected by schizophrenia) and three pairs of healthy twins (both of the siblings are healthy). We have prepared iPSC models for these 3 groups of patients with schizophrenia, unaffected co-twins, and the healthy twins. When the study started the co-twins were considered healthy and unaffected but both the co-twins were later diagnosed with a depressive disorder. The reprogrammed iPSCs were differentiated into hippocampal neurons to measure the neurophysiological abnormalities in the patients. We found that the neurons derived from the schizophrenia patients were less arborized, were hypoexcitable with immature spike features, and exhibited a significant reduction in synaptic activity with dysregulation in synapse-related genes. Interestingly, the neurons derived from the co-twin siblings who did not have schizophrenia formed another distinct group that was different from the neurons in the group of the affected twin siblings but also different from the neurons in the group of the control twins. Importantly, their synaptic activity was not affected. Our measurements that were obtained from schizophrenia patients and their monozygotic twin and compared also to control healthy twins point to hippocampal synaptic deficits as a central mechanism in schizophrenia.

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

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.

DNA methylation variation after a parenting program for child conduct problems: Findings from a randomized controlled trial

This study investigated associations of the Incredible Years (IY) parenting program with children's DNA methylation. Participants were 289 Dutch children aged 3-9 years (75% European ancestry, 48% female) with above-average conduct problems. Saliva was collected 2.5 years after families were randomized to IY or care as usual (CAU). Using an intention-to-treat approach, confirmatory multiple-regression analyses revealed no significant differences between the IY and CAU groups in children's methylation levels at the NR3C1 and FKBP5 genes. However, exploratory epigenome-wide analyses revealed nine differentially methylated regions between groups, coinciding with SLAMF1, MITF, FAM200B, PSD3, SNX31, and CELSR1. The study provides preliminary evidence for associations of IY with children's salivary methylation levels and highlights the need for further research into biological outcomes of parenting programs.

Cell-type deconvolution of bulk-blood RNA-seq reveals biological insights into neuropsychiatric disorders

Genome-wide association studies (GWASs) have uncovered susceptibility loci associated with psychiatric disorders such as bipolar disorder (BP) and schizophrenia (SCZ). However, most of these loci are in non-coding regions of the genome, and the causal mechanisms of the link between genetic variation and disease risk is unknown. Expression quantitative trait locus (eQTL) analysis of bulk tissue is a common approach used for deciphering underlying mechanisms, although this can obscure cell-type-specific signals and thus mask trait-relevant mechanisms. Although single-cell sequencing can be prohibitively expensive in large cohorts, computationally inferred cell-type proportions and cell-type gene expression estimates have the potential to overcome these problems and advance mechanistic studies. Using bulk RNA-seq from 1,730 samples derived from whole blood in a cohort ascertained from individuals with BP and SCZ, this study estimated cell-type proportions and their relation with disease status and medication. For each cell type, we found between 2,875 and 4,629 eGenes (genes with an associated eQTL), including 1,211 that are not found on the basis of bulk expression alone. We performed a colocalization test between cell-type eQTLs and various traits and identified hundreds of associations that occur between cell-type eQTLs and GWASs but that are not detected in bulk eQTLs. Finally, we investigated the effects of lithium use on the regulation of cell-type expression loci and found examples of genes that are differentially regulated according to lithium use. Our study suggests that applying computational methods to large bulk RNA-seq datasets of non-brain tissue can identify disease-relevant, cell-type-specific biology of psychiatric disorders and psychiatric medication.

The Cutting Edge of Epigenetic Clocks: In Search of Mechanisms Linking Aging and Mental Health

Individuals with psychiatric disorders are at increased risk of age-related diseases and early mortality. Recent studies demonstrate that this link between mental health and aging is reflected in epigenetic clocks, aging biomarkers based on DNA methylation. The reported relationships between epigenetic clocks and mental health are mostly correlational, and the mechanisms are poorly understood. Here, we review recent progress concerning the molecular and cellular processes underlying epigenetic clocks as well as novel technologies enabling further studies of the causes and consequences of epigenetic aging. We then review the current literature on how epigenetic clocks relate to specific aspects of mental health, such as stress, medications, substance use, health behaviors, and symptom clusters. We propose an integrated framework where mental health and epigenetic aging are each broken down into multiple distinct processes, which are then linked to each other, using stress and schizophrenia as examples. This framework incorporates the heterogeneity and complexity of both mental health conditions and aging, may help reconcile conflicting results, and provides a basis for further hypothesis-driven research in humans and model systems to investigate potentially causal mechanisms linking aging and mental health.

Negative parenting, epigenetic age, and psychological problems: prospective associations from adolescence to young adulthood

BACKGROUND

Epigenetic clocks are based on DNA methylation levels of several genomic loci and have been developed as indices of biological aging. Studies examining the effects of stressful environmental exposures have shown that stress is associated with differences between epigenetic age and chronological age (i.e., Epigenetic Age acceleration, EA). This pre-registered longitudinal study examined the long-term effects of negative parenting and psychological problems throughout adolescence (ages 13-17 years) on EA in late adolescence (age 17 years) and EA changes from late adolescence to young adulthood (age 25 years). Further, it examined how (change in) EA is related to changes in psychological problems from adolescence to young adulthood.

METHODS

We used data from a sample of 434 participants followed from age 13 to age 25, with saliva collected at ages 17 and 25. We estimated EA using four commonly used epigenetic clocks and analyzed the data using Structural Equation Modeling.

RESULTS

While negative parenting was not related to EA nor change in EA, (change in) EA was related to developmental indices such as externalizing problems and self-concept clarity.

CONCLUSIONS

Declining psychological well-being during young adulthood was preceded by EA.

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.

Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders

Posttraumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 novel). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (e.g., GRIA1, GRM8, CACNA1E ), developmental, axon guidance, and transcription factors (e.g., FOXP2, EFNA5, DCC ), synaptic structure and function genes (e.g., PCLO, NCAM1, PDE4B ), and endocrine or immune regulators (e.g., ESR1, TRAF3, TANK ). Additional top genes influence stress, immune, fear, and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

An international research agenda for clozapine-resistant schizophrenia

Treatment-resistant symptoms occur in about a third of patients with schizophrenia and are associated with a substantial reduction in their quality of life. The development of new treatment options for clozapine-resistant schizophrenia constitutes a crucial, unmet need in psychiatry. Additionally, an overview of past and possible future research avenues to optimise the early detection, diagnosis, and management of clozapine-resistant schizophrenia is unavailable. In this Health Policy, we discuss the ongoing challenges associated with clozapine-resistant schizophrenia faced by patients and health-care providers worldwide to improve the understanding of this condition. We then revisit several clozapine guidelines, the diagnostic tests and treatment options for clozapine-resistant schizophrenia, and currently applied research approaches in clozapine-resistant schizophrenia. We also suggest methodologies and targets for future research, divided into innovative nosology-oriented field trials (eg, examining dimensional symptom staging), translational approaches (eg, genetics), epidemiological research (eg, real-world studies), and interventional studies (eg, non-traditional trial designs incorporating lived experiences and caregivers' perspectives). Finally, we note that low-income and middle-income countries are under-represented in studies on clozapine-resistant schizophrenia and propose an agenda to guide multinational research on the cause and treatment of clozapine-resistant schizophrenia. We hope that this research agenda will empower better global representation of patients living with clozapine-resistant schizophrenia and ultimately improve their functional outcomes and quality of life.

Cell type deconvolution of bulk blood RNA-Seq to reveal biological insights of neuropsychiatric disorders

Genome-wide association studies (GWAS) have uncovered susceptibility loci associated with psychiatric disorders like bipolar disorder (BP) and schizophrenia (SCZ). However, most of these loci are in non-coding regions of the genome with unknown causal mechanisms of the link between genetic variation and disease risk. Expression quantitative trait loci (eQTL) analysis of bulk tissue is a common approach to decipher underlying mechanisms, though this can obscure cell-type specific signals thus masking trait-relevant mechanisms. While single-cell sequencing can be prohibitively expensive in large cohorts, computationally inferred cell type proportions and cell type gene expression estimates have the potential to overcome these problems and advance mechanistic studies. Using bulk RNA-Seq from 1,730 samples derived from whole blood in a cohort ascertained for individuals with BP and SCZ this study estimated cell type proportions and their relation with disease status and medication. We found between 2,875 and 4,629 eGenes for each cell type, including 1,211 eGenes that are not found using bulk expression alone. We performed a colocalization test between cell type eQTLs and various traits and identified hundreds of associations between cell type eQTLs and GWAS loci that are not detected in bulk eQTLs. Finally, we investigated the effects of lithium use on cell type expression regulation and found examples of genes that are differentially regulated dependent on lithium use. Our study suggests that computational methods can be applied to large bulk RNA-Seq datasets of non-brain tissue to identify disease-relevant, cell type specific biology of psychiatric disorders and psychiatric medication.

Independent contribution of polygenic risk for schizophrenia and cannabis use in predicting psychotic-like experiences in young adulthood: testing gene × environment moderation and mediation

BACKGROUND

It has not yet been determined if the commonly reported cannabis-psychosis association is limited to individuals with pre-existing genetic risk for psychotic disorders.

METHODS

We examined whether the relationship between polygenic risk score for schizophrenia (PRS-Sz) and psychotic-like experiences (PLEs), as measured by the Community Assessment of Psychic Experiences-42 (CAPE-42) questionnaire, is mediated or moderated by lifetime cannabis use at 16 years of age in 1740 of the individuals of the European IMAGEN cohort. Secondary analysis examined the relationships between lifetime cannabis use, PRS-Sz and the various sub-scales of the CAPE-42. Sensitivity analyses including covariates, including a PRS for cannabis use, were conducted and results were replicated using data from 1223 individuals in the Dutch Utrecht cannabis cohort.

RESULTS

PRS-Sz significantly predicted cannabis use (p = 0.027) and PLE (p = 0.004) in the IMAGEN cohort. In the full model, considering PRS-Sz and covariates, cannabis use was also significantly associated with PLE in IMAGEN (p = 0.007). Results remained consistent in the Utrecht cohort and through sensitivity analyses. Nevertheless, there was no evidence of a mediation or moderation effects.

CONCLUSIONS

These results suggest that cannabis use remains a risk factor for PLEs, over and above genetic vulnerability for schizophrenia. This research does not support the notion that the cannabis-psychosis link is limited to individuals who are genetically predisposed to psychosis and suggests a need for research focusing on cannabis-related processes in psychosis that cannot be explained by genetic vulnerability.

In vitro modeling of the neurobiological effects of glucocorticoids: A review

Hypothalamic-pituitary adrenal (HPA)axis dysregulation has long been implicated in stress-related disorders such as major depression and post-traumatic stress disorder. Glucocorticoids (GCs) are released from the adrenal glands as a result of HPA-axis activation. The release of GCs is implicated with several neurobiological changes that are associated with negative consequences of chronic stress and the onset and course of psychiatric disorders. Investigating the underlying neurobiological effects of GCs may help to better understand the pathophysiology of stress-related psychiatric disorders. GCs impact a plethora of neuronal processes at the genetic, epigenetic, cellular, and molecular levels. Given the scarcity and difficulty in accessing human brain samples, 2D and 3D in vitro neuronal cultures are becoming increasingly useful in studying GC effects. In this review, we provide an overview of in vitro studies investigating the effects of GCs on key neuronal processes such as proliferation and survival of progenitor cells, neurogenesis, synaptic plasticity, neuronal activity, inflammation, genetic vulnerability, and epigenetic alterations. Finally, we discuss the challenges in the field and offer suggestions for improving the use of in vitro models to investigate GC effects.

Childhood trauma is associated with reduced frontal gray matter volume: a large transdiagnostic structural MRI study

BACKGROUND

Childhood trauma increases risk for psychopathology and cognitive impairment. Prior research mainly focused on the hippocampus and amygdala in single diagnostic categories. However, other brain regions may be impacted by trauma as well, and effects may be independent of diagnosis. This cross-sectional study investigated cortical and subcortical gray matter volume in relation to childhood trauma severity.

METHODS

We included 554 participants: 250 bipolar-I patients, 84 schizophrenia-spectrum patients and 220 healthy individuals without a psychiatric history. Participants filled in the Childhood Trauma Questionnaire. Anatomical T1 MRI scans were acquired at 3T, regional brain morphology was assessed using Freesurfer.

RESULTS

In the total sample, trauma-related gray matter reductions were found in the frontal lobe (β = -0.049, p = 0.008; q = 0.048), this effect was driven by the right medial orbitofrontal, paracentral, superior frontal regions and the left precentral region. No trauma-related volume reductions were observed in any other (sub)cortical lobes nor the hippocampus or amygdala, trauma-by-group (i.e. both patient groups and healthy subjects) interaction effects were absent. A categorical approach confirmed a pattern of more pronounced frontal gray matter reductions in individuals reporting multiple forms of trauma and across quartiles of cumulative trauma scores. Similar dose-response patterns were revealed within the bipolar and healthy subgroups, but did not reach significance in schizophrenia-spectrum patients.

CONCLUSIONS

Findings show that childhood trauma is linked to frontal gray matter reductions, independent of psychiatric morbidity. Our results indicate that childhood trauma importantly contributes to the neurobiological changes commonly observed across psychiatric disorders. Frontal volume alterations may underpin affective and cognitive disturbances observed in trauma-exposed individuals.

Rare copy number variation in posttraumatic stress disorder

Posttraumatic stress disorder (PTSD) is a heritable (h2 = 24-71%) psychiatric illness. Copy number variation (CNV) is a form of rare genetic variation that has been implicated in the etiology of psychiatric disorders, but no large-scale investigation of CNV in PTSD has been performed. We present an association study of CNV burden and PTSD symptoms in a sample of 114,383 participants (13,036 cases and 101,347 controls) of European ancestry. CNVs were called using two calling algorithms and intersected to a consensus set. Quality control was performed to remove strong outlier samples. CNVs were examined for association with PTSD within each cohort using linear or logistic regression analysis adjusted for population structure and CNV quality metrics, then inverse variance weighted meta-analyzed across cohorts. We examined the genome-wide total span of CNVs, enrichment of CNVs within specified gene-sets, and CNVs overlapping individual genes and implicated neurodevelopmental regions. The total distance covered by deletions crossing over known neurodevelopmental CNV regions was significant (beta = 0.029, SE = 0.005, P = 6.3 × 10-8). The genome-wide neurodevelopmental CNV burden identified explains 0.034% of the variation in PTSD symptoms. The 15q11.2 BP1-BP2 microdeletion region was significantly associated with PTSD (beta = 0.0206, SE = 0.0056, P = 0.0002). No individual significant genes interrupted by CNV were identified. 22 gene pathways related to the function of the nervous system and brain were significant in pathway analysis (FDR q < 0.05), but these associations were not significant once NDD regions were removed. A larger sample size, better detection methods, and annotated resources of CNV are needed to explore this relationship further.

A computational solution for bolstering reliability of epigenetic clocks: Implications for clinical trials and longitudinal tracking

Epigenetic clocks are widely used aging biomarkers calculated from DNA methylation data, but this data can be surprisingly unreliable. Here we show technical noise produces deviations up to 9 years between replicates for six prominent epigenetic clocks, limiting their utility. We present a computational solution to bolster reliability, calculating principal components from CpG-level data as input for biological age prediction. Our retrained principal-component versions of six clocks show agreement between most replicates within 1.5 years, improved detection of clock associations and intervention effects, and reliable longitudinal trajectories in vivo and in vitro. This method entails only one additional step compared to traditional clocks, requires no replicates or prior knowledge of CpG reliabilities for training, and can be applied to any existing or future epigenetic biomarker. The high reliability of principal component-based clocks is critical for applications to personalized medicine, longitudinal tracking, in vitro studies, and clinical trials of aging interventions.

Clinical profiles of subsequent stages in bipolar disorder: Results from the Dutch Bipolar Cohort

INTRODUCTION

The manifestation of bipolar disorder (BD) is hypothesized to be determined by clinical characteristics such as familial loading, childhood abuse, age at onset, illness duration, comorbid psychiatric disorders, addiction, treatment resistance, and premorbid cognitive functioning. Which of these are associated with a more severe course and worse outcome is currently unknown. Our objective is to find a combination of clinical characteristics associated with advancement to subsequent stages in two clinical staging models for BD.

METHODS

Using cross-sectional data from the Dutch Bipolar Cohort, staging was applied to determine the progression of bipolar-I-disorder (BD-I; N = 1396). Model A is primarily defined by recurrence of mood episodes, ranging from prodromal to chronicity. Model B is defined by level of inter-episodic functioning, ranging from prodromal to inability to function autonomously. For both models, ordinal logistic regression was conducted to test which clinical characteristics are associated with subsequent stages.

RESULTS

For model A, familial loading, childhood abuse, earlier onset, longer illness duration, psychiatric comorbidity, and treatment resistance were all predictors for a higher stage in contrast to addiction and cognitive functioning. For model B, childhood abuse, psychiatric comorbidity, cognitive functioning, and treatment resistance were predictors for a more severe stage, whereas age at onset, illness duration, and addiction were not.

DISCUSSION/CONCLUSIONS

Differences in clinical characteristics across stages support the construct validity of both staging models. Characteristics associated with a higher stage largely overlapped across both models. This study is a first step toward determining different clinical profiles, with a corresponding course and outcome.

Childhood maltreatment mediates the effect of the genetic background on psychosis risk in young adults

Childhood maltreatment (CM) and genetic vulnerability are both risk factors for psychosis, but the relations between them are not fully understood. Guided by the recent identification of genetic risk to CM, this study investigates the hypothesis that genetic risk to schizophrenia also increases the risk of CM and thus impacts psychosis risk. The relationship between schizophrenia polygenetic risk, CM, and psychotic-like experiences (PLE) was investigated in participants from the Utrecht Cannabis Cohort (N = 1262) and replicated in the independent IMAGEN cohort (N = 1740). Schizophrenia polygenic risk score (SZ-PRS) were calculated from the most recent GWAS. The relationship between CM, PRS, and PLE was first investigated using multivariate linear regression. Next, mediation of CM in the pathway linking SZ-PRS and PLE was examined by structural equation modeling, while adjusting for a set of potential mediators including cannabis use, smoking, and neuroticism. In agreement with previous studies, PLE were strongly associated with SZ-PRS (B = 0.190, p = 0.009) and CM (B = 0.575, p < 0.001). Novel was that CM was also significantly associated with SZ-PRS (B = 0.171, p = 0.001), and substantially mediated the effects of SZ-PRS on PLE (proportion mediated = 29.9%, p = 0.001). In the replication cohort, the analyses yielded similar results, confirming equally strong mediation by CM (proportion mediated = 34.7%, p = 0.009). Our results suggest that CM acts as a mediator in the causal pathway linking SZ-PRS and psychosis risk. These findings open new perspectives on the relations between genetic and environmental risks and warrant further studies into potential interventions to reduce psychosis risk in vulnerable people.

Characterization of HIV-1 Infection in Microglia-Containing Human Cerebral Organoids

The achievement of an HIV cure is dependent on the eradication or permanent silencing of HIV-latent viral reservoirs, including the understudied central nervous system (CNS) reservoir. This requires a deep understanding of the molecular mechanisms of HIV’s entry into the CNS, latency establishment, persistence, and reversal. Therefore, representative CNS culture models that reflect the intercellular dynamics and pathophysiology of the human brain are urgently needed in order to study the CNS viral reservoir and HIV-induced neuropathogenesis. In this study, we characterized a human cerebral organoid model in which microglia grow intrinsically as a CNS culture model to study HIV infection in the CNS. We demonstrated that both cerebral organoids and isolated organoid-derived microglia (oMG), infected with replication-competent HIVbal reporter viruses, support productive HIV infection via the CCR5 co-receptor. Productive HIV infection was only observed in microglial cells. Fluorescence analysis revealed microglia as the only HIV target cell. Susceptibility to HIV infection was dependent on the co-expression of microglia-specific markers and the CD4 and CCR5 HIV receptors. Altogether, this model will be a valuable tool within the HIV research community to study HIV–CNS interactions, the underlying mechanisms of HIV-associated neurological disorders (HAND), and the efficacy of new therapeutic and curative strategies on the CNS viral reservoir.

Biomarkers in PTSD-susceptible and resistant veterans with war experience of more than ten years ago: FOCUS ON cortisol, thyroid hormones, testosterone and GABA

In the present study we measured the concentrations of cortisol, thyroid hormones, testosterone, and GABA (gamma aminobutyric acid) in am blood plasma samples of combatants with an at least 10 year history of military psychological trauma (N = 74) divided in groups that either suffer from post-traumatic stress disorder (PTSD) (N = 37) or are resistant (N = 37) as well as in a control group without traumatic experience in the anamnesis, (N = 34). PTSD symptoms were assessed using the Clinician-Administered PTSD Scale (CAPS). The results show that the am blood cortisol levels of individuals that were exposed to war zone experiences irrespective susceptibility for or resistance to PTSD were significantly higher than the values observed in the controls. Testosterone levels in PTSD patients differed neither from that observed in PTSD resistant nor control groups. In the resistant group testosterone levels were however significantly higher than in controls. The level of all thyroid hormones did not differ between the study groups. GABA level was significantly lower in the PTSD group compared with healthy controls. In the resistant group blood GABA levels were not significantly different from either PTSD patients or controls. In conclusion, the current data show that cortisol and to some extent testosterone may serve as biomarker of war zone stress per se, even if trauma was experienced at least ten years before, rather than of only PTSD or resistance to PTSD. GABA, in contrast, can be considered a potential marker of the protracted nature of PTSD.

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.

Enhancing Discovery of Genetic Variants for Posttraumatic Stress Disorder Through Integration of Quantitative Phenotypes and Trauma Exposure Information

BACKGROUND

Posttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).

METHODS

A GWAS on PTSD symptoms was performed in 51 cohorts followed by a fixed-effects meta-analysis (N = 182,199 European ancestry participants). A GWAS of LTE burden was performed in the UK Biobank cohort (N = 132,988). Genetic correlations were evaluated with linkage disequilibrium score regression. Multivariate analysis was performed using Multi-Trait Analysis of GWAS. Functional mapping and annotation of leading loci was performed with FUMA. Replication was evaluated using the Million Veteran Program GWAS of PTSD total symptoms.

RESULTS

GWASs of PTSD symptoms and LTE burden identified 5 and 6 independent genome-wide significant loci, respectively. There was a 72% genetic correlation between PTSD and LTE. PTSD and LTE showed largely similar patterns of genetic correlation with other traits, albeit with some distinctions. Adjusting PTSD for LTE reduced PTSD heritability by 31%. Multivariate analysis of PTSD and LTE increased the effective sample size of the PTSD GWAS by 20% and identified 4 additional loci. Four of these 9 PTSD loci were independently replicated in the Million Veteran Program.

CONCLUSIONS

Through using a quantitative trait measure of PTSD, we identified novel risk loci not previously identified using prior case-control analyses. PTSD and LTE have a high genetic overlap that can be leveraged to increase discovery power through multivariate methods.

Interrogating Associations Between Polygenic Liabilities and Electroconvulsive Therapy Effectiveness

BACKGROUND

Electroconvulsive therapy (ECT) is the most effective treatment for severe major depressive episodes (MDEs). Nonetheless, firmly established associations between ECT outcomes and biological variables are currently lacking. Polygenic risk scores (PRSs) carry clinical potential, but associations with treatment response in psychiatry are seldom reported. Here, we examined whether PRSs for major depressive disorder, schizophrenia (SCZ), cross-disorder, and pharmacological antidepressant response are associated with ECT effectiveness.

METHODS

A total of 288 patients with MDE from 3 countries were included. The main outcome was a change in the 17-item Hamilton Depression Rating Scale scores from before to after ECT treatment. Secondary outcomes were response and remission. Regression analyses with PRSs as independent variables and several covariates were performed. Explained variance (R²) at the optimal p-value threshold is reported.

RESULTS

In the 266 subjects passing quality control, the PRS-SCZ was positively associated with a larger Hamilton Depression Rating Scale decrease in linear regression (optimal p-value threshold = .05, R² = 6.94%, p < .0001), which was consistent across countries: Ireland (R² = 8.18%, p = .0013), Belgium (R² = 6.83%, p = .016), and the Netherlands (R² = 7.92%, p = .0077). The PRS-SCZ was also positively associated with remission (R² = 4.63%, p = .0018). Sensitivity and subgroup analyses, including in MDE without psychotic features (R² = 4.42%, p = .0024) and unipolar MDE only (R² = 9.08%, p < .0001), confirmed the results. The other PRSs were not associated with a change in the Hamilton Depression Rating Scale score at the predefined Bonferroni-corrected significance threshold.

CONCLUSIONS

A linear association between PRS-SCZ and ECT outcome was uncovered. Although it is too early to adopt PRSs in ECT clinical decision making, these findings strengthen the positioning of PRS-SCZ as relevant to treatment response in psychiatry.

Untargeted metabolic analysis in dried blood spots reveals metabolic signature in 22q11.2 deletion syndrome

The 22q11.2 deletion syndrome (22q11.2DS) is characterized by a well-defined microdeletion and is associated with increased risk of neurodevelopmental phenotypes including autism spectrum disorders (ASD) and intellectual impairment. The typically deleted region in 22q11.2DS contains multiple genes with the potential of altering metabolism. Deficits in metabolic processes during early brain development may help explain the increased prevalence of neurodevelopmental phenotypes seen in 22q11.2DS. However, relatively little is known about the metabolic impact of the 22q11.2 deletion, while such insight may lead to increased understanding of the etiology. We performed untargeted metabolic analysis in a large sample of dried blood spots derived from 49 22q11.2DS patients and 87 controls, to identify a metabolic signature for 22q11.2DS. We also examined trait-specific metabolomic patterns within 22q11.2DS patients, focusing on intelligence (intelligence quotient, IQ) and ASD. We used the Boruta algorithm to select metabolites distinguishing patients from controls, patients with ASD from patients without, and patients with an IQ score in the lowest range from patients with an IQ score in the highest range. The relevance of the selected metabolites was visualized with principal component score plots, after which random forest analysis and logistic regression were used to measure predictive performance of the selected metabolites. Analysis yielded a distinct metabolic signature for 22q11.2DS as compared to controls, and trait-specific (IQ and ASD) metabolomic patterns within 22q11.2DS patients. The metabolic characteristics of 22q11.2DS provide insights in biological mechanisms underlying the neurodevelopmental phenotype and may ultimately aid in identifying novel therapeutic targets for patients with developmental disorders.

Epigenome-wide meta-analysis of PTSD symptom severity in three military cohorts implicates DNA methylation changes in genes involved in immune system and oxidative stress

Epigenetic factors modify the effects of environmental factors on biological outcomes. Identification of epigenetic changes that associate with PTSD is therefore a crucial step in deciphering mechanisms of risk and resilience. In this study, our goal is to identify epigenetic signatures associated with PTSD symptom severity (PTSS) and changes in PTSS over time, using whole blood DNA methylation (DNAm) data (MethylationEPIC BeadChip) of military personnel prior to and following combat deployment. A total of 429 subjects (858 samples across 2 time points) from three male military cohorts were included in the analyses. We conducted two different meta-analyses to answer two different scientific questions: one to identify a DNAm profile of PTSS using a random effects model including both time points for each subject, and the other to identify a DNAm profile of change in PTSS conditioned on pre-deployment DNAm. Four CpGs near four genes (F2R, CNPY2, BAIAP2L1, and TBXAS1) and 88 differentially methylated regions (DMRs) were associated with PTSS. Change in PTSS after deployment was associated with 15 DMRs, of those 2 DMRs near OTUD5 and ELF4 were also associated with PTSS. Notably, three PTSS-associated CpGs near F2R, BAIAP2L1 and TBXAS1 also showed nominal evidence of association with change in PTSS. This study, which identifies PTSD-associated changes in genes involved in oxidative stress and immune system, provides novel evidence that epigenetic differences are associated with PTSS.

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.

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.

Shape and volume changes of the superior lateral ventricle after electroconvulsive therapy measured with ultra-high field MRI

The subventricular zone (SVZ) of the lateral ventricles harbors neuronal stem cells in adult mammals. Rodent studies report neurogenic effects in the SVZ of electroconvulsive stimulation. We hypothesize that if this finding translates to depressed patients undergoing electroconvulsive therapy (ECT), this would be reflected in shape changes at the SVZ. Using T1-weighted MR images acquired at ultra-high field strength (7T), the shape and volume of the ventricles were compared from pre to post ECT after 10 ECT sessions (in patients twice weekly) or 5 weeks apart (controls) using linear mixed models with age and gender as covariates. Ventricle shape significantly changed and volume significantly decreased over time in patients for the left ventricle, but not in controls. The decrease in volume of the ventricles was associated to a decrease in depression scores, and an increase in the left dentate gyrus, However, the shape changes of the ventricles were not restricted to the neurogenic niche in the lateral walls of the ventricles, providing no clear evidence for neurogenesis as sole explanation of volume changes in the ventricles after ECT.

A systematic review of biological, social and environmental factors associated with epigenetic clock acceleration

Aging involves a diverse set of biological changes accumulating over time that leads to increased risk of morbidity and mortality. Epigenetic clocks are now widely used to quantify biological aging, in order to investigate determinants that modify the rate of aging and to predict age-related outcomes. Numerous biological, social and environmental factors have been investigated for their relationship to epigenetic clock acceleration and deceleration. The aim of this review was to synthesize general trends concerning the associations between human epigenetic clocks and these investigated factors. We conducted a systematic review of all available literature and included 156 publications across 4 resource databases. We compiled a list of all presently existing blood-based epigenetic clocks. Subsequently, we created an extensive dataset of over 1300 study findings in which epigenetic clocks were utilized in blood tissue of human subjects to assess the relationship between these clocks and numeral environmental exposures and human traits. Statistical analysis was possible on 57 such relationships, measured across 4 different epigenetic clocks (Hannum, Horvath, Levine and GrimAge). We found that the Horvath, Hannum, Levine and GrimAge epigenetic clocks tend to agree in direction of effects, but vary in size. Body mass index, HIV infection, and male sex were significantly associated with acceleration of one or more epigenetic clocks. Acceleration of epigenetic clocks was also significantly related to mortality, cardiovascular disease, cancer and diabetes. Our findings provide a graphical and numerical synopsis of the past decade of epigenetic age estimation research and indicate areas where further attention could be focused in the coming years.

DNA methylation differences in cortical grey and white matter in schizophrenia

Aim: Identify grey- and white-matter-specific DNA-methylation differences between schizophrenia (SCZ) patients and controls in postmortem brain cortical tissue. Materials & methods: Grey and white matter were separated from postmortem brain tissue of the superior temporal and medial frontal gyrus from SCZ (n = 10) and control (n = 11) cases. Genome-wide DNA-methylation analysis was performed using the Infinium EPIC Methylation Array (Illumina, CA, USA). Results: Four differentially methylated regions associated with SCZ status and tissue type (grey vs white matter) were identified within or near KLF9, SFXN1, SPRED2 and ALS2CL genes. Gene-expression analysis showed differential expression of KLF9 and SFXN1 in SCZ. Conclusion: Our data show distinct differences in DNA methylation between grey and white matter that are unique to SCZ, providing new leads to unravel the pathogenesis of SCZ.

Distinct non-inflammatory signature of microglia in post-mortem brain tissue of patients with major depressive disorder

Findings from epidemiological studies, biomarker measurements and animal experiments suggest a role for aberrant immune processes in the pathogenesis of major depressive disorder (MDD). Microglia, the resident immune cells of the brain, are likely to play a key role in these processes. Previous post-mortem studies reported conflicting findings regarding microglial activation and an in-depth profiling of those cells in MDD is lacking. The aim of this study was therefore to characterize the phenotype and function of microglia in MDD. We isolated microglia from post-mortem brain tissue of patients with MDD (n = 13-19) and control donors (n = 12-25). Using flow cytometry and quantitative Polymerase Chain Reaction (qPCR), we measured protein and mRNA levels of a panel of microglial markers across four different brain regions (medial frontal gyrus, superior temporal gyrus, thalamus, and subventricular zone). In MDD cases, we found a significant upregulation of CX3CR1 and TMEM119 mRNA expression and a downregulation of CD163 mRNA expression and CD14 protein expression across the four brain regions. Expression levels of microglial activation markers, such as HLA-DRA, IL6, and IL1β, as well as the inflammatory responses to lipopolysaccharide and dexamethasone were unchanged. Our findings suggest that microglia enhance homeostatic functions in MDD but are not immune activated.

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.

Successful treatment of post-traumatic stress disorder reverses DNA methylation marks

Epigenetic mechanisms play a role in the detrimental effects of traumatic stress and the development of post-traumatic stress disorder (PTSD). However, it is unknown whether successful treatment of PTSD restores these epigenetic marks. This study investigated longitudinal changes of blood-based genome-wide DNA methylation levels in relation to trauma-focused psychotherapy for PTSD in soldiers that obtained remission (N = 21), non-remitted PTSD patients (N = 23), and trauma-exposed military controls (N = 23). In an independent prospective cohort, we then examined whether these DMRs were also relevant for the development of deployment-related PTSD (N = 85). Successful treatment of PTSD was accompanied by significant changes in DNA methylation at 12 differentially methylated regions (DMRs) in the genes: APOB, MUC4, EDN2, ZFP57, GPX6, CFAP45, AFF3, TP73, UBCLP1, RPL13P, and two intergenic regions (p values < 0.0001 were confirmed using permutation and sensitivity analyses). Of the 12 DMRs related to PTSD symptom reduction, consistent prospective evidence was found for ZFP57 methylation changes related to changing PTSD symptoms (B = -0.84, t = -2.49, p = 0.014). Increasing ZFP57 methylation related to PTSD symptom reduction was present over and above the relation with symptoms, suggesting that psychological treatments exert biological effects independent of symptom reduction. Together, these data provide longitudinal evidence that ZFP57 methylation is involved in both the development and successful treatment of deployment-related PTSD. This study is a first step to disentangle the interaction between psychological and biological systems to identify genomic regions relevant for the etiology and treatment of stress-related disorders such as PTSD.

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.

Epigenome-wide meta-analysis of PTSD across 10 military and civilian cohorts identifies methylation changes in AHRR

Epigenetic differences may help to distinguish between PTSD cases and trauma-exposed controls. Here, we describe the results of the largest DNA methylation meta-analysis of PTSD to date. Ten cohorts, military and civilian, contribute blood-derived DNA methylation data from 1,896 PTSD cases and trauma-exposed controls. Four CpG sites within the aryl-hydrocarbon receptor repressor (AHRR) associate with PTSD after adjustment for multiple comparisons, with lower DNA methylation in PTSD cases relative to controls. Although AHRR methylation is known to associate with smoking, the AHRR association with PTSD is most pronounced in non-smokers, suggesting the result was independent of smoking status. Evaluation of metabolomics data reveals that AHRR methylation associated with kynurenine levels, which are lower among subjects with PTSD. This study supports epigenetic differences in those with PTSD and suggests a role for decreased kynurenine as a contributor to immune dysregulation in PTSD.

Associations between the development of PTSD symptoms and longitudinal changes in the DNA methylome of deployed military servicemen: A comparison with polygenic risk scores

Objective

Method

Results

Conclusion

•

First study of methylation changes over a deployment period using three time-points.

•

Post-deployment methylation decreases in EP300 were associated with delayed onset PTSD.

•

Genetic risk was exclusively related to a non-delayed PTSD onset.

Methylation of oxytocin related genes and early life trauma together shape the N170 response to human faces

Childhood trauma fundamentally shapes social cognition and basic processing of social cues, which frequently cascade into adverse behavioral outcomes. Recent studies indicate that epigenetic changes in oxytocin functioning might contribute to these long-term effects, although a deeper understanding of the underlying mechanisms is still lacking. The electroencephalographic N170 response to faces might capture a neural response at the core of these interactive effects of oxytocin gene methylation and childhood adversity, given that this response is considered to reflect fundamental face processing, to be susceptible to oxytocin administration and also to be a biomarker of various psychiatric disorders. We assessed the N170 response to neutral faces in relation to participant's (81, women) recalled childhood trauma, methylation of their oxytocin structural (OXTg) and oxytocin receptor (OXTRg) genes, and endogenous levels of cortisol and testosterone. Additionally, we investigated the interactive effect of OXTg methylation and CTQ across three face sets of varying maturity. Methylation of OXTg relates to a weakened N170 response towards adults, children and infants. Moreover, methylation of both OXTRg and OXTg shaped the directionality of adversity effects, predicting a weakened N170 response in those with high methylation and hyper-vigilance with participants with low methylation. Our results are the first to relate OXT(R)g methylation to the N170 response. They shed light on biological processes linking childhood adversity and epigenetic marks to altered behavior and potentially psychopathologies.

Functional brain networks in the schizophrenia spectrum and bipolar disorder with psychosis

Psychotic experiences have been proposed to lie on a spectrum, ranging from subclinical experiences to treatment-resistant schizophrenia. We aimed to characterize functional connectivity and brain network characteristics in relation to the schizophrenia spectrum and bipolar disorder with psychosis to disentangle neural correlates to psychosis. Additionally, we studied antipsychotic medication and lithium effects on network characteristics. We analyzed functional connectivity strength and network topology in 487 resting-state functional MRI scans of individuals with schizophrenia spectrum disorder (SCZ), bipolar disorder with a history of psychotic experiences (BD), treatment-naïve subclinical psychosis (SCP), and healthy controls (HC). Since differences in connectivity strength may confound group comparisons of brain network topology, we analyzed characteristics of the minimum spanning tree (MST), a relatively unbiased backbone of the network. SCZ and SCP subjects had a lower connectivity strength than BD and HC individuals but showed no differences in network topology. In contrast, BD patients showed a less integrated network topology but no disturbances in connectivity strength. No differences in outcome measures were found between SCP and SCZ, or between BD patients that used antipsychotic medication or lithium and those that did not. We conclude that functional networks in patients prone to psychosis have different signatures for chronic SCZ patients and SCP compared to euthymic BD patients, with a limited role for medication. Connectivity strength effects may have confounded previous studies, as no functional network alterations were found in SCZ after strict correction for connectivity strength.

Higher cortisol levels may proceed a manic episode and are related to disease severity in patients with bipolar disorder

BACKGROUND

Regulation of the hypothalamic-pituitary-adrenal (HPA) axis is implicated in the pathogenesis of bipolar disorder (BD). However, the relationship between HPA-activity and disease severity is not fully elucidated. In this pilot study we aimed to explore the temporal relationship between HPA-activity and the risk of a manic episode in BD patients type I, by assessing long-term hair cortisol concentrations (HCC). Second, we explored the relation between HCC and the number of previous episodes.

METHODS

Hair samples were collected from 45 BD I patients in euthymic or manic state and compared to 17 controls. From each participant, two hair samples of 3 cm length were used to measure long-term cortisol, reflecting retrospect time frames of 1-3 months and 4-6 months respectively prior to sampling.

RESULTS

HCC in the BD group was slightly higher than in the control group in both hair segments (p = 0.049 and 0.03; after adjustment for age, sex, BMI and hair washing frequency p = 0.222 and 0.139). A significant peak in hair cortisol was observed prior to a manic episode (p = 0.036). Furthermore, we found a positive correlation between the number of mood episodes HCC (p = 0.03).

CONCLUSIONS

Our results indicate that long-term cortisol levels are slightly higher in BD, and in particular elevated in the months prior to a manic relapse. In addition HCC are positively associated with the number of previous mood episodes in the course of BD type I.

Cannabinoids and psychotic symptoms: A potential role for a genetic variant in the P2X purinoceptor 7 (P2RX7) gene

To investigate the biological mechanisms underlying the higher risk for psychosis in those that use cannabis, we conducted a genome-wide environment-interaction study (GWEIS). In a sample of individuals without a psychiatric disorder (N = 1262), we analyzed the interactions between regular cannabis use and genotype with psychotic-like experiences (PLE) as outcome. PLE were measured using the Community Assessment of Psychic Experiences (CAPE). The sample was enriched for those at the extremes of both cannabis use and PLE to increase power. A single nucleotide polymorphism in the P2RX7 gene (rs7958311) was associated with risk for a high level of psychotic experiences in regular cannabis users (p = 1.10 x10^-7) and in those with high levels of lifetime cannabis use (p = 4.5 × 10^-6). This interaction was replicated in individuals with high levels of lifetime cannabis use in the IMAGEN cohort (N = 1217, p = 0.020). Functional relevance of P2RX7 in cannabis users was suggested by in vitro experiments on activated monocytes. Exposure of these cells to tetrahydrocannabinol (THC) or cannabidiol (CBD) reduced the immunological response of the P2X7 receptor, which was dependent on the identified genetic variant. P2RX7 variants have been implicated in psychiatric disorders before and the P2X7 receptor is involved in pathways relevant to psychosis, such as neurotransmission, synaptic plasticity and immune regulation. We conclude that P2RX7 plays a role in vulnerability to develop psychotic symptoms when using cannabis and point to a new pathway that can potentially be targeted by newly developed P2X7 antagonists.

Schizophrenia and Epigenetic Aging Biomarkers: Increased Mortality, Reduced Cancer Risk, and Unique Clozapine Effects

BACKGROUND

Schizophrenia (SZ) is associated with increased all-cause mortality, smoking, and age-associated proteins, yet multiple previous studies found no association between SZ and biological age using Horvath's epigenetic clock, a well-established aging biomarker based on DNA methylation. However, numerous epigenetic clocks that may capture distinct aspects of aging have been developed. This study tested the hypothesis that altered aging in SZ manifests in these other clocks.

METHODS

We performed a comprehensive analysis of 14 epigenetic clocks categorized according to what they were trained to predict: chronological age, mortality, mitotic divisions, or telomere length. To understand the etiology of biological age differences, we also examined DNA methylation predictors of smoking, alcohol, body mass index, serum proteins, and cell proportions. We independently analyzed 3 publicly available multiethnic DNA methylation data sets from whole blood, a total of 567 SZ cases and 594 nonpsychiatric controls.

RESULTS

All data sets showed accelerations in SZ for the 3 mortality clocks up to 5 years, driven by smoking and elevated levels of 6 age-associated proteins. The 2 mitotic clocks were decelerated in SZ related to antitumor natural killer and CD8T cells, which may help explain conflicting reports about low cancer rates in epidemiological studies of SZ. One cohort with available medication data showed that clozapine is associated with male-specific decelerations up to 7 years in multiple chronological age clocks.

CONCLUSIONS

Our study demonstrates the utility of studying the various epigenetic clocks in tandem and highlights potential mechanisms by which mental illness influences long-term outcomes, including cancer and early mortality.

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.

Evolutionary modifications in human brain connectivity associated with schizophrenia

The genetic basis and human-specific character of schizophrenia has led to the hypothesis that human brain evolution may have played a role in the development of the disorder. We examined schizophrenia-related changes in brain connectivity in the context of evolutionary changes in human brain wiring by comparing in vivo neuroimaging data from humans and chimpanzees, one of our closest living evolutionary relatives and a species with which we share a very recent common ancestor. We contrasted the connectome layout between the chimpanzee and human brain and compared differences with the pattern of schizophrenia-related changes in brain connectivity as observed in patients. We show evidence of evolutionary modifications of human brain connectivity to significantly overlap with the cortical pattern of schizophrenia-related dysconnectivity (P < 0.001, permutation testing). We validated these effects in three additional, independent schizophrenia datasets. We further assessed the specificity of effects by examining brain dysconnectivity patterns in seven other psychiatric and neurological brain disorders (including, among others, major depressive disorder and obsessive-compulsive disorder, arguably characterized by behavioural symptoms that are less specific to humans), which showed no such associations with modifications of human brain connectivity. Comparisons of brain connectivity across humans, chimpanzee and macaques further suggest that features of connectivity that evolved in the human lineage showed the strongest association to the disorder, that is, brain circuits potentially related to human evolutionary specializations. Taken together, our findings suggest that human-specific features of connectome organization may be enriched for changes in brain connectivity related to schizophrenia. Modifications in human brain connectivity in service of higher order brain functions may have potentially also rendered the brain vulnerable to brain dysfunction.

An epigenome-wide association study of posttraumatic stress disorder in US veterans implicates several new DNA methylation loci

BACKGROUND

Previous studies using candidate gene and genome-wide approaches have identified epigenetic changes in DNA methylation (DNAm) associated with posttraumatic stress disorder (PTSD).

METHODS

In this study, we performed an EWAS of PTSD in a cohort of Veterans (n = 378 lifetime PTSD cases and 135 controls) from the Translational Research Center for TBI and Stress Disorders (TRACTS) cohort assessed using the Illumina EPIC Methylation BeadChip which assesses DNAm at more than 850,000 sites throughout the genome. Our model included covariates for ancestry, cell heterogeneity, sex, age, and a smoking score based on DNAm at 39 smoking-associated CpGs. We also examined in EPIC-based DNAm data generated from pre-frontal cortex (PFC) tissue from the National PTSD Brain Bank (n = 72).

RESULTS

The analysis of blood samples yielded one genome-wide significant association with PTSD at cg19534438 in the gene G0S2 (p = 1.19 × 10-7, padj = 0.048). This association was replicated in an independent PGC-PTSD-EWAS consortium meta-analysis of military cohorts (p = 0.0024). We also observed association with the smoking-related locus cg05575921 in AHRR despite inclusion of a methylation-based smoking score covariate (p = 9.16 × 10-6), which replicates a previously observed PGC-PTSD-EWAS association (Smith et al. 2019), and yields evidence consistent with a smoking-independent effect. The top 100 EWAS loci were then examined in the PFC data. One of the blood-based PTSD loci, cg04130728 in CHST11, which was in the top 10 loci in blood, but which was not genome-wide significant, was significantly associated with PTSD in brain tissue (in blood p = 1.19 × 10-5, padj = 0.60, in brain, p = 0.00032 with the same direction of effect). Gene set enrichment analysis of the top 500 EWAS loci yielded several significant overlapping GO terms involved in pathogen response, including "Response to lipopolysaccharide" (p = 6.97 × 10-6, padj = 0.042).

CONCLUSIONS

The cross replication observed in independent cohorts is evidence that DNA methylation in peripheral tissue can yield consistent and replicable PTSD associations, and our results also suggest that that some PTSD associations observed in peripheral tissue may mirror associations in the brain.

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.

Longitudinal epigenome-wide association studies of three male military cohorts reveal multiple CpG sites associated with post-traumatic stress disorder

BACKGROUND

Epigenetic mechanisms have been suggested to play a role in the development of post-traumatic stress disorder (PTSD). Here, blood-derived DNA methylation data (HumanMethylation450 BeadChip) collected prior to and following combat exposure in three cohorts of male military members were analyzed to assess whether DNA methylation profiles are associated with the development of PTSD. A total of 123 PTSD cases and 143 trauma-exposed controls were included in the analyses. The Psychiatric Genomics Consortium (PGC) PTSD EWAS QC pipeline was used on all cohorts, and results were combined using a sample size weighted meta-analysis in a two-stage design. In stage one, we jointly analyzed data of two new cohorts (N = 126 and 78) for gene discovery, and sought to replicate significant findings in a third, previously published cohort (N = 62) to assess the robustness of our results. In stage 2, we aimed at maximizing power for gene discovery by combining all three cohorts in a meta-analysis.

RESULTS

Stage 1 analyses identified four CpG sites in which, conditional on pre-deployment DNA methylation, post-deployment DNA methylation was significantly associated with PTSD status after epigenome-wide adjustment for multiple comparisons. The most significant (intergenic) CpG cg05656210 (p = 1.0 × 10-08) was located on 5q31 and significantly replicated in the third cohort. In addition, 19 differentially methylated regions (DMRs) were identified, but failed replication. Stage 2 analyses identified three epigenome-wide significant CpGs, the intergenic CpG cg05656210 and two additional CpGs located in MAD1L1 (cg12169700) and HEXDC (cg20756026). Interestingly, cg12169700 had an underlying single nucleotide polymorphism (SNP) which was located within the same LD block as a recently identified PTSD-associated SNP in MAD1L1. Stage 2 analyses further identified 12 significant differential methylated regions (DMRs), 1 of which was located in MAD1L1 and 4 were situated in the human leukocyte antigen (HLA) region.

CONCLUSIONS

This study suggests that the development of combat-related PTSD is associated with distinct methylation patterns in several genomic positions and regions. Our most prominent findings suggest the involvement of the immune system through the HLA region and HEXDC, and MAD1L1 which was previously associated with PTSD.

Multivariate genome-wide analysis of stress-related quantitative phenotypes

Exposure to traumatic stress increases the odds of developing a broad range of psychiatric conditions. Genetic studies targeting multiple stress-related quantitative phenotypes may shed light on mechanisms underlying vulnerability to psychopathology in the aftermath of stressful events. We applied a multivariate genome-wide association study (GWAS) to a unique military cohort (N = 583) in which we measured biochemical and behavioral phenotypes. The availability of pre- and post-deployment measurements allowed to capture changes in these phenotypes in response to stress. For genome-wide significant loci, we performed functional annotation, phenome-wide analysis and quasi-replication in PTSD case-control GWASs. We discovered one genetic variant reaching genome-wide significant association, surviving permutation and sensitivity analyses (rs10100651, p = 9.9 × 10^-9). Functional annotation prioritized the genes INTS8 and TP53INP1. A phenome-wide scan revealed a significant association of these same genes with sleeping problems, hypertension and subjective well-being. Finally, a targeted lookup revealed nominally significant association of rs10100651 in a PTSD case-control GWAS in the UK Biobank (p = 0.02). We provide comprehensive evidence from multiple resources hinting at a role of the highlighted genetic variant in the human stress response, marking the power of multivariate genome-wide analysis of quantitative measures in stress research. Future genetic and functional studies can target this locus to further assess its effects on stress mediation and its possible role in psychopathology or resilience.

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.

SCHIZOPHRENIA EPIGENETIC AGING PATTERNS REFLECT ALTERED MORTALITY AND CANCER RISKS

Schizophrenia (SZ) is associated with large increases in all-cause mortality, high smoking rates, and elevated levels of age-associated proteins—suggesting individuals with SZ may experience accelerated rates of biological aging. Yet surprisingly, multiple previous studies found no association between SZ and biological age using Horvath’s epigenetic clock, a well-recognized and validated biomarker of aging based on DNA methylation (DNAm) levels. However, numerous epigenetic clocks have been developed to date, many of which are better indicators of differential lifespan and healthspan than the original Horvath clock. Thus, we hypothesize that these epigenetic clocks may be better proxies for the presumed accelerated aging rate in SZ. Here we investigate 14 epigenetic clocks using three publicly available DNAm datasets from whole blood, comparing SZ to non-psychiatric controls (NPC). In all data sets, we find SZ age acceleration in three clocks previously shown to be most predictive of age-related morbidity and mortality risk. In contrast, two clocks developed to capture mitotic rate are decelerated in SZ, consistent with low cancer rates despite smoking observed in epidemiological studies of SZ. We use these clocks to investigate the determinants of altered aging in SZ, such as smoking, alcohol, BMI, age-associated proteins, blood cell composition, and psychotropic medications. Principal component analysis suggests mortality clock acceleration, mitotic clock deceleration, and medication effects are independent phenomena in SZ. Our study demonstrates the importance of studying the various epigenetic clocks in tandem and highlights their potential utility for understanding how mental illness influences long-term outcomes including cancer and early mortality.