Network analysis of depressive and anxiety symptoms in adolescents during and after the COVID-19 outbreak peak.

BACKGROUND: This study examined the extent to which the network structure of anxiety and depression among adolescents identified during the peak of the COVID-19 pandemic could be cross-validated in a sample of adolescents assessed after the COVID-19 peak.
METHODS: Two cross-sectional surveys were conducted between February 20 and 27, 2020 and between April 11 and 19, 2020, respectively. Depressive and anxiety symptoms were assessed using the 20-item Center for Epidemiological Studies-Depression and 7-item Generalized Anxiety Disorder, respectively. Anxiety-depression networks of the first and second assessments were estimated separately using a sparse Graphical Gaussian Model combined with the graphical least absolute shrinkage and selection operator method. A Network Comparison Test was conducted to assess differences between the two networks.
RESULTS: The most central symptoms in the first and second survey networks were Depressed affect and Nervousness. Compared with connections in the first survey network, connections in the second survey network analysis between Relax-Nervousness-Depressed affect-Interpersonal problems (diff, contrast: second survey-first survey. diff=-0.04, P = 0.04; diff=-0.03, P = 0.03; diff=-0.03, P = 0.04), and Irritability-Somatic complaints (diff=-0.04, P = 0.02) were weaker while connections of Somatic complaints-Nervousness (diff=0.05, P<0.001), Somatic complaints-Depressed affect (diff=0.03, P = 0.009), and Irritability-Control worry-Restlessness (diff=0.02, P = 0.03; diff=0.05, P = 0.02) were stronger.
CONCLUSIONS: Depressed affect emerged as a robust central symptom and bridge symptom across Anxiety-Depression networks. Considering the negative impact of depression and anxiety on daily life, timely interventions targeting depressed affect should be implemented to reduce the co-occurrence of anxious and depressive symptoms among adolescents during the COVID-19 pandemic.

Introduction

The coronavirus disease 2019 (COVID-19) pandemic is a global health emergency (Eurosurveillance Editorial Team., 2020; Mattila et al., 2021). As of early 2021, COVID-19 has spread to more than 200 countries and territories, with more than 100 million confirmed cases (WorldoMeters, 2021). To control rapid transmission, a range of strategies has been implemented and includes mass lockdowns, home quarantines and school closures, all of which could have a negative impact on mental health of certain subpopulations including adolescents and their families (Ghosh, 2020; Golberstein et al., 2020; Guessoum et al., 2020; Zhu et al., 2020). Previous studies found depressive and anxiety symptoms were common in adolescents during the COVID-19 pandemic (Qi et al., 2020; Zhou et al., 2020a, 2020b). Although they are found in separate diagnostic categories (Zbozinek et al., 2012), depression and anxiety often occur together, particularly in children and adolescents (Melton et al., 2016), in part, because they have some causes in common. Considering the co-occurrence of depression and anxiety in different populations during the COVID-19 pandemic as well as exploring the pattern and relationship between anxious and depressive symptoms during different phases of pandemics may aid in explaining their negative impact on daily life and functional outcomes, developing effective preventive measures, and facilitating timely treatments.

In traditional approaches to understanding the psychopathology of depression and anxiety, symptoms are secondary to an underlying common cause (Borsboom, 2008; Schmittmann et al., 2013). For instance, depression causes dysphoric feelings, insomnia symptoms and impaired cognition in a manner similar to infections as a cause of fever. However, many studies have found that individual symptoms of a particular psychiatric diagnosis may be related to different risk factors or causes (Fried, 2015; Fried and Nesse, 2014; Kapur et al., 2012; Mamun et al., 2020; Wasserman et al., 2021). Furthermore, differences between individual symptoms and their associations with typical symptoms of other syndromes are not usually emphasized in traditional theories of psychopathology (Beard et al., 2016; Contreras et al., 2019).

Network analysis is a new approach to investigating complex, dynamic relationships between individual psychiatric symptoms (Fisher et al., 2017). This method is based on the premise that a particular psychiatric disorder is an interacting cluster of symptoms (Jones et al., 2019), and emphasis on understanding the strength and nature of associations among symptoms (Beard et al., 2016; Rouquette et al., 2018). In network analysis, nodes reflect psychiatric symptoms and edges between nodes reflect relationships between symptoms, including the activation spread from one symptom to other symptoms through the network (Borsboom, 2017; Rouquette et al., 2018). Nodes can also function as bridge symptoms that transfer symptom activation from one disorder to another disorder. Therefore, network analysis has utility in elucidating connection patterns between individual psychiatric symptoms and also between psychiatric disorders (Borsboom and Cramer, 2013).

Several past network analysis studies have examined the comorbidity of depression and anxiety. For example, one such study on generalized anxiety disorder and major depressive disorder found that fatigue, restlessness and sleep disturbances were key symptoms associated with the comorbid depression and anxiety (Cramer et al., 2010). Another study that examined the association between anxiety as measured by the 7-item Generalized Anxiety Disorder Scale (GAD-7) and depression as measured by the Patient Health Questionnaire-9 (PHQ-9), found that sad mood and worry were the most prominent central symptoms (i.e., those are most highly connected to other symptoms) in a large psychiatric sample (Beard et al., 2016).

Some network analysis studies have been conducted during the COVID-19 pandemic. A study examining changes of depressive and anxiety symptoms between the COVID-19 outbreak peak and its aftermath among Chinese adults (Wang et al., 2020) indicated symptoms of restlessness and an inability to relax had high centrality during the outbreak but gradually remitted after the peak. Furthermore, symptoms of irritability and loss of energy were bridge symptoms (i.e., symptoms that connect one psychiatric syndrome/disorder to another and may be important as targets of treatments intended to reduce or prevent comorbid problems) after the pandemic peak. However, because clinical features of depression and anxiety can be quite different between adolescents and adults (Nardi et al., 2013; Pyramid Healthcare, 2018; Rice et al., 2019), there is no clear basis for assuming findings of network analyses in adults during the COVID-19 pandemic are also applicable to adolescents who experience these symptoms. Furthermore, it is not clear whether network findings based on only one assessment are replicable or can be generalized over time or within other samples drawn from the same population. Toward addressing this limitation, cross-validation procedures designed to evaluate the stability of initial network models across two samples are one strategy that can be incorporated to clarify more versus less stable features of models generated in network analysis studies.

In sum, although network analysis has aided in elucidating links between psychiatric symptoms in adult samples, little is known about trajectories of individual depressive and anxiety symptoms among adolescents. Therefore, we conducted this study to characterize depressive and anxiety symptom clusters using network analysis and explored the extent to which a network model generated in one large adolescent sample during the COVID-19 outbreak peak could be replicated within another large adolescent sample shortly after the peak.

Methods

Settings and participants

Two national, cross-sectional surveys were undertaken by the National Clinical Research Center for Mental Disorders, China. The first survey was conducted between February 20 and 27 (peak of the COVID-19 outbreak) during which adolescents were required to stay at home due to nationwide lockdown measures (Peng et al., 2020; Sanchez-Caballero et al., 2020). The second survey was conducted between April 11 and 19, 2020 after the COVID-19 outbreak peak had remitted and during which schools in some provinces in China reopened (Xinhuanet, 2020). To reduce the risk of contagion, face-to-face interviews were not adopted. Following other studies (Liu et al., 2020b; Wang et al., 2020), data were collected through the WeChat-based QuestionnaireStar application using a snowball sampling method. WeChat is a widely used social communication application with 1.2 billion users in China, and is employed in most secondary schools nationwide for the online teaching. To be eligible for the research, participants were to be: (1) secondary school students aged between 11 and 20 years and (2) living in China during the COVID-19 outbreak. All participants and their guardians were required to provide written informed consent. This study was approved by the Medical Ethical Committee of Beijing Anding Hospital, Capital Medical University, China.

Measurements

Depressive symptoms were assessed with the self-report Center for Epidemiological Studies Depression Scale (CES-D) (Radloff, 1977); the CES-D has 20 items covering four clusters of symptoms: depressed affect (items 3, 6, 9, 10, 14, 17 and 18), positive affect (items 4, 8, 12 and 16), somatic complaints (items 1, 2, 5, 7, 11, 13 and 20), and interpersonal problems (items 15 and 19) (Nguyen et al., 2004; Radloff, 1977). Each item is scored from 0 (not at all) to 3 (a lot); total scores range from 0 to 60, with higher scores indicating more severe depressive symptoms. The Chinese version of the CES-D has satisfactory psychometric properties and has been widely used in adolescents (Guo et al., 2016; Sun et al., 2017; Zhou et al., 2020a).

Anxiety symptoms were measured using the 7-item self-report Generalized Anxiety Disorder Scale (GAD-7) (Spitzer et al., 2006). Each item is scored from 0 (not at all) to 3 (nearly every day). Total scores range from 0 to 21 and higher total scores indicate more severe anxiety symptoms. The Chinese version of GAD-7 has been validated with satisfactory psychometric properties in Chinese samples (He et al., 2010; Tong et al., 2016). Item content of the two scales is listed in Table S1.

Network estimation

Network analyses were conducted using bootnet (Epskamp et al., 2018a) and qgraph (Epskamp et al., 2012) packages in the R program (version 3.6.3). In network parlance, symptoms were nodes, while correlations between individual symptoms were edges (Beard et al., 2016; Wang et al., 2020). In this study, 11 nodes/symptoms were included: (1) 4 clusters of depressive symptoms (depressive nodes) as measured by the CES-D (i.e., Depressed, Positive, Somatic, and Interpersonal) and (2) 7 anxiety symptoms (anxiety nodes) as measured by the GAD-7 (i.e., Nervousness, Control Worry, Worry A Lot, Relax, Restlessness, Irritability, and Afraid). Pearson correlation analyses were used to estimate strengths of association between nodes, with thicker edges indicating stronger relationships.

To examine relationships between symptoms within networks during the COVID-19 outbreak, network models from the two surveys were estimated separately using sparse Graphical Gaussian Models (GGM) combined with a graphical least absolute shrinkage and selection operator (LASSO) method (Friedman et al., 2008); model selection was based on the Extended Bayesian Information Criterion (EBIC) (Chen and Chen, 2008). The GGM method was applied to construct network models. However, the approach can produce a large number of spurious edges due to latent variables (Epskamp et al., 2018b; Montazeri et al., 2020). To reduce the number of spurious edges and improve interpretability of results, network models were regularized using LASSO, a well-established algorithm for regularization that reduces weak associations to zero by removing them from models as potentially “false positive” edges (Epskamp et al., 2018a; Friedman et al., 2008; Heeren et al., 2018). EBIC is a widely used approach to model selection based on a tuning parameter (Epskamp, 2016; van Borkulo et al., 2015) that is usually set between 0 and 0.5. If the tuning parameter is close to 0.5, EBIC generates a simple model containing fewer edges; in contrast, if the tuning parameter is close to 0, EBIC generates a model with more edges (Epskamp et al., 2018a; Heeren et al., 2018). Following previous studies (Beard et al., 2016; Wang et al., 2020), the ‘estimateNetwork’ function in bootnet (Epskamp et al., 2018a) was used, with ‘EBICglasso’ as the default method and 0.5 as the default tuning parameter.

To assess the importance of each node in the network, centrality indices including Strength, Betweenness and Closeness were calculated (Borgatti, 2005; Freeman, 1978; Opsahl et al., 2010; Pan and Liu, 2020), using the ‘centralityPlot’ function in qgraph package (Epskamp et al., 2012) and normalized (z-scored) values were plotted for each node. Strength is the sum of absolute edge weights of all direct connections between a specific node and other nodes, reflecting the importance of an individual symptom in the network. Betweenness represents the frequency of a node that lies on all the shortest paths between other nodes, reflecting the importance of a symptom as a “connecter” to other symptoms. Closeness refers to the reciprocal of the sum of the shortest path distance of a specific node to all other nodes in a network, indicating the strength of a symptom's indirect connections to other symptoms (Dalege et al., 2017; Wang et al., 2020). Furthermore, to identify individual symptoms that act as pathways in links between depressive and anxious symptoms, bridge centrality indices (e.g., bridge strength, bridge closeness, and bridge betweenness) were measured using the ‘bridge’ function in networktools (Jones, 2021). Bridge centrality measures the importance of a node in linking its network to other networks (Jones et al., 2019; Wang et al., 2020).

Estimation of network stability

Following previous studies (Beard et al., 2016; Belvederi Murri et al., 2020; Wang et al., 2020), the accuracy and stability of networks were tested using R package bootnet (Epskamp et al., 2018a). First, to estimate the accuracy of edge weights, non-parametric bootstrapping (1000 replicates, 8 cores) was performed to compute 95% confidence intervals (CI) of edge values. Second, bootstrapped difference tests were used to determine significant differences between edge weights and node strengths. Finally, to determine the stability of centrality indices, case-dropping subset bootstrapping (1000 replicates, 8 cores) was performed to compute correlation stability coefficients (CS). A host of correlation values was calculated between centrality indices in the original sample and subset centrality based on different subsets of original samples (e.g., 95% of the sample, 80%, 70%, …, 25%) (Epskamp et al., 2018a; Hevey, 2018). CS represents the maximum proportion of cases that can be dropped from the original sample; for instance, with a 95% probability, correlation coefficients between centrality indices in the original networks and centrality indices in case-subset networks should reach at least 0.7 (default) (Epskamp et al., 2018a). The CS should not be lower than 0.25, and is recommended to be above 0.5 (Epskamp et al., 2018a).

Comparison of networks between the first and the second surveys

To assess differences between first versus second surveys, a Network Comparison Test (NCT) was conducted via the R package NetworkComparisonTest (van Borkulo et al., 2017). NCT is a two-tailed permutation method to test differences between two networks based on several invariance measures (i.e., global strength, network structure, and edge strength) (Dalege et al., 2017; van Borkulo et al., 2017). In this procedure, NCTs were run 1000 times and mean differences in global strength (the sum of absolute edge weights of the two networks were compared (Opsahl et al., 2010)). In addition, significant differences between two networks were estimated separately using a P<0.05 (two-tailed) level of significance.

Sensitivity analyses

In order to test the robustness of results across samples, a sensitivity analysis was conducted by repeating the network analyses after excluding data from those who participated in both surveys (n = 6719) from the second survey sample (i.e., group sampled after the COVID-19 peak had remitted).

Results

Participants

A total of 20,158 adolescents from mainland China was included in the final analysis; 9,553 adolescents participated in the first survey, with a mean age of 15.11 1.94 years. In the second survey, 10,605 adolescents completed the assessment; their mean age was 15.44 1.74 years. Table 1 shows the basic demographic characteristics of participants in each sample.

Table 1

Demographic characteristics of survey samples.

	First Survey (N = 9553)		Second Survey (N = 10,605)		Second Survey1 (N = 3886)
Age (mean, SD)	15.11	1.94	15.44	1.74	15.79	1.60
Female (n,%)	4977	52.10	5708	53.80	2132	54.9
Junior secondary school	5459	57.10	4656	43.90	1262	32.5
CES-D (mean, SD)
Depressed affect	4.62	4.53	5.76	5.01	6.42	5.02
Positive affect	3.00	2.85	3.49	3.01	3.81	3.01
Somatic complaints	5.43	4.27	6.89	4.68	7.65	4.58
Interpersonal problems	1.01	1.48	1.23	1.56	1.41	1.60
GAD total score (mean, SD)	2.28	3.54	3.46	4.34	3.98	4.46

Note: 1 Based on a sensitivity analysis of the independent subgroup that completed only the second survey. CES-D, Center for Epidemiological Studies Depression Scale; GAD-7, 7-item Generalized Anxiety Disorder Scale. SD, standard deviation.

Network structures

As shown in Fig. 1 , the first survey network (during the COVID-19 outbreak peak) of 11 nodes revealed that 48 of 55 edges were estimated to be above zero. Similarly, the second survey network (after the outbreak peak) of 11 nodes indicated 47 of 55 edges were estimated to be above zero. In both the first and second sample networks, the strongest edges between nodes of depression were for Depressed affect and Somatic complaints as well as Depressive affect and Interpersonal problems while the two strongest edges between nodes of anxiety were for (1) Relax and Restlessness, and (2) Nervousness and Control Worry. Bootstrapped confidence intervals for the edges (Figure S1) and correlation matrices (Table S2 and S3) are shown in Supplementary materials.

Fig. 1

Symptom network of depressive and anxiety symptoms among adolescents during different stages of the COVID-19 outbreak. The first survey was conducted between February 20 and 27 (during the COVID-19 outbreak). The second survey was conducted between April 11 and 19 (after the COVID-19 outbreak). In the diagram symptom nodes with stronger connections are closer to each other. The orange nodes denote the CES-D factors; the blue nodes denote the GAD-7 items. The dark green lines represent positive correlations. The edge thickness represents the strength of the association between symptom nodes.

Fig. 2 shows centrality indices of all nodes as measured in the two surveys. Regarding closeness centrality in the first survey network, Nervousness was the strongest symptom connected with other nodes on the shortest path, followed by Irritability and Control Worry. Strength centrality data suggested that Depressed affect was the most central symptom, followed by anxiety symptoms of Relax and Control worry. Analyses of betweenness centrality identified Depressed affect as the node that most frequently lying on the shortest path between other nodes, followed by Nervousness, Restlessness and Irritability.

Fig. 2

Centrality measures of all symptoms within the network at different stages of the COVID-19 outbreak. The figure shows centrality measures (i.e., strength, betweenness, and closeness) of all symptoms within the network (z-scores).

Similarly, in the second survey network, Nervousness was cross-validated as the node with the highest closeness centrality, followed by Irritability and Control Worry. Depressed affect was the node with the highest strength centrality, followed by Control worry and Relax. Replicating the initial network model, Depressed affect also had the highest betweenness centrality, followed by Nervousness; however, unlike the initial network model, Control worry had the highest betweenness centrality. Bridge centrality indices underscored Depressed affect as the node with the highest bridge strength and bridge betweenness while Nervousness had the highest bridge closeness across networks representing adolescent samples assessed during two stages of the COVID-19 outbreak (Figure S2).

Network accuracy and stability

As shown in Fig. 3 , results of case-dropping subset bootstrapping indicated the strength, betweenness, closeness, bridge strength, bridge betweenness, and bridge closeness remained stable under various analysis conditions. CS coefficients for strength, betweenness and closeness were 0.75 in the first survey; 0.75, 0.67 and 0.75, respectively in the second survey. CS coefficients for bridge strength, bridge betweenness and bridge closeness were 0.75, 0.67 and 0.75, respectively in the first survey; 0.75, 0,67 and 0.59, respectively in the second survey.

Fig. 3

Stability of centrality indices by case dropping subset bootstrap. The x-axis represents the percentage of cases of original sample used at each step. The y-axis represents the average of correlations between the centrality indices from the original network and the centrality indices from the networks that were re-estimated after dropping increasing percentages of cases. Each line indicates the correlations of betweenness, strength, closeness, bridge betweenness, bridge strength, and bridge closeness, while areas indicate 95% CI.

Bootstrapped 95% CIs for edges at different stages of the COVID-19 outbreak were narrow and edge values were significantly higher than zero, indicating that the edges were stable (Figure S1). The bootstrapped difference test revealed that most of the edge weights and node strengths had statistically significant differences from one another in individual comparisons (Figure S3 and S4).

Network comparisons

Comparisons of the first versus the second surveys networks based on NCT results showed a significant difference in global strength (S = 0.09, P<0.001), but not network structure (M = 0.05, P = 0.22). Compared with the first survey network, the second survey network showed edge weight decreases for Relax- Nervousness – Depressed affect – Interpersonal problems (diff, contrast: second survey – first survey. diff=−0.04, P = 0.04; diff=−0.03, P = 0.03; diff=−0.03, P = 0.04), and Irritability – Somatic complaints (diff=−0.04, P = 0.02). Conversely, the second survey network revealed comparatively increased edge weights for Somatic complaints - Nervousness (diff=0.05, P<0.001), Somatic complaints – Depressed affect (diff=0.03, P = 0.009), and Irritability – Control worry – Restlessness (diff=0.02, P = 0.03; diff=0.05, P = 0.02) (Fig. 4 ).

Fig. 4

Significant different edge between the two stages of the COVID-19 outbreak. The blue nodes denote the GAD-7 items, and the orange nodes denote the CES-D factors. Compared with the first survey, the second survey network showed the edge weights of Relax- Nervousness - Depressed affect - Interpersonal problems, and Irritability - Somatic complaint were significant decreased (light green lines). In contrast, Somatic complaint- Nervous, Somatic complaint - Depressed affect, and Irritability - Control worry - Restlessness were significant increased (red lines).

Sensitivity analysis

After excluding adolescents who completed both surveys from the second survey (i.e., those who joined both the first and second surveys), 9,553 adolescents from the first survey and 3,886 adolescents from the second survey were included in the sensitivity analysis. Replicating primary results based on the original sample, the strongest edges between nodes of depression in the second survey network (N = 3,886) were the connections of Depressed affect with Somatic complaints and Interpersonal problems, while the two strongest edges between nodes of anxiety were connections between Relax and Restlessness, and between Nervousness and Control Worry (Figure S5). Depressed affect had the strongest strength centrality and betweenness centrality; Nervousness had the highest closeness centrality (Figure S6). Further, compared with the network based on the first survey, the network based on the independent sample that completed only the second survey showed edge weight decreases for Restlessness – Interpersonal problems – Depressed affect (diff=−0.03, P = 0.008; diff=−0.05, P = 0.001), and Irritability – Somatic complaints (diff=−0.05, P = 0.03). Conversely, the edge weight of Irritability – Control worry (diff=0.06, P = 0.002) increased in the network based on the sample that completed only the second survey (Figure S7).

Discussion

Although numerous network studies have documented relations between psychiatric symptoms within a single sample or time point during the COVID-19 outbreak, it is not clear whether observed models are replicable across two samples drawn from the same population and assessed during different phases of the pandemic. To the best of our knowledge, this was the first network analysis study to investigate the stability of depressive and anxious symptom networks using cross-validation procedures within national samples of adolescents during different stages of the COVID-19 outbreak. In general, network structures were stable across samples assessed during and shortly following the COVID-19 outbreak peak, suggesting that symptom relations observed in the initial network model are similar across Chinese adolescent samples assessed during different pandemic phases.

Depressed affect emerged as the most prominent central and bridge symptom for adolescents across survey networks during and following the COVID-19 outbreak peak in China. These results are consistent with previous network analysis findings showing that depressed affect is also the central node in depression networks of adults (Fried et al., 2016; van Borkulo et al., 2015). Still other network analysis research on depression and anxiety during the pandemic also found that Depressed affect was the central symptom in nursing students (Bai et al., 2021) and the general population (Hoffart et al., 2021). Depressed affect is a hallmark characteristic of depression (Beard et al., 2016) that refers to a cluster of negative mood experiences such as sadness, feeling “blue”, depressed mood, fear, loneliness, crying spells, and perceptions of the self as a failure (Danhauer et al., 2013). Due to the COVID-19 outbreak peak, adolescents in China could not pursue many routine activities outside of the home such as schooling, social communications and physical activities; presumably disrupted routines and isolation increased the likelihood of depression in this group (Courtney et al., 2020; Liu et al., 2020a). Furthermore, extended homestays due to long-term quarantine measures during the pandemic, may have contributed to loneliness, increased stress, and fear that could increase risk for negative mood experiences such as depression and anxiety (Courtney et al., 2020; Guessoum et al., 2020; Qi et al., 2020). Previous studies have found that depressive affect is the common risk factor underlying comorbid emotional disorders including depression and anxiety (Conway et al., 2017; Mineka et al., 1998; Naragon-Gainey et al., 2016; Prenoveau et al., 2010). Such evidence may help to could account for depressed affect as the most prominent symptom linking depressive and anxious symptoms across two samples in this study.

Also underscoring the replicability of observed network model relations across different samples of Chinese adolescents, Depressed affect had its strongest associations with other depression nodes reflecting somatic complaints and interpersonal problems across both surveys during the COVID-19 outbreak. Select past studies have found that both depressed affect and stressful events are associated with somatic symptoms (Bohman et al., 2018; Firat et al., 2018), while other studies have indicated depressed affect is also directly related to interpersonal problems (Majd Ara et al., 2017; Rice et al., 2017; Schwartz-Mette and Smith, 2018). In tandem with more severe depressed affect, adolescents might suffer from more frequent physical complaints, lower self-efficacy and reduced life contentment during the COVID-19 pandemic (Yun et al., 2021). Notwithstanding this evidence from other samples, the pandemic has resulted in increased restrictions related to engaging in routine physical activities and higher than typical levels of social isolation. Hence, these disruptions to daily life brought about by the COVID-19 pandemic may have contributed to more pronounced connections of depressed affect with somatic complaints and interpersonal problems across two different assessments of adolescents in this study.

The symptom of Nervousness presented high closeness and bridge closeness in both survey symptom networks, a finding that suggests Nervousness is more likely to affect and/or be affected by changes in other symptoms. Nervousness and anxiety are common among adolescents due to frequent exposure to negative events and feelings of insecurity (Kılınçel et al., 2021; Qi et al., 2020; Zhou et al., 2020b). This finding aligns with other recent China-based findings (Bai et al., 2021) identifying Nervousness as a bridge symptom between depression and anxiety during the COVID-19 pandemic among nursing students. Furthermore, connections between (1) Nervousness and Control Worry, and (2) Relax and Restlessness were relatively strong in both the first and second survey networks. Interactions between Nervousness and Control worries as well as contradictory experiences of Relaxation and Restlessness might have been strengthened among adolescents due to the widespread uncertainty and unpredictability resulting from the pandemic (Durodie, 2020; Orgilés et al., 2020).

Although the most prominent central symptoms, bridge symptoms and edges of the initial network model were cross-validated across samples and pandemic phases, select findings were not replicated across the samples. Specifically, connections of Relax – Nervousness – Depressed affect – Interpersonal problems, and Irritability – Somatic complaints were weaker in network model of the second sample. In addition, connections of Restlessness– Interpersonal problems– Depressed affect, and Irritability – Somatic complaints had significantly weaker connections in the second survey when the sensitivity analysis was run on the participant subgroup that had completed only the second survey (N = 3,886), As such, linkages of mood, communication and physical problems within anxiety-depression network models may be less stable or generalizable across pandemic phases and/or different groups of Chinese adolescents. It is not clear why these facets of the initial network model were less replicable. One hypothesis may be that changes in pandemic conditions contributed to changes in these relations. To elaborate, during the initial COVID-19 outbreak, many adolescents experienced heavy academic burdens and psychological distress, both of which had a negative impact on their physical and mental health in tandem with increased physical and social isolation (Jiao et al., 2020; Qi et al., 2020; Zhou et al., 2020b). After the peak of the COVID-19 outbreak, lockdown and quarantine measures were gradually lifted in most areas and schools reopened; these shifts toward regular routines may have may have contributed to more attenuated connections between these nodes within the symptom network following the pandemic peak (Fancourt et al., 2020; Lu et al., 2020; Wang et al., 2020).

Other relations from the initial network model that were not cross-validated included stronger connections between nodes of Depressed affect - Somatic complaints - Nervousness, and Irritability – Control worry – Restlessness observed in the network model generated for the second assessment of adolescents. Furthermore, the connection of Irritability – Control worry was significantly stronger for the second survey based on the sensitivity analysis of participants who completed only the second survey. Although it is not clear why these associations were not consistent in the post-peak survey sample, it should be noted that quarantine measures can have a lingering long-term negative impact on experiences of stress. For instance, following the COVID-19 outbreak peak, many adolescents had to re-establish rapport with teachers and classmates, adapt to the transition from online to classroom teaching, cope with heavy academic demands related to the end of an academic year, and live with the uncertainty of possible future waves of COVID-19 (Lee, 2020; Singh et al., 2020; Tasso et al., 2021); these adjustments could contribute to the persistence or increased connectivity between these nodes.

In relation to treatment implications, high overall depression and anxiety levels, highly consistent central symptoms, reliable bridge symptoms and stable connections between select nodes across anxiety-depression network models of two different evaluations of Chinese adolescents underscored the importance of monitoring adolescents’ long-term mental health across different phases of the COVID-19 pandemic as well as the potential utility of targeting particular symptoms, most notably in interventions tailored for the overall population of adolescents. Drawing upon recommendations from other network analysis research (Beard et al., 2016; Fried et al., 2016), the cross-validation of depressed affect both as a central symptom and a bridge symptom across both assessment samples suggests that interventions targeting experiences of sadness, depressed mood, nervousness, loneliness, and physical activity may be particularly beneficial in reducing distress and potential comorbidity of anxiety and depression among adolescents during different pandemic phases. For example, Cognitive Behavioral Therapy (CBT) strategies such as behavioral activation and cognitive restructuring may aid in reducing depressed affect and help to explain the efficacy of very brief CBT as an approach to reducing to depression in this population (Bjorgvinsson et al., 2014).

Strengths of this study included its use of a cross-validation approach to identify the most reliable features of anxiety-depression networks across two assessments of large national samples of adolescents, assessments at different phases of a pandemic, and use of network analysis to elucidate the nature and strength of associations between depressive and anxious symptoms of adolescents. In addition, the sensitivity analysis demonstrated that most of the main findings were robust and stable across independent samples of adolescents evaluated during different phases of the COVID-19 pandemic. The results may be useful in encouraging researchers and clinicians to consider the significance of individual symptoms of depression and anxiety across groups of adolescents drawn from the same population as well as their interrelations while questioning the assumption that all symptoms have an equal impact in the expression of these disturbances.

The main limitations should also be noted. First, data for each sample were cross-sectional in nature so they offered single snapshots of symptom networks for particular groups of adolescents at particular points during a pandemic rather than dynamic changes within a single group over different phases of the pandemic. Although our research helped to highlight anxiety-depression network features that are important across different groups of Chinese adolescents and stages of the COVID-19 pandemic, longitudinal extensions following adolescents over more extended periods are needed to evaluate whether and how the structure of adolescent depression-anxiety symptom networks may evolve over time within a single group of adolescents. Second, because results were based on self-report scales, replications based on interviewer-rated scales would be useful in bolstering the validity of results. Third, although most of the main results were replicated across two national samples, random selection was not employed to recruit participants, in part, because pandemic conditions did not allow for offline recruitment strategies. As such, findings may not be entirely representative of the adolescent population of China. Finally, it is not clear whether results from this survey generalize to adolescents in other countries or non-pandemic eras. This limitation provides another foundation for extensions.

In conclusion, network analyses of central and bridge symptoms of depression and anxiety found relatively stable structures of depressive and anxious symptoms across two evaluations of large national adolescent samples assessed during and after the peak of the COVID-19 outbreak in China. In addition, select relations between symptoms were not cross-validated across the two assessments, suggesting that some features of the network models are less stable and do not extend as well to the overall study population and/or different pandemic phases. The findings suggested that effective interventions targeting depressed affect in particular during both phases of the pandemic could be helpful reducing co-occurring depressive and anxious symptoms among Chinese adolescents.

Role of funding

The study was supported by the (81901368), Funding Support (ZYLX202128); the Beijing Municipal Administration of Hospitals Incubating Program (PX2020072) and the Clinical Medicine Development of special funding support (XMLX202128).

Author contributions

Study design: Rui Liu, Ling Zhang, Yu-Tao Xiang. Data collection, analysis and interpretation: Xu Chen, Han Qi, Rui Liu, Yuan Feng, Teris Cheung, Hui Lei. Drafting of the manuscript: Rui Liu, Han Qi, Zhaohui Su, Yu-Tao Xiang. Critical revision of the manuscript: Todd Jackson. Approval of the final version for publication: all co-authors.

Declaration of Competing Interest

The authors have no conflicts of interest to declare.