Brain Proteome-Wide Association Study Identifies Candidate Genes that Regulate Protein Abundance Associated with Post-Traumatic Stress Disorder

Proteome-wide association study using cis and trans variants and applied to blood cell and lipid-related traits in the Women's Health Initiative study

In most Proteome-Wide Association Studies (PWAS), variants near the protein-coding gene (±1 Mb), also known as cis single nucleotide polymorphisms (SNPs), are used to predict protein levels, which are then tested for association with phenotypes. However, proteins can be regulated through variants outside of the cis region. An intermediate GWAS step to identify protein quantitative trait loci (pQTL) allows for the inclusion of trans SNPs outside the cis region in protein-level prediction models. Here, we assess the prediction of 540 proteins in 1002 individuals from the Women's Health Initiative (WHI), split equally into a GWAS set, an elastic net training set, and a testing set. We compared the testing r2 between measured and predicted protein levels using this proposed approach, to the testing r2 using only cis SNPs. The two methods usually resulted in similar testing r2, but some proteins showed a significant increase in testing r2 with our method. For example, for cartilage acidic protein 1, the testing r2 increased from 0.101 to 0.351. We also demonstrate reproducible findings for predicted protein association with lipid and blood cell traits in WHI participants without proteomics data and in UK Biobank utilizing our PWAS weights.

Genetic Basis of Stress-Related Neuropsychiatric Disorders

Stress exposure is one of the major risk factors for the development of different psychiatric and neurodegenerative diseases, leading to various (mal)adaptive alternations in different parts of the brain, and affecting both brain structure and function [...].

TrkB-dependent regulation of molecular signaling across septal cell types

The lateral septum (LS), a GABAergic structure located in the basal forebrain, is implicated in social behavior, learning, and memory. We previously demonstrated that expression of tropomyosin kinase receptor B (TrkB) in LS neurons is required for social novelty recognition. To better understand molecular mechanisms by which TrkB signaling controls behavior, we locally knocked down TrkB in LS and used bulk RNA-sequencing to identify changes in gene expression downstream of TrkB. TrkB knockdown induces upregulation of genes associated with inflammation and immune responses, and downregulation of genes associated with synaptic signaling and plasticity. Next, we generated one of the first atlases of molecular profiles for LS cell types using single nucleus RNA-sequencing (snRNA-seq). We identified markers for the septum broadly, and the LS specifically, as well as for all neuronal cell types. We then investigated whether the differentially expressed genes (DEGs) induced by TrkB knockdown map to specific LS cell types. Enrichment testing identified that downregulated DEGs are broadly expressed across neuronal clusters. Enrichment analyses of these DEGs demonstrated that downregulated genes are uniquely expressed in the LS, and associated with either synaptic plasticity or neurodevelopmental disorders. Upregulated genes are enriched in LS microglia, associated with immune response and inflammation, and linked to both neurodegenerative disease and neuropsychiatric disorders. In addition, many of these genes are implicated in regulating social behaviors. In summary, the findings implicate TrkB signaling in the LS as a critical regulator of gene networks associated with psychiatric disorders that display social deficits, including schizophrenia and autism, and with neurodegenerative diseases, including Alzheimer's.

Proteome-wide mendelian randomization identifies causal plasma proteins in venous thromboembolism development

Genome-wide association studies (GWAS) have identified numerous risk loci for venous thromboembolism (VTE), but it is challenging to decipher the underlying mechanisms. We employed an integrative analytical pipeline to transform genetic associations to identify novel plasma proteins for VTE. Proteome-wide association studies (PWAS) were determined by functional summary-based imputation leveraging data from a genome-wide association analysis (14,429 VTE patients, 267,037 controls), blood proteomes (1348 cases), followed by Mendelian randomization, Bayesian colocalization, protein-protein interaction, and pathway enrichment analysis. Twenty genetically regulated circulating protein abundances (F2, F11, ABO, PLCG2, LRP4, PLEK, KLKB1, PROC, KNG1, THBS2, SERPINA1, RARRES2, CEL, GP6, SERPINE2, SERPINA10, OBP2B, EFEMP1, F5, and MSR1) were associated with VTE. Of these 13 proteins demonstrated Mendelian randomized correlations. Six proteins (F2, F11, PLEK, SERPINA1, RARRES2, and SERPINE2) had strong support in colocalization analysis. Utilizing multidimensional data, this study suggests PLEK, SERPINA1, and SERPINE2 as compelling proteins that may provide key hints for future research and possible diagnostic and therapeutic targets for VTE.

Peripheral signature of altered synaptic integrity in young onset cannabis use disorder: A proteomic study of circulating extracellular vesicles

BACKGROUND

Rates of Cannabis Use Disorder (CUD) are highest amongst young adults. Paucity of brain tissue samples limits the ability to examine the molecular basis of cannabis related neuropathology. Proteomic studies of neuron-derived extracellular vesicles (NDEs) isolated from the biofluids may reveal markers of neuropathology in CUD.

METHODS

NDEs were extracted using ExoSORT, an immunoaffinity method to enrich NDEs from plasma samples from patients with young onset CUD and matched controls. Differential proteomic profiles were explored with Label Free Quantification (LFQ) mass spectrometry. Selected proteins were validated using orthogonal methods.

RESULTS

A total of 231 (±10) proteins were identified in NDE preparations from CUD and controls of which 28 were differentially abundant between groups. The difference in abundance of properdin (CFP gene) was statistically significant. SHANK1 (SHANK1 gene), an adapter protein at the post-synaptic density, was nominally depleted in the CUD NDE preparations.

CONCLUSION

In this pilot study, we noted a decrease in SHANK1 protein, involved in the structural and functional integrity of glutamatergic post-synapse, a potential peripheral signature of CUD neuropathology. The study shows that LFQ mass spectrometry proteomic analysis of NDEs derived from plasma may yield important insights into the synaptic pathology associated with CUD.

Identifying Novel Drug Targets for Epilepsy Through a Brain Transcriptome-Wide Association Study and Protein-Wide Association Study with Chemical-Gene-Interaction Analysis

Epilepsy is a severe neurological condition affecting 50-65 million individuals worldwide that can lead to brain damage. Nevertheless, the etiology of epilepsy remains poorly understood. Meta-analyses of genome-wide association studies involving 15,212 epilepsy cases and 29,677 controls of the ILAE Consortium cohort were used to conduct transcriptome-wide association studies (TWAS) and protein-wide association studies (PWAS). Furthermore, a protein-protein interaction (PPI) network was generated using the STRING database, and significant epilepsy-susceptible genes were verified using chip data. Chemical-related gene set enrichment analysis (CGSEA) was performed to determine novel drug targets for epilepsy. TWAS analysis identified 21,170 genes, of which 58 were significant (TWASfdr < 0.05) in ten brain regions, and 16 differentially expressed genes were verified based on mRNA expression profiles. The PWAS identified 2249 genes, of which 2 were significant (PWASfdr < 0.05). Through chemical-gene set enrichment analysis, 287 environmental chemicals associated with epilepsy were identified. We identified five significant genes (WIPF1, IQSEC1, JAM2, ICAM3, and ZNF143) that had causal relationships with epilepsy. CGSEA identified 159 chemicals that were significantly correlated with epilepsy (Pcgsea < 0.05), such as pentobarbital, ketone bodies, and polychlorinated biphenyl. In summary, we performed TWAS, PWAS (for genetic factors), and CGSEA (for environmental factors) analyses and identified several epilepsy-associated genes and chemicals. The results of this study will contribute to our understanding of genetic and environmental factors for epilepsy and may predict novel drug targets.

TrkB-dependent regulation of molecular signaling across septal cell types

The lateral septum (LS), a GABAergic structure located in the basal forebrain, is implicated in social behavior, learning and memory. We previously demonstrated that expression of tropomyosin kinase receptor B (TrkB) in LS neurons is required for social novelty recognition. To better understand molecular mechanisms by which TrkB signaling controls behavior, we locally knocked down TrkB in LS and used bulk RNA-sequencing to identify changes in gene expression downstream of TrkB. TrkB knockdown induces upregulation of genes associated with inflammation and immune responses, and downregulation of genes associated with synaptic signaling and plasticity. Next, we generated one of the first atlases of molecular profiles for LS cell types using single nucleus RNA-sequencing (snRNA-seq). We identified markers for the septum broadly, and the LS specifically, as well as for all neuronal cell types. We then investigated whether the differentially expressed genes (DEGs) induced by TrkB knockdown map to specific LS cell types. Enrichment testing identified that downregulated DEGs are broadly expressed across neuronal clusters. Enrichment analyses of these DEGs demonstrated that downregulated genes are uniquely expressed in the LS, and associated with either synaptic plasticity or neurodevelopmental disorders. Upregulated genes are enriched in LS microglia, associated with immune response and inflammation, and linked to both neurodegenerative disease and neuropsychiatric disorders. In addition, many of these genes are implicated in regulating social behaviors. In summary, the findings implicate TrkB signaling in the LS as a critical regulator of gene networks associated with psychiatric disorders that display social deficits, including schizophrenia and autism, and with neurodegenerative diseases, including Alzheimer's.

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.