skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Compensating cortical thickness for cortical folding-related variation
Cortical thickness is a widely used biomarker of brain morphology and health, yet it is dependent on local cortical folding. Because gyral crowns are consistently thicker than sulcal fundi and cortical folds vary widely across individuals, these fluctuations introduce unmodeled nuisance variance that can obscure meaningful biological effects of interest. Previous global methods of folding compensation incompletely compensate for folding effects on cortical thickness. Spatial smoothing is commonly used to reduce these effects in the literature, but this markedly degrades spatial localization precision. To address these limitations, we developed a novel method for folding-compensated cortical thickness estimation that uses nonlinear local multiple regression with five folding measures to model and more completely remove folding-related variance from cortical thickness. This approach estimates what cortical thickness would have been in the absence of folding, yielding a more biologically interpretable measure of cortical architecture. We applied this new approach to data from the Young Adult Human Connectome Project (HCP-YA) and Aging Human Connectome Project (HCA), demonstrating substantial reductions in intra-areal and inter-individual variability, substantially increasing standardized effect sizes of age on cortical thickness (41% increase) while preserving neurobiologically expected patterns, and avoiding the loss of spatial precision that occurs with the spatial smoothing that has traditionally been used in the literature. The method has been integrated into the HCP pipelines, facilitating its widespread use. By attenuating folding-induced variability, this technique enhances cortical thickness as a structural phenotype and may support more accurate cortical parcellation, longitudinal tracking, and biomarker discovery in brain health and disease.  more » « less
Award ID(s):
2144412
PAR ID:
10664044
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
The MIT Press
Date Published:
Journal Name:
Imaging Neuroscience
Volume:
4
ISSN:
2837-6056
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract The cortical thickness is a characteristic biomarker for a wide variety of neurological disorders. While the structural organization of the cerebral cortex is tightly regulated and evolutionarily preserved, its thickness varies widely between 1.5 and 4.5 mm across the healthy adult human brain. It remains unclear whether these thickness variations are a cause or consequence of cortical development. Recent studies suggest that cortical thickness variations are primarily a result of genetic effects. Previous studies showed that a simple homogeneous bilayered system with a growing layer on an elastic substrate undergoes a unique symmetry breaking into a spatially heterogeneous system with discrete gyri and sulci. Here, we expand on that work to explore the evolution of cortical thickness variations over time to support our finding that cortical pattern formation and thickness variations can be explained – at least in part – by the physical forces that emerge during cortical folding. Strikingly, as growth progresses, the developing gyri universally thicken and the sulci thin, even in the complete absence of regional information. Using magnetic resonance images, we demonstrate that these naturally emerging thickness variations agree with the cortical folding pattern in n = 9 healthy adult human brains, in n = 564 healthy human brains ages 7–64, and in n = 73 infant brains scanned at birth, and at ages one and two. Additionally, we show that cortical organoids develop similar patterns throughout their growth. Our results suggest that genetic, geometric, and physical events during brain development are closely interrelated. Understanding regional and temporal variations in cortical thickness can provide insight into the evolution and causative factors of neurological disorders, inform the diagnosis of neurological conditions, and assess the efficacy of treatment options. 
    more » « less
  2. Across mammalia, brain morphology follows specific scaling patterns. Bigger bodies have bigger brains, with surface area outpacing volume growth, resulting in increased foldedness. We have recently studied scaling rules of cortical thickness, both local and global, finding that the cortical thickness difference between thick gyri and thin sulci also increases with brain size and foldedness. Here, we investigate early brain development in humans, using subjects from the Developing Human Connectome Project, scanned shortly after pre-term or full-term birth, yielding magnetic resonance images of the brain from 29 to 43 postmenstrual weeks. While the global cortical thickness does not change significantly during this development period, its distribution does, with sulci thinning, while gyri thickening. By comparing our results with our recent work on humans and 11 non-human primate species, we also compare the trajectories of primate evolution with human development, noticing that the 2 trends are distinct for volume, surface area, cortical thickness, and gyrification index. Finally, we introduce the global shape index as a proxy for gyrification index; while correlating very strongly with gyrification index, it offers the advantage of being calculated only from local quantities without generating a convex hull or alpha surface. 
    more » « less
  3. Abstract The relationship of human brain structure to cognitive function is complex, and how this relationship differs between childhood and adulthood is poorly understood. One strong hypothesis suggests the cognitive function of Fluid Intelligence (Gf) is dependent on prefrontal cortex and parietal cortex. In this work, we developed a novel graph convolutional neural networks (gCNNs) for the analysis of localized anatomic shape and prediction of Gf. Morphologic information of the cortical ribbons and subcortical structures was extracted from T1-weighted MRIs within two independent cohorts, the Adolescent Brain Cognitive Development Study (ABCD; age: 9.93 ± 0.62 years) of children and the Human Connectome Project (HCP; age: 28.81 ± 3.70 years). Prediction combining cortical and subcortical surfaces together yielded the highest accuracy of Gf for both ABCD (R = 0.314) and HCP datasets (R = 0.454), outperforming the state-of-the-art prediction of Gf from any other brain measures in the literature. Across both datasets, the morphology of the amygdala, hippocampus, and nucleus accumbens, along with temporal, parietal and cingulate cortex consistently drove the prediction of Gf, suggesting a significant reframing of the relationship between brain morphology and Gf to include systems involved with reward/aversion processing, judgment and decision-making, motivation, and emotion. 
    more » « less
  4. null (Ed.)
    Human cortex is patterned by a complex and interdigitated web of large-scale functional networks. Recent methodological breakthroughs reveal variation in the size, shape, and spatial topography of cortical networks across individuals. While spatial network organization emerges across development, is stable over time, and is predictive of behavior, it is not yet clear to what extent genetic factors underlie interindividual differences in network topography. Here, leveraging a nonlinear multidimensional estimation of heritability, we provide evidence that individual variability in the size and topographic organization of cortical networks are under genetic control. Using twin and family data from the Human Connectome Project ( n = 1,023), we find increased variability and reduced heritability in the size of heteromodal association networks ( h 2 : M = 0.34, SD = 0.070), relative to unimodal sensory/motor cortex ( h 2 : M = 0.40, SD = 0.097). We then demonstrate that the spatial layout of cortical networks is influenced by genetics, using our multidimensional estimation of heritability ( h 2 - multi; M = 0.14, SD = 0.015). However, topographic heritability did not differ between heteromodal and unimodal networks. Genetic factors had a regionally variable influence on brain organization, such that the heritability of network topography was greatest in prefrontal, precuneus, and posterior parietal cortex. Taken together, these data are consistent with relaxed genetic control of association cortices relative to primary sensory/motor regions and have implications for understanding population-level variability in brain functioning, guiding both individualized prediction and the interpretation of analyses that integrate genetics and neuroimaging. 
    more » « less
  5. Abstract The important mechanical parameters and their hierarchy in the growth and folding of the human brain have not been thoroughly understood. In this study, we developed a multiscale mechanical model to investigate how the interplay between initial geometrical undulations, differential tangential growth in the cortical plate, and axonal connectivity form and regulate the folding patterns of the human brain in a hierarchical order. To do so, different growth scenarios with bilayer spherical models that features initial undulations on the cortex and uniform or heterogeneous distribution of axonal fibers in the white matter were developed, statistically analyzed, and validated by the imaging observations. The results showed that the differential tangential growth is the inducer of cortical folding, and in a hierarchal order, high-amplitude initial undulations on the surface and axonal fibers in the substrate regulate the folding patterns and determine the location of gyri and sulci. The locations with dense axonal fibers after folding settle in gyri rather than sulci. The statistical results also indicated that there is a strong correlation between the location of positive (outward) and negative (inward) initial undulations and the locations of gyri and sulci after folding, respectively. In addition, the locations of 3-hinge gyral folds are strongly correlated with the initial positive undulations and locations of dense axonal fibers. As another finding, it was revealed that there is a correlation between the density of axonal fibers and local gyrification index, which has been observed in imaging studies but not yet fundamentally explained. This study is the first step in understanding the linkage between abnormal gyrification (surface morphology) and disruption in connectivity that has been observed in some brain disorders such as Autism Spectrum Disorder. Moreover, the findings of the study directly contribute to the concept of the regularity and variability of folding patterns in individual human brains. 
    more » « less