Hubness of strategic planning and sociality influences depressive mood and anxiety in College Population

Depressive mood and anxiety can reduce cognitive performance. Conversely, the presence of a biased cognitive tendency may serve as a trigger for depressive mood-anxiety. Previous studies have largely focused on group-wise correlations between clinical-neurocognitive variables. Using network analyses for intra-individual covariance, we sought to decipher the most influential clinical-neurocognitive hub in the differential severity of depressive-anxiety symptoms in a college population. Ninety college students were evaluated for depressive-anxiety symptoms, Minnesota multiphasic personality inventory-2(MMPI-2), and neuro-cognition. Weighted and undirected version of the intra-individual covariance networks, comprised of 18 clinical-neurocognitive variables satisfied small-worldness and modular organization in the sparsity range of K = 0.20–0.21. Furthermore, betweenness centrality of perseverative error for the Wisconsin card sorting test was reduced in more depressive individuals; higher anxiety was related to the increased betweenness centrality of MMPI-2 clinical scale 0(Si). Elevated edge-betweenness centrality of covariance between the MMPI-2 clinical scale 7(Pt) versus commission error of the continuous performance test predicted more anxiety higher than depressive mood. With intra-individual covariance network of clinical-neurocognitive variables, this study demonstrated critical drivers of depressive mood[attenuated influence of strategic planning] or anxiety[domination of social introversion/extroversion, in addition to the influence of compulsivity-impulsivity covariance as a shortcut component among various clinical-neurocognitive features].

Global network characteristics. In the optimal sparsity range of K = 0.20-0.21, for more than 95% of the participants, the weighted and undirected intra-individual covariance networks consisting of 18 clinical-neurocognitive variables satisfied the (1) small-world organization (sigma > 1), (2) modular organization (Q > 0.3; K = 0.10-0. 21), and showed (3) network connectedness [more than 80% (>15) of nodes from a total of 18 nodes were connected to at least one other nodes; K = 0.20-0.21]. Group-averaged mean and standard deviation for a total of four global network characteristics including normalized clustering coefficient (gamma), normalized characteristic path length (lambda), small-worldness (sigma) and modularity (Q) values across the searched sparsity range of K = 0.10-0.30 are illustrated in Fig. 1. T-score of MMPI-2 clinical scales: 1(Hs) 55. 8 11.4 T-score of MMPI-2 clinical scales: 2(D) 66  Betweenness centrality vs. severity of depressive mood or anxiety. The distribution of the rank-transformed betweenness centrality values per participant averaged in the K = 0.20-0.21 is shown in Fig. 2, drawn using distributionPlot.m (https://kr.mathworks.com/ matlabcentral/fileexchange/23661-violin -plots-for-plotting-multiple-distributions-distributionplot-m-) within Matlab R2014b software. When the rank-transformed betweenness centrality value is low, the original betweenness centrality value of a given node [clinical-neurocognitive measure] would be higher than other nodes in the intra-individual covariance network of that participant. Among these nodes, rank-transformed betweenness centrality of Wisconsin Card Sorting Test (WCST) perseverative error demonstrated statistically significant correlation with the total score of Beck Depression Inventory (BDI; Spearman's rho = 0.291, p = 0.005; Fig. 3(A)). Rank-transformed betweenness centrality of MMPI-2 clinical scale Si(0) demonstrated a statistically significant relationship with total score of Beck Anxiety Inventory (BAI; Spearman's rho = −0.321, p = 0.002; Fig. 3(B)). Otherwise, no statistically significant correlation was found (all p > 0.01).
Edge betweenness centrality vs. depressive mood-anxiety dynamic. The most influential top 24 edges [top 20% ranked edges across K = 0.20-0.21 in more than 45% of participants] are shown in Fig. 4 , most of the clinical-neurocognitive variables shared two or more influential edges with other variables. Of note, node MMPI-2: 4(Pd) had seven influential edges, including BAI total score, MMPI-2: 1(Hs) and 6(Pa), omission error during Continuous Performance Test (CPT), reaction time of Trail Making Test part A (TMT-A), perseverative error of WCST, and forward span score of digit span test. Next, node MMPI-2: 7(Pt) was connected with a total of five other clinical-neurocognitive variables including MMPI-2: 9(Ma), commission error and omission error of CPT, in addition to forward and backward span of the digit span test. Another node, named forward span of digit span test, shared five influential edges only with MMPI-2 clinical scales [4(Pd), 6(Pa), 7(Pt), and 0(Si)] and BAI total score.  Furthermore, among the 24 influential edges [reflecting the covariance between two different variables], there were eight influential edges mediating nodes [MMPI-2: 4(Pd) and 7(Pt)] that grasped more than four influential edges into another domain [connector between neurocognitive and clinical nodes]. Specifically, an edge reflecting the rank-transformed edge betweenness centrality of covariance between the MMPI-2: 7(Pt) and commission error of CPT demonstrated statistically significant correlation with difference between the 'z-score transformed BDI total score' and 'z-score transformed BAI total score' (Spearman's rho = 0.451, p = 8.2 × 10 -6 ; Fig. 5).

Discussion
To the best of the authors' knowledge, this is the first reported study to explore the intra-individual covariance-based hierarchical relationship among the clinical [BDI, BAI, MMPI-2 clinical scale]-neurocognitive  In the sparsity range of K = 0.20-0.21, weighted and undirected version of intra-individual covariance network of clinical-neurocognitive variables satisfied small-worldness and modular organization. Furthermore, using the regional network characteristics of betweenness centrality [=importance of specific node or edge as a mediator comprising the shortcut pathway among various nodes], this study uncovered a reduced influence of strategic planning-cognitive flexibility [=WCST perseverative error] on other clinical-neurocognitive characteristics in more depressive individuals; reports of higher anxiety levels were related to increased control of sociality trait [=MMPI-2 0(Si)] in the intra-individual covariance network. Furthermore, tight coupling of obsession-compulsivity [=MMPI-2 7(Pt)] with impulsivity [=CPT commission error] to mediate the covariance between other clinical-neurocognitive variables as a shortcut edge could predict reports of more anxiety symptoms than complaints of depressive mood.
Strategic planning (WCST perseverative error) against depressive mood. In this study, participants highly-ranked for the betweenness centrality of WCST perseverative error reported relatively mild depressive mood in terms of BDI total score (Fig. 2(A)). Appropriate strategic planning and flexible cognitive adaptation of hypotheses about current situation in response to the changing feedback would be crucial for the preservation of self-efficacy and timely performance in academic/occupational tasks [33][34][35] . Conversely, previous suffering or current experience of depressive symptoms has been found to coexist with executive functioning impairment [36][37][38] . When clinical-neurocognitive characteristics are primarily governed by strategic planning and cognitive flexibility [highly-ranked betweenness centrality of WCST perseverative error], stress resilience can be strengthened and vulnerability to psychopathology reduced [39][40][41] . With few therapeutic regimens proven to be effective for ameliorating depression-related executive dysfunction 42,43 , clinicians may inform depressive individuals about the enduring nature of cognitive vulnerability for depression while also referring to the therapeutic potential of cognitive training 44 or theta-burst transcranial magnetic stimulation 45 , among other treatments.
Influential sociality (MMPI-2 clinical scale Si(0)) and higher anxiety. This study demonstrated elevated betweenness centrality value (=higher rank-transformed centrality) of MMPI-2 clinical scale 0(Si) for participants with higher anxiety (Fig. 2(B)). This result is in line with the Heeren and McNally (2016), in which social introversion-related behavioural tendency such as avoidance and fear of social situations were highly ranked for centrality measurement in the inter-individual psychopathology network for social anxiety disorder. Social interaction style, as reflected in tasks such as following of others' eye gaze direction 46 or recognition of facial emotion 47,48 , greatly affects an individual's susceptibility to anxiety 49 . Moreover, when sociality factor measured using MMPI-2 0(Si) prevails in the hierarchical covariance network of clinical-neurocognitive features, this state might be likened to the trait of sensory processing sensitivity 50 specifically biased towards nonverbal social cues 51 . Given the consistency of sociality trait within individuals across time, psychoeducation for enduring trait anxiety combined with mindfulness-based stress reduction training could be helpful in reducing and controlling anxiety symptoms 52,53 . Tight compulsivity-impulsivity covariance as a shortcut to higher anxiety. In this study, participants who reported more severe levels of anxiety beyond the intensity of depressive mood, revealed highly ranked edge betweenness centrality for covariance between MMPI-2 clinical scale 7(Pt) and CPT commission error (Fig. 4). Although poor impulse control might be shown in a sub-population with a genetic loading for depression by mediation of the serotonin transporter 54 , depressive mood per se did not affect CPT performance 55 . Concurrent escalation of compulsivity and impulsivity have been related to substance use disorder in a college population 56 and suggests the highest clinical severity in patients with obsessive-compulsive disorder 57 . Importance of the compulsivity-impulsivity interplay, reflected as increased edge-betweenness centrality of 'MMPI-2: 7(Pt)-CPT: commission error' edge for participants who reported more severe anxiety beyond the intensity of depressive mood, has been also shown valid for anxiety disorders and attention-deficit/hyperactivity disorder 56 . When compulsivity or impulsivity becomes attenuated with symptom-specific cognitive behavioural therapy or brain circuitry-targeted repetitive transcranial magnetic stimulation 58 , the influence of compulsivity-impulsivity covariance as a shortcut mediator among the various clinical-neurocognitive variables might be reduced.
Limitations. This study has some limitations. First, participants did not undergo DSM-derived evaluation procedures using SCID-I or SCID-NP. Second, as participants were measured for their psychological characteristics at only one time-point, intra-individual comparison of covariance network characteristics varying with changed severity of depressive mood and/or anxiety were not possible. However, regardless of the psychiatric diagnosis, by deciphering the key driver of hierarchical dynamics among the clinical-neurocognitive variables per individual, this study has provided clinicians with a specific point of leverage for non-pharmacological treatment as an initial therapy.

Conclusions
Applying the format of an intra-individual covariance network into the hierarchical dynamic of clinical-neurocognitive variables per individual, this study demonstrated critical drivers of depressive mood [attenuated influence of strategic planning, measured using the WSCT perseverative error] or anxiety [domination of social introversion/extroversion reflected in the MMPI-2: 0(Si) in addition to the critical influence of covariance between compulsivity (MMPI-2: 7(Pt)) and impulsivity (CPT commission error) as a shortcut component among various clinical-neurocognitive features] in undergraduate and graduate students not receiving psychotropic medication. Further studies to explore intra-individual changes of hierarchical dynamic among the clinical-neurocognitive variables across multiple time-points would allow examination of the value of these components as effective points of therapeutic intervention.
Moreover, diverse domains of nonverbal neurocognitive performance were assessed with the omission error (the number of failures in responding to a target stimulus; related to inattentiveness) and commission error(responses to a non-target stimulus; measuring impulsivity) of CPT 63 , the reaction times of TMT-A (psychomotor speed) and TMT-B (visuomotor sequencing and set shifting)] 64 , the perseverative error of WCST (the number of erroneous sorting events that occurred when participant used the same principle applied during their previous sorts; concept formation, set-shifting and strategic planning) [65][66][67] and digit span from the Korean version of the Wechsler Adult Intelligence Scale [DS; forward span (sustained attention) and backward span (working memory)] 68 . The Institutional Review Board at Seoul National University Hospital approved the current study. Written informed consent was obtained from all subjects after the procedures had been fully explained, and all methods were performed in accordance with the relevant guidelines and regulations.

Construction of the intra-individual covariance network.
We adapted a total of 18 clinicalneurocognitive variables into an intra-individual covariance network. First, group-wise (n = 90) mean and standard deviation per variable was calculated 27,28 . To ensure that positive z-scores indicated a better condition, all of the z-values other than forward or backward scores of the digit span test were multiplied by -1. Second, all scores per participant were z-score transformed using the group-wise mean and standard deviation. Third, the difference (=covariance) between two z-transformed variables was converted into an edge of intra-individual covariance network of clinical-neurocognitive features using the formula ' ' , in which z(i, t) and z(j, t) denote z-transformed values of specific clinical-neurocognitive variables for participant t, to finally construct a total of 90 intra-individual covariance networks of clinical-neurocognitive features 30,31 . Network analyses: Global network characteristics. Global-regional characteristics of the intra-individual covariance networks of clinical-neurocognitive features (refer to the previous section) were explored with graph theory approach using the Brain Connectivity Toolbox (https://www.nitrc.org/projects/ bct/) in Matlab R2014b (https://kr.mathworks.com). First, each network matrix was trimmed to survive only top 10-30% ranked edges for edge strength [using 'threshold_proportional.m 32 ']. Second, to examine the presence of and to clearly identify the range of sparsity [K; the percentage of edge presence compared to the possible number of edges in a network], in which small-worldness and modular organization of given intra-individual covariance networks satisfied, four global network characteristics were calculated as below: (a) normalized clustering coefficient [gamma; was measured using 'clustering_coef_wu.m 69 , per node and were averaged over a total of 18 nodes, and finally were normalized using the same variable calculated from the 10,000 degree-constrained random networks produced from the original network with 'randmio_und.m 70  'distance_wei.m 32 ' , the characteristic path length among nodes in a given network were calculated using 'charpath. m 32 ' , and finally were normalized using the same variable calculated from the 10,000 random network alike in case of (a)], (C) small-worldness [sigma = gamma/lambda] 32 , and (d) modularity [Q; a statistic that quantifies the degree to which the network may be subdivided into more clearly delineated communities, averaged over 500 runs of trial using 'modularity_und.m 71,72 ' considering the heuristic nature of the algorithm used].
As a result, in the sparsity range of K = 0.20-0.21, more than 95% of intra-individual covariance networks satisfied the (1) small-world organization (sigma > 1), (2) modular organization (Q > 0.3; K = 0.10-0.21), and showed (3) network connectedness [more than 80% ( > 15) of nodes from a total of 18 nodes were connected to at least one other node; K = 0.20-0.21] 73 . Therefore, this sparsity range of K = 0.20-0.21 became the focus of the next exploration for regional network characteristics.
Network analyses: Regional network characteristics. Using the connection-length matrix, an inverse of original intra-individual covariance network matrix (K = 0.20-0.21), 'betweenness centrality [indicates the frequency at which a node is located in the path of the shortcut that connects two different nodes in a network; calculated using 'betweenness_wei.m'] 74 ' per node was calculated. In a scale-free network, distribution of the centrality value does not follow normal distribution. Thus, prior to correlation analyses with clinical characteristics (BDI total score and BAI total score), these betweenness centrality values were rank-transformed per participant using the 'tiedrank.m' function of Matlab R2014b.
Additionally, 'edge betweenness centrality [the fraction of a given edge being a component of shortest paths connecting different nodes in a network; measured using 'edge_betweenness_wei.m'] 74 , values were also gathered from the connection-length version of the original intra-individual covariance network matrix. After having been rank-transformed and being averaged over the sparsity range of K = 0.20-0.21 per participant, the most influential top 24 edges [top 20% ranked edges in more than 45% of participants] were retrieved 75,76 . Statistical analyses. Descriptive statistics of mean and standard deviation values were calculated using the 'mean.m' and 'std2.m' functions of Matlab R2014b, respectively. To explore the correlation between the rank-transformed betweenness centrality values and clinical symptom severity [BDI total score for depressive mood; BAI total score for anxiety], Spearman's correlation coefficients were calculated [using the 'spearman' option of the 'corr.m' function]; considering the exploratory nature of the analysis, the statistical threshold of significance was set at p < 0.01. Furthermore, to challenge the role of influential edges as a reflection of clinical symptom dynamics [reported severity of depressive mood versus anxiety] per participant, Spearman's correlation coefficients between the rank-transformed edge betweenness centrality values [the most influential top 24 edges; top 20% ranked edges in more than 45% of participants] versus difference between 'z-transformed BDI score' and 'z-transformed BAI score' were assessed (p < 0.002 = 0.05/24). Data Availability Statement. The authors will make materials, data and associated protocols promptly available to readers without undue qualifications in material transfer agreements.