Commonly preserved and species-specific gyral folding patterns across primate brains

Learning lifespan brain anatomical correspondence via cortical developmental continuity transfer

Identifying anatomical correspondences in the human brain throughout the lifespan is an essential prerequisite for studying brain development and aging. But given the tremendous individual variability in cortical folding patterns, the heterogeneity of different neurodevelopmental stages, and the scarce of neuroimaging data, it is difficult to infer reliable lifespan anatomical correspondence at finer scales. To solve this problem, in this work, we take the advantage of the developmental continuity of the cerebral cortex and propose a novel transfer learning strategy: the model is trained from scratch using the age group with the largest sample size, and then is transferred and adapted to the other groups following the cortical developmental trajectory. A novel loss function is designed to ensure that during the transfer process the common patterns will be extracted and preserved, while the group-specific new patterns will be captured. The proposed framework was evaluated using multiple datasets covering four lifespan age groups with 1,000+ brains (from 34 gestational weeks to young adult). Our experimental results show that: 1) the proposed transfer strategy can dramatically improve the model performance on populations (e.g., early neurodevelopment) with very limited number of training samples; and 2) with the transfer learning we are able to robustly infer the complicated many-to-many anatomical correspondences among different brains at different neurodevelopmental stages. (Code will be released soon: https://github.com/qidianzl/CDC-transfer).

Modeling refined differences of cortical folding patterns via spatial, morphological, and temporal fusion representations

The gyrus, a pivotal cortical folding pattern, is essential for integrating brain structure-function. This study focuses on 2-Hinge and 3-Hinge folds, characterized by the gyral convergence from various directions. Existing voxel-level studies may not adequately capture the precise spatial relationships within cortical folding patterns, especially when relying solely on local cortical characteristics due to their variable shapes and homogeneous frequency-specific features. To overcome these challenges, we introduced a novel model that combines spatial distribution, morphological structure, and functional magnetic resonance imaging data. We utilized spatio-morphological residual representations to enhance and extract subtle variations in cortical spatial distribution and morphological structure during blood oxygenation, integrating these with functional magnetic resonance imaging embeddings using self-attention for spatio-morphological-temporal representations. Testing these representations for identifying cortical folding patterns, including sulci, gyri, 2-Hinge, and 2-Hinge folds, and evaluating the impact of phenotypic data (e.g. stimulus) on recognition, our experimental results demonstrate the model's superior performance, revealing significant differences in cortical folding patterns under various stimulus. These differences are also evident in the characteristics of sulci and gyri folds between genders, with 3-Hinge showing more variations. Our findings indicate that our representations of cortical folding patterns could serve as biomarkers for understanding brain structure-function correlations.

Mechanical hierarchy in the formation and modulation of cortical folding patterns

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.

An integrated finite element method and machine learning algorithm for brain morphology prediction

The human brain development experiences a complex evolving cortical folding from a smooth surface to a convoluted ensemble of folds. Computational modeling of brain development has played an essential role in better understanding the process of cortical folding, but still leaves many questions to be answered. A major challenge faced by computational models is how to create massive brain developmental simulations with affordable computational sources to complement neuroimaging data and provide reliable predictions for brain folding. In this study, we leveraged the power of machine learning in data augmentation and prediction to develop a machine-learning-based finite element surrogate model to expedite brain computational simulations, predict brain folding morphology, and explore the underlying folding mechanism. To do so, massive finite element method (FEM) mechanical models were run to simulate brain development using the predefined brain patch growth models with adjustable surface curvature. Then, a GAN-based machine learning model was trained and validated with these produced computational data to predict brain folding morphology given a predefined initial configuration. The results indicate that the machine learning models can predict the complex morphology of folding patterns, including 3-hinge gyral folds. The close agreement between the folding patterns observed in FEM results and those predicted by machine learning models validate the feasibility of the proposed approach, offering a promising avenue to predict the brain development with given fetal brain configurations.

Gyral peaks and patterns in human brains

Cortical folding patterns are related to brain function, cognition, and behavior. Since the relationship has not been fully explained on a coarse scale, many efforts have been devoted to the identification of finer grained cortical landmarks, such as sulcal pits and gyral peaks, which were found to remain invariant across subjects and ages and the invariance may be related to gene mediated proto-map. However, gyral peaks were only investigated on macaque monkey brains, but not on human brains where the investigation is challenged due to high inter-individual variabilities. To this end, in this work, we successfully identified 96 gyral peaks both on the left and right hemispheres of human brains, respectively. These peaks are spatially consistent across individuals. Higher or sharper peaks are more consistent across subjects. Both structural and functional graph metrics of peaks are significantly different from other cortical regions, and more importantly, these nodal graph metrics are anti-correlated with the spatial consistency metrics within peaks. In addition, the distribution of peaks and various cortical anatomical, structural/functional connective features show hemispheric symmetry. These findings provide new clues to understanding the cortical landmarks, as well as their relationship with brain functions, cognition, behavior in both healthy and aberrant brains.

Cortex2vector: anatomical embedding of cortical folding patterns

Current brain mapping methods highly depend on the regularity, or commonality, of anatomical structure, by forcing the same atlas to be matched to different brains. As a result, individualized structural information can be overlooked. Recently, we conceptualized a new type of cortical folding pattern called the 3-hinge gyrus (3HG), which is defined as the conjunction of gyri coming from three directions. Many studies have confirmed that 3HGs are not only widely existing on different brains, but also possess both common and individual patterns. In this work, we put further effort, based on the identified 3HGs, to establish the correspondences of individual 3HGs. We developed a learning-based embedding framework to encode individual cortical folding patterns into a group of anatomically meaningful embedding vectors (cortex2vector). Each 3HG can be represented as a combination of these embedding vectors via a set of individual specific combining coefficients. In this way, the regularity of folding pattern is encoded into the embedding vectors, while the individual variations are preserved by the multi-hop combination coefficients. Results show that the learned embeddings can simultaneously encode the commonality and individuality of cortical folding patterns, as well as robustly infer the complicated many-to-many anatomical correspondences among different brains.

Gyral peaks: Novel gyral landmarks in developing macaque brains

Cerebral cortex development undergoes a variety of processes, which provide valuable information for the study of the developmental mechanism of cortical folding as well as its relationship to brain structural architectures and brain functions. Despite the variability in the anatomy-function relationship on the higher-order cortex, recent studies have succeeded in identifying typical cortical landmarks, such as sulcal pits, that bestow specific functional and cognitive patterns and remain invariant across subjects and ages with their invariance being related to a gene-mediated proto-map. Inspired by the success of these studies, we aim in this study at defining and identifying novel cortical landmarks, termed gyral peaks, which are the local highest foci on gyri. By analyzing data from 156 MRI scans of 32 macaque monkeys with the age spanned from 0 to 36 months, we identified 39 and 37 gyral peaks on the left and right hemispheres, respectively. Our investigation suggests that these gyral peaks are spatially consistent across individuals and relatively stable within the age range of this dataset. Moreover, compared with other gyri, gyral peaks have a thicker cortex, higher mean curvature, more pronounced hub-like features in structural connective networks, and are closer to the borders of structural connectivity-based cortical parcellations. The spatial distribution of gyral peaks was shown to correlate with that of other cortical landmarks, including sulcal pits. These results provide insights into the spatial arrangement and temporal development of gyral peaks as well as their relation to brain structure and function.

Cortical development coupling between surface area and sulcal depth on macaque brains

Postnatal development of cerebral cortex is associated with a variety of neuronal processes and is thus critical to development of brain function and cognition. Longitudinal changes of cortical morphology and topology, such as postnatal cortical thinning and flattening have been widely studied. However, thorough and systematic investigation of such cortical change, including how to quantify it from multiple spatial directions and how to relate it to surface topology, is rarely found. In this work, based on a longitudinal macaque neuroimaging dataset, we quantified local changes in gyral white matter's surface area and sulcal depth during early development. We also investigated how these two metrics are coupled and how this coupling is linked to cortical surface topology, underlying white matter, and positions of functional areas. Semi-parametric generalized additive models were adopted to quantify the longitudinal changes of surface area (A) and sulcal depth (D), and the coupling patterns between them. This resulted in four classes of regions, according to how they change compared with global change throughout early development: slower surface area change and slower sulcal depth change (slowA_slowD), slower surface area change and faster sulcal depth change (slowA_fastD), faster surface area change and slower sulcal depth change (fastA_slowD), and faster surface area change and faster sulcal depth change (fastA_fastD). We found that cortex-related metrics, including folding pattern and cortical thickness, vary along slowA_fastD-fastA_slowD axis, and structural connection-related metrics vary along fastA_fastD-slowA_slowD axis, with which brain functional sites align better. It is also found that cortical landmarks, including sulcal pits and gyral hinges, spatially reside on the borders of the four patterns. These findings shed new lights on the relationship between cortex development, surface topology, axonal wiring pattern and brain functions.

NAS-optimized topology-preserving transfer learning for differentiating cortical folding patterns

Increasing evidence suggests that cortical folding patterns of human cerebral cortex manifest overt structural and functional differences. However, for interpretability, few studies leverage advanced techniques (e.g., deep learning) to investigate the difference among cortical folds, resulting in more differences yet to be extensively explored. To this end, we proposed an effective topology-preserving transfer learning framework to differentiate cortical fMRI time series extracted from cortical folds. Our framework consists of three main parts: (1) Neural architecture search (NAS), which is used to devise a well-performing network structure based on an initialized hand-designed super-graph in an image dataset; (2) Topology-preserving transfer, which takes the model searched by NAS as the source network, keeping the topological connectivity in the network unchanged, while transforming all 2D operations including convolution and pooling into 1D, therefore resulting in a topology-preserving network, named TPNAS-Net; (3) Classification and correlation analysis, which involves using the TPNAS-Net to classify 1D cortical fMRI time series for each individual brain, and performing a group difference analysis between autism spectrum disorder (ASD) and healthy control (HC) and correlation analysis with clinical information (i.e., age). Extensive experiments on two ASD datasets obtain consistent results, demonstrating that the TPNAS-Net not only discriminates cortical folding patterns at high classification accuracy, but also captures subtle differences between ASD and HC (p-value = 0.042). In addition, we discover that there is a positive correlation between the classification accuracy and age in ASD (r = 0.39, p-value = 0.04). These findings together suggest that structural and functional differences in cortical folding patterns between ASD and HC may provide a potentially useful biomarker for the diagnosis of ASD.

Gyral Hinges Account for the Highest Cost and the Highest Communication Capacity in a Corticocortical Network

Prior studies reported the global structure of brain networks exhibits the "small-world" and "rich-world" attributes. However, the underlying structural and functional architecture highlighted by these graph theory findings hasn't been explicitly related to the morphology of the cortex. This could be attributed to the lower resolution of used folding patterns, such as gyro-sulcal patterns. By defining a novel gyral folding pattern, termed gyral hinge (GH), which is the conjunction of ordinary gyri from multiple directions, we found GHs possess the highest length and cost in the white matter fiber connective network, and the shortest paths in the network tend to travel through GHs in their middle part. Based on these findings, we would hypothesize GHs could reside in the centers of a network core, thereby accounting for the highest cost and the highest communication capacity in a corticocortical network. The following results further support our hypothesis: 1) GHs possess stronger functional network integration capacity. 2) Higher cost is found on the connection with GHs to hinges and GHs to GHs. 3) Moving GHs introduces higher extra network cost. Our findings and hypotheses could reveal a profound relationship among the cortical folding patterns, axonal wiring architectures, and brain functions.

Mechanism Exploration of 3-Hinge Gyral Formation and Pattern Recognition

The 3-hinge gyral folding is the conjunction of gyrus crest lines from three different orientations. Previous studies have not explored the possible mechanisms of formation of such 3-hinge gyri, which are preserved across species in primate brains. We develop a biomechanical model to mimic the formation of 3-hinge patterns on a real brain and determine how special types of 3-hinge patterns form in certain areas of the model. Our computational and experimental imaging results show that most tertiary convolutions and exact locations of 3-hinge patterns after growth and folding are unpredictable, but they help explain the consistency of locations and patterns of certain 3-hinge patterns. Growing fibers within the white matter is posited as a determining factor to affect the location and shape of these 3-hinge patterns. Even if the growing fibers do not exert strong enough forces to guide gyrification directly, they still may seed a heterogeneous growth profile that leads to the formation of 3-hinge patterns in specific locations. A minor difference in initial morphology between two growing model brains can lead to distinct numbers and locations of 3-hinge patterns after folding.

Fundamental functional differences between gyri and sulci: implications for brain function, cognition, and behavior

Folding of the cerebral cortex is a prominent characteristic of mammalian brains. Alterations or deficits in cortical folding are strongly correlated with abnormal brain function, cognition, and behavior. Therefore, a precise mapping between the anatomy and function of the brain is critical to our understanding of the mechanisms of brain structural architecture in both health and diseases. Gyri and sulci, the standard nomenclature for cortical anatomy, serve as building blocks to make up complex folding patterns, providing a window to decipher cortical anatomy and its relation with brain functions. Huge efforts have been devoted to this research topic from a variety of disciplines including genetics, cell biology, anatomy, neuroimaging, and neurology, as well as involving computational approaches based on machine learning and artificial intelligence algorithms. However, despite increasing progress, our understanding of the functional anatomy of gyro-sulcal patterns is still in its infancy. In this review, we present the current state of this field and provide our perspectives of the methodologies and conclusions concerning functional differentiation between gyri and sulci, as well as the supporting information from genetic, cell biology, and brain structure research. In particular, we will further present a proposed framework for attempting to interpret the dynamic mechanisms of the functional interplay between gyri and sulci. Hopefully, this review will provide a comprehensive summary of anatomo-functional relationships in the cortical gyro-sulcal system together with a consideration of how these contribute to brain function, cognition, and behavior, as well as to mental disorders.

Anatomical correlates of face patches in macaque inferotemporal cortex

Primate brains typically have regions within the ventral visual stream that are selectively responsive to faces. In macaques, these face patches are located in similar parts of inferotemporal cortex across individuals although correspondence with particular anatomical features has not been reported previously. Here, using high-resolution functional and anatomical imaging, we show that small "bumps," or buried gyri, along the lower bank of the superior temporal sulcus are predictive of the location of face-selective regions. Recordings from implanted multielectrode arrays verified that these bumps contain face-selective neurons. These bumps were present in monkeys raised without seeing faces and that lack face patches, indicating that these anatomical landmarks are predictive of, but not sufficient for, the presence of face selectivity. These bumps are found across primate species that span taxonomy lines, indicating common evolutionary developmental mechanisms. The bumps emerge during fetal development in macaques, indicating that they arise from general developmental mechanisms that result in the regularity of cortical folding of the entire brain.

Cortical 3-hinges could serve as hubs in cortico-cortical connective network

Mapping the relation between cortical convolution and structural/functional brain architectures could provide deep insights into the mechanisms of brain development, evolution and diseases. In our previous studies, we found a unique gyral folding pattern, termed a 3-hinge, which was defined as the conjunction of three gyral crests. The uniqueness of the 3-hinge was evidenced by its thicker cortex and stronger fiber connections than other gyral regions. However, the role that 3-hinges play in cortico-cortical connective architecture remains unclear. To this end, we conducted MRI studies by constructing structural cortico-cortical connective networks based on a fine-granular cortical parcellation, the parcels of which were automatically labeled as 3-hinge, 2-hinge (ordinary gyrus) or sulcus. On human brains, 3-hinges possess significantly higher degrees, strengths and betweennesses than 2-hinges, suggesting that 3-hinges could serve more like hubs in the cortico-cortical connective network. This hypothesis gains supports from human functional network analyses, in which 3-hinges are involved in more global functional networks than ordinary gyri. In addition, 3-hinges could serve as 'connector' hubs rather than 'provincial' hubs and they account for a dominant proportion of nodes in the high-level 'backbone' of the network. These structural results are reproduced on chimpanzee and macaque brains, while the roles of 3-hinges as hubs become more pronounced in higher order primates. Our new findings could provide a new window to the relation between cortical convolution, anatomical connection and brain function.

Identifying Cross-individual Correspondences of 3-hinge Gyri

Human brain alignment based on imaging data has long been an intriguing research topic. One of the challenges is the huge inter-individual variabilities, which are pronounced not only in cortical folding patterns, but also in the underlying structural and functional patterns. Also, it is still not fully understood how to link the cross-subject similarity of cortical folding patterns to the correspondences of structural brain wiring diagrams and brain functions. Recently, a specific cortical gyral folding pattern was identified, which is the conjunction of gyri from multiple directions and termed a "gyral hinge". These gyral hinges are characterized by the thickest cortices, the densest long-range fibers, and the most complex functional profiles in contrast to other gyri. In addition to their structural and functional importance, a small portion of 3-hinges found correspondences across subjects and even species by manual labeling. However, it is unclear if such cross-subject correspondences can be found for all 3-hinges, or if the correspondences are interpretable from structural and functional aspects. Given the huge variability of cortical folding patterns, we proposed a novel algorithm which jointly uses structural MRI-derived cortical folding patterns and diffusion-MRI-derived fiber shape features to estimate the correspondences. This algorithm was executed in a group-wise manner, whereby 3-hinges of all subjects were simultaneously aligned. The effectiveness of the algorithm was demonstrated by higher cross-subject 3-hinges' consistency with respect to structural and functional metrics, when compared with other methods. Our findings provide a novel approach to brain alignment and an insight to the linkage between cortical folding patterns and the underlying structural connective diagrams and brain functions.

Joint representation of connectome-scale structural and functional profiles for identification of consistent cortical landmarks in macaque brain

Discovery and representation of common structural and functional cortical architectures has been a significant yet challenging problem for years. Due to the remarkable variability of structural and functional cortical architectures in human brain, it is challenging to jointly represent a common cortical architecture which can comprehensively encode both structure and function characteristics. In order to better understand this challenge and considering that macaque monkey brain has much less variability in structure and function compared with human brain, in this paper, we propose a novel computational framework to apply our DICCCOL (Dense Individualized and Common Connectivity-based Cortical Landmarks) and HAFNI (Holistic Atlases of Functional Networks and Interactions) frameworks on macaque brains, in order to jointly represent structural and functional connectome-scale profiles for identification of a set of consistent and common cortical landmarks across different macaque brains based on multimodal DTI and resting state fMRI (rsfMRI) data. Experimental results demonstrate that 100 consistent and common cortical landmarks are successfully identified via the proposed framework, each of which has reasonably accurate anatomical, structural fiber connection pattern, and functional correspondences across different macaque brains. This set of 100 landmarks offer novel insights into the structural and functional cortical architectures in macaque brains.

Denser Growing Fiber Connections Induce 3-hinge Gyral Folding

Recent studies have shown that quantitative description of gyral shape patterns offers a novel window to examine the relationship between brain structure and function. Along this research line, this paper examines a unique and interesting type of cortical gyral region where 3 different gyral crests meet, termed 3-hinge gyral region. We extracted 3-hinge gyral regions in macaque/chimpanzee/human brains, quantified and compared the relevant DTI-derived fiber densities in 3-hinge and 2-hinge gyral regions. Our observations consistently showed that DTI-derived fiber densities in 3-hinge regions are much higher than those in 2-hinge regions. Therefore, we hypothesize that besides the cortical expansion, denser fiber connections can induce the formation of 3-hinge gyri. To examine the biomechanical basis of this hypothesis, we constructed a series of 3-dimensional finite element soft tissue models based on continuum growth theory to investigate fundamental biomechanical mechanisms of consistent 3-hinge gyri formation. Our computational simulation results consistently showed that during gyrification gyral regions with higher concentrations of growing axonal fibers tend to form 3-hinge gyri. Our integrative approach combining neuroimaging data analysis and computational modeling appears effective in probing a plausible theory of 3-hinge gyri formation and providing new insights into structural and functional cortical architectures and their relationship.

Exploring 3-hinge gyral folding patterns among HCP Q3 868 human subjects

Comparison and integration of neuroimaging data from different brains and populations are fundamental in neuroscience. Over the past decades, the neuroimaging field has largely depended on image registration to compare and integrate neuroimaging data from individuals in a common reference space, with a basic assumption that the brains are similar. However, the intrinsic neuroanatomical complexity and huge interindividual cortical folding variation remain underexplored. Here we focus on a specific cortical convolution pattern, termed 3-hinge gyral folding, which is the conjunction of gyri from multiple orientations and has unique and consistent anatomically, structurally, and functionally connective patterns across subjects. By developing a novel shape descriptor and a two-stage clustering pipeline, we devise an automatic method to identify 3-hinges in the Human Connectome Project Q3 868 human brains, and further parameterize the complexity of such a pattern and quantify its regularity and variation in terms of 3-hinge number, position, and morphology. Our results not only exhibit the huge interindividual variations, but also reveal regular relationship between gyral hinges and other factors, such as their locations and cortical morphologies. It is found that "line-shape" cortices have relatively more consistent 3-hinge shape pattern distributions, and certain types of 3-hinge patterns favor particular cortical morphologies. In addition, more 3-hinges are found on "line-shape" cortices while their numbers vary more across subjects than those on "non-line-shape" cortices. This study adds new insights into a better understanding of the regularity and variability of human brain anatomy, and their functional aspects.

Gyral net: A new representation of cortical folding organization

One distinct feature of the cerebral cortex is its convex (gyri) and concave (sulci) folding patterns. Due to the remarkable complexity and variability of gyral/sulcal shapes, it has been challenging to quantitatively model their organization patterns. Inspired by the observation that the lines of gyral crests can form a connected graph on each brain hemisphere, we propose a new representation of cortical gyri/sulci organization pattern - gyral net, which models cortical architecture from a graph perspective, starting with nodes and edges obtained from the reconstructed cortical surfaces. A novel computational framework is developed to efficiently and automatically construct gyral nets from surface meshes, and four measurements are devised to quantify the folding patterns. Using an MRI dataset for autism study as a test bed, we identified reduced local connectivity cost and increased curviness of gyral net bilaterally on the parietal lobe, occipital lobe, and temporal lobe in autistic patients. Additionally, we found that the cortical thickness and the gyral straightness of gyral joints are higher than the rest of gyral crest regions. The proposed representation offers a new tool for a comprehensive and reliable characterization of the cortical folding organization. This novel computational framework will enable large-scale analyses of cortical folding patterns in the future.

Connectome-scale functional intrinsic connectivity networks in macaques

There have been extensive studies of intrinsic connectivity networks (ICNs) in the human brains using resting-state functional magnetic resonance imaging (fMRI) in the literature. However, the functional organization of ICNs in macaque brains has been less explored so far, despite growing interests in the field. In this work, we propose a computational framework to identify connectome-scale group-wise consistent ICNs in macaques via sparse representation of whole-brain resting-state fMRI data. Experimental results demonstrate that 70 group-wise consistent ICNs are successfully identified in macaque brains via the proposed framework. These 70 ICNs are interpreted based on two publicly available parcellation maps of macaque brains and our work significantly expand currently known macaque ICNs already reported in the literature. In general, this set of connectome-scale group-wise consistent ICNs can potentially benefit a variety of studies in the neuroscience and brain-mapping fields, and they provide a foundation to better understand brain evolution in the future.

Radial Structure Scaffolds Convolution Patterns of Developing Cerebral Cortex

Commonly-preserved radial convolution is a prominent characteristic of the mammalian cerebral cortex. Endeavors from multiple disciplines have been devoted for decades to explore the causes for this enigmatic structure. However, the underlying mechanisms that lead to consistent cortical convolution patterns still remain poorly understood. In this work, inspired by prior studies, we propose and evaluate a plausible theory that radial convolution during the early development of the brain is sculptured by radial structures consisting of radial glial cells (RGCs) and maturing axons. Specifically, the regionally heterogeneous development and distribution of RGCs controlled by Trnp1 regulate the convex and concave convolution patterns (gyri and sulci) in the radial direction, while the interplay of RGCs' effects on convolution and axons regulates the convex (gyral) convolution patterns. This theory is assessed by observations and measurements in literature from multiple disciplines such as neurobiology, genetics, biomechanics, etc., at multiple scales to date. Particularly, this theory is further validated by multimodal imaging data analysis and computational simulations in this study. We offer a versatile and descriptive study model that can provide reasonable explanations of observations, experiments, and simulations of the characteristic mammalian cortical folding.