Transcriptome-wide association study identifies new susceptibility genes and pathways for depression

Sex-specific gene expression differences in the prefrontal cortex of major depressive disorder individuals

Major depressive disorder (MDD) is a leading global cause of disability, being more prevalent in females, possibly due to molecular and neuronal pathway differences between females and males. However, the connection between transcriptional changes and MDD remains unclear. We identified transcriptionally altered genes (TAGs) in MDD through gene and transcript expression analyses, focusing on sex-specific differences. Analyzing 263 brain samples from both sexes, we conducted differential gene expression, differential transcript expression, and differential transcript usage analyses, revealing 1169 unique TAGs, primarily in the prefrontal areas, with nearly half exhibiting transcript-level alterations. Females showed notable RNA splicing and export process disruptions in the orbitofrontal cortex, alongside altered DDX39B gene expression in five of the six brain regions in both sexes. Our findings suggest that disruptions in RNA processing pathways may play a vital role in MDD.

Differential interactions between gene expressions and stressors across the lifespan in major depressive disorder

BACKGROUND

Both genetic predispositions and exposures to stressors have collectively contributed to the development of major depressive disorder (MDD). To deep dive into their roles in MDD, our study aimed to examine which susceptible gene expression interacts with various dimensions of stressors in the MDD risk among a large population cohort.

METHODS

Data analyzed were from a longitudinal community-based cohort from Southwest Montreal, Canada (N = 1083). Latent profile models were used to identify distinct patterns of stressors for the study cohort. A transcriptome-wide association study (TWAS) method was performed to examine the interactive effects of three dimensions of stressors (threat, deprivation, and cumulative lifetime stress) and gene expression on the MDD risk in a total of 48 tissues from GTEx. Additional analyses were also conducted to further explore and specify these associations including colocalization, and fine-mapping analyses, in addition to enrichment analysis investigations based on TWAS.

RESULTS

We identified 3321 genes linked to MDD at the nominal p-value <0.05 and found that different patterns of stressors can amplify the genetic susceptibility to MDD. We also observed specific genes and pathways that interacted with deprivation and cumulative lifetime stressors, particularly in specific brain tissues including basal ganglia, prefrontal cortex, brain amygdala, brain cerebellum, brain cortex, and the whole blood. Colocalization analysis also identified these genes as having a high probability of sharing MDD causal variants.

LIMITATIONS

The study cohort was composed exclusively of individuals of Caucasians, which restricts the generalizability of the findings to other ethnic population groups.

CONCLUSIONS

The findings of the study unveiled significant interactions between potential tissue-specific gene expression × stressors in the MDD risk and shed light on the intricate etiological attributes of gene expression and specific stressors across the lifespan in MDD. These genetic and environmental attributes in MDD corroborate the vulnerability-stress theory and direct future stress research to have a closer examination of genetic predisposition and potential involvements of omics studies to specify the intricate relationships between genes and stressful environments.

Investigating the shared genetic architecture between depression and subcortical volumes

Depression, a widespread and highly heritable mental health condition, profoundly affects millions of individuals worldwide. Neuroimaging studies have consistently revealed volumetric abnormalities in subcortical structures associated with depression. However, the genetic underpinnings shared between depression and subcortical volumes remain inadequately understood. Here, we investigate the extent of polygenic overlap using the bivariate causal mixture model (MiXeR), leveraging summary statistics from the largest genome-wide association studies for depression (N = 674,452) and 14 subcortical volumetric phenotypes (N = 33,224). Additionally, we identify shared genomic loci through conditional/conjunctional FDR analyses. MiXeR shows that subcortical volumetric traits share a substantial proportion of genetic variants with depression, with 44 distinct shared loci identified by subsequent conjunctional FDR analysis. These shared loci are predominantly located in intronic regions (58.7%) and non-coding RNA intronic regions (25.4%). The 269 protein-coding genes mapped by these shared loci exhibit specific developmental trajectories, with the expression level of 55 genes linked to both depression and subcortical volumes, and 30 genes linked to cognitive abilities and behavioral symptoms. These findings highlight a shared genetic architecture between depression and subcortical volumetric phenotypes, enriching our understanding of the neurobiological underpinnings of depression.

A Single-Nucleus Transcriptome-Wide Association Study Implicates Novel Genes in Depression Pathogenesis

BACKGROUND

Depression, a common psychiatric illness and global public health problem, remains poorly understood across different life stages, which hampers the development of novel treatments.

METHODS

To identify new candidate genes for therapeutic development, we performed differential gene expression analysis of single-nucleus RNA sequencing data from the dorsolateral prefrontal cortex of older adults (n = 424) in relation to antemortem depressive symptoms. Additionally, we integrated genome-wide association study results for depression (n = 500,199) along with genetic tools for inferring the expression of 14,048 unique genes in 7 cell types and 52 cell subtypes to perform a transcriptome-wide association study of depression followed by Mendelian randomization.

RESULTS

Our single-nucleus transcriptome-wide association study analysis identified 68 candidate genes for depression and showed the greatest number being in excitatory and inhibitory neurons. Of the 68 genes, 53 were novel compared to previous studies. Notably, gene expression in different neuronal subtypes had varying effects on depression risk. Traits with high genetic correlations with depression, such as neuroticism, shared more transcriptome-wide association study genes than traits that were not highly correlated with depression. Complementing these analyses, differential gene expression analysis across 52 neocortical cell subtypes showed that genes such as KCNN2, SCAI, WASF3, and SOCS6 were associated with late-life depressive symptoms in specific cell subtypes.

CONCLUSIONS

These 2 sets of analyses illustrate the utility of large single-nucleus RNA sequencing data both to uncover genes whose expression is altered in specific cell subtypes in the context of depressive symptoms and to enhance the interpretation of well-powered genome-wide association studies so that we can prioritize specific susceptibility genes for further analysis and therapeutic development.

Mitochondrial proteases modulate mitochondrial stress signalling and cellular homeostasis in health and disease

Maintenance of mitochondrial homeostasis requires a plethora of coordinated quality control and adaptations' mechanisms in which mitochondrial proteases play a key role. Their activation or loss of function reverberate beyond local mitochondrial biochemical and metabolic remodelling into coordinated cellular pathways and stress responses that feedback onto the mitochondrial functionality and adaptability. Mitochondrial proteolysis modulates molecular and organellar quality control, metabolic adaptations, lipid homeostasis and regulates transcriptional stress responses. Defective mitochondrial proteolysis results in disease conditions most notably, mitochondrial diseases, neurodegeneration and cancer. Here, it will be discussed how mitochondrial proteases and mitochondria stress signalling impact cellular homeostasis and determine the cellular decision to survive or die, how these processes may impact disease etiopathology, and how modulation of proteolysis may offer novel therapeutic strategies.

Investigating neonatal health risk variables through cell-type specific methylome-wide association studies

Adverse neonatal outcomes are a prevailing risk factor for both short- and long-term mortality and morbidity in infants. Given the importance of these outcomes, refining their assessment is paramount for improving prevention and care. Here we aim to enhance the assessment of these often correlated and multifaceted neonatal outcomes. To achieve this, we employ factor analysis to identify common and unique effects and further confirm these effects using criterion-related validity testing. This validation leverages methylome-wide profiles from neonatal blood. Specifically, we investigate nine neonatal health risk variables, including gestational age, Apgar score, three indicators of body size, jaundice, birth diagnosis, maternal preeclampsia, and maternal age. The methylomic profiles used for this research capture data from nearly all 28 million methylation sites in human blood, derived from the blood spot collected from 333 neonates, within 72 h post-birth. Our factor analysis revealed two common factors, size factor, that captured the shared effects of weight, head size, height, and gestational age and disease factor capturing the orthogonal shared effects of gestational age, combined with jaundice and birth diagnosis. To minimize false positives in the validation studies, validation was limited to variables with significant cumulative association as estimated through an in-sample replication procedure. This screening resulted in that the two common factors and the unique effects for gestational age, jaundice and Apgar were further investigated with full-scale cell-type specific methylome-wide association analyses. Highly significant, cell-type specific, associations were detected for both common effect factors and for Apgar. Gene Ontology analyses revealed multiple significant biologically relevant terms for the five fully investigated neonatal health risk variables. Given the established links between adverse neonatal outcomes and both immediate and long-term health, the distinct factor effects (representing the common and unique effects of the risk variables) and their biological profiles confirmed in our work, suggest their potential role as clinical biomarkers for assessing health risks and enhancing personalized care.

A novel multi-omics mendelian randomization method for gene set enrichment and its application to psychiatric disorders

Genome-wide association studies (GWAS) of psychiatric disorders (PD) yield numerous loci with significant signals, but often do not implicate specific genes. Because GWAS risk loci are enriched in expression/protein/methylation quantitative loci (e/p/mQTL, hereafter xQTL), transcriptome/proteome/methylome-wide association studies (T/P/MWAS, hereafter XWAS) that integrate xQTL and GWAS information, can link GWAS signals to effects on specific genes. To further increase detection power, gene signals are aggregated within relevant gene sets (GS) by performing gene set enrichment (GSE) analyses. Often GSE methods test for enrichment of "signal" genes in curated GS while overlooking their linkage disequilibrium (LD) structure, allowing for the possibility of increased false positive rates. Moreover, no GSE tool uses xQTL information to perform mendelian randomization (MR) analysis. To make causal inference on association between PD and GS, we develop a novel MR GSE (MR-GSE) procedure. First, we generate a "synthetic" GWAS for each MSigDB GS by aggregating summary statistics for x-level (mRNA, protein or DNA methylation (DNAm) levels) from the largest xQTL studies available) of genes in a GS. Second, we use synthetic GS GWAS as exposure in a generalized summary-data-based-MR analysis of complex trait outcomes. We applied MR-GSE to GWAS of nine important PD. When applied to the underpowered opioid use disorder GWAS, none of the four analyses yielded any signals, which suggests a good control of false positive rates. For other PD, MR-GSE greatly increased the detection of GO terms signals (2,594) when compared to the commonly used (non-MR) GSE method (286). Some of the findings might be easier to adapt for treatment, e.g., our analyses suggest modest positive effects for supplementation with certain vitamins and/or omega-3 for schizophrenia, bipolar and major depression disorder patients. Similar to other MR methods, when applying MR-GSE researchers should be mindful of the confounding effects of horizontal pleiotropy on statistical inference.

Investigating the Mechanism of Chufan Yishen Formula in Treating Depression through Network Pharmacology and Experimental Verification

Objective: To investigate the antidepressant effect and potential mechanism of the Chufan Yishen Formula (CFYS) through network pharmacology, molecular docking, and experimental verification. Methods: The active ingredients and their target genes of CFYS were identified through Traditional Chinese Medicine Systems Pharmacology (TCMSP) and TCM-ID. We obtained the differentially expressed genes in patients with depression from the GEO database and screened out the genes intersecting with the target genes of CFYS to construct the PPI network. The key pathways were selected through STRING and KEGG. Then, molecular docking and experimental verification were performed. Results: A total of 113 effective components and 195 target genes were obtained. After intersecting the target genes with the differentially expressed genes in patients with depression, we obtained 37 differential target genes, among which HMOX1, VEGFA, etc., were the key genes. After enriching the differential target genes by KEGG, we found that the "chemical carcinogenesis-reactive oxygen species" pathway was the key pathway for the CFYS antidepressant effect. Besides, VEGFA might be a key marker for depression. Experimental verification found that CFYS could significantly improve the behavioral indicators of rats with depression models, including improving the antioxidant enzyme activity and increasing VEGFA levels. The results are consistent with the network pharmacology analysis. Conclusions: CFYS treatment for depression is a multicomponent, multitarget, and multipathway complex process, which may mainly exert an antidepressant effect by improving the neuron antioxidant stress response and regulating VEGFA levels.

Genome-wide meta-analysis, functional genomics and integrative analyses implicate new risk genes and therapeutic targets for anxiety disorders

Anxiety disorders are the most prevalent mental disorders. However, the genetic etiology of anxiety disorders remains largely unknown. Here we conducted a genome-wide meta-analysis on anxiety disorders by including 74,973 (28,392 proxy) cases and 400,243 (146,771 proxy) controls. We identified 14 risk loci, including 10 new associations near CNTNAP5, MAP2, RAB9BP1, BTN1A1, PRR16, PCLO, PTPRD, FARP1, CDH2 and RAB27B. Functional genomics and fine-mapping pinpointed the potential causal variants, and expression quantitative trait loci analysis revealed the potential target genes regulated by the risk variants. Integrative analyses, including transcriptome-wide association study, proteome-wide association study and colocalization analyses, prioritized potential causal genes (including CTNND1 and RAB27B). Evidence from multiple analyses revealed possibly causal genes, including RAB27B, BTN3A2, PCLO and CTNND1. Finally, we showed that Ctnnd1 knockdown affected dendritic spine density and resulted in anxiety-like behaviours in mice, revealing the potential role of CTNND1 in anxiety disorders. Our study identified new risk loci, potential causal variants and genes for anxiety disorders, providing insights into the genetic architecture of anxiety disorders and potential therapeutic targets.

A statistical method for image-mediated association studies discovers genes and pathways associated with four brain disorders

Brain imaging and genomics are critical tools enabling characterization of the genetic basis of brain disorders. However, imaging large cohorts is expensive and may be unavailable for legacy datasets used for genome-wide association studies (GWASs). Using an integrated feature selection/aggregation model, we developed an image-mediated association study (IMAS), which utilizes borrowed imaging/genomics data to conduct association mapping in legacy GWAS cohorts. By leveraging the UK Biobank image-derived phenotypes (IDPs), the IMAS discovered genetic bases underlying four neuropsychiatric disorders and verified them by analyzing annotations, pathways, and expression quantitative trait loci (eQTLs). A cerebellar-mediated mechanism was identified to be common to the four disorders. Simulations show that, if the goal is identifying genetic risk, our IMAS is more powerful than a hypothetical protocol in which the imaging results were available in the GWAS dataset. This implies the feasibility of reanalyzing legacy GWAS datasets without conducting additional imaging, yielding cost savings for integrated analysis of genetics and imaging.

Whole Blood-based Transcriptional Risk Score for Nonobese Type 2 Diabetes Predicts Dynamic Changes in Glucose Metabolism

CONTEXT

The performance of peripheral blood transcriptional markers in evaluating risk of type 2 diabetes (T2D) with normal body mass index (BMI) is unknown.

OBJECTIVE

We developed a whole blood-based transcriptional risk score (wb-TRS) for nonobese T2D and assessed its contributions on disease risk and dynamic changes in glucose metabolism.

METHODS

Using a community-based cohort with blood transcriptome data, we developed the wb-TRS in 1105 participants aged ≥40 years who maintained a normal BMI for up to 10 years, and we validated the wb-TRS in an external dataset. Potential biological significance was explored.

RESULTS

The wb-TRS included 144 gene transcripts. Compared to the lowest tertile, wb-TRS in tertile 3 was associated with 8.91-fold (95% CI, 3.53-22.5) higher risk and each 1-unit increment was associated with 2.63-fold (95% CI, 1.87-3.68) higher risk of nonobese T2D. Furthermore, baseline wb-TRS significantly associated with dynamic changes in average, daytime, nighttime, and 24-hour glucose, HbA1c values, and area under the curve of glucose measured by continuous glucose monitoring over 6 months of intervention. The wb-TRS improved the prediction performance for nonobese T2D, combined with fasting glucose, triglycerides, and demographic and anthropometric parameters. Multi-contrast gene set enrichment (Mitch) analysis implicated oxidative phosphorylation, mTORC1 signaling, and cholesterol metabolism involved in nonobese T2D pathogenesis.

CONCLUSION

A whole blood-based nonobese T2D-associated transcriptional risk score was validated to predict dynamic changes in glucose metabolism. These findings suggested several biological pathways involved in the pathogenesis of nonobese T2D.

Epigenetic modification related to cognitive changes during a cognitive training intervention in depression

BACKGROUND

DNA methylation as a biomarker is well suited to investigate dynamic processes, such as symptom improvement. For this study we focus on epigenomic state or trait markers as early signatures of cognitive improvement in individuals receiving a cognitive intervention. We performed a first epigenome-wide association study (EWAS) on patients with cognitive dysfunction in depression comparing those with vs without cognitive dysfunction and those cognitively improving vs non-improving following a cognitive intervention.

METHOD

Data from a randomized controlled trial (RCT) were used for this analysis, where cognitive function of 112 patients randomly assigned to a personalized cognitive intervention was compared to standard cognitive treatment. Cognition was measured for this study using the four cognitive tasks from the THINC-it battery. We compared individuals with cognitive impairment with individuals without cognitive impairment at baseline and after a cognitive intervention of 8 weeks. Blood for DNA methylation analysis (Illumina Infinium MethylationEPIC 850 k BeadChip) was collected at baseline and 8 weeks into the treatment. For the baseline analysis, after quality control, the final sample comprised 90 individuals, and analyses at week 8 were performed on 84 individuals. Data cleaning, quality control, and differential methylation analysis of DNA methylation data was performed using the RnBeads package (R). Analyses were corrected for gender, age, depression score (MADRS), reported years of education, height and weight, as well as surrogate variables estimated by the pipeline used. The within-individual paired longitudinal analysis was performed using Welch's t-test.

RESULTS

Analyses at baseline and at week 8 did not show any genome-wide significant CpGs (p < 5 × 10-8) comparing patients with and without cognitive impairment. The most significant result in the baseline analysis comparing the groups with and without cognitive impairment at baseline is located in an open Sea region with predominantly regulatory qualities (cg10962945; 6.61 × 10-7). The most significant CpG at 8 weeks was also located in open sea, though in exon 13 of the NTRK2-gene, linked to the BDNF pathway (cg13620631, 5.56 × 10-7). Finally, a within-individual paired longitudinal analysis with only patients that show improved cognitive function over time was performed, showing 65 CpGs that overlapped between the 1% most significant of this analysis and the 1% most significant CpGs from the cross-sectional analysis at 8 weeks.

CONCLUSION

Our result suggest that DNA methylation can be suitable to capture early signs of treatment response of a cognitive intervention in depression. In our layered approach we could capture dynamics that can help differentiate between biological trait and state markers of cognitive function in depression. Despite not being genome-wide significant, the CpG locations returned by our analysis comparing patients with and without cognitive impairment, are in line with prior knowledge on pathways and genes relevant for depression treatment and cognition.

TWAS revealed significant causal loci for milk production and its composition in Murrah buffaloes

Milk yield is the most complex trait in dairy animals, and mapping all causal variants even with smallest effect sizes has been difficult with the genome-wide association study (GWAS) sample sizes available in geographical regions with small livestock holdings such as Indian sub-continent. However, Transcriptome-wide association studies (TWAS) could serve as an alternate for fine mapping of expression quantitative trait loci (eQTLs). This is a maiden attempt to identify milk production and its composition related genes using TWAS in Murrah buffaloes (Bubalus bubalis). TWAS was conducted on a test (N = 136) set of Murrah buffaloes genotyped through ddRAD sequencing. Their gene expression level was predicted using reference (N = 8) animals having both genotype and mammary epithelial cell (MEC) transcriptome information. Gene expression prediction was performed using Elastic-Net and Dirichlet Process Regression (DPR) model with fivefold cross-validation and without any cross-validation. DPR model without cross-validation predicted 80.92% of the total genes in the test group of Murrah buffaloes which was highest compared to other methods. TWAS in test individuals based on predicted gene expression, identified a significant association of one unique gene for Fat%, and two for SNF% at Bonferroni corrected threshold. The false discovery rates (FDR) corrected P-values of the top ten SNPs identified through GWAS were comparatively higher than TWAS. Gene ontology of TWAS-identified genes was performed to understand the function of these genes, it was revealed that milk production and composition genes were mainly involved in Relaxin, AMPK, and JAK-STAT signaling pathway, along with CCRI, and several key metabolic processes. The present study indicates that TWAS offers a lower false discovery rate and higher significant hits than GWAS for milk production and its composition traits. Hence, it is concluded that TWAS can be effectively used to identify genes and cis-SNPs in a population, which can be used for fabricating a low-density genomic chip for predicting milk production in Murrah buffaloes.

ENAH-202 promotes cancer progression in oral squamous cell carcinoma by regulating ZNF502/VIM axis

BACKGROUND

We aimed to demonstrate the regulatory effect of long non-coding RNA (lncRNA) ENAH-202 on oral squamous cell carcinoma (OSCC) development as well as its molecular mechanism.

METHODS

We detected ENAH-202 expression in OSCC tissues and cell lines by quantitative real-time PCR (qPCR). The biological function of ENAH-202 was assessed in vitro and in vivo using CCK-8, colony formation assays, transwell assays, xenograft formation, and tail vein injection. The further molecular mechanism by which ENAH-202 promoted OSCC progression was identified using RNA pull-down, LS-MS/MS analysis, RNA immunoprecipitation (RIP), and chromatin immunoprecipitation (ChIP) assays.

RESULTS

ENAH-202 was significantly upregulated in OSCC tissues and cells. ENAH-202 promoted OSCC cell proliferation, migration, and invasion in vitro and in vivo. The expression of enabled homolog (ENAH) and epithelial-to-mesenchymal transition (EMT)-related proteins was changed with the expression of ENAH-202. Moreover, ENAH-202 promoted the transcription of Vimentin (VIM) by binding with ZNF502, which can help ENAH-202 promote OSCC progression.

CONCLUSIONS

ENAH-202 facilitated OSCC cell proliferation and metastasis by regulating ZNF502/VIM axis, which played an important role in OSCC progression.

Identification of region-specific splicing QTLs in human hippocampal tissue and its distinctive role in brain disorders

Alternative splicing (AS) regulation has an essential role in complex diseases. However, the AS profiles in the hippocampal (HIPPO) region of human brain are underexplored. Here, we investigated cis-acting sQTLs of HIPPO region in 264 samples and identified thousands of significant sQTLs. By enrichment analysis and functional characterization of these sQTLs, we found that the HIPPO sQTLs were enriched among histone-marked regions, transcription factors binding sites, RNA binding proteins sites, and brain disorders-associated loci. Comparative analyses with the dorsolateral prefrontal cortex revealed the importance of AS regulation in HIPPO (rg = 0.87). Furthermore, we performed a transcriptome-wide association study of Alzheimer's disease and identified 16 significant genes whose genetically regulated splicing levels may have a causal role in Alzheimer. Overall, our study improves our knowledge of the transcriptome gene regulation in the HIPPO region and provides novel insights into elucidating the pathogenesis of potential genes associated with brain disorders.

Integrative bioinformatics and artificial intelligence analyses of transcriptomics data identified genes associated with major depressive disorders including NRG1

Major depressive disorder (MDD) is a common mental disorder and is amongst the most prevalent psychiatric disorders. MDD remains challenging to diagnose and predict its onset due to its heterogeneous phenotype and complex etiology. Hence, early detection using diagnostic biomarkers is critical for rapid intervention. In this study, a mixture of AI and bioinformatics were used to mine transcriptomic data from publicly available datasets including 170 MDD patients and 121 healthy controls. Bioinformatics analysis using gene set enrichment analysis (GSEA) and machine learning (ML) algorithms were applied. The GSEA revealed that differentially expressed genes in MDD patients are mainly enriched in pathways related to immune response, inflammatory response, neurodegeneration pathways and cerebellar atrophy pathways. Feature selection methods and ML provided predicted models based on MDD-altered genes with ≥75% of accuracy. The integrative analysis between the bioinformatics and ML approaches identified ten key MDD-related biomarkers including NRG1, CEACAM8, CLEC12B, DEFA4, HP, LCN2, OLFM4, SERPING1, TCN1 and THBS1. Among them, NRG1, active in synaptic plasticity and neurotransmission, was the most robust and reliable to distinguish between MDD patients and healthy controls amongst independent external datasets consisting of a mixture of populations. Further evaluation using saliva samples from an independent cohort of MDD and healthy individuals confirmed the upregulation of NRG1 in patients with MDD compared to healthy controls. Functional mapping to the human brain regions showed NRG1 to have high expression in the main subcortical limbic brain regions implicated in depression. In conclusion, integrative bioinformatics and ML approaches identified putative non-invasive diagnostic MDD-related biomarkers panel for the onset of depression.

Whole-genome sequencing analysis of suicide deaths integrating brain-regulatory eQTLs data to identify risk loci and genes

Recent large-scale genome-wide association studies (GWAS) have started to identify potential genetic risk loci associated with risk of suicide; however, a large portion of suicide-associated genetic factors affecting gene expression remain elusive. Dysregulated gene expression, not assessed by GWAS, may play a significant role in increasing the risk of suicide death. We performed the first comprehensive genomic association analysis prioritizing brain expression quantitative trait loci (eQTLs) within regulatory regions in suicide deaths from the Utah Suicide Genetic Risk Study (USGRS). 440,324 brain-regulatory eQTLs were obtained by integrating brain eQTLs, histone modification ChIP-seq, ATAC-seq, DNase-seq, and Hi-C results from publicly available data. Subsequent genomic analyses were conducted in whole-genome sequencing (WGS) data from 986 suicide deaths of non-Finnish European (NFE) ancestry and 415 ancestrally matched controls. Additional independent USGRS suicide deaths with genotyping array data (n = 4657) and controls from the Genome Aggregation Database were explored for WGS result replication. One significant eQTL locus, rs926308 (p = 3.24e-06), was identified. The rs926308-T is associated with lower expression of RFPL3S, a gene important for neocortex development and implicated in arousal. Gene-based analyses performed using Sherlock Bayesian statistical integrative analysis also detected 20 genes with expression changes that may contribute to suicide risk. From analyzing publicly available transcriptomic data, ten of these genes have previous evidence of differential expression in suicide death or in psychiatric disorders that may be associated with suicide, including schizophrenia and autism (ZNF501, ZNF502, CNN3, IGF1R, KLHL36, NBL1, PDCD6IP, SNX19, BCAP29, and ARSA). Electronic health records (EHR) data was further merged to evaluate if there were clinically relevant subsets of suicide deaths associated with genetic variants. In summary, our study identified one risk locus and ten genes associated with suicide risk via gene expression, providing new insight into possible genetic and molecular mechanisms leading to suicide.

Investigation of the shared molecular mechanisms and hub genes between myocardial infarction and depression

Background

The pathogenesis of myocardial infarction complicating depression is still not fully understood. Bioinformatics is an effective method to study the shared pathogenesis of multiple diseases and has important application value in myocardial infarction complicating depression.

Methods

The differentially expressed genes (DEGs) between control group and myocardial infarction group (M-DEGs), control group and depression group (D-DEGs) were identified in the training set. M-DEGs and D-DEGs were intersected to obtain DEGs shared by the two diseases (S-DEGs). The GO, KEGG, GSEA and correlation analysis were conducted to analyze the function of DEGs. The biological function differences of myocardial infarction and depression were analyzed by GSVA and immune cell infiltration analysis. Four machine learning methods, nomogram, ROC analysis, calibration curve and decision curve were conducted to identify hub S-DEGs and predict depression risk. The unsupervised cluster analysis was constructed to identify myocardial infarction molecular subtype clusters based on hub S-DEGs. Finally, the value of these genes was verified in the validation set, and blood samples were collected for RT-qPCR experiments to further verify the changes in expression levels of these genes in myocardial infarction and depression.

Results

A total of 803 M-DEGs, 214 D-DEGs, 13 S-DEGs and 6 hub S-DEGs (CD24, CSTA, EXTL3, RPS7, SLC25A5 and ZMAT3) were obtained in the training set and they were all involved in immune inflammatory response. The GSVA and immune cell infiltration analysis results also suggested that immune inflammation may be the shared pathogenesis of myocardial infarction and depression. The diagnostic models based on 6 hub S-DEGs found that these genes showed satisfactory combined diagnostic performance for depression. Then, two molecular subtypes clusters of myocardial infarction were identified, many differences in immune inflammation related-biological functions were found between them, and the hub S-DEGs had satisfactory molecular subtypes identification performance. Finally, the analysis results of the validation set further confirmed the value of these hub genes, and the RT-qPCR results of blood samples further confirmed the expression levels of these hub genes in myocardial infarction and depression.

Conclusion

Immune inflammation may be the shared pathogenesis of myocardial infarction and depression. Meanwhile, hub S-DEGs may be potential biomarkers for the diagnosis and molecular subtype identification of myocardial infarction and depression.

Integrating genetics and transcriptomics to study major depressive disorder: a conceptual framework, bioinformatic approaches, and recent findings

Major depressive disorder (MDD) is a complex and heterogeneous psychiatric syndrome with genetic and environmental influences. In addition to neuroanatomical and circuit-level disturbances, dysregulation of the brain transcriptome is a key phenotypic signature of MDD. Postmortem brain gene expression data are uniquely valuable resources for identifying this signature and key genomic drivers in human depression; however, the scarcity of brain tissue limits our capacity to observe the dynamic transcriptional landscape of MDD. It is therefore crucial to explore and integrate depression and stress transcriptomic data from numerous, complementary perspectives to construct a richer understanding of the pathophysiology of depression. In this review, we discuss multiple approaches for exploring the brain transcriptome reflecting dynamic stages of MDD: predisposition, onset, and illness. We next highlight bioinformatic approaches for hypothesis-free, genome-wide analyses of genomic and transcriptomic data and their integration. Last, we summarize the findings of recent genetic and transcriptomic studies within this conceptual framework.

The Melanocortin System: A Promising Target for the Development of New Antidepressant Drugs

Major depression is one of the most prevalent mental disorders, causing significant human suffering and socioeconomic loss. Since conventional antidepressants are not sufficiently effective, there is an urgent need to develop new antidepressant medications. Despite marked advances in the neurobiology of depression, the etiology and pathophysiology of this disease remain poorly understood. Classical and newer hypotheses of depression suggest that an imbalance of brain monoamines, dysregulation of the hypothalamic-pituitary-adrenal axis (HPAA) and immune system, or impaired hippocampal neurogenesis and neurotrophic factors pathways are cause of depression. It is assumed that conventional antidepressants improve these closely related disturbances. The purpose of this review was to discuss the possibility of affecting these disturbances by targeting the melanocortin system, which includes adrenocorticotropic hormone-activated receptors and their peptide ligands (melanocortins). The melanocortin system is involved in the regulation of various processes in the brain and periphery. Melanocortins, including peripherally administered non-corticotropic agonists, regulate HPAA activity, exhibit anti-inflammatory effects, stimulate the levels of neurotrophic factors, and enhance hippocampal neurogenesis and neurotransmission. Therefore, endogenous melanocortins and their analogs are able to complexly affect the functioning of those body’s systems that are closely related to depression and the effects of antidepressants, thereby demonstrating a promising antidepressant potential.

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.

Depression recognition using a proposed speech chain model fusing speech production and perception features

BACKGROUND

Increasing depression patients puts great pressure on clinical diagnosis. Audio-based diagnosis is a helpful auxiliary tool for early mass screening. However, current methods consider only speech perception features, ignoring patients' vocal tract changes, which may partly result in the poor recognition.

METHODS

This work proposes a novel machine speech chain model for depression recognition (MSCDR) that can capture text-independent depressive speech representation from the speaker's mouth to the listener's ear to improve recognition performance. In the proposed MSCDR, linear predictive coding (LPC) and Mel-frequency cepstral coefficients (MFCC) features are extracted to describe the processes of speech generation and of speech perception, respectively. Then, a one-dimensional convolutional neural network and a long short-term memory network sequentially capture intra- and inter-segment dynamic depressive features for classification.

RESULTS

We tested the MSCDR on two public datasets with different languages and paradigms, namely, the Distress Analysis Interview Corpus-Wizard of Oz and the Multi-modal Open Dataset for Mental-disorder Analysis. The accuracy of the MSCDR on the two datasets was 0.77 and 0.86, and the average F1 score was 0.75 and 0.86, which were better than the other existing methods. This improvement reveals the complementarity of speech production and perception features in carrying depressive information.

LIMITATIONS

The sample size was relatively small, which may limit the application in clinical translation to some extent.

CONCLUSION

This experiment proves the good generalization ability and superiority of the proposed MSCDR and suggests that the vocal tract changes in patients with depression deserve attention for audio-based depression diagnosis.

Variable number tandem repeats (VNTRs) as modifiers of breast cancer risk in carriers of BRCA1 185delAG

Despite substantial efforts in identifying both rare and common variants affecting disease risk, in the majority of diseases, a large proportion of unexplained genetic risk remains. We propose that variable number tandem repeats (VNTRs) may explain a proportion of the missing genetic risk. Herein, in a pilot study with a retrospective cohort design, we tested whether VNTRs are causal modifiers of breast cancer risk in 347 female carriers of the BRCA1 185delAG pathogenic variant, an important group given their high risk of developing breast cancer. We performed targeted-capture to sequence VNTRs, called genotypes with adVNTR, tested the association of VNTRs and breast cancer risk using Cox regression models, and estimated the effect size using a retrospective likelihood approach. Of 303 VNTRs that passed quality control checks, 4 VNTRs were significantly associated with risk to develop breast cancer at false discovery rate [FDR] < 0.05 and an additional 4 VNTRs had FDR < 0.25. After determining the specific risk alleles, there was a significantly earlier age at diagnosis of breast cancer in carriers of the risk alleles compared to those without the risk alleles for seven of eight VNTRs. One example is a VNTR in exon 2 of LINC01973 with a per-allele hazard ratio of 1.58 (1.07-2.33) and 5.28 (2.79-9.99) for the homozygous risk-allele genotype. Results from this first systematic study of VNTRs demonstrate that VNTRs may explain a proportion of the unexplained genetic risk for breast cancer.

Transcriptome-wide association analyses identify an association between ARL14EP and polycystic ovary syndrome

Polycystic ovary syndrome (PCOS) is a common endocrine disorder, which is accompanied by a variety of comorbidities including metabolic, reproductive, and psychiatric disorders. Genome-wide association studies have identified several genetic variants that are associated with PCOS. However, these variants often occur outside of coding regions and require further investigation to understand their contribution to PCOS. A transcriptome-wide association study (TWAS) was performed to uncover heritable gene expression profiles that are associated with PCOS in two independent cohorts. Causal gene prioritization was subsequently performed and expression of genes prioritized through these analyses was examined in 49 PCOS patients and 30 controls. TWAS analyses revealed that increased expression of ARL14EP was significantly associated with PCOS risk in the discovery (P = 1.6 × 10^-6) and replication cohorts (P = 2.0 × 10^-13). Gene prioritization pipelines provided further evidence that ARL14EP is the most likely causal gene at this locus. ARL14EP gene expression was shown to be significantly different between PCOS cases and controls, after adjusting for body mass index, age and testosterone levels (P = 1.2 × 10^-13). This study has provided evidence for the role of ARL14EP in PCOS. Given that ARL14EP has been reported to play an important role in chromatin remodeling, variants affecting the expression of ARL14EP may also affect the expression of other genes that contribute to PCOS pathogenesis.

Multi-omics data integration methods and their applications in psychiatric disorders

To study mental illness and health, in the past researchers have often broken down their complexity into individual subsystems (e.g., genomics, transcriptomics, proteomics, clinical data) and explored the components independently. Technological advancements and decreasing costs of high throughput sequencing has led to an unprecedented increase in data generation. Furthermore, over the years it has become increasingly clear that these subsystems do not act in isolation but instead interact with each other to drive mental illness and health. Consequently, individual subsystems are now analysed jointly to promote a holistic understanding of the underlying biological complexity of health and disease. Complementing the increasing data availability, current research is geared towards developing novel methods that can efficiently combine the information rich multi-omics data to discover biologically meaningful biomarkers for diagnosis, treatment, and prognosis. However, clinical translation of the research is still challenging. In this review, we summarise conventional and state-of-the-art statistical and machine learning approaches for discovery of biomarker, diagnosis, as well as outcome and treatment response prediction through integrating multi-omics and clinical data. In addition, we describe the role of biological model systems and in silico multi-omics model designs in clinical translation of psychiatric research from bench to bedside. Finally, we discuss the current challenges and explore the application of multi-omics integration in future psychiatric research. The review provides a structured overview and latest updates in the field of multi-omics in psychiatry.

Personality traits as mediators in the association between SIRT1 rs12415800 polymorphism and depressive symptoms among Chinese college students

Background

Previous research has linked polymorphisms in the SIRT1 gene to depressive symptoms, particularly in Chinese individuals. However, it is not clear how personality traits may contribute to this association.

Methods

To explore the potential mediating effect of personality traits, we utilized a mediation model to examine the relationship between the SIRT1 rs12415800 polymorphism and depressive symptoms in 787 Chinese college students. Depressive symptoms were assessed using the Center for Epidemiologic Studies Depression (CES-D) scale, while personality traits were measured using the Big Five Inventory (BFI).

Results

Our analysis indicated a significant association between the SIRT1 rs12415800 polymorphism and depressive symptoms, with this relationship partially mediated by the personality traits of neuroticism and conscientiousness. Specifically, individuals who were heterozygous for the rs12415800 polymorphism and had higher levels of conscientiousness were less likely to experience depressive symptoms. Conversely, those who were homozygous for the rs12415800 polymorphism and had higher levels of neuroticism were more likely to experience depressive symptoms.

Conclusion

Our results suggest that personality traits, particularly neuroticism and conscientiousness, may play a critical role in the association between the SIRT1 rs12415800 polymorphism and depressive symptoms among Chinese college students. These findings highlight the importance of considering both genetic factors and personality traits when exploring the etiology of depressive symptoms in this population.

Benchmarking post-GWAS analysis tools in major depression: Challenges and implications

Our knowledge of complex disorders has increased in the last years thanks to the identification of genetic variants (GVs) significantly associated with disease phenotypes by genome-wide association studies (GWAS). However, we do not understand yet how these GVs functionally impact disease pathogenesis or their underlying biological mechanisms. Among the multiple post-GWAS methods available, fine-mapping and colocalization approaches are commonly used to identify causal GVs, meaning those with a biological effect on the trait, and their functional effects. Despite the variety of post-GWAS tools available, there is no guideline for method eligibility or validity, even though these methods work under different assumptions when accounting for linkage disequilibrium and integrating molecular annotation data. Moreover, there is no benchmarking of the available tools. In this context, we have applied two different fine-mapping and colocalization methods to the same GWAS on major depression (MD) and expression quantitative trait loci (eQTL) datasets. Our goal is to perform a systematic comparison of the results obtained by the different tools. To that end, we have evaluated their results at different levels: fine-mapped and colocalizing GVs, their target genes and tissue specificity according to gene expression information, as well as the biological processes in which they are involved. Our findings highlight the importance of fine-mapping as a key step for subsequent analysis. Notably, the colocalizing variants, altered genes and targeted tissues differed between methods, even regarding their biological implications. This contribution illustrates an important issue in post-GWAS analysis with relevant consequences on the use of GWAS results for elucidation of disease pathobiology, drug target prioritization and biomarker discovery.

Functional Genomics Analysis to Disentangle the Role of Genetic Variants in Major Depression

Understanding the molecular basis of major depression is critical for identifying new potential biomarkers and drug targets to alleviate its burden on society. Leveraging available GWAS data and functional genomic tools to assess regulatory variation could help explain the role of major depression-associated genetic variants in disease pathogenesis. We have conducted a fine-mapping analysis of genetic variants associated with major depression and applied a pipeline focused on gene expression regulation by using two complementary approaches: cis-eQTL colocalization analysis and alteration of transcription factor binding sites. The fine-mapping process uncovered putative causally associated variants whose proximal genes were linked with major depression pathophysiology. Four colocalizing genetic variants altered the expression of five genes, highlighting the role of SLC12A5 in neuronal chlorine homeostasis and MYRF in nervous system myelination and oligodendrocyte differentiation. The transcription factor binding analysis revealed the potential role of rs62259947 in modulating P4HTM expression by altering the YY1 binding site, altogether regulating hypoxia response. Overall, our pipeline could prioritize putative causal genetic variants in major depression. More importantly, it can be applied when only index genetic variants are available. Finally, the presented approach enabled the proposal of mechanistic hypotheses of these genetic variants and their role in disease pathogenesis.

dbBIP: a comprehensive bipolar disorder database for genetic research

Bipolar disorder (BIP) is one of the most common hereditary psychiatric disorders worldwide. Elucidating the genetic basis of BIP will play a pivotal role in mechanistic delineation. Genome-wide association studies (GWAS) have successfully reported multiple susceptibility loci conferring BIP risk, thus providing insight into the effects of its underlying pathobiology. However, difficulties remain in the extrication of important and biologically relevant data from genetic discoveries related to psychiatric disorders such as BIP. There is an urgent need for an integrated and comprehensive online database with unified access to genetic and multi-omics data for in-depth data mining. Here, we developed the dbBIP, a database for BIP genetic research based on published data. The dbBIP consists of several modules, i.e.: (i) single nucleotide polymorphism (SNP) module, containing large-scale GWAS genetic summary statistics and functional annotation information relevant to risk variants; (ii) gene module, containing BIP-related candidate risk genes from various sources and (iii) analysis module, providing a simple and user-friendly interface to analyze one's own data. We also conducted extensive analyses, including functional SNP annotation, integration (including summary-data-based Mendelian randomization and transcriptome-wide association studies), co-expression, gene expression, tissue expression, protein-protein interaction and brain expression quantitative trait loci analyses, thus shedding light on the genetic causes of BIP. Finally, we developed a graphical browser with powerful search tools to facilitate data navigation and access. The dbBIP provides a comprehensive resource for BIP genetic research as well as an integrated analysis platform for researchers and can be accessed online at http://dbbip.xialab.info. Database URL:  http://dbbip.xialab.info.

Transcriptome Profiling of the Dorsomedial Prefrontal Cortex in Suicide Victims

The default mode network (DMN) plays an outstanding role in psychiatric disorders. Still, gene expressional changes in its major component, the dorsomedial prefrontal cortex (DMPFC), have not been characterized. We used RNA sequencing in postmortem DMPFC samples to investigate suicide victims compared to control subjects. 1400 genes differed using log2FC > ±1 and adjusted p-value < 0.05 criteria between groups. Genes associated with depressive disorder, schizophrenia and impaired cognition were strongly overexpressed in top differentially expressed genes. Protein−protein interaction and co-expressional networks coupled with gene set enrichment analysis revealed that pathways related to cytokine receptor signaling were enriched in downregulated, while glutamatergic synaptic signaling upregulated genes in suicidal individuals. A validated differentially expressed gene, which is known to be associated with mGluR5, was the N-terminal EF-hand calcium-binding protein 2 (NECAB2). In situ hybridization histochemistry and immunohistochemistry proved that NECAB2 is expressed in two different types of inhibitory neurons located in layers II-IV and VI, respectively. Our results imply extensive gene expressional alterations in the DMPFC related to suicidal behavior. Some of these genes may contribute to the altered mental state and behavior of suicide victims.

Identifying causal genes for depression via integration of the proteome and transcriptome from brain and blood

Genome-wide association studies (GWASs) have identified numerous risk genes for depression. Nevertheless, genes crucial for understanding the molecular mechanisms of depression and effective antidepressant drug targets are largely unknown. Addressing this, we aimed to highlight potentially causal genes by systematically integrating the brain and blood protein and expression quantitative trait loci (QTL) data with a depression GWAS dataset via a statistical framework including Mendelian randomization (MR), Bayesian colocalization, and Steiger filtering analysis. In summary, we identified three candidate genes (TMEM106B, RAB27B, and GMPPB) based on brain data and two genes (TMEM106B and NEGR1) based on blood data with consistent robust evidence at both the protein and transcriptional levels. Furthermore, the protein-protein interaction (PPI) network provided new insights into the interaction between brain and blood in depression. Collectively, four genes (TMEM106B, RAB27B, GMPPB, and NEGR1) affect depression by influencing protein and gene expression level, which could guide future researches on candidate genes investigations in animal studies as well as prioritize antidepressant drug targets.

Transcriptome-based polygenic score links depression-related corticolimbic gene expression changes to sex-specific brain morphology and depression risk

Studies in post-mortem human brain tissue have associated major depressive disorder (MDD) with cortical transcriptomic changes, whose potential in vivo impact remains unexplored. To address this translational gap, we recently developed a transcriptome-based polygenic risk score (T-PRS) based on common functional variants capturing 'depression-like' shifts in cortical gene expression. Here, we used a non-clinical sample of young adults (n = 482, Duke Neurogenetics Study: 53% women; aged 19.8 ± 1.2 years) to map T-PRS onto brain morphology measures, including Freesurfer-derived subcortical volume, cortical thickness, surface area, and local gyrification index, as well as broad MDD risk, indexed by self-reported family history of depression. We conducted side-by-side comparisons with a PRS independently derived from a Psychiatric Genomics Consortium (PGC) MDD GWAS (PGC-PRS), and sought to link T-PRS with diagnosis and symptom severity directly in PGC-MDD participants (n = 29,340, 59% women; 12,923 MDD cases, 16,417 controls). T-PRS was associated with smaller amygdala volume in women (t = -3.478, p = 0.001) and lower prefrontal gyrification across sexes. In men, T-PRS was associated with hypergyrification in temporal and occipital regions. Prefrontal hypogyrification mediated a male-specific indirect link between T-PRS and familial depression (b = 0.005, p = 0.029). PGC-PRS was similarly associated with lower amygdala volume and cortical gyrification; however, both effects were male-specific and hypogyrification emerged in distinct parietal and temporo-occipital regions, unassociated with familial depression. In PGC-MDD, T-PRS did not predict diagnosis (OR = 1.007, 95% CI = [0.997-1.018]) but correlated with symptom severity in men (rho = 0.175, p = 7.957 × 10-4) in one cohort (N = 762, 48% men). Depression-like shifts in cortical gene expression have sex-specific effects on brain morphology and may contribute to broad depression vulnerability in men.

Proteome-wide Association Study Provides Insights Into the Genetic Component of Protein Abundance in Psychiatric Disorders

BACKGROUND

Genome-wide association studies have identified multiple risk variants for psychiatric disorders. Nevertheless, how the risk variants confer risk of psychiatric disorders remains largely unknown.

METHODS

We performed proteome-wide association studies to identify genes whose cis-regulated protein abundance change in the human brain were associated with psychiatric disorders.

RESULTS

By integrating genome-wide associations of four common psychiatric disorders and two independent brain proteomes (n = 376 and n = 152, respectively) from the dorsolateral prefrontal cortex, we identified 61 genes (including 48 genes for schizophrenia, 12 genes for bipolar disorder, 5 genes for depression, and 2 genes for attention-deficit/hyperactivity disorder) whose genetically regulated protein abundance levels were associated with risk of psychiatric disorders. Comparison with transcriptome-wide association studies identified 18 overlapping genes that showed significant associations with psychiatric disorders at both proteome-wide and transcriptome-wide levels, suggesting that genetic risk variants likely confer risk of psychiatric disorders by regulating messenger RNA expression and protein abundance of these genes.

CONCLUSIONS

Our study not only provides new insights into the genetic component of protein abundance in psychiatric disorders but also highlights several high-confidence risk proteins (including CNNM2 and CTNND1) for schizophrenia and depression. These high-confidence risk proteins represent promising therapeutic targets for future drug development.