Multi-omics analyses of cognitive traits and psychiatric disorders highlights brain-dependent mechanisms

Investigating causal relationships between gene expression and major depressive disorder via brain bulk-tissue and cell type-specific eQTL: A Mendelian randomization and Bayesian colocalization study

BACKGROUND

Major depressive disorder (MDD) is a highly prevalent psychiatric disorder with complex genetic underpinnings. While genome-wide association studies (GWAS) have identified multiple risk loci, pinpointing causal genes within the human brain remains challenging, particularly given the regulatory complexity across different cell types.

METHODS

We performed summary data-based MR (SMR) and Bayesian colocalization analyses by integrating bulk-tissue eQTL data from 888 individuals with single-cell eQTL datasets from 192 donors representing major brain cell types (excitatory and inhibitory neurons, astrocytes, microglia, oligodendrocytes, OPCs/COPs, endothelial cells, and pericytes). GWAS summary statistics for MDD (170,756 cases and 329,443 controls) were used to assess the causal impact of gene expression. Sensitivity analyses, including the heterogeneity in dependent instruments (HEIDI) test and Steiger filtering, ensured robust inference.

RESULTS

In bulk tissue analyses, five genes (BTN3A2, SLC12A5, AREL1, GMPPB, and ZNF660) emerged as having robust causal evidence for MDD, displaying consistent MR signals and strong colocalization. Cell type-specific analyses revealed additional candidate genes in excitatory neurons (FLOT1, AL450423.1), astrocytes (AL121821.1), and oligodendrocytes (YLPM1, COP1).

CONCLUSION

Our integrative approach reveals that causal gene expression profiles differ markedly between bulk-tissue and specific brain cell types, emphasizing cellular heterogeneity in MDD pathogenesis and informing precision therapeutic strategies. These findings underscore the necessity of considering cell type-specific gene regulation when developing therapeutic interventions for MDD.

Multiomics Analysis of the Molecular Response to Glucocorticoids: Insights Into Shared Genetic Risk From Psychiatric to Medical Disorders

BACKGROUND

Alterations in the effects of glucocorticoids have been implicated in mediating some of the negative health effects associated with chronic stress, including increased risk for psychiatric disorders and cardiovascular and metabolic diseases. In this study, we investigated how genetic variants influence gene expression and DNA methylation in response to glucocorticoid receptor (GR) activation and their association with disease risk.

METHODS

We measured DNA methylation (n = 199) and gene expression (n = 297) in peripheral blood before and after GR activation with dexamethasone, with matched genotype data available for all samples. A comprehensive molecular quantitative trait locus (QTL) analysis was conducted, mapping GR-response methylation (me)QTLs, GR-response expression (e)QTLs, and GR-response expression quantitative trait methylation (eQTMs). A multilevel network analysis was employed to map the complex relationships between the transcriptome, epigenome, and genetic variation.

RESULTS

We identified 3772 GR-response meCpGs corresponding to 104,828 local GR-response meQTLs that did not strongly overlap with baseline meQTLs. eQTM and eQTL analyses revealed distinct genetic influences on gene expression and DNA methylation. Multilevel network analysis uncovered GR-response network trio QTLs, characterized by SNP-CpG-transcript combinations where meQTLs act as both eQTLs and eQTMs. GR-response trio variants were enriched in a genome-wide association study for psychiatric, respiratory, autoimmune, and cardiovascular diseases and conferred a higher relative heritability per SNP than GR-response meQTL and baseline QTL SNPs.

CONCLUSIONS

Genetic variants modulating the molecular effects of glucocorticoids are associated with psychiatric as well as medical diseases and not uncovered in baseline QTL analyses.

MAAT: a new nonparametric Bayesian framework for incorporating multiple functional annotations in transcriptome-wide association studies

Transcriptome-wide association study (TWAS) has emerged as a powerful tool for translating the myriad variations identified by genome-wide association studies (GWAS) into regulated genes in the post-GWAS era. While integrating annotation information has been shown to enhance power, current annotation-assisted TWAS tools predominantly focus on epigenomic annotations. When including more annotations, the assumption of a positive correlation between annotation scores and SNPs' effect sizes, as adopted by current methods, often falls short. Here, we propose MAAT expanding the horizons of existing TWAS studies, generating a new model incorporating multiple annotations into TWAS and a new metric indicating the most important annotation.

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.

Brain transcriptome-wide association study implicates novel risk genes underlying schizophrenia risk

BACKGROUND

To identify risk genes whose expression are regulated by the reported risk variants and to explore the potential regulatory mechanism in schizophrenia (SCZ).

METHODS

We systematically integrated three independent brain expression quantitative traits (eQTLs) (CommonMind, GTEx, and BrainSeq Phase 2, a total of 1039 individuals) and GWAS data (56 418 cases and 78 818 controls), with the use of transcriptome-wide association study (TWAS). Diffusion magnetic resonance imaging was utilized to quantify the integrity of white matter bundles and determine whether polygenic risk of novel genes linked to brain structure was present in patients with first-episode antipsychotic SCZ.

RESULTS

TWAS showed that eight risk genes (CORO7, DDAH2, DDHD2, ELAC2, GLT8D1, PCDHA8, THOC7, and TYW5) reached transcriptome-wide significance (TWS) level. These findings were confirmed by an independent integrative approach (i.e. Sherlock). We further conducted conditional analyses and identified the potential risk genes that driven the TWAS association signal in each locus. Gene expression analysis showed that several TWS genes (including CORO7, DDAH2, DDHD2, ELAC2, GLT8D1, THOC7 and TYW5) were dysregulated in the dorsolateral prefrontal cortex of SCZ cases compared with controls. TWS genes were mainly expressed on the surface of glutamatergic neurons, GABAergic neurons, and microglia. Finally, SCZ cases had a substantially greater TWS genes-based polygenic risk (PRS) compared to controls, and we showed that fractional anisotropy of the cingulum-hippocampus mediates the influence of TWS genes PRS on SCZ.

CONCLUSIONS

Our findings identified novel SCZ risk genes and highlighted the importance of the TWS genes in frontal-limbic dysfunctions in SCZ, indicating possible therapeutic targets.

A single-nucleus transcriptome-wide association study implicates novel genes in depression pathogenesis

Background

Depression is a common psychiatric illness and global public health problem. However, our limited understanding of the biological basis of depression has hindered the development of novel treatments and interventions.

Methods

To identify new candidate genes for therapeutic development, we examined single-nucleus RNA sequencing (snucRNAseq) data from the dorsolateral prefrontal cortex (N=424) in relation to ante-mortem depressive symptoms. To complement these direct analyses, we also used genome-wide association study (GWAS) results for depression (N=500,199) along with genetic tools for inferring the expression of 22,159 genes in 7 cell types and 55 cell subtypes to perform transcriptome-wide association studies (TWAS) of depression followed by Mendelian randomization (MR).

Results

Our single-nucleus TWAS analysis identified 71 causal genes in depression that have a role in specific neocortical cell subtypes; 59 of 71 genes were novel compared to previous studies. Depression TWAS genes showed a cell type specific pattern, with the greatest enrichment being in both excitatory and inhibitory neurons as well as astrocytes. Gene expression in different neuron subtypes have different directions of effect on depression risk. Compared to lower genetically correlated traits (e.g. body mass index) with depression, higher correlated traits (e.g., neuroticism) have more common TWAS genes with depression. In parallel, we performed differential gene expression analysis in relation to depression in 55 cortical cell subtypes, and we found that genes such as ANKRD36, MADD, TAOK3, SCAI and CHUK are associated with depression in specific cell subtypes.

Conclusions

These two sets of analyses illustrate the utility of large snucRNAseq data to uncover both genes whose expression is altered in specific cell subtypes in the context of depression and to enhance the interpretation of well-powered GWAS so that we can prioritize specific susceptibility genes for further analysis and therapeutic development.

Gene co-expression architecture in peripheral blood in a cohort of remitted first-episode schizophrenia patients

A better understanding of schizophrenia subtypes is necessary to stratify the patients according to clinical attributes. To explore the genomic architecture of schizophrenia symptomatology, we analyzed blood co-expression modules and their association with clinical data from patients in remission after a first episode of schizophrenia. In total, 91 participants of the 2EPS project were included. Gene expression was assessed using the Clariom S Human Array. Weighted-gene co-expression network analysis (WGCNA) was applied to identify modules of co-expressed genes and to test its correlation with global functioning, clinical symptomatology, and premorbid adjustment. Among the 25 modules identified, six modules were significantly correlated with clinical data. These modules could be clustered in two groups according to their correlation with clinical data. Hub genes in each group showing overlap with risk genes for schizophrenia were enriched in biological processes related to metabolic processes, regulation of gene expression, cellular localization and protein transport, immune processes, and neurotrophin pathways. Our results indicate that modules with significant associations with clinical data showed overlap with gene sets previously identified in differential gene-expression analysis in brain, indicating that peripheral tissues could reveal pathogenic mechanisms. Hub genes involved in these modules revealed multiple signaling pathways previously related to schizophrenia, which may represent the complex interplay in the pathological mechanisms behind the disease. These genes could represent potential targets for the development of peripheral biomarkers underlying illness traits in clinical remission stages after a first episode of schizophrenia.

The emerging role of furin in neurodegenerative and neuropsychiatric diseases

Furin is an important mammalian proprotein convertase that catalyzes the proteolytic maturation of a variety of prohormones and proproteins in the secretory pathway. In the brain, the substrates of furin include the proproteins of growth factors, receptors and enzymes. Emerging evidence, such as reduced FURIN mRNA expression in the brains of Alzheimer's disease patients or schizophrenia patients, has implicated a crucial role of furin in the pathophysiology of neurodegenerative and neuropsychiatric diseases. Currently, compared to cancer and infectious diseases, the aberrant expression of furin and its pharmaceutical potentials in neurological diseases remain poorly understood. In this article, we provide an overview on the physiological roles of furin and its substrates in the brain, summarize the deregulation of furin expression and its effects in neurodegenerative and neuropsychiatric disorders, and discuss the implications and current approaches that target furin for therapeutic interventions. This review may expedite future studies to clarify the molecular mechanisms of furin deregulation and involvement in the pathogenesis of neurodegenerative and neuropsychiatric diseases, and to develop new diagnosis and treatment strategies for these diseases.

The human physiology of well-being: A systematic review on the association between neurotransmitters, hormones, inflammatory markers, the microbiome and well-being

To understand the pathways through which well-being contributes to health, we performed a systematic review according to the Preferred Reporting Items for Systematic Review and Meta-Analysis (PRISMA) guidelines on the association between well-being and physiological markers in four categories, neurotransmitters, hormones, inflammatory markers, and microbiome. We identified 91 studies. Neurotransmitter studies (k_{number of studies}=9) reported only a possible positive association between serotonin and well-being. For the hormone studies (k = 48), a lower momentary cortisol level was related to higher well-being (meta-analytic r = -0.06), and a steeper diurnal slope of cortisol levels. Inflammatory marker studies (k = 36) reported negative or non-significant relations with well-being, with meta-analytic estimates of respectively r = -0.07 and r = -0.05 for C-reactive protein and interleukin-6. Microbiome studies (k = 4) reported inconsistent associations between different bacteria abundance and well-being. The results indicate possible but small roles of serotonin, cortisol, and inflammatory markers in explaining differences in well-being. The inconsistent and limited results for other markers and microbiome require further research. Future directions for a complete picture of the physiological factors underlying well-being are proposed.

Genetics in psychiatry: Methods, clinical applications and future perspectives

Psychiatric disorders and related traits have a demonstrated genetic component, with heritability estimated by twin studies generally between 80% and 40%. Their pathogenesis is complex and multi-determined: environmental factors interact with a polygenic architecture, making difficult the development of models able to stratify patients or predict mental health outcomes. Despite this difficult challenge, relevant progress has been made in the field of psychiatric genetics in recent years. This review aims to present the main current methods in psychiatric genetics, their output, limitations, clinical applications, and possible future developments. Genome-wide association studies (GWASs) performed in increasingly large samples have led to the identification of replicated genetic loci associated with the risk of major psychiatric disorders, including schizophrenia and mood disorders. Statistical and biological approaches have been developed to improve our understanding of the etiopathogenetic mechanisms behind genome-wide significant associations, as well as for estimating the cumulative effect of risk variants at the individual level and the genetic overlap between different disorders, as pleiotropy is the rule rather than the exception. Clinical applications are available in the pharmacogenetics field. The main issues that remain to be addressed include improving ethnic diversity in genetic studies and the optimization of statistical power through methodological improvements, such as the definition of dimensional phenotypes with specific biological correlates and the integration of different types of omics data.

DNA methylation-based predictors of health: applications and statistical considerations

DNA methylation data have become a valuable source of information for biomarker development, because, unlike static genetic risk estimates, DNA methylation varies dynamically in relation to diverse exogenous and endogenous factors, including environmental risk factors and complex disease pathology. Reliable methods for genome-wide measurement at scale have led to the proliferation of epigenome-wide association studies and subsequently to the development of DNA methylation-based predictors across a wide range of health-related applications, from the identification of risk factors or exposures, such as age and smoking, to early detection of disease or progression in cancer, cardiovascular and neurological disease. This Review evaluates the progress of existing DNA methylation-based predictors, including the contribution of machine learning techniques, and assesses the uptake of key statistical best practices needed to ensure their reliable performance, such as data-driven feature selection, elimination of data leakage in performance estimates and use of generalizable, adequately powered training samples.

A comprehensive gene-centric pleiotropic association analysis for 14 psychiatric disorders with GWAS summary statistics

BACKGROUND

Recent genome-wide association studies (GWASs) have revealed the polygenic nature of psychiatric disorders and discovered a few of single-nucleotide polymorphisms (SNPs) associated with multiple psychiatric disorders. However, the extent and pattern of pleiotropy among distinct psychiatric disorders remain not completely clear.

METHODS

We analyzed 14 psychiatric disorders using summary statistics available from the largest GWASs by far. We first applied the cross-trait linkage disequilibrium score regression (LDSC) to estimate genetic correlation between disorders. Then, we performed a gene-based pleiotropy analysis by first aggregating a set of SNP-level associations into a single gene-level association signal using MAGMA. From a methodological perspective, we viewed the identification of pleiotropic associations across the entire genome as a high-dimensional problem of composite null hypothesis testing and utilized a novel method called PLACO for pleiotropy mapping. We ultimately implemented functional analysis for identified pleiotropic genes and used Mendelian randomization for detecting causal association between these disorders.

RESULTS

We confirmed extensive genetic correlation among psychiatric disorders, based on which these disorders can be grouped into three diverse categories. We detected a large number of pleiotropic genes including 5884 associations and 2424 unique genes and found that differentially expressed pleiotropic genes were significantly enriched in pancreas, liver, heart, and brain, and that the biological process of these genes was remarkably enriched in regulating neurodevelopment, neurogenesis, and neuron differentiation, offering substantial evidence supporting the validity of identified pleiotropic loci. We further demonstrated that among all the identified pleiotropic genes there were 342 unique ones linked with 6353 drugs with drug-gene interaction which can be classified into distinct types including inhibitor, agonist, blocker, antagonist, and modulator. We also revealed causal associations among psychiatric disorders, indicating that genetic overlap and causality commonly drove the observed co-existence of these disorders.

CONCLUSIONS

Our study is among the first large-scale effort to characterize gene-level pleiotropy among a greatly expanded set of psychiatric disorders and provides important insight into shared genetic etiology underlying these disorders. The findings would inform psychiatric nosology, identify potential neurobiological mechanisms predisposing to specific clinical presentations, and pave the way to effective drug targets for clinical treatment.