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

Systems biology dissection of PTSD and MDD across brain regions, cell types, and blood

The molecular pathology of stress-related disorders remains elusive. Our brain multiregion, multiomic study of posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) included the central nucleus of the amygdala, hippocampal dentate gyrus, and medial prefrontal cortex (mPFC). Genes and exons within the mPFC carried most disease signals replicated across two independent cohorts. Pathways pointed to immune function, neuronal and synaptic regulation, and stress hormones. Multiomic factor and gene network analyses provided the underlying genomic structure. Single nucleus RNA sequencing in dorsolateral PFC revealed dysregulated (stress-related) signals in neuronal and non-neuronal cell types. Analyses of brain-blood intersections in >50,000 UK Biobank participants were conducted along with fine-mapping of the results of PTSD and MDD genome-wide association studies to distinguish risk from disease processes. Our data suggest shared and distinct molecular pathology in both disorders and propose potential therapeutic targets and biomarkers.

Transcriptomic imputation of genetic risk variants uncovers novel whole-blood biomarkers of Parkinson's disease

Blood-based gene expression signatures could potentially be used as biomarkers for PD. However, it is unclear whether genetically-regulated transcriptomic signatures can provide novel gene candidates for use as PD biomarkers. We leveraged on the Genotype-Tissue Expression (GTEx) database to impute whole-blood transcriptomic expression using summary statistics of three large-scale PD GWAS. A random forest classifier was used with the consensus whole-blood imputed gene signature (IGS) to discriminate between cases and controls. Outcome measures included Area under the Curve (AUC) of Receiver Operating Characteristic (ROC) Curve. We demonstrated that the IGS (n = 37 genes) is conserved across PD GWAS studies and brain tissues. IGS discriminated between cases and controls in an independent whole-blood RNA-sequencing study (1176 PD, 254 prodromal, and 860 healthy controls) with mean AUC and accuracy of 64.8% and 69.4% for PD cohort, and 78.8% and 74% for prodromal cohort. PATL2 was the top-performing imputed gene in both PD and prodromal PD cohorts, whose classifier performance varied with biological sex (higher performance for males and females in the PD and prodromal PD, respectively). Single-cell RNA-sequencing studies (scRNA-seq) of healthy humans and PD patients found PATL2 to be enriched in terminal effector CD8+ and cytotoxic CD4+ cells, whose proportions are both increased in PD patients. We demonstrated the utility of GWAS transcriptomic imputation in identifying novel whole-blood transcriptomic signatures which could be leveraged upon for PD biomarker derivation. We identified PATL2 as a potential biomarker in both clinical and prodromic PD.

Detection of evolutionary conserved and accelerated genomic regions related to adaptation to thermal niches in Anolis lizards

Understanding the genetic basis for adapting to thermal environments is important due to serious effects of global warming on ectothermic species. Various genes associated with thermal adaptation in lizards have been identified mainly focusing on changes in gene expression or the detection of positively selected genes using coding regions. Only a few comprehensive genome-wide analyses have included noncoding regions. This study aimed to identify evolutionarily conserved and accelerated genomic regions using whole genomes of eight Anolis lizard species that have repeatedly adapted to similar thermal environments in multiple lineages. Evolutionarily conserved genomic regions were extracted as regions with overall sequence conservation (regions with fewer base substitutions) across all lineages compared with the neutral model. Genomic regions that underwent accelerated evolution in the lineage of interest were identified as those with more base substitutions in the target branch than in the entire background branch. Conserved elements across all branches were relatively abundant in "intergenic" genomic regions among noncoding regions. Accelerated regions (ARs) of each lineage contained a significantly greater proportion of noncoding RNA genes than the entire multiple alignment. Common genes containing ARs within 5 kb of their vicinity in lineages with similar thermal habitats were identified. Many genes associated with circadian rhythms and behavior were found in hot-open and cool-shaded habitat lineages. These genes might play a role in contributing to thermal adaptation and assist future studies examining the function of genes involved in thermal adaptation via genome editing.

Effects of sex and gender on the etiologies and presentation of select internalizing psychopathologies

The internalizing spectrum encompasses a subset of psychopathologies characterized by emotional liability, anhedonia, anxiousness, distress, and fear, and includes, among others, diagnoses of major depressive disorder (MDD), generalized anxiety disorder (GAD), and posttraumatic stress disorder (PTSD). In this review, we describe the vast body of work highlighting a role for sex and gender in the environment, symptom onset, genetic liability, and disorder progression and comorbidities of MDD, GAD, and PTSD. We also point the reader to different language used in diverse fields to describe sexual and gender minorities that may complicate the interpretation of emerging literature from the social sciences, psychiatric and psychological sciences, and genetics. Finally, we identify several gaps in knowledge that we hope serve as launch-points for expanding the scope of psychiatric studies beyond binarized sex-stratification. Despite being under-represented in genomics studies, placing emphasis on inclusion of sexual and gender diverse participants in these works will hopefully improve our understanding of disorder etiology using genetics as one tool to inform how biology (e.g., hormone concentration) and environmental variables (e.g., exposure to traumatic events) contribute to differences in symptom onset, pattern, and long-term trajectory.

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.

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.

Trauma Matters: Integrating Genetic and Environmental Components of PTSD

Trauma is ubiquitous, but only a subset of those who experience trauma will develop posttraumatic stress disorder (PTSD). In this review, it is argued that to determine who is at risk of developing PTSD, it is critical to examine the genetic etiology of the disorder and individual trauma profiles of those who are susceptible. First, the state of current PTSD genetic research is described, with a particular focus on studies that present evidence for trauma type specificity, or for differential genetic etiology according to gender or race. Next, approaches that leverage non-traditional phenotyping approaches are reviewed to identify PTSD-associated variants and biology, and the relative advantages and limitations inherent in these studies are reflected on. Finally, it is discussed how trauma might influence the heritability of PTSD, through type, risk factors, genetics, and associations with PTSD symptomology.

Identifying genetic loci that are associated with changes in gene expression in PTSD in a South African cohort

The molecular mechanisms underlying posttraumatic stress disorder (PTSD) are yet to be fully elucidated, especially in underrepresented population groups. Expression quantitative trait loci (eQTLs) are DNA sequence variants that influence gene expression, in a local (cis-) or distal (trans-) manner, and subsequently impact cellular, tissue, and system physiology. This study aims to identify genetic loci associated with gene expression changes in a South African PTSD cohort. Genome-wide genotype and RNA-sequencing data were obtained from 32 trauma-exposed controls and 35 PTSD cases of mixed-ancestry, as part of the SHARED ROOTS project. The first approach utilised 108 937 single-nucleotide polymorphisms (SNPs) (MAF > 10%) and 11 312 genes with Matrix eQTL to map potential eQTLs, while controlling for covariates as appropriate. The second analysis was focused on 5638 SNPs related to a previously calculated PTSD polygenic risk score for this cohort. SNP-gene pairs were considered eQTLs if they surpassed Bonferroni correction and had a false discovery rate <0.05. We did not identify eQTLs that significantly influenced gene expression in a PTSD-dependent manner. However, several known cis-eQTLs, independent of PTSD diagnosis, were observed. rs8521 (C > T) was associated with TAGLN and SIDT2 expression, and rs11085906 (C > T) was associated with ZNF333 expression. This exploratory study provides insight into the molecular mechanisms associated with PTSD in a non-European, admixed sample population. This study was limited by the cross-sectional design and insufficient statistical power. Overall, this study should encourage further multi-omics approaches towards investigating PTSD in diverse populations.

Treatment-associated mRNA co-expression changes in monocytes of patients with posttraumatic stress disorder

PTSD is a prevalent mental disorder that results from exposure to extreme and stressful life events and comes at high costs for both the individual and society. Therapeutic treatment presents the best way to deal with PTSD-the mechanisms underlying change after treatment, however, remain poorly understood. While stress and immune associated gene expression changes have been associated with PTSD development, studies investigating treatment effects at the molecular level so far tended to focus on DNA methylation. Here we use gene-network analysis on whole-transcriptome RNA-Seq data isolated from CD14⁺ monocytes of female PTSD patients (N = 51) to study pre-treatment signatures of therapy response and therapy-related changes at the level of gene expression. Patients who exhibited significant symptom improvement after therapy showed higher baseline expression in two modules involved in inflammatory processes (including notable examples IL1R2 and FKBP5) and blood coagulation. After therapy, expression of an inflammatory module was increased, and expression of a wound healing module was decreased. This supports findings reporting an association between PTSD and dysregulations of the inflammatory and the hemostatic system and mark both as potentially treatment sensitive.

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.

Stem Cell Models for Context-Specific Modeling in Psychiatric Disorders

Genome-wide association studies reveal the complex polygenic architecture underlying psychiatric disorder risk, but there is an unmet need to validate causal variants, resolve their target genes(s), and explore their functional impacts on disorder-related mechanisms. Disorder-associated loci regulate transcription of target genes in a cell type- and context-specific manner, which can be measured through expression quantitative trait loci. In this review, we discuss methods and insights from context-specific modeling of genetically and environmentally regulated expression. Human induced pluripotent stem cell-derived cell type and organoid models have uncovered context-specific psychiatric disorder associations by investigating tissue-, cell type-, sex-, age-, and stressor-specific genetic regulation of expression. Techniques such as massively parallel reporter assays and pooled CRISPR (clustered regularly interspaced short palindromic repeats) screens make it possible to functionally fine-map genome-wide association study loci and validate their target genes at scale. Integration of disorder-associated contexts with these patient-specific human induced pluripotent stem cell models makes it possible to uncover gene by environment interactions that mediate disorder risk, which will ultimately improve our ability to diagnose and treat psychiatric disorders.

Transcriptional signal and cell specificity of genes related to cortical structural differences of post-traumatic stress disorder

Due to the diversity of traumatic events, the diagnosis of Post-traumatic Stress Disorder is heterogeneous. The pathogenesis has been explored in the fields of brain imaging and genomics separately, but the results are inconsistent. Previous research evidenced that there existed structural differences between PTSD and healthy controls in multiple brain regions. This study further looked into the differences of brain structure in PTSD at the whole brain level and analyzed the difference-related genomes. The brain structure imaging data of 36 patients and 32 healthy controls were taken as morphological indexes. Partial least squares regression and transcriptome data were used to extract genomes related to structural differences. Additional data sets were used to study transcription characteristics of genome. Morphological differences were found in cingulate gyrus between patients and control group. Differentially expressed genes related to Morphometric similarity networks difference space were also observed. The obtained genes (i.e., RORA, PRKG1 and FKBP5) were proved to be related to the disorder with no significant correlation with other mental illnesses. In the subsequent cell type analysis, astrocytes, excitatory neurons and inhibitory neurons were evidenced to have the most significant correlation with these genes. This study found morphologically different brain regions related to PTSD. The related genome transcription analysis connects the structural differences and molecular mechanisms.

Polygenic risk of Social-isolation and its influence on social behavior, psychosis, depression and autism spectrum disorder

Social-isolation has been linked to a range of psychiatric issues, but the behavioral component that drives it is not well understood. Here, a GWAS is carried out to identify genetic variants which contribute to Social-isolation behaviors in up to 449,609 participants from the UK Biobank. 17 loci were identified at genome-wide significance, contributing to a 4% SNP heritability estimate. Using the Social-isolation GWAS, polygenic risk scores (PRS) were derived in ALSPAC, an independent, developmental cohort, and used to test for association with friendship quality. At age 18, friendship scores were associated with the Social-isolation PRS, demonstrating that the genetic factors are able to predict related social traits. LD score regression using the GWAS demonstrated genetic correlation with autism spectrum disorder, schizophrenia, and major depressive disorder. However, no evidence of causality was found using a conservative Mendelian randomization approach other than that of autism spectrum disorder on Social-isolation. Our results show that Social-isolation has a small heritable component which may drive those behaviors which is associated genetically with other social traits such as friendship satisfaction as well as psychiatric disorders.

Aberrant white matter microstructure evaluation by automated fiber quantification in typhoon-related post-traumatic stress disorder

Super typhoons can lead to post-traumatic stress disorder (PTSD), which can adversely affect a person's mental health after a disaster. Neuroimaging studies suggest that patients with PTSD may have post-exposure abnormalities of the white matter. However, little is known about these defects, if they are localized to specific regions of the white matter fibers, or whether they may be potential biomarkers for PTSD. Typhoon survivors with PTSD (n = 27), trauma-exposed controls (TEC) (n = 33), and healthy controls (HCs) (n = 30) were enrolled. We used automated fiber quantification (AFQ) to process the participants' DTI and compared diffusion metrics among the three groups. To evaluate diagnostic value, we used support vector machine (SVM) and a random forest (RF) classifier to build a machine learning model. White matter fiber segmentation between the three groups was found to be statistically significant for the fractional anisotropy (FA) value of the right anterior thalamic radiation (ATR) (26-50 nodes) and right uncinate fasciculus (UF) (60-72 nodes) (FDR correction, p < 0.05). By analyzing the characteristics of the machine learning model, the two most important variables were the right ATR and right UF for differentiating PTSD and trauma-exposed controls (TEC) from the healthy controls (HC). In addition, the left cingulum cingulate and left UF were the most critical variables in the differentiation of PTSD and TEC. AFQ with machine learning can localize abnormalities in specific regions of white matter fibers. These regions may be used as a diagnostic biomarker for PTSD.

A translational genomics approach identifies IL10RB as the top candidate gene target for COVID-19 susceptibility

Recent efforts have identified genetic loci that are associated with coronavirus disease 2019 (COVID-19) infection rates and disease outcome severity. Translating these genetic findings into druggable genes that reduce COVID-19 host susceptibility is a critical next step. Using a translational genomics approach that integrates COVID-19 genetic susceptibility variants, multi-tissue genetically regulated gene expression (GReX), and perturbagen signatures, we identified IL10RB as the top candidate gene target for COVID-19 host susceptibility. In a series of validation steps, we show that predicted GReX upregulation of IL10RB and higher IL10RB expression in COVID-19 patient blood is associated with worse COVID-19 outcomes and that in vitro IL10RB overexpression is associated with increased viral load and activation of disease-relevant molecular pathways.

Integrating human brain proteomes with genome-wide association data implicates novel proteins in post-traumatic stress disorder

Genome-wide association studies (GWAS) have identified several risk loci for post-traumatic stress disorder (PTSD); however, how they confer PTSD risk remains unclear. We aimed to identify genes that confer PTSD risk through their effects on brain protein abundance to provide new insights into PTSD pathogenesis. To that end, we integrated human brain proteomes with PTSD GWAS results to perform a proteome-wide association study (PWAS) of PTSD, followed by Mendelian randomization, using a discovery and confirmatory study design. Brain proteomes (N = 525) were profiled from the dorsolateral prefrontal cortex using mass spectrometry. The Million Veteran Program (MVP) PTSD GWAS (n = 186,689) was used for the discovery PWAS, and the Psychiatric Genomics Consortium PTSD GWAS (n = 174,659) was used for the confirmatory PWAS. To understand whether genes identified at the protein-level were also evident at the transcript-level, we performed a transcriptome-wide association study (TWAS) using human brain transcriptomes (N = 888) and the MVP PTSD GWAS results. We identified 11 genes that contribute to PTSD pathogenesis via their respective cis-regulated brain protein abundance. Seven of 11 genes (64%) replicated in the confirmatory PWAS and 4 of 11 also had their cis-regulated brain mRNA levels associated with PTSD. High confidence level was assigned to 9 of 11 genes after considering evidence from the confirmatory PWAS and TWAS. Most of the identified genes are expressed in other PTSD-relevant brain regions and several are preferentially expressed in excitatory neurons, astrocytes, and oligodendrocyte precursor cells. These genes are novel, promising targets for mechanistic and therapeutic studies to find new treatments for PTSD.

Didymin Ameliorates Liver Fibrosis by Alleviating Endoplasmic Reticulum Stress and Glycerophospholipid Metabolism: Based on Transcriptomics and Metabolomics

Introduction

Origanum vulgare L. is a traditional Chinese herb, having a strong hepatoprotective effect. In our previous experiments, we have isolated an ingredient from this herb and identified it as didymin. This study aimed to investigate the effects and underlying mechanisms of didymin on liver injury and fibrosis, elucidating whether it was the pharmacodynamic material basis of Origanum vulgare L.

Methods

Mice were injected with CCl₄ for 10 weeks to induce liver fibrosis, followed by didymin treatment for 6 weeks. Then, biochemical analysis and histopathological examinations were conducted to evaluate the therapeutic effects of didymin in alleviating fibrosis. Next, the possible mechanisms of didymin were predicted by transcriptomics and then verified by the multiple relevant examinations.

Results

The pharmacodynamic experiments indicated that didymin significantly attenuated CCl₄-induced hepatic injury and fibrogenesis, as evidenced by the ameliorative pathological tissue, low transaminase activity, and decreased collagen accumulation. Interestingly, the transcriptome analysis predicted that the potential targets were likely to be endoplasmic reticulum stress (ERS), inflammation, apoptosis, and metabolic pathways. And the predictions were then verified by the following examinations: (1) didymin significantly inhibited ERS by regulating the ATF6, IRE1α, and PERK pathways; (2) didymin markedly alleviated hepatocyte apoptosis by restoring the expression of Bcl-2 and caspase families, as well as the mitochondrial dysfunction; (3) didymin significantly decreased the production of the pro-inflammatory cytokines (IL-1β and IL-6); (4) didymin inhibited the glycerophospholipid metabolism pathway by decreasing the synthesis of phosphatidylethanolamines and phosphatidylcholines.

Conclusion

Our findings demonstrate that didymin can ameliorate liver fibrosis, which is mainly attributed to the inhibition of ERS, inflammation, and glycerophospholipid metabolism.

Post-traumatic stress disorder: clinical and translational neuroscience from cells to circuits

Post-traumatic stress disorder (PTSD) is a maladaptive and debilitating psychiatric disorder, characterized by re-experiencing, avoidance, negative emotions and thoughts, and hyperarousal in the months and years following exposure to severe trauma. PTSD has a prevalence of approximately 6-8% in the general population, although this can increase to 25% among groups who have experienced severe psychological trauma, such as combat veterans, refugees and victims of assault. The risk of developing PTSD in the aftermath of severe trauma is determined by multiple factors, including genetics - at least 30-40% of the risk of PTSD is heritable - and past history, for example, prior adult and childhood trauma. Many of the primary symptoms of PTSD, including hyperarousal and sleep dysregulation, are increasingly understood through translational neuroscience. In addition, a large amount of evidence suggests that PTSD can be viewed, at least in part, as a disorder that involves dysregulation of normal fear processes. The neural circuitry underlying fear and threat-related behaviour and learning in mammals, including the amygdala-hippocampus-medial prefrontal cortex circuit, is among the most well-understood in behavioural neuroscience. Furthermore, the study of threat-responding and its underlying circuitry has led to rapid progress in understanding learning and memory processes. By combining molecular-genetic approaches with a translational, mechanistic knowledge of fear circuitry, transformational advances in the conceptual framework, diagnosis and treatment of PTSD are possible. In this Review, we describe the clinical features and current treatments for PTSD, examine the neurobiology of symptom domains, highlight genomic advances and discuss translational approaches to understanding mechanisms and identifying new treatments and interventions for this devastating syndrome.

DNA methylation of Nuclear Factor of Activated T Cells 1 mediates the prospective relation between exposure to different traumatic event types and post-traumatic stress disorder

The mechanisms through which exposure to differing trauma types become biologically embedded to shape the risk for post-traumatic stress disorder (PTSD) is unclear. DNA methylation (5-mC), particularly in stress-relevant genes, may play a role in this relationship. Here, we conducted path analysis using generalized structural equation modeling to investigate whether blood-derived 5-mC in Nuclear Factor of Activated T Cells 1 (NFATC1) mediates the prospective association between each of five different trauma types ("assaultive violence", "other injury or shocking experience", "learning of trauma to loved one", "sudden, unexpected death of a close friend or relative", and "other") and lifetime PTSD. All five trauma types were significantly associated with reduced methylation at NFATC1 CpG site, cg17057218. Two of the five trauma types were significantly associated with increased methylation at NFATC1 CpG site, cg22324981. Moreover, methylation at cg17057218 significantly mediated 21-32% of the total effect for four of the five trauma types, while methylation at cg22324981 mediated 27-40% of the total effect for two of the five trauma types. These CpG sites were differentially associated with transcription factor binding sites and chromatin state signatures. NFATC1 5-mC may be a potential mechanism in the relationship between some trauma types and prospective risk for PTSD.

Altered gene expression and PTSD symptom dimensions in World Trade Center responders

Despite experiencing a significant trauma, only a subset of World Trade Center (WTC) rescue and recovery workers developed posttraumatic stress disorder (PTSD). Identification of biomarkers is critical to the development of targeted interventions for treating disaster responders and potentially preventing the development of PTSD in this population. Analysis of gene expression from these individuals can help in identifying biomarkers of PTSD. We established a well-phenotyped sample of 371 WTC responders, recruited from a longitudinal WTC responder cohort using stratified random sampling, by obtaining blood, self-reported and clinical interview data. Using bulk RNA-sequencing from whole blood, we examined the association between gene expression and WTC-related PTSD symptom severity on (i) highest lifetime Clinician-Administered PTSD Scale (CAPS) score, (ii) past-month CAPS score, and (iii) PTSD symptom dimensions using a 5-factor model of re-experiencing, avoidance, emotional numbing, dysphoric arousal and anxious arousal symptoms. We corrected for sex, age, genotype-derived principal components and surrogate variables. Finally, we performed a meta-analysis with existing PTSD studies (total N = 1016), using case/control status as the predictor and correcting for these variables. We identified 66 genes significantly associated with total highest lifetime CAPS score (FDR-corrected p < 0.05), and 31 genes associated with total past-month CAPS score. Our more granular analyses of PTSD symptom dimensions identified additional genes that did not reach statistical significance in our analyses with total CAPS scores. In particular, we identified 82 genes significantly associated with lifetime anxious arousal symptoms. Several genes significantly associated with multiple PTSD symptom dimensions and total lifetime CAPS score (SERPINA1, RPS6KA1, and STAT3) have been previously associated with PTSD. Geneset enrichment of these findings has identified pathways significant in metabolism, immune signaling, other psychiatric disorders, neurological signaling, and cellular structure. Our meta-analysis revealed 10 genes that reached genome-wide significance, all of which were downregulated in cases compared to controls (CIRBP, TMSB10, FCGRT, CLIC1, RPS6KB2, HNRNPUL1, ALDOA, NACA, ZNF429 and COPE). Additionally, cellular deconvolution highlighted an enrichment in CD4 T cells and eosinophils in responders with PTSD compared to controls. The distinction in significant genes between total lifetime CAPS score and the anxious arousal symptom dimension of PTSD highlights a potential biological difference in the mechanism underlying the heterogeneity of the PTSD phenotype. Future studies should be clear about methods used to analyze PTSD status, as phenotypes based on PTSD symptom dimensions may yield different gene sets than combined CAPS score analysis. Potential biomarkers implicated from our meta-analysis may help improve therapeutic target development for PTSD.

Mineralocorticoid receptor and glucocorticoid receptor work alone and together in cell-type-specific manner: Implications for resilience prediction and targeted therapy

‘You can't roll the clock back and reverse the effects of experiences' Bruce McEwen used to say when explaining how allostasis labels the adaptive process. Here we will for once roll the clock back to the times that the science of the glucocorticoid hormone was honored with a Nobel prize and highlight the discovery of their receptors in the hippocampus as inroad to its current status as master regulator in control of stress coping and adaptation. Glucocorticoids operate in concert with numerous neurotransmitters, neuropeptides, and other hormones with the aim to facilitate processing of information in the neurocircuitry of stress, from anticipation and perception of a novel experience to behavioral adaptation and memory storage. This action, exerted by the glucocorticoids, is guided by two complementary receptor systems, mineralocorticoid receptors (MR) and glucocorticoid receptors (GR), that need to be balanced for a healthy stress response pattern. Here we discuss the cellular, neuroendocrine, and behavioral studies underlying the MR:GR balance concept, highlight the relevance of hypothalamic-pituitary-adrenal (HPA) -axis patterns and note the limited understanding yet of sexual dimorphism in glucocorticoid actions. We conclude with the prospect that (i) genetically and epigenetically regulated receptor variants dictate cell-type-specific transcriptome signatures of stress-related neuropsychiatric symptoms and (ii) selective receptor modulators are becoming available for more targeted treatment. These two new developments may help to ‘restart the clock’ with the prospect to support resilience.

Enhancing Discovery of Genetic Variants for Posttraumatic Stress Disorder Through Integration of Quantitative Phenotypes and Trauma Exposure Information

BACKGROUND

Posttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).

METHODS

A GWAS on PTSD symptoms was performed in 51 cohorts followed by a fixed-effects meta-analysis (N = 182,199 European ancestry participants). A GWAS of LTE burden was performed in the UK Biobank cohort (N = 132,988). Genetic correlations were evaluated with linkage disequilibrium score regression. Multivariate analysis was performed using Multi-Trait Analysis of GWAS. Functional mapping and annotation of leading loci was performed with FUMA. Replication was evaluated using the Million Veteran Program GWAS of PTSD total symptoms.

RESULTS

GWASs of PTSD symptoms and LTE burden identified 5 and 6 independent genome-wide significant loci, respectively. There was a 72% genetic correlation between PTSD and LTE. PTSD and LTE showed largely similar patterns of genetic correlation with other traits, albeit with some distinctions. Adjusting PTSD for LTE reduced PTSD heritability by 31%. Multivariate analysis of PTSD and LTE increased the effective sample size of the PTSD GWAS by 20% and identified 4 additional loci. Four of these 9 PTSD loci were independently replicated in the Million Veteran Program.

CONCLUSIONS

Through using a quantitative trait measure of PTSD, we identified novel risk loci not previously identified using prior case-control analyses. PTSD and LTE have a high genetic overlap that can be leveraged to increase discovery power through multivariate methods.

Posttraumatic Stress Disorder Brain Transcriptomics: Convergent Genomic Signatures Across Biological Sex

While a definitive understanding of the molecular pathology of posttraumatic stress disorder (PTSD) is far from a current reality, it has become increasingly clear that many of the molecular effects of PTSD are sex specific. Women are twice as likely as men to develop PTSD after a traumatic event, and neurobiological evidence suggests that there are structural differences between the brains of males versus females with PTSD. Recent advances in genomic technologies have begun to shed light on the sex-specific molecular determinants of PTSD, which seem to be governed predominantly by dysfunction of GABAergic (gamma-aminobutyric acidergic) signaling and immune function. We review the current state of the field of PTSD genomics focusing on the effect of sex. We provide an overview of difference in heritability of PTSD based on sex, how difference in gene regulation based on sex impacts the PTSD brain, and what is known about genomic regulation that is dysregulated in specific cell types in PTSD.

Contributions of PTSD polygenic risk and environmental stress to suicidality in preadolescents

Suicidal ideation and attempts (i.e., suicidality) are complex behaviors driven by environmental stress, genetic susceptibility, and their interaction. Preadolescent suicidality is a major health problem with rising rates, yet its underlying biology is understudied. Here we studied effects of genetic stress susceptibility, approximated by the polygenic risk score (PRS) for post-traumatic-stress-disorder (PTSD), on preadolescent suicidality in participants from the Adolescent Brain Cognitive Development (ABCD) Study®. We further evaluated PTSD-PRS effects on suicidality in the presence of environmental stressors that are established suicide risk factors. Analyses included both European and African ancestry participants using PRS calculated based on summary statistics from ancestry-specific genome-wide association studies. In European ancestry participants (N = 4,619, n = 378 suicidal), PTSD-PRS was associated with preadolescent suicidality (odds ratio [OR] = 1.12, 95%CI 1-1.25, p = 0.038). Results in African ancestry participants (N = 1,334, n = 130 suicidal) showed a similar direction but were not statistically significant (OR = 1.21, 95%CI 0.93-1.57, p = 0.153). Sensitivity analyses using non-psychiatric polygenic score for height and using cross-ancestry PTSD-PRS did not reveal any association with suicidality, supporting the specificity of the association of ancestry-specific PTSD-PRS with suicidality. Environmental stressors were robustly associated with suicidality across ancestries with moderate effect size for negative life events and family conflict (OR 1.27-1.6); and with large effect size (OR ∼ 4) for sexual-orientation discrimination. When combined with environmental factors, PTSD-PRS showed marginal additive effects in explaining variability in suicidality, with no evidence for G × E interaction. Results support use of cross-phenotype PRS, specifically stress-susceptibility, as a genetic marker for suicidality risk early in the lifespan.

Transcriptomics of the depressed and PTSD brain

Stress is the response of an organism to demands for change, yet excessive or chronic stress contributes to nearly all psychiatric disorders. The advent of high-throughput transcriptomic methods such as single cell RNA sequencing poses new opportunities to understand the neurobiology of stress, yet substantial barriers to understanding stress remain. Stress adaptation is an organismal survival mechanism conserved across all organisms, yet there is an infinity of potential stressful experiences. Unraveling shared and separate transcriptional programs for adapting to stressful experience remains a challenge, despite methodological and analytic advances. Here we review the state of the field focusing on the technologies used to study the transcriptome for the stress neurobiologist, and also attempt to identify central questions about the heterogeneity of stress for those applying transcriptomic approaches. We further explore how postmortem transcriptome studies aided by preclinical animal models are converging on common molecular pathways for adaptation to aversive experience. Finally, we discuss approaches to integrate large genomic datasets with human neuroimaging and other datasets.

From genetics to systems biology of stress-related mental disorders

Many individuals will be exposed to some form of traumatic stress in their lifetime which, in turn, increases the likelihood of developing stress-related disorders such as post-traumatic stress disorder (PTSD), major depressive disorder (MDD) and anxiety disorders (ANX). The development of these disorders is also influenced by genetics and have heritability estimates ranging between ∼30 and 70%. In this review, we provide an overview of the findings of genome-wide association studies for PTSD, depression and ANX, and we observe a clear genetic overlap between these three diagnostic categories. We go on to highlight the results from transcriptomic and epigenomic studies, and, given the multifactorial nature of stress-related disorders, we provide an overview of the gene-environment studies that have been conducted to date. Finally, we discuss systems biology approaches that are now seeing wider utility in determining a more holistic view of these complex disorders.

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.

Posttraumatic stress disorder: from gene discovery to disease biology

Posttraumatic stress disorder (PTSD) is a complex mental disorder afflicting approximately 7% of the population. The diverse number of traumatic events and the wide array of symptom combinations leading to PTSD diagnosis contribute substantial heterogeneity to studies of the disorder. Genomic and complimentary-omic investigations have rapidly increased our understanding of the heritable risk for PTSD. In this review, we emphasize the contributions of genome-wide association, epigenome-wide association, transcriptomic, and neuroimaging studies to our understanding of PTSD etiology. We also discuss the shared risk between PTSD and other complex traits derived from studies of causal inference, co-expression, and brain morphological similarities. The investigations completed so far converge on stark contrasts in PTSD risk between sexes, partially attributed to sex-specific prevalence of traumatic experiences with high conditional risk of PTSD. To further understand PTSD biology, future studies should focus on detecting risk for PTSD while accounting for substantial cohort-level heterogeneity (e.g. civilian v. combat-exposed PTSD cases or PTSD risk among cases exposed to specific traumas), expanding ancestral diversity among study cohorts, and remaining cognizant of how these data influence social stigma associated with certain traumatic events among underrepresented minorities and/or high-risk populations.

Gene Expression Analysis in Three Posttraumatic Stress Disorder Cohorts Implicates Inflammation and Innate Immunity Pathways and Uncovers Shared Genetic Risk With Major Depressive Disorder

Posttraumatic stress disorder (PTSD) is a complex psychiatric disorder that can develop following exposure to traumatic events. The Psychiatric Genomics Consortium PTSD group (PGC-PTSD) has collected over 20,000 multi-ethnic PTSD cases and controls and has identified both genetic and epigenetic factors associated with PTSD risk. To further investigate biological correlates of PTSD risk, we examined three PGC-PTSD cohorts comprising 977 subjects to identify differentially expressed genes among PTSD cases and controls. Whole blood gene expression was quantified with the HumanHT-12 v4 Expression BeadChip for 726 OEF/OIF veterans from the Veterans Affairs (VA) Mental Illness Research Education and Clinical Center (MIRECC), 155 samples from the Injury and Traumatic Stress (INTRuST) Clinical Consortium, and 96 Australian Vietnam War veterans. Differential gene expression analysis was performed in each cohort separately followed by meta-analysis. In the largest cohort, we performed co-expression analysis to identify modules of genes that are associated with PTSD and MDD. We then conducted expression quantitative trait loci (eQTL) analysis and assessed the presence of eQTL interactions involving PTSD and major depressive disorder (MDD). Finally, we utilized PTSD and MDD GWAS summary statistics to identify regions that colocalize with eQTLs. Although not surpassing correction for multiple testing, the most differentially expressed genes in meta-analysis were interleukin-1 beta (IL1B), a pro-inflammatory cytokine previously associated with PTSD, and integrin-linked kinase (ILK), which is highly expressed in brain and can rescue dysregulated hippocampal neurogenesis and memory deficits. Pathway analysis revealed enrichment of toll-like receptor (TLR) and interleukin-1 receptor genes, which are integral to cellular innate immune response. Co-expression analysis identified four modules of genes associated with PTSD, two of which are also associated with MDD, demonstrating common biological pathways underlying the two conditions. Lastly, we identified four genes (UBA7, HLA-F, HSPA1B, and RERE) with high probability of a shared causal eQTL variant with PTSD and/or MDD GWAS variants, thereby providing a potential mechanism by which the GWAS variant contributes to disease risk. In summary, we provide additional evidence for genes and pathways previously reported and identified plausible novel candidates for PTSD. These data provide further insight into genetic factors and pathways involved in PTSD, as well as potential regions of pleiotropy between PTSD and MDD.

Unveiling the Pathogenesis of Psychiatric Disorders Using Network Models

Psychiatric disorders are complex brain disorders with a high degree of genetic heterogeneity, affecting millions of people worldwide. Despite advances in psychiatric genetics, the underlying pathogenic mechanisms of psychiatric disorders are still largely elusive, which impedes the development of novel rational therapies. There has been accumulating evidence suggesting that the genetics of complex disorders can be viewed through an omnigenic lens, which involves contextualizing genes in highly interconnected networks. Thus, applying network-based multi-omics integration methods could cast new light on the pathophysiology of psychiatric disorders. In this review, we first provide an overview of the recent advances in psychiatric genetics and highlight gaps in translating molecular associations into mechanistic insights. We then present an overview of network methodologies and review previous applications of network methods in the study of psychiatric disorders. Lastly, we describe the potential of such methodologies within a multi-tissue, multi-omics approach, and summarize the future directions in adopting diverse network approaches.

Putative Blood Somatic Mutations in Post-Traumatic Stress Disorder-Symptomatic Soldiers: High Impact of Cytoskeletal and Inflammatory Proteins

Background: We recently discovered autism/intellectual disability somatic mutations in postmortem brains, presenting higher frequency in Alzheimer’s disease subjects, compared with the controls. We further revealed high impact cytoskeletal gene mutations, coupled with potential cytoskeleton-targeted repair mechanisms. Objective: The current study was aimed at further discerning if somatic mutations in brain diseases are presented only in the most affected tissue (the brain), or if blood samples phenocopy the brain, toward potential diagnostics. Methods: Variant calling analyses on an RNA-seq database including peripheral blood samples from 85 soldiers (58 controls and 27 with symptoms of post-traumatic stress disorder, PTSD) was performed. Results: High (e.g., protein truncating) as well as moderate impact (e.g., single amino acid change) germline and putative somatic mutations in thousands of genes were found. Further crossing the mutated genes with autism, intellectual disability, cytoskeleton, inflammation, and DNA repair databases, identified the highest number of cytoskeletal-mutated genes (187 high and 442 moderate impact). Most of the mutated genes were shared and only when crossed with the inflammation database, more putative high impact mutated genes specific to the PTSD-symptom cohorts versus the controls (14 versus 13) were revealed, highlighting tumor necrosis factor specifically in the PTSD-symptom cohorts. Conclusion: With microtubules and neuro-immune interactions playing essential roles in brain neuroprotection and Alzheimer-related neurodegeneration, the current mutation discoveries contribute to mechanistic understanding of PTSD and brain protection, as well as provide future diagnostics toward personalized military deployment strategies and drug design.

Genome-wide association analyses of post-traumatic stress disorder and its symptom subdomains in the Million Veteran Program

We conducted genome-wide association analyses of over 250,000 participants of European (EUR) and African (AFR) ancestry from the Million Veteran Program using electronic health record-validated post-traumatic stress disorder (PTSD) diagnosis and quantitative symptom phenotypes. Applying genome-wide multiple testing correction, we identified three significant loci in European case-control analyses and 15 loci in quantitative symptom analyses. Genomic structural equation modeling indicated tight coherence of a PTSD symptom factor that shares genetic variance with a distinct internalizing (mood-anxiety-neuroticism) factor. Partitioned heritability indicated enrichment in several cortical and subcortical regions, and imputed genetically regulated gene expression in these regions was used to identify potential drug repositioning candidates. These results validate the biological coherence of the PTSD syndrome, inform its relationship to comorbid anxiety and depressive disorders and provide new considerations for treatment.

Massively parallel techniques for cataloguing the regulome of the human brain

Complex brain disorders are highly heritable and arise from a complex polygenic risk architecture. Many disease-associated loci are found in non-coding regions that house regulatory elements. These elements influence the transcription of target genes-many of which demonstrate cell-type-specific expression patterns-and thereby affect phenotypically relevant molecular pathways. Thus, cell-type-specificity must be considered when prioritizing candidate risk loci, variants and target genes. This Review discusses the use of high-throughput assays in human induced pluripotent stem cell-based neurodevelopmental models to probe genetic risk in a cell-type- and patient-specific manner. The application of massively parallel reporter assays in human induced pluripotent stem cells can characterize the human regulome and test the transcriptional responses of putative regulatory elements. Parallel CRISPR-based screens can further functionally dissect this genetic regulatory architecture. The integration of these emerging technologies could decode genetic risk into medically actionable information, thereby improving genetic diagnosis and identifying novel points of therapeutic intervention.

Exploring the impact of trauma type and extent of exposure on posttraumatic alterations in 5-HT1A expression

The long-term behavioral, psychological, and neurobiological effects of exposure to potentially traumatic events vary within the human population. Studies conducted on trauma-exposed human subjects suggest that differences in trauma type and extent of exposure combine to affect development, maintenance, and treatment of a variety of psychiatric syndromes. The serotonin 1-A receptor (5-HT1A) is an inhibitory G protein-coupled serotonin receptor encoded by the HTR1A gene that plays a role in regulating serotonin release, physiological stress responding, and emotional behavior. Studies from the preclinical and human literature suggest that dysfunctional expression of 5-HT1A is associated with a multitude of psychiatric symptoms commonly seen in trauma-exposed individuals. Here, we synthesize the literature, including numerous preclinical studies, examining differences in alterations in 5-HT1A expression following trauma exposure. Collectively, these findings suggest that the impact of trauma exposure on 5-HT1A expression is dependent, in part, on trauma type and extent of exposure. Furthermore, preclinical and human studies suggest that this observation likely applies to additional molecular targets and may help explain variation in trauma-induced changes in behavior and treatment responsivity. In order to understand the neurobiological impact of trauma, including the impact on 5-HT1A expression, it is crucial to consider both trauma type and extent of exposure.