Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Psychol., 11 August 2025

Sec. Health Psychology

Volume 16 - 2025 | https://doi.org/10.3389/fpsyg.2025.1572011

This article is part of the Research TopicBridging the Gap: Physical Manifestations and Mental Health in AdolescentsView all 8 articles

Exploring the co-occurrence of depression and anxiety symptoms among Chinese vocational high school adolescents through a network approach

  • 1School of Psychology, Shenzhen University, Shenzhen, China
  • 2The Shenzhen Humanities and Social Sciences Key Research Bases of the Center for Mental Health, Shenzhen, China
  • 3Department of Psychology, School of Sociology and Psychology, Central University of Finance and Economics, Beijing, China
  • 4Department of Psychology, Faculty of Social Sciences, University of Macau, Macau SAR, China
  • 5Experimental Middle School to Sichuan Institute of Education Sciences, Chengdu, Sichuan, China

Depression and anxiety, prevalent among adolescents, show higher rates in females than males. However, little is known about symptom-level interactions and gender-specific network characteristics in secondary vocational students. This study adopted network analysis, including Gaussian graph models and directed acyclic graphs, to explore symptom-level relationships between depression and anxiety in secondary vocational students, identifying bridge and core symptoms, potential causation, and gender-specific network characteristics. Involving 7,471 students in a cross-sectional study, we utilized the Patient Health Questionnaire 8-item (PHQ-8) and Generalized Anxiety Disorders 7-item (GAD-7) for symptom measurement. Key findings include: bridge symptom was “sad mood (PHQ2).” Core symptoms were “fatigue (PHQ4),” “loss of control (GAD2),” “trouble relaxing (GAD4),” “sad mood (PHQ2),” “irritability (GAD6),” “excessive worry (GAD3),” and “feeling guilty (PHQ6).” There were numerous cross-disorder links between anxiety and depressive symptoms, predominantly from anxiety symptoms to depressive symptoms. “Irritability (GAD6)” was the initial symptom in the directed network. Females showed higher prevalence of depression and anxiety, and network global strength. This study firstly highlights the complex interplay between depressive and anxiety symptoms in secondary vocational students, with notable gender differences. The insights offer valuable guidance for developing targeted interventions and support strategies for this group.

1 Introduction

The largest vocational education system in the world has been built in China, and it starts in vocational high school education. After completing 9 years of compulsory education, students with better academic performance usually attend regular high schools, while those with lower performance attend vocational high schools (Liu et al., 2022). Of note, almost equal numbers of adolescents enter the two education systems. Despite being at the same stage development, vocational high school students often face stigmatization and derogatory labeling regarding their personality traits, behavior patterns, and lower social status (Spruyt et al., 2015). While academic performance differences may partially explain educational tracking, pre-existing mental health challenges could also contribute to lower academic achievement, complicating causal interpretations (Coledam et al., 2022).

Notably, these rates in vocational high school students reached 46.5% for depression and 58.7% for anxiety (Yu et al., 2017), higher than the prevalence of depression and anxiety among regular high school students stands at 28.0% and 26.3% (Yu et al., 2022). Depression and anxiety are not only prevalent but also commonly co-occur, presenting a significant challenge in adolescents’ mental health (Garber and Weersing, 2010). When suffering from either depression and anxiety, the incidence of comorbidity is 10%–75% (Cummings et al., 2014). This comorbidity complicates treatment (Huang et al., 2019), extends disease duration and recovery period (Hung et al., 2020), and leads to more severe outcomes, including higher rates of disability (Zhou et al., 2017). Meanwhile, the relationship between depression and anxiety is intricate and bidirectional, with each potentially predicting the other (Jacobson and Newman, 2017). Therefore, beyond concentrating solely on depression or anxiety, it is important to consider the prevention and intervention strategies for their comorbidity.

A promising framework for investigating depression, anxiety, and their comorbidity is the network approach, which has been widely applied in recent years. Network analysis conceptualizes psychiatric disorders as emerging from dynamic interactions among their constituent symptoms. Within this framework, individual symptoms are represented as nodes in a network, while their statistical associations are depicted as edges (Borsboom and Cramer, 2013). In the current study, this methodology revealed critical structural features of symptom co-occurrence, including the identification of bridge symptoms (e.g., sad mood) that connect distinct symptom clusters, and core symptoms (e.g., fatigue and irritability) that exhibit high centrality within the network. These findings underscore the utility of network analysis in delineating symptom-level mechanisms underlying the comorbidity of depression and anxiety. Network analysis is a promising approach to Fig out important symptoms for targeted interventions. In the context of comorbidity, bridge symptoms and core symptoms are considered important. Symptoms that link disorders are termed “bridge symptoms,” identified by bridge centrality measures, play a critical role in the overlap between disorders (Jones et al., 2021). This insight leads to more targeted and effective interventions, as it highlights pivotal symptoms that drive comorbidity between disorders (Jones et al., 2021). Crucially, bridge symptoms differ across populations and clinical contexts. For instance, in a Gaussian graphical model (GGM) analysis of Chinese adolescents, guilt (PHQ6) was identified as a bridge symptom specifically linking depressive symptoms (e.g., fatigue, sad mood) to anxiety symptoms (e.g., excessive worry, irritability) (Cai et al., 2022). In contrast, anhedonia (PHQ1) emerged as a bridge symptom connecting depressive and generalized anxiety disorder symptoms in Bayesian network analyses of adult populations (Bai et al., 2021; Garabiles et al., 2019). Regarding core symptoms—those with high connectivity to multiple nodes within the network—prior research utilizing undirected network models has demonstrated their utility as ideal therapeutic targets. For example, interventions targeting highly central symptoms (e.g., insomnia, psychomotor agitation) have shown cascading benefits in reducing comorbid symptom severity (Borsboom et al., 2011). However, these studies relied largely on Gaussian Graphical Model (GGM), which lacks the ability to establish causality (Briganti et al., 2023). It remains unclear whether these important symptoms cause the disorder or are a consequence of it.

Bayesian networks structured as directed acyclic graph (DAG) can be used to estimate the activation sequence of symptoms in the network. DAG structures the relationships between variables within a probabilistic framework, reflecting the principles of causal reasoning (Pearl and Mackenzie, 2018; Pearl, 2009). Within a DAG, the joint probability distribution defines the intricate interplay among variables, with each arrow symbolizing a potential causal link. In psychopathology, DAG can elucidate potential causal relationships between symptoms, where each directed edge implies causation from one symptom to another, offering a directional perspective that undirected networks lack (Briganti et al., 2023). For example, DAG analysis has shed light on the causal role of re-experiencing symptoms in the comorbidity of depression and PTSD (Lazarov et al., 2020). However, it appears that this approach has not yet been applied to study the comorbidity of depression and anxiety.

While network analysis has been increasingly applied to identify symptom-level intervention targets in adolescent populations (Fang et al., 2024), no studies to date have employed this methodology to investigate depression and anxiety comorbidity among Chinese vocational high school students. Prior studies have typically employed undirected network models, which provide insights into the associative patterns of symptoms but do not address causality. Moreover, while research has indicated that bridge symptoms and core symptoms differ across populations, there is a lack of study specifically concerning the networks of depression and anxiety among vocational high school students, whether through undirected or directed approach. This study aims to bridge these gaps, This study identifies critical symptoms by employing Gaussian Graphical Models (GGM) and directed acyclic graphs (DAG) to elucidate undirected and directed symptom networks, respectively.

2 Materials and methods

2.1 Participants and procedures

A cross-sectional study was carried out in Chongqing, China, from 11 March to 17 April 2021, as part of a large mental health screening focusing on depression and anxiety in vocational high school students. Two schools were randomly chosen, and all their first-year students participated. Exclusion criteria were incomplete questionnaire data and failing an attention check (e.g., not selecting the specified option on control questions to filter out inattentive respondents). Of 8,923 submitted questionnaires, 7,471 were valid (83.73%).

The survey utilized “Questionnaire Star,” a popular online tool in China. Students and their guardians were informed about the study’s aim, methods, risks, and benefits before participation. Digital consent was obtained from participants and their guardians, following ethical research standards. The survey took 20–30 min, was conducted in computer classrooms, with instructors facilitating and ensuring access. They were assured that participation was voluntary and anonymous, only those who provided consent and signed the informed consent form proceeded to complete the questionnaire, and data would be kept confidentially and used for research only. The study was approved by the Medical Ethics Committee of Shenzhen University School of Medicine (Granting number: 2020005).

2.2 Measures

Depressive symptoms were measured using the Patient Health Questionnaire (PHQ-8), consisting of eight items each representing a depressive symptom. These are rated on a four-point Likert scale from 0 (not at all) to 3 (nearly every day), with the total score indicating depression severity over the past 2 weeks (Kroenke et al., 2009). The reliability and validity of the PHQ-8 has been established in prior research (Kroenke et al., 2009). Although specific validation for the Chinese PHQ-8 is limited, the Chinese version of the PHQ-9, including an additional item assessing suicidal ideation, showed good reliability and validity (Wang et al., 2014). In this study, the PHQ-8 showed excellent internal consistency (Cronbach’s alpha = 0.886) and construct validity, confirmed by confirmatory factor analysis with satisfactory fit indices (χ2/df = 27.165, CFI = 0.980, TLI = 0.972, RMSEA = 0.059, SRMR = 0.022).

Anxiety symptoms were measured using the Generalized Anxiety Disorder 7-item Scale (GAD-7). This instrument evaluates seven symptoms of anxiety disorders on a four-point Likert scale, ranging from 0 (not at all) to 3 (nearly every day), with the total score reflecting the severity of anxiety symptoms over the past 2 weeks. The Chinese version has been validated in Chinese adolescent (Sun et al., 2021). In the current study, the GAD-7 also demonstrated excellent internal consistency (Cronbach’s alpha = 0.909) and construct validity, confirmed by confirmatory factor analysis with strong fit indices (χ2/df = 28.122, CFI = 0.982, TLI = 0.970, RMSEA = 0.060, SRMR = 0.021).

2.3 Statistic analysis

Descriptive statistics were analyzed using the “bruceR (version 2023.9)” package (Bao, 2023). Our network analysis included a Gaussian Graphical Model (GGM) for undirected networks and a Directed Acyclic Graph (DAG) for directed networks. In GGM, we followed the guideline proposed by Epskamp and Fried (2018a,b), including network estimation and visualization, and calculation of network centrality, network stability and accuracy. Bridge centrality in GGM was measured using Jones et al. (2021) bridge strength method. For the DAG analysis, we adhered to Briganti et al. (2023) tutorial and referenced McNally et al. (2017). Gender differences in network structures were examined using the Network Comparison Test (NCT) developed by van Borkulo et al. (2022).

2.3.1 GGM network

Network estimation and visualization. We modeled symptoms as nodes and their partial correlations as edges using Spearman’s matrices. Considering our large dataset with 35 variables, we chose unregularized models, as recommended by Williams and Rast (2020). We identified the GGM using the ggmModSelect algorithm, which selects the best GGM without regularization, employing a repetitive method and the extended Bayesian information criterion (Williams and Rast, 2020). The analysis was conducted using the R packages “bootnet (version 1.5)” and “qgraph (version1.9.8)” (Epskamp and Fried, 2018b; Epskamp et al., 2012). For network visualization, we applied the Fruchterman and Reingold (1991) algorithm, clustering related nodes and placing less connected ones peripherally.

Node’s bridge centrality and centrality. We used the “networkTools” and “qgraph” R packages to calculate bridge and strength centralities, selecting bridge strength and strength as node attributes (Epskamp et al., 2012; Jones and Jones, 2017). These robust metrics, highlighted in prior studies, quantify node importance (Bringmann et al., 2019; Opsahl et al., 2010). Bridge strength centrality sums absolute edge weights linking nodes across symptom clusters or disorders, identifying key “bridges” (Jones et al., 2021). Strength centrality, adding absolute weights of all connected edges, indicates a node’s influence on network connectivity (Opsahl et al., 2010). These measure identifies bridge symptoms that potentially serve as critical connectors within the network, reflecting their role in the co-occurrence and interaction of depressive and anxiety symptoms.

Network stability and accuracy estimation. To confirm the accuracy and stability of network, we used the bootnet package in R, following the procedures by Epskamp and Fried (2018a,b). In this context, “accuracy” refers to the precision of the estimated edge weights in the network—that is, the extent to which these estimates approximate their true values. We evaluated accuracy by generating 1,000 bootstrap samples to compute 95% confidence intervals (CIs) for the edge weights; narrower CIs indicate higher accuracy. Similarly, “stability” refers to the consistency of the network structure (particularly the centrality indices) when the data is subjected to variations or when subsets of the data are used. We assessed stability through a case-dropping bootstrap method, calculating the CS-coefficient from 1,000 samples. The CS-coefficient (ideally above 0.25) reflects the proportion of the sample that can be dropped while still retaining a correlation of at least 0.7 between the full dataset and the subset estimates of the centrality indices (Epskamp and Fried, 2018a; Epskamp and Fried, 2018b). Lastly, we applied a bootstrap difference test (1,000 samples, alpha = 0.05) to assess network edge and node differences, comparing bootstrapped CIs of edge weights and centrality indices (Epskamp and Fried, 2018a; Epskamp and Fried, 2018b).

2.3.2 Bayesian network analysis

To explore potential causal relationships in the symptom network, we used Directed Acyclic Graph (DAG) to estimate the Bayesian network with the bnlearn package in R 4.2.3 (Scutari, 2009). DAGs help infer potential causal pathways by representing conditional dependencies between variables. We employed a bootstrapped Hill-Climbing (HC) algorithm, following Briganti et al. (2023). The HC algorithm, a score-based method, iteratively modifies edges to maximize a score while keeping the graph acyclic (Russell and Norvig, 2010). Our analysis was strengthened by creating 1,000 bootstrap samples and configuring multiple restarts to avoid local optima, following McNally et al. (2017). After bootstrapping, we synthesized an averaged network from all samples, evaluated using the Bayesian Information Criterion (BIC) for model fit. We also computed arc strengths using the BIC-Gaussian score method to determine the certainty of each directed edge. The final DAG, visualized using the bnlearn package, showed elliptical nodes and edges weighted by arc strengths, highlighting the most reliable causal relationships between symptoms.

2.3.3 Network comparisons between gender

To compare the network structures between different groups, we employed the “Network Comparison Test (NCT)” using the corresponding R package (van Borkulo et al., 2022). NCT enables statistical comparisons of networks using metrics like global strength (total of all absolute edge weights), individual edge weights, and centrality measures such as strength. This allows NCT to identify significant differences in overall network connectivity and the roles/importance of specific nodes and edges. Considering extensive research has consistently shown that females are nearly twice as likely as males to experience depression in their lifetime (Kuehner, 2017), we applied NCT to contrast networks from female and male groups. Thus, we could know whether the observed gender differences extend beyond mere prevalence rates and also reflect differences in the interrelationships among individual symptoms.

3 Results

3.1 Descriptive statistics

Among 7,471 participants (average age 16.37 ± 0.94), 59.29% were males (N = 4,428, average age = 16.39 ± 0.89), and 40.71% were females (N = 3,043, average age = 16.35 ± 0.92). The average scores for the PHQ-8 and GAD-7, along with the prevalence of depression and anxiety defined by scores of 10 or higher (Kroenke et al., 2009; Spitzer et al., 2006), are detailed in Table 1. It is noteworthy that the prevalence and total score of depression and anxiety were significantly higher among female participants than male counterparts (ps < 0.001). Moreover, females scored significantly higher than males in depressive symptoms and anxiety symptoms, as shown in Supplementary Table 1. The comorbidity rate was 11.31% in the total sample, with no significant difference between genders.

TABLE 1
www.frontiersin.org

Table 1. Participant demographics and descriptive statistics of depression and anxiety.

3.2 GGM network analysis of depressive symptoms and anxiety symptoms

Network estimation and visualization. Figure 1 displays the GGM network for depressive and anxiety symptoms. Out of 105 potential node pairs, 77 showed associations, leading to a network density of 0.733 and an average edge weight of 0.09. Edge weight details are in the Supplementary Table 2. The network visualization shows clusters of depressive and anxiety symptoms, indicating their inter-connectivity. All network edges were positive, suggesting a positive correlation between these symptoms.

FIGURE 1
Network diagram depicting relationships between anxiety (GAD) and depression (PHQ) symptoms. Orange nodes represent anxiety symptoms, and blue nodes represent depression symptoms. Lines indicate connections, with varying thickness to show strength. Legends specify symptoms for both categories.

Figure 1. Network structure of depressive and anxiety symptoms. (1) Nodes are color-coded by DSM-5 symptom clusters: orange for anxiety symptoms, blue for depressive symptoms. (2) Edge colors indicate correlations: purple for positive, red for negative, with thickness denoting the strength of pairwise correlations – thicker edges signify stronger correlations.

Node’s Bridge Centrality and Centrality. Figure 2a shows node bridge strengths, and Supplementary Figure 1 presents results from the non-parametric bootstrap difference test. With the highest bridge strength and significant differences than others (P < 0.05, Supplementary Figure 1), sad mood was the bridge symptom of the network. Figure 2b illustrates node strengths, with Fatigue having the highest. Additionally, Supplementary Figure 2’s non-parametric bootstrap difference test indicates that loss of control, trouble relaxing, sad mood, irritability, excessive worry, feeling guilty have significantly higher strengths than other symptoms, these symptoms are thus considered central within the network.

FIGURE 2
Two line graphs labeled (a) and (b) compare variables related to “Bridge Strength” and “Strength”. Both graphs plot different variables: PHQ9, PHQ7, PHQ3, and others, on the vertical axis against a horizontal scale from -1 to 2. In graph (a), the variables are more evenly distributed across the scale, while in graph (b), they cluster closer to the center.

Figure 2. The variables are more evenly distributed across the scale, while in graph (b), they cluster closer to the center. (a) Bridge strength (Z-standardized); (b) Strength (Z-standardized). GAD1, nervousness; GAD2, loss of control; GAD3, excessive worry; GAD4, trouble relaxing; GAD5, restlessness; GAD6, irritability; GAD7, feeling afraid; PHQ1, anhedonia; PHQ2, sad mood; PHQ3, sleep problems; PHQ4, fatigue; PHQ5, appetite changes; PHQ6, guilt; PHQ7, difficulty concentrating; PHQ8-motor.

Network stability and accuracy. All results are in Supplementary Figure 3, both the nodes and edges of the network have high stability and accuracy. The results of the Subsetting Bootstrap are shown in Supplementary Figure 3a, and the stability coefficients (CS-Coefficient) of node bridge strength and strength are all 0.75. The results of the edge weight bootstrap procedure are shown in Supplementary Figure 3b. The edge weights in the current sample are the same as the bootstrap sample, which means the network has high stability.

3.3 Bayesian network analysis of depressive symptoms and anxiety symptoms

The DAG analysis results, shown in Figures 3, 4, map potential causal links between depressive and anxiety symptoms. Figure 3 offers an overview of all symptom interactions, and Figure 4 focuses on connections between depressive and anxiety symptoms. In Figure 3, irritability emerges as the initial node, highlighting its potential role in triggering or maintaining other symptoms, thus being a critical target for early intervention. Figure 4 shows a dominant directional flow from anxiety to depressive symptoms, with 31 edges from anxiety to depressive symptoms, suggesting anxiety’s role in developing depressive symptoms. Conversely, the transition from depressive to anxiety symptoms is less evident, with only four edges, mainly from sad mood and one from guilty. This asymmetry indicates that while depressive symptoms can lead to anxiety, this is less common than the reverse.

FIGURE 3
Network diagram illustrating connections between anxiety and depression indicators. Anxiety indicators (GAD1-GAD7) include nervousness, worry control, and fear. Depression indicators (PHQ1-PHQ8) include anhedonia, sad mood, and fatigue. Lines of varying thickness show relationships between items. Arrows indicate the direction of influence between nodes. Legend lists specific symptoms for each node.

Figure 3. Network diagram illustrating connections between anxiety and depression indicators. Arrows indicate the direction of influence between nodes. Legend lists specific symptoms for each node.

FIGURE 4
Network diagram showing relationships between anxiety and depression symptoms. Anxiety nodes (GAD1-7) are orange; depression nodes (PHQ1-8) are blue. Arrows indicate the direction of influence between nodes. Legend lists specific symptoms for each node.

Figure 4. Subgraph of directed edges between depressive and anxiety symptoms. (1) A Subgraph of Directed acyclic graphs (DAG) network nodes represent symptoms and anxiety symptoms and edges represent directed connections between them. (2) The thickness of the arrows represents directional probability.

3.4 Network comparison test between gender

The networks showed notable gender differences in global strength, with females’ network exhibiting a higher global strength (7.04) than the males’ network (6.90; P < 0.01), which suggests a density connections of symptoms among females (Borsboom and Cramer, 2013; Cramer et al., 2010). Specifically, the symptom of anhedonia presented a significantly higher strength in the females’ network than males’ network (diff_PHQ1 = −0.17, P < 0.05), pointing to its greater relevance or impact among females.

4 Discussion

As the first study to apply network analysis specifically to depressive and anxiety symptoms among Chinese vocational high school students, this research systematically identified critical nodes within both directed (Bayesian-directed acyclic graph, DAG) and undirected (Gaussian graphical model, GGM) symptom networks. Our approach revealed one bridge symptom and seven core symptoms, explored potential directional relationships between depressive and anxiety symptoms, and examined gender differences in both prevalence and network structure. Furthermore, we confirmed the stability and accuracy of our network models, ensuring robust conclusions. We found one bridge symptom and seven core symptoms, uncovered bidirectional associations between symptoms and gender differences in GGM network structure. These findings provide novel insights from a network perspective. To validate these exploratory findings, future studies should prioritize longitudinal designs to track symptom dynamics and experimental interventions targeting identified core/bridge symptoms.

A key finding of our study is the mutual influence between depressive and anxiety symptoms, with a more pronounced trend from anxiety symptoms to depressive symptoms. This pattern suggests that anxiety may be an earlier and more influential factor in the development of comorbidity. It aligns with existing literature indicating that anxiety can predict depression, with the predictive effect of anxiety on depression being stronger in the short term and influenced by the baseline level of anxiety (Jacobson and Newman, 2017), potentially acting as a catalyst for the development of depressive symptoms. Our study conducted following a new COVID-19 outbreaks, students faced the uncertainty and threat brought about by the pandemic may lead to increased anxiety symptoms in individuals (Bendau et al., 2021), and the absence of effective coping strategies lead to an exacerbation of anxiety and depressive symptoms (Miola et al., 2023). Thus, the development and implementation of positive coping strategies are crucial in mitigating the severity of these symptoms and promoting mental health resilience.

Notably, irritability emerged as the initial node in the DAG network, influencing all other symptoms and highlighting its importance for preventive interventions. Irritability also emerged as a core symptom in Chinese female nursing students and a bridge symptom in epilepsy patients (He et al., 2023; Ren et al., 2021). A recent meta-analysis of 121 studies confirmed irritability as a transdiagnostic status across youth mental-health disorders (Kim and Asbury, 2020). Cross-lagged-panel work further indicates that elevations in irritability at Time 1 precipitate increases in worry and somatic tension six months later, but not vice-versa (Chin et al., 2025). Considering our research carried out followed a new COVID-19 outbreaks, with a significant reductions in social activities (Losiewicz et al., 2023), and heightened perceptions of risk (Webb, 2021). These threat factors, uncertainty, and isolation during COVID-19 pandemics, have been consistently associated with heightened stress responses and the emergence of anxiety symptoms (Buheji et al., 2020). Because our DAG implies top-down influence, targeting irritability may yield a cascade effect. Irritability should therefore be prioritized for early identification and intervention to prevent further mental health issues in vocational high school students. Youth-specific RCT already showed medium-to-large reductions in global impairment after irritability-focused program—an exposure-based CBT protocol. Adapting these brief (6–8 session) formats to homeroom periods could provide cost-effective, population-level prevention. We recommend incorporating a single-item irritability screener into routine school health checks and intervention in frustration-tolerance skills (Naim et al., 2024).

Importantly, sad mood exerts a significant influence on symptoms transition across depression and anxiety. Identified as a bridge symptom in the GGM network, consistent with various studies (He et al., 2023; Ren et al., 2021; Wang et al., 2020), and was found to demonstrate potential directional relationships with both anxiety and depression symptoms in the DAG network for the first time. This bi-directionality indicates a complex interplay where sad mood not only results from but also contributes to the exacerbation of other symptoms within the network. Our results corroborate Clark and Watson’s (1991) tripartite model of anxiety and depression, which posits negative affectivity, like sad mood and despair, as a shared characteristic of both conditions. The significant role of these symptoms in cross-disorder symptom transmission becomes evident due to their shared nature. Building on Kaiser et al. (2021) recommendation for targeted interventions on bridge symptoms, and giving the stronger evidence supporting the role of sad mood in the network, we advocate for therapeutic interventions focused on this symptom to potentially reduce the occurrence and progression of comorbid depression and anxiety.

Moreover, fatigue, loss of control, trouble relaxing, sad mood, irritability, excessive worry, and guilt were the core symptoms. Notably, some researchers found similar core symptoms (e.g., low energy and trouble relaxing) in university students, but with physical activity and stress management as key protective factors—a contrast highlighting the need for age-specific interventions (Sun et al., 2024). Regarding depressive symptoms, while fatigue and sad mood are well-recognized core symptoms (Beard et al., 2016; He et al., 2023; Kaiser et al., 2021), the prominence of guilt is a novel finding, possibly reflecting the cultural and educational context in East Asian Confucian societies, where academic failure often leads to feelings of inadequacy or guilt (Fwu et al., 2017). Regarding anxiety symptoms, their identification as core may be attributed to the timing of our study during the COVID-19 pandemic. Another network analysis conducted during the COVID-19 pandemic also highlighted the central role of anxiety symptoms in comorbid presentations. Loss of control was also a core symptom in nursing students during the COVID-19 (Bai et al., 2021), irritability and trouble relaxing showed high centrality during the after peak stage of Covid-19 pandemic (Naim et al., 2024). These findings emphasize the significant impact of anxiety symptoms on mental health during the pandemic and suggest targeting these symptoms for prevention and intervention in similar public health crises.

Finally, the study found significant gender differences in the prevalence of depression and anxiety, network global strength, and the strength of anhedonia, all of which were higher in females. This suggests more severe symptoms and a more connected network in the female. Such gender disparity aligns with literature indicating females experience more severe depression and anxiety challenges (Bang et al., 2020). The emergence of this gender difference is consistent with Sun et al. (2024)’s findings in primary/middle school students, in which screen time had stronger associations with negative emotions in females. This suggests that gender-specific vulnerability patterns may persist across developmental stages, though the contributing factors may differ. The stronger symptom connection network may be one of the possible reasons to explain why the depression and anxiety prevalence in females is higher than in males. According to Borsboom (2017), in tightly connected networks, symptoms can sustain activity even after initial triggers subside, leading to disorder states. Thus, females, with a denser network, are more susceptible susceptible to comorbid depression and anxiety. Additionally, anhedonia showed significantly higher strength in females, but literature on its gender-specific impact is mixed (Bolton et al., 2018; Chan et al., 2012), indicating a need for further research. In general, our study not only identifies a higher prevalence of comorbid depression and anxiety in females but also, from a network analysis perspective, sheds light on the potential underlying mechanisms.

This study utilized network analysis to examine depression and anxiety symptom patterns and gender differences in vocational high school students. In doing so, we have provided insights into symptom-level understanding of depression and anxiety in this large but overlooked target group, and provided recommendations for subsequent prevention and intervention. We discovered potential directional relationships in the spread of symptoms across diseases, as well as key symptoms within the network structure. Based on these findings and relevant theories, we proposed directions for future interpretation and practical application based on relevant theories and our findings, and verified the existence of gender differences, providing an explanation from a network perspective. Firstly, we recognized anxiety symptoms, especially irritability as early indicators for comorbidity onset and progression, advocating for prompt preventive measures. Concurrently, we suggested therapeutic interventions focused on sad mood to effectively prevent further development and exacerbation of comorbid depression and anxiety. Lastly, considering the higher prevalence and susceptibility in females, we recommend prioritizing mental health resources for female students in settings with limited resources.

The network analysis and sample used in this study to reveal patterns of depression and anxiety symptoms have several limitations. First, the study’s cross-sectional nature means the DAG network only suggests potential directional relationships. Time-series or longitudinal studies are needed for deeper insights into symptom activation spread and symptom causality (Forbes et al., 2017). Second, only DSM-5 symptoms of anxiety and depression were included. Non-DSM-5 symptoms could also be crucial in these networks (Fried et al., 2016). This may limit the comprehensiveness of the identified symptom interactions. Future studies could expand symptom inventories to capture broader clinical profiles. Future research should consider a broader range of symptoms to more accurately explain comorbidity. Third, the present study is limited to first-year vocational high school students in China, without comparing their networks to those of students from other grades or regular high schools, which limits the understanding of both commonalities and differences across different student groups. Additionally, the regional focus on Chongqing may constrain the generalizability of findings to other socioeconomic or cultural contexts in China. Future studies could expand to multi-grade cohorts and diverse school types to validate these patterns. Network connectivity differences may exist across different samples or populations (Schultz et al., 2019). Fourth, while our study provides novel network insights into depression and anxiety among vocational students, we did not assess developmental factors such as identity formation or hormonal changes that may influence symptom networks across different adolescent populations (Sun et al., 2024; Sun et al., 2025). Fifth, approximately 16% of responses were ultimately excluded—primarily due to failed validity checks or incomplete data. These exclusions were necessary to ensure the reliability of the network analysis; however, we acknowledge that students experiencing greater psychological distress (e.g., severe anxiety or depression) may have been more likely to provide invalid or incomplete responses. Further research comparing these groups could lead to more tailored interventions. Including clinical samples in future studies could also help verify or challenge these findings.

Despite its limitations, the present study yields critical insights into the depression and anxiety among secondary vocational school students. These implications are crucial for both educational policy and mental health intervention strategies. Firstly, early psychological assessment is crucial, with irritability and sad mood as key indicators. Systematic screening for these symptoms in schools can enable early intervention, potentially altering comorbidity trajectories. Secondly, the pronounced gender disparity in comorbidity, particularly among female students, underscores the need for gender-specific psychological support strategies, such as tailored support groups, counseling, and psycho-educational programs. However, future research should prioritize elucidating the mechanisms underlying this disparity rather than solely detecting it. Potential contributors warrant systematic investigation. Longitudinal designs integrating multi-modal data could disentangle whether observed network density differences reflect biological predispositions, contextual stressors, or their interaction. Such mechanistic insights are critical for developing interventions that address root causes rather than symptomatic manifestations of the gender gap. Third, our findings are consistent with recent studies demonstrating that Chinese students continued to experience significant mental health challenges even after the relaxation of COVID-19 restrictions (Cai et al., 2025). These persistent effects, particularly among secondary school students, suggest that the psychological impacts of the pandemic have persisted long after the acute crisis period, reinforcing the continued relevance of our findings in the post-pandemic educational context. Lastly, integrating mental health literacy into the vocational curriculum, with a focus on psychological resilience and access to support, is essential. This approach aims to foster a mentally resilient and informed student body.

5 Conclusion

The network analysis identified sad mood (PHQ2) as the primary bridge symptom linking depression and anxiety, while fatigue (PHQ4) and six other symptoms (loss of control, trouble relaxing, irritability, excessive worry, guilt, and difficulty concentrating) emerged as core nodes with high centrality. Directed acyclic graph analysis further revealed that anxiety symptoms (e.g., irritability) predominantly activated depressive symptoms, with irritability (GAD6) serving as the initial trigger in the symptom cascade. Notably, female participants exhibited stronger global network connectivity, suggesting denser symptom interactions that may underlie their higher prevalence of comorbid depression and anxiety. These findings underscore the urgency of gender-sensitive interventions targeting bridge and core symptoms—particularly irritability and sad mood—to disrupt symptom propagation pathways in this vulnerable population.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: if data sets are required, contact the corresponding author.

Ethics statement

The studies involving humans were approved by the Medical Ethics Committee of Shenzhen University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.

Author contributions

YH: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing. DY: Conceptualization, Data curation, Writing – original draft, Writing – review and editing. KL: Conceptualization, Writing – review and editing. YZ: Conceptualization, Data curation, Writing – original draft. HZ: Funding acquisition, Project administration, Writing – review and editing. XC: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by Humanities and Social Sciences of Ministry of Education Planning Foundation (Grant 23YJA190002), Guangdong Philosophy and Social Science Planning 2023 Co-construction Project (Grant GD23XXL17), Revitalization Project of Humanities and Social Sciences, Shenzhen University (Grant WKZX0311), Youth Research Project of Shenzhen University 2035 Outstanding Research Plan (Grant ZYQN2311), and National Social Science Foundation (Grant 24BSH048).

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.

Generative AI statement

The authors declare that no Generative AI was used in the creation of this manuscript.

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/fpsyg.2025.1572011/full#supplementary-material

References

Bai, W., Cai, H., Liu, S., Chen, X., Sha, S., Cheung, T., et al. (2021). Anxiety and depressive symptoms in college students during the late stage of the COVID-19 outbreak: A network approach. Transl. Psychiatry 11:638. doi: 10.1038/s41398-021-01738-4

PubMed Abstract | Crossref Full Text | Google Scholar

Bang, H., Won, D., and Park, S. (2020). School engagement, self-esteem, and depression of adolescents: The role of sport participation and volunteering activity and gender differences. Child. Youth Serv. Rev. 113:105012. doi: 10.1016/j.childyouth.2020.105012

Crossref Full Text | Google Scholar

Bao, H. W. S. (2023). bruceR: Broadly useful convenient and efficient R functions. Available online at: https://psychbruce.github.io/bruceR/ (accessed November 26, 2023).

Google Scholar

Beard, C., Millner, A. J., Forgeard, M. J., Fried, E. I., Hsu, K. J., Treadway, M. T., et al. (2016). Network analysis of depression and anxiety symptom relationships in a psychiatric sample. Psychol. Med. 46, 3359–3369. doi: 10.1017/S0033291716002300

PubMed Abstract | Crossref Full Text | Google Scholar

Bendau, A., Petzold, M. B., Wyka, S., Pyrkosch, L., Plag, J., and Ströhle, A. (2021). [Anxiety in times of COVID-19 and other health crises]. Nervenarzt 92, 417–425. doi: 10.1007/s00115-020-01030-8

PubMed Abstract | Crossref Full Text | Google Scholar

Bolton, J. L., Molet, J., Regev, L., Chen, Y., Rismanchi, N., Haddad, E., et al. (2018). Anhedonia following early-life adversity involves aberrant interaction of reward and anxiety circuits and is reversed by partial silencing of amygdala corticotropin-releasing hormone gene. Biol. Psychiatry 83, 137–147. doi: 10.1016/j.biopsych.2017.08.023

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Borsboom, D., Cramer, A. O., Schmittmann, V. D., Epskamp, S., and Waldorp, L. J. (2011). The small world of psychopathology. PLoS One 6:e27407. doi: 10.1371/journal.pone.0027407

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Briganti, G., Scutari, M., and McNally, R. J. (2023). A tutorial on bayesian networks for psychopathology researchers. Psychol. Methods 28, 947–961. doi: 10.1037/met0000479

PubMed Abstract | Crossref Full Text | Google Scholar

Bringmann, L. F., Elmer, T., Epskamp, S., Krause, R. W., Schoch, D., Wichers, M., et al. (2019). What do centrality measures measure in psychological networks? J. Abnorm. Psychol. 128, 892–903. doi: 10.1037/abn0000446

PubMed Abstract | Crossref Full Text | Google Scholar

Buheji, M., Jahrami, H., and Sabah Dhahi, A. (2020). Minimising stress exposure during pandemics similar to COVID-19. IJPBS 10, 9–16. doi: 10.5923/j.ijpbs.20201001.02

PubMed Abstract | Crossref Full Text | Google Scholar

Cai, H., Bai, W., Liu, H., Chen, X., Qi, H., Liu, R., et al. (2022). Network analysis of depressive and anxiety symptoms in adolescents during the later stage of the COVID-19 pandemic. Transl. Psychiatry 12:98. doi: 10.1038/s41398-022-01838-9

PubMed Abstract | Crossref Full Text | Google Scholar

Cai, J., Wang, Y., Wang, C., Deng, Z. Y., Mu, Y. F., Deng, A. P., et al. (2025). Mental health problems and associated factors of students at different learning stages during and after COVID-19 pandemic in Sichuan, China. BMC Psychiatry 25:144. doi: 10.1186/s12888-025-06588-8

PubMed Abstract | Crossref Full Text | Google Scholar

Chan, R. C., Wang, Y., Yan, C., Zhao, Q., McGrath, J., Hsi, X., et al. (2012). A study of trait anhedonia in non-clinical Chinese samples: Evidence from the chapman scales for physical and social anhedonia. PLoS One 7:e34275. doi: 10.1371/journal.pone.0034275

PubMed Abstract | Crossref Full Text | Google Scholar

Chin, M., Robson, D. A., Woodbridge, H., and Hawes, D. J. (2025). Irritability as a transdiagnostic construct across childhood and adolescence: A systematic review and meta-analysis. Clin. Child Fam. Psychol. Rev. 28, 101–124. doi: 10.1007/s10567-024-00512-4

PubMed Abstract | Crossref Full Text | Google Scholar

Clark, L. A., and Watson, D. (1991). Tripartite model of anxiety and depression: Psychometric evidence and taxonomic implications. J. Abnorm. Psychol. 100, 316–336. doi: 10.1037//0021-843x.100.3.316

PubMed Abstract | Crossref Full Text | Google Scholar

Coledam, D. H. C., Frotta, B. M., and Ré, A. H. N. (2022). General versus vocational education in high school: Cross-sectional associations with student’s health. J. Sch. Health 92, 570–580. doi: 10.1111/josh.13165

PubMed Abstract | Crossref Full Text | Google Scholar

Cramer, A. O., Waldorp, L. J., van der Maas, H. L., and Borsboom, D. (2010). Comorbidity: A network perspective. Behav. Brain Sci. 33:137–150; discussion 150–193. doi: 10.1017/S0140525X09991567.

PubMed Abstract | Crossref Full Text | Google Scholar

Cummings, C. M., Caporino, N. E., and Kendall, P. C. (2014). Comorbidity of anxiety and depression in children and adolescents: 20 years after. Psychol. Bull. 140, 816–845. doi: 10.1037/a0034733

PubMed Abstract | Crossref Full Text | Google Scholar

Epskamp, S., Cramer, A. O. J., Waldorp, L. J., Schmittmann, V. D., and Borsboom, D. (2012). qgraph: Network visualizations of relationships in psychometric data. J. Stat. Soft. 48, 1–18. doi: 10.18637/jss.v048.i04

Crossref Full Text | Google Scholar

Epskamp, S., and Fried, E. I. (2018a). A tutorial on regularized partial correlation networks. Psychol. Methods 23, 617–634. doi: 10.1037/met0000167

PubMed Abstract | Crossref Full Text | Google Scholar

Epskamp, S., and Fried, E. I. (2018b). bootnet: Bootstrap methods for various network estimation routines (R package version 1.1.0). Vienna: R Foundation for Statistical Computing. Available online at: https://CRAN.R-project.org/package=bootnet

Google Scholar

Fang, G., Wang, Y., Yuan, H., Yan, N., and Zhi, S. (2024). Unraveling the core symptoms of mental health in senior grade three students- a network analysis. Front. Psychiatry 15:1364334. doi: 10.3389/fpsyt.2024.1364334

PubMed Abstract | Crossref Full Text | Google Scholar

Forbes, M. K., Wright, A. G. C., Markon, K. E., and Krueger, R. F. (2017). Evidence that psychopathology symptom networks have limited replicability. J. Abnorm. Psychol. 126, 969–988. doi: 10.1037/abn0000276

PubMed Abstract | Crossref Full Text | Google Scholar

Fried, E. I., Epskamp, S., Nesse, R. M., Tuerlinckx, F., and Borsboom, D. (2016). What are ‘good’ depression symptoms? Comparing the centrality of DSM and non-DSM symptoms of depression in a network analysis. J. Affect. Disord. 189, 314–320. doi: 10.1016/j.jad.2015.09.005

PubMed Abstract | Crossref Full Text | Google Scholar

Fruchterman, T. M. J., and Reingold, E. M. (1991). Graph drawing by force-directed placement. Softw. Pract. Exp. 21, 1129–1164. doi: 10.1002/spe.4380211102

Crossref Full Text | Google Scholar

Fwu, B. J., Wang, H. H., Chen, S. W., and Wei, C. F. (2017). ‘Feeling bad’or ‘being bad?’The trapping effect of effort in academic failure in a Confucian cultural context. Educ. Psychol. 37, 506–519. doi: 10.1080/01443410.2016.1152355

Crossref Full Text | Google Scholar

Garabiles, M. R., Lao, C. K., Xiong, Y., and Hall, B. J. (2019). Exploring comorbidity between anxiety and depression among migrant Filipino domestic workers: A network approach. J. Affect. Disord. 250, 85–93. doi: 10.1016/j.jad.2019.02.062

PubMed Abstract | Crossref Full Text | Google Scholar

Garber, J., and Weersing, V. R. (2010). Comorbidity of anxiety and depression in youth: Implications for treatment and prevention. Clin. Psychol. 17, 293–306. doi: 10.1111/j.1468-2850.2010.01221.x

PubMed Abstract | Crossref Full Text | Google Scholar

He, Y., Wu, C., LeMoult, J., Huang, J., Zhao, Y., Liang, K., et al. (2023). Exploring symptom-level associations between anxiety and depression across developmental stages of adolescence: A network analysis approach. BMC Psychiatry 23:941. doi: 10.1186/s12888-023-05449-6

PubMed Abstract | Crossref Full Text | Google Scholar

Huang, Y., Wang, Y., Wang, H., Liu, Z., Yu, X., Yan, J., et al. (2019). Prevalence of mental disorders in China: A cross-sectional epidemiological study. Lancet Psychiatry 6, 211–224. doi: 10.1016/S2215-0366(18)30511-X

PubMed Abstract | Crossref Full Text | Google Scholar

Hung, C. I., Liu, C. Y., Yang, C. H., and Gan, S. T. (2020). Comorbidity with more anxiety disorders associated with a poorer prognosis persisting at the 10-year follow-up among patients with major depressive disorder. J. Affect. Disord. 260, 97–104. doi: 10.1016/j.jad.2019.08.085

PubMed Abstract | Crossref Full Text | Google Scholar

Jacobson, N. C., and Newman, M. G. (2017). Anxiety and depression as bidirectional risk factors for one another: A meta-analysis of longitudinal studies. Psychol. Bull. 143, 1155–1200. doi: 10.1037/bul0000111

PubMed Abstract | Crossref Full Text | Google Scholar

Jones, P. J., Ma, R., and McNally, R. J. (2021). Bridge centrality: A network approach to understanding comorbidity. Multiv. Behav. Res. 56, 353–367. doi: 10.1080/00273171.2019.1614898

PubMed Abstract | Crossref Full Text | Google Scholar

Jones, P., and Jones, M. P. (2017). Package ‘networktools.’. Available online at: https://cran.rproject.org/web/packages/networktools/networktools.pdf (accessed June, 2018).

Google Scholar

Kaiser, T., Herzog, P., Voderholzer, U., and Brakemeier, E. L. (2021). Unraveling the comorbidity of depression and anxiety in a large inpatient sample: Network analysis to examine bridge symptoms. Dep. Anxiety 38, 307–317. doi: 10.1002/da.23136

PubMed Abstract | Crossref Full Text | Google Scholar

Kim, L. E., and Asbury, K. (2020). ‘Like a rug had been pulled from under you’: The impact of COVID-19 on teachers in England during the first six weeks of the UK lockdown. Br. J. Educ. Psychol. 90, 1062–1083. doi: 10.1111/bjep.12381

PubMed Abstract | Crossref Full Text | Google Scholar

Kroenke, K., Strine, T. W., Spitzer, R. L., Williams, J. B., Berry, J. T., and Mokdad, A. H. (2009). The PHQ-8 as a measure of current depression in the general population. J. Affect. Disord. 114, 163–173. doi: 10.1016/j.jad.2008.06.026

PubMed Abstract | Crossref Full Text | Google Scholar

Kuehner, C. (2017). Why is depression more common among women than among men? Lancet Psychiatry 4, 146–158. doi: 10.1016/S2215-0366(16)30263-2

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, J., Teng, Z., Chen, Z., Wei, Z., Zou, T., Qin, Y., et al. (2022). Exploring the associations between behavioral health risk factors, abnormal eating attitudes and socio-demographic factors among Chinese youth: Survey of 7,984 vocational high school students in Hunan in 2020. Front. Psychiatry 13:1000821. doi: 10.3389/fpsyt.2022.1000821

PubMed Abstract | Crossref Full Text | Google Scholar

Losiewicz, O. M., Metts, A. V., Zinbarg, R. E., Hammen, C., and Craske, M. G. (2023). Examining the indirect contributions of irritability and chronic interpersonal stress on symptoms of anxiety and depression in adolescents. J. Affect. Disord. 329, 350–358. doi: 10.1016/j.jad.2023.02.125

PubMed Abstract | Crossref Full Text | Google Scholar

McNally, R. J., Mair, P., Mugno, B. L., and Riemann, B. C. (2017). Co-morbid obsessive-compulsive disorder and depression: A Bayesian network approach. Psychol. Med. 47, 1204–1214. doi: 10.1017/S0033291716003287

PubMed Abstract | Crossref Full Text | Google Scholar

Miola, A., Caiolo, S., Pontoni, G., Pozzan, E., Moriglia, C., Simionato, F., et al. (2023). Anxiety and depression during the second wave of the COVID-19 pandemic: The role of coping strategies. Int. J. Environ. Res. Public Health 20:2974. doi: 10.3390/ijerph20042974

PubMed Abstract | Crossref Full Text | Google Scholar

Naim, R., Dombek, K., German, R. E., Haller, S. P., Kircanski, K., and Brotman, M. A. (2024). An exposure-based cognitive-behavioral therapy for youth with severe irritability: Feasibility and preliminary efficacy. J. Clin. Child Adolesc. Psychol. 53, 260–276. doi: 10.1080/15374416.2023.2264385

PubMed Abstract | Crossref Full Text | Google Scholar

Opsahl, T., Agneessens, F., and Skvoretz, J. (2010). Node centrality in weighted networks: Generalizing degree and shortest paths. Soc. Netw. 32, 245–251. doi: 10.1016/j.socnet.2010.03.006

Crossref Full Text | Google Scholar

Pearl, J. (2009). Causality: Models, reasoning, and inference, 2nd Edn. Cambridge, MA: Cambridge University Press.

Google Scholar

Pearl, J., and Mackenzie, D. (2018). The book of why: The new science of cause and effect. New York, NY: Basic books.

Google Scholar

Ren, L., Wang, Y., Wu, L., Wei, Z., Cui, L. B., Wei, X., et al. (2021). Network structure of depression and anxiety symptoms in Chinese female nursing students. BMC Psychiatry 21:279. doi: 10.1186/s12888-021-03276-1

PubMed Abstract | Crossref Full Text | Google Scholar

Russell, S. J., and Norvig, P. (2010). Artificial intelligence: A modern approach, 3rd Edn. Upper Saddle River, NJ: Prentice Hall.

Google Scholar

Schultz, D. H., Ito, T., Solomyak, L. I., Chen, R. H., Mill, R. D., Anticevic, A., et al. (2019). Global connectivity of the fronto-parietal cognitive control network is related to depression symptoms in the general population. Netw. Neurosci. 3, 107–123. doi: 10.1162/netn_a_00056

PubMed Abstract | Crossref Full Text | Google Scholar

Scutari, M. (2009). Learning bayesian networks with the bnlearn R package. arXiv [Preprint] doi: 10.48550/arXiv.0908.3817

Crossref Full Text | Google Scholar

Spitzer, R. L., Kroenke, K., Williams, J. B., and Löwe, B. (2006). A brief measure for assessing generalized anxiety disorder: The GAD-7. Arch. Intern. Med. 166, 1092–1097. doi: 10.1001/archinte.166.10.1092

PubMed Abstract | Crossref Full Text | Google Scholar

Spruyt, B., Van Droogenbroeck, F., and Kavadias, D. (2015). Educational tracking and sense of futility: A matter of stigma consciousness? Oxf. Rev. Educ. 41, 747–765. doi: 10.1080/03054985.2015.1117971

Crossref Full Text | Google Scholar

Sun, C., Hu, B., Zhang, P., Wu, T., Wang, L., Xue, W., et al. (2025). Negative emotions among Chinese primary and middle school students and association with screen time: A network analysis. J. Affect. Disord. 370, 373–380. doi: 10.1016/j.jad.2024.11.030

PubMed Abstract | Crossref Full Text | Google Scholar

Sun, C., Zhu, Z., Zhang, P., Wang, L., Zhang, Q., Guo, Y., et al. (2024). Exploring the interconnections of anxiety, depression, sleep problems and health-promoting lifestyles among Chinese university students: A comprehensive network approach. Front. Psychiatry 15:1402680. doi: 10.3389/fpsyt.2024.1402680

PubMed Abstract | Crossref Full Text | Google Scholar

Sun, J., Liang, K., Chi, X., and Chen, S. (2021). Psychometric properties of the generalized anxiety disorder scale-7 item (GAD-7) in a large sample of chinese adolescents. Healthcare 9:1709. doi: 10.3390/healthcare9121709

PubMed Abstract | Crossref Full Text | Google Scholar

van Borkulo, C. D., van Bork, R., Boschloo, L., Kossakowski, J. J., Tio, P., Schoevers, R. A., et al. (2022). Comparing network structures on three aspects: A permutation test. Psychol. Methods 28, 1273–1285. doi: 10.1037/met0000476

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, W., Bian, Q., Zhao, Y., Li, X., Wang, W., Du, J., et al. (2014). Reliability and validity of the Chinese version of the Patient Health Questionnaire (PHQ-9) in the general population. Gen. Hosp. Psychiatry 36, 539–544. doi: 10.1016/j.genhosppsych.2014.05.021

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y., Hu, Z., Feng, Y., Wilson, A., and Chen, R. (2020). Changes in network centrality of psychopathology symptoms between the COVID-19 outbreak and after peak. Mol. Psychiatry 25, 3140–3149. doi: 10.1038/s41380-020-00881-6

PubMed Abstract | Crossref Full Text | Google Scholar

Webb, L. (2021). COVID-19 lockdown: A perfect storm for older people’s mental health. J. Psychiatr. Ment. Health Nurs. 28, 300–300. doi: 10.1111/jpm.12644

PubMed Abstract | Crossref Full Text | Google Scholar

Williams, D. R., and Rast, P. (2020). Back to the basics: Rethinking partial correlation network methodology. Br. J. Math. Stat. Psychol. 73, 187–212. doi: 10.1111/bmsp.12173

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, X., Zhang, Y., and Yu, G. (2022). Prevalence of mental health problems among senior high school students in mainland of China from 2010 to 2020: A meta-analysis. Adv. Psychol. Sci. 30, 978–990. doi: 10.3724/SP.J.1042.2022.00978

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, X. D., Yu, J. C., Wu, Q. F., Chen, J. Y., Wang, Y. C., Yan, D., et al. (2017). [The relationship among depression, anxiety, stress and addictive substance use behavior in 5 935 secondary vocational students]. Zhonghua Yu Fang Yi Xue Za Zhi 51, 226–231. doi: 10.3760/cma.j.issn.0253-9624.2017.03.007

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, Y., Cao, Z., Yang, M., Xi, X., Guo, Y., Fang, M., et al. (2017). Comorbid generalized anxiety disorder and its association with quality of life in patients with major depressive disorder. Sci. Rep. 7:40511. doi: 10.1038/srep40511

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: vocational high school students, depression, anxiety, gender differences, network analysis, Bayesian networks

Citation: He Y, Yang D, Liang K, Zhang Y, Zhang H and Chi X (2025) Exploring the co-occurrence of depression and anxiety symptoms among Chinese vocational high school adolescents through a network approach. Front. Psychol. 16:1572011. doi: 10.3389/fpsyg.2025.1572011

Received: 06 February 2025; Accepted: 16 July 2025;
Published: 11 August 2025.

Edited by:

Małgorzata Stępień-Nycz, Jagiellonian University, Poland

Reviewed by:

Marika Orlandi, Neurological Institute Foundation Casimiro Mondino (IRCCS), Italy
Peijia Zhang, Henan Provincial People’s Hospital, China

Copyright © 2025 He, Yang, Liang, Zhang, Zhang and Chi. 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: Xinli Chi, eGlubGljaGlAMTI2LmNvbQ==; Hongchuan Zhang, aG9uZ2NodWFuLnpoYW5nQGN1ZmUuZWR1LmNu

These authors have contributed equally to this work and share first authorship

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.