Advertisement
Archival Report|Articles in Press

Effects of the COVID-19 Pandemic on Mental Health and Brain Maturation in Adolescents: Implications for Analyzing Longitudinal Data

Open AccessPublished:December 01, 2022DOI:https://doi.org/10.1016/j.bpsgos.2022.11.002

      Abstract

      Background

      The COVID-19 pandemic has caused significant stress and disruption for young people, likely leading to alterations in their mental health and neurodevelopment. In this context, it is not clear whether youths who lived through the pandemic and its shutdowns are comparable psychobiologically to their age- and sex-matched peers assessed before the pandemic. This question is particularly important for researchers who are analyzing longitudinal data that span the pandemic.

      Methods

      We compared carefully matched youths assessed before the pandemic (n = 81) and after the pandemic-related shutdowns ended (n = 82).

      Results

      We found that youths assessed after the pandemic shutdowns had more severe internalizing mental health problems, reduced cortical thickness, larger hippocampal and amygdala volume, and more advanced brain age.

      Conclusions

      The COVID-19 pandemic not only appears to have led to poorer mental health and accelerated brain aging in adolescents, but it also poses significant challenges to researchers analyzing data from longitudinal studies of normative development that were interrupted by the pandemic.

      Keywords

      The COVID-19 pandemic has been a generation-defining event and a major source of adversity. Given the shelter-in-place orders in spring 2020 that led to school closures, academic disruptions, social restrictions, and reduced access to school-based mental health services (
      • Golberstein E.
      • Wen H.
      • Miller B.F.
      Coronavirus disease 2019 (COVID-19) and mental health for children and adolescents.
      ), the pandemic appears to have been particularly difficult for children and adolescents (
      • Jiao W.Y.
      • Wang L.N.
      • Liu J.
      • Fang S.F.
      • Jiao F.Y.
      • Pettoello-Mantovani M.
      • Somekh E.
      Behavioral and emotional disorders in children during the COVID-19 epidemic.
      ,
      • Xie X.
      • Xue Q.
      • Zhou Y.
      • Zhu K.
      • Liu Q.
      • Zhang J.
      • Song R.X.
      Mental health status among children in home confinement during the coronavirus disease 2019 outbreak in Hubei Province, China.
      ,
      • Zhou S.J.
      • Zhang L.G.
      • Wang L.L.
      • Guo Z.C.
      • Wang J.Q.
      • Chen J.C.
      • et al.
      Prevalence and socio-demographic correlates of psychological health problems in Chinese adolescents during the outbreak of COVID-19.
      ). In fact, a recent meta-analysis found that the prevalence of internalizing symptoms in youths has doubled since the onset of the COVID-19 pandemic (
      • Racine N.
      • McArthur B.A.
      • Cooke J.E.
      • Eirich R.
      • Zhu J.
      • Madigan S.N.
      Global prevalence of depressive and anxiety symptoms in children and adolescents during COVID-19: A meta-analysis.
      ). Despite this alarming statistic, however, the potential implications of the pandemic for children’s neurodevelopment have not been delineated.
      Research conducted before the pandemic found that exposure to early life adversity, including violence, neglect, and family dysfunction, is associated not only with poorer mental health, but also with maladaptive neurodevelopmental outcomes that indicate accelerated brain maturation or aging (
      • Colich N.L.
      • Rosen M.L.
      • Williams E.S.
      • McLaughlin K.A.
      Biological aging in childhood and adolescence following experiences of threat and deprivation: A systematic review and meta-analysis.
      ). For example, cortical thickness, which decreases with age (
      • Wierenga L.M.
      • Langen M.
      • Oranje B.
      • Durston S.
      Unique developmental trajectories of cortical thickness and surface area.
      ), is further reduced in youths with a history of early adversity (
      • Colich N.L.
      • Rosen M.L.
      • Williams E.S.
      • McLaughlin K.A.
      Biological aging in childhood and adolescence following experiences of threat and deprivation: A systematic review and meta-analysis.
      ). Recently, researchers have used machine learning algorithms to predict individuals’ ages from their neuroanatomical features (
      • Cole J.H.
      • Franke K.
      Predicting age using neuroimaging: Innovative brain ageing biomarkers.
      ). In adolescents, exposure to adversity has been associated with a brain age gap estimation (BrainAGE) suggestive of accelerated aging (i.e., having a predicted brain age older than one’s chronological age) (
      • Drobinin V.
      • Van Gestel H.
      • Helmick C.A.
      • Schmidt M.H.
      • Bowen C.V.
      • Uher R.
      The developmental brain age is associated with adversity, depression, and functional outcomes among adolescents.
      ). As a result of social isolation and distancing during the shutdown, virtually all youth experienced adversity in the form of significant departures from their normal routines. In addition, financial strain, threats to physical health, and exposure to increased familial violence were alarmingly common during the pandemic (
      • Prime H.
      • Wade M.
      • Browne D.T.
      Risk and resilience in family well-being during the COVID-19 pandemic.
      ,
      • Francisco R.
      • Pedro M.
      • Delvecchio E.
      • Espada J.P.
      • Morales A.
      • Mazzeschi C.
      • Orgilés M.
      Psychological symptoms and behavioral changes in children and adolescents during the early phase of COVID-19 quarantine in three European countries.
      ). If the pandemic has adversely affected adolescents’ mental health and neurodevelopment, such that adolescents who are assessed now differ in significant respects from their age- and sex-matched peers who were assessed before the pandemic, researchers must give serious consideration to how they accurately analyze and interpret longitudinal developmental data that span years on both sides of this extraordinary event.
      In this study, we matched a group of adolescents who experienced the pandemic shutdown (peri-COVID group) with a group of adolescents, matched on age, sex, puberty, exposure to early life stress, and socioeconomic status, who underwent the same assessment before the pandemic (pre-COVID group). We expected that compared with the pre-COVID group, the peri-COVID group would report more severe mental health problems and have older, or more mature, brains.

      Methods and Materials

      Participants

      Participants in this study were 163 adolescents (103 females) living in the San Francisco Bay Area who were participating in a larger longitudinal study assessing the effects of early life stress on psychobiology across puberty (N = 214) (
      • Chahal R.
      • Miller J.G.
      • Yuan J.P.
      • Buthmann J.L.
      • Gotlib I.H.
      An exploration of dimensions of early adversity and the development of functional brain network connectivity during adolescence: Implications for trajectories of internalizing symptoms.
      ,
      • King L.S.
      • Humphreys K.L.
      • Camacho M.C.
      • Gotlib I.H.
      A person-centered approach to the assessment of early life stress: Associations with the volume of stress-sensitive brain regions in early adolescence.
      ,
      • Chahal R.
      • Ho T.C.
      • Miller J.G.
      • Borchers L.R.
      • Gotlib I.H.
      Sex-specific vulnerability to depressive symptoms across adolescence and during the COVID-19 pandemic: The role of the cingulum bundle.
      ). Exclusion criteria were postpubertal status, nonfluency in English, inability to undergo magnetic resonance imaging (MRI), and history of neurological disorder or major medical illness. Participants were invited to return for follow-up assessments approximately every 2 years; however, the approximately 1-year-long COVID-19 pandemic shutdown beginning in March 2020 interrupted participants’ in-person assessments [see (
      • Miller J.G.
      • Ho T.C.
      • Kirshenbaum J.S.
      • Chahal R.
      • Gifuni A.J.
      • Gotlib I.H.
      Testing a developmental model of positive parenting, amygdala–subgenual anterior cingulate cortex connectivity, and depressive symptoms in adolescents before and during the COVID-19 pandemic.
      ) for more details]. All participants and their legal guardians gave informed assent and consent, respectively, and were compensated for their time. All study procedures were approved by the Stanford University Institutional Review Board.
      From this larger cohort, we constructed 2 matched subsamples using data collected either before the pandemic (from November 2016 to November 2019; pre-COVID group, n = 81) or during the pandemic but following the end of the Bay Area shutdown (from October 2020 to March 2022; peri-COVID group, n = 82). We constructed these subsamples to maximize group sizes and to match the 2 groups on sex, age, pubertal status, race/ethnicity, parental education, annual household income, and severity of early life stress based on panel ratings of participants’ responses to interview (
      • Chahal R.
      • Miller J.G.
      • Yuan J.P.
      • Buthmann J.L.
      • Gotlib I.H.
      An exploration of dimensions of early adversity and the development of functional brain network connectivity during adolescence: Implications for trajectories of internalizing symptoms.
      ,
      • King L.S.
      • Humphreys K.L.
      • Camacho M.C.
      • Gotlib I.H.
      A person-centered approach to the assessment of early life stress: Associations with the volume of stress-sensitive brain regions in early adolescence.
      ). Specifically, we attempted to match the peri-COVID participants with pre-COVID participants with respect to age and sex as closely as possible at the group level. Not all peri-COVID participants could be matched to pre-COVID participants given their older age, and not all pre-COVID participants were needed to be matched to the smaller peri-COVID group (and the youngest pre-COVID participants were too young to be matched to the peri-COVID participants). For the mental health symptoms sample of 81 pre-COVID and 82 peri-COVID participants (see below), we excluded from analyses 50 pre-COVID and 12 peri-COVID participants who could not be appropriately matched. For the neuroimaging sample, we were able to match by age and sex 64 of the 104 participants who underwent scanning in the peri-COVID period with 64 pre-COVID participants.

      Mental Health Symptoms

      Participants self-reported their depressive symptoms using the 10-item version of the Children’s Depression Inventory (
      • Kovacs M.
      Children’s Depression Inventory (CDI): Technical Manual Update.
      ). This widely used reliable measure (
      • Saylor C.F.
      • Finch A.J.
      • Spirito A.
      • Bennett B.
      The Children’s Depression Inventory: A systematic evaluation of psychometric properties.
      ) has been shown to have convergent validity with clinician ratings of depression symptoms and diagnosis (
      • Timbremont B.
      • Braet C.
      • Dreessen L.
      Assessing depression in youth: Relation between the Children’s Depression Inventory and a structured interview.
      ). We assessed anxiety symptoms using total score of the Social Anxiety and Physical Symptoms subscales of the Multidimensional Anxiety Scale for Children (MASC) (
      • March J.S.
      • Parker J.D.
      • Sullivan K.
      • Stallings P.
      • Conners C.K.
      The Multidimensional Anxiety Scale for Children (MASC): Factor structure, reliability, and validity.
      ). The full MASC assesses a wide range of anxiety symptoms, including symptoms that are not as relevant for the age range of our participants (e.g., separation anxiety). For this reason and to reduce participant burden, we administered only the Social Anxiety and Physical Symptoms subscales of the MASC for this study; therefore, the MASC total score in this study reflects the sum of these 2 subscales. Finally, we assessed internalizing and externalizing symptoms using the validated subscales of the Youth Self-Report version of the Child Behavior Checklist (
      • Achenbach T.M.
      Manual for the Youth Self-Report and 1991 Profile.
      ).

      Neuroimaging

      A subset of the participants (matched n = 64 per group) completed a T1-weighted MRI scan at the Center for Cognitive and Neurobiological Imaging at Stanford University. All participants in the pre-COVID group completed their scans using a 3T Discovery MR750 (GE Medical Systems). As of March 16, 2020, the Discovery MR750 was upgraded to an ultra-high performance system. Thus, all peri-COVID participants were scanned on the upgraded scanner. Participants in both groups were scanned using a 32-channel head coil (Nova Medical, Inc.). Prior work suggests that FreeSurfer-based cortical thickness and subcortical measures are highly reliable across scanner upgrades (
      • Jovicich J.
      • Czanner S.
      • Han X.
      • Salat D.
      • van der Kouwe A.
      • Quinn B.
      • et al.
      MRI-derived measurements of human subcortical, ventricular and intracranial brain volumes: Reliability effects of scan sessions, acquisition sequences, data analyses, scanner upgrade, scanner vendors and field strengths.
      ,
      • Han X.
      • Jovicich J.
      • Salat D.
      • van der Kouwe A.
      • Quinn B.
      • Czanner S.
      • et al.
      Reliability of MRI-derived measurements of human cerebral cortical thickness: The effects of field strength, scanner upgrade and manufacturer.
      ,
      • Brown E.M.
      • Pierce M.E.
      • Clark D.C.
      • Fischl B.R.
      • Iglesias J.E.
      • Milberg W.P.
      • et al.
      Test-retest reliability of FreeSurfer automated hippocampal subfield segmentation within and across scanners.
      ). For example, Han et al. (
      • Han X.
      • Jovicich J.
      • Salat D.
      • van der Kouwe A.
      • Quinn B.
      • Czanner S.
      • et al.
      Reliability of MRI-derived measurements of human cerebral cortical thickness: The effects of field strength, scanner upgrade and manufacturer.
      ) did not find evidence that scanner upgrades introduce bias for cortical thickness measures, and Brown et al. (
      • Brown E.M.
      • Pierce M.E.
      • Clark D.C.
      • Fischl B.R.
      • Iglesias J.E.
      • Milberg W.P.
      • et al.
      Test-retest reliability of FreeSurfer automated hippocampal subfield segmentation within and across scanners.
      ) found that hippocampal measures are reliable across scanners. In addition, we conducted analyses with our own data to assess potential differences in T1-weighted image quality related to scanner upgrade. Specifically, in a subset of 31 participants with imaging data before and after the scanner upgrade, we tested within-participant changes in gray–white matter contrast-to-noise ratio using FreeSurfer’s mri_cnr quality metric command. We did not find significant differences in contrast-to-noise ratio from pre- to post-upgrade in either the left (t30 = 0.81, p = .425) or the right (t30 = 0.66, p = .513) hemisphere. Thus, the scanner upgrade does not appear to have introduced a systematic bias in image quality. Whole-brain T1-weighted images were collected for all participants using the following spoiled gradient echo pulse sequence: 186 sagittal slices, repetition time/echo time/inversion time = 6.24/2.34/450 ms, flip angle = 12°, voxel size = 0.9 mm × 0.9 mm × 0.9 mm, scan duration = 315 seconds. The spoiled gradient echo sequence was repeated up to 2 additional times if the first acquisition did not yield clear images. For each participant with multiple acquisitions, the single spoiled gradient echo image with the clearest structural boundaries (i.e., that was free from motion or other artifacts) was used for further analysis.

      Segmentation of Cortical and Subcortical Regions

      We used FreeSurfer version 6.0 (https://surfer.nmr.mgh.harvard.edu/) recon-all function to automatically skull strip and segment cortical and subcortical volumes from the T1-weighted structural images (
      • 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.
      ), which has been shown to have acceptable scan-rescan reliability (
      • Jovicich J.
      • Czanner S.
      • Han X.
      • Salat D.
      • van der Kouwe A.
      • Quinn B.
      • et al.
      MRI-derived measurements of human subcortical, ventricular and intracranial brain volumes: Reliability effects of scan sessions, acquisition sequences, data analyses, scanner upgrade, scanner vendors and field strengths.
      ) and comparable accuracy to manual labeling techniques (
      • 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.
      ,
      • Fischl B.
      • Dale A.M.
      Measuring the thickness of the human cerebral cortex from magnetic resonance images.
      ,
      • Morey R.A.
      • Petty C.M.
      • Xu Y.
      • Hayes J.P.
      • Wagner 2nd, H.R.
      • Lewis D.V.
      • et al.
      A comparison of automated segmentation and manual tracing for quantifying hippocampal and amygdala volumes.
      ). We implemented structural image processing protocols established by ENIGMA (https://enigma.ini.usc.edu/protocols/imaging-protocols/) to extract and perform quality assurance checks on the cortical thickness and subcortical volume estimates from the FreeSurfer outputs. Using FreeView image viewer, all cortical and subcortical outputs were visually inspected to quality check for processing and segmentation errors. As previously described (
      • Ho T.C.
      • Teresi G.I.
      • Ojha A.
      • Walker J.C.
      • Kirshenbaum J.S.
      • Singh M.K.
      • Gotlib I.H.
      Smaller caudate gray matter volume is associated with greater implicit suicidal ideation in depressed adolescents.
      ,
      • Ho T.C.
      • Cichocki A.C.
      • Gifuni A.J.
      • Camacho M.C.
      • Ordaz S.J.
      • Singh M.K.
      • Gotlib I.H.
      Reduced dorsal striatal gray matter volume predicts implicit suicidal ideation in adolescents.
      ), we converted gray matter volumes from each hemisphere into z scores; volumes with z scores greater than 2.5 or less than −2.5 were visually examined again for accuracy, and any segmentations that failed any of these steps were removed from final analyses. We focused on mean cortical thickness (average of cortical thickness values across individual regions as defined by the Desikan-Killiany atlas) (
      • Desikan R.S.
      • Ségonne F.
      • Fischl B.
      • Quinn B.T.
      • Dickerson B.C.
      • Blacker D.
      • et al.
      An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest.
      ) and unstandardized residuals of subcortical volumes regressed on total intracranial volume.

      Brain Age Gap Estimates

      Based on cortical and subcortical features, we computed BrainAGE values for male and female participants using sex-specific machine learning–based models developed by the ENIGMA-Brain Age working group (
      • Han L.K.
      • Dinga R.
      • Hahn T.
      • Ching C.R.
      • Eyler L.T.
      • Aftanas L.
      • et al.
      Brain aging in major depressive disorder: Results from the ENIGMA major depressive disorder working group.
      ). These models use data from 14 subcortical gray matter volumes, 2 lateral ventricles, 68 cortical thickness measures, 68 surface area measures, and total intracranial volume to predict chronological age (i.e., predicted brain age). We computed brain age gap estimates by subtracting chronological age from predicted brain age. Given that BrainAGE values are often overestimated in younger individuals and underestimated in older individuals, Le et al. (
      • Le T.T.
      • Kuplicki R.T.
      • McKinney B.A.
      • Yeh H.W.
      • Thompson W.K.
      • Paulus M.P.
      Tulsa 1000 Investigators
      A nonlinear simulation framework supports adjusting for age when analyzing BrainAGE.
      ) proposed adjusting for chronological age in analyses of BrainAGE. Therefore, we regressed gap estimates onto chronological age and used the unstandardized residuals as the BrainAGE outcome variable in our statistical analyses.

      Statistical Approach

      All statistical analyses were conducted using R version 4.0.2 (

      R Core Team (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available at: https://www.R-project.org/. Accessed November 7, 2022.

      ). To examine whether adolescents who experienced the pandemic differed from their pre-pandemic peers, we conducted between-group tests on measures of internalizing and externalizing symptoms, cortical thickness, and subcortical volume (regions of interest were bilateral amygdala, hippocampus, and nucleus accumbens). For analyses of mental health problems, we first conducted a one-way multiple analysis of variance test to examine whether there were group differences in overall mental health scores across measures. We conducted follow-up independent sample t tests to examine whether the pre-COVID and peri-COVID groups differed in specific aspects of mental health as assessed by different measures. We repeated these steps for analyses of brain metrics. Given our expectations based on recent work suggesting that mental health problems have increased during the pandemic (
      • Racine N.
      • McArthur B.A.
      • Cooke J.E.
      • Eirich R.
      • Zhu J.
      • Madigan S.N.
      Global prevalence of depressive and anxiety symptoms in children and adolescents during COVID-19: A meta-analysis.
      ), we used one-tailed hypothesis tests for follow-up analyses of mental health outcomes. We used two-tailed hypothesis tests for follow-up analyses of brain outcomes.

      Results

      Participant Characteristics

      The demographic and clinical characteristics of the pre- and peri-COVID subgroups of participants are presented in Table 1. Participants’ parents reported on their annual household income, from which we computed an income-to-needs ratio by dividing the midpoint of their reported income bin by the low-income value for Santa Clara County. Importantly, this calculation considers the number of people in the home and the time period in which the study occurred [https://www.huduser.gov/portal/datasets/il/il2017/2017summary.odn; (
      • King L.S.
      • Dennis E.L.
      • Humphreys K.L.
      • Thompson P.M.
      • Gotlib I.H.
      Cross-sectional and longitudinal associations of family income-to-needs ratio with cortical and subcortical brain volume in adolescent boys and girls.
      )]. Attesting to the success of our careful matching procedure, there were no significant group differences in participant characteristics between the pre-COVID and peri-COVID subgroups for either the mental health or the brain samples (all individual ps > .06).
      Table 1Participant Characteristics
      Mental Health Subsample
      Subsample of participants who completed the measures of mental health.
      Brain Subsample
      Subsample of participants who also successfully completed the neuroimaging protocol.
      VariablePre-COVID, n = 81Peri-COVID, n = 82Pre-COVID, n = 64Peri-COVID, n = 64
      Sex, Female51 (63%)52 (63%)34 (53%)34 (53%)
      Age, Years15.87 (1.14)16.17 (0.93)16.08 (0.90)16.43 (1.19)
      Race/Ethnicity
       Asian/Asian American13 (16%)11 (13%)7 (11%)9 (14%)
       Biracial9 (11%)21 (26%)11 (17%)12 (19%)
       Black/African American6 (7%)7 (9%)3 (5%)7 (11%)
       Hispanic/Latinx5 (6%)7 (9%)4 (6%)6 (9%)
       Other
      Primarily participants who identified as multiracial.
      5 (6%)3 (4%)4 (6%)4 (6%)
       White43 (53%)33 (40%)35 (55%)26 (41%)
      Income-to-Needs Ratio1.37 (0.53)1.27 (0.54)1.30 (0.59)1.31 (0.53)
      Early Life Stress6.60 (4.89)6.31 (4.97)7.30 (5.88)5.84 (4.68)
      Parental Education
       No GED/high school diploma0 (0%)1 (1%)0 (0%)1 (2%)
       GED/high school diploma0 (0%)4 (5%)0 (0%)0 (0%)
       Some college10 (12%)15 (18%)4 (6%)10 (16%)
       2-Year college degree5 (6%)7 (9%)5 (8%)5 (8%)
       4-Year college degree29 (36%)25 (30%)24 (38%)27 (42%)
       Master’s degree30 (37%)20 (24%)20 (31%)17 (27%)
       Professional degree1 (1%)5 (6%)3 (5%)3 (5%)
       Doctorate4 (5%)2 (2%)2 (3%)0 (0%)
       Not reported2 (2%)3 (4%)6 (9%)1 (2%)
      COVID-19 Impact
       Individual diagnosisN/A1 (1%)N/A1 (1%)
       Household diagnosisN/A3 (4%)N/A2 (3%)
       Financial strainN/A13 (16%)N/A11 (17%)
       Job lossN/A7 (9%)N/A8 (13%)
      Values are mean (SD) or n (%).
      GED, General Educational Development; N/A, not applicable.
      a Subsample of participants who completed the measures of mental health.
      b Subsample of participants who also successfully completed the neuroimaging protocol.
      c Primarily participants who identified as multiracial.

      Mental Health

      Group differences on the mental health measures are presented in Figure 1. A one-way multivariate analysis of variance indicated that the pre-COVID (n = 81) and peri-COVID (n = 82) groups differed significantly in their self-reported mental health difficulties (F4,158 = 2.67, p = .034). Follow-up t tests showed that the peri-COVID group reported more severe symptoms of anxiety (t161 = 3.15, p < .001; Cohen’s d = 0.49), depression (t161 = 1.92, p = .029; d = 0.30), and internalizing problems (t161 = 1.77, p = .039; d = 0.28); the 2 groups did not differ in externalizing problems (t161 = 1.25, p = .108).
      Figure thumbnail gr1
      Figure 1Group differences on the Children’s Depression Inventory (CDI), Multidimensional Anxiety Scale for Children (MASC) (sum of the Social Anxiety and Physical Symptom subscales), and Youth Self-Report (YSR) internalizing and externalizing. ∗p < .05, ∗∗∗p < .001. n.s., not significant.

      Brain Metrics

      Group differences in cortical thickness, subcortical volumes, and BrainAGE are shown in Figure 2. A one-way multivariate analysis of variance conducted on all the brain metrics yielded a significant difference between the pre-COVID (n = 61–64) and peri-COVID (n = 63–64) groups (F5,116 = 7.13, p < .001). Follow-up tests indicated that the peri-COVID group had reduced bilateral cortical thickness (t122 = 3.67, p < .001; d = 0.66) and, controlling for intracranial volume, larger bilateral hippocampal volume (t125 = 3.56, p < .001; d = 0.63) and bilateral amygdala volume (t125 = 2.01, p = .047; d = 0.36); the 2 groups did not differ in bilateral nucleus accumbens volume (t125 = 0.68, p = .248; d = 0.12). Finally, despite the fact that the 2 groups were matched on age and other relevant demographic characteristics, adolescents in the peri-COVID group had an older BrainAGE than their peers who were assessed before the pandemic (t125 = 2.31, p = .022; d = 0.41).
      Figure thumbnail gr2
      Figure 2Raw data are plotted for visualization. Significance levels are based on group differences in subcortical volumes (in mm) adjusted for intracranial volume, cortical thickness, and brain age gap estimation (BrainAGE) adjusted for chronological age. ∗p < .05, ∗∗∗p < .001. n.s., not significant.

      Interval Between COVID-19 Shutdown and Peri-COVID Assessments

      Finally, given the possibility that participants’ mental health difficulties and their brain metrics increased with the duration of the pandemic, we examined our clinical functioning and brain metrics as a function of time since the Bay Area shelter-in-place orders were initiated (March 17, 2020). The peri-COVID participants completed measures of clinical functioning between January 10, 2021, and September 30, 2021, and MRI scans between October 13, 2020, and March 22, 2022. Within the peri-COVID group, we examined associations between the number of days from the start of shelter-in-place orders to the dates that participants completed measures of psychopathology (mean [SD] = 346.49 [131.70] days; range, 133–720 days). There were no significant associations between this interval and participants’ scores on the measures of depression (r80 = 0.01, p = .901), anxiety (r80 = −0.06, p = .544), internalizing symptoms (r80 = 0.07, p = .506), or externalizing symptoms (r80 = 0.00, p = .980). We repeated these analyses for the brain metrics (mean [SD] interval = 379.00 [119.24] days; range, 210–735 days). Again, there were no significant associations between the interval and residuals of amygdala volume (r62 = 0.01, p = .935), hippocampal volume (r62 = 0.15, p = .245), nucleus accumbens volume (r62 = 0.05, p = .681), mean cortical thickness (r61 = 0.04, p = .303), or residuals of BrainAGE (r62 = −0.09, p = .459).

      Discussion

      In addition to replicating prior findings that the pandemic has adversely affected the mental health of young people (
      • Racine N.
      • McArthur B.A.
      • Cooke J.E.
      • Eirich R.
      • Zhu J.
      • Madigan S.N.
      Global prevalence of depressive and anxiety symptoms in children and adolescents during COVID-19: A meta-analysis.
      ), we found that adolescents assessed during the pandemic have neuroanatomical features that are more typical of individuals who are older or who experienced significant adversity in childhood. Compared with carefully matched peers assessed before the pandemic, adolescents assessed during the pandemic showed signs of advanced cortical thinning and had larger bilateral hippocampal and amygdala volumes. Given that volume in these structures typically increases over adolescence (
      • Coupé P.
      • Catheline G.
      • Lanuza E.
      • Manjón J.V.
      Alzheimer’s Disease Neuroimaging Initiative
      Towards a unified analysis of brain maturation and aging across the entire lifespan: A MRI analysis.
      ), these neural alterations may reflect accelerated brain maturation in the context of the pandemic. Indeed, adolescents assessed during the pandemic also had larger positive brain age gap estimates, indicative of older-appearing brains.
      It appears, therefore, that the pandemic not only has adversely affected mental health of adolescents, but also has accelerated their brain maturation. These findings have critical implications for researchers who are conducting longitudinal studies that were interrupted due to pandemic-related shutdowns. In our own longitudinal study, we had been assessing a sample of approximately 200 adolescents at each of 4 time points, at 2-year intervals, to examine the effects of early adversity on trajectories of neurodevelopment and clinical symptoms. At the time of the shutdown, we were two thirds of the way through the third assessment, when our participants were 13 to 17 years of age. We had originally planned to simply use participants’ age in analyzing trajectories from our 4 time points of data. Although some participants would have had a longer interval than others between assessments that bracketed the shutdown, we would control statistically for those differences. It is important to recognize that this analytic approach assumes that, for example, 16-year-olds who were assessed after the shutdown ended are equivalent in their clinical functioning and neurodevelopment to 16-year-olds who were assessed before the pandemic and would simply be grouped together. Our results suggest that this assumption is not correct. Rather, the pandemic appears to have altered adolescent mental health and neurodevelopment, at least in the short term, which will present a challenge for researchers in analyzing longitudinal data from studies of normative development that were interrupted by the pandemic.
      In order to not confound age-related changes in brain maturation with experiences and consequences of the COVID-19 pandemic, some researchers, including our group, have used a dummy-coded variable to control statistically for whether participants were assessed before or during the pandemic [e.g. (
      • Ho T.C.
      • Kulla A.
      • Teresi G.I.
      • Sisk L.M.
      • Rosenberg-Hasson Y.
      • Maecker H.T.
      • Gotlib I.H.
      Inflammatory cytokines and callosal white matter microstructure in adolescents.
      )]. Nevertheless, restrictions around COVID-19 are constantly changing; therefore, additional measures may need to be used as covariates, including the interval between shelter-in-place orders and time of assessment as well as the nature and severity of the individual’s stress and experience during the pandemic (e.g., COVID-19 infection, upheaval in living situation, financial strain).
      We should note that our sample is of relatively high socioeconomic status and represents the racial/ethnic composition of the San Francisco Bay Area. Researchers have reported that sample composition influences age-related effects on brain structure (
      • LeWinn K.Z.
      • Sheridan M.A.
      • Keyes K.M.
      • Hamilton A.
      • McLaughlin K.A.
      Sample composition alters associations between age and brain structure.
      ) and, more specifically, that the psychosocial and health consequences of the pandemic have been more severe among individuals from socially marginalized groups [e.g., lower socioeconomic status (
      • Lee H.
      • Singh G.K.
      Monthly trends in self-reported health status and depression by race/ethnicity and socioeconomic status during the COVID-19 pandemic, United States, April 2020–May 2021.
      ,
      • Lopez 3rd, L.
      • Hart 3rd, L.H.
      • Katz M.H.
      Racial and ethnic health disparities related to COVID-19.
      ,
      • Mena G.E.
      • Martinez P.P.
      • Mahmud A.S.
      • Marquet P.A.
      • Buckee C.O.
      • Santillana M.
      Socioeconomic status determines COVID-19 incidence and related mortality in Santiago, Chile.
      )]. Therefore, it is important that investigators examine the effects of the COVID-19 pandemic on psychopathology and brain metrics in more diverse samples of adolescents that are representative of the broader population.
      Another critical task for future research is to determine whether these alterations are temporary effects of the pandemic or stable changes that will characterize the current generation of youths. If these changes are found to be enduring, accounting for and interpreting data acquired during this period will require additional attention and consideration. For example, as more researchers publish data concerning normative developmental trajectories of MRI-derived anatomical features [e.g. (
      • Ball G.
      • Beare R.
      • Seal M.L.
      Charting shared developmental trajectories of cortical thickness and structural connectivity in childhood and adolescence.
      )], it will be possible to compare COVID-19–impacted neurodevelopmental trajectories with normative trajectories and, indeed, to compute COVID-19–adjusted metrics of brain maturation. Regardless, however, we emphasize that it is important that we continue to follow and assess individuals who were recruited and assessed before the pandemic; this type of research offers the strongest possibility for us to examine the effects of a major stressor experienced on a global scale.

      Acknowledgments and Disclosures

      This work was supported by the National Institutes of Health (Grant Nos. R37MH101495 [to IHG] and K01MH117442 [to TCH]).
      The National Institutes of Health played no role in the study design, data collection, analysis, or preparation, review, or approval of the manuscript.
      IHG designed the study and helped to write the manuscript; JGM, LRB, SMC, and TCH analyzed the data and helped to write the manuscript; SMC, LAC, and JMG helped to collect the data and write the manuscript.
      We thank the youths and their families for participating in this research.
      The authors report no biomedical financial interests or potential conflicts of interest.

      Supplementary Material

      References

        • Golberstein E.
        • Wen H.
        • Miller B.F.
        Coronavirus disease 2019 (COVID-19) and mental health for children and adolescents.
        JAMA Pediatr. 2020; 174: 819-820
        • Jiao W.Y.
        • Wang L.N.
        • Liu J.
        • Fang S.F.
        • Jiao F.Y.
        • Pettoello-Mantovani M.
        • Somekh E.
        Behavioral and emotional disorders in children during the COVID-19 epidemic.
        J Pediatr. 2020; 221: 264-266.e1
        • Xie X.
        • Xue Q.
        • Zhou Y.
        • Zhu K.
        • Liu Q.
        • Zhang J.
        • Song R.X.
        Mental health status among children in home confinement during the coronavirus disease 2019 outbreak in Hubei Province, China.
        JAMA Pediatr. 2020; 174: 898-900
        • Zhou S.J.
        • Zhang L.G.
        • Wang L.L.
        • Guo Z.C.
        • Wang J.Q.
        • Chen J.C.
        • et al.
        Prevalence and socio-demographic correlates of psychological health problems in Chinese adolescents during the outbreak of COVID-19.
        Eur Child Adolesc Psychiatry. 2020; 29: 749-758
        • Racine N.
        • McArthur B.A.
        • Cooke J.E.
        • Eirich R.
        • Zhu J.
        • Madigan S.N.
        Global prevalence of depressive and anxiety symptoms in children and adolescents during COVID-19: A meta-analysis.
        JAMA Pediatr. 2021; 175: 1142-1150
        • Colich N.L.
        • Rosen M.L.
        • Williams E.S.
        • McLaughlin K.A.
        Biological aging in childhood and adolescence following experiences of threat and deprivation: A systematic review and meta-analysis.
        Psychol Bull. 2020; 146: 721-764
        • Wierenga L.M.
        • Langen M.
        • Oranje B.
        • Durston S.
        Unique developmental trajectories of cortical thickness and surface area.
        Neuroimage. 2014; 87: 120-126
        • Cole J.H.
        • Franke K.
        Predicting age using neuroimaging: Innovative brain ageing biomarkers.
        Trends Neurosci. 2017; 40: 681-690
        • Drobinin V.
        • Van Gestel H.
        • Helmick C.A.
        • Schmidt M.H.
        • Bowen C.V.
        • Uher R.
        The developmental brain age is associated with adversity, depression, and functional outcomes among adolescents.
        Biol Psychiatry Cogn Neurosci Neuroimaging. 2022; 7: 406-414
        • Prime H.
        • Wade M.
        • Browne D.T.
        Risk and resilience in family well-being during the COVID-19 pandemic.
        Am Psychol. 2020; 75: 631-643
        • Francisco R.
        • Pedro M.
        • Delvecchio E.
        • Espada J.P.
        • Morales A.
        • Mazzeschi C.
        • Orgilés M.
        Psychological symptoms and behavioral changes in children and adolescents during the early phase of COVID-19 quarantine in three European countries.
        Front Psychiatry. 2020; 11570164
        • Chahal R.
        • Miller J.G.
        • Yuan J.P.
        • Buthmann J.L.
        • Gotlib I.H.
        An exploration of dimensions of early adversity and the development of functional brain network connectivity during adolescence: Implications for trajectories of internalizing symptoms.
        Dev Psychopathol. 2022; 34: 557-571
        • King L.S.
        • Humphreys K.L.
        • Camacho M.C.
        • Gotlib I.H.
        A person-centered approach to the assessment of early life stress: Associations with the volume of stress-sensitive brain regions in early adolescence.
        Dev Psychopathol. 2019; 31: 643-655
        • Chahal R.
        • Ho T.C.
        • Miller J.G.
        • Borchers L.R.
        • Gotlib I.H.
        Sex-specific vulnerability to depressive symptoms across adolescence and during the COVID-19 pandemic: The role of the cingulum bundle.
        JCPP Adv. 2022; 2e12061
        • Miller J.G.
        • Ho T.C.
        • Kirshenbaum J.S.
        • Chahal R.
        • Gifuni A.J.
        • Gotlib I.H.
        Testing a developmental model of positive parenting, amygdala–subgenual anterior cingulate cortex connectivity, and depressive symptoms in adolescents before and during the COVID-19 pandemic.
        Biol Psychiatry Glob Open Sci. 2021; 1: 291-299
        • Kovacs M.
        Children’s Depression Inventory (CDI): Technical Manual Update.
        Multi-Health Systems, Inc, North Tonawanda, NY2003
        • Saylor C.F.
        • Finch A.J.
        • Spirito A.
        • Bennett B.
        The Children’s Depression Inventory: A systematic evaluation of psychometric properties.
        J Consult Clin Psychol. 1984; 52: 955-967
        • Timbremont B.
        • Braet C.
        • Dreessen L.
        Assessing depression in youth: Relation between the Children’s Depression Inventory and a structured interview.
        J Clin Child Adolesc Psychol. 2004; 33: 149-157
        • March J.S.
        • Parker J.D.
        • Sullivan K.
        • Stallings P.
        • Conners C.K.
        The Multidimensional Anxiety Scale for Children (MASC): Factor structure, reliability, and validity.
        J Am Acad Child Adolesc Psychiatry. 1997; 36: 554-565
        • Achenbach T.M.
        Manual for the Youth Self-Report and 1991 Profile.
        University of Vermont Department of Psychiatry, Burlington, VT1991
        • Jovicich J.
        • Czanner S.
        • Han X.
        • Salat D.
        • van der Kouwe A.
        • Quinn B.
        • et al.
        MRI-derived measurements of human subcortical, ventricular and intracranial brain volumes: Reliability effects of scan sessions, acquisition sequences, data analyses, scanner upgrade, scanner vendors and field strengths.
        Neuroimage. 2009; 46: 177-192
        • Han X.
        • Jovicich J.
        • Salat D.
        • van der Kouwe A.
        • Quinn B.
        • Czanner S.
        • et al.
        Reliability of MRI-derived measurements of human cerebral cortical thickness: The effects of field strength, scanner upgrade and manufacturer.
        Neuroimage. 2006; 32: 180-194
        • Brown E.M.
        • Pierce M.E.
        • Clark D.C.
        • Fischl B.R.
        • Iglesias J.E.
        • Milberg W.P.
        • et al.
        Test-retest reliability of FreeSurfer automated hippocampal subfield segmentation within and across scanners.
        Neuroimage. 2020; 210116563
        • 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
        • Fischl B.
        • Dale A.M.
        Measuring the thickness of the human cerebral cortex from magnetic resonance images.
        Proc Natl Acad Sci U S A. 2000; 97: 11050-11055
        • Morey R.A.
        • Petty C.M.
        • Xu Y.
        • Hayes J.P.
        • Wagner 2nd, H.R.
        • Lewis D.V.
        • et al.
        A comparison of automated segmentation and manual tracing for quantifying hippocampal and amygdala volumes.
        Neuroimage. 2009; 45: 855-866
        • Ho T.C.
        • Teresi G.I.
        • Ojha A.
        • Walker J.C.
        • Kirshenbaum J.S.
        • Singh M.K.
        • Gotlib I.H.
        Smaller caudate gray matter volume is associated with greater implicit suicidal ideation in depressed adolescents.
        J Affect Disord. 2021; 278: 650-657
        • Ho T.C.
        • Cichocki A.C.
        • Gifuni A.J.
        • Camacho M.C.
        • Ordaz S.J.
        • Singh M.K.
        • Gotlib I.H.
        Reduced dorsal striatal gray matter volume predicts implicit suicidal ideation in adolescents.
        Soc Cogn Affect Neurosci. 2018; 13: 1215-1224
        • Desikan R.S.
        • Ségonne F.
        • Fischl B.
        • Quinn B.T.
        • Dickerson B.C.
        • Blacker D.
        • et al.
        An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest.
        Neuroimage. 2006; 31: 968-980
        • Han L.K.
        • Dinga R.
        • Hahn T.
        • Ching C.R.
        • Eyler L.T.
        • Aftanas L.
        • et al.
        Brain aging in major depressive disorder: Results from the ENIGMA major depressive disorder working group.
        Mol Psychiatry. 2021; 26: 5124-5139
        • Le T.T.
        • Kuplicki R.T.
        • McKinney B.A.
        • Yeh H.W.
        • Thompson W.K.
        • Paulus M.P.
        • Tulsa 1000 Investigators
        A nonlinear simulation framework supports adjusting for age when analyzing BrainAGE.
        Front Aging Neurosci. 2018; 10317
      1. R Core Team (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available at: https://www.R-project.org/. Accessed November 7, 2022.

        • King L.S.
        • Dennis E.L.
        • Humphreys K.L.
        • Thompson P.M.
        • Gotlib I.H.
        Cross-sectional and longitudinal associations of family income-to-needs ratio with cortical and subcortical brain volume in adolescent boys and girls.
        Dev Cogn Neurosci. 2020; 44100796
        • Coupé P.
        • Catheline G.
        • Lanuza E.
        • Manjón J.V.
        • Alzheimer’s Disease Neuroimaging Initiative
        Towards a unified analysis of brain maturation and aging across the entire lifespan: A MRI analysis.
        Hum Brain Mapp. 2017; 38: 5501-5518
        • Ho T.C.
        • Kulla A.
        • Teresi G.I.
        • Sisk L.M.
        • Rosenberg-Hasson Y.
        • Maecker H.T.
        • Gotlib I.H.
        Inflammatory cytokines and callosal white matter microstructure in adolescents.
        Brain Behav Immun. 2022; 100: 321-331
        • LeWinn K.Z.
        • Sheridan M.A.
        • Keyes K.M.
        • Hamilton A.
        • McLaughlin K.A.
        Sample composition alters associations between age and brain structure.
        Nat Commun. 2017; 8: 874
        • Lee H.
        • Singh G.K.
        Monthly trends in self-reported health status and depression by race/ethnicity and socioeconomic status during the COVID-19 pandemic, United States, April 2020–May 2021.
        Ann Epidemiol. 2021; 63: 52-62
        • Lopez 3rd, L.
        • Hart 3rd, L.H.
        • Katz M.H.
        Racial and ethnic health disparities related to COVID-19.
        JAMA. 2021; 325: 719-720
        • Mena G.E.
        • Martinez P.P.
        • Mahmud A.S.
        • Marquet P.A.
        • Buckee C.O.
        • Santillana M.
        Socioeconomic status determines COVID-19 incidence and related mortality in Santiago, Chile.
        Science. 2021; 372eabg5298
        • Ball G.
        • Beare R.
        • Seal M.L.
        Charting shared developmental trajectories of cortical thickness and structural connectivity in childhood and adolescence.
        Hum Brain Mapp. 2019; 40: 4630-4644