The impact of chronic pain on brain gene expression.

Background: Chronic pain affects one fifth of American adults, contributing significant public health burden. Chronic pain mechanisms can be further understood through investigating brain gene expression. Methods: We tested differentially expressed genes (DEGs) in chronic pain, migraine, lifetime fentanyl and oxymorphone use, and with chronic pain genetic risk in four brain regions (dACC, DLPFC, MeA, BLA) and imputed cell type expression data from 304 postmortem donors. We compared findings across traits and with independent transcriptomics resources, and performed gene-set enrichment. Results: We identified two chronic pain DEGs: B4GALT and VEGFB in bulk dACC. We found over 2000 (primarily BLA microglia) chronic pain cell type DEGs. Findings were enriched for mouse microglia pain genes, and for hypoxia and immune response. Cross-trait DEG overlap was minimal. Conclusions: Chronic pain-associated gene expression is heterogeneous across cell type, largely distinct from that in pain-related traits, and shows BLA microglia are a key cell type.

Tuning parameters for polygenic risk score methods using GWAS summary statistics from training data.

Various polygenic risk scores (PRS) methods have been proposed to combine the estimated effects of single nucleotide polymorphisms (SNPs) to predict genetic risks for common diseases, using data collected from genome-wide association studies (GWAS). Some methods require external individual-level GWAS dataset for parameter tuning, posing privacy and security-related concerns. Leaving out partial data for parameter tuning can also reduce model prediction accuracy. In this article, we propose PRStuning, a method that tunes parameters for different PRS methods using GWAS summary statistics from the training data. PRStuning predicts the PRS performance with different parameters, and then selects the best-performing parameters. Because directly using training data effects tends to overestimate the performance in the testing data, we adopt an empirical Bayes approach to shrinking the predicted performance in accordance with the genetic architecture of the disease. Extensive simulations and real data applications demonstrate PRStuning's accuracy across PRS methods and parameters.

Uncommon Protein-Coding Variants Associated With Suicide Attempt in a Diverse Sample of U.S. Army Soldiers.

BACKGROUND: Suicide is a societal and public health concern of global scale. Identifying genetic risk factors for suicide attempt can characterize underlying biology and enable early interventions to prevent deaths. Recent studies have described common genetic variants for suicide-related behaviors. Here, we advance this search for genetic risk by analyzing the association between suicide attempt and uncommon variation exome-wide in a large, ancestrally diverse sample. METHODS: We sequenced whole genomes of 13,584 soldiers from the Army STARRS (Army Study to Assess Risk and Resilience in Servicemembers), including 979 individuals with a history of suicide attempt. Uncommon, nonsilent protein-coding variants were analyzed exome-wide for association with suicide attempt using gene-collapsed and single-variant analyses. RESULTS: We identified 19 genes with variants enriched in individuals with history of suicide attempt, either through gene-collapsed or single-variant analysis (Bonferroni padjusted < .05). These genes were CIB2, MLF1, HERC1, YWHAE, RCN2, VWA5B1, ATAD3A, NACA, EP400, ZNF585A, LYST, RC3H2, PSD3, STARD9, SGMS1, ACTR6, RGS7BP, DIRAS2, and KRTAP10-1. Most genes had variants across multiple genomic ancestry groups. Seventeen of these genes were expressed in healthy brain tissue, with 9 genes expressed at the highest levels in the brain versus other tissues. Brains from individuals deceased from suicide aberrantly expressed RGS7BP (padjusted = .035) in addition to nominally significant genes including YWHAE and ACTR6, all of which have reported associations with other mental disorders. CONCLUSIONS: These results advance the molecular characterization of suicide attempt behavior and support the utility of whole-genome sequencing for complementing the findings of genome-wide association studies in suicide research.

Novel hippocampal genes involved in enhanced susceptibility to chronic pain-induced behavioral emotionality.

Altered mood and psychiatric disorders are commonly associated with chronic pain conditions; however, brain mechanisms linking pain and comorbid clinical depression are still largely unknown. In this study, we aimed to identify whether key genes/cellular mechanisms underlie susceptibility/resiliency to development of depressive-like behaviors during chronic pain state. Genome-wide RNA-seq analysis was used to examine the transcriptomic profile of the hippocampus, a limbic brain region that regulates mood and stress responses, from male rats exposed to chronic inflammatory pain. Pain-exposed animals were separated into either 'resilient' or 'susceptible' to development of enhanced behavioral emotionality based on behavioral testing. RNA-seq bioinformatic analysis, followed by validation using qPCR, revealed dysregulation of hippocampal genes involved in neuroinflammation, cell cycle/neurogenesis and blood-brain barrier integrity. Specifically, ADAM Metallopeptidase Domain 8 (Adam8) and Aurora Kinase B (Aurkb), genes with functional roles in activation of the NLRP3 inflammasome and microgliosis, respectively, were significantly upregulated in the hippocampus of 'susceptible' animals expressing increased behavioral emotionality. In addition, genes associated with blood-brain barrier integrity, such as the Claudin 4 (Cldn4), a tight junction protein and a known marker of astrocyte activation, were also significantly dysregulated between 'resilient' or 'susceptible' pain groups. Furthermore, differentially expressed genes (DEGs) were further characterized in rodents stress models to determine whether their hippocampal dysregulation is driven by common stress responses vs. affective pain processing. Altogether these results continue to strengthen the connection between dysregulation of hippocampal genes involved in neuroinflammatory and neurodegenerative processes with increased behavioral emotionality often expressed in chronic pain state.

Ketamine and the neurobiology of depression: Toward next-generation rapid-acting antidepressant treatments.

Ketamine has emerged as a transformative and mechanistically novel pharmacotherapy for depression. Its rapid onset of action, efficacy for treatment-resistant symptoms, and protection against relapse distinguish it from prior antidepressants. Its discovery emerged from a reconceptualization of the neurobiology of depression and, in turn, insights from the elaboration of its mechanisms of action inform studies of the pathophysiology of depression and related disorders. It has been 25 y since we first presented our ketamine findings in depression. Thus, it is timely for this review to consider what we have learned from studies of ketamine and to suggest future directions for the optimization of rapid-acting antidepressant treatment.

iHerd: an integrative hierarchical graph representation learning framework to quantify network changes and prioritize risk genes in disease.

Different genes form complex networks within cells to carry out critical cellular functions, while network alterations in this process can potentially introduce downstream transcriptome perturbations and phenotypic variations. Therefore, developing efficient and interpretable methods to quantify network changes and pinpoint driver genes across conditions is crucial. We propose a hierarchical graph representation learning method, called iHerd. Given a set of networks, iHerd first hierarchically generates a series of coarsened sub-graphs in a data-driven manner, representing network modules at different resolutions (e.g., the level of signaling pathways). Then, it sequentially learns low-dimensional node representations at all hierarchical levels via efficient graph embedding. Lastly, iHerd projects separate gene embeddings onto the same latent space in its graph alignment module to calculate a rewiring index for driver gene prioritization. To demonstrate its effectiveness, we applied iHerd on a tumor-to-normal GRN rewiring analysis and cell-type-specific GCN analysis using single-cell multiome data of the brain. We showed that iHerd can effectively pinpoint novel and well-known risk genes in different diseases. Distinct from existing models, iHerd's graph coarsening for hierarchical learning allows us to successfully classify network driver genes into early and late divergent genes (EDGs and LDGs), emphasizing genes with extensive network changes across and within signaling pathway levels. This unique approach for driver gene classification can provide us with deeper molecular insights. The code is freely available at https://github.com/aicb-ZhangLabs/iHerd. All other relevant data are within the manuscript and supporting information files.

Tuning Parameters for Polygenic Risk Score Methods Using GWAS Summary Statistics from Training Data.

Predicting genetic risks for common diseases may improve their prevention and early treatment. In recent years, various additive-model-based polygenic risk scores (PRS) methods have been proposed to combine the estimated effects of single nucleotide polymorphisms (SNPs) using data collected from genome-wide association studies (GWAS). Some of these methods require access to another external individual-level GWAS dataset to tune the hyperparameters, which can be difficult because of privacy and security-related concerns. Additionally, leaving out partial data for hyperparameter tuning can reduce the predictive accuracy of the constructed PRS model. In this article, we propose a novel method, called PRStuning, to automatically tune hyperparameters for different PRS methods using only GWAS summary statistics from the training data. The core idea is to first predict the performance of the PRS method with different parameter values, and then select the parameters with the best prediction performance. Because directly using the effects observed from the training data tends to overestimate the performance in the testing data (a phenomenon known as overfitting), we adopt an empirical Bayes approach to shrinking the predicted performance in accordance with the estimated genetic architecture of the disease. Results from extensive simulations and real data applications demonstrate that PRStuning can accurately predict the PRS performance across PRS methods and parameters, and it can help select the best-performing parameters.

Single-Nucleus Transcriptome Profiling of Dorsolateral Prefrontal Cortex: Mechanistic Roles for Neuronal Gene Expression, Including the 17q21.31 Locus, in PTSD Stress Response.

OBJECTIVE: Multidisciplinary studies of posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) implicate the dorsolateral prefrontal cortex (DLPFC) in disease risk and pathophysiology. Postmortem brain studies have relied on bulk-tissue RNA sequencing (RNA-seq), but single-cell RNA-seq is needed to dissect cell-type-specific mechanisms. The authors conducted the first single-nucleus RNA-seq postmortem brain study in PTSD to elucidate disease transcriptomic pathology with cell-type-specific resolution. METHOD: Profiling of 32 DLPFC samples from 11 individuals with PTSD, 10 with MDD, and 11 control subjects was conducted (∼415K nuclei; >13K cells per sample). A replication sample included 15 DLPFC samples (∼160K nuclei; >11K cells per sample). RESULTS: Differential gene expression analyses identified significant single-nucleus RNA-seq differentially expressed genes (snDEGs) in excitatory (EX) and inhibitory (IN) neurons and astrocytes, but not in other cell types or bulk tissue. MDD samples had more false discovery rate-corrected significant snDEGs, and PTSD samples had a greater replication rate. In EX and IN neurons, biological pathways that were differentially enriched in PTSD compared with MDD included glucocorticoid signaling. Furthermore, glucocorticoid signaling in induced pluripotent stem cell (iPSC)-derived cortical neurons demonstrated greater relevance in PTSD and opposite direction of regulation compared with MDD, especially in EX neurons. Many snDEGs were from the 17q21.31 locus and are particularly interesting given causal roles in disease pathogenesis and DLPFC-based neuroimaging (PTSD: ARL17B, LINC02210-CRHR1, and LRRC37A2; MDD: LRRC37A and LRP4), while others were regulated by glucocorticoids in iPSC-derived neurons (PTSD: SLC16A6, TAF1C; MDD: CDH3). CONCLUSIONS: The study findings point to cell-type-specific mechanisms of brain stress response in PTSD and MDD, highlighting the importance of examining cell-type-specific gene expression and indicating promising novel biomarkers and therapeutic targets.

Profiling neuronal methylome and hydroxymethylome of opioid use disorder in the human orbitofrontal cortex.

Opioid use disorder (OUD) is influenced by genetic and environmental factors. While recent research suggests epigenetic disturbances in OUD, this is mostly limited to DNA methylation (5mC). DNA hydroxymethylation (5hmC) has been widely understudied. We conducted a multi-omics profiling of OUD in a male cohort, integrating neuronal-specific 5mC and 5hmC as well as gene expression profiles from human postmortem orbitofrontal cortex (OUD = 12; non-OUD = 26). Single locus methylomic analysis and co-methylation analysis showed a higher number of OUD-associated genes and gene networks for 5hmC compared to 5mC; these were enriched for GPCR, Wnt, neurogenesis, and opioid signaling. 5hmC marks also showed a higher correlation with gene expression patterns and enriched for GWAS of psychiatric traits. Drug interaction analysis revealed interactions with opioid-related drugs, some used as OUD treatments. Our multi-omics findings suggest an important role of 5hmC and reveal loci epigenetically dysregulated in OFC neurons of individuals with OUD.

Decoding the role of transcriptomic clocks in the human prefrontal cortex.

Aging is a complex process with interindividual variability, which can be measured by aging biological clocks. Aging clocks are machine-learning algorithms guided by biological information and associated with mortality risk and a wide range of health outcomes. One of these aging clocks are transcriptomic clocks, which uses gene expression data to predict biological age; however, their functional role is unknown. Here, we profiled two transcriptomic clocks (RNAAgeCalc and knowledge-based deep neural network clock) in a large dataset of human postmortem prefrontal cortex (PFC) samples. We identified that deep-learning transcriptomic clock outperforms RNAAgeCalc to predict transcriptomic age in the human PFC. We identified associations of transcriptomic clocks with psychiatric-related traits. Further, we applied system biology algorithms to identify common gene networks among both clocks and performed pathways enrichment analyses to assess its functionality and prioritize genes involved in the aging processes. Identified gene networks showed enrichment for diseases of signal transduction by growth factor receptors and second messenger pathways. We also observed enrichment of genome-wide signals of mental and physical health outcomes and identified genes previously associated with human brain aging. Our findings suggest a link between transcriptomic aging and health disorders, including psychiatric traits. Further, it reveals functional genes within the human PFC that may play an important role in aging and health risk.

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.

Genetic liability to suicidal thoughts and behaviors and risk of suicide attempt in US military veterans: moderating effects of cumulative trauma burden.

BACKGROUND: Little is known about environmental factors that may influence associations between genetic liability to suicidality and suicidal behavior. METHODS: This study examined whether a suicidality polygenic risk score (PRS) derived from a large genome-wide association study (N = 122,935) was associated with suicide attempts in a population-based sample of European-American US military veterans (N = 1664; 92.5% male), and whether cumulative lifetime trauma exposure moderated this association. RESULTS: Eighty-five veterans (weighted 6.3%) reported a history of suicide attempt. After adjusting for sociodemographic and psychiatric characteristics, suicidality PRS was associated with lifetime suicide attempt (odds ratio 2.65; 95% CI 1.37-5.11). A significant suicidality PRS-by-trauma exposure interaction emerged, such that veterans with higher levels of suicidality PRS and greater trauma burden had the highest probability of lifetime suicide attempt (16.6%), whereas the probability of attempts was substantially lower among those with high suicidality PRS and low trauma exposure (1.4%). The PRS-by-trauma interaction effect was enriched for genes implicated in cellular and developmental processes, and nervous system development, with variants annotated to the DAB2 and SPNS2 genes, which are implicated in inflammatory processes. Drug repurposing analyses revealed upregulation of suicide gene-sets in the context of medrysone, a drug targeting chronic inflammation, and clofibrate, a triacylglyceride level lowering agent. CONCLUSION: Results suggest that genetic liability to suicidality is associated with increased risk of suicide attempt among veterans, particularly in the presence of high levels of cumulative trauma exposure. Additional research is warranted to investigate whether incorporation of genomic information may improve suicide prediction models.

Common genetic variation impacts stress response in the brain.

To explain why individuals exposed to identical stressors experience divergent clinical outcomes, we determine how molecular encoding of stress modifies genetic risk for brain disorders. Analysis of post-mortem brain (n=304) revealed 8557 stress-interactive expression quantitative trait loci (eQTLs) that dysregulate expression of 915 eGenes in response to stress, and lie in stress-related transcription factor binding sites. Response to stress is robust across experimental paradigms: up to 50% of stress-interactive eGenes validate in glucocorticoid treated hiPSC-derived neurons (n=39 donors). Stress-interactive eGenes show brain region- and cell type-specificity, and, in post-mortem brain, implicate glial and endothelial mechanisms. Stress dysregulates long-term expression of disorder risk genes in a genotype-dependent manner; stress-interactive transcriptomic imputation uncovered 139 novel genes conferring brain disorder risk only in the context of traumatic stress. Molecular stress-encoding explains individualized responses to traumatic stress; incorporating trauma into genomic studies of brain disorders is likely to improve diagnosis, prognosis, and drug discovery.

The long noncoding RNA FEDORA is a cell type- and sex-specific regulator of depression.

Women suffer from depression at twice the rate of men, but the underlying molecular mechanisms are poorly understood. Here, we identify marked baseline sex differences in the expression of long noncoding RNAs (lncRNAs), a class of regulatory transcripts, in human postmortem brain tissue that are profoundly lost in depression. One such human lncRNA, RP11-298D21.1 (which we termed FEDORA), is enriched in oligodendrocytes and neurons and up-regulated in the prefrontal cortex (PFC) of depressed females only. We found that virally expressing FEDORA selectively either in neurons or in oligodendrocytes of PFC promoted depression-like behavioral abnormalities in female mice only, changes associated with cell type-specific regulation of synaptic properties, myelin thickness, and gene expression. We also found that blood FEDORA levels have diagnostic implications for depressed women and are associated with clinical response to ketamine. These findings demonstrate the important role played by lncRNAs, and FEDORA in particular, in shaping the sex-specific landscape of the brain and contributing to sex differences in depression.

Venus: An efficient virus infection detection and fusion site discovery method using single-cell and bulk RNA-seq data.

Early and accurate detection of viruses in clinical and environmental samples is essential for effective public healthcare, treatment, and therapeutics. While PCR detects potential pathogens with high sensitivity, it is difficult to scale and requires knowledge of the exact sequence of the pathogen. With the advent of next-gen single-cell sequencing, it is now possible to scrutinize viral transcriptomics at the finest possible resolution-cells. This newfound ability to investigate individual cells opens new avenues to understand viral pathophysiology with unprecedented resolution. To leverage this ability, we propose an efficient and accurate computational pipeline, named Venus, for virus detection and integration site discovery in both single-cell and bulk-tissue RNA-seq data. Specifically, Venus addresses two main questions: whether a tissue/cell type is infected by viruses or a virus of interest? And if infected, whether and where has the virus inserted itself into the human genome? Our analysis can be broken into two parts-validation and discovery. Firstly, for validation, we applied Venus on well-studied viral datasets, such as HBV- hepatocellular carcinoma and HIV-infection treated with antiretroviral therapy. Secondly, for discovery, we analyzed datasets such as HIV-infected neurological patients and deeply sequenced T-cells. We detected viral transcripts in the novel target of the brain and high-confidence integration sites in immune cells. In conclusion, here we describe Venus, a publicly available software which we believe will be a valuable virus investigation tool for the scientific community at large.

Molecular and cellular evolution of the primate dorsolateral prefrontal cortex.

The granular dorsolateral prefrontal cortex (dlPFC) is an evolutionary specialization of primates that is centrally involved in cognition. We assessed more than 600,000 single-nucleus transcriptomes from adult human, chimpanzee, macaque, and marmoset dlPFC. Although most cell subtypes defined transcriptomically are conserved, we detected several that exist only in a subset of species as well as substantial species-specific molecular differences across homologous neuronal, glial, and non-neural subtypes. The latter are exemplified by human-specific switching between expression of the neuropeptide somatostatin and tyrosine hydroxylase, the rate-limiting enzyme in dopamine production in certain interneurons. The above molecular differences are also illustrated by expression of the neuropsychiatric risk gene FOXP2, which is human-specific in microglia and primate-specific in layer 4 granular neurons. We generated a comprehensive survey of the dlPFC cellular repertoire and its shared and divergent features in anthropoid primates.

Translator: A Transfer Learning Approach to Facilitate Single-Cell ATAC-Seq Data Analysis from Reference Dataset.

Recent advances in single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) have allowed simultaneous epigenetic profiling over thousands of individual cells to dissect the cellular heterogeneity and elucidate regulatory mechanisms at the finest possible resolution. However, scATAC-seq is challenging to model computationally due to the ultra-high dimensionality, low signal-to-noise ratio, complex feature interactions, and high vulnerability to various confounding factors. In this study, we present Translator, an efficient transfer learning approach to capture generalizable chromatin interactions from high-quality (HQ) reference scATAC-seq data to obtain robust cell representations in low-to-moderate quality target scATAC-seq data. We applied Translator on various simulated and real scATAC-seq datasets and demonstrated that Translator could learn more biologically meaningful cell representations than other methods by incorporating information learned from the reference data, thus facilitating various downstream analyses such as clustering and motif enrichment measurements. Moreover, Translator's block-wise deep learning framework can handle nonlinear relationships with restricted connections using fewer parameters to boost computational efficiency through Graphics Processing Unit (GPU) parallelism. Finally, we have implemented Translator as a free software package available for the community to leverage large-scale, HQ reference data to study target scATAC-seq data.

InsuLock: A Weakly Supervised Learning Approach for Accurate Insulator Prediction, and Variant Impact Quantification.

Mapping chromatin insulator loops is crucial to investigating genome evolution, elucidating critical biological functions, and ultimately quantifying variant impact in diseases. However, chromatin conformation profiling assays are usually expensive, time-consuming, and may report fuzzy insulator annotations with low resolution. Therefore, we propose a weakly supervised deep learning method, InsuLock, to address these challenges. Specifically, InsuLock first utilizes a Siamese neural network to predict the existence of insulators within a given region (up to 2000 bp). Then, it uses an object detection module for precise insulator boundary localization via gradient-weighted class activation mapping (~40 bp resolution). Finally, it quantifies variant impacts by comparing the insulator score differences between the wild-type and mutant alleles. We applied InsuLock on various bulk and single-cell datasets for performance testing and benchmarking. We showed that it outperformed existing methods with an AUROC of ~0.96 and condensed insulator annotations to ~2.5% of their original size while still demonstrating higher conservation scores and better motif enrichments. Finally, we utilized InsuLock to make cell-type-specific variant impacts from brain scATAC-seq data and identified a schizophrenia GWAS variant disrupting an insulator loop proximal to a known risk gene, indicating a possible new mechanism of action for the disease.

Positive modulation of N-methyl-D-aspartate receptors in the mPFC reduces the spontaneous recovery of fear.

N-methyl-D-aspartate receptor (NMDAR) modulators have recently received increased attention as potential therapeutics for posttraumatic stress disorder (PTSD). Here, we tested a novel NMDAR-positive modulator, NYX-783, in the following two rodent models of PTSD: an auditory fear-conditioning model and a single-prolonged stress (SPS) model. We examined the ability of NYX-783 to reduce subsequent fear-based behaviors by measuring enhanced fear extinction and reduced spontaneous recovery (spontaneous return of fear) in male mice. NYX-783 administration significantly reduced spontaneous recovery in both PTSD models and enhanced fear extinction in the SPS model. Furthermore, NYX-783 increased the NMDA-induced inward currents of excitatory and inhibitory neurons in the infralimbic medial prefrontal cortex (IL mPFC) and that the GluN2B subunit of NMDARs on pyramidal neurons in the IL mPFC is required for its effect on spontaneous recovery. The downstream expression of brain-derived neurotrophic factor was required for NYX-783 to achieve its behavioral effect. These results elucidate the cellular targets of NYX-783 and the molecular mechanisms underlying the inhibition of spontaneous recovery. These preclinical findings support the hypothesis that NYX-783 may have therapeutic potential for PTSD treatment and may be particularly useful for inhibiting spontaneous recovery.