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.

Gene expression in the dorsolateral and ventromedial prefrontal cortices implicates immune-related gene networks in PTSD

Studies evaluating neuroimaging, genetically predicted gene expression, and pre-clinical genetic models of PTSD, have identified PTSD-related abnormalities in the prefrontal cortex (PFC) of the brain, particularly in dorsolateral and ventromedial PFC (dlPFC and vmPFC). In this study, RNA sequencing was used to examine gene expression in the dlPFC and vmPFC using tissue from the VA National PTSD Brain Bank in donors with histories of PTSD with or without depression (dlPFC n = 38, vmPFC n = 35), depression cases without PTSD (n = 32), and psychopathology-free controls (dlPFC n = 24, vmPFC n = 20). Analyses compared PTSD cases to controls. Follow-up analyses contrasted depression cases to controls. Twenty-one genes were differentially expressed in PTSD after strict multiple testing correction. PTSD-associated genes with roles in learning and memory (FOS, NR4A1), immune regulation (CFH, KPNA1) and myelination (MBP, MOBP, ERMN) were identified. PTSD-associated genes partially overlapped depression-associated genes. Co-expression network analyses identified PTSD-associated networks enriched for immune-related genes across the two brain regions. However, the immune-related genes and association patterns were distinct. The immune gene IL1B was significantly associated with PTSD in candidate-gene analysis and was an upstream regulator of PTSD-associated genes in both regions. There was evidence of replication of dlPFC associations in an independent cohort from a recent study, and a strong correlation between the dlPFC PTSD effect sizes for significant genes in the two studies (r = 0.66, p < 2.2 × 10⁻¹⁶). In conclusion, this study identified several novel PTSD-associated genes and brain region specific PTSD-associated immune-related networks.

Klotho, PTSD, and advanced epigenetic age in cortical tissue

This study examined the klotho (KL) longevity gene polymorphism rs9315202 and psychopathology, including posttraumatic stress disorder (PTSD), depression, and alcohol-use disorders, in association with advanced epigenetic age in three postmortem cortical tissue regions: dorsolateral and ventromedial prefrontal cortices and motor cortex. Using data from the VA National PTSD Brain Bank (n = 117), we found that rs9315202 interacted with PTSD to predict advanced epigenetic age in motor cortex among the subset of relatively older (>=45 years), white non-Hispanic decedents (corrected p = 0.014, n = 42). An evaluation of 211 additional common KL variants revealed that only variants in linkage disequilibrium with rs9315202 showed similarly high levels of significance. Alcohol abuse was nominally associated with advanced epigenetic age in motor cortex (p = 0.039, n = 114). The rs9315202 SNP interacted with PTSD to predict decreased KL expression via DNAm age residuals in motor cortex among older white non-Hispanics decedents (indirect β = -0.198, p = 0.027). Finally, in dual-luciferase enhancer reporter system experiments, we found that inserting the minor allele of rs9315202 in a human kidney cell line HK-2 genomic DNA resulted in a change in KL transcriptional activities, likely operating via long noncoding RNA in this region. This was the first study to examine multiple forms of psychopathology in association with advanced DNA methylation age across several brain regions, to extend work concerning the association between rs9315202 and advanced epigenetic to brain tissue, and to identify the effects of rs9315202 on KL gene expression. KL augmentation holds promise as a therapeutic intervention to slow the pace of cellular aging, disease onset, and neuropathology, particularly in older, stressed populations.

PTSD and the klotho longevity gene: Evaluation of longitudinal effects on inflammation via DNA methylation

BACKGROUND

Longevity gene klotho (KL) is associated with age-related phenotypes including lifespan, cardiometabolic disorders, cognition, and brain morphology, in part, by conferring protection against inflammation. We hypothesized that the KL/inflammation association might be altered in the presence of psychiatric stress and operate via epigenetic pathways. We examined KL polymorphisms, and their interaction with posttraumatic stress disorder (PTSD) symptoms, in association with KL DNA methylation in blood. We further examined KL DNA methylation as a predictor of longitudinal changes in a peripheral biomarker of inflammation (C-reactive protein; CRP).

METHODS

The sample comprised 309 white non-Hispanic military veterans (93.5 % male; mean age: 32 years, range: 19-65; 30 % PTSD per structured diagnostic interview); 111 were reassessed approximately two years later.

RESULTS

Analyses revealed a methylation quantitative trait locus at rs9527025 (C370S, previously implicated in numerous studies of aging) in association with a Cytosine-phosphate-Guanine site (cg00129557; B = -.65, p =  1.29 X 10^-20), located within a DNase hypersensitivity site in the body of KL. There was also a rs9527025 x PTSD severity interaction (B = .004, p = .035) on methylation at this locus such that the minor allele was associated with reduced cg00129557 methylation in individuals with few or no PTSD symptoms while this effect was attenuated in those with elevated levels of PTSD. Path models revealed that methylation at cg00129557 was inversely associated with CRP over time (B = -.14, p = .005), controlling for baseline CRP. There was also an indirect effect of rs9527025 X PTSD on subsequent CRP via cg00129557 methylation (indirect B = -.002, p = .033).

CONCLUSIONS

Results contribute to our understanding of the epigenetic correlates of inflammation in PTSD and suggest that KL methylation may be a mechanism by which KL genotype confers risk vs. resilience to accelerated aging in those experiencing traumatic stress.

Analysis of Genetically Regulated Gene Expression Identifies a Prefrontal PTSD Gene, SNRNP35, Specific to Military Cohorts

To reveal post-traumatic stress disorder (PTSD) genetic risk influences on tissue-specific gene expression, we use brain and non-brain transcriptomic imputation. We impute genetically regulated gene expression (GReX) in 29,539 PTSD cases and 166,145 controls from 70 ancestry-specific cohorts and identify 18 significant GReX-PTSD associations corresponding to specific tissue-gene pairs. The results suggest substantial genetic heterogeneity based on ancestry, cohort type (military versus civilian), and sex. Two study-wide significant PTSD associations are identified in European and military European cohorts; ZNF140 is predicted to be upregulated in whole blood, and SNRNP35 is predicted to be downregulated in dorsolateral prefrontal cortex, respectively. In peripheral leukocytes from 175 marines, the observed PTSD differential gene expression correlates with the predicted differences for these individuals, and deployment stress produces glucocorticoid-regulated expression changes that include downregulation of both ZNF140 and SNRNP35. SNRNP35 knockdown in cells validates its functional role in U12-intron splicing. Finally, exogenous glucocorticoids in mice downregulate prefrontal Snrnp35 expression.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

The PPM1F gene moderates the association between PTSD and cortical thickness

BACKGROUND

Evidence suggests that single nucleotide polymorphisms (SNPs) in genes involved in serotonergic signaling and stress response pathways moderate associations between PTSD and cortical thickness. This study examined a genetic regulator of these pathways, the PPM1F gene, which has also been implicated in mechanisms of stress responding and is differentially expressed in individuals with comorbid PTSD and depression compared to controls.

METHODS

Drawing from a sample of 240 white non-Hispanic trauma-exposed veterans, we tested 18 SNPs spanning the PPM1F gene for association with PTSD and cortical thickness.

RESULTS

Analyses revealed six PPM1F SNPs that moderated associations between PTSD symptom severity and cortical thickness of bilateral superior frontal and orbitofrontal regions as well as the right pars triangularis (all corrected p's < 0.05) such that greater PTSD severity was related to reduced cortical thickness as a function of genotype. A whole-cortex vertex-wise analysis using the most associated SNP (rs9610608) revealed this effect to be localized to a cluster in the right superior frontal gyrus (cluster-corrected p < 0.02).

LIMITATIONS

Limitations of this study include the small sample size and that the sample was all-white, non-Hispanic predominately male veterans.

CONCLUSIONS

These results extend prior work linking PPM1F to PTSD and suggest that variants in this gene may have bearing on the neural integrity of the prefrontal cortex (PFC).

Investigation of bidirectional longitudinal associations between advanced epigenetic age and peripheral biomarkers of inflammation and metabolic syndrome

Epigenetic age estimations based on DNA methylation (DNAm) can predict human chronological age with a high level of accuracy. These DNAm age algorithms can also be used to index advanced cellular age, when estimated DNAm age exceeds chronological age. Advanced DNAm age has been associated with several diseases and metabolic and inflammatory pathology, but the causal direction of this association is unclear. The goal of this study was to examine potential bidirectional associations between advanced epigenetic age and metabolic and inflammatory markers over time in a longitudinal cohort of 179 veterans with a high prevalence of posttraumatic stress disorder (PTSD) who were assessed over the course of two years. Analyses focused on two commonly investigated metrics of advanced DNAm age derived from the Horvath (developed across multiple tissue types) and Hannum (developed in whole blood) DNAm age algorithms. Results of cross-lagged panel models revealed that advanced Hannum DNAm age at Time 1 (T1) was associated with increased (i.e., accounting for T1 levels) metabolic syndrome (MetS) severity at Time 2 (T2; p = < 0.001). This association was specific to worsening lipid panels and indicators of abdominal obesity (p = 0.001). In contrast, no baseline measures of inflammation or metabolic pathology were associated with changes in advanced epigenetic age over time. No associations emerged between advanced Horvath DNAm age and any of the examined biological parameters. Results suggest that advanced epigenetic age, when measured using an algorithm developed in whole blood, may be a prognostic marker of pathological metabolic processes. This carries implications for understanding pathways linking advanced epigenetic age to morbidity and mortality.

Posttraumatic psychopathology and the pace of the epigenetic clock: a longitudinal investigation

BACKGROUND

Posttraumatic stress disorder (PTSD) and stress/trauma exposure are cross-sectionally associated with advanced DNA methylation age relative to chronological age. However, longitudinal inquiry and examination of associations between advanced DNA methylation age and a broader range of psychiatric disorders is lacking. The aim of this study was to examine if PTSD, depression, generalized anxiety, and alcohol-use disorders predicted acceleration of DNA methylation age over time (i.e. an increasing pace, or rate of advancement, of the epigenetic clock).

METHODS

Genome-wide DNA methylation and a comprehensive set of psychiatric symptoms and diagnoses were assessed in 179 Iraq/Afghanistan war veterans who completed two assessments over the course of approximately 2 years. Two DNA methylation age indices (Horvath and Hannum), each a weighted index of an array of genome-wide DNA methylation probes, were quantified. The pace of the epigenetic clock was operationalized as change in DNA methylation age as a function of time between assessments.

RESULTS

Analyses revealed that alcohol-use disorders (p = 0.001) and PTSD avoidance and numbing symptoms (p = 0.02) at Time 1 were associated with an increasing pace of the epigenetic clock over time, per the Horvath (but not the Hannum) index of cellular aging.

CONCLUSIONS

This is the first study to suggest that posttraumatic psychopathology is longitudinally associated with a quickened pace of the epigenetic clock. Results raise the possibility that accelerated cellular aging is a common biological consequence of stress-related psychopathology, which carries implications for identifying mechanisms of stress-related cellular aging and developing interventions to slow its pace.

DNA methylation correlates of PTSD: Recent findings and technical challenges

There is increasing evidence that epigenetic factors play a critical role in posttraumatic stress disorder (PTSD), by mediating the impact of environmental exposures to trauma on the regulation of gene expression. DNA methylation is one epigenetic process that has been highly studied in PTSD. This review will begin by providing an overview of DNA methylation (DNAm) methods, and will then highlight two major biological systems that have been identified in the epigenetic regulation in PTSD: (a) the immune system and (b) the stress response system. In addition to candidate gene approaches, we will review novel strategies to study epigenome-wide PTSD-related effects, including epigenome-wide algorithms that distill information from many loci into a single summary score (e.g., measures of "epigenetic age" which have been associated with PTSD). This review will also cover recent epigenome wide association studies (EWAS) of PTSD, and biological pathway models used to identify gene sets enriched in PTSD. Finally, we address technical and methodological advances and challenges to the field, and highlight exciting directions for future research.

Reversing Behavioral, Neuroanatomical, and Germline Influences of Intergenerational Stress

BACKGROUND

Stressors affect populations exposed to them as well as offspring. Strategies preventing the intergenerational propagation of effects of stress would benefit public health. Olfactory cue-based fear conditioning provides a framework to address this issue.

METHODS

We 1) exposed adult male mice to an odor, acetophenone (Ace) or Lyral (parental generation [F0]-Exposed), 2) trained mice to associate these odors with mild foot shocks (F0-Trained), and 3) trained mice to associate these odors with mild foot shocks and then extinguished their fear toward these odors with odor-only presentations (F0-Extinguished). We then examined sensitivity of future generation (F1) offspring to these odors, expression of M71 odorant (Ace-responsive) and MOR23 odorant (Lyral-responsive) receptor-expressing cell populations in F1 offspring, and DNA methylation at genes encoding the Ace- (Olfr151, Olfr160) and Lyral- (Olfr16) responsive receptors in F0 sperm.

RESULTS

Extinguishing fear toward Ace or Lyral of F0 male mice (F0-Extinguished) that had been fear conditioned with Ace or Lyral, respectively, results in F1-Extinguished offspring that do not demonstrate behavioral sensitivity to Ace or Lyral, respectively, and do not have enhanced representation for M71 or MOR23 odorant receptors in the olfactory system, as is observed in F1-Trained-Ace or F1-Trained-Lyral cohorts, respectively. The promoters of genes encoding Olfr151 and Olfr160 receptors are less methylated in F0-Trained-Ace sperm compared with F0-Exposed-Ace sperm. The Olfr16 promoter is less methylated in F0-Trained-Lyral sperm compared with F0-Exposed-Lyral sperm, and F0-Extinguished-Lyral sperm have methylation levels comparable to F0-Exposed-Lyral sperm.

CONCLUSIONS

Our study demonstrates the potential of using extinction-based behavioral strategies to reverse influences of parental stress in offspring and in the parental germline.

Reduced interleukin 1A gene expression in the dorsolateral prefrontal cortex of individuals with PTSD and depression

The inflammatory system has been implicated in the pathophysiology of a variety of psychiatric conditions. Individuals with PTSD, depression, and other fear- and anxiety-related disorders exhibit alterations in peripheral circulating inflammatory markers, suggesting dysregulation of the inflammatory system. The relationship between inflammation and PTSD has been investigated almost exclusively in the periphery, and has not been extensively explored in human postmortem brain tissue. Interleukins (ILs) represent a subtype of cytokines and are key signaling proteins in the immune and inflammatory systems. Based on prior research implicating IL signaling in PTSD and depression, we performed a preliminary investigation of IL gene expression in a region of the cortex involved in emotion regulation and PTSD, the dorsolateral prefrontal cortex (dlPFC), using tissue from the newly established VA National PTSD Brain Bank. Gene expression analyses were conducted on post-mortem tissue from the dlPFC from 50 donors: 13 controls, 12 PTSD cases, and 25 depressed cases. RNA was extracted from frozen dlPFC tissue, reverse transcribed to cDNA, and quantitative polymerase chain reaction (qPCR) was performed to assess gene expression of IL1A, IL1B, IL6, IL8, IL10, IL13, and IL15. We found a multiple-testing corrected significant decrease in IL1A expression in the dlPFC for PTSD and depression cases compared to controls (p < 0.005) with age at death, sex, race and RNA integrity number (RIN) included as covariates. To our knowledge this finding is the first demonstration of altered IL expression in brain tissue from deceased individuals with histories of PTSD and/or depression.

Cell-type-specific interrogation of CeA Drd2 neurons to identify targets for pharmacological modulation of fear extinction

Behavioral and molecular characterization of cell-type-specific populations governing fear learning and behavior is a promising avenue for the rational identification of potential therapeutics for fear-related disorders. Examining cell-type-specific changes in neuronal translation following fear learning allows for targeted pharmacological intervention during fear extinction learning, mirroring possible treatment strategies in humans. Here we identify the central amygdala (CeA) Drd2-expressing population as a novel fear-supporting neuronal population that is molecularly distinct from other, previously identified, fear-supporting CeA populations. Sequencing of actively translating transcripts of Drd2 neurons using translating ribosome affinity purification (TRAP) technology identifies mRNAs that are differentially regulated following fear learning. Differentially expressed transcripts with potentially targetable gene products include Npy5r, Rxrg, Adora2a, Sst5r, Fgf3, Erbb4, Fkbp14, Dlk1, and Ssh3. Direct pharmacological manipulation of NPY5R, RXR, and ADORA2A confirms the importance of this cell population and these cell-type-specific receptors in fear behavior. Furthermore, these findings validate the use of functionally identified specific cell populations to predict novel pharmacological targets for the modulation of emotional learning.

Traumatic Stress and Accelerated Cellular Aging: From Epigenetics to Cardiometabolic Disease

PURPOSE OF REVIEW

The aim of this paper is to review the recent literature on traumatic stress-related accelerated aging, including a focus on cellular mechanisms and biomarkers of cellular aging and on the clinical manifestations of accelerated biological aging.

RECENT FINDINGS

Multiple lines of research converge to suggest that PTSD is associated with accelerated aging in the epigenome, and the immune and inflammation systems, and this may be reflected in premature onset of cardiometabolic and cardiovascular disease. The current state of research paves the way for future work focused on identifying the peripheral and central biological mechanisms linking traumatic stress to accelerated biological aging and medical morbidity, with an emphasis on processes involved in inflammation, immune functioning, oxidative stress, autonomic arousal, and stress response. Ultimately, such work could help reduce the pace of biological aging and improve health and wellness.

Genetic approaches for the study of PTSD: Advances and challenges

Post-traumatic stress disorder (PTSD) is a highly debilitating stress and anxiety-related disorder that occurs in response to specific trauma or abuse. Genetic risk factors may account for up to 30-40% of the heritability of PTSD. Understanding the gene pathways that are associated with PTSD, and how those genes interact with the fear and stress circuitry to mediate risk and resilience for PTSD will enable the development of targeted therapies to prevent the occurrence of or decrease the severity of this complex multi-gene disorder. This review will summarize recent research on genetic approaches to understanding PTSD risk and resilience in human populations, including candidate genes and their epigenetic modifications, genome-wide association studies and neural imaging genetics approaches. Despite challenges faced within this field of study such as inconsistent results and replications, genetic approaches still offer exciting opportunities for the identification and development of novel therapeutic targets and therapies in the future.

Bridging the Gap: Towards a cell-type specific understanding of neural circuits underlying fear behaviors

Fear and anxiety-related disorders are remarkably common and debilitating, and are often characterized by dysregulated fear responses. Rodent models of fear learning and memory have taken great strides towards elucidating the specific neuronal circuitries underlying the learning of fear responses. The present review addresses recent research utilizing optogenetic approaches to parse circuitries underlying fear behaviors. It also highlights the powerful advances made when optogenetic techniques are utilized in a genetically defined, cell-type specific, manner. The application of next-generation genetic and sequencing approaches in a cell-type specific context will be essential for a mechanistic understanding of the neural circuitry underlying fear behavior and for the rational design of targeted, circuit specific, pharmacologic interventions for the treatment and prevention of fear-related disorders.

From the neurobiology of extinction to improved clinical treatments

The neural circuitry underlying the fear response is extremely well conserved across mammalian species, which has allowed for the rapid translation of research findings in rodent models of fear to therapeutic interventions in human populations. Many aspects of exposure-based psychotherapy treatments in humans, which are widely used in the treatment of PTSD, panic disorder, phobias, and other anxiety disorders, are closely paralleled by extinction training in rodent fear conditioning models. Here, we discuss how the neural circuitry of fear learning and extinction in rodent animal models may be used to understand the underlying neural circuitry of fear-related disorders, such as PTSD in humans. We examine the factors that contribute to the pathology and development of PTSD. Next, we will review how fear is measured in animal models using classical Pavlovian fear conditioning paradigms, as well as brain regions such as the amygdala, which are involved in the fear response across species. Finally, we highlight the following three systems involved in the extinction of fear, all of which represent promising avenues for therapeutic interventions in the clinic: (1) the role of the glutamatergic N-methyl-d-aspartate (NMDA) receptor, (2) the role of the brain-derived neurotrophic factor (BDNF)-tyrosine kinase B (TrkB) induced signaling pathway, and (3) the role of the renin-angiotensin system. The modulation of pathways underlying fear learning and extinction, such as the ones presented in this review, in combination with extinction-based exposure therapy, represents promising avenues for therapeutic intervention in the treatment of human fear related disorders.