Feature similarity gradients detect alterations in the neonatal cortex associated with preterm birth.

The early life environment programmes cortical architecture and cognition across the life course. A measure of cortical organisation that integrates information from multimodal MRI and is unbound by arbitrary parcellations has proven elusive, which hampers efforts to uncover the perinatal origins of cortical health. Here, we use the Vogt-Bailey index to provide a fine-grained description of regional homogeneities and sharp variations in cortical microstructure based on feature gradients, and we investigate the impact of being born preterm on cortical development at term-equivalent age. Compared with term-born controls, preterm infants have a homogeneous microstructure in temporal and occipital lobes, and the medial parietal, cingulate, and frontal cortices, compared with term infants. These observations replicated across two independent datasets and were robust to differences that remain in the data after matching samples and alignment of processing and quality control strategies. We conclude that cortical microstructural architecture is altered in preterm infants in a spatially distributed rather than localised fashion.

BOUNTI: Brain vOlumetry and aUtomated parcellatioN for 3D feTal MRI.

Fetal MRI is widely used for quantitative brain volumetry studies. However, currently, there is a lack of universally accepted protocols for fetal brain parcellation and segmentation. Published clinical studies tend to use different segmentation approaches that also reportedly require significant amounts of time-consuming manual refinement. In this work, we propose to address this challenge by developing a new robust deep learning-based fetal brain segmentation pipeline for 3D T2w motion corrected brain images. At first, we defined a new refined brain tissue parcellation protocol with 19 regions-of-interest using the new fetal brain MRI atlas from the Developing Human Connectome Project. This protocol design was based on evidence from histological brain atlases, clear visibility of the structures in individual subject 3D T2w images and the clinical relevance to quantitative studies. It was then used as a basis for developing an automated deep learning brain tissue parcellation pipeline trained on 360 fetal MRI datasets with different acquisition parameters using semi-supervised approach with manually refined labels propagated from the atlas. The pipeline demonstrated robust performance for different acquisition protocols and GA ranges. Analysis of tissue volumetry for 390 normal participants (21-38 weeks gestational age range), scanned with three different acquisition protocols, did not reveal significant differences for major structures in the growth charts. Only minor errors were present in < 15% of cases thus significantly reducing the need for manual refinement. In addition, quantitative comparison between 65 fetuses with ventriculomegaly and 60 normal control cases were in agreement with the findings reported in our earlier work based on manual segmentations. These preliminary results support the feasibility of the proposed atlas-based deep learning approach for large-scale volumetric analysis. The created fetal brain volumetry centiles and a docker with the proposed pipeline are publicly available online at https://hub.docker.com/r/fetalsvrtk/segmentation (tag brain_bounti_tissue).

Structural and functional asymmetry of the neonatal cerebral cortex.

Features of brain asymmetry have been implicated in a broad range of cognitive processes; however, their origins are still poorly understood. Here we investigated cortical asymmetries in 442 healthy term-born neonates using structural and functional magnetic resonance images from the Developing Human Connectome Project. Our results demonstrate that the neonatal cortex is markedly asymmetric in both structure and function. Cortical asymmetries observed in the term cohort were contextualized in two ways: by comparing them against cortical asymmetries observed in 103 preterm neonates scanned at term-equivalent age, and by comparing structural asymmetries against those observed in 1,110 healthy young adults from the Human Connectome Project. While associations with preterm birth and biological sex were minimal, significant differences exist between birth and adulthood.

The Multiscale Surface Vision Transformer.

Surface meshes are a favoured domain for representing structural and functional information on the human cortex, but their complex topology and geometry pose significant challenges for deep learning analysis. While Transformers have excelled as domain-agnostic architectures for sequence-to-sequence learning, notably for structures where the translation of the convolution operation is non-trivial, the quadratic cost of the self-attention operation remains an obstacle for many dense prediction tasks. Inspired by some of the latest advances in hierarchical modelling with vision transformers, we introduce the Multiscale Surface Vision Transformer (MS-SiT) as a backbone architecture for surface deep learning. The self-attention mechanism is applied within local-mesh-windows to allow for high-resolution sampling of the underlying data, while a shifted-window strategy improves the sharing of information between windows. Neighbouring patches are successively merged, allowing the MS-SiT to learn hierarchical representations suitable for any prediction task. Results demonstrate that the MS-SiT outperforms existing surface deep learning methods for neonatal phenotyping prediction tasks using the Developing Human Connectome Project (dHCP) dataset. Furthermore, building the MS-SiT backbone into a U-shaped architecture for surface segmentation demonstrates competitive results on cortical parcellation using the UK Biobank (UKB) and manually-annotated MindBoggle datasets. Code and trained models are publicly available at https://github.com/metrics-lab/surface-vision-transformers.

Telomere length and brain imaging phenotypes in UK Biobank.

Telomeres form protective caps at the ends of chromosomes, and their attrition is a marker of biological aging. Short telomeres are associated with an increased risk of neurological and psychiatric disorders including dementia. The mechanism underlying this risk is unclear, and may involve brain structure and function. However, the relationship between telomere length and neuroimaging markers is poorly characterized. Here we show that leucocyte telomere length (LTL) is associated with multi-modal MRI phenotypes in 31,661 UK Biobank participants. Longer LTL is associated with: i) larger global and subcortical grey matter volumes including the hippocampus, ii) lower T1-weighted grey-white tissue contrast in sensory cortices, iii) white-matter microstructure measures in corpus callosum and association fibres, iv) lower volume of white matter hyperintensities, and v) lower basal ganglia iron. Longer LTL was protective against certain related clinical manifestations, namely all-cause dementia (HR 0.93, 95% CI: 0.91-0.96), but not stroke or Parkinson's disease. LTL is associated with multiple MRI endophenotypes of neurodegenerative disease, suggesting a pathway by which longer LTL may confer protective against dementia.

ICAM-Reg: Interpretable Classification and Regression With Feature Attribution for Mapping Neurological Phenotypes in Individual Scans.

An important goal of medical imaging is to be able to precisely detect patterns of disease specific to individual scans; however, this is challenged in brain imaging by the degree of heterogeneity of shape and appearance. Traditional methods, based on image registration, historically fail to detect variable features of disease, as they utilise population-based analyses, suited primarily to studying group-average effects. In this paper we therefore take advantage of recent developments in generative deep learning to develop a method for simultaneous classification, or regression, and feature attribution (FA). Specifically, we explore the use of a VAE-GAN (variational autoencoder - general adversarial network) for translation called ICAM, to explicitly disentangle class relevant features, from background confounds, for improved interpretability and regression of neurological phenotypes. We validate our method on the tasks of Mini-Mental State Examination (MMSE) cognitive test score prediction for the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort, as well as brain age prediction, for both neurodevelopment and neurodegeneration, using the developing Human Connectome Project (dHCP) and UK Biobank datasets. We show that the generated FA maps can be used to explain outlier predictions and demonstrate that the inclusion of a regression module improves the disentanglement of the latent space. Our code is freely available on GitHub https://github.com/CherBass/ICAM.

CortexODE: Learning Cortical Surface Reconstruction by Neural ODEs.

We present CortexODE, a deep learning framework for cortical surface reconstruction. CortexODE leverages neural ordinary differential equations (ODEs) to deform an input surface into a target shape by learning a diffeomorphic flow. The trajectories of the points on the surface are modeled as ODEs, where the derivatives of their coordinates are parameterized via a learnable Lipschitz-continuous deformation network. This provides theoretical guarantees for the prevention of self-intersections. CortexODE can be integrated to an automatic learning-based pipeline, which reconstructs cortical surfaces efficiently in less than 5 seconds. The pipeline utilizes a 3D U-Net to predict a white matter segmentation from brain Magnetic Resonance Imaging (MRI) scans, and further generates a signed distance function that represents an initial surface. Fast topology correction is introduced to guarantee homeomorphism to a sphere. Following the isosurface extraction step, two CortexODE models are trained to deform the initial surface to white matter and pial surfaces respectively. The proposed pipeline is evaluated on large-scale neuroimage datasets in various age groups including neonates (25-45 weeks), young adults (22-36 years) and elderly subjects (55-90 years). Our experiments demonstrate that the CortexODE-based pipeline can achieve less than 0.2mm average geometric error while being orders of magnitude faster compared to conventional processing pipelines.

A bipolar taxonomy of adult human brain sulcal morphology related to timing of fetal sulcation and trans-sulcal gene expression gradients.

We developed a computational pipeline (now provided as a resource) for measuring morphological similarity between cortical surface sulci to construct a sulcal phenotype network (SPN) from each magnetic resonance imaging (MRI) scan in an adult cohort (N=34,725; 45-82 years). Networks estimated from pairwise similarities of 40 sulci on 5 morphological metrics comprised two clusters of sulci, represented also by the bipolar distribution of sulci on a linear-to-complex dimension. Linear sulci were more heritable and typically located in unimodal cortex; complex sulci were less heritable and typically located in heteromodal cortex. Aligning these results with an independent fetal brain MRI cohort (N=228; 21-36 gestational weeks), we found that linear sulci formed earlier, and the earliest and latest-forming sulci had the least between-adult variation. Using high-resolution maps of cortical gene expression, we found that linear sulcation is mechanistically underpinned by trans-sulcal gene expression gradients enriched for developmental processes.

Neonatal multi-modal cortical profiles predict 18-month developmental outcomes.

Developmental delays in infanthood often persist, turning into life-long difficulties, and coming at great cost for the individual and community. By examining the developing brain and its relation to developmental outcomes we can start to elucidate how the emergence of brain circuits is manifested in variability of infant motor, cognitive and behavioural capacities. In this study, we examined if cortical structural covariance at birth, indexing coordinated development, is related to later infant behaviour. We included 193 healthy term-born infants from the Developing Human Connectome Project (dHCP). An individual cortical connectivity matrix derived from morphological and microstructural features was computed for each subject (morphometric similarity networks, MSNs) and was used as input for the prediction of behavioural scores at 18 months using Connectome-Based Predictive Modeling (CPM). Neonatal MSNs successfully predicted social-emotional performance. Predictive edges were distributed between and within known functional cortical divisions with a specific important role for primary and posterior cortical regions. These results reveal that multi-modal neonatal cortical profiles showing coordinated maturation are related to developmental outcomes and that network organization at birth provides an early infrastructure for future functional skills.

Preterm birth alters the development of cortical microstructure and morphology at term-equivalent age.

INTRODUCTION: The dynamic nature and complexity of the cellular events that take place during the last trimester of pregnancy make the developing cortex particularly vulnerable to perturbations. Abrupt interruption to normal gestation can lead to significant deviations to many of these processes, resulting in atypical trajectory of cortical maturation in preterm birth survivors. METHODS: We sought to first map typical cortical micro- and macrostructure development using invivo MRI in a large sample of healthy term-born infants scanned after birth (n = 259). Then we offer a comprehensive characterization of the cortical consequences of preterm birth in 76 preterm infants scanned at term-equivalent age (37-44 weeks postmenstrual age). We describe the group-average atypicality, the heterogeneity across individual preterm infants, and relate individual deviations from normative development to age at birth and neurodevelopment at 18 months. RESULTS: In the term-born neonatal brain, we observed heterogeneous and regionally specific associations between age at scan and measures of cortical morphology and microstructure, including rapid surface expansion, greater cortical thickness, lower cortical anisotropy and higher neurite orientation dispersion. By term-equivalent age, preterm infants had on average increased cortical tissue water content and reduced neurite density index in the posterior parts of the cortex, and greater cortical thickness anteriorly compared to term-born infants. While individual preterm infants were more likely to show extreme deviations (over 3.1 standard deviations) from normative cortical maturation compared to term-born infants, these extreme deviations were highly variable and showed very little spatial overlap between individuals. Measures of regional cortical development were associated with age at birth, but not with neurodevelopment at 18 months. CONCLUSION: We showed that preterm birth alters cortical micro- and macrostructural maturation near the time of full-term birth. Deviations from normative development were highly variable between individual preterm infants.

A survey on active learning and human-in-the-loop deep learning for medical image analysis.

Fully automatic deep learning has become the state-of-the-art technique for many tasks including image acquisition, analysis and interpretation, and for the extraction of clinically useful information for computer-aided detection, diagnosis, treatment planning, intervention and therapy. However, the unique challenges posed by medical image analysis suggest that retaining a human end-user in any deep learning enabled system will be beneficial. In this review we investigate the role that humans might play in the development and deployment of deep learning enabled diagnostic applications and focus on techniques that will retain a significant input from a human end user. Human-in-the-Loop computing is an area that we see as increasingly important in future research due to the safety-critical nature of working in the medical domain. We evaluate four key areas that we consider vital for deep learning in the clinical practice: (1) Active Learning to choose the best data to annotate for optimal model performance; (2) Interaction with model outputs - using iterative feedback to steer models to optima for a given prediction and offering meaningful ways to interpret and respond to predictions; (3) Practical considerations - developing full scale applications and the key considerations that need to be made before deployment; (4) Future Prospective and Unanswered Questions - knowledge gaps and related research fields that will benefit human-in-the-loop computing as they evolve. We offer our opinions on the most promising directions of research and how various aspects of each area might be unified towards common goals.

Phenotyping the Preterm Brain: Characterizing Individual Deviations From Normative Volumetric Development in Two Large Infant Cohorts.

The diverse cerebral consequences of preterm birth create significant challenges for understanding pathogenesis or predicting later outcome. Instead of focusing on describing effects common to the group, comparing individual infants against robust normative data offers a powerful alternative to study brain maturation. Here we used Gaussian process regression to create normative curves characterizing brain volumetric development in 274 term-born infants, modeling for age at scan and sex. We then compared 89 preterm infants scanned at term-equivalent age with these normative charts, relating individual deviations from typical volumetric development to perinatal risk factors and later neurocognitive scores. To test generalizability, we used a second independent dataset comprising of 253 preterm infants scanned using different acquisition parameters and scanner. We describe rapid, nonuniform brain growth during the neonatal period. In both preterm cohorts, cerebral atypicalities were widespread, often multiple, and varied highly between individuals. Deviations from normative development were associated with respiratory support, nutrition, birth weight, and later neurocognition, demonstrating their clinical relevance. Group-level understanding of the preterm brain disguises a large degree of individual differences. We provide a method and normative dataset that offer a more precise characterization of the cerebral consequences of preterm birth by profiling the individual neonatal brain.

Cortical morphology at birth reflects spatiotemporal patterns of gene expression in the fetal human brain.

Interruption to gestation through preterm birth can significantly impact cortical development and have long-lasting adverse effects on neurodevelopmental outcome. We compared cortical morphology captured by high-resolution, multimodal magnetic resonance imaging (MRI) in n = 292 healthy newborn infants (mean age at birth = 39.9 weeks) with regional patterns of gene expression in the fetal cortex across gestation (n = 156 samples from 16 brains, aged 12 to 37 postconceptional weeks [pcw]). We tested the hypothesis that noninvasive measures of cortical structure at birth mirror areal differences in cortical gene expression across gestation, and in a cohort of n = 64 preterm infants (mean age at birth = 32.0 weeks), we tested whether cortical alterations observed after preterm birth were associated with altered gene expression in specific developmental cell populations. Neonatal cortical structure was aligned to differential patterns of cell-specific gene expression in the fetal cortex. Principal component analysis (PCA) of 6 measures of cortical morphology and microstructure showed that cortical regions were ordered along a principal axis, with primary cortex clearly separated from heteromodal cortex. This axis was correlated with estimated tissue maturity, indexed by differential expression of genes expressed by progenitor cells and neurons, and engaged in stem cell differentiation, neuron migration, and forebrain development. Preterm birth was associated with altered regional MRI metrics and patterns of differential gene expression in glial cell populations. The spatial patterning of gene expression in the developing cortex was thus mirrored by regional variation in cortical morphology and microstructure at term, and this was disrupted by preterm birth. This work provides a framework to link molecular mechanisms to noninvasive measures of cortical development in early life and highlights novel pathways to injury in neonatal populations at increased risk of neurodevelopmental disorder.

The developing Human Connectome Project (dHCP) automated resting-state functional processing framework for newborn infants.

The developing Human Connectome Project (dHCP) aims to create a detailed 4-dimensional connectome of early life spanning 20-45 weeks post-menstrual age. This is being achieved through the acquisition of multi-modal MRI data from over 1000 in- and ex-utero subjects combined with the development of optimised pre-processing pipelines. In this paper we present an automated and robust pipeline to minimally pre-process highly confounded neonatal resting-state fMRI data, robustly, with low failure rates and high quality-assurance. The pipeline has been designed to specifically address the challenges that neonatal data presents including low and variable contrast and high levels of head motion. We provide a detailed description and evaluation of the pipeline which includes integrated slice-to-volume motion correction and dynamic susceptibility distortion correction, a robust multimodal registration approach, bespoke ICA-based denoising, and an automated QC framework. We assess these components on a large cohort of dHCP subjects and demonstrate that processing refinements integrated into the pipeline provide substantial reduction in movement related distortions, resulting in significant improvements in SNR, and detection of high quality RSNs from neonates.

Development of Microstructural and Morphological Cortical Profiles in the Neonatal Brain.

Interruptions to neurodevelopment during the perinatal period may have long-lasting consequences. However, to be able to investigate deviations in the foundation of proper connectivity and functional circuits, we need a measure of how this architecture evolves in the typically developing brain. To this end, in a cohort of 241 term-born infants, we used magnetic resonance imaging to estimate cortical profiles based on morphometry and microstructure over the perinatal period (37-44 weeks postmenstrual age, PMA). Using the covariance of these profiles as a measure of inter-areal network similarity (morphometric similarity networks; MSN), we clustered these networks into distinct modules. The resulting modules were consistent and symmetric, and corresponded to known functional distinctions, including sensory-motor, limbic, and association regions, and were spatially mapped onto known cytoarchitectonic tissue classes. Posterior regions became more morphometrically similar with increasing age, while peri-cingulate and medial temporal regions became more dissimilar. Network strength was associated with age: Within-network similarity increased over age suggesting emerging network distinction. These changes in cortical network architecture over an 8-week period are consistent with, and likely underpin, the highly dynamic processes occurring during this critical period. The resulting cortical profiles might provide normative reference to investigate atypical early brain development.

Cross-species cortical alignment identifies different types of anatomical reorganization in the primate temporal lobe.

Evolutionary adaptations of temporo-parietal cortex are considered to be a critical specialization of the human brain. Cortical adaptations, however, can affect different aspects of brain architecture, including local expansion of the cortical sheet or changes in connectivity between cortical areas. We distinguish different types of changes in brain architecture using a computational neuroanatomy approach. We investigate the extent to which between-species alignment, based on cortical myelin, can predict changes in connectivity patterns across macaque, chimpanzee, and human. We show that expansion and relocation of brain areas can predict terminations of several white matter tracts in temporo-parietal cortex, including the middle and superior longitudinal fasciculus, but not the arcuate fasciculus. This demonstrates that the arcuate fasciculus underwent additional evolutionary modifications affecting the temporal lobe connectivity pattern. This approach can flexibly be extended to include other features of cortical organization and other species, allowing direct tests of comparative hypotheses of brain organization.

How did language evolve? Since the human lineage diverged from that of the other great apes millions of years ago, changes in the brain have given rise to behaviors that are unique to humans, such as language. Some of these changes involved alterations in the size and relative positions of brain areas, while others required changes in the connections between those regions. But did these changes occur independently, or can the changes observed in one actually explain the changes we see in the other? One way to answer this question is to use neuroimaging to compare the brains of related species, using different techniques to examine different aspects of brain structure. Imaging a fatty substance called myelin, for example, can produce maps showing the size and position of brain areas. Measuring how easily water molecules diffuse through brain tissue, by contrast, provides information about connections between areas. Eichert et al. performed both types of imaging in macaques and healthy human volunteers, and compared the results to existing data from chimpanzees. Computer simulations were used to manipulate the myelin-based images so that equivalent brain areas in each species occupied the same positions. In most cases, the distortions ­ or 'warping' ­ needed to superimpose brain regions on top of one another also predicted the differences between species in the connections between those regions. This suggests that movement of brain regions over the course of evolution explain the differences previously observed in brain connectivity. But there was one notable exception, namely a bundle of fibers with a key role in language called the arcuate fasciculus. This structure follows a slightly different route through the brain in humans compared to chimpanzees and macaques. Eichert et al. show that this difference cannot be explained solely by changes in the positions of brain regions. Instead, the arcuate fasciculus underwent additional changes in its course, which may have contributed to the evolution of language. The framework developed by Eichert et al. can be used to study evolution in many different species. Interspecies comparisons can provide clues to how brain structure and activity relate to each other and to behavior, and this knowledge could ultimately help to understand and treat brain disorders.

Modelling brain development to detect white matter injury in term and preterm born neonates.

Premature birth occurs during a period of rapid brain growth. In this context, interpreting clinical neuroimaging can be complicated by the typical changes in brain contrast, size and gyrification occurring in the background to any pathology. To model and describe this evolving background in brain shape and contrast, we used a Bayesian regression technique, Gaussian process regression, adapted to multiple correlated outputs. Using MRI, we simultaneously estimated brain tissue intensity on T1- and T2-weighted scans as well as local tissue shape in a large cohort of 408 neonates scanned cross-sectionally across the perinatal period. The resulting model provided a continuous estimate of brain shape and intensity, appropriate to age at scan, degree of prematurity and sex. Next, we investigated the clinical utility of this model to detect focal white matter injury. In individual neonates, we calculated deviations of a neonate's observed MRI from that predicted by the model to detect punctate white matter lesions with very good accuracy (area under the curve > 0.95). To investigate longitudinal consistency of the model, we calculated model deviations in 46 neonates who were scanned on a second occasion. These infants' voxelwise deviations from the model could be used to identify them from the other 408 images in 83% (T2-weighted) and 76% (T1-weighted) of cases, indicating an anatomical fingerprint. Our approach provides accurate estimates of non-linear changes in brain tissue intensity and shape with clear potential for radiological use.

Extending the Human Connectome Project across ages: Imaging protocols for the Lifespan Development and Aging projects.

The Human Connectome Projects in Development (HCP-D) and Aging (HCP-A) are two large-scale brain imaging studies that will extend the recently completed HCP Young-Adult (HCP-YA) project to nearly the full lifespan, collecting structural, resting-state fMRI, task-fMRI, diffusion, and perfusion MRI in participants from 5 to 100+ years of age. HCP-D is enrolling 1300+ healthy children, adolescents, and young adults (ages 5-21), and HCP-A is enrolling 1200+ healthy adults (ages 36-100+), with each study collecting longitudinal data in a subset of individuals at particular age ranges. The imaging protocols of the HCP-D and HCP-A studies are very similar, differing primarily in the selection of different task-fMRI paradigms. We strove to harmonize the imaging protocol to the greatest extent feasible with the completed HCP-YA (1200+ participants, aged 22-35), but some imaging-related changes were motivated or necessitated by hardware changes, the need to reduce the total amount of scanning per participant, and/or the additional challenges of working with young and elderly populations. Here, we provide an overview of the common HCP-D/A imaging protocol including data and rationales for protocol decisions and changes relative to HCP-YA. The result will be a large, rich, multi-modal, and freely available set of consistently acquired data for use by the scientific community to investigate and define normative developmental and aging related changes in the healthy human brain.

Construction of a neonatal cortical surface atlas using Multimodal Surface Matching in the Developing Human Connectome Project.

We propose a method for constructing a spatio-temporal cortical surface atlas of neonatal brains aged between 36 and 44 weeks of post-menstrual age (PMA) at the time of scan. The data were acquired as part of the Developing Human Connectome Project (dHCP), and the constructed surface atlases are publicly available. The method is based on a spherical registration approach: Multimodal Surface Matching (MSM), using cortical folding for driving the alignment. Templates have been generated for the anatomical cortical surface and for the cortical feature maps: sulcal depth, curvature, thickness, T1w/T2w myelin maps and cortical regions. To achieve this, cortical surfaces from 270 infants were first projected onto the sphere. Templates were then generated in two stages: first, a reference space was initialised via affine alignment to a group average adult template. Following this, templates were iteratively refined through repeated alignment of individuals to the template space until the variability of the average feature sets converged. Finally, bias towards the adult reference was removed by applying the inverse of the average affine transformations on the template and de-drifting the template. We used temporal adaptive kernel regression to produce age-dependant atlases for 9 weeks (36-44 weeks PMA). The generated templates capture expected patterns of cortical development including an increase in gyrification as well as an increase in thickness and T1w/T2w myelination with increasing age.

Dynamic patterns of cortical expansion during folding of the preterm human brain.

During the third trimester of human brain development, the cerebral cortex undergoes dramatic surface expansion and folding. Physical models suggest that relatively rapid growth of the cortical gray matter helps drive this folding, and structural data suggest that growth may vary in both space (by region on the cortical surface) and time. In this study, we propose a unique method to estimate local growth from sequential cortical reconstructions. Using anatomically constrained multimodal surface matching (aMSM), we obtain accurate, physically guided point correspondence between younger and older cortical reconstructions of the same individual. From each pair of surfaces, we calculate continuous, smooth maps of cortical expansion with unprecedented precision. By considering 30 preterm infants scanned two to four times during the period of rapid cortical expansion (28-38 wk postmenstrual age), we observe significant regional differences in growth across the cortical surface that are consistent with the emergence of new folds. Furthermore, these growth patterns shift over the course of development, with noninjured subjects following a highly consistent trajectory. This information provides a detailed picture of dynamic changes in cortical growth, connecting what is known about patterns of development at the microscopic (cellular) and macroscopic (folding) scales. Since our method provides specific growth maps for individual brains, we are also able to detect alterations due to injury. This fully automated surface analysis, based on tools freely available to the brain-mapping community, may also serve as a useful approach for future studies of abnormal growth due to genetic disorders, injury, or other environmental variables.