Skip to main content

ORIGINAL RESEARCH article

Front. Public Health, 04 May 2023
Sec. Public Mental Health
This article is part of the Research Topic Challenges and Emerging Issues on Firefighter’s Toxic Chemical Exposure: Smoke Chemicals, Contaminated PPE, and Off-gassing View all 9 articles

Post-traumatic stress disorder and depressive symptoms among firefighters: a network analysis

Peng ChengPeng Cheng1Lirong WangLirong Wang2Ying ZhouYing Zhou1Wenjing MaWenjing Ma1Guangju ZhaoGuangju Zhao1Li ZhangLi Zhang1Weihui Li
Weihui Li1*
  • 1Department of Psychiatry, and National Clinical Research Center for Mental Disorders, The Second Xiangya Hospital of Central South University, Changsha, Hunan, China
  • 2Xiangya School of Medicine, Xiangya Hospital, Central South University, Changsha, China

Background: Firefighters, as first responders with a high risk of occupational exposure to traumatic events and heavy working stress, have a high prevalence of PTSD symptoms and depressive symptoms. But no previous studies analyzed the relationships and hierarchies of PTSD and depressive symptoms among firefighters. Network analysis is a novel and effective method for investigating the complex interactions of mental disorders at the symptom level and providing a new understanding of psychopathology. The current study was designed to characterize the PTSD and depressive symptoms network structure in the Chinese firefighters.

Method: The Primary Care PTSD Screen for DSM-5 (PC-PTSD-5) and the Self-Rating Depression Scale (SDS) were applied to assess PTSD and depressive symptoms, respectively. The network structure of PTSD and depressive symptoms was characterized using “expected influence (EI)” and “bridge EI” as centrality indices. The Walktrap algorithm was conducted to identify communities in the PTSD and depressive symptoms network. Finally, Network accuracy and stability were examined using the Bootstrapped test and the case-dropping procedure.

Results: A total of 1,768 firefighters were enrolled in our research. Network analysis revealed that the relationship between PTSD symptoms, “Flashback” and “Avoidance,” was the strongest. “Life emptiness” was the most central symptom with the highest EI in the PTSD and depression network model. Followed by “Fatigue” and “Interest loss.” Bridge symptoms connecting PTSD and depressive symptoms in our study were “Numb,” “High alertness,” “Sad mood,” and “Compunction and blame,” successively. The data-driven community detection suggested the differences in PTSD symptoms in the clustering process. The reliability of the network was approved by both stability and accuracy tests.

Conclusion: To the best of our knowledge, the current study first demonstrated the network structure of PTSD and depressive symptoms among Chinese firefighters, identifying the central and bridge symptoms. Targeting interventions to the symptoms mentioned above may effectively treat firefighters suffering from PTSD and depressive symptoms.

1. Introduction

Post-traumatic stress disorder (PTSD) is a chronic severe mental disorder after traumatic exposure, which is composed of four core symptoms, intrusion symptoms associated with the traumatic events, avoidance of traumatic stimuli, negative alterations in cognitions and mood, and marked alterations in arousal and reactivity (1). Given that over 70% of people may experience a traumatic event in their lifetime (2), previous studies indicate that people with PTSD symptoms are more likely to have other psychiatric symptom comorbidities (3), especially depressive symptoms (4). The high rate of symptom co-occurrence of PTSD and depression is demonstrated in populations with different categories (natural disasters, combat violence, the break of an intimate relationship, etc.) of traumatic exposures (5–8). Either PTSD symptoms or depressive symptoms are associated with higher final mortality and morbidity of physical symptoms [musculoskeletal pain, cardio-respiratory (CR) symptoms, gastrointestinal (GI) symptoms, etc.] (9, 10). Moreover, emerging research suggests that the co-occurrence of PTSD and MDD possibly has more severe adverse effects on mental health and life quality than either of them alone (11).

Firefighters, first responders with a high risk of occupational exposure to traumatic events and heavy working stress (12), are the population with a high prevalence of PTSD symptoms and depressive symptoms. A systematic review suggests that the mean prevalence of PTSD among firefighters is 12.3%, and the mean prevalence of depression is 18.7%. Both of them are higher than those of the community population (13). Since stressful and frightening working situations are almost inevitable for firefighters, the potential risk of a series of negative results caused by PTSD and depression is extremely serious, including alcohol abuse, occupation burnout, and even suicidal ideation or attempt (14–16). Besides, from the clinical perspective, there is a substantial amount of symptoms overlap between PTSD and depression, like sleep disturbance, inattention, avoidance, and withdrawal, etc. Diagnostic confusion brought by the similarity of symptoms is common, especially when clinicians lack information of traumatic event history (17). The concealment of psycho-trauma related information is not unusual in the Asian population, as high levels of stigma are general in the Asian culture and mental disorders are regarded as a sign of weakness (18). Therefore, we consider that it is meaningful to figure out the potential relationships and significance hierarchy among symptoms of PTSD and depression, which can provide references for accurate psychological intervention for firefighters.

In recent years, network analysis has become a widely applied method for exploring the associations among individual psychiatric symptoms and ascertaining their relative significance in psychopathology (19). Network analysis posits that psychiatric disorders/syndromes can be modeled as phenomena arising as a complex network of interacting and mutually reinforcing symptoms (20). Thus, specific symptoms act as positive roles in activating and maintaining mental disorders instead of passive indicators of mental disorders (21). Compared with the traditional statistical analysis method mainly investigate the unidirectional association (e.g., structural equation modeling, regression analysis, etc.) (22, 23), which demands a hypothesis of the potential relationships among symptoms, network analysis can identify the dynamic and reciprocal associations between various symptoms (24, 25).

Network analysis can provide the corresponding centrality and predictability index for each symptom node to test its significance and controllability in the whole network (26, 27). The centrality index can be used to ascertain central symptoms which contribute to the majority of the network stability, for further understanding mechanisms about the onset and maintenance of a disorder or syndrome and potential targets for clinical interventions (24). Moreover, when individuals suffer from different psychiatric disorders, some specific symptoms of one disorder may increase the risk of the other disorder, which is considered as bridge symptoms in the network. The bridge symptoms of the network play essential roles in maintaining and developing comorbidities, suggesting references for clinical prevention and treatments (28).

Previous studies of the network of the comorbidity of PTSD and depression were mainly conducted in veterans (29, 30), though heterogeneity exists between different studies, indicating higher centralities of symptoms “Flashback,” “Getting upset by trauma reminders,” “concentration difficulties,” and “anhedonia.” Another network analysis enrolling male participants with domestic violence indicates that symptoms “Feelings of worthlessness” and “Avoiding internal reminders of the traumatic experiences” are the most central symptoms (31). To the best of our knowledge, similar studies involving the comorbidity of PTSD and depression are rarely conducted in firefighter samples. Concerned about the difference in exposure to traumatic events and the potential effect of stigma rooted in Asian culture on mental symptoms, we speculated that features of the network of comorbidity of PTSD and depressive symptoms in Chinese firefighters might differ from those of other kind populations in the background of western culture.

As far as we know, the current study is the first study applying the network analysis method to investigate the network structure of comorbidities of PTSD symptoms and depressive symptoms among Chinese firefighters, which inspired us to conduct this study to fill the gap. The purpose of our study was to examine the associations between PTSD and depressive symptoms among Chinese firefighters, then to figure out the center symptoms and bridge symptoms of this PTSD-depressive symptoms network. Our results can provide theoretical references for establishing precise psychological interventions and normalized mental health care for firefighters.

2. Materials and methods

2.1. Participants and study design

Our study was a cross-sectional survey conducted in March 2021 in Changsha, Hunan Province, China. The Web-based questionnaire was applied in our research for higher data collection efficiency and lower risk of transmission of COVID-19. The questionnaire was set to be qualified for uploading when all items were filled correctly and checked by software, which could improve the reliability of raw data. After being approved by the administration of Hunan Fire Brigade, electronic questionnaires set by us were distributed by the management department of Hunan Fire Brigade. Inclusion criteria were (1) age above 18 years old, (2) certificated firefighters in fire stations, and (3) understanding the purpose and content of this research. Exclusion criteria were (1) not being at work due to any reasons, (2) not able to finish the whole questionnaire due to health reasons, and (3) without any experience of formal emergency mission. This study was approved by the Ethic Committee of the Second Xiangya Hospital of Central South University. All the participants have provided electronic written informed consent.

2.2. Measures

The primary care PTSD screen for DSM-5 (PC-PTSD-5) was used for assessing PTSD symptoms. PC-PTSD-5 is a five item self-report screening scale, which is set based on the diagnostic criteria of PTSD in DSM-5 (32). Items are scored dichotomously as either zero or one (0 = No; 1 = Yes). Higher scores mean more severe PTSD symptoms. The reliability of the Chinese version of PC-PTSD-5 in Chinese was validated and the cutoff score of potential PTSD was set as two points (33).

Depressive symptoms were assessed by the Self-Rating Depression Scale (SDS), which was compiled by Zung et al. (34). SDS is a classical and widely applied four-point Likert scale for depression in China (35), scored from one (a little of the time) to four (most of the time). Higher scores mean more severe depressive symptoms. According to the SDS Chinese manual (36), individuals with a total score of 40 or more were considered to have potential depression.

2.3. Statistical analysis

All the data were analyzed by the R program (37). The process of network analysis was divided into five domains, including network estimation and visualization, Exploratory community analysis, centrality and predictability analysis, network accuracy and stability, and network comparison.

2.3.1. Network estimation and visualization

Partial correlation analysis was applied to analyze the association between each pairwise nodes for controlling the confounding effects of other nodes in the network. The least absolute shrinkage and selection operator (LASSO) algorithm was used in the regularization process to shrink all edges in the network and set small correlations to zero, which enabled as few nodes and edges as possible in the network (38). Meanwhile, the extended Bayesian Information Criteria (EBIC) were also adopted to obtain a sparse and interpretable network model (39). The turning parameter was set for 0.5, which was widely applied for controlling spurious correlations in the network estimation (40).

Fruchterman-Reingold algorithm (41) was used to visualize the network. Nodes with stronger and more frequent associations with other nodes were placed closer together and more centralized in the network.

According to the parlance of network analysis, each item represented a symptom node, and each edge from one node to another was indicated as the association between two nodes. Thicker and more saturated edges meant stronger relationships. The color of edges represented the correlation, generally green for positive and red for negative (42). Concerned that the primary purpose of the current study was to investigate the inducing relationships between PTSD and depressive symptoms, and also for the simplicity and better understanding of the network structure, only positive relationships were depicted in the visualization of network. The possible negative associations among nodes were retained in the correlation matrix for the integrity of data. The R packages bootnet (Version 1.5) (43) and qgraph (Version 1.9.2) (42) were used to estimate and visualize the network.

2.3.2. Exploratory community analysis

Concerned that actual dimensions of symptoms in the PTSD and depression network might differed from initial hypothetical dimensions, we conducted the data-driven method, community detection, to identify the final dimensions of symptoms in the network. Meanwhile, the centralities (both original and bridge) of the final community detection grouping were calculated to be compared with the centralities of the original hypothetical grouping.

The Walktrap algorithm, based on the principle that adjacent nodes tend to belong to the same community (44), was applied to identify communities in the network model (45). This algorithm has been proven that performs well on psychological networks (46, 47). The function cluster_walktrap, available in the R package igraph was used to detect communities in the current study.

2.3.3. Centrality and predictability analysis

The centrality index, expected influence (EI), was calculated by the R package qgraph (1.9.2) (42) to quantify the importance of each symptom node in the network model. Symptom node EI refers to the sum of the value of all edges connecting to a given node, including both positive and negative values (48). Since there were not only positive associations but also negative associations in the network of our study, EI is more appropriate than other centrality indexes. Moreover, the bridge EI was analyzed to evaluate the significance of a node in connecting external symptom dimensions, by using the function bridge of the R package networktools (version 1.4.0) (49).

Predictability, an index suggesting the extent to which a node was predicted by its neighboring nodes (39), was computed by the function predict of the R package mgm (version 1.2–12) (50). Predictability represents the controllability of a node (27). In the picture of network model, the ring area around each node illustrated the value of predictability.

2.3.4. Network accuracy and stability

A total of three steps were conducted by the R package bootnet (version 1.5) (43) to assess the accuracy and stability of the network model in our research. Firstly, the accuracy of edge weights was estimated by constructing a 95% confidence interval (CI) with non-parametric bootstrapping method (51), which was visualized with a line diagram. A narrower bootstrapped CI of edge weights suggested a more precise estimation of the edges (52).

Secondly, the stability of symptom node EI was tested by calculating the correlation stability coefficient (CS-C), using a case-dropping bootstrap procedure (53). In this procedure, increasing percentage of cases was dropped from the dataset, and the EI indexes were re-estimated. If the EI of symptoms nodes did not change significantly after excluding a subset of the sample, the network structure could be considered stable. CS-C means the maximum proportion of samples that could be dropped while the EI correlation between the networks of original sample and case-dropping subsets was at least 0.7 with a 95% probability (43). In general, the value of CS-C needs to be above 0.25 and is preferably above 0.5 (43).

Thirdly, the significant differences between edge weights and node EI centralities were computed using the bootstrapped difference tests based on 95% CIs, which suggested that there were statistical differences between two edges weights or two EI centralities if zero was not included in the CIs (43).

2.3.5. Network comparison

Only-child status was collected as demographic information in the current study, according to which network comparison analysis was conducted to investigate the difference in network structure between only-child firefighters and non-only-child firefighters. China is a country with a huge population, family planning policy was conducted by the government so far, to control the rapid increase of population. Concerned that previous studies proved the effect of only-child status on mental health (54–56), we speculated that only-child and non-only-child firefighters might have a different network structure from each other.

The function network comparison test (NCT) of the R package NetworkComparisonTest (57) was applied to compare the network structures. The principle of NCT was a permutation-based test, which randomly regroups participants from each sub-network repeatedly (1,000 times) and then calculates the differences between networks. A total of three invariance measures were examined of sub-networks of genders, including global strength, network structure, and edge strength (57). The global strength tested the difference in overall network connectivity between sub-networks. The network structure examined the overall difference between all the possible edges between two sub-networks. Only when a statistically significant difference was found in either global strength or network structure comparisons, edge strength was further calculated to investigate the possible difference in each edge between the two networks. Moreover, the original centralities of nodes between sub-networks were also tested for the completeness of analysis.

3. Results

3.1. Descriptive statistics

A total of 1,781 firefighters were included in our study. For the demographic information, all of the participants were male, 624 (35.0%) of which were only-child and 1,157 were non-only-child. The average age of sample was 26.60 [standard deviation (SD) = 4.82]. The average total scores were 0.34 (SD = 0.83) for PC-PTSD-5 and 36.21 (SD = 9.34) for SDS. According to the cutoff scores applied in the current research, the proportions of firefighters with potential PTSD and potential depression were 8.7 and 36.3%, respectively. The mean score and abbreviation for each symptom node are shown in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Descriptive statistics of PTSD and depressive symptom nodes.

3.2. Network structure

The network of PTSD and depressive symptoms is depicted in Figure 1. The detailed Correlation matrix is represented in Supplementary Table 1.

FIGURE 1
www.frontiersin.org

Figure 1. Network structure of post-traumatic stress disorder (PTSD) and depressive symptoms in Chinese firefighters. Symptom nodes with stronger associations are placed closer to each other. The dark green lines represent positive correlations. The red lines represent negative correlations. The line thickness represents the strength of the connection between symptom nodes. Ring-shaped pie charts illustrate predictability (e.g., a fully filled dark ring means that 100% of the symptom’s variance can be explained by its intercorrelations with the other symptom nodes in the network).

A total of 300 edges [25 × (25–1)/2] were initially estimated, of which 140 edges were non-zero weights and entered further analysis. For the connection between symptoms, “Flashback” (PCPTSD1) and “Avoidance” (PCPTSD2) were the strongest edges (weight = 0.46), followed by the edges between “Hopeless” (SDS14) and “Life emptiness” (SDS18) (weight = 0.33), between “Unworthiness” (SDS17) and “Life emptiness” (SDS18) (weight = 0.30), and between “Confusion” (SDS11) and “Ability decline” (SDS12) (weight = 0.30).

In terms of the predictability, the node “Life emptiness” had the highest predictability value (71.1%), while the mean level of predictability of nodes was 38.3 ± 18.2%. These results meant this symptom node was the most central symptom in the PTSD-depression interactive network but also the node highest-predicted by surrounding nodes.

3.3. Community and centrality

The Walktrap algorithm detected three communities in the comorbidity of PTSD and depressive symptoms (Figure 2). “Flashback” (PCPTSD1) and “Avoidance” (PCPTSD2) remained in the community of PTSD (Community 2), while “High alertness” (PCPTSD3), “Numb feeling” (PCPTSD4), and “Compunction and blame” (PCPTSD5) were included in the Community 1 with depressive symptoms. The Community 3 was entirely generated by the rest depressive symptoms.

FIGURE 2
www.frontiersin.org

Figure 2. Communities detected by the Walktrap algorithm in the network of PTSD and depressive symptoms.

All centralities (original EIs and bridge EIs) based on the theoretical model and data-driven detected model were both calculated for comparison. Details are shown in Table 2. As the original EI, measuring the significance of a certain node in the whole network, was not affected by the dimension distribution. Thus, original EIs based on the theoretical model were as same as the data-driven detected model, which are shown in Figure 3. The depressive symptom, “Life emptiness” (SAS18), had the highest EI value, indicating that this symptom was the most central symptom in the network model, which was followed by the symptom “Fatigue” (SDS 10), “Interest loss” (SDS 20), and “Hopeless” (SDS 14).

TABLE 2
www.frontiersin.org

Table 2. Original centralities and bridge centralities of PTSD and depressive nodes in the network.

FIGURE 3
www.frontiersin.org

Figure 3. Expected influences of PTSD and depressive symptom nodes in the network (shown as standardized values z scores).

Regarding the bridge EI, in the theoretical model, “Numb feeling” (PC-PTSD-4) had the highest bridge EI, followed by “High alertness” (PC-PTSD-3), “Sad mood” (SDS1), and “Compunction and blame” (PC-PTSD-5). Details are shown in Figure 4.

FIGURE 4
www.frontiersin.org

Figure 4. Bridge expected influences of PTSD and depressive symptom nodes in the theoretical model (shown as standardized values z scores).

In the data-driven model (detected by the Walktrap algorithm), the bridge EIs of symptom nodes were calculated for each community (Figure 5). In community 1, “Mood circadian rhythm” (SDS2) was the node with the highest bridge EI, being significant in connecting external nodes. In community 2, “Flashback” (PCPTSD1) was the node with the highest bridge EI. In community 3, “Life emptiness” (SAS18) had the highest bridge EI.

FIGURE 5
www.frontiersin.org

Figure 5. Bridge expected influences of PTSD and depressive symptom nodes in the data-driven model (shown as standardized values z scores).

3.4. Network accuracy and stability

For the accuracy of edge weight, as shown in Supplementary Figure 1, bootstrap 95% CIs were narrow, indicating that the estimations of edges were reliable. As presented in Figure 6, the case-dropping bootstrap procedure indicates that the CS-C of EI was 0.75, reflecting that the network kept being stable after dropping 75% of the sample. In the Bootstrapped difference tests for the node EIs (Figure 7), central symptom nodes differed from most other nodes. And the Bootstrapped difference tests for edge weights (Supplementary Figure 2) also showed that most comparisons between edges were statistically different, especially the strongest edges.

FIGURE 6
www.frontiersin.org

Figure 6. Stability of expected influence using case-dropping bootstrap. The x-axis represents the percentage of cases in the original sample used at each step. The y-axis represents the average of correlations between the expected influence in the original network and the expected influence in the networks that were re-estimated after dropping increasing percentages of cases.

FIGURE 7
www.frontiersin.org

Figure 7. Estimation of node expected influence difference by bootstrapped difference test. Gray boxes indicate nodes that do not significantly differ from one another. Black boxes represent nodes that differ significantly from one another. Values in box on the diagonal represent the expected influence of each node.

3.5. Network comparison

The network structures for both only-child firefighter sub-network and non-only child firefighter sub-network are, respectively, shown in Supplementary Figures 3, 4. The network comparison test (NCT) did not find statistically significant difference in the global strength invariance between only-child firefighters and non-only-child firefighters (only-child: 10.497 vs. non-only-child: 9.688, Test statistic S = 0.809, p = 0.125). Meanwhile, there was no significant difference found in network structure invariance (Test statistic M = 0.175, p = 0.311). No significant difference of original centrality was detected in the comparison of only-child sub-network and non-only child sub-network (Test statistic C = 0.133–0.468, p = 0.20–0.99). Details of original centrality comparison are shown in Supplementary Table 2.

4. Discussion

To the best of our knowledge, the current study is the first study focusing on characterizing the network structure of PTSD and depressive symptoms among the firefighters in China. For the associations between symptom nodes, there were no cross-dimension edges in the strong edges. The strongest edge was the connection between “Flashback” and “Avoidance” within the PTSD sub-network.

Compared with a previous study conducted in veterans in which the edge between PTSD and depressive symptoms was the strongest in the whole network (29), our results suggested that the association between “Flashback” and “Avoidance” is the strongest connection, even in the symptom network simultaneously including PTSD and depressive symptoms. We speculated that our results revealed the effect of PTSD symptoms in the conjunct network was more limited in firefighters. Thus, the risk of inducing depressive symptoms caused by traumatic-related symptoms was lower.

According to the node EI and predictability results, “Life emptiness” was the most central symptom with the highest EI and the symptom node with the highest predictability. Feeling of emptiness, which could be described as a sense of deadness or absence of inner feelings, is a complexed and negative emotion, including physical component, aloneness component and component involving personal unfulfillment, and lack of purpose (58, 59). Feeling of emptiness was observed in people with depression (60) and PTSD (61), thus, it should be recognized as a transdiagnostic construct. Additionally, previous studies demonstrated that the comorbidity of PTSD and depression could be explained by similar psychopathology symptom dimensions across these two illnesses. And PTSD and depression might present a same traumatic stress construct with common related factors, especially in chronic cases (4). Based previous results, the current study further demonstrated that “feeling of emptiness” was the most significant factor in activating and stabilizing the PTSD-depression psychopathology network. Concerned the high incidence of PTSD and depression comorbidity (17, 62, 63), it is meaningful to figure out the core symptom node in the PTSD-depression network and manipulate accurate intervention aiming “feeling of emptiness” could simultaneously alleviate both PTSD and depressive symptoms.

As for the bridge symptoms, from the perspective of the association between PTSD symptoms and depressive symptoms (i.e., the theoretical model), the most influential bridge symptom node was the “Numb feeling.” Numb feeling was considered as one of the hallmark symptoms of PTSD, including markedly diminished interest in significant activities (Criterion C-4), feelings of detachment or estrangement from others (C-5), and restricted range of affect (C-6), which were represented in the diagnostic criteria of PTSD of DSM-IV (64). Meanwhile, anhedonia is a symptomatic deficit in positive affect, involving the loss of enjoyment in pleasurable activities or the loss of desire to engage in them. This is similar to the Criterion C-4 criterion of numb feeling. Anhedonia was considered to be related to the core symptom of depression and was so common that approximately one-third of depressed individuals have clinically significant anhedonia symptoms (65). Based on the previous studies and the results of our study, we speculated that our results demonstrated the “Numb feeling” acted as “mediating node” between the PTSD symptom cluster and the depressive symptom cluster in the Chinese firefighters, a population with a high risk of traumatic exposure. Additionally, an emotional dysregulation model of PTSD illustrated by Litz et al. (66) held that the “emotionally numb” of individuals with PTSD was actually a result of hyperresponsivity to negative emotional stimulation, individuals with PTSD required more intense positive stimulation to access pleasure. Thus, our results verified this theoretical model among Chinese firefighters from the network analysis perspective. Concerned depressive symptoms usually followed or co-occurred with the PTSD symptoms, instead of existing alone in the context of the traumatic exposure, psychological intervention targeting at “Numb feeling” might control the comorbidity of PTSD and depression.

Notably, the bridge symptoms in the network were also supported by the results of data-driven model. “High alertness,” “Numb feeling,” and “Compunction and blame,” three PTSD symptom nodes with relatively high bridge EI were divided into a community with some depressive symptom nodes. While “Flashback” and “Avoidance” composed a community without depressive nodes. Our results suggested that among PTSD symptom nodes, “Flashback” and “Avoidance” were the two traumatic-event related but less associated with depressive symptoms, while “High alertness,” “Numb feeling,” and “Compunction and blame” were highly related to the depression. The results of community detection could act as a theoretical reference for the clinical diagnosis. Concerning the similarity of clinical symptoms between PTSD and depression (1), if firefighters were only observed with PTSD symptoms mentioned above which were closely related to depressive symptoms, the possibility of single depression diagnosis, instead of PTSD, should not be overlooked.

Additionally, “Mood circadian rhythm” and “Life emptiness” were the depressive nodes with the highest bridge EI in the respective community. As “Life emptiness” was also the node with the highest original EI, the data-driven results again proved its significance in the PTSD and depressive symptom network. Whereas, “Mood circadian rhythm” had a relatively low original EI but the highest bridge EI. Combined with the location of “Mood circadian rhythm” node in the total network (mainly connecting Community 2 and Community 3), we considered that this situation reflected the circadian rhythm of a series of depressive symptoms in Community 3. Given that the node with high bridge EI play a critical role in inducing the external community symptoms, and that Chinese firefighters are a high-risk, high-stress occupation with a lack of work routines, individualized psychological interventions based on the circadian rhythms of depressive symptoms are essential. Meanwhile, ensuring their healthy life routines can promote the recovery of depressive symptoms.

5. Limitations

Notwithstanding our study firstly demonstrated the PTSD and depression symptom network structure of Chinese firefighters, there are still some limitations need to be clarified. First, the current study was cross-sectional, thus the dynamic alterations and causality among symptom nodes could not be investigated. Second, firefighters recruited in our study were all male, which might limit the promotion of our results among female firefighters. As Cao et al. (67) demonstrated that the sex difference in symptoms connectivity of the PTSD symptom network among adolescents, we speculated that our results need to be further verified in female firefighters. Third, the symptom networks were generated specific to self-reported assessments applied in our study; therefore, the network structure was possibly affected by recall bias and participants’ social desirability. Additionally, self-reported assessments were applied rather than clinical interview; thus, atypical clinical features of PTSD and/or depression might be ignored.

6. Conclusion

To the best of our knowledge, the current study firstly demonstrated the network structure of PTSD and depressive symptoms among Chinese firefighters, identifying the central symptoms (i.e., Life emptiness, Fatigue, and Hopeless) and bridge symptoms (i.e., Numb feeling, Sad mood, and Compunction and blame). Meanwhile, the association between two clusters of symptoms was further clarified by the community detection analysis, indicating that core PTSD symptoms were “Flashback” and “Avoidance,” whereas “High alertness,” “Numb feeling,” and “Compunction and blame” were more relevant to depressive symptoms. Our work could act as an impetus for future studies investigating symptom networks of firefighters or other populations with high occupational risk (e.g., policemen, soldiers, etc.) as a method to verify the core symptoms that generalize across various populations with different culture backgrounds and those which are specific for certain populations. Moreover, central and bridge symptoms demonstrated by the current study can be possible targets for clinical intervention to treat firefighters who are suffering from PTSD and depressive symptoms.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethic Committee of the Second Xiangya Hospital of Central South University. The patients/participants provided their written informed consent to participate in this study.

Author contributions

PC and LW: literature review and manuscript drafting. YZ, WM, and GZ: data acquisition. LZ and WL: study design and manuscript revision. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by Natural Science Foundation of China (No. 82171518 to LZ), the Natural Science Foundation of Hunan Province, China (No. 2020JJ5844 to LZ), and Technological Innovation Foundation of Central South University (No. 502801002 to WL).

Acknowledgments

We sincerely thank all participants for their time and patience in this study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2023.1096771/full#supplementary-material

References

1. American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders (DSM-5) Arlington, VA: American Psychiatric Pub (2013).

Google Scholar

2. Benjet, C, Bromet, E, Karam, EG, Kessler, RC, McLaughlin, KA, Ruscio, AM, et al. The epidemiology of traumatic event exposure worldwide: results from the world mental health survey consortium. Psychol Med. (2016) 46:327–43. doi: 10.1017/s0033291715001981

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Sareen, J, Cox, BJ, Stein, MB, Afifi, TO, Fleet, C, and Asmundson, GJ. Physical and mental comorbidity, disability, and suicidal behavior associated with posttraumatic stress disorder in a large community sample. Psychosom Med. (2007) 69:242–8. doi: 10.1097/PSY.0b013e31803146d8

PubMed Abstract | CrossRef Full Text | Google Scholar

4. O'Donnell, ML, Creamer, M, and Pattison, P. Posttraumatic stress disorder and depression following trauma: understanding comorbidity. Am J Psychiatry. (2004) 161:1390–6. doi: 10.1176/appi.ajp.161.8.1390

CrossRef Full Text | Google Scholar

5. Hankin, CS, Spiro, A 3rd, Miller, DR, and Kazis, L. Mental disorders and mental health treatment among U.S. Department of Veterans Affairs outpatients: the veterans health study. Am J Psychiatry. (1999) 156:1924–30. doi: 10.1176/ajp.156.12.1924

CrossRef Full Text | Google Scholar

6. Hsu, CC, Chong, MY, Yang, P, and Yen, CF. Posttraumatic stress disorder among adolescent earthquake victims in Taiwan. J Am Acad Child Adolesc Psychiatry. (2002) 41:875–81. doi: 10.1097/00004583-200207000-00022

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Nixon, RD, Resick, PA, and Nishith, P. An exploration of comorbid depression among female victims of intimate partner violence with posttraumatic stress disorder. J Affect Disord. (2004) 82:315–20. doi: 10.1016/j.jad.2004.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Norris, FH, Murphy, AD, Baker, CK, and Perilla, JL. Postdisaster PTSD over four waves of a panel study of Mexico's 1999 flood. J Trauma Stress. (2004) 17:283–92. doi: 10.1023/B:JOTS.0000038476.87634.9b

CrossRef Full Text | Google Scholar

9. Lépine, JP, and Briley, M. The increasing burden of depression. Neuropsychiatr Dis Treat. (2011) 7:3–7. doi: 10.2147/ndt.S19617

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Pacella, ML, Hruska, B, and Delahanty, DL. The physical health consequences of PTSD and PTSD symptoms: a meta-analytic review. J Anxiety Disord. (2013) 27:33–46. doi: 10.1016/j.janxdis.2012.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Flory, JD, and Yehuda, R. Comorbidity between post-traumatic stress disorder and major depressive disorder: alternative explanations and treatment considerations. Dialogues Clin Neurosci. (2015) 17:141–50. doi: 10.31887/DCNS.2015.17.2/jflory

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Van Eerd, D, Irvin, E, Harbin, S, Mahood, Q, and Tiong, M. Occupational exposure and post-traumatic stress disorder: a rapid review. Work. (2021) 68:721–31. doi: 10.3233/wor-203406

CrossRef Full Text | Google Scholar

13. Wagner, SL, White, N, Randall, C, Regehr, C, White, M, Alden, LE, et al. Mental disorders in firefighters following large-scale disaster. Disaster Med Public Health Prep. (2021) 15:504–17. doi: 10.1017/dmp.2020.61

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Gulliver, SB, Zimering, RT, Dobani, F, Pennington, ML, Morissette, SB, Kamholz, BW, et al. Alcohol use and mental health symptoms in female firefighter recruits. Occup Med (Lond). (2019) 69:625–31. doi: 10.1093/occmed/kqaa015

CrossRef Full Text | Google Scholar

15. Katsavouni, F, Bebetsos, E, Malliou, P, and Beneka, A. The relationship between burnout, PTSD symptoms and injuries in firefighters. Occup Med (Lond). (2016) 66:32–7. doi: 10.1093/occmed/kqv144

CrossRef Full Text | Google Scholar

16. Martin, CE, Vujanovic, AA, Paulus, DJ, Bartlett, B, Gallagher, MW, and Tran, JK. Alcohol use and suicidality in firefighters: associations with depressive symptoms and posttraumatic stress. Compr Psychiatry. (2017) 74:44–52. doi: 10.1016/j.comppsych.2017.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Brady, KT, Killeen, TK, Brewerton, T, and Lucerini, S. Comorbidity of psychiatric disorders and posttraumatic stress disorder. J Clin Psychiatry. (2000) 61:22–32.

Google Scholar

18. Karasz, A. Cultural differences in conceptual models of depression. Soc Sci Med. (2005) 60:1625–35. doi: 10.1016/j.socscimed.2004.08.011

CrossRef Full Text | Google Scholar

19. Borsboom, D, and Cramer, AO. Network analysis: an integrative approach to the structure of psychopathology. Annu Rev Clin Psychol. (2013) 9:91–121. doi: 10.1146/annurev-clinpsy-050212-185608

CrossRef Full Text | Google Scholar

20. Borsboom, D. A network theory of mental disorders. World Psychiatry. (2017) 16:5–13. doi: 10.1002/wps.20375

PubMed Abstract | CrossRef Full Text | Google Scholar

21. McNally, RJ. Can network analysis transform psychopathology? Behav Res Ther. (2016) 86:95–104. doi: 10.1016/j.brat.2016.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Carrà, G, Crocamo, C, Angermeyer, M, Brugha, T, Toumi, M, and Bebbington, P. Positive and negative symptoms in schizophrenia: a longitudinal analysis using latent variable structural equation modelling. Schizophr Res. (2019) 204:58–64. doi: 10.1016/j.schres.2018.08.018

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Fowler, D, Hodgekins, J, Garety, P, Freeman, D, Kuipers, E, Dunn, G, et al. Negative cognition, depressed mood, and paranoia: a longitudinal pathway analysis using structural equation modeling. Schizophr Bull. (2012) 38:1063–73. doi: 10.1093/schbul/sbr019

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Belvederi Murri, M, Amore, M, Respino, M, and Alexopoulos, GS. The symptom network structure of depressive symptoms in late-life: results from a European population study. Mol Psychiatry. (2020) 25:1447–56. doi: 10.1038/s41380-018-0232-0

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Herniman, SE, Phillips, LJ, Wood, SJ, Cotton, SM, Liemburg, EJ, and Allott, KA. Interrelationships between depressive symptoms and positive and negative symptoms of recent onset schizophrenia spectrum disorders: a network analytical approach. J Psychiatr Res. (2021) 140:373–80. doi: 10.1016/j.jpsychires.2021.05.038

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Contreras, A, Nieto, I, Valiente, C, Espinosa, R, and Vazquez, C. The study of psychopathology from the network analysis perspective: a systematic review. Psychother Psychosom. (2019) 88:71–83. doi: 10.1159/000497425

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Haslbeck, JMB, and Fried, EI. How predictable are symptoms in psychopathological networks? A reanalysis of 18 published datasets. Psychol Med. (2017) 47:2767–76. doi: 10.1017/s0033291717001258

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Jones, PJ, Ma, R, and McNally, RJ. Bridge centrality: a network approach to understanding comorbidity. Multivar Behav Res. (2021) 56:353–67. doi: 10.1080/00273171.2019.1614898

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Duek, O, Spiller, TR, Pietrzak, RH, Fried, EI, and Harpaz-Rotem, I. Network analysis of PTSD and depressive symptoms in 158,139 treatment-seeking veterans with PTSD. Depress Anxiety. (2021) 38:554–62. doi: 10.1002/da.23112

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Lazarov, A, Suarez-Jimenez, B, Levi, O, Coppersmith, DDL, Lubin, G, Pine, DS, et al. Symptom structure of PTSD and co-morbid depressive symptoms - a network analysis of combat veteran patients. Psychol Med. (2020) 50:2154–70. doi: 10.1017/s0033291719002034

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Gilbar, O. Examining the boundaries between ICD-11 PTSD/CPTSD and depression and anxiety symptoms: a network analysis perspective. J Affect Disord. (2020) 262:429–39. doi: 10.1016/j.jad.2019.11.060

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Prins, A, Bovin, MJ, Smolenski, DJ, Marx, BP, Kimerling, R, Jenkins-Guarnieri, MA, et al. The primary care PTSD screen for DSM-5 (PC-PTSD-5): development and evaluation within a veteran primary care sample. J Gen Intern Med. (2016) 31:1206–11. doi: 10.1007/s11606-016-3703-5

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Cheng, P, Jasinski, N, Zheng, W, Yadava, A, Wang, L, Li, L, et al. Psychometric properties of the primary care PTSD screen for DSM-5: findings from family members of Chinese healthcare workers during the outbreak of COVID-19. Front Psychol. (2021) 12:695678. doi: 10.3389/fpsyt.2021.695678

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Zung, WW. A Self-Rating Depression Scale. Arch Gen Psychiatry. (1965) 12:63–70. doi: 10.1001/archpsyc.1965.01720310065008

CrossRef Full Text | Google Scholar

35. Wang, X, Wang, X, and Ma, H. Manual of mental health rating scale. Chin Ment Health J. (1999) 13:31–5.

Google Scholar

36. Guo, NF. China National Vocational Qualification Training Course-Psychological Counselor. Beijing, China: Nationalities Press (2012).

Google Scholar

37. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing (2022).

Google Scholar

38. Friedman, J, Hastie, T, and Tibshirani, R. Sparse inverse covariance estimation with the graphical lasso. Biostatistics. (2008) 9:432–41. doi: 10.1093/biostatistics/kxm045

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Haslbeck, JMB, and Waldorp, LJ. How well do network models predict observations? On the importance of predictability in network models. Behav Res Methods. (2018) 50:853–61. doi: 10.3758/s13428-017-0910-x

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Isvoranu, AM, van Borkulo, CD, Boyette, LL, Wigman, JT, Vinkers, CH, and Borsboom, D. A network approach to psychosis: pathways between childhood trauma and psychotic symptoms. Schizophr Bull. (2017) 43:187–96. doi: 10.1093/schbul/sbw055

CrossRef Full Text | Google Scholar

41. Fruchterman, TMJ, and Reingold, EM. Graph drawing by force-directed placement. Softw Pract Exp. (1991) 21:1129–64. doi: 10.1002/spe.4380211102

CrossRef Full Text | Google Scholar

42. Epskamp, S, Cramer, AOJ, Waldorp, LJ, Schmittmann, VD, and Borsboom, D. Qgraph: network visualizations of relationships in psychometric data. J Stat Softw. (2012) 48:1–8. doi: 10.18637/jss.v048.i04

CrossRef Full Text | Google Scholar

43. Epskamp, S, Borsboom, D, and Fried, EI. Estimating psychological networks and their accuracy: a tutorial paper. Behav Res Methods. (2018) 50:195–212. doi: 10.3758/s13428-017-0862-1

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Yang, Z, Algesheimer, R, and Tessone, CJ. A comparative analysis of community detection algorithms on artificial networks. Sci Rep. (2016) 6:30750. doi: 10.1038/srep30750

CrossRef Full Text | Google Scholar

45. Briganti, G, and Linkowski, P. Item and domain network structures of the resilience scale for adults in 675 university students. Epidemiol Psychiatr Sci. (2019) 29:e33. doi: 10.1017/s2045796019000222

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Gates, KM, Henry, T, Steinley, D, and Fair, DA. A Monte Carlo evaluation of weighted community detection algorithms. Front Neuroinform. (2016) 10:45. doi: 10.3389/fninf.2016.00045

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Golino, HF, and Epskamp, S. Exploratory graph analysis: a new approach for estimating the number of dimensions in psychological research. PLoS One. (2017) 12:e0174035. doi: 10.1371/journal.pone.0174035

CrossRef Full Text | Google Scholar

48. Robinaugh, DJ, Millner, AJ, and McNally, RJ. Identifying highly influential nodes in the complicated grief network. J Abnorm Psychol. (2016) 125:747–57. doi: 10.1037/abn0000181

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Payton, J. (2021). "Networktools: Tools for identifying important nodes in networks".

Google Scholar

50. Haslbeck, J, and Waldorp, L. Mgm: estimating time-varying mixed graphical models in high-dimensional data. J Stat Softw. (2020) 93:1–46. doi: 10.18637/jss.v093.i08

CrossRef Full Text | Google Scholar

51. Sherman, M. Book review of: M. R. Chernick, bootstrap methods: a guide for practitioners and researchers. J Am Stat Assoc. (2008) 103:1705–28. doi: 10.1198/016214508000000887

CrossRef Full Text | Google Scholar

52. Marchetti, I. Hopelessness: a network analysis. Cogn Ther Res. (2019) 43:611–9. doi: 10.1007/s10608-018-9981-y

CrossRef Full Text | Google Scholar

53. Costenbader, E, and Valente, TW. The stability of centrality measures when networks are sampled. Soc Networks. (2003) 25:283–307. doi: 10.1016/S0378-8733(03)00012-1

CrossRef Full Text | Google Scholar

54. Chow, E, and Zhao, S. The one-child policy and parent-child relationships: a comparison of one-child with multiple-child families in China. Int J Sociol Soc Policy. (1996) 16:35–62. doi: 10.1108/eb013285

CrossRef Full Text | Google Scholar

55. Mancillas, A. Challenging the stereotypes about only children: a review of the literature and implications for practice. J Couns Dev. (2006) 84:268–75. doi: 10.1002/j.1556-6678.2006.tb00405.x

CrossRef Full Text | Google Scholar

56. Yang, J, Hou, X, Wei, D, Wang, K, Li, Y, and Qiu, J. Only-child and non-only-child exhibit differences in creativity and agreeableness: evidence from behavioral and anatomical structural studies. Brain Imag Behav. (2017) 11:493–502. doi: 10.1007/s11682-016-9530-9

PubMed Abstract | CrossRef Full Text | Google Scholar

57. van Borkulo, C, van Bork, R, Boschloo, L, Kossakowski, J, Tio, P, Schoevers, R, et al. Comparing network structures on three aspects: a permutation test. Psychol Methods. (2021). doi: 10.1037/met0000476

CrossRef Full Text | Google Scholar

58. D'Agostino, A, Pepi, R, Rossi Monti, M, and Starcevic, V. The feeling of emptiness: a review of a complex subjective experience. Harv Rev Psychiatry. (2020) 28:287–95. doi: 10.1097/hrp.0000000000000269

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Zandersen, M, and Parnas, J. Identity disturbance, feelings of emptiness, and the boundaries of the schizophrenia Spectrum. Schizophr Bull. (2019) 45:106–13. doi: 10.1093/schbul/sbx183

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Rhodes, J, Hackney, S, and Smith, J. Emptiness, engulfment, and life struggle: an interpretative phenomenological analysis of chronic depression. J Constr Psychol. (2018) 32:390–407. doi: 10.1080/10720537.2018.1515046

CrossRef Full Text | Google Scholar

61. Darves-Bornoz, JM, Lépine, JP, Choquet, M, Berger, C, Degiovanni, A, and Gaillard, P. Predictive factors of chronic post-traumatic stress disorder in rape victims. Eur Psychiatry. (1998) 13:281–7. doi: 10.1016/s0924-9338(98)80045-x

CrossRef Full Text | Google Scholar

62. Creamer, M, Burgess, P, and McFarlane, AC. Post-traumatic stress disorder: findings from the Australian National Survey of mental health and well-being. Psychol Med. (2001) 31:1237–47. doi: 10.1017/s0033291701004287

CrossRef Full Text | Google Scholar

63. Kessler, RC, Sonnega, A, Bromet, E, Hughes, M, and Nelson, CB. Posttraumatic stress disorder in the National Comorbidity Survey. Arch Gen Psychiatry. (1995) 52:1048–60. doi: 10.1001/archpsyc.1995.03950240066012

CrossRef Full Text | Google Scholar

64. American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders. 4th ed. Washington DC: American Psychiatric Association (1994).

Google Scholar

65. Pelizza, L, and Ferrari, A. Anhedonia in schizophrenia and major depression: state or trait? Ann General Psychiatry. (2009) 8:22. doi: 10.1186/1744-859x-8-22

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Litz, BT, and Gray, MJ. Emotional numbing in posttraumatic stress disorder: current and future research directions. Aust N Z J Psychiatry. (2002) 36:198–204. doi: 10.1046/j.1440-1614.2002.01002.x

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Cao, X, Wang, L, Cao, C, Fang, R, Chen, C, Hall, BJ, et al. Sex differences in global and local connectivity of adolescent posttraumatic stress disorder symptoms. J Child Psychol Psychiatry. (2019) 60:216–24. doi: 10.1111/jcpp.12963

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: post-traumatic stress disorder, depression, network analysis, firefighter, mental health

Citation: Cheng P, Wang L, Zhou Y, Ma W, Zhao G, Zhang L and Li W (2023) Post-traumatic stress disorder and depressive symptoms among firefighters: a network analysis. Front. Public Health 11:1096771. doi: 10.3389/fpubh.2023.1096771

Received: 12 November 2022; Accepted: 04 April 2023;
Published: 04 May 2023.

Edited by:

Changiz Mohiyeddini, Oakland University William Beaumont School of Medicine, United States

Reviewed by:

Cristian Ramos-Vera, Cesar Vallejo University, Peru
Shengdong Chen, Qufu Normal University, China

Copyright © 2023 Cheng, Wang, Zhou, Ma, Zhao, Zhang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Weihui Li, weihui_li@csu.edu.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.