Abstract:Cortical folding is a critical process during brain development, resulting in morphologies that are both consistent and distinct between individuals and species. While earlier studies have highlighted important aspects of cortical folding, most existing computational models, based on the differential growth theory, fall short of explaining why folds tend to appear in particular locations. The axon tension hypothesis may provide insight into this conundrum; however, there has been significant controversy about a potential role of axonal tension during the gyrification. The common opinion in the field is that axonal tension is inadequate to drive gyrification, but we currently run the risk of discarding this hypothesis without comprehensively studying the role of axonal tension. Here we propose a novel bi-layered finite element model incorporating the two theories, including characteristic axonal tension in the subcortex and differential cortical growth. We show that axon tension can serve as a perturbation sufficient to trigger buckling in simulations; similarly to other types of perturbations, the natural stability behavior of the system tends to determine some characteristics of the folding morphology (e.g. the wavelength) while the perturbation determines the location of folds. Certain geometries, however, can interact or compete with the natural stability of the system to change the wavelength. When multiple perturbations are present, they similarly compete with each other. We found that an axon bundle of reasonable size will overpower up to a 5\% thickness perturbation (typical in the literature) and determine fold placement. Finally, when multiple axon tracts are present, even a slight difference in axon stiffness, representing the heterogeneity of axonal connections, is enough to significantly change the folding pattern. While the simulations presented here are a very simple representation of white matter connectivity, our findings point to urgent future research on the role of axon connectivity in cortical folding.

Role of Axonal Fibers in the Cortical Folding Patterns: A Tale of Variability and Regularity

Role of mechanical factors in cortical folding development.

Mechanical Hierarchy in the Formation and Modulation of Cortical Folding Patterns.

Cortical Folding Pattern and its Consistency Induced by Biological Growth

Axonal tension contributes to consistent fold placement

Exploring the role of different cell types on cortical folding in the developing human brain through computational modeling

A Computational Model of Cerebral Cortex Folding

Relationship between Cortical Folding Pattern and Brain Network Characteristics

Joint Analysis of Cortical Folding Pattern, Thickness and Fiber Density Using Structural and DTI Data

The role of thickness inhomogeneities in hierarchical cortical folding

On early brain folding patterns using biomechanical growth modeling

Role of Data-driven Regional Growth Model in Shaping Brain Folding Patterns

The influence of biophysical parameters in a biomechanical model of cortical folding patterns

Axonal Fiber Terminations Concentrate on Gyri

Coevolution of Gyral Folding and Structural Connection Patterns in Primate Brains.

Evolutionarily-preserved Consistent Gyral Folding Patterns Across Primate Brains

Brain Cortical Functional Gradients Predict Cortical Folding Patterns via Attention Mesh Convolution

Global Perturbation of Initial Geometry in a Biomechanical Model of Cortical Morphogenesis

Numerical investigation of biomechanically coupled growth in cortical folding

Gene regulatory landscape of cerebral cortex folding

A mechanical method of cerebral cortical folding development based on thermal expansion