Abstract: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.

Decreased mononuclear cell NR3C1 SKA2 and FKPB5 expression levels among adult survivors of suicide bombing terror attacks in childhood are associated with the development of PTSD

Reduced Cingulate Gyrus Volume Associated with Enhanced Cortisol Awakening Response in Young Healthy Adults Reporting Childhood Trauma

Cell Type-Specific Gene Expression Patterns Associated with Posttraumatic Stress Disorder in World Trade Center Responders

- Inflammatory gene expression differences observed in depression treatment response groups: exploring the role of childhood trauma

Increased methylation of glucocorticoid receptor gene (NR3C1) in adults with a history of childhood maltreatment: a link with the severity and type of trauma

Transcriptome-wide association study of post-trauma symptom trajectories identified GRIN3B as a potential biomarker for PTSD development

Epigenetic regulation of the glucocorticoid receptor in human brain associates with childhood abuse

Early Childhood Traumas Interact with Variation in P2RX7 Gene Implicated in Neuroinflammation on Markers of Current Suicide Risk

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

NR3C1 and NR3C2 Genes Increase the Risk of Suicide Attempt in Psychiatric Disorder Patients with History of Childhood Trauma

Modeling gene × environment interactions in PTSD using human neurons reveals diagnosis-specific glucocorticoid-induced gene expression

Differential transcriptional response following glucocorticoid activation in cultured blood immune cells: a novel approach to PTSD biomarker development

Microglia-mediated neuroimmune suppression in PTSD is associated with anhedonia

Decreased Functional Connectivity of Hippocampal Subregions and Methylation of the NR3C1 Gene in Han Chinese Adults Who Lost Their Only Child

The relationship of maternal and child methylation of the glucocorticoid receptor NR3C1 during early childhood and subsequent child psychopathology at school-age in the context of maternal interpersonal violence-related post-traumatic stress disorder

PTSD is associated with accelerated transcriptional aging in World Trade Center responders

GWAS in interaction with childhood traumas implicates novel variants and genes previously associated with suicide-related factors in the background of suicidal ideation

Circulating cell-free mitochondrial DNA levels and glucocorticoid sensitivity in a cohort of male veterans with and without combat-related PTSD

Methylation and expression of glucocorticoid receptor exon-1 variants and FKBP5 in teenage suicide-completers

Epigenetic modification of the glucocorticoid receptor gene is linked to traumatic memory and post-traumatic stress disorder risk in genocide survivors

Two HPA axis genes, CRHBP and FKBP5, interact with childhood trauma to increase the risk for suicidal behavior