A new method to measure cortical growth in the developing brain

Folding of the cerebral cortex is a critical phase of brain development in higher mammals but the biomechanics of folding remain incompletely understood. During folding, the growth of the cortical surface is heterogeneous and anisotropic. We developed and applied a new technique to measure spatial and directional variations in surface growth from longitudinal magnetic resonance imaging (MRI) studies of a single animal or human subject. MRI provides high resolution 3D image volumes of the brain at different stages of development. Surface representations of the cerebral cortex are obtained by segmentation of these volumes. Estimation of local surface growth between two times requires establishment of a point-to-point correspondence ("registration") between surfaces measured at those times. Here we present a novel approach for the registration of two surfaces in which an energy function is minimized by solving a partial differential equation on a spherical surface. The energy function includes a strain-energy term due to distortion and an "error energy" term due to mismatch between surface features. This algorithm, implemented with the finite element method, brings surface features into approximate alignment while minimizing deformation in regions without explicit matching criteria. The method was validated by application to three simulated test cases and applied to characterize growth of the ferret cortex during folding. Cortical surfaces were created from MRI data acquired in vivo at 14 days, 21 days, and 28 days of life. Deformation gradient and Lagrangian strain tensors describe the kinematics of growth over this interval. These quantitative results illuminate the spatial, temporal, and directional patterns of growth during cortical folding.     

Cerebral cortex expansion and folding: what have we learned?

One of the most prominent features of the human brain is the fabulous size of the cerebral cortex and its intricate folding. Cortical folding takes place during embryonic development and is important to optimize the functional organization and wiring of the brain, as well as to allow fitting a large cortex in a limited cranial volume. Pathological alterations in size or folding of the human cortex lead to severe intellectual disability and intractable epilepsy. Hence, cortical expansion and folding are viewed as key processes in mammalian brain development and evolution, ultimately leading to increased intellectual performance and, eventually, to the emergence of human cognition. Here, we provide an overview and discuss some of the most significant advances in our understanding of cortical expansion and folding over the last decades. These include discoveries in multiple and diverse disciplines, from cellular and molecular mechanisms regulating cortical development and neurogenesis, genetic mechanisms defining the patterns of cortical folds, the biomechanics of cortical growth and buckling, lessons from human disease, and how genetic evolution steered cortical size and folding during mammalian evolution .     

Mechanical role of a growing solid tumor on cortical folding

Cortical folding, or convolution of the brain, is a vital process in mammals that causes the brain to have a wrinkled appearance. The existence of different types of prenatal solid tumors may alter this complex phenomenon and cause severe brain disorders. Here we interpret the effects of a growing solid tumor on the cortical folding in the fetal brain by virtue of theoretical analyses and computational modeling. The developing fetal brain is modeled as a simple, double-layered, and soft structure with an outer cortex and an inner core, in combination with a circular tumor model imbedded in the structure to investigate the developmental mechanism of cortical convolution. Analytical approaches offer introductory insight into the deformation field and stress distribution of a developing brain. After the onset of instability, analytical approaches fail to capture complex secondary evolution patterns, therefore a series of non-linear finite element simulations are carried out to study the crease formation and the influence from a growing solid tumor inside the structure. Parametric studies show the dependency of the cortical folding pattern on the size, location, and growth speed of a solid tumor in fetal brain. It is noteworthy to mention that there is a critical distance from the cortex/core interface where the growing tumor shows its pronounced effect on the cortical convolution, and that a growing tumor decreases the gyrification index of cortical convolution while its stiffness does not have a profound effect on the gyrification process.     

Emerging roles of neural stem cells in cerebral cortex development and evolution

Expansion and folding of the cerebral cortex are landmark features of mammalian brain evolution, which are recapitulated during embryonic development. Neural stem cells and their derived germinal cells are coordinated during cerebral cortex development to produce the appropriate amounts and types of neurons. This process is further complicated in gyrencephalic species, where newborn neurons must disperse in the tangential axis to expand the cerebral cortex in surface area. Here, we review advances that have been made over the last decade in understanding the nature and diversity of telencephalic neural stem cells and their roles in cortical development, and we discuss recent progress on how newly identified types of cortical progenitor cell populations may have evolved to drive the expansion and folding of the mammalian cerebral cortex.     

综述: 大脑新皮层和神经发育疾病

哺乳动物进化过程中,大脑皮层逐渐增大增厚和脑容量增大,从而构成了脑神经环路复杂性的细胞生物学基础.皮层出现皱褶是非人类灵长类演化的重要特征.成体人脑大约由近860多亿个神经细胞组成,其中,在人脑神经发生高峰,每小时有近400多万个兴奋性神经细胞产生.如此高速的神经生成过程需要精确的细胞与分子调控机制.本文主要讨论调控大脑皮层增大增厚的细胞与分子机制和相关的脑发育疾病.     

A computational model of hemodynamic parameters in cortical capillary networks

The analysis of hemodynamic parameters and functional reactivity of cerebral capillaries is still controversial. To assess the hemodynamic parameters in the cortical capillary network, a generic model was created using 2D voronoi tessellation in which each edge represents a capillary segment. This method is capable of creating an appropriate generic model of cerebral capillary network relating to each part of the brain cortex because the geometric model is able to vary the capillary density. The modeling presented here is based on morphometric parameters extracted from physiological data of the human cortex. The pertinent hemodynamic parameters were obtained by numerical simulation based on effective blood viscosity as a function of hematocrit and microvessel diameter, phase separation and plasma skimming effects. The hemodynamic parameters of capillary networks with two different densities (consistent with the variation of the morphometric data in the human cortical capillary network) were analyzed. The results show pertinent hemodynamic parameters for each model. The heterogeneity (coefficient variation) and the mean value of hematocrits, flow rates and velocities of the both network models were specified. The distributions of blood flow throughout the both models seem to confirm the hypothesis in which all capillaries in a cortical network are recruited at rest (normal condition). The results also demonstrate a discrepancy of the network resistance between two models, which are derived from the difference in the number density of capillary segments between the models.     

Dopamine D1-stimulated adenylyl cyclase activity in cerebral cortex of autopsied human brain.

Although the cerebral cortical dopamine D(1) receptor is considered to play a role in normal and abnormal brain function, little information is available on its characteristics in human brain. We compared dopamine-stimulated adenylyl cyclase (AC) activity in homogenates of cerebral cortex (frontal, temporal, parietal, occipital and cingulate cortex) of autopsied brain of neurologically normal subjects to that in striatum. Cerebral cortical AC activity was modestly and dose-dependently stimulated by dopamine (maximal 20-30%) with low microM EC50s and such stimulation was inhibited by the selective dopamine D1 receptor antagonist SCH23390. The magnitude of the maximal stimulation by dopamine was similar in autopsied and biopsied cerebral cortex. The extent of maximal stimulation was similar to that in dopamine-rich striatum (caudate, putamen and nucleus accumbens), despite much lower density of dopamine D1 receptors in cerebral cortex vs. striatum. The EC50 for dopamine stimulation in cerebral cortex (approximately 1 microM) was lower than that for caudate and putamen (approximately 3 microM). No detectable dopamine stimulation was observed in cerebellar cortex, thalamus or hippocampus. Dopamine stimulation in both cerebral cortex and striatum was independent of calcium activation. We conclude that dopamine stimulated AC can be measured in cerebral cortex of human brain allowing for the possibility that this process can be examined in human brain disorders in which dopaminergic abnormalities are suspected.     

Change in the Cortical Complexity of Spinocerebellar Ataxia Type 3 Appears Earlier than Clinical Symptoms

Patients with spinocerebellar ataxia type 3 (SCA3) have exhibited cerebral cortical involvement and various mental deficits in previous studies. Clinically, conventional measurements, such as the Mini-Mental State Examination (MMSE) and electroencephalography (EEG), are insensitive to cerebral cortical involvement and mental deficits associated with SCA3, particularly at the early stage of the disease. We applied a three-dimensional fractal dimension (3D-FD) method, which can be used to quantify the shape complexity of cortical folding, in assessing cortical degeneration. We evaluated 48 genetically confirmed SCA3 patients by employing clinical scales and magnetic resonance imaging and using 50 healthy participants as a control group. According to the Scale for the Assessment and Rating of Ataxia (SARA), the SCA3 patients were diagnosed with cortical dysfunction in the cerebellar cortex; however, no significant difference in the cerebral cortex was observed according to the patients’ MMSE ratings. Using the 3D-FD method, we determined that cortical involvement was more extensive than involvement of traditional olivopontocerebellar regions and the corticocerebellar system. Moreover, the significant correlation between decreased 3D-FD values and disease duration may indicate atrophy of the cerebellar cortex and cerebral cortex in SCA3 patients. The change of the cerebral complexity in the SCA3 patients can be detected throughout the disease duration, especially it becomes substantial at the late stage of the disease. Furthermore, we determined that atrophy of the cerebral cortex may occur earlier than changes in MMSE scores and EEG signals.     

On the Possible Shapes of the Brain

The human brain is unique among primates in its complexity and variability. Here I argue that this variability is, however, strongly constrained by developmental processes common to all mammals. Comparative analyses of grey and white matter volume, cortical surface area and cortical folding show that the rostro?Ccaudal axis of the central nervous system is a main direction along which mammalian neuroanatomical diversity is organised. Phylogenetically, rostral structures are often disproportionately larger and more differentiated in large mammals compared with small ones. Ontogenetically, caudal structures differentiate earlier but show less variation among species than rostral structures, which differentiate later and for a longer period. Theoretical considerations suggest that growth oriented along the rostro?Ccaudal axis should produce non-linear differences in white matter volume and cortical folding. Growth appears then as a fundamental parameter to understand mammalian neuroanatomical variability, whose effects should be common to all species. This seems to be indeed the case for humans: the volume of different brain structures as well as changes in the extension and folding of the cerebral cortex resemble the trends observed across mammals. A strong global pattern of coordinated variability emerges, where differences in total brain volume are non-linearly related to local neuroanatomical changes. Finally, I review evidence suggesting that the changes related to this global pattern of variability may have an influence on the organisation of behaviour, modulating the development of certain cognitive traits or even affecting the susceptibility to psychiatric disorders.     

Organizational levels of the cerebral cortex: An integrated model

We propose a theoretical model of the cerebral cortex which is based on its cellular components and integrates its different levels of organization: (1) cells have general adaptive and memorization properties; (2) cortical columns are repetitive interneuronal circuits which determine an adaptive processing specific to the cerebral cortex; (3) cortical maps effect selective combinations which are very efficient to learn basic behaviourial adaptations such as invariant recognition of forms, visually-guided hand movements, or execution of structured motor programs; (4) the network between cortical areas has a global architecture which integrates successive learning experiences into coherent functions such as the human language.     

Role of mechanical factors in the morphology of the primate cerebral cortex

The convoluted cortex of primates is instantly recognizable in its principal morphologic features, yet puzzling in its complex finer structure. Various hypotheses have been proposed about the mechanisms of its formation. Based on the analysis of databases of quantitative architectonic and connection data for primate prefrontal cortices, we offer support for the hypothesis that tension exerted by corticocortical connections is a significant factor in shaping the cerebral cortical landscape. Moreover, forces generated by cortical folding influence laminar morphology, and appear to have a previously unsuspected impact on cellular migration during cortical development. The evidence for a significant role of mechanical factors in cortical morphology opens the possibility of constructing computational models of cortical develoment based on physical principles. Such models are particularly relevant for understanding the relationship of cortical morphology to the connectivity of normal brains, and structurally altered brains in diseases of developmental origin, such as schizophrenia and autism.     

Differential involvement of cell cycle reactivation between striatal and cortical neurons in cell death induced by 3-nitropropionic acid

Recent evidence suggests that unscheduled cell cycle activity leads to neuronal cell death. 3-Nitropropionic acid (3-NP) is an irreversible inhibitor of succinate dehydrogenase and induces cell death in both striatum and cerebral cortex. Here we analyzed the involvement of aberrant cell cycle progression in 3-NP-induced cell death in these brain regions. 3-NP reduced the level of cyclin-dependent kinase inhibitor p27 in striatum but not in cerebral cortex. 3-NP also induced phosphorylation of retinoblastoma protein, a marker of cell cycle progression at late G(1) phase, only in striatum. Pharmacological experiments revealed that cyclin-dependent kinase activity and N-methyl-d-aspartate (NMDA) receptor were cooperatively involved in cell death by 3-NP in striatal neurons, whereas only NMDA receptor was involved in 3-NP-induced neurotoxicity in cortical neurons. Death of striatal neurons was preceded by elevation of somatic Ca(2+) and activation of calpain, a Ca(2+)-dependent protease. Both striatal p27 down-regulation and cell death provoked by 3-NP were dependent on calpain activity. Moreover, transfection of p27 small interfering RNA reduced striatal cell viability. In cortical neurons, however, there was no change in somatic Ca(2+) and calpain activity by 3-NP, and calpain inhibitors were not protective. These results suggest that 3-NP induces aberrant cell cycle progression and neuronal cell death via p27 down-regulation by calpain in striatum but not in the cerebral cortex. This is the first report for differential involvement of cell cycle reactivation in different brain regions and lightens the mechanism for region-selective vulnerability in human disease, including Huntington disease.     

3D model-based approach to identification of laminar structures of the cerebral cortex: Application to Brodmann areas 17 and 18

The human cerebral cortex may be subdivided into architectonic fields according to variations within its laminar structure. Studies have shown correspondences between the locations of functional activation foci and architectonic regions. In order to perform accurate localization of functional activation foci to architectonic regions, a parcellation algorithm capable of segmenting architectonic regions on in vivo imaging datasets is required. This paper presents a novel 3D model-based approach to directly detect cortical layers and classify architectonic fields. The column-like structure of the cortex is modeled using a Laplace equation method which generates a collection of intensity profiles that span the cortical mantle. Bayesian evidence for intensity profile elements belonging to hyper- or hypo-intense bands, which represent cell or myelin poor or rich layers in imaging data, is gathered. A non-isotropic Markov Random Field model is used to encourage contiguous bands as well as a penalty term that completes bands across highly curved cortical regions where neighbouring evidence for banding is strong. This algorithm is validated on a 3D histological dataset of a macaque brain with visible layering at intermediate resolution between high-resolution MRI and histology. The algorithm detects the myelin-rich Stria of Gennari and uses this as the basis for finding the Brodmann Area 17/18 boundary.     

Specific modification of heparan sulphate is required for normal cerebral cortical development

Proteoglycans are cell surface and extracellular matrix molecules to which long, unbranched glycosaminoglycan side chains are attached. Heparan sulphate, a type of glycosaminoglycan chain, has been proposed as a co-factor necessary for signalling by a range of growth factors. Here we provide evidence that loss of 2-O-sulphation in heparan sulphate leads to a significant reduction in cell proliferation in the developing cerebral cortex. The gene encoding heparan sulphate 2-sulphotransferase (Hs2st) is expressed in embryonic cortex and histological analysis of mice homozygous for a null mutation in Hs2st indicated a reduction in the thickness of the embryonic cerebral cortex. Using 5′-bromodeoxyuridine (BrdU) incorporation assays we found a reduction of approximately 40% in labelling indices of cortical precursor cells at E12. Comparison of the fates of cortical cells born on E13 and E15 in Hs2st^−/− mutant and wildtype littermate embryos revealed no differences in the pattern of cell migration. Our findings suggest a critical role for 2-O-sulphation of heparan sulphate proteoglycan (HSPG) in regulating cell proliferation during development of the cerebral cortex, perhaps through the modulation of cellular responses to growth factor signalling.     

A geometric model for the cortical folding pattern of simple folded brains

A model for the folding of the cerebral cortex is described which effects a direct link between the geometry of the brain and its folding pattern. Using some simple approximations to the brain's geometry the model is able to predict folding patterns. As the model is based on an assumption of uniformity of the cortex, the success of its application tends to lend support to homogeneous mechanical models of folding.     

Dopamine D4-like binding sites labeled by [3H]nemonapride include substantial serotonin 5-HT2A receptors in primate cerebral cortex

Dopamine D4-like binding sites are abundant in human cerebral cortex as detected by [3H]nemonapride. The extremely low density of D4 mRNA in human cerebral cortex is inconsistent with the high amount of D4-like binding sites. To investigate the nature of the D4-like receptors, [3H]nemonapride binding sites in the nonhuman primate cerebral cortex were characterized. Although [3H]nemonapride binding sites were D4-like, displaceable by clozapine but not raclopride, [3H]nemonapride binding was not displaced by selective D4 antagonists but was displaced by the selective 5-HT2A antagonist MDL100907. Using [3H]ketanserin as a 5-HT2A ligand, nemonapride showed high affinity for monkey (Ki = 10.4 nM) and cloned human (Ki = 9.4 nM) 5-HT2A receptors, while its affinity for rat receptors was lower (Ki = 140 nM). The present study demonstrates that cerebral cortical D4-like binding sites labeled by [3H]nemonapride in nonhuman primates consist of a very small portion of D4, but a substantial portion of 5-HT2A receptors. The unexpectedly high affinity of nemonapride for primate 5-HT2A receptor suggests reconsidering previous data from other studies using [3H]nemonapride, particularly those on D4-like receptors.     

Brain structure and spatial sensitivity profile assessing by near‐infrared spectroscopy modeling based on 3D MRI data

The goal of this study is to prove that the light propagation in the head by used the 3‐D optical model from in vivo MRI data set can also provide significant characteristics on the spatial sensitivity of cerebral cortex folding geometry based on Monte Carlo simulation. Thus, we proposed a MRI based approach for 3‐D brain modeling of near‐infrared spectroscopy (NIRS). In the results, the spatial sensitivity profile of the cerebral cortex folding geometry and the arrangement of source‐detector separation have being necessarily considered for applications of functional NIRS. The optimal choice of source‐detector separation is suggested within 3–3.5 cm by the received intensity with different source‐detector separations and the ratio of received light from the gray and white matter layer is greater than 50%. Additionally, this study has demonstrated the capability of NIRS in not only assessing the functional but also detecting the structural change of the brain by taking advantage of the low scattering and absorption coefficients observed in CSF of sagittal view. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Quantifying Agreement between Anatomical and Functional Interhemispheric Correspondences in the Resting Brain

The human brain is composed of two broadly symmetric cerebral hemispheres, with an abundance of reciprocal anatomical connections between homotopic locations. However, to date, studies of hemispheric symmetries have not identified correspondency precisely due to variable cortical folding patterns. Here we present a method to establish accurate correspondency using position on the unfolded cortical surface relative to gyral and sulcal landmarks. The landmark method is shown to outperform the method of reversing standard volume coordinates, and it is used to quantify the functional symmetry in resting fMRI data throughout the cortex. Resting brain activity was found to be maximally correlated with locations less than 1 cm away on the cortical surface from the corresponding anatomical location in nearly half of the cortex. While select locations exhibited asymmetric patterns, precise symmetric relationships were found to be the norm, with fine-grained symmetric functional maps demonstrated in motor, occipital, and inferior frontal cortex.     

Orchestrated neuronal migration and cortical folding: A computational and experimental study

Brain development involves precisely orchestrated genetic, biochemical, and mechanical events. At the cellular level, neuronal proliferation in the innermost zone of the brain followed by migration towards the outermost layer results in a rapid increase in brain surface area, outpacing the volumetric growth of the brain, and forming the highly folded cortex. This work aims to provide mechanistic insights into the process of brain development and cortical folding using a biomechanical model that couples cell division and migration with volumetric growth. Unlike phenomenological growth models, our model tracks the spatio-temporal development of cohorts of neurons born at different times, with each cohort modeled separately as an advection-diffusion process and the total cell density determining the extent of volume growth. We numerically implement our model in Abaqus/Standard (2020) by writing user-defined element (UEL) subroutines. For model calibration, we apply in utero electroporation (IUE) to ferret brains to visualize and track cohorts of neurons born at different stages of embryonic development. Our calibrated simulations of cortical folding align qualitatively with the ferret experiments. We have made our experimental data and finite-element implementation available online to offer other researchers a modeling platform for future study of neurological disorders associated with atypical neurodevelopment and cortical malformations.     

Regulation of cerebral cortex size and folding by expansion of basal progenitors

Size and folding of the cerebral cortex increased massively during mammalian evolution leading to the current diversity of brain morphologies. Various subtypes of neural stem and progenitor cells have been proposed to contribute differently in regulating thickness or folding of the cerebral cortex during development, but their specific roles have not been demonstrated. We report that the controlled expansion of unipotent basal progenitors in mouse embryos led to megalencephaly, with increased surface area of the cerebral cortex, but not to cortical folding. In contrast, expansion of multipotent basal progenitors in the naturally gyrencephalic ferret was sufficient to drive the formation of additional folds and fissures. In both models, changes occurred while preserving a structurally normal, six‐layered cortex. Our results are the first experimental demonstration of specific and distinct roles for basal progenitor subtypes in regulating cerebral cortex size and folding during development underlying the superior intellectual capability acquired by higher mammals during evolution.