Mapping the transcriptomics landscape of post-traumatic stress disorder symptom dimensions in World Trade Center responders

Gene expression has provided promising insights into the pathophysiology of post-traumatic stress disorder (PTSD); however, specific regulatory transcriptomic mechanisms remain unknown. The present study addressed this limitation by performing transcriptome-wide RNA-Seq of whole-blood samples from 226 World Trade Center responders. The investigation focused on differential expression (DE) at the gene, isoform, and for the first time, alternative splicing (AS) levels associated with the symptoms of PTSD: total burden, re-experiencing, avoidance, numbing, and hyperarousal subdimensions. These symptoms were associated with 76, 1, 48, 15, and 49 DE genes, respectively (FDR < 0.05). Moreover, they were associated with 103, 11, 0, 43, and 32 AS events. Avoidance differed the most from other dimensions with respect to DE genes and AS events. Gene set enrichment analysis (GSEA) identified pathways involved in inflammatory and metabolic processes, which may have implications in the treatment of PTSD. Overall, the findings shed a novel light on the wide range of transcriptomic alterations associated with PTSD at the gene and AS levels. The results of DE analysis associated with PTSD subdimensions highlights the importance of studying PTSD symptom heterogeneity.