Study Population
Longitudinal data were collected as a part of the ongoing longitudinal ABCD Study, which enrolled 11,876 children at ages 9–10 years across 21 study sites. Study enrollment criteria included age (≤ 10-years-old at initial visit) and English language proficiency. Exclusion criteria included major medical or neurological conditions, history of traumatic brain injury, diagnosis of schizophrenia, moderate/severe autism spectrum disorder, intellectual disability, alcohol/substance use disorder, premature birth (gestational age < 28 weeks), low birthweight (< 1200 g), and contraindications to MRI scanning.68 All study procedures were approved by the centralized institutional review board at the University of California San Diego; each study site also obtained approval from their institutional review boards. Participants provided written assent and legal guardians provided written consent.
We used a subset of data from the ABCD Study, including magnetic resonance imaging (MRI) from the baseline and/or year-2 follow-up study visits and measures of participants’ age, sex at birth, and sociodemographics held constant from the baseline assessment. Only high-quality imaging scans completed before March 1, 2020 were included to remove potential confounding effects of stress inherent to the COVID-19 pandemic. We filtered for valid air pollution estimates (see quality control details below), and randomly selected one subject per family to reduce the number of hierarchical levels, uneven by study design (i.e., the number of both siblings and twins vary by site). Our final sample included 8,182 subjects across 21 study sites. Of these, 3,679 (45%) had two time points of high-quality DWI data, while 4503 (55%) had one DWI time point, either from the baseline or 2-year follow-up visit (see details below; Table 1). All data used here were obtained from ABCD’s 4.0 data release (http://dx.doi.org/10.15154/1523041).
Ambient Air Pollution Estimates
Annual ambient air pollution concentration for PM2.5, NO2, and O3 were assigned to primary residential addresses of each child as previously described.69 Briefly, daily estimates were derived at a 1-km2 resolution using hybrid spatiotemporal models, utilizing satellite-based aerosol optical depth models, land-use regression, and chemical transport models,69–71 and averaged over the 2016 calendar year, corresponding with enrollment for the baseline assessment. One-year annual average concentrations were then assigned to primary residential address at the baseline assessment when children were aged 9–10 years. PM2.5 was positively correlated with NO2 (r = 0.21, p = 3.44e-81) and negatively correlated with O3 (r = -0.19, p = 1.78e-64); there was no correlation between NO2 and O3 (r = -0.02, p = 0.12).
Diffusion Weighted Imaging (DWI): Acquisition, Processing, and Quality Control
A harmonized neuroimaging protocol was utilized across sites, given the differences in scanner manufacturer (3T Siemens, Phillips, or GE). The multi-shell DWI acquisition included a voxel size of 1.7 mm isotropic, implemented multiband EPI72,73 with slice acceleration factor 3, and included a fieldmap scan for B0 distortion correction. Seven b = 0 frames and 96 total diffusion directions at 4 b-values (6 with b = 500 s/mm2, 15 with b = 1000 s/mm2, 15 with b = 2000 s/mm2, and 60 with b = 3000 s/mm2) were collected.74 All images underwent distortion, bias field, and motion correction, and manual and automated quality control.74 After preprocessing, white matter tracts were identified using the probabilistic atlas AtlasTrack.75 Only images without clinically significant incidental findings (mrif_score = 1 or 2) that passed all ABCD quality-control parameters (imgincl_dmri_include = 1) were included in analysis.
Restriction Spectrum Imaging (RSI)
Restriction spectrum imaging (RSI) utilizes all 96 directions in ABCD’s multi-shell acquisition protocol.76 RSI provides detailed information regarding both the extracellular and intracellular compartments of tissue within the brain.20 RSI model outputs are normalized measures, unitless on a scale of 0–1. We focused on restricted (intracellular) normalized isotropic signal fraction (RNI) and restricted normalized directional signal fraction (RND) white matter fiber tract ROIs created with AtlasTrack.75 We explored all tracts excluding summary tracts (14 in the left hemisphere, 14 in the right hemisphere, and 3 spanning both hemispheres), including the right and left fornix, cingulate cingulum, parahippocampal cingulum, corticospinal tract, anterior thalamic radiations, uncinate fasciculi, inferior longitudinal fasciculi, inferior fronto-occipital fasciculi, temporal and parietal superior longitudinal fasciculi, frontal and parietal superior corticostriate, striatal to inferior frontal cortex, and inferior frontal to superior frontal cortex, as well as the forceps major, minor, and corpus callosum.
Predictors
Predictors were chosen using a directed acyclic graph and included demographic and socioeconomic variables: race/ethnicity (race_ethnicity variable with the following categories: White, Black, Hispanic, Asian, or Other), annual household income (USD; >100K, 50-100K, < 50K, or Don’t Know/Refuse to Answer), and highest household education (Post-Graduate, Bachelor, Some College, High School Diploma/GED, or < High School Diploma). Pollution levels are higher in minority communities and those from disadvantaged social status backgrounds due to structural racism and class bias increasing the likely proximity of these communities to major sources of pollution in the U.S. 77,78 Census Tract Urban Classification (Rural, Urban Clusters, or Urbanized Area) was included as air pollution levels vary by degree of urbanicity. We also included the subject-specific precision variable, handedness (right, left, or mixed) and MRI-related precision variables such as scanner manufacturer (Siemens, Philips, GE) to account for differences in both scanner hardware and software, average frame displacement (mm) to account for head motion, and tract volume. Lastly, due to the potential acute differences in seasonality of pollutant concentrations at the time of each visit, we included the meteorological season of the MRI scan date as an additional time-varying variable.
Statistical Analyses
We used hierarchical linear mixed-effect models, as implemented in lme4::lmer() in R statistical software (Version 4.1.2.).79 We first tested a developmental model, examining the main effects of age and sex, as well as an age-by-sex interaction term. Next, we tested the longitudinal change of RSI outcomes in the context of specific pollutants (PM2.5, O3, and NO2), opting for sex-stratified models over including pollutant-by-age-by-sex or pollutant-by-sex interaction terms; the outcomes and a number of predictors in the model have demonstrated sex-specific effects and the inclusion of an interaction term would likely introduce bias.24 In each model, subjects nested within ABCD sites were modeled as random effects, to account for the multi-level data structure. Age was centered on the lowest age within our sample (107 months), resulting in a scaled age score of 0 equivalent to 8.9 years. We accounted for non-linearity of age by utilizing a piecewise linear spline model, placing a knot at median age of 127 months. This two-piece linear spline model was parameterized to include an overall effect of age and an age-deviation (ageD) term. To assess how air pollution exposures affect age-related maturation in WM microstructure, we have included an interaction of pollutants with the age-specific spline terms. Specifically, each model included the pollutant of interest (PM2.5, O3, or NO2), age, ageD, interactions of the specific pollutant with age and ageD, and all predictors discussed above. To account for co-exposure of the three criteria pollutants, we additionally controlled for the other two pollutants not included in the age-by-pollutant interaction term of interest. For example, we included NO2 and O3 as covariates in the model testing the effects of PM2.5, age, and age-by-PM2.5 (plus ageD, ageD-by-PM2.5, and previously mentioned predictors) on RSI outcomes. Parameters of interest included the fixed effects of the pollutant on attainted WM microstructure at age 9 (i.e., scaled age score of 0), age (i.e., time), and the age-by-pollutant interaction term to investigate how WM maturation may be altered by air pollution exposure. To account for multiple comparisons due to modeling RNI and RND for all available white matter tracts, we performed a false discovery rate (FDR) correction for 62 tests (31 total white matter tracts across both sexes).