Abstract
The ongoing NewTics study examines children who have had tics for less than 9 months (NT group), a population on which little research exists. Here, we further investigate relationships between subcortical shape and tic symptom outcomes. 187 children were assessed at baseline and a 12-month follow-up: 88 with NT, 60 tic-free healthy controls (HC), and 39 with chronic tic disorder or Tourette syndrome (TS), using T1-weighted MRI and total tic scores (TTS) from the Yale Global Tic Severity Scale to evaluate symptom change. Subcortical surface maps were generated using FreeSurfer-initialized large deformation diffeomorphic metric mapping, and linear regression models were constructed to correlate structural shapes with TTS while accounting for co-variates, with relationships mapped onto structure surfaces. We found the NT group to have a larger right hippocampus compared to healthy controls. Subsequent analyses including clinical symptoms revealed a significant correlation between the future worsening of tic symptoms and a larger pallidum and thalamus at baseline. Surface maps illustrate distinct patterns of inward deformation (localized volume loss) in the putamen and outward deformation (localized volume gain) in the thalamus for the NT group compared to healthy controls. We also found distinct patterns of outward deformation in almost all studied structures when comparing the TS group to healthy controls. In the significant vertices of this comparison, the caudate further exhibited an overall trend of outward deformation compared to the average template in the TS group compared to both the NT group and controls. When comparing the NT and TS groups, the NT group showed consistent outward deformation in the caudate, accumbens, putamen, and thalamus. Since the NT group has had tics only for a few months, we can rule out the possibility that these subcortical volume differences are caused by living with tics for years; they are more likely related to the cause of tics. These observations constitute some of the first prognostic biomarkers for tic disorders and suggest localized circuitry that may be associated with outcome of tic disorders.
Introduction
Chronic tic disorders (CTD) were once thought to be rare but are now known to be relatively common [1]. Tics are sudden, repetitive, nonrhythmic movements or vocalizations such as blinks or grunting [2]. Transient tics affect at least 20% of children, though only about 3% of all children have tics for a full year, the requirement to diagnose a chronic tic disorder or Tourette syndrome (TS) [3, 4]. When tics are present but have not yet persisted for one year since onset, Provisional Tic Disorder is diagnosed [1]. Efforts to identify biomarkers for tics and study the pathophysiology of tic disorders have recently been increasing, although our understanding is still limited [5, 6].
Studies exploring differences in subcortical structure and function have often been contradictory, with some finding no significant differences in basal ganglia volumes or shape between children with TS and matched control children [6]. Two groups found increased putamen volume in TS compared to HC, but a larger study found decreased volume [7-9]. A large study of basal ganglia volume in vivo found the caudate to be 4.9% smaller in the TS group [10]. Smaller studies found lower caudate volume [11, 12], and another large study identified no change [13].
Previous such studies have also primarily focused only on TS (diagnosed chronic tics) and control samples. Therefore, we cannot determine whether the identified differences reflect an underlying cause of tics or secondary changes due to prolonged tic presence. Examining participants at the onset of tic symptoms will more likely lead to identifying biomarkers related to the primary cause of tics.
Thus, the NewTics study examined children who had experienced tics for less than 9 months (median 4 months; new tics, or NT) [14]. Little research exists on this population, and there are even fewer results on predictive outcome analyses, which have also been contradictory [4]. The NewTics study aimed to see whether features including subcortical structures measured shortly after tic onset can predict symptom severity 12 months after tic onset [14]. A previous volumetric MRI analysis using data from 65 children with NT found that striatal volumes did not predict outcome, but a larger hippocampus at baseline predicted worse severity at follow-up [15].
In the present study, we further investigated neurobiological characteristics and predictors of tic disorders by examining relationships of the shape of these subcortical structures with tic symptom outcomes (NT) and with diagnosis (NT, TS, and tic-free controls). Three-dimensional surface analysis can detect subtler or more localized volumetric changes that are not revealed in whole-structure, scalar volumetric analysis [16]. Using whole-structure volume estimates alone may yield false negative findings by overlooking local deformities in shape. Previously, diffeomorphic mapping of structural magnetic resonance imaging (MRI) has successfully mapped pathological biomarker patterns onto surface-based representations of anatomical structures [17].
We predicted that baseline volumes would differ across NT, TS, and control groups, and that subcortical shape would demonstrate distinct patterns of shape deformation in tic disorders. We further predicted that we would find distinct regions of shape deformation in subcortical structures at baseline that predict clinical outcome in terms of tic severity changes after one year. We also tested whether shape deformation analyses would confirm the previous finding that hippocampal volume predicted symptom severity outcome, using a 3D method in an expanded sample.
Materials and Methods
Subjects and data collection
Subjects
The sample consisted of 187 children across 3 groups: children examined within 9 months after tic onset (median 3.5 months; new tic, or NT), tic-free children with no parental or sibling history of tics (healthy controls, HC), and children who at the time of screening already have TS/CTD (TS group) [14]. NT subjects returned for clinical evaluation at the one-year anniversary of the best estimate of tic onset. For TS and HC subjects, follow-up visits occurred as near as possible to the same time after screening as it did for their matched NT subject, based on age, sex, and handedness.
Enrollment criteria
Subjects were limited to ages 5-15 at enrollment to reduce sample variance (changed to ages 5-10 after enrolling only 2 participants over age 10 in the first year or so). NT children had current tic symptoms, with the first tic starting less than 9 months before enrollment. Exclusion criteria are detailed elsewhere [14]. The TS group included children who met DSM-5 criteria for TS/CTD at enrollment. The exclusion criteria used reflect those of the NT group. Control children had no current or past tics.
Clinical data collection
For participants exhibiting tic symptoms, a best-estimate date of onset was taken as described in [18]. Total tic scores from the Yale Global Tic Severity Scale (YGTSS), reflecting current tic severity, were determined during a neurological and psychiatric examination performed by Dr. Black for NT and TS subjects. Additional assessments were done at the time of screening [14]. The YGTSS is a tool used to quantify the severity of tic symptoms in children and adolescents. The Total Tic Severity Score (TTS) comprises half of the YGTSS score and has a range of 0-50; a higher score indicates more severe tic symptoms [19].
Imaging data collection
All subjects at entry were assessed through 1 mm3 T1-weighted magnetization prepared rapid gradient echo (MPRAGE) images and T2-weighted images using three different scanners across the data acquisition period. Details of scan parameters are given in Kim et al (2020) [15]. About half of the scans (the newer ones) were acquired with a prospective motion correction sequence (vNavs) that substantially reduced head motion artifact [20].
Additional measures were taken prior to scanning to reduce motion effects on images, including the use of a mock scanner, as well as an informational video and a game for children to practice holding still during scanning, adapted from a previous study in children [21].
Data processing / Image processing
Data processing began with the use of the FreeSurfer (version 6.0.0) software’s probabilistic voxel-based classification, which provided initial subcortical segmentations [22]. Surfaces of the hippocampus, amygdala, basal ganglia (caudate, putamen, pallidum, and nucleus accumbens), and thalamus were automatically generated for each participant using multi-atlas FreeSurfer-initialized large deformation diffeomorphic metric mapping (FS+LDDMM), which utilizes automated brain segmentations based on multiple template images and allows for image alignment and intensity normalization to produce smooth transformations for each region of interest [22, 23]. Combining maps from multiple atlases that best match an individual’s scan features has shown improved segmentation accuracy and reduced biases [24]. An experienced rater (the first author) inspected the final surfaces and made minor manual edits on the initial segmentation on poor surface maps, affecting maps of 6 subjects. The edited segmentations were then re-processed via LDDMM to yield accurate maps to be included in subsequent surface analyses.
Surfaces were then scaled with a scale factor calculated for each subject using the population TIV (estimated total intracranial volume) over the individual’s TIV, as well as the voxel resolution of their scan. To calculate the (one-dimensional) scale factor, we used the formula: (Population TIV / Individual TIV)? * Voxel Resolution.
Local shape variation for each participant was calculated from the population average of all participants by quantifying the vertex-to-vertex perpendicular change between surfaces, which were assigned a positive (outward variation from the population average) or negative (inward variation from the population average) value [22]. Subcortical volumes for each participant were determined using the volume enclosed within the surfaces. FreeSurfer also reported estimated total intracranial volume (TIV) to be used as a covariate in surface analyses, as TIV estimated by FreeSurfer segmentation has been shown to influence subcortical volumes [23].
Statistical analysis
TTS: Paired T-tests were conducted to determine significance of TTS changes from baseline to 12-months within each group, as well as between the NT and TS groups.
TIV: We first conducted a one-way ANOVA on TIV (dependent variable) to determine group (independent variable) effects with post hoc Tukey HSD tests for pairwise differences, with and without age and sex as covariates. The tests showed significant group effects (see results); therefore, TIV was used as an additional covariate in all subsequent analyses.
Scanner type: We first conducted a chi-square test that determined that neither scanner types nor race differed between groups. Additionally, we conducted one-way ANOVA tests to determine the effects of scanner type on structural volumes. We found that the subcortical structures did not differ in volume based on the scanner type, and subsequent ANCOVA and post hoc tests found that group differences found in TIV between NT and healthy controls still occur after controlling for the scanner type. Thus, we did not use scanner as a covariate in subsequent analyses.
Race: The tests done for scanner type (see above) were applied to race as well, and similar results were found. We thus also did not use race as a covariate in subsequent analyses.
Subcortical volumes: Using the surface volumes for each subcortical structure, we first found estimated marginal means and standard errors for each group.
ANCOVA and post hoc tests on subcortical volumes (dependent variables) were then conducted to compare group (independent variable) differences of baseline structural volumes, using age and sex as covariates. All structures were examined by combining left and right subcortical structures, since we did not have a lateralized hypothesis. If groups were found to differ significantly for a specific structure, further analyses were conducted to examine left and right structures separately. Analyses done in R used version 4.0.5, with psych and ggplot2 packages [25, 26].
We then performed a partial correlation analysis using baseline structural volume to predict 12-month TTS in the NT group while controlling for screen TTS, age, and sex.
Shape: For group comparisons, ANCOVAs were conducted to compare pairwise group differences of baseline surface shape (i.e., NT vs. HC, TS vs. HC, and TS vs. NT). All models included covariates for age and sex. Surface comparisons were conducted using SurfStat implemented in MATLAB [27]. We applied random field theory (RFT) to identify significant clusters of vertices at the family-wise error rate (FWER) of p<0.05 within each subcortical structure, to account for the multiple comparisons inherent in surface maps [28]. Group differences were visualized as a color map displayed on the overall average surface [22].
We then extracted the significant vertices found in the TS-HC surface comparison in applicable structures and found the mean values of the deformation from the average template of these vertices for each subject. Additionally, significant vertices from the NT-HC surface comparisons were studied to run partial correlation analyses on each subject’s mean deformation value and their 12-month clinical score (TTS), while controlling for screen TTS, age, and sex. Subsequently, we repeated this process using significant vertices we extracted from the NT-TS surface comparisons.
Results
Subject demographics and clinical features
In total, we analyzed 187 participants: 60 HC (44M / 16F), 88 NT (63M / 25F), and 39 TS (28M / 11F) (Table 1). Sex, mean age (all near 8 years old), scanner type used, and race did not differ significantly between groups.
Baseline total intracranial volume
The one-way ANOVA on TIV with age and sex as covariates revealed significant group differences [F(2,182)=3.316, p=0.004; see Fig. 1]. Post hoc Tukey HSD tests revealed that the NT group had smaller TIV than controls at baseline (NT=1484 ± 154 cm3 and HC= 1596 ± 141 cm3; p=0.03). Group differences in TIV between NT and healthy controls remained significant after controlling for both scanner type and race independently. Additionally, the TIV differed significantly across scanner types (p=0.04) and race (p=0.04). Since TIV differed significantly among groups, we henceforth used baseline TIV as a covariate in subsequent analyses.
Baseline group comparisons of subcortical structural volumes
Table 2 summarizes the estimated marginal means and standard errors of the subcortical volume for each group accounting for age and sex, as well as statistical comparisons for each group. Here, we describe the results for each structure in detail.
Hippocampus
There were no statistically significant differences (though very near significance) among groups when examining the entire hippocampus. However, differences between the NT group and healthy controls were found when isolating the right hippocampus. Patients in NT had on average an 8.5% larger hippo-campus compared to those in the HC group.
Amygdala, caudate, accumbens, putamen, pallidum, and thalamus
No group differences were found among HC, NT, and TS groups.
Prediction of 12-month TTS from baseline structural volume
Longitudinal TTS analyses included 80 NT subjects. Over the course of the 12 months between the baseline and second visits, the average TTS decreased significantly from 16.96 ± 5.64 to 14.16 ± 6.94 (t = 3.82, df = 79, p-value = 0.0003).
Table 3 summarizes the partial correlation values between structural baseline volumes and TTS changes, while controlling for screen TTS, age, and sex. Here, we describe the results for each structure in detail.
Hippocampus, amygdala, caudate, accumbens, and putamen
Scaled baseline structural volumes did not predict TTS changes in the NT group.
Pallidum
We found a significant positive correlation between baseline pallidal volume and TTS change from baseline to 12-month measurements (r=0.24, p=0.04). In other words, a larger pallidum at baseline was significantly correlated with the worsening of tic symptoms.
Thalamus
We found a significant positive correlation between baseline thalamic volume and TTS change from baseline to 12-month measurements (r=0.23, p=0.05). In other words, a larger thalamus at baseline was significantly correlated with the worsening of tic symptoms.
Baseline group comparisons of subcortical structural shape
When comparing the NT group to healthy controls, two structures were found to have significant regions of deformation as described below.
Putamen. (Figure 2 Panel A). Inferiorly, a small region of inward deformation is seen in the right medial putamen.
Thalamus. (Figure 2 Panel B). We found a region of outward deformation in the medial parts of the left thalamus inferiorly.
Hippocampus, amygdala, caudate, accumbens, and pallidum. No group differences were found among HC, NT, and TS groups.
When comparing the NT group to TS children, six structures were found to have significant regions of deformation as described below.
Hippocampus. (Figure 3 Panel A). Superiorly, regions of outward deformation are concentrated towards the lateral-posterior parts of both the left and right hippo-campus.
Amygdala. No group differences were found among HC, NT, and TS groups.
Caudate. (Figure 3 Panel B). Superiorly, large regions of outward deformation can be seen along the lateral edge of both the left and right caudate. Inferiorly, regions of outward deformation are seen in the medial anterior parts of the left and right caudate. Additionally, a small region of inward deformation exists along the lateral edge of the left caudate inferiorly.
Accumbens. (Figure 3 Panel C).
Putamen. (Figure 3 Panel D). Inferiorly, a small region of inward deformation is seen in the right medial putamen.
Thalamus. (Figure 3 Panel E). We found a region of outward deformation in the medial parts of the left thalamus inferiorly.
Pallidum. (Figure 3 Panel F). We found regions of outward deformation in both the left and right pallidum inferiorly and superiorly.
When comparing the TS group to healthy controls, four structures were found to have significant regions of outward deformation. In other words, areas of outward deformation described below illustrate localized volume gain in the specified regions in TS.
Hippocampus, amygdala, and pallidum. No group differences were found among HC, NT, and TS groups.
Caudate. (Figure 3 Panel A). Inferiorly, two small regions of outward deformation is seen in the both the left and right medial caudate.
Accumbens. (Figure 3 Panel B). Superiorly and inferiorly, regions of outward deformation can be seen in the right accumbens.
Putamen. (Figure 3 Panel C). Inferiorly, small regions of outward deformation is seen in the left and right medial anterior putamen.
Thalamus. (Figure 3 Panel D). We found regions of outward deformation towards the medial ends of the left and right thalamus.
Table 4 summarizes the mean values of the deformation from the average template for each group within the significant vertices found in the TS-HC surface comparison, as well as statistical comparisons for each group. Here, we describe the results for each structure in detail.
Caudate. We found an overall trend of outward deformation compared to the average template in the TS group compared to both the HC and NT groups.
Accumbens. There were no statistically significant differences between groups. Putamen. There were no statistically significant differences between groups.
Thalamus. There were no statistically significant differences between groups.
All significant vertices. There were no statistically significant differences between groups.
Prediction of 12-month TTS from baseline structural shape
When comparing the NT group with healthy controls, we found that 18 out of the 13638 vertices on the putamen surface had a significant deformation value. On the thalamus, 45 out of 9580 vertices were significant.
Table 5 summarizes the correlation and p-values between each NT subject’s mean deformation value in the significant vertices of the NT-HC comparison and their 12-month TTS, while controlling for screen TTS, age, and sex. We did not find a significant correlation between mean deformation values and 12-month TTS.
When comparing the NT group with the TS group, amounts of significant vertices on structure surfaces were found as follows:183/13222 significant vertices on the hippo-campus surface; 3734/24744 significant vertices on the caudate surface; 2279/5804 significant vertices on the nucleus accumbens surface; 1030/5398 significant vertices on the pallidum surface; 3345/13638 significant vertices on the putamen surface; 1956/9580 significant vertices on the thalamus surface.
NT-TS Partial Correlations between Deformation and 12-month TTS
Table 6 summarizes the correlation and p-values between each NT subject’s mean deformation value in the significant vertices of the NT-TS comparison and their 12-month TTS, while controlling for screen TTS, age, and sex. We did not find a significant correlation between mean deformation values and 12-month TTS.
Discussion
In this study, we focus on the first year of tic development (NT group) and identified several baseline subcortical volume and shape characteristics related to tic symptoms, including some that predict clinical tic outcome 3 to 12 (mean 7.9) months later.
Baseline volume analyses showed that the right hippocampus was on average 8.5% larger in NT children compared to healthy controls. We also found a significant correlation between the worsening of tic symptoms and both a larger pallidum and thalamus at baseline.
Surface analyses demonstrated distinct patterns of subcortical surface deformation in several structures across all group comparisons. When comparing the NT group to healthy controls, we found that the putamen exhibited inward deformation (i.e., localized volume loss), while the thalamus showed outward deformation (localized volume gain). Upon comparing the TS group to healthy controls, we found that structures consistently showed patterns of outward deformation in the TS group, found significantly in almost all of the studied structures, including the hippocampus, caudate, accumbens, putamen, thalamus, and pallidum. Similarly, we also found that this pattern was present when comparing the NT group to the TS group, with the NT group showing consistent patterns of outward deformation in the caudate, accumbens, putamen, and thalamus.
In the significant vertices found in the TS-HC surface comparison of the caudate, we found an overall trend of outward deformation compared to the average template in the TS group compared to both the HC and NT groups. When examining the relationship between each NT subject’s mean deformation value and their 12-month TTS, we did not find a significant correlation in neither the NT-HC significant vertices nor the NT-TS significant vertices.
Functional implications
Hippocampus
Children in the NT group had larger right hippocampal volume compared to healthy controls. Further, the NT group also exhibited localized volume increases compared to the TS group. The hippocampus plays a role in memory consolidation in both the cognitive and motor domains [29]. Specifically, the hippocampus supports visuospatial ability and memory. In a study on patients with amnestic mild cognitive impairment, a larger right hippocampus was significantly associated with higher visuospatial memory scores on the Rey’s complex figure test, which measures visuospatial ability and memory [30]. Furthermore, in a study on the effects of traumatic brain injury in the hippocampus, damage to the right hippocampus correlated with lower memory scores on the Rey’s complex figure test [31]. Our conjecture is that damage here impairs visuospatial ability and memory, and abnormal enlargement may conversely lead to abnormal preservation of visuospatial memory. There is evidence that in children with TS, those with more persistent motor memory show more severe tics, as they were found to take longer to unlearn a previously learned motor pattern [32]. As we saw that the NT group had a larger right hippocampus at baseline compared to healthy controls, perhaps this region may be implicated in tic development and/or persistence.
Basal Ganglia (accumbens, pallidum, putamen, and caudate)
Through structural volume analyses, we saw that the worsening of tic symptoms was associated with a larger baseline pallidum. When analyzing shape differences, we found significant patterns of localized volume increases in NT children compared to TS children in all four structures. However, the putamen in the NT group exhibited localized volume loss when compared to healthy controls, but this deformation was concentrated towards the posterior end. Additionally, the TS group exhibited localized volume increases in the caudate, accumbens, and putamen when compared to healthy controls. In the significant vertices found in the TS-HC surface comparison of the caudate, we found an overall trend of outward deformation compared to the average template in the TS group compared to both the HC and NT groups.
The basal ganglia are subcortical nuclei that, along with their associated connections, have been a large focus of TS/CTD research. They are involved in motor control in other movement disorders such as Parkinson’s and Huntington’s disease. Cortico-striatal-thalamo-cortical (CSTC) circuits are involved in inhibitory control and habit formation, both to which are affected in TS [33, 34]. This understanding may relate to the shape and volume differences we found in basal ganglia structures and the pallidum’s association with observed changes in tic symptoms.
Thalamus
In the thalamus, there were no significant volumetric differences between groups. When analyzing surfaces, however, we found that the thalamus displayed localized volume increases in NT children when compared to both TS children and healthy controls. Previous volume and surface morphology studies have found greater grey matter volume and outward deformation in the thalamus of TS children compared to healthy controls, which is again seen here in our surface comparison between TS and healthy children [13, 35]. Further, we found a significant correlation between a larger thalamus at baseline and the worsening of tic symptoms. The position of the thalamus in CSTC circuits, mediating output from motor-related basal ganglia and motor areas of the cerebral cortex, supports its involvement in motor control [36]. Perhaps this region may be implicated in tic development and/or persistence.
We found the NT group to have a larger right hippocampus at baseline compared to healthy controls. In a prior analysis of a smaller subset of the present sample using different methods, a larger hippocampus at baseline predicted worse severity at follow-up [15]. Since the NT group has had tics for only a few months, the subcortical volume differences in this group are especially meaningful as they are more likely to be related to the cause of tics rather than a result of living with chronic tics.
Our results are also consistent with the results of a meta-analysis on task-based fMRI studies in patients with TS [37]. In the thalamus, the meta-analysis showed a clear overlap of all the conditions involving various aspects of voluntary motor execution, response inhibition, and tic generation. The pallidum and thalamus both showed more consistent activation in free-to-tic conditions. This result is consistent with our findings of larger thalamic and pallidal volume at baseline predicting the worsening of tic symptoms. The meta-analysis also identified a positive correlation between tic severity and BOLD activity in the thalamus and putamen [37]. Similarly, in the putamen and thalamus, we found that the NT group exhibited localized volume gain compared to the TS group, and the TS group showed localized volume gain compared to the healthy controls. It seems that these regions are implicated in both tic severity and persistence.
The primary limitation of this project is the cross-sectional nature of the results, as we are collecting more 12-month scan and TTS data. Future studies shall include longitudinal and within-subject analyses using both baseline and 12-month scans, in combination with clinical scores. Nevertheless, our study is, to our knowledge, the first imaging study of Provisional Tic Disorder.
These new findings have potential clinical relevance, and continuing to investigate neuroanatomical characteristics in TS/CTD may further provide insight into prognostic biomarkers. In many children with Provisional Tic Disorder, tics improve within the first year, often to the point of clinical insignificance. Understanding the mechanisms related to these outcomes may provide clinical insight into the pathophysiological traits in tic disorders and thus potentially lead to improved treatment approaches. In a previous memory study, multiple-session stimulation increased functional connectivity among cortical-hippocampal network regions and simultaneously improved associative memory performance [38]. With a deeper understanding of neural networks related to tic symptoms, we may be able to similarly provide targets that could be manipulated to prevent worsening of tic symptoms.
Data Availability
All data from this project is being archived at the NIMH Data Archive (https://nda.nih.gov, study ID C2692).
Author Contributions
Conceptualization, K.J.B., B.L.S. and L.W.; methodology, L.W.; formal analysis, T.C., L.W., J.D.; investigation, K.J.B., D.J.G., B.L.S., and S.K.; data curation, K.J.B. and S.K.; writing—original draft preparation, T.C.; writing—review and editing, all authors; visualization, T.C.; project administration, K.J.B.; funding acquisition, K.J.B., D.J.G. and B.L.S. All authors have read and agreed to the published version of the manuscript.
Funding
Research reported in this publication was supported by the National Institutes of Health: National Institute of Mental Health under award number K24MH087913 to K.J.B.; K01MH104592 to D.J.G.; R01MH104030 to K.J.B. and B.L.S.; R01MH118217 to D.J.G.; National institute of Neurological Disorders and Stroke R21NS091635 to B.L.S. and K.J.B.; the Washington University Institute of Clinical and Translational Sciences grants UL1RR024992 and UL1TR000448; the Eunice Kennedy Shriver National Institute of Child Health and Human Development of the National Institutes of Health under Award Number U54HD087011 to the Intellectual and Developmental Disabilities Research Center at Washington University, and K23DC006638 to J.E.C. Lieu; the Mallinckrodt Institute of Radiology MIR-IDDRC Pilot Study Fund; the McDonnell Center for Systems Neuro-science; the National Institute of Biomedical Imaging and Bioengineering grant R01 EB020062 to L.W.; the National Science Foundation grants 1734853 and 1636893 to L.W. The studies presented in this work were carried out in part in the East Building MR Facility of the Washington University Medical Center. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health or of the MIR.
Institutional Review Board Statement
Protocols 201109157 and 201707059 were approved by the Washington University Human Research Protection Office (IRB).
Informed Consent Statement
Each child assented, and a parent (guardian) gave informed consent. Data shared from other projects were shared after appropriate human subjects review and consent.
Data Availability Statement
All data from this project is being archived at the NIMH Data Archive (https://nda.nih.gov, study ID C2692).
Conflicts of Interest
The authors declare no conflict of interest.
Footnotes
tiffanieche2019{at}u.northwestern.edu; Lei.Wang{at}osumc.edu
soyoung{at}brainkim.com
deannagreene{at}ucsd.edu
ashley.heywood{at}northwestern.edu
jmding{at}wustl.edu
schlaggar{at}kennedykrieger.org
kevin{at}wustl.edu
Additional subject data included; Results and figures updated