Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE
PreviousNext
Research Articles, Development/Plasticity/Repair

Graded Variation in T1w/T2w Ratio during Adolescence: Measurement, Caveats, and Implications for Development of Cortical Myelin

Graham L. Baum, John C. Flournoy, Matthew F. Glasser, Michael P. Harms, Patrick Mair, Ashley F. P. Sanders, Deanna M. Barch, Randy L. Buckner, Susan Bookheimer, Mirella Dapretto, Stephen Smith, Kathleen M. Thomas, Essa Yacoub, David C. Van Essen and Leah H. Somerville
Journal of Neuroscience 20 July 2022, 42 (29) 5681-5694; https://doi.org/10.1523/JNEUROSCI.2380-21.2022
Graham L. Baum
1Department of Psychology and Center for Brain Science, Harvard University, Cambridge, MA, USA, 02138
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John C. Flournoy
1Department of Psychology and Center for Brain Science, Harvard University, Cambridge, MA, USA, 02138
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Matthew F. Glasser
2Department of Radiology, Washington University School of Medicine, St. Louis, MO, USA, 63110
3Department of Neuroscience, Washington University School of Medicine, St. Louis, MO, USA, 63110
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael P. Harms
4Department of Psychiatry, Washington University School of Medicine, St. Louis, MO, USA, 63110
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Patrick Mair
1Department of Psychology and Center for Brain Science, Harvard University, Cambridge, MA, USA, 02138
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ashley F. P. Sanders
4Department of Psychiatry, Washington University School of Medicine, St. Louis, MO, USA, 63110
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ashley F. P. Sanders
Deanna M. Barch
4Department of Psychiatry, Washington University School of Medicine, St. Louis, MO, USA, 63110
5Department of Psychological and Brain Sciences, Washington University, St. Louis, MO, USA, MO 63130
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Deanna M. Barch
Randy L. Buckner
1Department of Psychology and Center for Brain Science, Harvard University, Cambridge, MA, USA, 02138
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Susan Bookheimer
6Department of Psychiatry and Biobehavioral Sciences, David Geffen School of Medicine at UCLA, University of California, Los Angeles, CA, USA, 90095
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mirella Dapretto
6Department of Psychiatry and Biobehavioral Sciences, David Geffen School of Medicine at UCLA, University of California, Los Angeles, CA, USA, 90095
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stephen Smith
7Wellcome Centre for Integrative Neuroimaging, University of Oxford, Oxford, United Kingdom, OX3 9DU
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kathleen M. Thomas
8Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, MN, USA, 55455
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Essa Yacoub
8Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, MN, USA, 55455
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David C. Van Essen
3Department of Neuroscience, Washington University School of Medicine, St. Louis, MO, USA, 63110
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for David C. Van Essen
Leah H. Somerville
1Department of Psychology and Center for Brain Science, Harvard University, Cambridge, MA, USA, 02138
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Adolescence is characterized by the maturation of cortical microstructure and connectivity supporting complex cognition and behavior. Axonal myelination influences brain connectivity during development by enhancing neural signaling speed and inhibiting plasticity. However, the maturational timing of cortical myelination during human adolescence remains poorly understood. Here, we take advantage of recent advances in high-resolution cortical T1w/T2w mapping methods, including principled correction of B1+ transmit field effects, using data from the Human Connectome Project in Development (HCP-D; N = 628, ages 8–21). We characterize microstructural changes relevant to myelination by estimating age-related differences in T1w/T2w throughout the cerebral neocortex from childhood to early adulthood. We apply Bayesian spline models and clustering analysis to demonstrate graded variation in age-dependent cortical T1w/T2w differences that are correlated with the sensorimotor-association (S-A) axis of cortical organization reported by others. In sensorimotor areas, T1w/T2w ratio measures start at high levels at early ages, increase at a fast pace, and decelerate at later ages (18–21). In intermediate multimodal areas along the S-A axis, T1w/T2w starts at intermediate levels and increases linearly at an intermediate pace. In transmodal/paralimbic association areas, T1w/T2w starts at low levels and increases linearly at the slowest pace. These data provide evidence for graded variation of the T1w/T2w ratio along the S-A axis that may reflect cortical myelination changes during adolescence underlying the development of complex information processing and psychological functioning. We discuss the implications of these results as well as caveats in interpreting magnetic resonance imaging (MRI)-based estimates of myelination.

SIGNIFICANCE STATEMENT Myelin is a lipid membrane that is essential to healthy brain function. Myelin wraps axons to increase neural signaling speed, enabling complex neuronal functioning underlying learning and cognition. Here, we characterize the developmental timing of myelination across the cerebral cortex during adolescence using a noninvasive proxy measure, T1w/T2w mapping. Our results provide new evidence demonstrating graded variation across the cortex in the timing of T1w/T2w changes during adolescence, with rapid T1w/T2w increases in lower-order sensory areas and gradual T1w/T2w increases in higher-order association areas. This spatial pattern of microstructural brain development closely parallels the sensorimotor-to-association axis of cortical organization and plasticity during ontogeny.

  • adolescence
  • development
  • MRI
  • myelin
  • T1w/T2w

Introduction

Adolescence, the transition from childhood to adulthood, is characterized by the refinement and stabilization of neural circuits supporting the dynamic control of attention and behavior (Larsen and Luna, 2018). Axonal myelination in the cerebral cortex enhances signaling speed (McDougall et al., 2018) and regulates synapse formation during development to support the reliable instantiation of adaptive behavior (Mount and Monje, 2017). While the role of cortical myelination in regulating plasticity during sensitive periods of neurodevelopment has been demonstrated in animal models (Takesian and Hensch, 2013), there are few in vivo studies characterizing the maturational timing of cortical myelination during sensitive periods of human development.

The cerebral cortex follows an overarching organizational axis that is anchored on one end by sensorimotor areas underlying perception and movement, and at the opposing end by heteromodal and paralimbic association areas involved with cognitive control and socioemotional processing (Huntenburg et al., 2018). The archetypal sensorimotor-association (S-A) axis of cortical organization was derived in a recent study by computing a composite ranking of cortical maps that capture feature variation in cortical microstructure, functional connectivity, metabolism, gene expression, and evolutionary expansion of the cortical mantle (Sydnor et al., 2021). Further, the spatiotemporal pattern of cortical development may also progress along the S-A axis, such that the early maturation of sensorimotor systems scaffolds the development of association areas that integrate more complex cognitive processes. Both myelin-sensitive imaging and postmortem histology have provided evidence for continuous remodeling of intracortical myelin content into adulthood (R.A. Hill et al., 2018; Hughes et al., 2018), but it remains unclear whether cortical areas exhibit graded differences from childhood through early adulthood in the developmental timing of myelination according to their position along the S-A axis.

Advances in noninvasive magnetic resonance imaging (MRI) methods have facilitated in vivo investigations of myelin content and tissue changes across the human lifespan (Edwards et al., 2018). The ratio of T1-weighted and T2-weighted MRI images (T1w/T2w) provides an indirect estimate of cortical myelin that correlates with histologic myelin stains and quantitative MRI (e.g., magnetization transfer, R1, R2) indices of cortical myelin content (Glasser and Van Essen, 2011; Glasser et al., 2014, 2021; Shams et al., 2019). T1w/T2w mapping has been increasingly used in recent years because it is based on widely available MRI sequences that are routinely collected in brain imaging studies at high resolution, at clinical field strengths, and with short acquisition times. However, prior studies evaluating age-related differences in T1w/T2w (Grydeland et al., 2019; Norbom et al., 2020) have not accounted for residual radio frequency transmit field (B1+) biases, which impact T1w/T2w estimates (as it does quantitative relaxometry measures) in an age-dependent manner through their relationship to body size (Glasser et al., 2021; MacLennan et al., 2022). Because studies that neglected B1+ bias field effects may be prone to spurious results, it is critical to characterize age-related differences in T1w/T2w after mitigating these confounds.

Here, we apply high-resolution cortical T1w/T2w mapping methods from the Human Connectome Project in Development (HCP-D; N = 628, ages 8–21) to characterize the maturational timing of cortical myelination from childhood through early adulthood. We test whether cortical areas exhibit graded differences in the rate of T1w/T2w development according to their position along the S-A axis. We use a novel approach for mitigating B1+ confound effects in T1w/T2w mapping that enables us to detect unbiased age-related differences in T1w/T2w (Glasser et al., 2021). We characterize properties of age-related differences in T1w/T2w across the cerebral cortex, including estimates of the variance in T1w/T2w explained by age, rate of change, age of peak growth, and degree of nonlinearity in T1w/T2w development from 8 to 21 years old. These data provide evidence for graded variation in age-dependent T1w/T2w differences along the S-A axis that may reflect variation in the rate of cortical myelination during adolescence and undergird the development of complex information processing.

Materials and Methods

Participants

Neuroimaging datasets included 628 typically developing participants aged 8–21 years (53.5% female; Fig. 1A) who were part of the HCP-D. The HCP-D is a large cross-sectional and longitudinal study aiming to characterize brain connectivity development in a sample approximating the demographics of youth in the United States with respect to race, ethnicity, and socioeconomic status (Harms et al., 2018; Somerville et al., 2018; Elam et al., 2021). Participants in this sample were recruited across four sites: Harvard University, University of California-Los Angeles, University of Minnesota, and Washington University in St. Louis. Exclusion criteria for recruitment included (1) premature birth (<37 weeks gestation); (2) serious neurologic condition (e.g., stroke, cerebral palsy); (3) serious endocrine condition (e.g., precocious puberty, untreated growth hormone deficiency); (4) long-term use of immunosuppressants or steroids; (5) history of serious head injury; (6) hospitalization >2 d for certain physical or psychiatric conditions or substance use; (7) treatment >12 months for psychiatric conditions; (8) claustrophobia; or (9) pregnancy or other contraindications for MRI. Participants provided written informed consent and assent, and parents of participants under 18 years provided written informed consent for their child's participation. All procedures were approved by a central Institutional Review Board administered at Washington University in St. Louis (IRB #201603135).

This study focuses on cross-sectional data drawn from the publicly available HCP-D “Release 2.0” dataset in the National Institute of Mental Health Data Archive (NDA; N = 652). From this dataset, we excluded 21 participants under the age of 8 because of insufficient data available for 5–7 year olds at the time of this study. We excluded two additional participants because of face-masking errors during image reconstruction, and we excluded one additional participant with poor T1w/T2w data quality indicated by visible artifacts and substantial outliers in T1w/T2w maps, yielding our final sample of 628 youth.

Myelin measurement: caveats and assumptions

MRI signal contrasts in cortical gray matter reflect a composite of tissue properties (Fukunaga et al., 2010; Stüber et al., 2014). Changes in the MRI signal may reflect changes in iron, myelin, lipids, water, and local cell density (Fukunaga et al., 2010; Stüber et al., 2014; Edwards et al., 2018). Several noninvasive MRI methods have been developed in recent years to estimate cortical myelin content, including quantitative T1, quantitative T2*, quantitative magnetization transfer, T2-based myelin water fraction methods (Carey et al., 2018; Edwards et al., 2018), and T1w/T2w ratio image intensities. These approaches are correlated with one another and with histologic myelin stains in the cerebral cortex (Glasser and Van Essen, 2011; Glasser et al., 2014, 2021), but the relative sensitivity to the various sources of MRI tissue contrast differ across these methods. Importantly, differences in these measures over the course of development may reflect changes in multiple tissue properties (e.g., not just myelin) and the properties that are changing may differ from those that primarily drive the overall image contrast (e.g., the known and well validated differences across the cortex between heavily myelinated and lightly myelinated cortical areas). Thus, age-dependent differences in cortical myelin measures should be interpreted cautiously as microstructural development likely reflecting myelination, while acknowledging the complexity of cortical tissue changes during development and the limitations of indirect proxy measures.

Image acquisition

High-resolution T1w MRI images were acquired on a 3T Siemens Prisma with a 32 channel head coil using a 3D multiecho MPRAGE sequence [Mugler and Brookeman, 1990; van der Kouwe et al., 2008; 0.8-mm isotropic voxels, TR/TI = 2500/1000 ms, TE = 1.8/3.6/5.4/7.2 ms, flip angle = 8°, in-plane (iPAT) acceleration factor of 2, TA = 8:22, up to 30 reacquired TRs]. Structural T2w images were acquired at 0.8 mm isotropic using the variable-flip-angle turbo-spin-echo 3D SPACE sequence (Mugler et al., 2000; TR/TE = 3200/564 ms; same in-plane acceleration, TA = 6:35, up to 25 reacquired TRs).

Both “PreScan Normalized” (Siemens' approach for removing the receive coil intensity profile) and non-normalized reconstructions were generated at the scanner. The former were used for image quality review at the scanner, while the latter were used as the inputs for subsequent processing (consistent with the use of non-normalized reconstructions as the input for the processing of the HCP Young Adult acquisitions; Van Essen et al., 2013). Both versions (of the T1w) image were used for estimating the B1– receive field to correct for effects of subject motion between the T1w and T2w images. Only the first two echoes of the T1w image were used in processing because of artifacts in the later echoes that affected surface reconstructions and T1w/T2w maps (Elam et al., 2021).

Volumetric navigators (vNavs) were embedded in the T1w and T2w sequences for prospective motion correction and for selective reacquisition of the lines in k-space that were heavily corrupted by subject motion (Tisdall et al., 2012). Real-time motion correction can substantially reduce bias in brain morphometry analyses, where motion might induce measurable morphometric differences (Reuter et al., 2015; Tisdall et al., 2016). If the T1w or T2w structural scans were nonetheless deemed to be of poor quality at the time of acquisition, they were reacquired (typically immediately in the same imaging session, although sometimes in a different session). Only the single pair of T1w and T2w scans rated the highest in quality from a single session were used in subsequent processing. For further neuroimaging protocol description, see Harms et al. (2018). Additionally, 2-mm isotropic gradient echo (GRE) and spin echo (SE) images were acquired and used for computing the pseudo-transmit field described below (Glasser et al., 2021).

Image processing

The structural MRI data were analyzed using the HCP Pipelines (Glasser et al., 2013) version 4.0.0, instantiated into the QuNex container environment (qunex.yale.edu), and the data were released as a part of the Lifespan HCP Release 2.0 in the NDA (https://nda.nih.gov/). Briefly, the T1w and T2w volumes were processed through the PreFreeSurfer pipeline, which included gradient nonlinearity distortion correction, initial brain-extraction, and rigid registration into an anterior/posterior-commissure aligned “native” space, registration of the T2w volume to the T1w volume using boundary-based registration (BBR; Greve and Fischl, 2009), correction for the receiver coil bias field based on the smoothed square root of the product of the T1w and T2w images, and registration of the structural images to MNI space. Next, the FreeSurfer pipeline (v6.0.0; Dale et al., 1999; Fischl, 2012) optimized for use with high-spatial resolution (>1 mm isotropic) structural images was used for computing the “white” and “pial” surfaces, including use of the T2w volume to optimize the pial surface placement. Lastly, the PostFreeSurfer pipeline produced cortical surface models in GIFTI format and surface-related data in CIFTI format, and each subject's cortical surface was then registered to a common 32k_FS_LR mesh using “MSMAll” areal-feature-based cortical surface registration, which is a multimodal registration constrained by cortical T1w/T2w maps and resting-state network maps (Robinson et al., 2014, 2018; Glasser et al., 2016b).

Following cortical surface reconstruction, a single experienced individual performed a “SurfaceQC” review of the white and gray matter surface placement, informed by the T1w/T2w ratio maps (Glasser and Van Essen, 2011; Elam et al., 2021). Participants with more than minor (focal) issues were flagged for possible future editing and excluded from the cohort analyzed for the current study. This “SurfaceQC” review of the HCP-D data revealed some degradation of the accuracy of surface placement relative to expectations established by the HCP Young Adult project, which we subsequently traced to artifacts in the longer echoes. Therefore, to reduce the prevalence of surface segmentation errors in this developmental sample, we used the mean of the shortest two echoes (i.e., excluded the longest two of four echoes) as the T1w input to the HCP Pipelines (Elam et al., 2021).

T1w/T2w processing

T1w/T2w maps were generated as a cortical ribbon volume and a surface map using the methods described previously (Glasser and Van Essen, 2011; Glasser et al., 2013), together with more recent improvements (Glasser et al., 2014, 2021). Briefly, the T1w/T2w ratio in the voxels between the white and pial surfaces is mapped to the surface mesh in a way that emphasizes the middle layers and underemphasizes voxels near pial and gray/white surface to reduce partial volume effects (Glasser et al., 2013). Division of the T1w image by the T2w image mathematically cancels the signal intensity bias related to the sensitivity profile of the radio frequency receiver coils, which is the same in both images in the absence of subject head motion. Taking the ratio will also increase the contrast related to myelin content since both input images show myelin-related contrast (Glasser and Van Essen, 2011), inverted in the T2w image relative to the T1w image. Individual T1w/T2w maps were parcellated according to the HCP-multimodal atlas (Glasser et al., 2016a) using wb_command -cifti-parcellate in Connectome Workbench v1.4.2 (Marcus et al., 2011). This computes the average value of all vertices within a parcel and is thus a neuroanatomically constrained form of smoothing. T1w/T2w units are arbitrary, representing relative estimates of intracortical myelin content that are comparable within a consistently acquired study. Further, since the scaling of T1w and T2w images depends on the scanner and acquisition parameters, T1w/T2w units should not be directly compared across studies without appropriate harmonization.

B1+ transmit field correction of T1w/T2w myelin maps

T1w/T2w maps were initially developed for neuroanatomical analyses such as constraining the anatomic boundaries of cortical areas within individuals (Glasser et al., 2016a) and residual biases from the B1+ radio frequency transmit field were mitigated using the MyelinMap_BC approach (Glasser et al., 2013); however, this approach also removes genuine low spatial frequency individual differences in the T1w/T2w ratio. Given the increasing interest in using T1w/T2w maps to perform statistical comparisons across individuals and groups with other variables of interest, such as age (Kwon et al., 2018; Grydeland et al., 2019; Norbom et al., 2020), a new approach to B1+ transmit field bias correction was developed (Glasser et al., 2021). B1+ biases are influenced by the loading of the body coil (influenced by both the participant's head and body). Thus, variables such as head size, body size, and body mass index (BMI), which are often correlated with age and other variables of interest, can modulate B1+ bias and may lead to potentially spurious results in cross-subject statistical analyses of T1w/T2w maps.

In this study, we applied a novel, empirically validated “pseudo-transmit field” correction to mitigate B1+ bias in individual T1w/T2w maps, thereby reducing potentially spurious age-related differences in T1w/T2w values (Glasser et al., 2021). The B1+ correction approach in this study relies on computing a pseudo-transmit field based on the average across phase encoding directions of GRE images divided by SE images (mimicking the GRE T1w scan divided by SE T2w scan). First, a reference T1w/T2w map was generated at the group level by finding the scaling between the group average pseudo-transmit field and group average T1w/T2w map that minimizes the correlated left-right differences between the two maps (i.e., the clearly spurious left-right asymmetries). This reference group map was used to correct the individual maps. For the individual correction, the pseudo-transmit map was scaled to minimize the correlated differences between the individual's T1w/T2w map and the reference T1w/T2w map and the pseudo-transmit map (which includes all differences, not simply left-right ones, and is more robust at the individual level). Before estimating this correction, any residual B1– effects because of subject head motion between the T1w and T2w images were also removed using the scanner-computed B1– receive field. The pseudo-transmit field requires regularization by thresholding regions of T2*-related signal loss combined with spatial smoothing (with compensation for intensity changes induced by smoothing); it is then scaled to equal 1 at the value where the GRE/SE ratio corresponds to the flip angle prescribed by the scanner, a reference value that is determined at the group level. For more methodological details and validation of the correction, see Glasser et al. (2021), which demonstrates that this correction approach eliminates spurious relationships between the T1w/T2w ratio and B1+ bias modulated by body size.

Statistical analysis

Bayesian generalized additive models were fitted using the brms (Bürkner, 2018) interface to the Stan modeling language (Stan Development Team, 2021). Models were fitted separately to each parcellated cortical area to characterize different properties of age-related changes in T1w/T2w across the cerebral cortex. This approach samples from the posterior smooth function of T1w/T2w conditional on age and covariates, representing age-related changes in T1w/T2w, and the posterior derivative of the smooth function on age, representing the rate of change across the age range. Importantly, Bayesian thin-plate spline models allow us to estimate both linear and nonlinear age-related changes in T1w/T2w without specifying a linear or polynomial function a priori (Fahrmeir and Kneib, 2011; Hastie and Tibshirani, 2017; Wood, 2017). All statistical analyses were conducted in R version 3.5.1 (R Core Team, 2018).

We evaluated age-related change in T1w/T2w for each cortical parcel while controlling for participant sex, scanner, and several covariates related to the B1+ transmit field correction including the scanner transmit voltage, the mean of the pseudotransmit map, four regularization parameters [(1) T2* dropout threshold; (2) smoothing FWHM; (3) correction factor for smoothing's effect on the pseudotransmit field's intensities; and (4) the slope parameter of the correction], and a corrected T1w/T2w lateral ventricular CSF regressor (after excluding partial volume voxels and CSF flow effects) that are described in detail previously (Glasser et al., 2021). The brms syntax for the full Bayesian model for each cortical parcel was expressed as follows:

brms(T1w/T2w ∼ s(Age, k = −1, bs = 'tp') + Sex + Scanner + Transmit_Voltage + Mean_Pseudo_Transmit_Map + T2_Dropout_Threshold + Smoothing_FWHM_mm + Smoothing_Correction_Parameter + Correction_Slope + Corrected_CSF_Regressor)

Thin plate regression splines were used for the smoothing basis, and k was left at the default value of −1 which entails setting the basis dimension to 10. We retained the default priors set by brms; for regression coefficients this is a flat prior over the reals; for the intercept this is a Student's t distribution with df = 3, µ set at the mean of the outcome, and σ = 2.5; for the SD of the splines and SD of the error term this is a Student's t distribution with df = 3, µ = 0, and σ = 2.5. We ran four chains with 4500 total iterations (2000 warm-up) for each chain, yielding 10,000 posterior draws in total. The posterior smooth function for each cortical parcel was used to estimate other properties of age-related change in T1w/T2w, as described below.

Measurements of age-related change in T1w/T2w

Partial R2 values for age splines in each regional model were calculated by taking the difference between full models including an age spline and reduced models without age. These R2 values reflect the amount of variance in cortical T1w/T2w explained by age.

Annualized rate of change (AROC) in T1w/T2w

The annualized estimated rate of change in T1w/T2w was calculated for each cortical parcel by taking the difference between posterior samples for T1w/T2w at age 21 and age 8, and then dividing that difference by the number of years in the age range of the sample. This yielded a posterior distribution for this difference which we use to compute a point estimate using the median and 95% credible interval. The estimated rate of change in T1w/T2w was based on cross-sectional differences and does not reflect direct longitudinal estimates.

Linearity of age-related changes in T1w/T2w

Given prior work showing variation in the linearity or shape of neurodevelopmental changes during adolescence (Shaw et al., 2008), we used two complementary approaches to evaluate the linearity of age-related differences in T1w/T2w for each cortical parcel. Both approaches were dependent on the Bayesian spline models described above, which estimate both linear and nonlinear age-related differences in T1w/T2w without specifying a linear or polynomial function a priori (Wood, 2017). Therefore, to estimate the linearity of age curves regardless of functional form, we defined a measure of linearity as the mean absolute posterior second derivative of the smooth function of T1w/T2w on age.

First, the mean absolute posterior second derivative of the smooth function of T1w/T2w on age was calculated as a continuous measure of linearity for each cortical parcel using the curvish package in R (Flournoy, 2021). This measure reflects the rate of change in the slope of T1w/T2w age curves (e.g., how the rate of T1w/T2w change varies from 8 to 21 years).

Second, we conducted a separate, more conservative test evaluating whether a (nonlinear) spline model or a linear model of age-related changes in T1w/T2w provides greater out-of-sample predictive performance based on the leave-one-out-information criterion (LOOIC; Vehtari et al., 2017). Specifically, the smooth function from a generalized additive model was identified as nonlinear when the LOOIC difference between a Bayesian linear and spline model was at least 2.0 times the approximate standard error of its estimate (Bürkner, 2018). This is approximately analogous to a two-sided Z-test. The LOOIC indexes the out-of-sample predictive performance of a model. Therefore, cortical parcels for which the generalized additive model fits better are those in which we have credible evidence that this more complex, nonlinear form has greater out-of-sample predictive performance.

Age of peak and slowing T1w/T2w growth

The posterior age of peak growth (i.e., steepest slope) in T1w/T2w was estimated using the curvish package in R (Flournoy, 2021). Using the posterior of the smooth function and the method of finite differences at 100 equally spaced points, we evaluate the first derivative. Then, for each posterior sample we find the age at which the derivative is maximized, yielding a posterior density reflecting the probability of the slope being maximized at any given age.

To identify windows in development where the rate of age-related increases in T1w/T2w credibly slowed down, we computed the difference between the posterior derivative of the smooth at each evaluation point and the posterior derivative at the age of maximum median derivative. Regions of this 95% credible interval that did not include zero were used to mark regions of the smooth that were credibly less steep than the slope at the point of greatest maturation.

Identifying cortical areas with similar rate and shape of age-dependent T1w/T2w changes using clustering analysis

Based on prior literature characterizing variation in the maturational timing of cortical neurodevelopment across brain areas (Shaw et al., 2008; Somerville, 2016; Sydnor et al., 2021), we aimed to identify data-driven clusters of cortical parcels with similar T1w/T2w smooth functions on age. Further, rather than clustering parcels based on a single scalar property, we were interested in clustering cortical parcels based on latent dimensions of age-dependent T1w/T2w changes that capture subtle variation in the rate and shape of cortical development. To this end, we applied a functional latent mixture model using the funHDDC package in R (Bouveyron and Jacques, 2011), which consistently outperforms alternative clustering algorithms applied on spline coefficients (Bouveyron and Jacques, 2011; Schmutz et al., 2020).

First, for Bayesian models in each cortical parcel, we subtracted the posterior T1w/T2w estimates at age 8 from predicted T1w/T2w values at 500 age bins distributed equally between 8 and 21 years. This effectively removed intercept differences (i.e., starting T1w/T2w level) across cortical parcels, which were not of interest in this analysis. Second, we converted a brain region by age matrix into a functional data object for input into the funHDDC algorithm. Third, we ran funHDDC on the functional T1w/T2w trajectories, allowing the optimal number of clusters to vary freely between 2 and 7. Fourth and finally, we used the funHDDC::slopeHeuristic function to identify the most parsimonious model solution based on penalized log-likelihood (Schmutz and Bouveyron, 2021).

Sensitivity analyses

Cortical thickness

Prior studies have shown that T1w/T2w and cortical thickness measures are inversely correlated across the cerebral cortex (Glasser and Van Essen, 2011; Glasser et al., 2014; Natu et al., 2019). To demonstrate that our results were specific to age-related differences in T1w/T2w and are not attributable to concurrent changes in cortical thickness, we fit Bayesian spline models for each cortical parcel while controlling for estimates of cortical thickness (i.e., “*_V1_MR.thickness.32k_fs_LR.dscalar.nii”). We then calculated the spatial correlation between parcellated cortical maps of the AROC in T1w/T2w with and without controlling for cortical thickness.

Head motion

Head motion during structural MRI acquisition can impact image quality, inducing spurious artifacts in image intensity and morphometric measurements (Tisdall et al., 2016; Savalia et al., 2017). Age-related differences in motion (younger individuals commonly move more than older individuals) are a common concern when evaluating age-related changes in brain imaging measures (Satterthwaite et al., 2013; Reuter et al., 2015; Baum et al., 2018). Thus, we undertook sensitivity analyses in which participants with the highest motion estimates were excluded to test whether the primary results were evident in this lower-motion sample. To do so, we estimated the AROC in T1w/T2w after excluding 114 participants who exceeded the number of allowed TR reacquisitions for either the T1w (up to 30 allowed) or T2w (up to 25 allowed) acquisition, and thus had a considerable degree of ongoing head motion that was not fully corrected. For these participants, the number of k-space lines that still exceeded the motion threshold after the allowed number of reacquisitions is unknown, and thus the quality of these images is degraded relative to what could have been achieved if there was no limit placed on the number of reacquisitions (Tisdall et al., 2012). We then calculated the spatial correlation between parcellated cortical maps of the AROC in T1w/T2w with and without including participants with excess head motion.

Linking T1w/T2w development with the S-A axis of cortical organization

To evaluate whether cortical parcels exhibited differences in age-related changes in T1w/T2w according to their position along the S-A axis, we calculated the Spearman correlation between parcellated maps of cortical T1w/T2w development and S-A axis rankings. We used a spatial permutation test (described in detail below) to assess the significance of spatial correlations between each measure of age-related change in T1w/T2w described above (partial R2, rate of change, age of peak growth, degree of nonlinearity) and the S-A axis of cortical organization.

Briefly, S-A rankings were derived from cortical maps of ten fundamental brain features exhibiting systematic variation between lower-order sensorimotor areas and higher-order association areas (Sydnor et al., 2021). These cortical maps measured the anatomic hierarchy from T1w/T2w mapping (Glasser and Van Essen, 2011), the functional hierarchy quantified by the principal gradient of functional connectivity (Margulies et al., 2016), the evolutionary hierarchy quantified by macaque-to-human cortical expansion (J. Hill et al., 2010), allometric scaling quantified as the relative extent of areal scaling with overall brain size (Reardon et al., 2018), aerobic glycolysis quantified from positron emission tomography measures of oxygen consumption and glucose use (Vaishnavi et al., 2010), cerebral blood flow quantified by arterial spin labeling (Satterthwaite et al., 2014), gene expression quantified by the first principal component of brain-expressed genes (Burt et al., 2018), the first principal component of NeuroSynth meta-analytic decodings (Yarkoni et al., 2011), externopyramidization quantified as the ratio of supragranular pyramidal neuron soma size to infragranular pyramidal neuron soma size (Paquola et al., 2020b), and cortical thickness quantified from structural MRI (Human Connectome Project S1200 data). See Sydnor et al. (2021) for further details.

Spatial permutation test

The statistical significance of spatial correlations (Pspin) between parcellated cortical maps of age-related change in T1w/T2w and the archetypal S-A axis were assessed using a conservative parcel-based spatial permutation test that preserves the spatial covariance structure of cortical maps (Alexander-Bloch et al., 2018), as implemented by Váša et al. (2018). This permutation procedure generates a distribution of 10,000 spatial null maps randomly rotated on the spherical cortical surface to retain spatial contiguity and hemispheric symmetry of the original cortical maps. For each parcellated map showing age-related change in T1w/T2w (i.e., R2, rate of change, age of peak growth, nonlinearity, latent subgroups), we generated 10,000 randomly rotated “null” maps on the cortical surface. We then estimated the null distribution of Spearman correlation coefficients between the archetypal S-A axis map and the randomly rotated maps of T1w/T2w development. Permutation-based p-values (Pspin) were calculated based on the number of times the empirical correlation coefficient was higher than the permuted correlation coefficient.

False-positive error rate correction

When evaluating statistical significance of spatial correlations between cortical maps of age-related change in T1w/T2w and the archetypal S-A axis, the spatial permutation (“spin”) test provides family-wise control of multiple comparisons (Type I error) by accounting for the spatial and contralateral dependence of areal measurements in cortical maps (Alexander-Bloch et al., 2018). Further, we used the Holm correction (Holm, 1979), which is more powerful than Bonferroni and valid under arbitrary assumptions, to adjust the p-values across the 5 tests that are used to test the hypothesis that T1w/T2w development is linked with the S-A axis of cortical organization. The same correction was applied across all 15 spatial correlation tests between pairs of parcellated maps showing different properties of T1w/T2w development and the S-A axis (α = 0.05; Table 1). For descriptive purposes, we dichotomize parcels into those with or without a more complex functional form, and those with or without evident windows where the rate of age-related T1w/T2w growth credibly slows down, comprising 360 tests for each set. We do not attempt to control the family-wise false positive rate for the linearity decisions, or probability of making a sign error (Gelman and Carlin, 2014) when identifying windows where the rate of age-related T1w/T2w growth credibly slows down, as we report these as hypothesis-generating descriptions for which a higher Type I error rate (and lower Type II error rate) is preferable.

Data access and availability

Minimally preprocessed structural neuroimaging data from the HCP-D is available for download through the NIMH Data Archive (https://nda.nih.gov/). Parcellated cortical maps presented in this study are available for download through the BALSA database (https://balsa.wustl.edu/study/P2DmK).

Results

We analyzed high-resolution cortical T1w/T2w maps from 628 participants (ages 8–21 years old; Fig. 1A,B) to characterize the development of cortical T1w/T2w during youth after correcting for transmit field biases (Glasser et al., 2021). The posterior smooth function of T1w/T2w on age and the posterior derivative of the age function were estimated by fitting Bayesian thin-plate spline models across 180 cortical parcels in each hemisphere, thereby obtaining measures of the overall effect size, rate, and nonlinearity of age-related change in T1w/T2w across the neocortical sheet. For example, we illustrate age-related differences in T1w/T2w estimated for the left primary motor cortex using the posterior smooth function (Fig. 1C) and the posterior derivative (Fig. 1D) of T1w/T2w on age.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Charting T1w/T2w development during youth. A, Age/sex histogram of the current sample of 628 youth who completed structural neuroimaging as part of the HCP-D. B, T1w/T2w maps were parcellated using the HCP multimodal atlas (Glasser et al., 2016a) and averaged across participants. T1w/T2w units are arbitrary, representing relative estimates of intracortical myelin content that are comparable within a consistently acquired study. T1w/T2w units of 1.4 and 1.9 correspond to the second and 98th percentiles in this dataset. Data from the left primary motor cortex (highlighted with black border) are shown in panels C, D. C, We fit Bayesian generalized additive models with thin-plate splines to estimate different properties of age-related change in cortical T1w/T2w. Specifically, we estimated the posterior smooth function of T1w/T2w on age for each of the 360 cortical parcels (data shown for the left primary motor cortex, highlighted in panel B). The shaded area represents the 95% credible interval of the posterior smooth function. Navy blue segments of the posterior smooth function indicate the slope is credibly as steep as the maximum slope, while yellow segments of the posterior smooth function indicate the slope is credibly less steep than the maximum slope (see panel D). D, Bayesian generalized additive models were fitted to estimate the posterior derivative of the smooth function of T1w/T2w on age to characterize the rate of change, with higher values indicating a steeper slope of change per unit age and 0 representing a flat slope (i.e., no age-related change). The shaded area represents the 95% credible interval of posterior derivative estimates. The vertical line marks the age of the maximum median derivative (steepest slope). To identify windows where the rate of age-related T1w/T2w growth credibly slows down, we computed the difference between the posterior derivative of the smooth and the posterior derivative at the age of the maximum median derivative; regions of this 95% credible interval (data not shown) that did not include zero were used to mark regions of the smooth that were credibly less steep than the slope at the point of greatest maturation. Data are shown for an exemplar cortical parcel to illustrate our analytical approach.

For a descriptive visualization of cortical T1w/T2w development, we averaged data into three age groups (8–10 years approximating childhood, 14–16 years approximating mid-adolescence, and 19–21 years approximating emerging adulthood) to illustrate the spatial patterning of age-related differences in T1w/T2w from childhood to adulthood. The pattern of cortical T1w/T2w differences during adolescent development reflects the S-A axis reported by Sydnor et al. (2021; Fig. 2, inset right). The S-A hierarchy was evident in 8–10 year olds (n = 145; Fig. 2, left), with high T1w/T2w in sensory cortex and low T1w/T2w in association cortex. This hierarchy was present across all age groups, suggesting that the S-A hierarchy in cortical T1w/T2w is reinforced from childhood through adulthood. Relative to younger participants, 14–16 year olds (n = 154; Fig. 2, middle) had higher T1w/T2w across sensorimotor and intermediate association areas of the cerebral cortex. In 19–21 years olds, we observed a spatial pattern of cortical T1w/T2w extending from primary sensory to multimodal association cortex (n = 120; Fig. 2, right).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Waves of cortical T1w/T2w maturation from childhood through early adulthood. Parcellated T1w/T2w maps were averaged across participants within three age groups to illustrate the spatial patterning of age-related increases in T1w/T2w during youth. T1w/T2w units of 1.4 and 1.9 correspond to the second and 98th percentiles in this dataset. Age groups include 8–10 year olds (left; n = 145), 14–16 year olds (middle; n = 154), and 19–21 year olds (right; n = 120). Age-related increases in cortical T1w/T2w were observed across the cortical sheet, reinforcing a S-A hierarchy in cortical microstructure. The archetypal S-A axis (inset right) was adapted with permission from Sydnor et al. (2021).

Next, we evaluated the topography of the effect size and rate of age-related change in T1w/T2w from childhood through adulthood. In a multiple linear regression, all cortical parcels showed significant age-related differences in T1w/T2w after family-wise error correction (p < 0.05) with the exception of the bilateral subgenual cingulate (Brodmann area 25), which suffers from low SNR because of proximity of the ventricles and medial wall. Moreover, we found systematic variation in the amount of variance in T1w/T2w explained by age. Specifically, age explained a greater proportion of variance in T1w/T2w in sensorimotor areas such as right area V2 (R2 = 0.25; Fig. 3B) compared with association zones where age explained little variance in T1w/T2w, such as the left anterior cingulate cortex (area p32; R2 = 0.028; Fig. 3C). Figure 3D shows the spatial correlation between the variance in T1w/T2w explained by age (R2) and the cortical parcel's position along the archetypal S-A axis of cortical organization (r = −0.65, Pspin < 0.001).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Effect size estimates for regional Bayesian models of T1w/T2w development. A, Partial R2 values for age splines were estimated for each cortical parcel, representing the proportion of variance in T1w/T2w differences explained by age. B, Age-related increases in T1w/T2w were relatively strong in sensorimotor areas such as the right V2 (outlined and labeled in panel A), where age explained 25% of variance in T1w/T2w. C, By contrast, age-related increases in T1w/T2w were relatively weak in heteromodal and paralimbic association areas such as the left anterior cingulate cortex (area p32; outlined and labeled in panel A), where age explained 2.8% of variance in T1w/T2w. D, The variance in T1w/T2w explained by age (R2) was correlated with the cortical area's position along the archetypal S-A axis of cortical organization. In panel D, the blue data point corresponds to the exemplar sensorimotor parcel (V2) and the yellow data point corresponds to the exemplar association parcel (ACC) highlighted in panels A–C. ACC = anterior cingulate cortex; S-A = sensorimotor-association. Pspin is the permutation-based p-value calculated from a conservative parcel-based spatial permutation (“spin”) test. Data are shown for exemplar sensorimotor and association parcels to illustrate differences in T1w/T2w development. Models were estimated independently for all 180 cortical parcels in each hemisphere (360 in total).

We also observed spatial variation in the AROC in T1w/T2w estimated for each cortical parcel (Fig. 4A). Sensorimotor areas had relatively high AROC in T1w/T2w. Figure 4B shows the posterior smooth for the right motor cortex (area 4), an exemplar sensorimotor area with relatively high AROC in T1w/T2w (AROC = 0.025). Figure 4C shows the posterior smooth for the left anterior cingulate cortex (area p32), an exemplar association area with relatively low AROC in T1w/T2w (AROC = 0.005). From childhood through early adulthood, spatial variation in the AROC in T1w/T2w reflected S-A hierarchy (r = −0.65, Pspin < 0.001; Fig. 4D), with heavily myelinated sensorimotor areas showing rapid T1w/T2w increase relative to heteromodal and paralimbic association areas, which had slower age-related increases in T1w/T2w.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Annualized rate of change in T1w/T2w during youth. A, Annualized rate of change (AROC) in T1w/T2w was estimated for each cortical parcel. B, The AROC was relatively high in sensorimotor areas such as the right primary motor cortex (area 4; highlighted in panel A, right). C, In contrast, the AROC in T1w/T2w was relatively low in prefrontal and paralimbic association areas such as the left medial prefrontal cortex (area 9m; highlighted in panel A, left). D, The AROC in T1w/T2w was correlated with the cortical parcel's position along the archetypal S-A axis of cortical organization. The blue data point corresponds to the exemplar sensorimotor parcel (primary motor cortex) and the yellow data point corresponds to the exemplar association parcel (medial PFC) highlighted in panels A–C. Shaded areas in the plots of panels B, C represent the 95% credible interval of the posterior smooth function on age. Segments of the age curve where the slope is credibly less steep than the maximum slope (indicating a credibly reduced rate of change) are highlighted in yellow. PFC = prefrontal cortex; S-A = sensorimotor-association. Pspin is the permutation-based p-value calculated from a conservative parcel-based spatial permutation (“spin”) test. Data are shown for exemplar sensorimotor and association parcels to illustrate differences in T1w/T2w development. Models were estimated independently for all 180 cortical parcels in each hemisphere (360 in total).

Sensitivity analyses revealed high correspondence between the original parcellated cortical map of AROC in T1w/T2w and the parcellated map when controlling for cortical thickness in each cortical area (r = 0.96, Pspin < 0.001). We also observed a highly consistent spatial pattern of T1w/T2w development when excluding participants with excess head motion (n = 451; r = 0.97, Pspin < 0.001).

Next, we evaluated whether cortical parcels varied systematically in the linearity of age-related increases in T1w/T2w, as nonlinear T1w/T2w smooth functions might reflect diminishing age-related changes or a developmental inflection point. Figure 5A shows the linearity of age-related differences in T1w/T2w using the mean absolute posterior second derivative, which captures change in the slope over time. Figure 5B shows the posterior smooth function of T1w/T2w on age for the left V2, an exemplar sensorimotor area with relatively nonlinear age-related increases in T1w/T2w (mean absolute second derivative = 0.004). Figure 5C shows the posterior smooth function of T1w/T2w on age for the right anterior ventral insula, an exemplar association area with relatively linear age-related increases in T1w/T2w (mean absolute second derivative = 0.002). Spatial variation in the linearity of age-related increases in T1w/T2w was significantly correlated with the cortical parcel's position along the S-A axis (r = −0.45, Pspin < 0.001; Fig. 5D), with sensorimotor areas showing relatively nonlinear T1w/T2w increase relative to heteromodal and paralimbic association areas, which had relatively linear T1w/T2w increase with age. The LOOIC difference between linear and nonlinear spline models indicated that 63 of 360 cortical parcels (17.5%) had better out-of-sample predictive performance when including the spline. We found that 116 of 360 (32.2%) cortical parcels had segments of the age curve that were credibly less steep than the maximum slope, indicating a credible decrease in the slope of T1w/T2w maturation over time. This deceleration of age-related increases in T1w/T2w occurred primarily between 18 and 21 years old.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Nonlinearity of age-related changes in T1w/T2w during youth. A, The nonlinearity of age-related changes in T1w/T2w was estimated for each cortical parcel using the mean absolute posterior second derivative, where higher values indicate more nonlinear growth. B, Sensorimotor areas such as the left V2 (highlighted in panel A, left) had relatively nonlinear T1w/T2w growth, with the rate of change decreasing credibly by 17.5 years old (indicated by yellow segment). C, Cortical parcels in heteromodal and paralimbic association areas such as the right anterior ventral insula (highlighted in panel A, right) had relatively linear growth in T1w/T2w (constant slope). D, The nonlinearity of age-related changes in T1w/T2w was correlated with the cortical parcel's position along the archetypal S-A axis of cortical organization (Sydnor et al., 2021). Sensorimotor areas exhibited significantly higher nonlinearity in age-related increases in T1w/T2w compared with areas of association cortex. In panel D, the blue data point corresponds to the exemplar sensorimotor parcel (V2) and the yellow data point corresponds to the exemplar association parcel (anterior insula) highlighted in panels A–C. Shaded areas in the plots of panels B, C represent the 95% credible interval of the posterior smooth. Segments of the posterior smooth where the slope is credibly less steep than the maximum slope (indicating a credibly reduced rate of change) are highlighted in yellow. Pspin is the permutation-based p-value calculated from a conservative parcel-based spatial permutation (“spin”) test. Data are shown for exemplar sensorimotor and association parcels to illustrate differences in T1w/T2w development. Models were estimated independently for all 180 cortical parcels in each hemisphere (360 in total).

We further characterized cortical variation in the maturational timing of peak age-related increases in T1w/T2w based on the median posterior age where the first derivative was maximized (i.e., age of steepest slope; Fig. 6A). Figure 6B shows the shape and posterior age of peak T1w/T2w growth for the left V1, an exemplar sensorimotor area with a relatively early age of peak growth in T1w/T2w (median age = 9.6 years). Figure 6C shows the shape and posterior age of peak T1w/T2w growth for the right prefrontal area 8c, an exemplar association parcel with a relatively late age of peak growth in T1w/T2w (median age = 15.0 years). Spatial variation in the age of peak T1w/T2w growth from 8 to 21 years old closely followed the S-A axis of cortical organization (r = 0.60, Pspin < 0.001; Fig. 6D), with sensorimotor areas showing relatively early peak growth in T1w/T2w compared with heteromodal and paralimbic association areas, which had relatively late peak T1w/T2w growth. In general, areas with more linear T1w/T2w growth had greater uncertainty in the age of maximum slope, as illustrated by the width of the posterior distribution (Fig. 6B,C, upper inset). This may account for why the scatter (deviation from the best-fitting line) increases with S-A ranking in Figure 6D.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Age of peak T1w/T2w growth during youth. A, The age of peak growth (i.e., steepest slope) in T1w/T2w myelin was estimated for each cortical parcel as the median posterior age where the slope is at its maximum. B, Sensorimotor areas such as left V1 (highlighted in panel A) had relatively early age of peak T1w/T2w myelin growth (median age = 9.6 years). C, Association areas such as the right prefrontal cortex (area 8c; highlighted in panel A) had relatively later age of peak T1w/T2w myelin growth (median age = 15.0 years). D, The age of peak T1w/T2w myelin growth was correlated with the cortical parcel's position along the archetypal S-A axis of cortical organization. Sensorimotor areas exhibited a significantly earlier age of peak growth compared with areas of association cortex. In panel D, the blue data point corresponds to the exemplar sensorimotor parcel (V1) and the yellow data point corresponds to the exemplar association parcel (PFC) highlighted in panels A–C. Segments of the age curve where the slope is credibly less steep than the maximum slope (indicating a credibly reduced rate of change) are highlighted in yellow. Insets at the top of panels B, C show the posterior density distribution of the age of maximum slope for exemplar areas. The shaded blue area represents the 95% credible interval of the posterior distribution; black vertical lines mark the median of the posterior distribution (values projected on the cortical surface in panel A). PFC = prefrontal cortex. Pspin is the permutation-based p-value calculated from a conservative parcel-based spatial permutation (“spin”) test. Data are shown for exemplar cortical parcels to illustrate differences in T1w/T2w development. Models were estimated independently for all 180 cortical parcels in each hemisphere (360 in total).

We assessed whether groups of cortical parcels demonstrated similar temporal aspects of age-related change in T1w/T2w, as such groups of cortical parcels might reflect the development of coordinated patterns of brain connectivity. A functional latent mixture model was used to identify groups of cortical parcels exhibiting similar age-related patterns of T1w/T2w growth during youth. The most parsimonious model identified three latent clusters or groups of cortical parcels with graded differences in the rate of age-related changes in T1w/T2w (Fig. 7A). The spatial patterning of the functional latent clusters closely recapitulated the S-A axis (Sydnor et al., 2021; F = 166.6, p < 1.0 × 10−10): sensorimotor areas at one end of the S-A axis were grouped in latent cluster 1, intermediate association areas were grouped in latent cluster 2, and heteromodal/paralimbic association areas at the opposite end of the S-A axis were grouped in latent cluster three based on variation in the developmental timing of T1w/T2w changes (Fig. 7B). Cluster 1 (Fig. 7C, magenta) primarily included sensorimotor parcels, which had the highest rate of cortical T1w/T2w growth. Cluster 2 (Fig. 7C, orange) included heteromodal association parcels showing an intermediate rate of cortical T1w/T2w growth. Cluster 3 (Fig. 7C, yellow) included paralimbic association parcels in medial prefrontal and insular cortex, which had the lowest rate of cortical T1w/T2w growth.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Graded variation across data-driven clusters of cortical T1w/T2w development. A functional latent mixture model was used to identify groups of cortical parcels exhibiting similar age-related patterns of T1w/T2w growth during youth. A, The most parsimonious model identified three latent clusters or groups of cortical parcels with graded differences in the rate of age-related changes in T1w/T2w. The inset shows the penalized log-likelihood for each model solution from k = 2 to k = 7 clusters. The best fitting three-cluster solution is highlighted in red. B, Functional latent clusters closely aligned with the S-A axis of cortical organization. Sensorimotor areas, primarily in cluster 1 (magenta), had the highest rate of cortical T1w/T2w growth. We observed graded decreases in the rate of cortical T1w/T2w growth in cluster 2 (orange), which included heteromodal association parcels in frontoparietal and temporal cortex, and further in cluster 3 (yellow), which included paralimbic association parcels in the medial prefrontal and insular cortices. C, Regional cluster assignments based on T1w/T2w development reflect the S-A axis. L-L = log-likelihood.

We observed similar spatial patterning across each measurement of age-related changes in T1w/T2w. Table 1 shows the Spearman correlation coefficients and permutation-based p-values (Pspin) for each pairwise correlation among measurements. The high correspondence between these measures suggests that different temporal properties of T1w/T2w changes are yoked during adolescence and depend on the cortical parcel's position along the S-A axis.

View this table:
  • View inline
  • View popup
Table 1.

Correlated properties of T1w/T2w development and cortical S-A organization

Discussion

Here, we characterize the developmental timing of microstructural changes across the cerebral cortex during adolescence using recent advances in high-resolution T1w/T2w mapping applied to data from the HCP-D. This study is the first to evaluate the maturational timing of estimated cortical myelination in youth using B1+ transmit field-corrected T1w/T2w myelin maps, providing unbiased estimates of age-related differences in cortical microstructure. We demonstrate graded variation across the cortex in the timing of T1w/T2w changes during youth, with rapid age-related increases in sensorimotor cortex, and more gradual age-related increases in association cortex. This spatial pattern of microstructural brain development recapitulates the sensorimotor-to-association axis of cortical organization and plasticity during ontogeny. Moreover, our findings provide new evidence for protracted T1w/T2w changes in distributed association zones during youth, which could reflect ongoing cortical myelination underlying complex information processing and psychological functioning.

Graded variation in cortical T1w/T2w development during youth reflects the S-A axis of cortical organization

We observed graded variation in the rate and timing of age-related increases in cortical T1w/T2w from 8 to 21 years old. For each cortical parcel, Bayesian spline models provided posterior estimates of the AROC, the age of peak T1w/T2w growth, and the degree of nonlinearity in age-related changes. Each of these temporal properties of T1/T2w development followed a cortical topography reflecting S-A organization (Sydnor et al., 2021). Sensorimotor areas at one end of the S-A axis had a high rate of T1w/T2w growth, with early peak growth and relatively nonlinear age-related increases in T1w/T2w. Heteromodal and paralimbic association areas at the opposing end of the S-A axis had a slower rate of T1w/T2w changes, with late peak growth and relatively linear age-related increases in T1w/T2w. This coupling between measurements of age-related change in T1w/T2w suggest that the microstructural differentiation of cortical areas follows distinct developmental programs depending on the area's position along the S-A axis of cortical organization.

A functional latent mixture modeling procedure provided data-driven validation of the systematic variation we observed in T1w/T2w development. The most parsimonious clustering solution identified three clusters of cortical parcels based on the rate of age-dependent differences in T1w/T2w. This approach confirmed that the rate of cortical T1w/T2w growth during youth parallels the areal position along the S-A axis.

Our results provide new in vivo evidence aligned with histologic studies demonstrating topographic variation in the rate of cortical gray matter myelination during childhood and adulthood (Flechsig, 1901; Yakovlev and Lecours, 1967), underscoring a neurodevelopmental S-A axis. Graded variation in the rate of cortical myelination could reflect cascading sensitive windows for refining increasingly complex cognitive functions (Takesian and Hensch, 2013). Our findings are also consistent with comparative postmortem studies demonstrating prolonged myelination in human association cortex relative to chimpanzees (Miller et al., 2012), which could reflect an uniquely extended window of ongoing plasticity in cortical circuits underlying complex socioemotional and cognitive processing (Larsen and Luna, 2018).

Our main findings also converge with prior neuroimaging studies evaluating age-related differences in cortical myelin content using quantitative MRI methods such as magnetization transfer (MT) imaging, which demonstrate relatively late myelination in association cortices during adolescence into the third decade of life (Whitaker et al., 2016; Ziegler et al., 2019; Paquola et al., 2020a). Our findings are also consistent with a recent study evaluating age-related changes in T1w/T2w from 3–21 years old, which identified global increases in T1w/T2w across the neocortex (Norbom et al., 2020). Further, a lifespan study of 8–85 year olds characterized an early age of peak T1w/T2w growth in sensorimotor areas, and a later age of peak T1w/T2w growth in association areas (Grydeland et al., 2019). Despite their strengths, these studies did not control for B1+ transmit field bias, which impacts T1w/T2w estimates in a manner that is systematically related to age, sex, and body size (Glasser et al., 2021; MacLennan et al., 2022). Our study extends prior work by rigorously characterizing the pace and timing of cortical microstructural development in youth, using submillimeter resolution, B1+ bias-corrected T1w/T2w mapping, and mapping data to a multimodal cortical parcellation in a large, demographically diverse sample.

Spatial pattern of age-related increases in T1w/T2w remain consistent while controlling for cortical thickness

Cortical thickness and myelin content are systematically related, where thinner areas in sensorimotor cortex have higher myelin content, and thicker areas in paralimbic cortex have lower myelin content (but with exceptions such as the thick, heavily myelinated primary motor cortex; Glasser et al., 2014). Previous studies have observed variable rates of cortical thinning across the cortex during early childhood and adolescence (Shaw et al., 2008; Vandekar et al., 2015). Further, age-related changes in cortical thickness have been difficult to distinguish from changes in cortical myelin content under certain circumstances (Natu et al., 2019). In our study, age-related differences in T1w/T2w were not driven by cortical thinning, as we found highly consistent topography of age-related differences in T1w/T2w after controlling for cortical thickness in each parcel. This indicates that cortical thinning and T1w/T2w differences during youth are dissociable developmental processes.

Neurobiological mechanisms: the role of cortical myelin in regulating plasticity

While both myelin-sensitive imaging and postmortem histology have provided evidence for continuous increases in intracortical myelin content into adulthood (R.A. Hill et al., 2018; Hughes et al., 2018; Grydeland et al., 2019), the neurobiological mechanisms promoting cortical myelination throughout the lifespan remain largely unclear. Studies in animal models suggest that cortical myelination provides a mechanism to continuously tune the firing properties of neural circuits to support adaptive behavior within an individual's environment (Makinodan et al., 2012; Mount and Monje, 2017). Axonal myelination enhances neural signaling speed (McDougall et al., 2018) and plays an important role in regulating the timing of sensitive periods in neurocognitive development by providing structural constraints that limit local plasticity (Takesian and Hensch, 2013). Research using animal models has provided strong evidence for the role of myelin-related proteins such as Nogo-A in reducing plasticity by inhibiting neurite outgrowth and dendritic arborization in cortical circuits (Chen et al., 2000; McGee et al., 2005).

A role for cortical myelination in regulating plasticity during early sensitive periods for sensory cortex development has been demonstrated in animal models (Takesian and Hensch, 2013; Toyoizumi et al., 2013). However, the role of cortical myelination in regulating higher-order cognitive development in humans remains poorly understood. Recent work suggests that a large proportion of neocortical myelin ensheathes axons of parvalbumin-positive (PV+) fast-spiking inhibitory neurons (Micheva et al., 2016). Further, the myelination of PV+ interneurons is critical for maintaining excitation-inhibition (E/I) balance in maturing cortical sensory circuits (Benamer et al., 2020). The maturation of PV+ interneurons and E/I balance in neural circuits is a crucial factor regulating the opening and closure of sensitive periods in neurodevelopment (Takesian and Hensch, 2013; Toyoizumi et al., 2013; Takesian et al., 2018). Moreover, we speculate that cortical myelin development could be linked to organizational gradients of PV gene expression and the modulation of inhibitory circuitry throughout the cerebral cortex (Burt et al., 2018).

Limitations

Cortical MR signals are sensitive to many tissue properties including iron, myelin, cell density, and water content, and contrast images such as T1w/T2w generally represent a mix of these properties that can vary nonlinearly across laminar architecture (Carey et al., 2018; Glasser et al., 2022). Further, age-related differences in T1w/T2w could reflect a combination of these cortical tissue changes. For example, in addition to increased intracortical myelin content, age-related increases in T1w/T2w could also reflect an increase in tissue iron concentration (Larsen et al., 2020), the reduction in unmyelinated tissue processes such as dendrites, which could result in higher apparent myelin concentration per unit volume (expected to be associated with cortical thinning), or the dehydration of tissues (as occurs in early development).

Future directions

While this study characterized age-related changes in cortical T1w/T2w during youth aged 8–21 years, the cross-sectional design precluded inferences on intraindividual changes in cortical T1w/T2w. Although the HCP-D study contains a longitudinal component (Somerville et al., 2018), it was not feasible to incorporate the longitudinal measurements in this investigation because of ongoing data collection and processing. A more comprehensive lifespan approach could help identify milestones of cortical myelin development and senescence, expanding on previous work (Grydeland et al., 2019). The current study provides an important reference dataset characterizing adolescent cortical microstructural development that can serve as a benchmark for replication using the longitudinal subset of the HCP-D sample or other developmental longitudinal datasets such as ABCD (Casey et al., 2018). Such efforts would extend our understanding of cortical myelin development to allow us to evaluate the temporal relationships between cortical myelination and improving cognitive abilities. Future work could also investigate how socioeconomic factors, exposure to early life adversity, or other environmental contexts impact the rate of cortical myelination during youth, as well as the multifaceted relationships between developing cortical microstructure, connectivity, and psychological maturation.

Footnotes

  • This work was supported by National Institute of Mental Health (NIMH) Grants R01MH129493, R24MH108315, R24MH122820, U01MH109589, and U01MH109589-S1; the 14 National Institutes of Health (NIH) institutes and centers that support the NIH Blueprint for Neuroscience Research; the McDonnell Center for Systems Neuroscience at Washington University; the Office of the Provost at Washington University; and a Seed Grant from the Harvard Brain Institute Bipolar Disorder Seed Grant Program. Portions of this research were carried out at the Harvard Center for Brain Science using instrumentation supported by the NIH Shared Instrumentation Grant Program S10OD020039. This work was also supported by NIMH grants R01 MH60974 and the National Institute of Biomedical Imaging and Bioengineering (NIBIB) grant T32 EB021955. We thank Melanie Grad-Freilich for assistance with data collection and data access, Erin Reid for assistance with neuroimaging data quality assurance, and Valerie Sydnor for helpful discussion and sharing data related to the sensorimotor-association axis.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Graham L. Baum at gbaum{at}fas.harvard.edu

SfN exclusive license.

References

  1. ↵
    1. Alexander-Bloch AF,
    2. Shou H,
    3. Liu S,
    4. Satterthwaite TD,
    5. Glahn DC,
    6. Shinohara RT,
    7. Vandekar SN,
    8. Raznahan A
    (2018) On testing for spatial correspondence between maps of human brain structure and function. Neuroimage 178:540–551. doi:10.1016/j.neuroimage.2018.05.070 pmid:29860082
    OpenUrlCrossRefPubMed
  2. ↵
    1. Baum GL,
    2. Roalf DR,
    3. Cook PA,
    4. Ciric R,
    5. Rosen AFG,
    6. Xia C,
    7. Elliott MA,
    8. Ruparel K,
    9. Verma R,
    10. Tunç B,
    11. Gur RC,
    12. Gur RE,
    13. Bassett DS,
    14. Satterthwaite TD
    (2018) The impact of in-scanner head motion on structural connectivity derived from diffusion MRI. Neuroimage 173:275–286. doi:10.1016/j.neuroimage.2018.02.041 pmid:29486323
    OpenUrlCrossRefPubMed
  3. ↵
    1. Benamer N,
    2. Vidal M,
    3. Balia M,
    4. Angulo MC
    (2020) Myelination of parvalbumin interneurons shapes the function of cortical sensory inhibitory circuits. Nat Commun 11:5151. pmid:33051462
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bouveyron C,
    2. Jacques J
    (2011) Model-based clustering of time series in group-specific functional subspaces. Adv Data Anal Classif 5:281–300. doi:10.1007/s11634-011-0095-6
    OpenUrlCrossRef
  5. ↵
    1. Bürkner P-C
    (2018) Advanced Bayesian multilevel modeling with the R package brms. R J 10:395–411. doi:10.32614/RJ-2018-017
    OpenUrlCrossRef
  6. ↵
    1. Burt JB,
    2. Demirtaş M,
    3. Eckner WJ,
    4. Navejar NM,
    5. Ji JL,
    6. Martin WJ,
    7. Bernacchia A,
    8. Anticevic A,
    9. Murray JD
    (2018) Hierarchy of transcriptomic specialization across human cortex captured by structural neuroimaging topography. Nat Neurosci 21:1251–1259. doi:10.1038/s41593-018-0195-0 pmid:30082915
    OpenUrlCrossRefPubMed
  7. ↵
    1. Carey D,
    2. Caprini F,
    3. Allen M,
    4. Lutti A,
    5. Weiskopf N,
    6. Rees G,
    7. Callaghan MF,
    8. Dick F
    (2018) Quantitative MRI provides markers of intra-, inter-regional, and age-related differences in young adult cortical microstructure. Neuroimage 182:429–440. doi:10.1016/j.neuroimage.2017.11.066 pmid:29203455
    OpenUrlCrossRefPubMed
  8. ↵
    1. Casey BJ, et al
    . (2018) The adolescent brain cognitive development (ABCD) study: imaging acquisition across 21 sites. Dev Cogn Neurosci 32:43–54. doi:10.1016/j.dcn.2018.03.001 pmid:29567376
    OpenUrlCrossRefPubMed
  9. ↵
    1. Chen MS,
    2. Huber AB,
    3. van der Haar ME,
    4. Frank M,
    5. Schnell L,
    6. Spillmann AA,
    7. Christ F,
    8. Schwab ME
    (2000) Nogo-A is a myelin-associated neurite outgrowth inhibitor and an antigen for monoclonal antibody IN-1. Nature 403:434–439. doi:10.1038/35000219 pmid:10667796
    OpenUrlCrossRefPubMed
  10. ↵
    1. Dale AM,
    2. Fischl B,
    3. Sereno MI
    (1999) Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage 9:179–194. doi:10.1006/nimg.1998.0395 pmid:9931268
    OpenUrlCrossRefPubMed
  11. ↵
    1. Edwards LJ,
    2. Kirilina E,
    3. Mohammadi S,
    4. Weiskopf N
    (2018) Microstructural imaging of human neocortex in vivo. Neuroimage 182:184–206. doi:10.1016/j.neuroimage.2018.02.055 pmid:29588229
    OpenUrlCrossRefPubMed
  12. ↵
    1. Elam JS,
    2. Glasser MF,
    3. Harms MP,
    4. Sotiropoulos SN,
    5. Andersson JLR,
    6. Burgess GC,
    7. Curtiss SW,
    8. Oostenveld R,
    9. Larson-Prior LJ,
    10. Schoffelen J-M,
    11. Hodge MR,
    12. Cler EA,
    13. Marcus DM,
    14. Barch DM,
    15. Yacoub E,
    16. Smith SM,
    17. Ugurbil K,
    18. Van Essen DC
    (2021) The Human Connectome Project: a retrospective. Neuroimage 244:118543. doi:10.1016/j.neuroimage.2021.118543 pmid:34508893
    OpenUrlCrossRefPubMed
  13. ↵
    1. Fahrmeir L,
    2. Kneib T
    (2011) Bayesian smoothing and regression for longitudinal, spatial and event history data. Oxford: Oxford University Press.
  14. ↵
    1. Fischl B
    (2012) FreeSurfer. Neuroimage 62:774–781. doi:10.1016/j.neuroimage.2012.01.021 pmid:22248573
    OpenUrlCrossRefPubMed
  15. ↵
    1. Flechsig P
    (1901) Developmental (myelogenetic) localisation of the cerebral cortex in the human subject. Lancet 2:1027–1030.
    OpenUrl
  16. ↵
    1. Flournoy J
    (2021) jflournoy/curvish. Zenodo. Available at https://doi.org/10.5281/zenodo.5718790.
  17. ↵
    1. Fukunaga M,
    2. Li T-Q,
    3. van Gelderen P,
    4. de Zwart JA,
    5. Shmueli K,
    6. Yao B,
    7. Lee J,
    8. Maric D,
    9. Aronova MA,
    10. Zhang G,
    11. Leapman RD,
    12. Schenck JF,
    13. Merkle H,
    14. Duyn JH
    (2010) Layer-specific variation of iron content in cerebral cortex as a source of MRI contrast. Proc Natl Acad Sci U S A 107:3834–3839. doi:10.1073/pnas.0911177107 pmid:20133720
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Gelman A,
    2. Carlin J
    (2014) Beyond power calculations: assessing type S (sign) and type M (magnitude) errors. Perspect Psychol Sci 9:641–651. doi:10.1177/1745691614551642 pmid:26186114
    OpenUrlCrossRefPubMed
  19. ↵
    1. Glasser MF,
    2. Van Essen DC
    (2011) Mapping human cortical areas in vivo based on myelin content as revealed by T1- and T2-weighted MRI. J Neurosci 31:11597–11616. doi:10.1523/JNEUROSCI.2180-11.2011 pmid:21832190
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Glasser MF,
    2. Sotiropoulos SN,
    3. Wilson JA,
    4. Coalson TS,
    5. Fischl B,
    6. Andersson JL,
    7. Xu J,
    8. Jbabdi S,
    9. Webster M,
    10. Polimeni JR,
    11. Van Essen DC,
    12. Jenkinson M
    ; WU-Minn HCP Consortium (2013) The minimal preprocessing pipelines for the Human Connectome Project. Neuroimage 80:105–124. doi:10.1016/j.neuroimage.2013.04.127 pmid:23668970
    OpenUrlCrossRefPubMed
  21. ↵
    1. Glasser MF,
    2. Goyal MS,
    3. Preuss TM,
    4. Raichle ME,
    5. Van Essen DC
    (2014) Trends and properties of human cerebral cortex: correlations with cortical myelin content. Neuroimage 93:165–175. doi:10.1016/j.neuroimage.2013.03.060 pmid:23567887
    OpenUrlCrossRefPubMed
  22. ↵
    1. Glasser MF,
    2. Coalson TS,
    3. Robinson EC,
    4. Hacker CD,
    5. Harwell J,
    6. Yacoub E,
    7. Ugurbil K,
    8. Andersson J,
    9. Beckmann CF,
    10. Jenkinson M,
    11. Smith SM,
    12. Van Essen DC
    (2016a) A multi-modal parcellation of human cerebral cortex. Nature 536:171–178. doi:10.1038/nature18933 pmid:27437579
    OpenUrlCrossRefPubMed
  23. ↵
    1. Glasser MF,
    2. Smith SM,
    3. Marcus DS,
    4. Andersson J,
    5. Auerbach EJ,
    6. Behrens TEJ,
    7. Coalson TS,
    8. Harms MP,
    9. Jenkinson M,
    10. Moeller S,
    11. Robinson EC,
    12. Sotiropoulos SN,
    13. Xu J,
    14. Yacoub E,
    15. Ugurbil K,
    16. Van Essen DC
    (2016b) The Human Connectome Project's neuroimaging approach. Nat Neurosci 19:1175–1187. doi:10.1038/nn.4361 pmid:27571196
    OpenUrlCrossRefPubMed
  24. ↵
    1. Glasser MF,
    2. Coalson TS,
    3. Harms MP,
    4. Baum GL,
    5. Autio JA,
    6. Auerbach EJ,
    7. Xu J,
    8. Greve DN,
    9. Yacoub E,
    10. Essen DCV,
    11. Bock NA,
    12. Hayashi T
    (2022) Transmit field bias correction of T1w/T2w myelin maps. Neuroimage. doi:10.1016/j.neuroimage.2022.119360.
    OpenUrlCrossRef
  25. ↵
    1. Greve DN,
    2. Fischl B
    (2009) Accurate and robust brain image alignment using boundary-based registration. Neuroimage 48:63–72. doi:10.1016/j.neuroimage.2009.06.060 pmid:19573611
    OpenUrlCrossRefPubMed
  26. ↵
    1. Grydeland H,
    2. Vértes PE,
    3. Váša F,
    4. Romero-Garcia R,
    5. Whitaker K,
    6. Alexander-Bloch AF,
    7. Bjørnerud A,
    8. Patel AX,
    9. Sederevičius D,
    10. Tamnes CK,
    11. Westlye LT,
    12. White SR,
    13. Walhovd KB,
    14. Fjell AM,
    15. Bullmore ET
    (2019) Waves of maturation and senescence in micro-structural MRI markers of human cortical myelination over the lifespan. Cereb Cortex 29:1369–1381. doi:10.1093/cercor/bhy330 pmid:30590439
    OpenUrlCrossRefPubMed
  27. ↵
    1. Harms MP, et al
    . (2018) Extending the Human Connectome project across ages: imaging protocols for the Lifespan Development and Aging projects. Neuroimage 183:972–984. doi:10.1016/j.neuroimage.2018.09.060 pmid:30261308
    OpenUrlCrossRefPubMed
  28. ↵
    1. Hastie TJ,
    2. Tibshirani RJ
    (2017) Generalized additive models. Boca Raton: Routledge.
  29. ↵
    1. Hill J,
    2. Inder T,
    3. Neil J,
    4. Dierker D,
    5. Harwell J,
    6. Essen DV
    (2010) Similar patterns of cortical expansion during human development and evolution. Proc Natl Acad Sci U S A 107:13135–13140. doi:10.1073/pnas.1001229107 pmid:20624964
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Hill RA,
    2. Li AM,
    3. Grutzendler J
    (2018) Lifelong cortical myelin plasticity and age-related degeneration in the live mammalian brain. Nat Neurosci 21:683–695. doi:10.1038/s41593-018-0120-6 pmid:29556031
    OpenUrlCrossRefPubMed
  31. ↵
    1. Holm S
    (1979) A simple sequentially rejective multiple test procedure. Scand J Stat 6:65–70.
    OpenUrl
  32. ↵
    1. Hughes EG,
    2. Orthmann-Murphy JL,
    3. Langseth AJ,
    4. Bergles DE
    (2018) Myelin remodeling through experience-dependent oligodendrogenesis in the adult somatosensory cortex. Nat Neurosci 21:696–706. doi:10.1038/s41593-018-0121-5 pmid:29556025
    OpenUrlCrossRefPubMed
  33. ↵
    1. Huntenburg JM,
    2. Bazin P-L,
    3. Margulies DS
    (2018) Large-scale gradients in human cortical organization. Trends Cogn Sci 22:21–31. doi:10.1016/j.tics.2017.11.002 pmid:29203085
    OpenUrlCrossRefPubMed
  34. ↵
    1. Kwon D,
    2. Pfefferbaum A,
    3. Sullivan EV,
    4. Pohl KM
    (2018) Regional growth trajectories of cortical myelination in adolescents and young adults: longitudinal validation and functional correlates. Brain Imaging Behav 14:242–266.
    OpenUrl
  35. ↵
    1. Larsen B,
    2. Luna B
    (2018) Adolescence as a neurobiological critical period for the development of higher-order cognition. Neurosci Biobehav Rev 94:179–195. doi:10.1016/j.neubiorev.2018.09.005 pmid:30201220
    OpenUrlCrossRefPubMed
  36. ↵
    1. Larsen B,
    2. Bourque J,
    3. Moore TM,
    4. Adebimpe A,
    5. Calkins ME,
    6. Elliott MA,
    7. Gur RC,
    8. Gur RE,
    9. Moberg PJ,
    10. Roalf DR,
    11. Ruparel K,
    12. Turetsky BI,
    13. Vandekar SN,
    14. Wolf DH,
    15. Shinohara RT,
    16. Satterthwaite TD
    (2020) Longitudinal development of brain iron is linked to cognition in youth. J Neurosci 40:1810–1818. doi:10.1523/JNEUROSCI.2434-19.2020 pmid:31988059
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. MacLennan T,
    2. Seres P,
    3. Rickard J,
    4. Stolz E,
    5. Beaulieu C,
    6. Wilman AH
    (2022) Characterization of B1+ field variation in brain at 3 T using 385 healthy individuals across the lifespan. Magn Reson Med 87:960–971. doi:10.1002/mrm.29011
    OpenUrlCrossRef
  38. ↵
    1. Makinodan M,
    2. Rosen KM,
    3. Ito S,
    4. Corfas G
    (2012) A critical period for social experience–dependent oligodendrocyte maturation and myelination. Science 337:1357–1360. doi:10.1126/science.1220845 pmid:22984073
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Marcus DS,
    2. Harwell J,
    3. Olsen T,
    4. Hodge M,
    5. Glasser MF,
    6. Prior F,
    7. Jenkinson M,
    8. Laumann T,
    9. Curtiss SW,
    10. Van Essen DC
    (2011) Informatics and data mining tools and strategies for the human connectome project. Front Neuroinformatics 5:4.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Margulies DS,
    2. Ghosh SS,
    3. Goulas A,
    4. Falkiewicz M,
    5. Huntenburg JM,
    6. Langs G,
    7. Bezgin G,
    8. Eickhoff SB,
    9. Castellanos FX,
    10. Petrides M,
    11. Jefferies E,
    12. Smallwood J
    (2016) Situating the default-mode network along a principal gradient of macroscale cortical organization. Proc Natl Acad Sci U S A 113:12574–12579. doi:10.1073/pnas.1608282113 pmid:27791099
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. McDougall S,
    2. Riad WV,
    3. Silva-Gotay A,
    4. Tavares ER,
    5. Harpalani D,
    6. Li GL,
    7. Richardson HN
    (2018) Myelination of axons corresponds with faster transmission speed in the prefrontal cortex of developing male rats. eNeuro 5:ENEURO.0203-18.2018. doi:10.1523/ENEURO.0203-18.2018
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. McGee AW,
    2. Yang Y,
    3. Fischer QS,
    4. Daw NW,
    5. Strittmatter SM
    (2005) Experience-driven plasticity of visual cortex limited by myelin and Nogo receptor. Science 309:2222–2226. doi:10.1126/science.1114362 pmid:16195464
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Micheva KD,
    2. Wolman D,
    3. Mensh BD,
    4. Pax E,
    5. Buchanan J,
    6. Smith SJ,
    7. Bock DD
    (2016) A large fraction of neocortical myelin ensheathes axons of local inhibitory neurons. Elife 5:e15784. doi:10.7554/eLife.15784
    OpenUrlCrossRefPubMed
  44. ↵
    1. Miller DJ,
    2. Duka T,
    3. Stimpson CD,
    4. Schapiro SJ,
    5. Baze WB,
    6. McArthur MJ,
    7. Fobbs AJ,
    8. Sousa AMM,
    9. Šestan N,
    10. Wildman DE,
    11. Lipovich L,
    12. Kuzawa CW,
    13. Hof PR,
    14. Sherwood CC
    (2012) Prolonged myelination in human neocortical evolution. Proc Natl Acad Sci U S A 109:16480–16485. doi:10.1073/pnas.1117943109 pmid:23012402
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Mount CW,
    2. Monje M
    (2017) Wrapped to adapt: experience-dependent myelination. Neuron 95:743–756. doi:10.1016/j.neuron.2017.07.009 pmid:28817797
    OpenUrlCrossRefPubMed
  46. ↵
    1. Mugler JP,
    2. Brookeman JR
    (1990) Three-dimensional magnetization-prepared rapid gradient-echo imaging (3D MP RAGE). Magn Reson Med 15:152–157. doi:10.1002/mrm.1910150117 pmid:2374495
    OpenUrlCrossRefPubMed
  47. ↵
    1. Mugler JP,
    2. Bao S,
    3. Mulkern RV,
    4. Guttmann CRG,
    5. Robertson RL,
    6. Jolesz FA,
    7. Brookeman JR
    (2000) Optimized single-slab three-dimensional spin-echo MR imaging of the brain. Radiology 216:891–899. doi:10.1148/radiology.216.3.r00au46891 pmid:10966728
    OpenUrlCrossRefPubMed
  48. ↵
    1. Natu VS,
    2. Gomez J,
    3. Barnett M,
    4. Jeska B,
    5. Kirilina E,
    6. Jaeger C,
    7. Zhen Z,
    8. Cox S,
    9. Weiner KS,
    10. Weiskopf N,
    11. Grill-Spector K
    (2019) Apparent thinning of human visual cortex during childhood is associated with myelination. Proc Natl Acad Sci U S A 116:20750–20759. doi:10.1073/pnas.1904931116 pmid:31548375
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Norbom LB,
    2. Rokicki J,
    3. Alnæs D,
    4. Kaufmann T,
    5. Doan NT,
    6. Andreassen OA,
    7. Westlye LT,
    8. Tamnes CK
    (2020) Maturation of cortical microstructure and cognitive development in childhood and adolescence: a T1w/T2w ratio MRI study. Hum Brain Mapp 41:4676–4690. doi:10.1002/hbm.25149 pmid:32744409
    OpenUrlCrossRefPubMed
  50. ↵
    1. Paquola C,
    2. Bethlehem RA,
    3. Seidlitz J,
    4. Wagstyl K,
    5. Romero-Garcia R,
    6. Whitaker KJ,
    7. Vos de Wael R,
    8. Williams GB
    ; NSPN Consortium, Vértes PE, Margulies DS, Bernhardt B, Bullmore ET (2020a) Shifts in myeloarchitecture characterise adolescent development of cortical gradients. Elife 8:e50482. doi:10.7554/eLife.50482
    OpenUrlCrossRef
  51. ↵
    1. Paquola C,
    2. Seidlitz J,
    3. Benkarim O,
    4. Royer J,
    5. Klimes P,
    6. Bethlehem RAI,
    7. Larivière S,
    8. Vos de Wael R,
    9. Rodríguez-Cruces R,
    10. Hall JA,
    11. Frauscher B,
    12. Smallwood J,
    13. Bernhardt BC
    (2020b) A multi-scale cortical wiring space links cellular architecture and functional dynamics in the human brain. PLoS Biol 18:e3000979. doi:10.1371/journal.pbio.3000979 pmid:33253185
    OpenUrlCrossRefPubMed
  52. ↵
    R Core Team (2018) R: a language and environment for statistical computing, version 3.5.1. Vienna: R Foundation for Statistical Computing.
  53. ↵
    1. Reardon PK,
    2. Seidlitz J,
    3. Vandekar S,
    4. Liu S,
    5. Patel R,
    6. Park MTM,
    7. Alexander-Bloch A,
    8. Clasen LS,
    9. Blumenthal JD,
    10. Lalonde FM,
    11. Giedd JN,
    12. Gur RC,
    13. Gur RE,
    14. Lerch JP,
    15. Chakravarty MM,
    16. Satterthwaite TD,
    17. Shinohara RT,
    18. Raznahan A
    (2018) Normative brain size variation and brain shape diversity in humans. Science 360:1222–1227. doi:10.1126/science.aar2578 pmid:29853553
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Reuter M,
    2. Tisdall MD,
    3. Qureshi A,
    4. Buckner RL,
    5. van der Kouwe AJW,
    6. Fischl B
    (2015) Head motion during MRI acquisition reduces gray matter volume and thickness estimates. Neuroimage 107:107–115. doi:10.1016/j.neuroimage.2014.12.006 pmid:25498430
    OpenUrlCrossRefPubMed
  55. ↵
    1. Robinson EC,
    2. Jbabdi S,
    3. Glasser MF,
    4. Andersson J,
    5. Burgess GC,
    6. Harms MP,
    7. Smith SM,
    8. Van Essen DC,
    9. Jenkinson M
    (2014) MSM: a new flexible framework for multimodal surface matching. Neuroimage 100:414–426. doi:10.1016/j.neuroimage.2014.05.069 pmid:24939340
    OpenUrlCrossRefPubMed
  56. ↵
    1. Robinson EC, et al
    . (2018) Multimodal surface matching with higher-order smoothness constraints. Neuroimage 167:453–465. doi:10.1016/j.neuroimage.2017.10.037 pmid:29100940
    OpenUrlCrossRefPubMed
  57. ↵
    1. Satterthwaite TD,
    2. Wolf DH,
    3. Ruparel K,
    4. Erus G,
    5. Elliott MA,
    6. Eickhoff SB,
    7. Gennatas ED,
    8. Jackson C,
    9. Prabhakaran K,
    10. Smith A,
    11. Hakonarson H,
    12. Verma R,
    13. Davatzikos C,
    14. Gur RE,
    15. Gur RC
    (2013) Heterogeneous impact of motion on fundamental patterns of developmental changes in functional connectivity during youth. Neuroimage 83:45–57. doi:10.1016/j.neuroimage.2013.06.045 pmid:23792981
    OpenUrlCrossRefPubMed
  58. ↵
    1. Satterthwaite TD,
    2. Shinohara RT,
    3. Wolf DH,
    4. Hopson RD,
    5. Elliott MA,
    6. Vandekar SN,
    7. Ruparel K,
    8. Calkins ME,
    9. Roalf DR,
    10. Gennatas ED,
    11. Jackson C,
    12. Erus G,
    13. Prabhakaran K,
    14. Davatzikos C,
    15. Detre JA,
    16. Hakonarson H,
    17. Gur RC,
    18. Gur RE
    (2014) Impact of puberty on the evolution of cerebral perfusion during adolescence. Proc Natl Acad Sci U S A 111:8643–8648. doi:10.1073/pnas.1400178111 pmid:24912164
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Savalia NK,
    2. Agres PF,
    3. Chan MY,
    4. Feczko EJ,
    5. Kennedy KM,
    6. Wig GS
    (2017) Motion-related artifacts in structural brain images revealed with independent estimates of in-scanner head motion. Hum Brain Mapp 38:472–492. doi:10.1002/hbm.23397 pmid:27634551
    OpenUrlCrossRefPubMed
  60. ↵
    1. Schmutz A,
    2. Bouveyron C
    (2021) funHDDC: univariate and multivariate model-based clustering in group-specific functional subspaces. https://CRAN.R-project.org/package=funHDDC.
  61. ↵
    1. Schmutz A,
    2. Jacques J,
    3. Bouveyron C,
    4. Chèze L,
    5. Martin P
    (2020) Clustering multivariate functional data in group-specific functional subspaces. Comput Stat 35:1101–1131. doi:10.1007/s00180-020-00958-4
    OpenUrlCrossRef
  62. ↵
    1. Shams Z,
    2. Norris DG,
    3. Marques JP
    (2019) A comparison of in vivo MRI based cortical myelin mapping using T1w/T2w and R1 mapping at 3T. PLoS One 14:e0218089. doi:10.1371/journal.pone.0218089 pmid:31269041
    OpenUrlCrossRefPubMed
  63. ↵
    1. Shaw P,
    2. Kabani NJ,
    3. Lerch JP,
    4. Eckstrand K,
    5. Lenroot R,
    6. Gogtay N,
    7. Greenstein D,
    8. Clasen L,
    9. Evans A,
    10. Rapoport JL,
    11. Giedd JN,
    12. Wise SP
    (2008) Neurodevelopmental trajectories of the human cerebral cortex. J Neurosci 28:3586–3594. doi:10.1523/JNEUROSCI.5309-07.2008 pmid:18385317
    OpenUrlAbstract/FREE Full Text
  64. ↵
    1. Somerville LH
    (2016) Searching for signatures of brain maturity: what are we searching for? Neuron 92:1164–1167. doi:10.1016/j.neuron.2016.10.059 pmid:28009272
    OpenUrlCrossRefPubMed
  65. ↵
    1. Somerville LH,
    2. Bookheimer SY,
    3. Buckner RL,
    4. Burgess GC,
    5. Curtiss SW,
    6. Dapretto M,
    7. Elam JS,
    8. Gaffrey MS,
    9. Harms MP,
    10. Hodge C,
    11. Kandala S,
    12. Kastman EK,
    13. Nichols TE,
    14. Schlaggar BL,
    15. Smith SM,
    16. Thomas KM,
    17. Yacoub E,
    18. Van Essen DC,
    19. Barch DM
    (2018) The lifespan human connectome project in development: a large-scale study of brain connectivity development in 5–21 year olds. Neuroimage 183:456–468. doi:10.1016/j.neuroimage.2018.08.050 pmid:30142446
    OpenUrlCrossRefPubMed
  66. ↵
    Stan Development Team (2021) Stan modeling language users guide and reference manual, v2.28. https://mc-stan.org.
  67. ↵
    1. Stüber C,
    2. Morawski M,
    3. Schäfer A,
    4. Labadie C,
    5. Wähnert M,
    6. Leuze C,
    7. Streicher M,
    8. Barapatre N,
    9. Reimann K,
    10. Geyer S,
    11. Spemann D,
    12. Turner R
    (2014) Myelin and iron concentration in the human brain: a quantitative study of MRI contrast. Neuroimage 93:95–106. doi:10.1016/j.neuroimage.2014.02.026
    OpenUrlCrossRefPubMed
  68. ↵
    1. Sydnor VJ,
    2. Larsen B,
    3. Bassett DS,
    4. Alexander-Bloch A,
    5. Fair DA,
    6. Liston C,
    7. Mackey AP,
    8. Milham MP,
    9. Pines A,
    10. Roalf DR,
    11. Seidlitz J,
    12. Xu T,
    13. Raznahan A,
    14. Satterthwaite TD
    (2021) Neurodevelopment of the association cortices: patterns, mechanisms, and implications for psychopathology. Neuron 109:2820–2846. doi:10.1016/j.neuron.2021.06.016 pmid:34270921
    OpenUrlCrossRefPubMed
  69. ↵
    1. Takesian AE,
    2. Hensch TK
    (2013) Balancing plasticity/stability across brain development. Prog Brain Res 207:3–34. pmid:24309249
    OpenUrlCrossRefPubMed
  70. ↵
    1. Takesian AE,
    2. Bogart LJ,
    3. Lichtman JW,
    4. Hensch TK
    (2018) Inhibitory circuit gating of auditory critical-period plasticity. Nat Neurosci 21:218–227. doi:10.1038/s41593-017-0064-2 pmid:29358666
    OpenUrlCrossRefPubMed
  71. ↵
    1. Tisdall MD,
    2. Hess AT,
    3. Reuter M,
    4. Meintjes EM,
    5. Fischl B,
    6. van der Kouwe AJW
    (2012) Volumetric navigators for prospective motion correction and selective reacquisition in neuroanatomical MRI. Magn Reson Med 68:389–399. doi:10.1002/mrm.23228 pmid:22213578
    OpenUrlCrossRefPubMed
  72. ↵
    1. Tisdall MD,
    2. Reuter M,
    3. Qureshi A,
    4. Buckner RL,
    5. Fischl B,
    6. van der Kouwe AJW
    (2016) Prospective motion correction with volumetric navigators (vNavs) reduces the bias and variance in brain morphometry induced by subject motion. Neuroimage 127:11–22. doi:10.1016/j.neuroimage.2015.11.054 pmid:26654788
    OpenUrlCrossRefPubMed
  73. ↵
    1. Toyoizumi T,
    2. Miyamoto H,
    3. Yazaki-Sugiyama Y,
    4. Atapour N,
    5. Hensch TK,
    6. Miller KD
    (2013) A theory of the transition to critical period plasticity: inhibition selectively suppresses spontaneous activity. Neuron 80:51–63. doi:10.1016/j.neuron.2013.07.022 pmid:24094102
    OpenUrlCrossRefPubMed
  74. ↵
    1. Vaishnavi SN,
    2. Vlassenko AG,
    3. Rundle MM,
    4. Snyder AZ,
    5. Mintun MA,
    6. Raichle ME
    (2010) Regional aerobic glycolysis in the human brain. Proc Natl Acad Sci U S A 107:17757–17762. doi:10.1073/pnas.1010459107 pmid:20837536
    OpenUrlAbstract/FREE Full Text
  75. ↵
    1. van der Kouwe AJW,
    2. Benner T,
    3. Salat DH,
    4. Fischl B
    (2008) Brain morphometry with multiecho MPRAGE. Neuroimage 40:559–569. doi:10.1016/j.neuroimage.2007.12.025 pmid:18242102
    OpenUrlCrossRefPubMed
  76. ↵
    1. Van Essen DC,
    2. Smith SM,
    3. Barch DM,
    4. Behrens TEJ,
    5. Yacoub E,
    6. Ugurbil K
    ; WU-Minn HCP Consortium (2013) The WU-Minn Human Connectome Project: an overview. Neuroimage 80:62–79. doi:10.1016/j.neuroimage.2013.05.041 pmid:23684880
    OpenUrlCrossRefPubMed
  77. ↵
    1. Vandekar SN,
    2. Shinohara RT,
    3. Raznahan A,
    4. Roalf DR,
    5. Ross M,
    6. DeLeo N,
    7. Ruparel K,
    8. Verma R,
    9. Wolf DH,
    10. Gur RC,
    11. Gur RE,
    12. Satterthwaite TD
    (2015) Topologically dissociable patterns of development of the human cerebral cortex. J Neurosci 35:599–609. doi:10.1523/JNEUROSCI.3628-14.2015 pmid:25589754
    OpenUrlAbstract/FREE Full Text
  78. ↵
    1. Váša F,
    2. Seidlitz J,
    3. Romero-Garcia R,
    4. Whitaker KJ,
    5. Rosenthal G,
    6. Vértes PE,
    7. Shinn M,
    8. Alexander-Bloch A,
    9. Fonagy P,
    10. Dolan RJ,
    11. Jones PB,
    12. Goodyer IM
    ; NSPN Consortium, Sporns O, Bullmore ET (2018) Adolescent tuning of association cortex in human structural brain networks. Cereb Cortex 28:281–294. doi:10.1093/cercor/bhx249 pmid:29088339
    OpenUrlCrossRefPubMed
  79. ↵
    1. Vehtari A,
    2. Gelman A,
    3. Gabry J
    (2017) Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC. Stat Comput 27:1413–1432. doi:10.1007/s11222-016-9696-4
    OpenUrlCrossRef
  80. ↵
    1. Whitaker KJ,
    2. Vértes PE,
    3. Romero-Garcia R,
    4. Váša F,
    5. Moutoussis M,
    6. Prabhu G,
    7. Weiskopf N,
    8. Callaghan MF,
    9. Wagstyl K,
    10. Rittman T,
    11. Tait R,
    12. Ooi C,
    13. Suckling J,
    14. Inkster B,
    15. Fonagy P,
    16. Dolan RJ,
    17. Jones PB,
    18. Goodyer IM
    ; NSPN Consortium, Bullmore ET (2016) Adolescence is associated with genomically patterned consolidation of the hubs of the human brain connectome. Proc Natl Acad Sci U S A 113:9105–9110. doi:10.1073/pnas.1601745113 pmid:27457931
    OpenUrlAbstract/FREE Full Text
  81. ↵
    1. Wood SN
    (2017) Generalized additive models: an introduction with R, Ed 2. Boca Raton: CRC.
  82. ↵
    1. Yakovlev P,
    2. Lecours A
    (1967) The myelogenetic cycles of regional maturation of the brain. In: Resional development of the brain in early life (Minkowski A, eds), pp 3–70. Oxford: Blackwell.
  83. ↵
    1. Yarkoni T,
    2. Poldrack RA,
    3. Nichols TE,
    4. Van Essen DC,
    5. Wager TD
    (2011) Large-scale automated synthesis of human functional neuroimaging data. Nat Methods 8:665–670. doi:10.1038/nmeth.1635 pmid:21706013
    OpenUrlCrossRefPubMed
  84. ↵
    1. Ziegler G,
    2. Hauser TU,
    3. Moutoussis M,
    4. Bullmore ET,
    5. Goodyer IM,
    6. Fonagy P,
    7. Jones PB,
    8. Lindenberger U,
    9. Dolan RJ
    ; NSPN Consortium (2019) Compulsivity and impulsivity traits linked to attenuated developmental frontostriatal myelination trajectories. Nat Neurosci 22:992–999. doi:10.1038/s41593-019-0394-3 pmid:31086316
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 42 (29)
Journal of Neuroscience
Vol. 42, Issue 29
20 Jul 2022
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Masthead (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Graded Variation in T1w/T2w Ratio during Adolescence: Measurement, Caveats, and Implications for Development of Cortical Myelin
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Graded Variation in T1w/T2w Ratio during Adolescence: Measurement, Caveats, and Implications for Development of Cortical Myelin
Graham L. Baum, John C. Flournoy, Matthew F. Glasser, Michael P. Harms, Patrick Mair, Ashley F. P. Sanders, Deanna M. Barch, Randy L. Buckner, Susan Bookheimer, Mirella Dapretto, Stephen Smith, Kathleen M. Thomas, Essa Yacoub, David C. Van Essen, Leah H. Somerville
Journal of Neuroscience 20 July 2022, 42 (29) 5681-5694; DOI: 10.1523/JNEUROSCI.2380-21.2022

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Graded Variation in T1w/T2w Ratio during Adolescence: Measurement, Caveats, and Implications for Development of Cortical Myelin
Graham L. Baum, John C. Flournoy, Matthew F. Glasser, Michael P. Harms, Patrick Mair, Ashley F. P. Sanders, Deanna M. Barch, Randy L. Buckner, Susan Bookheimer, Mirella Dapretto, Stephen Smith, Kathleen M. Thomas, Essa Yacoub, David C. Van Essen, Leah H. Somerville
Journal of Neuroscience 20 July 2022, 42 (29) 5681-5694; DOI: 10.1523/JNEUROSCI.2380-21.2022
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • adolescence
  • development
  • MRI
  • myelin
  • T1w/T2w

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Regional Excitatory–Inhibitory Balance Relates to Self-Reference Effect on Recollection via the Precuneus/Posterior Cingulate Cortex–Medial Prefrontal Cortex Connectivity
  • Attention Alters Population Spatial Frequency Tuning
  • Sequence Diversity and Encoded Enzymatic Differences of Monocistronic L1 ORF2 mRNA Variants in the Aged Normal and Alzheimer's Disease Brain
Show more Research Articles

Development/Plasticity/Repair

  • Change of spiny neuron structure in the basal ganglia song circuit and its regulation by miR-9 during song development
  • Stereotyped Spatiotemporal Dynamics of Spontaneous Activity in Visual Cortex Prior to Eye Opening
  • The epigenetic reader PHF23 is required for embryonic neurogenesis
Show more Development/Plasticity/Repair
  • Home
  • Alerts
  • Follow SFN on BlueSky
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Notice
  • Contact
  • Accessibility
(JNeurosci logo)
(SfN logo)

Copyright © 2025 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.