Brain-age prediction: Systematic evaluation of site effects, and sample age range and size.

Structural neuroimaging data have been used to compute an estimate of the biological age of the brain (brain-age) which has been associated with other biologically and behaviorally meaningful measures of brain development and aging. The ongoing research interest in brain-age has highlighted the need for robust and publicly available brain-age models pre-trained on data from large samples of healthy individuals. To address this need we have previously released a developmental brain-age model. Here we expand this work to develop, empirically validate, and disseminate a pre-trained brain-age model to cover most of the human lifespan. To achieve this, we selected the best-performing model after systematically examining the impact of seven site harmonization strategies, age range, and sample size on brain-age prediction in a discovery sample of brain morphometric measures from 35,683 healthy individuals (age range: 5-90 years; 53.59% female). The pre-trained models were tested for cross-dataset generalizability in an independent sample comprising 2101 healthy individuals (age range: 8-80 years; 55.35% female) and for longitudinal consistency in a further sample comprising 377 healthy individuals (age range: 9-25 years; 49.87% female). This empirical examination yielded the following findings: (1) the accuracy of age prediction from morphometry data was higher when no site harmonization was applied; (2) dividing the discovery sample into two age-bins (5-40 and 40-90 years) provided a better balance between model accuracy and explained age variance than other alternatives; (3) model accuracy for brain-age prediction plateaued at a sample size exceeding 1600 participants. These findings have been incorporated into CentileBrain (https://centilebrain.org/#/brainAGE2), an open-science, web-based platform for individualized neuroimaging metrics.

Precision Assessment of Real-World Associations Between Stress and Sleep Duration Using Actigraphy Data Collected Continuously for an Academic Year: Individual-Level Modeling Study.

BACKGROUND: Heightened stress and insufficient sleep are common in the transition to college, often co-occur, and have both been linked to negative health outcomes. A challenge concerns disentangling whether perceived stress precedes or succeeds changes in sleep. These day-to-day associations may vary across individuals, but short study periods and group-level analyses in prior research may have obscured person-specific phenotypes. OBJECTIVE: This study aims to obtain stable estimates of lead-lag associations between perceived stress and objective sleep duration in the individual, unbiased by the group, by developing an individual-level linear model that can leverage intensive longitudinal data while remaining parsimonious. METHODS: In total, 55 college students (n=6, 11% second-year students and n=49, 89% first-year students) volunteered to provide daily self-reports of perceived stress via a smartphone app and wore an actigraphy wristband for the estimation of daily sleep duration continuously throughout the academic year (median usable daily observations per participant: 178, IQR 65.5). The individual-level linear model, developed in a Bayesian framework, included the predictor and outcome of interest and a covariate for the day of the week to account for weekly patterns. We validated the model on the cohort of second-year students (n=6, used as a pilot sample) by applying it to variables expected to correlate positively within individuals: objective sleep duration and self-reported sleep quality. The model was then applied to the fully independent target sample of first-year students (n=49) for the examination of bidirectional associations between daily stress levels and sleep duration. RESULTS: Proof-of-concept analyses captured expected associations between objective sleep duration and subjective sleep quality in every pilot participant. Target analyses revealed negative associations between sleep duration and perceived stress in most of the participants (45/49, 92%), but their temporal association varied. Of the 49 participants, 19 (39%) showed a significant association (probability of direction>0.975): 8 (16%) showed elevated stress in the day associated with shorter sleep later that night, 5 (10%) showed shorter sleep associated with elevated stress the next day, and 6 (12%) showed both directions of association. Of note, when analyzed using a group-based multilevel model, individual estimates were systematically attenuated, and some even reversed sign. CONCLUSIONS: The dynamic interplay of stress and sleep in daily life is likely person specific. Paired with intensive longitudinal data, our individual-level linear model provides a precision framework for the estimation of stable real-world behavioral and psychological dynamics and may support the personalized prioritization of intervention targets for health and well-being.

Organization of the human cerebral cortex estimated within individuals: networks, global topography, and function.

The cerebral cortex is populated by specialized regions that are organized into networks. Here we estimated networks from functional MRI (fMRI) data in intensively sampled participants. The procedure was developed in two participants (scanned 31 times) and then prospectively applied to 15 participants (scanned 8-11 times). Analysis of the networks revealed a global organization. Locally organized first-order sensory and motor networks were surrounded by spatially adjacent second-order networks that linked to distant regions. Third-order networks possessed regions distributed widely throughout association cortex. Regions of distinct third-order networks displayed side-by-side juxtapositions with a pattern that repeated across multiple cortical zones. We refer to these as supra-areal association megaclusters (SAAMs). Within each SAAM, two candidate control regions were adjacent to three separate domain-specialized regions. Response properties were explored with task data. The somatomotor and visual networks responded to body movements and visual stimulation, respectively. Second-order networks responded to transients in an oddball detection task, consistent with a role in orienting to salient events. The third-order networks, including distinct regions within each SAAM, showed two levels of functional specialization. Regions linked to candidate control networks responded to working memory load across multiple stimulus domains. The remaining regions dissociated across language, social, and spatial/episodic processing domains. These results suggest that progressively higher-order networks nest outward from primary sensory and motor cortices. Within the apex zones of association cortex, there is specialization that repeatedly divides domain-flexible from domain-specialized regions. We discuss implications of these findings, including how repeating organizational motifs may emerge during development.NEW & NOTEWORTHY The organization of cerebral networks was estimated within individuals with intensive, repeat sampling of fMRI data. A hierarchical organization emerged in each individual that delineated first-, second-, and third-order cortical networks. Regions of distinct third-order association networks consistently exhibited side-by-side juxtapositions that repeated across multiple cortical zones, with clear and robust functional specialization among the embedded regions.

Human striatal association megaclusters.

The striatum receives projections from multiple regions of the cerebral cortex consistent with the role of the basal ganglia in diverse motor, affective, and cognitive functions. Within the striatum, the caudate receives projections from association cortex, including multiple distinct regions of prefrontal cortex. Building on recent insights about the details of how juxtaposed cortical networks are specialized for distinct aspects of higher-order cognition, we revisited caudate organization using within-individual precision neuroimaging initially in two intensively scanned individuals (each scanned 31 times). Results revealed that the caudate has side-by-side regions that are coupled to at least five distinct distributed association networks, paralleling the organization observed in the cerebral cortex. We refer to these spatial groupings of regions as striatal association megaclusters. Correlation maps from closely juxtaposed seed regions placed within the megaclusters recapitulated the five distinct cortical networks, including their multiple spatially distributed regions. Striatal association megaclusters were explored in 15 additional participants (each scanned at least 8 times), finding that their presence generalizes to new participants. Analysis of the laterality of the regions within the megaclusters further revealed that they possess asymmetries paralleling their cortical counterparts. For example, caudate regions linked to the language network were left lateralized. These results extend the general notion of parallel specialized basal ganglia circuits with the additional discovery that, even within the caudate, there is fine-grained separation of multiple distinct higher-order networks that reflects the organization and lateralization found in the cerebral cortex.NEW & NOTEWORTHY An individualized precision neuroimaging approach reveals juxtaposed zones of the caudate that are coupled with five distinct networks in association cortex. The organization of these caudate zones recapitulates organization observed in the cerebral cortex and extends the notion of specialized basal ganglia circuits.

Normative modelling of brain morphometry across the lifespan with CentileBrain: algorithm benchmarking and model optimisation.

The value of normative models in research and clinical practice relies on their robustness and a systematic comparison of different modelling algorithms and parameters; however, this has not been done to date. We aimed to identify the optimal approach for normative modelling of brain morphometric data through systematic empirical benchmarking, by quantifying the accuracy of different algorithms and identifying parameters that optimised model performance. We developed this framework with regional morphometric data from 37 407 healthy individuals (53% female and 47% male; aged 3-90 years) from 87 datasets from Europe, Australia, the USA, South Africa, and east Asia following a comparative evaluation of eight algorithms and multiple covariate combinations pertaining to image acquisition and quality, parcellation software versions, global neuroimaging measures, and longitudinal stability. The multivariate fractional polynomial regression (MFPR) emerged as the preferred algorithm, optimised with non-linear polynomials for age and linear effects of global measures as covariates. The MFPR models showed excellent accuracy across the lifespan and within distinct age-bins and longitudinal stability over a 2-year period. The performance of all MFPR models plateaued at sample sizes exceeding 3000 study participants. This model can inform about the biological and behavioural implications of deviations from typical age-related neuroanatomical changes and support future study designs. The model and scripts described here are freely available through CentileBrain.

A hierarchical atlas of the human cerebellum for functional precision mapping.

The human cerebellum is activated by a wide variety of cognitive and motor tasks. Previous functional atlases have relied on single task-based or resting-state fMRI datasets. Here, we present a functional atlas that integrates information from 7 large-scale datasets, outperforming existing group atlasses. The new atlas has three further advantages: First, the atlas allows for precision mapping in individuals: The integration of the probabilistic group atlas with an individual localizer scan results in a marked improvement in prediction of individual boundaries. Second, we provide both asymmetric and symmetric versions of the atlas. The symmetric version, which is obtained by constraining the boundaries to be the same across hemispheres, is especially useful in studying functional lateralization. Finally, the regions are hierarchically organized across 3 levels, allowing analyses at the appropriate level of granularity. Overall, the new atlas is an important resource for the study of the interdigitated functional organization of the human cerebellum in health and disease.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers.

BACKGROUND: Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure. METHODS: Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference. RESULTS: For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness. CONCLUSIONS: We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Normative Modeling of Brain Morphometry Across the Lifespan Using CentileBrain: Algorithm Benchmarking and Model Optimization.

We present an empirically benchmarked framework for sex-specific normative modeling of brain morphometry that can inform about the biological and behavioral significance of deviations from typical age-related neuroanatomical changes and support future study designs. This framework was developed using regional morphometric data from 37,407 healthy individuals (53% female; aged 3-90 years) following a comparative evaluation of eight algorithms and multiple covariate combinations pertaining to image acquisition and quality, parcellation software versions, global neuroimaging measures, and longitudinal stability. The Multivariate Factorial Polynomial Regression (MFPR) emerged as the preferred algorithm optimized using nonlinear polynomials for age and linear effects of global measures as covariates. The MFPR models showed excellent accuracy across the lifespan and within distinct age-bins, and longitudinal stability over a 2-year period. The performance of all MFPR models plateaued at sample sizes exceeding 3,000 study participants. The model and scripts described here are freely available through CentileBrain (https://centilebrain.org/).

Side-by-side regions in dorsolateral prefrontal cortex estimated within the individual respond differentially to domain-specific and domain-flexible processes.

A recurring debate concerns whether regions of primate prefrontal cortex (PFC) support domain-flexible or domain-specific processes. Here we tested the hypothesis with functional MRI (fMRI) that side-by-side PFC regions, within distinct parallel association networks, differentially support domain-flexible and domain-specialized processing. Individuals (N = 9) were intensively sampled, and all effects were estimated within their own idiosyncratic anatomy. Within each individual, we identified PFC regions linked to distinct networks, including a dorsolateral PFC (DLPFC) region coupled to the medial temporal lobe (MTL) and an extended region associated with the canonical multiple-demand network. We further identified an inferior PFC region coupled to the language network. Exploration in separate task data, collected within the same individuals, revealed a robust functional triple dissociation. The DLPFC region linked to the MTL was recruited during remembering and imagining the future, distinct from juxtaposed regions that were modulated in a domain-flexible manner during working memory. The inferior PFC region linked to the language network was recruited during sentence processing. Detailed analysis of the trial-level responses further revealed that the DLPFC region linked to the MTL specifically tracked processes associated with scene construction. These results suggest that the DLPFC possesses a domain-specialized region that is small and easily confused with nearby (larger) regions associated with cognitive control. The newly described region is domain specialized for functions traditionally associated with the MTL. We discuss the implications of these findings in relation to convergent anatomical analysis in the monkey.NEW & NOTEWORTHY Competing hypotheses link regions of prefrontal cortex (PFC) to domain-flexible or domain-specific processes. Here, using a precision neuroimaging approach, we identify a domain-specialized region in dorsolateral PFC, coupled to the medial temporal lobe and recruited for scene construction. This region is juxtaposed to, but distinct from, broader PFC regions recruited flexibly for cognitive control. Region distinctions align with broader network differences, suggesting that PFC regions gain dissociable processing properties via segregated anatomical projections.

Human Striatal Association Megaclusters.

The striatum receives projections from multiple regions of the cerebral cortex consistent with its role in diverse motor, affective, and cognitive functions. Supporting cognitive functions, the caudate receives projections from cortical association regions. Building on recent insights about the details of how multiple cortical networks are specialized for distinct aspects of higher-order cognition, we revisited caudate organization using within-individual precision neuroimaging (n=2, each participant scanned 31 times). Detailed analysis revealed that the caudate has side-by-side zones that are coupled to at least Give distinct distributed association networks, paralleling the specialization observed in the cerebral cortex. Examining correlation maps from closely juxtaposed seed regions in the caudate recapitulated the Give distinct cerebral networks including their multiple spatially distributed regions. These results extend the general notion of parallel specialized basal ganglia circuits, with the additional discovery that even within the caudate, there is Gine-grained separation of multiple distinct higher-order networks.

Within-Individual Organization of the Human Cerebral Cortex: Networks, Global Topography, and Function.

The human cerebral cortex is populated by specialized regions that are organized into networks. Here we estimated networks using a Multi-Session Hierarchical Bayesian Model (MS-HBM) applied to intensively sampled within-individual functional MRI (fMRI) data. The network estimation procedure was initially developed and tested in two participants (each scanned 31 times) and then prospectively applied to 15 new participants (each scanned 8 to 11 times). Detailed analysis of the networks revealed a global organization. Locally organized first-order sensory and motor networks were surrounded by spatially adjacent second-order networks that also linked to distant regions. Third-order networks each possessed regions distributed widely throughout association cortex. Moreover, regions of distinct third-order networks displayed side-by-side juxtapositions with a pattern that repeated similarly across multiple cortical zones. We refer to these as Supra-Areal Association Megaclusters (SAAMs). Within each SAAM, two candidate control regions were typically adjacent to three separate domain-specialized regions. Independent task data were analyzed to explore functional response properties. The somatomotor and visual first-order networks responded to body movements and visual stimulation, respectively. A subset of the second-order networks responded to transients in an oddball detection task, consistent with a role in orienting to salient or novel events. The third-order networks, including distinct regions within each SAAM, showed two levels of functional specialization. Regions linked to candidate control networks responded to working memory load across multiple stimulus domains. The remaining regions within each SAAM did not track working memory load but rather dissociated across language, social, and spatial / episodic processing domains. These results support a model of the cerebral cortex in which progressively higher-order networks nest outwards from primary sensory and motor cortices. Within the apex zones of association cortex there is specialization of large-scale networks that divides domain-flexible from domain-specialized regions repeatedly across parietal, temporal, and prefrontal cortices. We discuss implications of these findings including how repeating organizational motifs may emerge during development.

Brain morphometry in older adults with and without dementia using extremely rapid structural scans.

T1-weighted structural MRI is widely used to measure brain morphometry (e.g., cortical thickness and subcortical volumes). Accelerated scans as fast as one minute or less are now available but it is unclear if they are adequate for quantitative morphometry. Here we compared the measurement properties of a widely adopted 1.0 mm resolution scan from the Alzheimer's Disease Neuroimaging Initiative (ADNI = 5'12'') with two variants of highly accelerated 1.0 mm scans (compressed-sensing, CSx6 = 1'12''; and wave-controlled aliasing in parallel imaging, WAVEx9 = 1'09'') in a test-retest study of 37 older adults aged 54 to 86 (including 19 individuals diagnosed with a neurodegenerative dementia). Rapid scans produced highly reliable morphometric measures that largely matched the quality of morphometrics derived from the ADNI scan. Regions of lower reliability and relative divergence between ADNI and rapid scan alternatives tended to occur in midline regions and regions with susceptibility-induced artifacts. Critically, the rapid scans yielded morphometric measures similar to the ADNI scan in regions of high atrophy. The results converge to suggest that, for many current uses, extremely rapid scans can replace longer scans. As a final test, we explored the possibility of a 0'49'' 1.2 mm CSx6 structural scan, which also showed promise. Rapid structural scans may benefit MRI studies by shortening the scan session and reducing cost, minimizing opportunity for movement, creating room for additional scan sequences, and allowing for the repetition of structural scans to increase precision of the estimates.

Functional specialization of parallel distributed networks revealed by analysis of trial-to-trial variation in processing demands.

Multiple large-scale networks populate human association cortex. Here, we explored the functional properties of these networks by exploiting trial-to-trial variation in component-processing demands. In two behavioral studies (n = 136 and n = 238), participants quantified strategies used to solve individual task trials that spanned remembering, imagining future scenarios, and various control trials. These trials were also all scanned in an independent sample of functional MRI participants (n = 10), each with sufficient data to precisely define within-individual networks. Stable latent factors varied across trials and correlated with trial-level functional responses selectively across networks. One network linked to parahippocampal cortex, labeled Default Network A (DN-A), tracked scene construction, including for control trials that possessed minimal episodic memory demands. To the degree, a trial encouraged participants to construct a mental scene with imagery and awareness about spatial locations of objects or places, the response in DN-A increased. The juxtaposed Default Network B (DN-B) showed no such response but varied in relation to social processing demands. Another adjacent network, labeled Frontoparietal Network B (FPN-B), robustly correlated with trial difficulty. These results support that DN-A and DN-B are specialized networks differentially supporting information processing within spatial and social domains. Both networks are dissociable from a closely juxtaposed domain-general control network that tracks cognitive effort.NEW & NOTEWORTHY Tasks shown to differentially recruit parallel association networks are multifaceted, leaving open questions about network processes. Here, examining trial-to-trial network response properties in relation to trial traits reveals new insights into network functions. In particular, processes linked to scene construction selectively recruit a distributed network with links to parahippocampal and retrosplenial cortices, including during trials designed not to rely on the personal past. Adjacent networks show distinct patterns, providing novel evidence of functional specialization.

Specialization of the Human Hippocampal Long Axis Revisited.

The hippocampus possesses anatomical differences along its long axis. Here the functional specialization of the human hippocampal long axis was explored using network-anchored precision functional MRI (N = 11) paired with behavioral analyses (N=266). Functional connectivity analyses demonstrated that the anterior hippocampus was preferentially correlated with a cerebral network associated with remembering, while the posterior hippocampus was correlated with a distinct network associated with behavioral salience. Seed regions placed within the hippocampus recapitulated the distinct cerebral networks. Functional characterization using task data within the same intensively sampled individuals discovered a functional double dissociation between the anterior and posterior hippocampal regions. The anterior hippocampal region was sensitive to remembering and imagining the future, specifically tracking the process of scene construction, while the posterior hippocampal region displayed transient responses to targets in an oddball detection task and to transitions between task blocks. These findings suggest specialization along the long axis of the hippocampus with differential responses reflecting the functional properties of the partner cerebral networks.

Within-Individual Organization of the Human Cognitive Cerebellum: Evidence for Closely Juxtaposed, Functionally Specialized Regions.

The human cerebellum possesses multiple regions linked to cerebral association cortex. Here we mapped the cerebellum using precision functional MRI within individual participants (N=15), first estimating regions using connectivity and then prospectively testing functional properties using independent task data. Network estimates in all participants revealed a Crus I / II cerebellar megacluster of five higher-order association networks often with multiple, discontinuous regions for the same network. Seed regions placed within the megaclusters, including the disjointed regions, yielded spatially selective networks in the cerebral cortex. Compelling evidence for functional specialization within the cerebellar megaclusters emerged from the task responses. Reflecting functional distinctions found in the cerebrum, domain-flexible cerebellar regions involved in cognitive control dissociated from distinct domain-specialized regions with differential responses to language, social, and spatial / episodic task demands. These findings provide a clear demonstration that the cerebellum encompasses multiple zones dedicated to cognition, featuring juxtaposed regions specialized for distinct processing domains.

Precision Brain Morphometry Using Cluster Scanning.

Measurement error limits the statistical power to detect group differences and longitudinal change in structural MRI morphometric measures (e.g., hippocampal volume, prefrontal thickness). Recent advances in scan acceleration enable extremely fast T1-weighted scans (~1 minute) to achieve morphometric errors that are close to the errors in longer traditional scans. As acceleration allows multiple scans to be acquired in rapid succession, it becomes possible to pool estimates to increase measurement precision, a strategy known as "cluster scanning." Here we explored brain morphometry using cluster scanning in a test-retest study of 40 individuals (12 younger adults, 18 cognitively unimpaired older adults, and 10 adults diagnosed with mild cognitive impairment or Alzheimer's Dementia). Morphometric errors from a single compressed sensing (CS) 1.0mm scan with 6x acceleration (CSx6) were, on average, 12% larger than a traditional scan using the Alzheimer's Disease Neuroimaging Initiative (ADNI) protocol. Pooled estimates from four clustered CSx6 acquisitions led to errors that were 34% smaller than ADNI despite having a shorter total acquisition time. Given a fixed amount of time, a gain in measurement precision can thus be achieved by acquiring multiple rapid scans instead of a single traditional scan. Errors were further reduced when estimates were pooled from eight CSx6 scans (51% smaller than ADNI). Neither pooling across a break nor pooling across multiple scan resolutions boosted this benefit. We discuss the potential of cluster scanning to improve morphometric precision, boost statistical power, and produce more sensitive disease progression biomarkers.

A third somatomotor representation in the human cerebellum.

Seminal neurophysiological studies in the 1940s discovered two somatomotor maps in the cerebellum-an inverted anterior lobe map and an upright posterior lobe map. Both maps have been confirmed in the human using noninvasive neuroimaging with additional hints of a third map within and near to the cerebellar vermis. Here, we sought direct evidence for the third somatomotor map by using intensive, repeated functional MRI (fMRI) scanning of individuals performing movements across multiple body parts (tongue, hands, glutes, and feet). An initial discovery sample (n = 4, 4 sessions per individual including 576 separate blocks of body movements) yielded evidence for the two established cerebellar somatomotor maps, as well as evidence for a third discontinuous foot representation within the vermis. When the left versus right foot movements were directly contrasted, the third representation could be clearly distinguished from the second representation in multiple individuals. Functional connectivity from seed regions in the third somatomotor representation confirmed anatomically specific connectivity with the cerebral cortex, paralleling the patterns observed for the two well-established maps. All results were prospectively replicated in an independent dataset with new individuals (n = 4). These collective findings provide direct support for a third somatomotor representation in the vermis of the cerebellum that may be part of a third map. We discuss the relations of this candidate third map to the broader topography of the cerebellum as well as its implications for understanding the specific organization of the human cerebellar vermis where distinct zones appear functionally specialized for somatomotor and visual domains.NEW & NOTEWORTHY A third somatomotor representation exists in the vermis of the human cerebellum. Evidence for this elusive representation arises specifically from mapping the foot. Separate foot representations distinguish the third from the nearby second somatomotor representation. A third somatomotor representation in the posterior vermis supports a large-scale organization hypothesis in which the cerebellum possesses three sets of roughly homotopic representations of the full cerebrum.

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

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.