Adolescent Brain Maturation and Cortical Folding: Evidence for Reductions in Gyrification (2014)

PLoS One. 2014; 9(1): e84914.

Published online Jan 15, 2014. doi:  10.1371/journal.pone.0084914
PMCID: PMC3893168
Maurice Ptito, Editor


Evidence from anatomical and functional imaging studies have highlighted major modifications of cortical circuits during adolescence. These include reductions of gray matter (GM), increases in the myelination of cortico-cortical connections and changes in the architecture of large-scale cortical networks. It is currently unclear, however, how the ongoing developmental processes impact upon the folding of the cerebral cortex and how changes in gyrification relate to maturation of GM/WM-volume, thickness and surface area. In the current study, we acquired high-resolution (3 Tesla) magnetic resonance imaging (MRI) data from 79 healthy subjects (34 males and 45 females) between the ages of 12 and 23 years and performed whole brain analysis of cortical folding patterns with the gyrification index (GI). In addition to GI-values, we obtained estimates of cortical thickness, surface area, GM and white matter (WM) volume which permitted correlations with changes in gyrification. Our data show pronounced and widespread reductions in GI-values during adolescence in several cortical regions which include precentral, temporal and frontal areas. Decreases in gyrification overlap only partially with changes in the thickness, volume and surface of GM and were characterized overall by a linear developmental trajectory. Our data suggest that the observed reductions in GI-values represent an additional, important modification of the cerebral cortex during late brain maturation which may be related to cognitive development.


A large body of work during last two decades has highlighted the importance of adolescence for the continued maturation of cortical circuits [1][3]. Starting with the observation of Huttenlocher [4] of marked decreases in the number of synaptic contacts, magnetic resonance imaging (MRI) studies have disclosed pronounced reductions in the volume and thickness of gray matter (GM) [5], [6]. In contrast, the amount of white matter (WM) has been shown to increase as a result of improved myelination of cortico-cortical connections [7][10]. More recent research has indicated that modifications in GM/WM extend into the third decade of life [11], [12] and involve changes in the large-scale organization of anatomical and functional networks [13]. These findings provided novel insights into the importance of adolescence as a critical period of human brain development which may also hold important clues for the emergence of psychiatric disorders, such as schizophrenia, which typically manifest during the transition from adolescence to adulthood [14], [15].

While modifications in the volume of GM/WM have been extensively characterized, relatively little evidence exists on maturational changes in the folding of the cortical surface. The cerebral cortex in humans has as one of its distinguishing characteristics a highly convoluted folding pattern that leads to a significantly increased cortical surface. For example, the surface area of the human cortex is per average ten times larger than that of the macaques monkey but only twice as thick [16]. The increased cortical surface in humans may be related to the emergence of higher cognitive functions because of the large number of neurons and cortico-cortical connections that can be accommodated.

There is evidence that the cortical folding pattern is subject to developmental changes. After 5 months in utero, cortical folds appear and continue to develop at least into the first postpartum year [17]. During early childhood, the degree of gyrification further increases and has so far been assumed to stabilize thereafter. Post-mortem analyses by Armstrong et al. [18], however, observed a significant overshot in cortical folding until the first year followed by a reduction until adulthood.

This finding is supported by recent MRI-studies which have investigated GI-values during brain maturation. Raznahan et al. [19] demonstrated a global decrease in gyrification during adolescence. More recently, Mutlu et al. [20] showed that GI-values declined between 6–29 years of age in frontal and parietal cortices which is consistent with data from Su and colleagues [21] who applied a novel approach of gyrification measurement towards a small sample of children and adolescents. Finally, data by Hogstrom et al. [22] suggest that modfications in gyrification continue until old age.

In the present study, we sought to comprehensively characterize the development of gyrification during adolescence through investigating whole-brain GI-values in MRI-data. In addition, we obtained GM-parameters (cortical thickness, volume and surface area) as well as WM-volume estimates to determine the relationship between age-dependent changes in gyrification and GM/WM parameters. Our results show widespread reduction in GI-values which occur in overlapping but also distinct areas of GM-change, such as in precentral, temporal and frontal regions, which highlight the ongoing anatomical modification of the cerebral cortex during adolescence.

Materials and Methods


85 right-handed participants (36 males and 49 females) between the ages of 12 and 23 years were recruited from local high-schools and the Goethe University Frankfurt and were screened for the presence of psychiatric disorders, neurological illness and substance abuse. Written informed consent was obtained from all participants. For participants younger than 18 years, written consent was given by their parents. The Hamburger-Wechsler intelligence testing battery (HAWI-E/K) [23], [24] was performed. Six participants were excluded cause of lacking or incomplete MRI-data. The study was approved by the ethics board of the Goethe-University Frankfurt.

MR Data Acquisition

Structural magnetic resonance images were obtained with a 3-Tesla Siemens Trio scanner (Siemens, Erlangen, Germany), using a CP head coil for RF transmission and signal reception. We used a T1-weighted three-dimensional (3D) Magnetization Prepared Rapid Acquisition Gradient Echo (MPRAGE) sequence with the following parameters: time repetition (TR): 2250 ms., time echo (TE): 2.6 ms., field of view (FOV): 256×256 mm3, slices: 176 and a voxel size of 1×1×1.1 mm3.

Surface reconstruction

MRI-data were processed with the surface and volume pipeline of the FreeSurfer-software version 5.1.0 ( [25], [26] and estimates of cortical thickness, GM- and WM- volume, cortical surface area, the 3-D local gyrification Index (lGI) and estimated intracranial volume (eTIV) were obtained. The standard FreeSurfer pipeline was followed and automatically reconstructed surfaces were inspected for accuracy and if necessary, manual interventions using FreeSurfer correction tools were used.

Pre-processing included Talairach transformation, motion correction, intensity normalization, non brain tissue removing, segmentation and tessellation of the gray and white matter boundary, automatically topology correction and surface deformation and is described in more detail elsewhere [25], [27][29]. In addition, a spherical atlas registration, inflation and a gyral/sulcal based parcellation of the cortical surface was performed for inter-individual analyses which yielded 33 cortical areas per hemisphere [30].

Cortical thickness, Cortical Surface Area and GM-volume

Cortical thickness was measured as the distance between the WM- boundary and the GM-matter surface at each point (vertex) on the tesselated surface [27]. Cortical surface area maps were generated through area estimations of each triangle in a standardized surface tessellation [31]. Area estimations were mapped back to the individual cortical space by means of a spherical atlas registration [32]. This yielded vertex-by-vertex estimates of the relative areal expansion or compression [33]. Estimates of GM-volume were derived from cortical thickness measures and the area around the corresponding vertex on the cortical surface [34].

3-D local gyrification index (lGI)

A 3-D lGI was computed [35] which has been employed in previous MR-studies [36], [37]. In brief, the lGI involves a 3-D reconstruction of the cortical surface where the degree of gyrification is defined as the amount of cortex surface buried within the sulcal folds compared with the amount of visible cortex in circular regions of interest [38]. In the first step, a triangulated outer surface which tightly wraps the pial surface was created through a morphological closing procedure. After converting the pial mesh into a binary volume, we used a diameter of 15 mm to close the main sulci for for generating the sphere [35]. For creating the circular region of interest (ROI), we choose a radius of 25 mm to include more than one sulcus to obtain an optimal resolution [38]. The initial lGI-values of a vertex were defined as the ratio between the surface of the outer ROI and the surface on the pial surface. For statistical comparisons, the outer lGI-values were mapped back to the individual coordinate system which reduced interindividual sulcal misalignment [35].


The regional WM-volume below parcellated cortical GM-regions were estimated. Each white matter voxel was labeled to the nearest cortical GM-voxel with a distance limit of 5 mm resulting in 33 WM-volumes of the corresponding 33 gyral labeled GM areas [39] which has been used in previous studies [9], [40].

Estimated intracranial volume (eTIV)

The estimated intracranial Volume (eTIV) in the FreeSurfer pipeline was derived from an atlas normalization procedure. Through the Atlas Scaling Factor (ASF), which represents a volume-scaling factor to match an individual to an atlas target, calculations of each eTIV were performed [41].

Statistical Analysis

The analyses steps are summarized in Figure 1. Surfaces of the right and left hemispheres of all 79 participants were averaged and individual surfaces were resampled into the average spherical coordinate system. To increase the signal to noise ratio, we used 20 mm full-width at half maximum (FWHM) smoothing for the estimation of cortical thickness, GM- volume and cortical surface area and 5 mm FWHM for the lGI.

Figure 1 

Analyses Steps for the lGI-Values and Correlations with Anatomical Parameters (GM/WM-volume, Cortical Surface Area and Cortical Thickness).

In the first step, we investigated whole-brain lGI-values, cortical thickness, cortical surface area and volume of GM in a vertex-by-vertex analysis. A General Linear Model (GLM) was employed to analyze the effect of age on the different anatomical parameters (lGI, cortical thickness, cortical surface area and GM-volume). All analyses were performed while controlling for the effects of gender and eTIV. We employed a false discovery rate approach (FDR) [42] to correct for multiple comparisons with a criterion for cortical thickness, surface area and GM-volume of q 0.05 and q 0.005 for lGI estimates. Different statistical thresholds were chosen because of the widespread, age-dependent changes in lGI-values compared to cortical thickness, cortical surface area and GM-volume. In addition we analyzed age2 and age3 effects for all anatomical parameters which were controlled for the influence of age, gender and eTIV.

To obtain estimates of area size, we selected vertices with the largest lGI-values and their corresponding Talairach coordinates and applied the automatic mri_surfcluster function in FreeSurfer ( In addition, Cohen’s d [43] was obtained for brain areas with the largest age-dependent changes through the comparison between the mean values in the youngest (age: 12–14, n = 13) and oldest participant group (age: 21–23, n = 18). Effects sizes are reported in figure legends.

In a second step, we examined Pearson correlation coefficients between age-dependent lGI-effects and changes in cortical thickness, cortical surface area and GM/WM-volume. To include WM-volume data, parcellation based regional analyses were performed. Four vertices from the vertex-by-vertex analyses per hemisphere with pronounced age-lGI effects (statistical threshold p<10−4) were assigned to FreeSurfers gyral based areas [30] and for the corresponding labels mean cortical thickness, GM/WM-volume and cortical surface area were extracted.


Vertex-by-vertex analyses of age-dependent changes in lGI

lGI-values decreased with age in 12 clusters in the left and 10 clusters in the right hemisphere (FDR at 0.005) (Figure 2 and ​and3,3, Table 1). Brain areas with the largest lGI-reductions were localized to the left precentral (area size = 22211.63 mm2, p = 10−8.42, BA 6 and 7), left superior-frontal (area size = 3804.76 mm2, p = 10−5.69, BA 10), left inferior-temporal (area size = 2477.53 mm2, p = 10−4.61, BA 19, 20 and 37), left lateral-orbitofrontal (area size = 1834.36 mm2, p = 10−4.45, BA 47 and 11) and right precentral cortex (area size = 12152.39 mm2, p = 10−7.47, BA 6 and 7), right pars triangularis (area size = 271.76 mm2, p = 10−4.57, BA 10 and 46), right rostral-middlefrontal (area size = 1200.69 mm2, p = 10−4.57, BA 9) and superior parietal (area size = 1834.36 mm2, p = 10−4.26, BA 19 and 39). No significant effects of gender were found for changes in lGI-values at a FDR at 0.005 and age-related reductions in gyrification followed nonlinear (cubic) trajectories (Figure 3).

Figure 2 

Whole-Brain Analyses of the Local Gyrification Index (lGI) during Adolescence.
Figure 3 

Scatter plots for the nine brain areas with significant correlations between age and lGI-values.
Table 1 

Age-Related Decreases in Gyrification.

Vertex-by-vertex analyses of age-dependent changes in Cortical Thickness, GM-Volume and Cortical Surface Area

Cortical thickness decreased most prominently in the superior-frontal (area size = 2608.63 mm2, p = 10−7.13, BA 6, 8 and 9) and rostral-middle-frontal (area size = 12859.08 mm2, p = 10−6.08, BA 11, 44, 45 and 46) cortices in the left hemisphere and in the precentral cluster in the right hemisphere (area size = 14735.38 mm2, p = 10−6.16, BA 6, 44 and 45) (Figure 4). The cortical thickness decrease can be described by a cubic trajectory (R2 = 0.191 for left rostral-middle-frontal, R2 = 0.126 for left superior-frontal and R2  = 0.134 for right precentral clusters). Moreover, we found age-dependent, bilateral decreases in GM-volume which were localized to the superior-frontal (area size = 45212.15 mm2, p = 10−7.60, BA 6, 8 and 9) lobe in the left hemisphere and to the pars orbitalis (area size = 19200.11 mm2, p = 10−6.68, BA 44, 45 and 47) and to the inferior-parietal (area size = 16614.72 mm2, p = 10−5.03 BA 19 and 39) lobe of the right hemisphere (Figure 4). GM-volume reductions followed cubic trajectories (R2 = 0.132 for left superior-frontal, R2 = 0.185 for right pars orbitalis and R2 = 0.204 for right inferior parietal clusters).

Figure 4 

Comparison of Age-Related Changes between GM-Volume, Cortical Thickness, Cortical Surface Area and Gyrification.

For surface area, we found a significant reduction in precentral (area size = 2296.99 mm2, p = 10−9.64, BA 4), caudal middle frontal (area size = 609.mm2, p = 10−6.03, BA 6) and supramarginal (area size = 1647.24 mm2, p = 10−4.88, BA 22) clusters in the left hemisphere. Surface area decreased in the right hemisphere most prominently in precentral (area size = 1371.37 mm2, p = 10−6.34, BA 4), inferior parietal (area size = 1248.36 mm2, p = 10−5.99, BA 7) and superior parietal (area size = 652.77 mm2, p = 10−4.11, BA 7) cortices (Figure 4). Reductions in surface area were best described by a cubic trajectory (R2 = 0.095 for left precentral, R2 = 0.026 left caudal-middle frontal, R2 = 0.024 left supramarginal, R2 = 0.116 right hemisphere, R2 = 0.156 right superior-parietal and R2  = 0.046 for right precentral clusters). No significant effects of gender were found for changes in cortical thickness, GM-volume and surface area at a FDR at 0.005

Correlations between Gyrification, Cortical Thickness, Surface Area and GM/WM-Volume

To test for relationships between lGI-values and changes in GM/WM, 8 areas with the largest age-dependent changes in gyrification were selected and lGI-values were correlated with cortical thickness, cortical surface area and GM/WM-Volume (Figure 5, Table 2). We found large and positive correlations between cortical surface area and GM-volume with lGI-values. Such a relationship was not found for correlations between cortical thickness and lGI-estimates. Increased WM-volume also showed a significant albeit weaker relationship than GM-volume and surface area with enhanced gyrification in several frontal regions and in parietal cortex.

Figure 5 

Based on the FreeSurfers Desikan labeling, eight regions of interest (ROI’s) were selected to analyze the relationships between lGI, Cortical Thickness, GM-volume, Cortical Surface Area and WM-volume.
Table 2 

Correlations between mean lGI-Values with Thickness, WM-, GM-Volume and Surface Area.

Non-Linear Relationships between Changes in Anatomical Parameters and Age: A Vertex-by-Vertex Analyses


We found 16 (left hemisphere) and 7 Clusters (hemisphere) where age2 and lGI were negatively correlated (Figure S1). The strongest age 2 effects on lGI were localized in left superior-frontal (area size = 2147.01 mm2, p = 10−5.48, BA 8, 9 and 10), left superior-parietal (area size = 5233.35 mm2, p = 10−4.51, BA 1, 2, 3, and 4) and left pericalcarine (area size = 243.34 mm2, p = 10−3.80, BA 17) clusters. For the right hemisphere, effects were observed in a precentral region (area size = 1165.59 mm2, p = 10−4.81, BA 1, 2, 3, 4, and 6), postcentral (area size = 465.07 mm2, p = 10−3.53, BA 1, 2 and 3) and in superiorfrontal cortices (area size = 330.55 mm2, p = 10−3.48, BA 8).

Cubic effects of age on lGI were found in 18 (left hemisphere) and 7 Clusters (right hemisphere). Regions with the strongest cubic effects were localized in a large superior-frontal (area size = 5598.96 mm2, p = 10−6.54, BA 8, 9, 10, 11, 45, 46 and 47), superior-parietal (area size = 11513.02 mm2, p = 10−6.11, BA 1, 2, 3, 4, 5, 6, 7, 8 and 9) and pericalcarine (area size = 292.35 mm2, p = 10−3.73, BA 17) cluster for the left hemisphere. In the right hemisphere, strongest cubic age and lGI relations were found in a precentral (area size = 5862.33 mm2, p = 10−5.52, BA 6, 4, 5, and 7), caudal-middlefrontal (area size = 503.66 mm2, p = 10−3.56, BA 8 and 9) and middle-temporal cluster (area size = 152.44 mm2, p = 10−2.98, BA 21).


Age2 effects on GMV were confined to the left hemisphere (Figure S2). Strongest effects were seen in extended parts of the pars opercularis (area size = 630.89 mm2, p = 10−4.35, BA 13, 44 and 45), paracentral (area size = 495.23 mm2, p = 10−4.11, BA 4, 6 and 31) and inferior-parietal (area size = 144.45 mm2, p = 10−3.71, BA 39 and 22) cortices.

Cubic age effects on GMV were located in 3 cortices in the left hemisphere. One cluster in the posterior parts of the gyrus cinguli (area size = 175.00 mm2, p = 10−4.55, BA 31), a part of the gyrus inferior frontalis-pars opercularis- (area size = 124.78 mm2, p = 10−4.25, BA 44) and the banks of superior temporal sulcus (area size = 7.12 mm2, p = 10−3.61, BA 39) were characterized by a significant age3 and lGI relationship (Figure S2).

CT/SA: No significant age2/age3 effects we found for CT and SA.


The results of our study highlight widespread changes in the gyrification-pattern of the cerebral cortex during adolescence. Previous post-mortem [18] and MRI-studies [19][21] indicated a decrease of lGI-values during later developmental periods but the extent of change, the brain regions involved and the relationship with concurrent anatomical process have remained unclear. Cortical areas which were characterized by the strongest reductions in lGI-values were precentral, temporal and frontal regions. These brain areas overlapped only partially with regions characterized by changes in GM and effect sizes were in the range and above for cortical thickness and GM-volume, suggesting that the observed modifications in gyrification represent an additional, important modification of the cerebral cortex during adolescence.

Cortical Regions of IGl-Changes

The largest cortical region characterized by reductions in gyrification was a cluster in the precentral cortex which included BA 3, 6 and 7. In comparison, changes in the thickness and volume of GM were focused over frontal (BA 8 and 9) and temporal (BA 20 and 21) cortices, which is consistent with data from previous longitudinal studies [6] but overlapped only partially with decreased lGI-values.

Although the precentral cluster, which extended to pre-/post-central gyrus, supramarginal gyrus as well as to the superior parietal cortex, has been less consistently involved in adolescent brain maturation, there is evidence to suggest that these brain areas may be related to ongoing changes in cognition and behavior. A recent study by Ramsden et al. [44] demonstrated that fluctuations in intelligence during adolescence are closely related to GM-changes in left motor speech regions. Similarly, there is ongoing improvement in motor cortex as revealed through studies with transcranial magnetic stimulation (TMS) [45] and EEG [46]. Finally, BA 7 is critical for the developmental of cortical networks underlying higher cognitive functions during adolescence, such as working memory (WM), because BOLD-activity in the superior parietal cortex shows substantial developmental increases during the manipulation of WM-items [47].

A second region of pronounced changes in IGl-values was the frontal cortex which has been consistently linked to changes in anatomy and behavior during adolescence. In the present study, decreased lGI-values were found in the frontal pole (BA 10), orbitofrontal cortex (BA 11) and the inferior frontal gyrus (BA 47). A large body of work has indicated that these regions are centrally involved in the behavioral modifications during adolescence, such as the improvements in cognitive inhibition [48], risk-taking [49] and mentalizing [50].

Finally, substantial reductions in gyrification were found in a cluster corresponding to BA 19, 20 and 37 which comprises early visual areas and cortical regions dedicated to object recognition. In addition to modifications in higher cognitive functions, adolescence is also associated with improvements in neural oscillations elicited by simple and complex visual stimuli [51], [52] as well as with maturation of object processing in the ventral stream [53].

Strong quadratic effects of age on lGI were found in left superior-frontal (BA 8, 9 and 10) and righthemispheric frontal (BA 8) clusters, which is in line with a previous study by (Hogstrom et al. [22]. Cubic age-lGI relationships are localized in left superior-frontal (BA 8, 9, 10, 11, 45, 46 and 47), superior-parietal (BA 1, 2, 3, 4, 5, 6, 7, 8 and 9), right caudal-middlefrontal (BA 8 and 9) and middle-temporal (BA 21) areas.

The current data thus provide a novel perspective on regions involved in gyrification development during adolescence which overall are characterized by a linear developmental trajectory with some regions showing curvilinear and cubic effects. Previous studies with smaller sample sizes [20], [21] identified predominantly changes in GI-values in temporal, parietal and frontal regions. In addition, Mutlu and colleagues [20] observed a steeper lGI decrease with age in males than females in prefrontal regions which was not confirmed by the present study.

Development of Cortical Folding during Adolescence: Relationship with GM/WM-change

Several mechanisms have been proposed for the changes in gyrification during development [54]. Van Essen [55] suggested that the folding pattern of the cerebral cortex can be explained by the mechanical tension along axons. According to this theory, the formation of gyri is the result of mechanical forces between densely linked regions as tension pulls strongly interconnected regions together. In addition, alternative accounts emphasized the role of differential growth between inner and outer cortical layers [17]. Finally, there is evidence that cortical folding is under genetic control [56] and that sex-differences exist in the mature cortex [57].

While the current study does not allow insights into the mechanisms underlying the reductions in gyrification during adolescence, comparison with changes in GM- and WM-parameters may be important for the question whether the observed changes in cortical folding are influenced by ongoing anatomical modifications. An important finding of the current study is that the reductions in lGI-values occur in cortical regions which are largely distinct from reductions in the volume and thickness of GM. Correlations between lGI-values in regions which were characterized by pronounced age-dependent decreases and GM/WM-parameters suggest, however, that the degree of cortical folding is nonetheless related to GM-volume and surface area. Specifically, we observed a positive relationship between increased lGI-values with surface area and volume of GM. Interestingly, this was not the case for the thickness of GM. Finally, WM-volume also contributed to higher lGI-values in 5 out of 7 cortical regions.

Gyrification, Behavior and Psychopathology

Despite the widespread reductions in cortical folding during adolescence and the large effect sizes associated with decreased lGI-values, the implications for changes in cognition and behaviour during adolescence remain to be established. Previous research has indicated that individual differences in cortical folding in frontal regions influence executive processes in adults [58] and behavioral modifications, such as meditation [59], impact on gyrification, suggesting a role of cortical folding in cognition and experience-dependent plasticity.

Furthermore, there is a large body of evidence that gyrification patterns are associated with psychopathology which underlines the potential importance of understanding developmental changes in gyrification and the relationship to cognition and behavior. Several neurodevelopmental disorders, such as Williams Syndrome (WS) and Autism Spectrum Disorders (ASDs), are associated with abnormal cortical folding patterns. Specifically, participants with WS are characterized by reductions in the depth of sulci in parieto-occipital regions which are prominently involved in the visuo-constructive deficits [60]. In contrast, gyrification patterns in ASDs are characterized by increased folding relative to normally developing children [61].

Schizophrenia is a severe psychiatric disorder with a typical onset during the transition from adolescence to adulthood which also involves aberrant gyrification. Post-mortem [62] and MRI-studies [63], [64] observed an increase in cortical folding, especially in the prefrontal cortex, which furthermore is predictive for the development of schizophrenia in at-risk subjects [65]. More recently, folding defects have also been shown to predict poor treatment response in first-episode psychosis [66].

Because our data strongly suggest that cortical folding undergoes major modifications during adolescence, one possibility is that in addition to early neurodevelopmental influences, abnormal brain development during adolescence contributes to the aberrant anatomy of the neocortex and the manifestation of cognitive dysfunctions and clinical symptoms.


The findings support the view that the adolescence involves fundamental changes in the architecture of the cerebral cortex. Specifically, we can show that cortical folding patterns undergo pronounced change which involves a reduction in gyrification across large areas of the cerebral cortex, in particular in precentral, frontal and temporal regions. Future studies need to establish the functional relevance of these modification for concurrent changes in behavior, cognition and physiology through correlations with neuropsychological data and functional brain imaging methods, such as fMRI and MEG.

Supporting Information

Figure S1

Nonlinear age effects on the local gyrification index (lGI) in a whole brain, vertex-by-vertex analyses projected onto an average template brain. Top Row: Age2 effects are illustrated for the left hemisphere (left) and right hemisphere (right) from lateral and medial views. Bottom Row: Correlations between age3 and lGI are shown for the left (left) and right hemisphere (right) from lateral and medial views. Blue colors indicate a significant decrease of lGI-values with increasing age, whereas warmer colors are coded for an increase in lGI. All analyses were performed by controlling for the effects of gender, eTIV and age (linear). Note: No significant correlations between age3 and lGI were found by controlling for the effects of gender, eTIV, age (linear) and age2.


Figure S2

Nonlinear age effects on GMV in a whole brain, vertex-by-vertex analyses projected onto an average template brain. Left: Age2 effects on GMV for the left hemisphere from lateral and medial view. Right: Effects of age3 are illustrated for the left hemisphere from lateral and medial view. Blue colors indicate a significant decrease of GMV with increasing age, whereas warmer colors are coded for an increase in GMV. All analyses were performed by controlling for the effects of gender, eTIV and age (linear). Note: No significant correlations between age3 and GMV were found by controlling for the effects of gender, eTIV, age (linear) and age2.



We would like to thank Sandra Anti for help with MRI data acquisition.

Funding Statement

This work was supported by the Max Planck Society (P. J. Uhlhaas) and by National Research Foundation of Korea funded by the Ministry of Education, Science and Technology (R32-10142, C.E. Han). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Blakemore SJ (2012) Imaging brain development: the adolescent brain. Neuroimage 61: 397–406. [PubMed]
2. Galvan A, Van Leijenhorst L, McGlennen KM (2012) Considerations for imaging the adolescent brain. Dev Cogn Neurosci 2: 293–302. [PubMed]
3. Giedd JN, Rapoport JL (2010) Structural MRI of pediatric brain development: what have we learned and where are we going? Neuron 67: 728–734. [PMC free article] [PubMed]
4. Huttenlocher PR (1984) Synapse elimination and plasticity in developing human cerebral cortex. Am J Ment Defic 88: 488–496. [PubMed]
5. Giedd JN, Jeffries NO, Blumenthal J, Castellanos FX, Vaituzis AC, et al. (1999) Childhood-onset schizophrenia: progressive brain changes during adolescence. Biol Psychiatry 46: 892–898. [PubMed]
6. Gogtay N, Giedd JN, Lusk L, Hayashi KM, Greenstein D, et al. (2004) Dynamic mapping of human cortical development during childhood through early adulthood. Proc Natl Acad Sci U S A 101: 8174–8179. [PMC free article] [PubMed]
7. Paus T (2010) Growth of white matter in the adolescent brain: myelin or axon? Brain Cogn 72: 26–35. [PubMed]
8. Paus T, Zijdenbos A, Worsley K, Collins DL, Blumenthal J, et al. (1999) Structural maturation of neural pathways in children and adolescents: in vivo study. Science 283: 1908–1911. [PubMed]
9. Tamnes CK, Ostby Y, Fjell AM, Westlye LT, Due-Tonnessen P, et al. (2010) Brain maturation in adolescence and young adulthood: regional age-related changes in cortical thickness and white matter volume and microstructure. Cereb Cortex 20: 534–548. [PubMed]
10. Colby JB, Van Horn JD, Sowell ER (2011) Quantitative in vivo evidence for broad regional gradients in the timing of white matter maturation during adolescence. Neuroimage 54: 25–31. [PMC free article] [PubMed]
11. Petanjek Z, Judas M, Simic G, Rasin MR, Uylings HB, et al. (2011) Extraordinary neoteny of synaptic spines in the human prefrontal cortex. Proc Natl Acad Sci U S A 108: 13281–13286. [PMC free article] [PubMed]
12. Lebel C, Beaulieu C (2011) Longitudinal development of human brain wiring continues from childhood into adulthood. J Neurosci 31: 10937–10947. [PubMed]
13. Raznahan A, Lerch JP, Lee N, Greenstein D, Wallace GL, et al. (2011) Patterns of coordinated anatomical change in human cortical development: a longitudinal neuroimaging study of maturational coupling. Neuron 72: 873–884. [PubMed]
14. Uhlhaas PJ, Singer W (2011) The development of neural synchrony and large-scale cortical networks during adolescence: relevance for the pathophysiology of schizophrenia and neurodevelopmental hypothesis. Schizophr Bull 37: 514–523. [PMC free article] [PubMed]
15. Paus T, Keshavan M, Giedd JN (2008) Why do many psychiatric disorders emerge during adolescence? Nat Rev Neurosci 9: 947–957. [PMC free article] [PubMed]
16. Rakic P (1995) A small step for the cell, a giant leap for mankind: a hypothesis of neocortical expansion during evolution. Trends Neurosci 18: 383–388. [PubMed]
17. Caviness VS Jr (1975) Mechanical model of brain convolutional development. Science 189: 18–21. [PubMed]
18. Armstrong E, Schleicher A, Omran H, Curtis M, Zilles K (1995) The ontogeny of human gyrification. Cereb Cortex 5: 56–63. [PubMed]
19. Raznahan A, Shaw P, Lalonde F, Stockman M, Wallace GL, et al. (2011) How does your cortex grow? J Neurosci 31: 7174–7177. [PMC free article] [PubMed]
20. Mutlu AK, Schneider M, Debbane M, Badoud D, Eliez S, et al. (2013) Sex differences in thickness, and folding developments throughout the cortex. Neuroimage 82: 200–207. [PubMed]
21. Su S, White T, Schmidt M, Kao CY, Sapiro G (2013) Geometric computation of human gyrification indexes from magnetic resonance images. Hum Brain Mapp 34: 1230–1244. [PubMed]
22. Hogstrom LJ, Westlye LT, Walhovd KB, Fjell AM (2012) The Structure of the Cerebral Cortex Across Adult Life: Age-Related Patterns of Surface Area, Thickness, and Gyrification. Cereb Cortex . [PubMed]
23. Petermann F, Petermann U (2010) HAWIK-IV. Bern: Huber.
24. Tewes U (1991) HAWIE-R. Hamburg-Wechsler-Intelligenztest für Erwachsene. Bern: Huber.
25. Dale AM, Fischl B, Sereno MI (1999) Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage 9: 179–194. [PubMed]
26. Fischl B, van der Kouwe A, Destrieux C, Halgren E, Segonne F, et al. (2004) Automatically parcellating the human cerebral cortex. Cereb Cortex 14: 11–22. [PubMed]
27. Fischl B, Dale AM (2000) Measuring the thickness of the human cerebral cortex from magnetic resonance images. Proc Natl Acad Sci U S A 97: 11050–11055. [PMC free article] [PubMed]
28. Fischl B, Sereno MI, Dale AM (1999) Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system. Neuroimage 9: 195–207. [PubMed]
29. Fischl B, Liu A, Dale AM (2001) Automated manifold surgery: constructing geometrically accurate and topologically correct models of the human cerebral cortex. IEEE Trans Med Imaging 20: 70–80. [PubMed]
30. Desikan RS, Segonne F, Fischl B, Quinn BT, Dickerson BC, et al. (2006) An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31: 968–980. [PubMed]
31. Joyner AH, J CR, Bloss CS, Bakken TE, Rimol LM, et al. (2009) A common MECP2 haplotype associates with reduced cortical surface area in humans in two independent populations. Proc Natl Acad Sci U S A 106: 15483–15488. [PMC free article] [PubMed]
32. Bakken TE, Roddey JC, Djurovic S, Akshoomoff N, Amaral DG, et al. (2012) Association of common genetic variants in GPCPD1 with scaling of visual cortical surface area in humans. Proc Natl Acad Sci U S A 109: 3985–3990. [PMC free article] [PubMed]
33. Rimol LM, Agartz I, Djurovic S, Brown AA, Roddey JC, et al. (2010) Sex-dependent association of common variants of microcephaly genes with brain structure. Proc Natl Acad Sci U S A 107: 384–388. [PMC free article] [PubMed]
34. Rimol LM, Nesvag R, Hagler DJ Jr, Bergmann O, Fennema-Notestine C, et al. (2012) Cortical volume, surface area, and thickness in schizophrenia and bipolar disorder. Biol Psychiatry 71: 552–560. [PubMed]
35. Schaer M, Cuadra MB, Tamarit L, Lazeyras F, Eliez S, et al. (2008) A surface-based approach to quantify local cortical gyrification. IEEE Trans Med Imaging 27: 161–170. [PubMed]
36. Palaniyappan L, Mallikarjun P, Joseph V, White TP, Liddle PF (2011) Folding of the prefrontal cortex in schizophrenia: regional differences in gyrification. Biol Psychiatry 69: 974–979. [PubMed]
37. Schaer M, Glaser B, Cuadra MB, Debbane M, Thiran JP, et al. (2009) Congenital heart disease affects local gyrification in 22q11.2 deletion syndrome. Dev Med Child Neurol 51: 746–753. [PubMed]
38. Schaer M, Cuadra MB, Schmansky N, Fischl B, Thiran JP, et al. (2012) How to measure cortical folding from MR images: a step-by-step tutorial to compute local gyrification index. J Vis Exp e3417. [PMC free article] [PubMed]
39. Fjell AM, Westlye LT, Greve DN, Fischl B, Benner T, et al. (2008) The relationship between diffusion tensor imaging and volumetry as measures of white matter properties. Neuroimage 42: 1654–1668. [PMC free article] [PubMed]
40. Salat DH, Greve DN, Pacheco JL, Quinn BT, Helmer KG, et al. (2009) Regional white matter volume differences in nondemented aging and Alzheimer’s disease. Neuroimage 44: 1247–1258. [PMC free article] [PubMed]
41. Buckner RL, Head D, Parker J, Fotenos AF, Marcus D, et al. (2004) A unified approach for morphometric and functional data analysis in young, old, and demented adults using automated atlas-based head size normalization: reliability and validation against manual measurement of total intracranial volume. Neuroimage 23: 724–738. [PubMed]
42. Genovese CR, Lazar NA, Nichols T (2002) Thresholding of statistical maps in functional neuroimaging using the false discovery rate. Neuroimage 15: 870–878. [PubMed]
43. Cohen J (1988) Statistical power analysis for the behavioral sciences. Hillsdale, NJ Lawrence Earlbaum Associates.
44. Ramsden S, Richardson FM, Josse G, Thomas MSC, Ellis C, et al. (2011) Verbal and non-verbal intelligence changes in the teenage brain. Nature 479: 113–116. [PMC free article] [PubMed]
45. Garvey MA, Ziemann U, Bartko JJ, Denckla MB, Barker CA, et al. (2003) Cortical correlates of neuromotor development in healthy children. Clin Neurophysiol 114: 1662–1670. [PubMed]
46. Farmer SF, Gibbs J, Halliday DM, Harrison LM, James LM, et al. (2007) Changes in EMG coherence between long and short thumb abductor muscles during human development. J Physiol 579: 389–402. [PMC free article] [PubMed]
47. Crone EA, Wendelken C, Donohue S, van Leijenhorst L, Bunge SA (2006) Neurocognitive development of the ability to manipulate information in working memory. Proc Natl Acad Sci U S A 103: 9315–9320. [PMC free article] [PubMed]
48. Rubia K, Smith AB, Taylor E, Brammer M (2007) Linear age-correlated functional development of right inferior fronto-striato-cerebellar networks during response inhibition and anterior cingulate during error-related processes. Hum Brain Mapp 28: 1163–1177. [PubMed]
49. Galvan A, Hare TA, Parra CE, Penn J, Voss H, et al. (2006) Earlier development of the accumbens relative to orbitofrontal cortex might underlie risk-taking behavior in adolescents. J Neurosci 26: 6885–6892. [PubMed]
50. Blakemore SJ (2008) Development of the social brain during adolescence. Q J Exp Psychol (Hove) 61: 40–49. [PubMed]
51. Werkle-Bergner M, Shing YL, Muller V, Li SC, Lindenberger U (2009) EEG gamma-band synchronization in visual coding from childhood to old age: evidence from evoked power and inter-trial phase locking. Clin Neurophysiol 120: 1291–1302. [PubMed]
52. Uhlhaas PJ, Roux F, Singer W, Haenschel C, Sireteanu R, et al. (2009) The development of neural synchrony reflects late maturation and restructuring of functional networks in humans. Proc Natl Acad Sci U S A 106: 9866–9871. [PMC free article] [PubMed]
53. Golarai G, Ghahremani DG, Whitfield-Gabrieli S, Reiss A, Eberhardt JL, et al. (2007) Differential development of high-level visual cortex correlates with category-specific recognition memory. Nat Neurosci 10: 512–522. [PMC free article] [PubMed]
54. Zilles K, Palomero-Gallagher N, Amunts K (2013) Development of cortical folding during evolution and ontogeny. Trends Neurosci 36: 275–284. [PubMed]
55. Van Essen DC (1997) A tension-based theory of morphogenesis and compact wiring in the central nervous system. Nature 385: 313–318. [PubMed]
56. Rogers J, Kochunov P, Zilles K, Shelledy W, Lancaster J, et al. (2010) On the genetic architecture of cortical folding and brain volume in primates. Neuroimage 53: 1103–1108. [PMC free article] [PubMed]
57. Luders E, Narr KL, Thompson PM, Rex DE, Jancke L, et al. (2004) Gender differences in cortical complexity. Nat Neurosci 7: 799–800. [PubMed]
58. Fornito A, Yucel M, Wood S, Stuart GW, Buchanan JA, et al. (2004) Individual differences in anterior cingulate/paracingulate morphology are related to executive functions in healthy males. Cereb Cortex 14: 424–431. [PubMed]
59. Luders E, Kurth F, Mayer EA, Toga AW, Narr KL, et al. (2012) The unique brain anatomy of meditation practitioners: alterations in cortical gyrification. Front Hum Neurosci 6: 34. [PMC free article] [PubMed]
60. Kippenhan JS, Olsen RK, Mervis CB, Morris CA, Kohn P, et al. (2005) Genetic contributions to human gyrification: sulcal morphometry in Williams syndrome. J Neurosci 25: 7840–7846. [PubMed]
61. Jou RJ, Minshew NJ, Keshavan MS, Hardan AY (2010) Cortical gyrification in autistic and Asperger disorders: a preliminary magnetic resonance imaging study. J Child Neurol 25: 1462–1467. [PMC free article] [PubMed]
62. Vogeley K, Schneider-Axmann T, Pfeiffer U, Tepest R, Bayer TA, et al. (2000) Disturbed gyrification of the prefrontal region in male schizophrenic patients: A morphometric postmortem study. Am J Psychiatry 157: 34–39. [PubMed]
63. Kulynych JJ, Luevano LF, Jones DW, Weinberger DR (1997) Cortical abnormality in schizophrenia: an in vivo application of the gyrification index. Biol Psychiatry 41: 995–999. [PubMed]
64. Palaniyappan L, Liddle PF (2012) Aberrant cortical gyrification in schizophrenia: a surface-based morphometry study. J Psychiatry Neurosci 37: 399–406. [PMC free article] [PubMed]
65. Harris JM, Whalley H, Yates S, Miller P, Johnstone EC, et al. (2004) Abnormal cortical folding in high-risk individuals: a predictor of the development of schizophrenia? Biol Psychiatry 56: 182–189. [PubMed]
66. Palaniyappan L, Marques TR, Taylor H, Handley R, Mondelli V, et al. (2013) Cortical folding defects as markers of poor treatment response in first-episode psychosis. JAMA Psychiatry 70: 1031–1040. [PubMed]