Advertisement
Archival Report|Articles in Press

Development of Neural Mechanisms Underlying Threat Processing: Associations With Childhood Social Reticence and Adolescent Anxiety

Open AccessPublished:February 02, 2023DOI:https://doi.org/10.1016/j.bpsgos.2023.01.008

      Abstract

      Background

      Social reticence in early childhood is characterized by shy and anxiously avoidant behavior, and it confers risk for pediatric anxiety disorders later in development. Aberrant threat processing may play a critical role in this association between early reticent behavior and later psychopathology. The goal of this longitudinal study is to characterize developmental trajectories of neural mechanisms underlying threat processing and relate these trajectories to associations between early-childhood social reticence and adolescent anxiety.

      Methods

      In this 16-year longitudinal study, social reticence was assessed from 2 to 7 years of age; anxiety symptoms and neural mechanisms during the dot-probe task were assessed at 10, 13, and 16 years of age. The sample included 144 participants: 71 children provided data at age 10 (43 girls, meanage = 10.62), 85 at age 13 (46 girls, meanage = 13.25), and 74 at age 16 (36 girls, meanage = 16.27).

      Results

      A significant interaction manifested among social reticence, anxiety symptoms, and time, on functional connectivity between the left amygdala and the left dorsolateral prefrontal cortex, voxelwise p < .001, clusterwise familywise error p < .05. Children with high social reticence showed a negative association between amygdala–dorsolateral prefrontal cortex connectivity and anxiety symptoms with age, compared to children with low social reticence, suggesting distinct neurodevelopmental pathways to anxiety.

      Conclusions

      These findings were present across all conditions, suggesting task-general effects in potential threat processing. Additionally, the timing of these neurodevelopmental pathways differed for children with high versus low social reticence, which could affect the timing of effective preventive interventions.

      Keywords

      Early-childhood inhibited behavior is an important risk factor for the development of pediatric anxiety disorders, which typically emerge in adolescence (
      • Clauss J.A.
      • Blackford J.U.
      Behavioral inhibition and risk for developing social anxiety disorder: A meta-analytic study.
      ,
      • Sandstrom A.
      • Uher R.
      • Pavlova B.
      Prospective association between childhood behavioral inhibition and anxiety: A meta-analysis.
      ). However, not all children with inhibited behavior will eventually develop an anxiety disorder. Therefore, discovering factors associated with early behavior and later anxiety is of clinical importance. Previous studies have proposed that threat processing is associated with early-childhood behavior and anxiety (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ,
      • Hardee J.E.
      • Benson B.E.
      • Bar-Haim Y.
      • Mogg K.
      • Bradley B.P.
      • Chen G.
      • et al.
      Patterns of neural connectivity during an attention bias task moderate associations between early childhood temperament and internalizing symptoms in young adulthood.
      ). In this 16-year longitudinal study with three waves of neuroimaging data assessing threat processing in a sample well characterized for early-childhood behavior, we extend this work on neurodevelopmental trajectories of pediatric anxiety disorders.
      Early-childhood inhibited behavior has been defined in several ways (
      • Chronis-Tuscano A.
      • Danko C.M.
      • Rubin K.H.
      • Coplan R.J.
      • Novick D.R.
      Future directions for research on early intervention for young children at risk for social anxiety.
      ,
      • Rubin K.H.
      • Coplan R.J.
      • Bowker J.C.
      Social withdrawal in childhood.
      ). Our laboratory has defined behavioral inhibition as a young child’s response to novelty and uncertainty in a laboratory setting (
      • Fox N.A.
      • Henderson H.A.
      • Rubin K.H.
      • Calkins S.D.
      • Schmidt L.A.
      Continuity and discontinuity of behavioral inhibition and exuberance: Psychophysiological and behavioral influences across the first four years of life.
      ). The child is observed responding to an unfamiliar adult and to novel objects while the child’s caregiver is in the room. Social reticence is defined as shy, anxious, and avoidant behavior observed during interactions with either familiar or unfamiliar same-age peers (
      • Degnan K.A.
      • Almas A.N.
      • Henderson H.A.
      • Hane A.A.
      • Walker O.L.
      • Fox N.A.
      Longitudinal trajectories of social reticence with unfamiliar peers across early childhood.
      ). Social reticence is associated with social competence in middle childhood (
      • Nelson E.E.
      • Jarcho J.M.
      • Guyer A.E.
      Social re-orientation and brain development: An expanded and updated view.
      ) and is a risk factor for anxiety symptoms (
      • Clarkson T.
      • Eaton N.R.
      • Nelson E.E.
      • Fox N.A.
      • Leibenluft E.
      • Pine D.S.
      • et al.
      Early childhood social reticence and neural response to peers in preadolescence predict social anxiety symptoms in midadolescence.
      ). While behavioral inhibition predicts social reticence in the current sample, we include children who were not assessed for behavioral inhibition and so focus exclusively on social reticence.
      Because not all children with social reticence develop anxiety, it is important to study associations with other risk factors. Neural mechanisms underlying threat processing may affect associations between early-childhood behavior and anxiety (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ,
      • Hardee J.E.
      • Benson B.E.
      • Bar-Haim Y.
      • Mogg K.
      • Bradley B.P.
      • Chen G.
      • et al.
      Patterns of neural connectivity during an attention bias task moderate associations between early childhood temperament and internalizing symptoms in young adulthood.
      ). Moreover, threat bias, the tendency to overly attend to threatening stimuli, relates to anxiety in children and adults (
      • Bar-Haim Y.
      • Lamy D.
      • Pergamin L.
      • Bakermans-Kranenburg M.J.
      • van IJzendoorn M.H.
      Threat-related attentional bias in anxious and nonanxious individuals: A meta-analytic study.
      ,
      • Abend R.
      • de Voogd L.
      • Salemink E.
      • Wiers R.W.
      • Pérez-Edgar K.
      • Fitzgerald A.
      • et al.
      Association between attention bias to threat and anxiety symptoms in children and adolescents.
      ). Previous studies on the neural mechanisms underlying threat bias suggest that amygdala–prefrontal cortex (PFC) connectivity is stable over a 9-week period and that these measures may be more reliable than behavioral bias indices (
      • White L.K.
      • Britton J.C.
      • Sequeira S.
      • Ronkin E.G.
      • Chen G.
      • Bar-Haim Y.
      • et al.
      Behavioral and neural stability of attention bias to threat in healthy adolescents.
      ). While some studies focus on brain activity (
      • Auday E.S.
      • Taber-Thomas B.C.
      • Pérez-Edgar K.E.
      Neural correlates of attention bias to masked facial threat cues: Examining children at-risk for social anxiety disorder.
      ), most focus on amygdala connectivity during the dot-probe task. For example, children with anxiety disorders showed increased positive amygdala-insula connectivity while maintaining attention on the location of threats, whereas healthy control subjects showed increased positive amygdala-insula connectivity while shifting attention away from the location of threats (
      • White L.K.
      • Sequeira S.
      • Britton J.C.
      • Brotman M.A.
      • Gold A.L.
      • Berman E.
      • et al.
      Complementary features of attention bias modification therapy and cognitive-behavioral therapy in pediatric anxiety disorders [published correction appears in Am J Psychiatry 2018;175:83].
      ). In adults, a history of childhood behavioral inhibition was related to increased negative amygdala–dorsolateral PFC (dlPFC) connectivity while both maintaining and shifting attention from threat-related stimuli (
      • Hardee J.E.
      • Benson B.E.
      • Bar-Haim Y.
      • Mogg K.
      • Bradley B.P.
      • Chen G.
      • et al.
      Patterns of neural connectivity during an attention bias task moderate associations between early childhood temperament and internalizing symptoms in young adulthood.
      ). Finally, in socially anxious adults, avoidant orientation and slow disengagement during the dot-probe task were related to amygdala–superior temporal sulcus connectivity (
      • Evans T.C.
      • Bar-Haim Y.
      • Fox N.A.
      • Pine D.S.
      • Britton J.C.
      Neural mechanisms underlying heterogeneous expression of threat-related attention in social anxiety.
      ). These studies all focused on amygdala connectivity at one time point, whereas Abend et al. (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ) included two time points. They showed that behavioral inhibition was related to distinct neurodevelopmental pathways leading to pediatric anxiety symptoms: children with higher levels of behavioral inhibition showed a negative association between amygdala-dlPFC connectivity when maintaining attention to threat and anxiety symptoms with age, whereas children with lower levels of behavioral inhibition showed a positive association with age. The authors concluded that children with higher levels of behavioral inhibition might have an early-emerging deficiency in the capacity to regulate attention capture by threats (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ). In that study, threat processing and anxiety symptoms were measured at 10 and 13 years of age. Many anxiety symptoms typically develop in early adolescence (
      • Kessler R.C.
      • Berglund P.
      • Demler O.
      • Jin R.
      • Merikangas K.R.
      • Walters E.E.
      Lifetime prevalence and age-of-onset distributions of DSM-IV disorders in the national comorbidity survey replication [published correction appears in Arch Gen Psychiatry 2005;62:768].
      ), and amygdala connectivity with other brain regions continues to develop during adolescence (
      • Gee D.G.
      • Humphreys K.L.
      • Flannery J.
      • Goff B.
      • Telzer E.H.
      • Shapiro M.
      • et al.
      A developmental shift from positive to negative connectivity in human amygdala-prefrontal circuitry.
      ,
      • Casey B.J.
      • Glatt C.E.
      • Lee F.S.
      Treating the developing versus developed brain: Translating preclinical mouse and human studies.
      ). Hence, the current study extends these prior findings to late adolescence.
      Here, we used a sample that was previous published (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ) and added a late-adolescent (age 16) data point, a time when anxiety often manifests. We focused on social reticence rather than behavioral inhibition, which allowed us to increase our sample size by adding a second, unselected sample to the sample used in the prior work. The goal of this longitudinal study is to characterize the developmental trajectories of neural correlates underlying threat processing and relate these correlates to the association between early-childhood social reticence and adolescent anxiety. Social reticence was assessed repeatedly between 2 and 7 years of age; threat processing and anxiety symptoms were assessed at 10, 13, and 16 years of age. Based on Abend et al. (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ) in an overlapping sample, we hypothesized distinct neurodevelopmental pathways to pediatric anxiety based on social reticence: children with high social reticence were expected to show a negative association between amygdala-dlPFC connectivity when maintaining attention to threat and anxiety symptoms with age, whereas children with low social reticence were expected to show a positive association with age.

      Methods and Materials

      Participants

      Children were recruited from a longitudinal study on behavioral inhibition and early-childhood reticence. The current study includes data from 2 samples: the first sample (n = 291) was selected at 4 months of age based on reactivity to novelty (
      • Hane A.A.
      • Fox N.A.
      • Henderson H.A.
      • Marshall P.J.
      Behavioral reactivity and approach-withdrawal bias in infancy.
      ); the second sample (n = 384) was recruited from the community at 2 years of age, unselected on the basis of reactivity or other traits (
      • Jarcho J.M.
      • Davis M.M.
      • Shechner T.
      • Degnan K.A.
      • Henderson H.A.
      • Stoddard J.
      • et al.
      Early-childhood social reticence predicts brain function in preadolescent youths during distinct forms of peer evaluation.
      ,
      • Michalska K.J.
      • Feldman J.S.
      • Ivie E.J.
      • Shechner T.
      • Sequeira S.
      • Averbeck B.
      • et al.
      Early-childhood social reticence predicts SCR-BOLD coupling during fear extinction recall in preadolescent youth.
      ). Both samples were followed until 16 years of age. Data from sample 1 have previously been reported, testing similar research questions in relation to behavioral inhibition but with fewer data points (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ). The data from these samples were combined, and sample was included as a covariate in all analyses. Study procedures were approved by the National Institute of Mental Health and University of Maryland-College Park Institutional Review Boards. Informed consent and assent were obtained from parents and children, respectively.
      Threat processing was assessed in the magnetic resonance imaging (MRI) scanner at 10, 13, and 16 years of age. Functional MRI (fMRI) data were provided by 183 children for at least one time point. Ten children who provided data at age 10 were excluded from analyses (1 for subthreshold accuracy [below 70%] on the task, 1 for aborting the task, 3 for excessive head motion during the anatomical scan, 2 for excessive head motion during the functional scan [see fMRI individual-level analysis: activation], 2 for missing social reticence scores, and 1 for missing anxiety symptom scores) (see Figures S1 and S2 for an overview per sample). Fifteen children who provided data at age 13 were excluded (2 for subthreshold accuracy, 1 for aborting the task, 1 for technical issues during data collection, and 11 for missing social reticence scores). Sixteen children who provided data at age 16 were excluded (6 for subthreshold accuracy, 2 for excessive head motion during anatomical scan, 1 for aborting the task, 1 for psychotropic medication use, and 6 for missing social reticence scores). The final sample consisted of data provided by 144 participants (Table 1 and Figure S3 for a correlation matrix of the study variables). Nineteen children provided data at all three time points
      Children with fMRI data at all three time points did not differ in social reticence and anxiety symptoms from children with fMRI data at 1 or 2 time points, Fs < 3.23, ps > .08.
      and 48 children provided data at two time points. Linear mixed-effects modeling accounted for missing data in this longitudinal design, resulting in more reliable effect estimates than complete-case analysis (
      • Donders A.R.T.
      • van der Heijden G.J.
      • Stijnen T.
      • Moons K.G.
      Review: A gentle introduction to imputation of missing values.
      ,
      • Chen G.
      • Saad Z.S.
      • Britton J.C.
      • Pine D.S.
      • Cox R.W.
      Linear mixed-effects modeling approach to FMRI group analysis.
      ,
      • Matta T.H.
      • Flournoy J.C.
      • Byrne M.L.
      Making an unknown unknown a known unknown: Missing data in longitudinal neuroimaging studies.
      ). Four children were Asian (2.78%), 24 were Black or African American (16.67%), 5 were Hispanic or Latino (3.47%), 16 were multiracial (11.11%), and 83 were White (57.64%). The race or ethnicity of 12 participants was unknown (8.33%). fMRI data from 61 children from sample 1 at age 10 and from 64 children from sample 1 at age 13 have been previously published, testing similar research questions on behavioral inhibition (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ).
      Table 1Overview of the Demographics, Clinical Indices, and Dot-Probe Performance in the Current Sample
      MeasuresAge 10Age 13Age 16
      Demographics
      Total, n (n female)71 (43 F)85 (46 F)74 (36 F)
       Sample 159 (35 F)65 (38 F)43 (24 F)
       Sample 212 (8 F)20 (8 F)31 (12 F)
      Age, Years, Mean (SD) [Range]10.62 (0.46) [9.50 to 11.67]13.25 (0.76) [12.23 to 15.49]16.27 (0.68) [12.92 to 17.69]
      IQ
      IQ and social reticence were measured at one time point. The means seem to change over time, because of the different number of children with usable functional magnetic resonance imaging data at each time point.
      ,
      IQ scores were missing for 6 participants.
      , Mean (SD) [Range]
      115.17 (11.97) [88 to 138]116.01 (13.33) [82 to 138]115.58 (13.99) [77 to 141]
      Clinical Indices
      Current Anxiety Disorder, n10103
      Social Reticence
      IQ and social reticence were measured at one time point. The means seem to change over time, because of the different number of children with usable functional magnetic resonance imaging data at each time point.
      , Mean (SD) [Range]
      −0.02 (0.45) [−0.57 to 1.77]0.03 (0.53) [−0.62 to 2.21]−0.04 (0.47) [−0.83 to 2.21]
      Anxiety Symptoms, Mean (SD) [Range]17.35 (9.16) [2 to 42]10.23 (7.17) [0 to 32]9.15 (7.72) [0 to 29.5]
      Social Anxiety, Mean (SD) [Range]3.95 (2.77) [0 to 12]3.94 (2.86) [0 to 12]3.32 (2.87) [0 to 12]
      Dot-Probe Performance
      Accuracy, %, Mean (SD) [Range]89.61% (7.53%) [69.20% to 99.00%]91.86% (5.78%) [73.80% to 100%]92.98% (4.40%) [80.40% to 99.40%]
      Threat Bias, ms, Mean (SD) [Range]2.43 (31.99) [−73.42 to 79.25]6.62 (33.04) [−55.89 to 106.44]7.61 (21.47) [−48.85 to 96.01]
      Happy Bias, ms, Mean (SD) [Range]7.06 (30.48) [−54.42 to 79.48]5.60 (30.60) [−68.91 to 89.80]0.78 (23.29) [−77.73 to 56.93]
      ABV Threat, Mean (SD) [Range]0.05 (0.02) [0.01 to 0.11]0.05 (0.02) [0.01 to 0.14]0.06 (0.02) [0.03 to 0.14]
      ABV Happy, Mean (SD) [Range]0.05 (0.02) [0.02 to 0.11]0.05 (0.02) [0.02 to 0.10]0.06 (0.02) [0.03 to 0.15]
      ABV, attention bias variability; F, female.
      a IQ and social reticence were measured at one time point. The means seem to change over time, because of the different number of children with usable functional magnetic resonance imaging data at each time point.
      b IQ scores were missing for 6 participants.

      Social Reticence

      Social reticence was observed by independent raters during free play, cleanup, and social problem-solving interactions with unfamiliar age- and sex-matched peers at 2, 3, 4, 5, and 7 years of age [see Supplemental Methods and Degnan et al. (
      • Degnan K.A.
      • Almas A.N.
      • Henderson H.A.
      • Hane A.A.
      • Walker O.L.
      • Fox N.A.
      Longitudinal trajectories of social reticence with unfamiliar peers across early childhood.
      ,
      • Degnan K.A.
      • Hane A.A.
      • Henderson H.A.
      • Moas O.L.
      • Reeb-Sutherland B.C.
      • Fox N.A.
      Longitudinal stability of temperamental exuberance and social-emotional outcomes in early childhood.
      )]. Children from sample 1 were randomly paired with children from sample 2 for these interactions.
      Dyadic nesting was not taken into account in calculating the composite score because children were paired with different children at all time points, which would average out the effects of partner.
      All behavior scores were standardized within time point and averaged together to create a social reticence composite (Table S1A, B for the correlations among time points
      It should be noted that these correlations are modest and that some of the correlations are lower in the current sample than in the full sample.
      ). Social reticence was correlated with behavioral inhibition, r = 0.32, p < .001.

      Anxiety Symptoms

      Anxiety symptoms were assessed within 6 weeks of each scan using the Screen for Child Anxiety Related Emotional Disorders (SCARED), a reliable 41-item child- and parent-report measure of anxiety symptomology (Figure S4) (
      • Birmaher B.
      • Khetarpal S.
      • Brent D.
      • Cully M.
      • Balach L.
      • Kaufman J.
      • Neer S.M.
      The screen for child anxiety related emotional disorders (SCARED): Scale construction and psychometric characteristics.
      ,
      • Muris P.
      • Dreessen L.
      • Bögels S.
      • Weckx M.
      • van Melick M.
      A questionnaire for screening a broad range of DSM-defined anxiety disorder symptoms in clinically referred children and adolescents.
      ). Internal consistency of the SCARED is good for both self-report and parent-report (Cronbach’s α > 0.74) (
      • Birmaher B.
      • Khetarpal S.
      • Brent D.
      • Cully M.
      • Balach L.
      • Kaufman J.
      • Neer S.M.
      The screen for child anxiety related emotional disorders (SCARED): Scale construction and psychometric characteristics.
      ). Participants missing more than 20% of questions would be excluded (n = 0); missing data were replaced with the mean of the available items. A total score of 25 or higher has been suggested to indicate clinically significant anxiety (
      • Rappaport B.I.
      • Pagliaccio D.
      • Pine D.S.
      • Klein D.N.
      • Jarcho J.M.
      Discriminant validity, diagnostic utility, and parent-child agreement on the screen for child anxiety related emotional disorders (SCARED) in treatment- and non-treatment-seeking youth.
      ). Total scores from each parent and youth pair (r = 0.37, p < .001) were averaged into a mean anxiety score per participant for each time point (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ) to mitigate informant differences (
      • Behrens B.
      • Swetlitz C.
      • Pine D.S.
      • Pagliaccio D.
      The Screen for Child Anxiety Related Emotional Disorders (SCARED): Informant discrepancy, measurement invariance, and test-retest reliability.
      ). Anxiety scores were mean centered at each time point.

      Dot-Probe Task

      Threat-related attention was assessed using the dot-probe task (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ,
      • Bar-Haim Y.
      • Lamy D.
      • Pergamin L.
      • Bakermans-Kranenburg M.J.
      • van IJzendoorn M.H.
      Threat-related attentional bias in anxious and nonanxious individuals: A meta-analytic study.
      ,
      • White L.K.
      • Sequeira S.
      • Britton J.C.
      • Brotman M.A.
      • Gold A.L.
      • Berman E.
      • et al.
      Complementary features of attention bias modification therapy and cognitive-behavioral therapy in pediatric anxiety disorders [published correction appears in Am J Psychiatry 2018;175:83].
      ,
      • Macleod C.
      • Mathews A.
      • Tata P.
      Attentional bias in emotional disorders.
      ). Participants completed the same version of the fMRI dot-probe task at each time point (Figure 1; Supplemental Methods). This task included 5 trial types: angry-congruent (48 trials), angry-incongruent (48 trials), happy-congruent (48 trials), happy-incongruent (48 trials), and neutral-neutral (96 trials). Trials occurred across 4 blocks of 4 minutes and 15 seconds; neural and behavioral responses to each trial type were recorded. The task was programmed and administered using E-Prime (Psychology Software Tools) and included faces from the validated NimStim set (
      • Tottenham N.
      • Tanaka J.W.
      • Leon A.C.
      • McCarry T.
      • Nurse M.
      • Hare T.A.
      • et al.
      The NimStim set of facial expressions: Judgments from untrained research participants.
      ). In line with prior work, trials were included in analyses only if reaction time (RT) was between 150 and 2000 ms and <2.5 standard deviations from the child’s mean RT and if the probe type was correctly identified. Children with mean accuracy below 70% at 1 time point were excluded from that time point (n = 9) (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ,
      • White L.K.
      • Sequeira S.
      • Britton J.C.
      • Brotman M.A.
      • Gold A.L.
      • Berman E.
      • et al.
      Complementary features of attention bias modification therapy and cognitive-behavioral therapy in pediatric anxiety disorders [published correction appears in Am J Psychiatry 2018;175:83].
      ).
      Figure thumbnail gr1
      Figure 1Overview of the dot-probe task. Reproduced with permission from Harrewijn et al. (
      • Harrewijn A.
      • Abend R.
      • Linke J.
      • Brotman M.A.
      • Fox N.A.
      • Leibenluft E.
      • et al.
      Combining fMRI during resting state and an attention bias task in children.
      ).
      Behavioral analyses included two indices of attentional processes: attention bias scores and attention bias variability (ABV) (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ,
      • Bar-Haim Y.
      • Lamy D.
      • Pergamin L.
      • Bakermans-Kranenburg M.J.
      • van IJzendoorn M.H.
      Threat-related attentional bias in anxious and nonanxious individuals: A meta-analytic study.
      ,
      • Shechner T.
      • Bar-Haim Y.
      Threat monitoring and attention-bias modification in anxiety and stress-related disorders.
      ). Threat bias scores were computed by subtracting the mean RT in threat-congruent trials from the mean RT in threat-incongruent trials, with positive scores (i.e., faster responses to threat-congruent trials) reflecting attention bias toward threat (
      • White L.K.
      • Sequeira S.
      • Britton J.C.
      • Brotman M.A.
      • Gold A.L.
      • Berman E.
      • et al.
      Complementary features of attention bias modification therapy and cognitive-behavioral therapy in pediatric anxiety disorders [published correction appears in Am J Psychiatry 2018;175:83].
      ,
      • Badura-Brack A.S.
      • Naim R.
      • Ryan T.J.
      • Levy O.
      • Abend R.
      • Khanna M.M.
      • et al.
      Effect of attention training on attention bias variability and PTSD symptoms: Randomized controlled trials in Israeli and U.S. combat veterans.
      ). Happy bias scores were computed with the same procedure for happy trials. Both angry and happy stimuli were considered to test for bias specificity to threat.
      ABV measures fluctuations in attention allocation and has been suggested to reflect impairments in attentional control with relation to anxiety (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ,
      • Iacoviello B.M.
      • Wu G.
      • Abend R.
      • Murrough J.W.
      • Feder A.
      • Fruchter E.
      • et al.
      Attention bias variability and symptoms of posttraumatic stress disorder.
      ). Using a moving-window algorithm, threat bias scores were calculated (as described above) for 10 successive angry-neutral trials. To control for associations between mean and variance, the standard deviation of these scores was calculated and then divided by the participant’s overall mean RT (
      • Badura-Brack A.S.
      • Naim R.
      • Ryan T.J.
      • Levy O.
      • Abend R.
      • Khanna M.M.
      • et al.
      Effect of attention training on attention bias variability and PTSD symptoms: Randomized controlled trials in Israeli and U.S. combat veterans.
      ,
      • Naim R.
      • Abend R.
      • Wald I.
      • Eldar S.
      • Levi O.
      • Fruchter E.
      • et al.
      Threat-related attention bias variability and posttraumatic stress.
      ). The same procedure was applied to happy-neutral trials to calculate happy ABV scores. Table S2 shows the intraclass correlations across time points for the behavioral measures.

      Behavioral Data Analysis

      To test associations among social reticence, anxiety, and attention bias at all time points, we ran two linear mixed-effects models in R (nlme package) (
      • Pinheiro J.
      • Bates D.
      • DebRoy S.
      • Sarkar D.
      • Heisterkamp S.
      • et al.
      EISPACK Authors
      Package ‘nlme’: Linear and nonlinear mixed effects models.
      ) with attention bias or ABV scores as the dependent variable. Anxiety symptoms (SCARED scores), social reticence, task condition (threat and happy), and time (age 10, 13, and 16 years) were included as independent variables; sample (1 or 2) and sex (male or female) were included as covariates; subject was modeled as random effect. All statistical tests were two sided; significance threshold was set to α < 0.05. Follow-up paired t tests were conducted for significant effects using the phia (
      • De Rosario-Marinez H.
      • Fox J.
      • CoreTeam
      Package ‘phia’: Post-hoc interaction analysis..
      ) and emmeans (
      • Lenth R.V.
      • Buerkner P.
      • Giné-Vázquez I.
      • Herve M.
      • Jung M.
      • Love J.
      • et al.
      emmeans: Estimated marginal means, aka least-squares means..
      ) packages.

      fMRI Individual-Level Analysis

      fMRI data were collected during the dot-probe task on two identical 3T MR750 General Electric scanners (Waukesha) with identical 32-channel head coils and scan parameters. fMRI data were analyzed with AFNI version 20.2 (
      • Cox R.W.
      AFNI: Software for analysis and visualization of functional magnetic resonance neuroimages.
      ,
      • Cox R.W.
      • Hyde J.S.
      Software tools for analysis and visualization of fMRI data.
      ). Data preprocessing included the following steps: removing the first 4 pre–steady-state volumes, despiking, slice timing correction, aligning functional scans to structural scans, nonlinear registration to standard Talairach space (TT_N27 template), volume registration, spatial smoothing with a 6 mm (full width at half maximum) Gaussian kernel, and scaling (allowing interpretation of effect estimates values as percent signal change relative to the mean). We used generalized psychophysiological interaction (gPPI) (
      • McLaren D.G.
      • Ries M.L.
      • Xu G.
      • Johnson S.C.
      A generalized form of context-dependent psychophysiological interactions (gPPI): A comparison to standard approaches.
      ) analysis to assess task-specific differences in amygdala functional connectivity. We used FreeSurfer (
      • Fischl B.
      • Salat D.H.
      • Busa E.
      • Albert M.
      • Dieterich M.
      • Haselgrove C.
      • et al.
      Whole brain segmentation: Automated labeling of neuroanatomical structures in the human brain.
      ) to extract subject-specific amygdala seeds. Then, we ran a whole-brain individual-level general linear model including correct responses in the 5 conditions, incorrect responses, 6 motion regressors (displacement in x, y, and z axes, rotational movement of roll, pitch, and yaw), time series in the amygdala seed, and the PPI terms (products of the detrended and demeaned seed and the 5 task condition regressors). gPPI analysis was run for left and right amygdala separately.

      fMRI Group-Level Analysis (gPPI)

      To test associations between social reticence, anxiety, and task-specific functional connectivity at all time points, we ran a linear mixed-effects model using AFNI’s 3dLMEr (
      • Chen G.
      • Saad Z.S.
      • Britton J.C.
      • Pine D.S.
      • Cox R.W.
      Linear mixed-effects modeling approach to FMRI group analysis.
      ) with the output from the individual-level gPPI as the dependent variable. Anxiety symptoms (total SCARED scores), social reticence, task condition (angry-congruent, angry-incongruent, happy-congruent, happy-incongruent, and neutral-neutral), and time (age 10, 13, and 16 years)
      Time point (instead of actual age) was included because we do not assume similar age effects across the whole range (e.g., between 10 and 13 years and 13 and 16 years).
      were included as independent variables; recruitment source (sample 1 or 2) and sex (female or male) were included as covariates; subject was modeled as random effect. AFNI’s 3dClustSim (
      • Cox R.W.
      • Chen G.
      • Glen D.R.
      • Reynolds R.C.
      • Taylor P.A.
      FMRI clustering in AFNI: False-positive rates redux.
      ), which assumes a non-Gaussian autocorrelation smoothing function (
      • Cox R.W.
      • Chen G.
      • Glen D.R.
      • Reynolds R.C.
      • Taylor P.A.
      FMRI clustering in AFNI: False-positive rates redux.
      ) in light of Eklund et al. (
      • Eklund A.
      • Nichols T.E.
      • Knutsson H.
      Cluster failure: Why fMRI inferences for spatial extent have inflated false-positive rates [published correction appears in Proc Natl Acad Sci U S A 2016;113:E4929].
      ), was run on all data and showed a minimum size of 45 contiguous voxels (Nearest Neighbor = 2) for a voxelwise threshold of p < .001 (two sided) and a clusterwise α < 0.05, based on 10,000 Monte-Carlo simulations. For significant effects, we extracted individual-level betas for post hoc analysis using linear mixed-effects models in R (nlme package) (
      • Pinheiro J.
      • Bates D.
      • DebRoy S.
      • Sarkar D.
      • Heisterkamp S.
      • et al.
      EISPACK Authors
      Package ‘nlme’: Linear and nonlinear mixed effects models.
      ).

      Results

      Behavioral Analyses

      No significant main or interaction effects manifested for social reticence, anxiety, task condition, or time on attention bias scores, p > .10 (Table S3). There was a significant main effect of time on ABV scores, with children showing greater ABV at age 16 versus at age 10, b = 0.01, p < .001 and at age 16 versus at age 13, b = 0.01, p < .001 (Figure 2). There were no other significant main or interaction effects, ps > .19 (Table S4).
      Figure thumbnail gr2
      Figure 2Effect of time on attention bias variability (ABV), averaged across angry-neutral and happy-neutral trials. ∗∗∗p < .001.

      Task-Specific Functional Connectivity

      Left Amygdala

      No brain regions showed a significant 4-way interaction among social reticence, time point, anxiety symptoms, and condition in left amygdala functional connectivity. However, there was a significant 3-way interaction among social reticence, anxiety symptoms, and time with the left dlPFC (coordinates: 46.2, −11.2, 38.8; k = 62) and the right occipital cortex (coordinates: −43.8, 58.8, 31.2; k = 53) (voxelwise threshold of p < .001 and a clusterwise familywise error p < .05). Results for the lower-order interactions and main effects are provided in Table S5.
      To decompose these interaction effects, follow-up analyses were conducted separately for children with high versus low levels of social reticence (median split) (Figure 3A; see Figure S5 for Johnson-Neyman plots). For children with high social reticence, simple slope analysis indicated that anxiety and amygdala-dlPFC connectivity were positively associated at 10 years of age (b = 0.02; 95% CI, 0.0004–0.03), not associated at 13 years of age (b = −0.01; 95% CI, −0.03 to 0.002), and negatively associated at 16 years of age (b = −0.2; 95% CI, −0.03 to −0.002). These associations were different between time points 10 and 13 and 10 and 16, b10–13yrs = 0.03, b10–16yrs = 0.03, χ21s > 11.63, ps < .001, but not between time points 13 and 16, b13–16yrs = 0.004, χ21 = 0.18, p = .67. A different trend emerged for children with low social reticence, whereby anxiety and amygdala-dlPFC connectivity were not associated at 10 years of age (b = −0.002; 95% CI, −0.02 to 0.01), 13 years of age (b = −0.007; 95% CI, −0.02 to 0.01), and 16 years of age (b = 0.01; 95% CI, −0.003 to 0.03). These associations were different between time points 13 and 16, b13–16yrs = −0.02, χ21 = 4.83, p = .03, but not between time points 10 and 13 or 10 and 16, b10–13yrs = 0.004, b10–16yrs = −0.01, χ21s < 3.04, ps > .08.
      Figure thumbnail gr3
      Figure 3Three-way interaction among social reticence (SR), anxiety symptoms, and time point on functional connectivity of the left amygdala with left dorsolateral prefrontal cortex (dlPFC) (A) and right occipital cortex (B). SR was dichotomized (by median split) for follow-up analyses.
      Follow-up analyses for amygdala–occipital cortex connectivity (Figure 3B; Figure S6 for Johnson-Neyman plots) indicated that for children with high social reticence, the relationship between anxiety and functional connectivity was positive at 10 years of age (b = 0.01; 95% CI, 0.001–0.02), negative at 13 years of age (b = −0.01; 95% CI, −0.02 to −0.0003) and not associated at 16 years of age (b = −0.01; 95% CI, −0.02 to 0.01). These associations were different between time points 10 and 13 and 10 and 16, b10–13yrs = 0.03, b10–16yrs = 0.02, χ21s > 6.75, ps < .01, but not between time points 13 and 16, b13–16yrs = −0.01, χ21 = 0.57, p = .45. In contrast, for children with low social reticence, the relationship between anxiety and functional connectivity was not significant at 10 years of age (b = −0.01; 95% CI, −0.02 to 0.002), 13 years of age (b = −0.01; 95% CI, −0.02 to 0.004), and 16 years of age (b = 0.01; 95% CI, −0.002 to 0.02). These associations were different between time points 10 and 16 and 13 and 16, b10–16yrs = −0.02, b13–16yrs = −0.02, χ21s > 6.36, ps < .05, but not between time points 10 and 13, b10–13yrs = −0.002, χ21 = 0.10, p = .75.

      Right Amygdala

      No brain regions showed significant 4-way interactions among social reticence, time point, anxiety symptoms, and condition in functional connectivity with the right amygdala. Unlike for the left amygdala, no 3-way interaction occurred among social reticence, anxiety symptoms, and time point for the dlPFC. There was a 3-way interaction in the cerebellum, middle occipital cortex, superior temporal gyrus, and thalamus, but these will not be discussed here because we hypothesized findings in the dlPFC. Results for the lower-order interactions and main effects are shown in Table S6.

      Discussion

      The goal of this study was to characterize the developmental trajectories of neural mechanisms underlying threat processing and relate these trajectories to associations between early-childhood social reticence and adolescent anxiety. We continued the study by Abend et al. (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ) by including an adolescent time point and a second unselected sample and examining social reticence rather than behavioral inhibition. Attention bias and ABV were not associated with social reticence and anxiety symptoms. Children with high social reticence showed a negative association between amygdala-dlPFC connectivity and anxiety symptoms with age, compared to children with low social reticence. However, this amygdala-dlPFC connectivity was in response to potential threat processing in general and not specific to threat bias.
      Differential amygdala-PFC connectivity has been found in studies examining threat processing and anxiety (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ,
      • Hardee J.E.
      • Benson B.E.
      • Bar-Haim Y.
      • Mogg K.
      • Bradley B.P.
      • Chen G.
      • et al.
      Patterns of neural connectivity during an attention bias task moderate associations between early childhood temperament and internalizing symptoms in young adulthood.
      ,
      • White L.K.
      • Sequeira S.
      • Britton J.C.
      • Brotman M.A.
      • Gold A.L.
      • Berman E.
      • et al.
      Complementary features of attention bias modification therapy and cognitive-behavioral therapy in pediatric anxiety disorders [published correction appears in Am J Psychiatry 2018;175:83].
      ). Amygdala-PFC connectivity supports emotion regulation and develops in adolescence (
      • Gee D.G.
      • Humphreys K.L.
      • Flannery J.
      • Goff B.
      • Telzer E.H.
      • Shapiro M.
      • et al.
      A developmental shift from positive to negative connectivity in human amygdala-prefrontal circuitry.
      ). The negative association between the amygdala and the dlPFC might suggest less efficient emotion regulation during threat processing. Here, we show that the combination of high social reticence in early childhood and strong negative amygdala-PFC connectivity at 13 and 16 years of age increases risk for adolescent anxiety symptoms. Understanding relationships between socially reticent behavior and brain function could inform targeted prevention. This study also found that overall response to potential threats, rather than responses to threats in association with attention demands, relates to anxiety symptoms. This could suggest that future studies attempting to identify risk might focus on overall levels of potential threat responsivity.
      Similar neurodevelopmental pathways were reported by Abend et al. (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ) in an overlapping sample. Extending these past findings into late adolescence builds on findings showing risk for anxiety disorders to change past 13 years of age (
      • Kessler R.C.
      • Berglund P.
      • Demler O.
      • Jin R.
      • Merikangas K.R.
      • Walters E.E.
      Lifetime prevalence and age-of-onset distributions of DSM-IV disorders in the national comorbidity survey replication [published correction appears in Arch Gen Psychiatry 2005;62:768].
      ), a time when amygdala connectivity also changes (
      • Gee D.G.
      • Humphreys K.L.
      • Flannery J.
      • Goff B.
      • Telzer E.H.
      • Shapiro M.
      • et al.
      A developmental shift from positive to negative connectivity in human amygdala-prefrontal circuitry.
      ,
      • Casey B.J.
      • Glatt C.E.
      • Lee F.S.
      Treating the developing versus developed brain: Translating preclinical mouse and human studies.
      ). We showed that the timing of neurodevelopmental pathways might differ for children with high versus low early social reticence. Children with high social reticence showed an early change (from 10 to 13 years of age) from a negative to a positive association between amygdala-dlPFC connectivity and anxiety; in contrast, children with low social reticence showed a later change (from 13 to 16 years of age) in the association between amygdala-dlPFC connectivity and anxiety. Because high social reticence places children at risk for anxiety, the current findings could inform timing of preventive interventions in these children. More specifically, preventive interventions might work better when started before this early change at 13 years of age.
      Beyond amygdala-dlPFC connectivity, other findings manifested for amygdala–occipital cortex connectivity. This has not been reported previously because previous studies have focused exclusively on connectivity between the amygdala and the PFC and/or the insula (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ,
      • Hardee J.E.
      • Benson B.E.
      • Bar-Haim Y.
      • Mogg K.
      • Bradley B.P.
      • Chen G.
      • et al.
      Patterns of neural connectivity during an attention bias task moderate associations between early childhood temperament and internalizing symptoms in young adulthood.
      ,
      • White L.K.
      • Sequeira S.
      • Britton J.C.
      • Brotman M.A.
      • Gold A.L.
      • Berman E.
      • et al.
      Complementary features of attention bias modification therapy and cognitive-behavioral therapy in pediatric anxiety disorders [published correction appears in Am J Psychiatry 2018;175:83].
      ,
      • Evans T.C.
      • Bar-Haim Y.
      • Fox N.A.
      • Pine D.S.
      • Britton J.C.
      Neural mechanisms underlying heterogeneous expression of threat-related attention in social anxiety.
      ). Future studies reporting whole-brain findings are necessary to reveal whether these effects are chance findings or more stable patterns in the dot-probe task.
      The direction of this connectivity (positive vs. negative correlation) differs between studies. Moreover, other studies [including Abend et al. (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ) in an overlapping sample] have found connectivity with the right, instead of the left, amygdala. However, it should be noted that the literature on laterality of the amygdala in emotion processing is unclear (
      • Ocklenburg S.
      • Peterburs J.
      • Mundorf A.
      Hemispheric asymmetries in the amygdala: A comparative primer.
      ,
      • McMenamin B.W.
      • Marsolek C.J.
      Can theories of visual representation help to explain asymmetries in amygdala function?.
      ). Finally, the associations between amygdala-PFC connectivity, behavior, and anxiety are complex. Some have found significant amygdala-PFC connectivity only when maintaining attention to threat [in an overlapping sample (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      )] (
      • White L.K.
      • Sequeira S.
      • Britton J.C.
      • Brotman M.A.
      • Gold A.L.
      • Berman E.
      • et al.
      Complementary features of attention bias modification therapy and cognitive-behavioral therapy in pediatric anxiety disorders [published correction appears in Am J Psychiatry 2018;175:83].
      ), whereas others have found this across all task conditions (S.P. Haller et al., unpublished data, May 2022). More broadly, reliability in fMRI paradigms is lower for contrasts involving conditions with similar demands (e.g., angry vs. happy faces) than for contrasts involving distinct conditions (e.g., faces vs. objects) (
      • Haller S.P.
      • Chen G.
      • Kitt E.R.
      • Smith A.R.
      • Stoddard J.
      • Abend R.
      • et al.
      Reliability of task-evoked neural activation during face-emotion paradigms: Effects of scanner and psychological processes.
      ). We found amygdala-dlPFC connectivity across task conditions (i.e., no 4-way interaction), possibly reflecting individual differences in potential threat processing generally, rather than specific aspects of threat processing. Such a task-general effect is in line with the absence of behavioral findings in our study and with findings on differential amygdala-PFC connectivity during resting state in children with anxiety disorders (
      • Birn R.M.
      • Shackman A.J.
      • Oler J.A.
      • Williams L.E.
      • McFarlin D.R.
      • Rogers G.M.
      • et al.
      Evolutionarily conserved prefrontal-amygdalar dysfunction in early-life anxiety.
      ). These task-general effects can be studied in two different ways in task-based fMRI. First, functional connectivity can be computed independently of task events. This task-general connectivity is more reflective of individual differences in behavior than functional connectivity during resting state (
      • Finn E.S.
      • Shen X.
      • Scheinost D.
      • Rosenberg M.D.
      • Huang J.
      • Chun M.M.
      • et al.
      Functional connectome fingerprinting: Identifying individuals using patterns of brain connectivity.
      ,
      • Greene A.S.
      • Gao S.
      • Scheinost D.
      • Constable R.T.
      Task-induced brain state manipulation improves prediction of individual traits.
      ) but has not yet been applied to the dot-probe task. Second, the similarity between task-general connectivity and resting-state connectivity is also a measure of task-general effects and is related to task performance (
      • Schultz D.H.
      • Cole M.W.
      Higher intelligence is associated with less task-related brain network reconfiguration.
      ). In the dot-probe task, task performance is related to the similarity between task-general and resting-state functional connectivity (
      • Harrewijn A.
      • Abend R.
      • Linke J.
      • Brotman M.A.
      • Fox N.A.
      • Leibenluft E.
      • et al.
      Combining fMRI during resting state and an attention bias task in children.
      ). These studies suggest that measures other than task-specific connectivity may be more predictive of individual differences in the dot-probe task.
      The dot-probe task rightfully has been criticized due to low reliability for the behaviors it engages (
      • White L.K.
      • Britton J.C.
      • Sequeira S.
      • Ronkin E.G.
      • Chen G.
      • Bar-Haim Y.
      • et al.
      Behavioral and neural stability of attention bias to threat in healthy adolescents.
      ,
      • MacLeod C.
      • Grafton B.
      • Notebaert L.
      Anxiety-linked attentional bias: Is it reliable?.
      ,
      • Rodebaugh T.L.
      • Scullin R.B.
      • Langer J.K.
      • Dixon D.J.
      • Huppert J.D.
      • Bernstein A.
      • et al.
      Unreliability as a threat to understanding psychopathology: The cautionary tale of attentional bias [published correction appears in J Abnorm Psychol 2016;125:878].
      ,
      • Price R.B.
      • Kuckertz J.M.
      • Siegle G.J.
      • Ladouceur C.D.
      • Silk J.S.
      • Ryan N.D.
      • et al.
      Empirical recommendations for improving the stability of the dot-probe task in clinical research [published correction appears in Psychol Assess 2015;27:376].
      ), consistent with lack of behavioral effects in the current study. Moreover, the findings with ABV as reported in Abend et al. (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ) in an overlapping sample were not found here. The association between behavioral inhibition, anxiety, and fluctuations in attention allocation was found at 13 and not at 10 years of age. The lack of findings with ABV in the current study could suggest that these associations are less sensitive to developmental effects beyond age 13, which would explain why there was no interaction between social reticence, anxiety, and time. Even though other measures during the dot-probe task, such as computational (
      • Evans T.C.
      • Bar-Haim Y.
      • Fox N.A.
      • Pine D.S.
      • Britton J.C.
      Neural mechanisms underlying heterogeneous expression of threat-related attention in social anxiety.
      ,
      • Price R.B.
      • Brown V.
      • Siegle G.J.
      Computational modeling applied to the dot-probe task yields improved reliability and mechanistic insights.
      ) and neuroimaging (
      • White L.K.
      • Britton J.C.
      • Sequeira S.
      • Ronkin E.G.
      • Chen G.
      • Bar-Haim Y.
      • et al.
      Behavioral and neural stability of attention bias to threat in healthy adolescents.
      ) measures, are more reliable, it is problematic to relate these other measures to unreliable behavioral measures. One possible solution is to relate neuroimaging findings to real-life measures of threat processing. For example, decreased amygdala-dlPFC connectivity during the dot-probe task was associated with using more distraction after negative events in real life (
      • Price R.B.
      • Allen K.B.
      • Silk J.S.
      • Ladouceur C.D.
      • Ryan N.D.
      • Dahl R.E.
      • et al.
      Vigilance in the laboratory predicts avoidance in the real world: A dimensional analysis of neural, behavioral, and ecological momentary data in anxious youth.
      ). Moreover, increased amygdala–anterior PFC connectivity while receiving social feedback was associated with a real-life attention bias toward a potentially critical judge (
      • Sequeira S.L.
      • Rosen D.K.
      • Silk J.S.
      • Hutchinson E.
      • Allen K.B.
      • Jones N.P.
      • et al.
      Don’t judge me!”: Links between in vivo attention bias toward a potentially critical judge and fronto-amygdala functional connectivity during rejection in adolescent girls.
      ). These studies are examples of how neuroimaging findings could be related to real-life processes to enhance ecologic validity.
      To our knowledge, this was the first study to include three waves of longitudinal neuroimaging data in a sample well characterized for social reticence, continuing previous research (
      • Abend R.
      • Swetlitz C.
      • White L.K.
      • Shechner T.
      • Bar-Haim Y.
      • Filippi C.
      • et al.
      Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
      ) by focusing on social reticence, including data from late adolescence, and adding a second unselected sample. A few limitations should be noted. First, children were recruited from two samples, one group was selected at 4 months of age based on reactivity to novelty and one unselected group was recruited from the community at 2 years of age. However, children from both samples only differed in age, and we accounted for sample in all analyses. Second, attrition from inclusion at 4 months or 2 years of age to inclusion in the MRI scans was high. Third, this was not a clinical sample, so anxiety symptoms were relatively low. Fourth, adult faces were used in the dot-probe task. It is important to study whether these effects are similar when adolescent faces are used. Fifth, puberty was assessed, but the analyses would have been underpowered if puberty were included. It should be noted that puberty could also be related to the change from negative to positive connectivity-anxiety association from 10 to 13 years of age. Sixth, it is not possible to estimate effect sizes in linear mixed-effects models with fMRI data. It should be noted that effect sizes in fMRI research in general are typically no larger than medium. This might limit the actionable conclusions drawn from this study, so these findings should be replicated in a larger, more diverse sample.
      To conclude, we found distinct neurodevelopmental pathways to pediatric anxiety based on social reticence: children with high social reticence showed a negative association between amygdala-dlPFC connectivity and anxiety symptoms with age compared to children with low social reticence. These patterns of functional connectivity were present across all task conditions, early in children with high social reticence, and later in children with low social reticence.

      Acknowledgments and Disclosures

      This research was supported by the Intramural Research Program of the National Institute of Mental Health, United States (Grant Nos. ZIA-MH-002782 and NCT00018057 [to DSP]) and by the National Institute of Mental Health, United States (Grant Nos. R37HD17899 and U01MH093349 [to NAF]). The funding sources did not have a role in the writing of the manuscript or the decision to submit it for publication.
      The authors report no biomedical financial interests or potential conflicts of interest.

      Supplementary Material

      References

        • Clauss J.A.
        • Blackford J.U.
        Behavioral inhibition and risk for developing social anxiety disorder: A meta-analytic study.
        J Am Acad Child Adolesc Psychiatry. 2012; 51: 1066-1075.e1
        • Sandstrom A.
        • Uher R.
        • Pavlova B.
        Prospective association between childhood behavioral inhibition and anxiety: A meta-analysis.
        J Abnorm Child Psychol. 2020; 48: 57-66
        • Abend R.
        • Swetlitz C.
        • White L.K.
        • Shechner T.
        • Bar-Haim Y.
        • Filippi C.
        • et al.
        Levels of early-childhood behavioral inhibition predict distinct neurodevelopmental pathways to pediatric anxiety.
        Psychol Med. 2020; 50: 96-106
        • Hardee J.E.
        • Benson B.E.
        • Bar-Haim Y.
        • Mogg K.
        • Bradley B.P.
        • Chen G.
        • et al.
        Patterns of neural connectivity during an attention bias task moderate associations between early childhood temperament and internalizing symptoms in young adulthood.
        Biol Psychiatry. 2013; 74: 273-279
        • Chronis-Tuscano A.
        • Danko C.M.
        • Rubin K.H.
        • Coplan R.J.
        • Novick D.R.
        Future directions for research on early intervention for young children at risk for social anxiety.
        J Clin Child Adolesc Psychol. 2018; 47: 655-667
        • Rubin K.H.
        • Coplan R.J.
        • Bowker J.C.
        Social withdrawal in childhood.
        Annu Rev Psychol. 2009; 60: 141-171
        • Fox N.A.
        • Henderson H.A.
        • Rubin K.H.
        • Calkins S.D.
        • Schmidt L.A.
        Continuity and discontinuity of behavioral inhibition and exuberance: Psychophysiological and behavioral influences across the first four years of life.
        Child Dev. 2001; 72: 1-21
        • Degnan K.A.
        • Almas A.N.
        • Henderson H.A.
        • Hane A.A.
        • Walker O.L.
        • Fox N.A.
        Longitudinal trajectories of social reticence with unfamiliar peers across early childhood.
        Dev Psychol. 2014; 50: 2311-2323
        • Nelson E.E.
        • Jarcho J.M.
        • Guyer A.E.
        Social re-orientation and brain development: An expanded and updated view.
        Dev Cogn Neurosci. 2016; 17: 118-127
        • Clarkson T.
        • Eaton N.R.
        • Nelson E.E.
        • Fox N.A.
        • Leibenluft E.
        • Pine D.S.
        • et al.
        Early childhood social reticence and neural response to peers in preadolescence predict social anxiety symptoms in midadolescence.
        Depress Anxiety. 2019; 36: 676-689
        • Bar-Haim Y.
        • Lamy D.
        • Pergamin L.
        • Bakermans-Kranenburg M.J.
        • van IJzendoorn M.H.
        Threat-related attentional bias in anxious and nonanxious individuals: A meta-analytic study.
        Psychol Bull. 2007; 133: 1-24
        • Abend R.
        • de Voogd L.
        • Salemink E.
        • Wiers R.W.
        • Pérez-Edgar K.
        • Fitzgerald A.
        • et al.
        Association between attention bias to threat and anxiety symptoms in children and adolescents.
        Depress Anxiety. 2018; 35: 229-238
        • White L.K.
        • Britton J.C.
        • Sequeira S.
        • Ronkin E.G.
        • Chen G.
        • Bar-Haim Y.
        • et al.
        Behavioral and neural stability of attention bias to threat in healthy adolescents.
        Neuroimage. 2016; 136: 84-93
        • Auday E.S.
        • Taber-Thomas B.C.
        • Pérez-Edgar K.E.
        Neural correlates of attention bias to masked facial threat cues: Examining children at-risk for social anxiety disorder.
        Neuroimage Clin. 2018; 19: 202-212
        • White L.K.
        • Sequeira S.
        • Britton J.C.
        • Brotman M.A.
        • Gold A.L.
        • Berman E.
        • et al.
        Complementary features of attention bias modification therapy and cognitive-behavioral therapy in pediatric anxiety disorders [published correction appears in Am J Psychiatry 2018;175:83].
        Am J Psychiatry. 2017; 174: 775-784
        • Evans T.C.
        • Bar-Haim Y.
        • Fox N.A.
        • Pine D.S.
        • Britton J.C.
        Neural mechanisms underlying heterogeneous expression of threat-related attention in social anxiety.
        Behav Res Ther. 2020; 132103657
        • Kessler R.C.
        • Berglund P.
        • Demler O.
        • Jin R.
        • Merikangas K.R.
        • Walters E.E.
        Lifetime prevalence and age-of-onset distributions of DSM-IV disorders in the national comorbidity survey replication [published correction appears in Arch Gen Psychiatry 2005;62:768].
        Arch Gen Psychiatry. 2005; 62: 593-602
        • Gee D.G.
        • Humphreys K.L.
        • Flannery J.
        • Goff B.
        • Telzer E.H.
        • Shapiro M.
        • et al.
        A developmental shift from positive to negative connectivity in human amygdala-prefrontal circuitry.
        J Neurosci. 2013; 33: 4584-4593
        • Casey B.J.
        • Glatt C.E.
        • Lee F.S.
        Treating the developing versus developed brain: Translating preclinical mouse and human studies.
        Neuron. 2015; 86: 1358-1368
        • Hane A.A.
        • Fox N.A.
        • Henderson H.A.
        • Marshall P.J.
        Behavioral reactivity and approach-withdrawal bias in infancy.
        Dev Psychol. 2008; 44: 1491-1496
        • Jarcho J.M.
        • Davis M.M.
        • Shechner T.
        • Degnan K.A.
        • Henderson H.A.
        • Stoddard J.
        • et al.
        Early-childhood social reticence predicts brain function in preadolescent youths during distinct forms of peer evaluation.
        Psychol Sci. 2016; 27: 821-835
        • Michalska K.J.
        • Feldman J.S.
        • Ivie E.J.
        • Shechner T.
        • Sequeira S.
        • Averbeck B.
        • et al.
        Early-childhood social reticence predicts SCR-BOLD coupling during fear extinction recall in preadolescent youth.
        Dev Cogn Neurosci. 2019; 36100605
        • Donders A.R.T.
        • van der Heijden G.J.
        • Stijnen T.
        • Moons K.G.
        Review: A gentle introduction to imputation of missing values.
        J Clin Epidemiol. 2006; 59: 1087-1091
        • Chen G.
        • Saad Z.S.
        • Britton J.C.
        • Pine D.S.
        • Cox R.W.
        Linear mixed-effects modeling approach to FMRI group analysis.
        Neuroimage. 2013; 73: 176-190
        • Matta T.H.
        • Flournoy J.C.
        • Byrne M.L.
        Making an unknown unknown a known unknown: Missing data in longitudinal neuroimaging studies.
        Dev Cogn Neurosci. 2018; 33: 83-98
        • Degnan K.A.
        • Hane A.A.
        • Henderson H.A.
        • Moas O.L.
        • Reeb-Sutherland B.C.
        • Fox N.A.
        Longitudinal stability of temperamental exuberance and social-emotional outcomes in early childhood.
        Dev Psychol. 2011; 47: 765-780
        • Birmaher B.
        • Khetarpal S.
        • Brent D.
        • Cully M.
        • Balach L.
        • Kaufman J.
        • Neer S.M.
        The screen for child anxiety related emotional disorders (SCARED): Scale construction and psychometric characteristics.
        J Am Acad Child Adolesc Psychiatry. 1997; 36: 545-553
        • Muris P.
        • Dreessen L.
        • Bögels S.
        • Weckx M.
        • van Melick M.
        A questionnaire for screening a broad range of DSM-defined anxiety disorder symptoms in clinically referred children and adolescents.
        J Child Psychol Psychiatry. 2004; 45: 813-820
        • Rappaport B.I.
        • Pagliaccio D.
        • Pine D.S.
        • Klein D.N.
        • Jarcho J.M.
        Discriminant validity, diagnostic utility, and parent-child agreement on the screen for child anxiety related emotional disorders (SCARED) in treatment- and non-treatment-seeking youth.
        J Anxiety Disord. 2017; 51: 22-31
        • Behrens B.
        • Swetlitz C.
        • Pine D.S.
        • Pagliaccio D.
        The Screen for Child Anxiety Related Emotional Disorders (SCARED): Informant discrepancy, measurement invariance, and test-retest reliability.
        Child Psychiatry Hum Dev. 2019; 50: 473-482
        • Macleod C.
        • Mathews A.
        • Tata P.
        Attentional bias in emotional disorders.
        J Abnorm Psychol. 1986; 95: 15-20
        • Tottenham N.
        • Tanaka J.W.
        • Leon A.C.
        • McCarry T.
        • Nurse M.
        • Hare T.A.
        • et al.
        The NimStim set of facial expressions: Judgments from untrained research participants.
        Psychiatry Res. 2009; 168: 242-249
        • Shechner T.
        • Bar-Haim Y.
        Threat monitoring and attention-bias modification in anxiety and stress-related disorders.
        Curr Dir Psychol Sci. 2016; 25: 431-437
        • Badura-Brack A.S.
        • Naim R.
        • Ryan T.J.
        • Levy O.
        • Abend R.
        • Khanna M.M.
        • et al.
        Effect of attention training on attention bias variability and PTSD symptoms: Randomized controlled trials in Israeli and U.S. combat veterans.
        Am J Psychiatry. 2015; 172: 1233-1241
        • Iacoviello B.M.
        • Wu G.
        • Abend R.
        • Murrough J.W.
        • Feder A.
        • Fruchter E.
        • et al.
        Attention bias variability and symptoms of posttraumatic stress disorder.
        J Trauma Stress. 2014; 27: 232-239
        • Naim R.
        • Abend R.
        • Wald I.
        • Eldar S.
        • Levi O.
        • Fruchter E.
        • et al.
        Threat-related attention bias variability and posttraumatic stress.
        Am J Psychiatry. 2015; 172: 1242-1250
        • Pinheiro J.
        • Bates D.
        • DebRoy S.
        • Sarkar D.
        • Heisterkamp S.
        • et al.
        • EISPACK Authors
        Package ‘nlme’: Linear and nonlinear mixed effects models.
        (Available at:)
        • De Rosario-Marinez H.
        • Fox J.
        • CoreTeam
        Package ‘phia’: Post-hoc interaction analysis..
        (Available at:)
        • Lenth R.V.
        • Buerkner P.
        • Giné-Vázquez I.
        • Herve M.
        • Jung M.
        • Love J.
        • et al.
        emmeans: Estimated marginal means, aka least-squares means..
        (Available at:)
        • Cox R.W.
        AFNI: Software for analysis and visualization of functional magnetic resonance neuroimages.
        Comput Biomed Res. 1996; 29: 162-173
        • Cox R.W.
        • Hyde J.S.
        Software tools for analysis and visualization of fMRI data.
        Nmr Biomed. 1997; 10: 171-178
        • McLaren D.G.
        • Ries M.L.
        • Xu G.
        • Johnson S.C.
        A generalized form of context-dependent psychophysiological interactions (gPPI): A comparison to standard approaches.
        Neuroimage. 2012; 61: 1277-1286
        • Fischl B.
        • Salat D.H.
        • Busa E.
        • Albert M.
        • Dieterich M.
        • Haselgrove C.
        • et al.
        Whole brain segmentation: Automated labeling of neuroanatomical structures in the human brain.
        Neuron. 2002; 33: 341-355
        • Cox R.W.
        • Chen G.
        • Glen D.R.
        • Reynolds R.C.
        • Taylor P.A.
        FMRI clustering in AFNI: False-positive rates redux.
        Brain Connect. 2017; 7: 152-171
        • Eklund A.
        • Nichols T.E.
        • Knutsson H.
        Cluster failure: Why fMRI inferences for spatial extent have inflated false-positive rates [published correction appears in Proc Natl Acad Sci U S A 2016;113:E4929].
        Proc Natl Acad Sci U S A. 2016; 113: 7900-7905
        • Ocklenburg S.
        • Peterburs J.
        • Mundorf A.
        Hemispheric asymmetries in the amygdala: A comparative primer.
        Prog Neurobiol. 2022; 214102283
        • McMenamin B.W.
        • Marsolek C.J.
        Can theories of visual representation help to explain asymmetries in amygdala function?.
        Cogn Affect Behav Neurosci. 2013; 13: 211-224
        • Haller S.P.
        • Chen G.
        • Kitt E.R.
        • Smith A.R.
        • Stoddard J.
        • Abend R.
        • et al.
        Reliability of task-evoked neural activation during face-emotion paradigms: Effects of scanner and psychological processes.
        Hum Brain Mapp. 2022; 43: 2109-2120
        • Birn R.M.
        • Shackman A.J.
        • Oler J.A.
        • Williams L.E.
        • McFarlin D.R.
        • Rogers G.M.
        • et al.
        Evolutionarily conserved prefrontal-amygdalar dysfunction in early-life anxiety.
        Mol Psychiatry. 2014; 19: 915-922
        • Finn E.S.
        • Shen X.
        • Scheinost D.
        • Rosenberg M.D.
        • Huang J.
        • Chun M.M.
        • et al.
        Functional connectome fingerprinting: Identifying individuals using patterns of brain connectivity.
        Nat Neurosci. 2015; 18: 1664-1671
        • Greene A.S.
        • Gao S.
        • Scheinost D.
        • Constable R.T.
        Task-induced brain state manipulation improves prediction of individual traits.
        Nat Commun. 2018; 9: 2807
        • Schultz D.H.
        • Cole M.W.
        Higher intelligence is associated with less task-related brain network reconfiguration.
        J Neurosci. 2016; 36: 8551-8561
        • Harrewijn A.
        • Abend R.
        • Linke J.
        • Brotman M.A.
        • Fox N.A.
        • Leibenluft E.
        • et al.
        Combining fMRI during resting state and an attention bias task in children.
        Neuroimage. 2020; 205116301
        • MacLeod C.
        • Grafton B.
        • Notebaert L.
        Anxiety-linked attentional bias: Is it reliable?.
        Annual Review of Clinical Psychology. 2019; 15: 529-554
        • Rodebaugh T.L.
        • Scullin R.B.
        • Langer J.K.
        • Dixon D.J.
        • Huppert J.D.
        • Bernstein A.
        • et al.
        Unreliability as a threat to understanding psychopathology: The cautionary tale of attentional bias [published correction appears in J Abnorm Psychol 2016;125:878].
        J Abnorm Psychol. 2016; 125: 840-851
        • Price R.B.
        • Kuckertz J.M.
        • Siegle G.J.
        • Ladouceur C.D.
        • Silk J.S.
        • Ryan N.D.
        • et al.
        Empirical recommendations for improving the stability of the dot-probe task in clinical research [published correction appears in Psychol Assess 2015;27:376].
        Psychol Assess. 2015; 27: 365-376
        • Price R.B.
        • Brown V.
        • Siegle G.J.
        Computational modeling applied to the dot-probe task yields improved reliability and mechanistic insights.
        Biol Psychiatry. 2019; 85: 606-612
        • Price R.B.
        • Allen K.B.
        • Silk J.S.
        • Ladouceur C.D.
        • Ryan N.D.
        • Dahl R.E.
        • et al.
        Vigilance in the laboratory predicts avoidance in the real world: A dimensional analysis of neural, behavioral, and ecological momentary data in anxious youth.
        Dev Cogn Neurosci. 2016; 19: 128-136
        • Sequeira S.L.
        • Rosen D.K.
        • Silk J.S.
        • Hutchinson E.
        • Allen K.B.
        • Jones N.P.
        • et al.
        Don’t judge me!”: Links between in vivo attention bias toward a potentially critical judge and fronto-amygdala functional connectivity during rejection in adolescent girls.
        Dev Cogn Neurosci. 2021; 49100960