A functional model of cortical gyri and sulci

Cortical Morphological Networks Differ Between Gyri and Sulci

This study explored how the human cortical folding pattern composed of convex gyri and concave sulci affected single-subject morphological brain networks, which are becoming an important method for studying the human brain connectome. We found that gyri-gyri networks exhibited higher morphological similarity, lower small-world parameters, and lower long-term test-retest reliability than sulci-sulci networks for cortical thickness- and gyrification index-based networks, while opposite patterns were observed for fractal dimension-based networks. Further behavioral association analysis revealed that gyri-gyri networks and connections between gyral and sulcal regions significantly explained inter-individual variance in Cognition and Motor domains for fractal dimension- and sulcal depth-based networks. Finally, the clinical application showed that only sulci-sulci networks exhibited morphological similarity reductions in major depressive disorder for cortical thickness-, fractal dimension-, and gyrification index-based networks. Taken together, these findings provide novel insights into the constraint of the cortical folding pattern to the network organization of the human brain.

Anatomy-Guided Spatio-Temporal Graph Convolutional Networks (AG-STGCNs) for Modeling Functional Connectivity Between Gyri and Sulci Across Multiple Task Domains

The cerebral cortex is folded as gyri and sulci, which provide the foundation to unveil anatomo-functional relationship of brain. Previous studies have extensively demonstrated that gyri and sulci exhibit intrinsic functional difference, which is further supported by morphological, genetic, and structural evidences. Therefore, systematically investigating the gyro-sulcal (G-S) functional difference can help deeply understand the functional mechanism of brain. By integrating functional magnetic resonance imaging (fMRI) with advanced deep learning models, recent studies have unveiled the temporal difference in functional activity between gyri and sulci. However, the potential difference of functional connectivity, which represents functional dependency between gyri and sulci, is much unknown. Moreover, the regularity and variability of the G-S functional connectivity difference across multiple task domains remains to be explored. To address the two concerns, this study developed new anatomy-guided spatio-temporal graph convolutional networks (AG-STGCNs) to investigate the regularity and variability of functional connectivity differences between gyri and sulci across multiple task domains. Based on 830 subjects with seven different task-based and one resting state fMRI (rs-fMRI) datasets from the public Human Connectome Project (HCP), we consistently found that there are significant differences of functional connectivity between gyral and sulcal regions within task domains compared with resting state (RS). Furthermore, there is considerable variability of such functional connectivity and information flow between gyri and sulci across different task domains, which are correlated with individual cognitive behaviors. Our study helps better understand the functional segregation of gyri and sulci within task domains as well as the anatomo-functional-behavioral relationship of the human brain.

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.

Hierarchical functional differences between gyri and sulci at different scales

Gyri and sulci are 2 fundamental cortical folding patterns of the human brain. Recent studies have suggested that gyri and sulci may play different functional roles given their structural and functional heterogeneity. However, our understanding of the functional differences between gyri and sulci remains limited due to several factors. Firstly, previous studies have typically focused on either the spatial or temporal domain, neglecting the inherently spatiotemporal nature of brain functions. Secondly, analyses have often been restricted to either local or global scales, leaving the question of hierarchical functional differences unresolved. Lastly, there has been a lack of appropriate analytical tools for interpreting the hierarchical spatiotemporal features that could provide insights into these differences. To overcome these limitations, in this paper, we proposed a novel hierarchical interpretable autoencoder (HIAE) to explore the hierarchical functional difference between gyri and sulci. Central to our approach is its capability to extract hierarchical features via a deep convolutional autoencoder and then to map these features into an embedding vector using a carefully designed feature interpreter. This process transforms the features into interpretable spatiotemporal patterns, which are pivotal in investigating the functional disparities between gyri and sulci. We evaluate the proposed framework on Human Connectome Project task functional magnetic resonance imaging dataset. The experiments demonstrate that the HIAE model can effectively extract and interpret hierarchical spatiotemporal features that are neuroscientifically meaningful. The analyses based on the interpreted features suggest that gyri are more globally activated, whereas sulci are more locally activated, demonstrating a distinct transition in activation patterns as the scale shifts from local to global. Overall, our study provides novel insights into the brain's anatomy-function relationship.

Frequency-specific functional difference between gyri and sulci in naturalistic paradigm fMRI

Disentangling functional difference between cortical folding patterns of gyri and sulci provides novel insights into the relationship between brain structure and function. Previous studies using resting-state functional magnetic resonance imaging (rsfMRI) have revealed that sulcal signals exhibit stronger high-frequency but weaker low-frequency components compared to gyral ones, suggesting that gyri may serve as functional integration centers while sulci are segregated local processing units. In this study, we utilize naturalistic paradigm fMRI (nfMRI) to explore the functional difference between gyri and sulci as it has proven to record stronger functional integrations compared to rsfMRI. We adopt a convolutional neural network (CNN) to classify gyral and sulcal fMRI signals in the whole brain (the global model) and within functional brain networks (the local models). The frequency-specific difference between gyri and sulci is then inferred from the power spectral density (PSD) profiles of the learned filters in the CNN model. Our experimental results show that nfMRI shows higher gyral-sulcal PSD contrast effect sizes in the global model compared to rsfMRI. In the local models, the effect sizes are either increased or decreased depending on frequency bands and functional complexity of the FBNs. This study highlights the advantages of nfMRI in depicting the functional difference between gyri and sulci, and provides novel insights into unraveling the relationship between brain structure and function.

Regularity and variability of functional brain connectivity characteristics between gyri and sulci under naturalistic stimulus

The human cerebral cortex is folded into two fundamentally anatomical units: gyri and sulci. Previous studies have demonstrated the genetical, structural, and functional differences between gyri and sulci, providing a unique perspective for revealing the relationship among brain function, cognition, and behavior. While previous studies mainly focus on the functional differences between gyri and sulci under resting or task-evoked state, such characteristics under naturalistic stimulus (NS) which reflects real-world dynamic environments are largely unknown. To address this question, this study systematically investigates spatio-temporal functional connectivity (FC) characteristics between gyri and sulci under NS using a spatio-temporal graph convolutional network model. Based on the public Human Connectome Project dataset of 174 subjects with four different runs of both movie-watching NS and resting state 7T functional MRI data, we successfully identify unique FC features under NS, which are mainly involved in visual, auditory, emotional and cognitive control, and achieve high discriminative accuracy 93.06 % to resting state. Moreover, gyral regions as well as gyro-gyral connections consistently participate more as functional information exchange hubs than sulcal ones among these networks. This study provides novel insights into the functional brain mechanism under NS and lays a solid foundation for accurately mapping the brain anatomy-function relationship.

Fractal Neurodynamics

The neuronal ongoing electrical activity in the brain network, the neurodynamics, reflects the structure and functionality of generating neuronal pools. The activity of neurons due to their excitatory and inhibitory projections is associated with specific brain functions. Here, the purpose was to investigate if the local ongoing electrical activity exhibits its characteristic spectral and fractal features in wakefulness and sleep across and within subjects. Moreover, we aimed to show that measures typical of complex systems catch physiological features missed by linear spectral analyses. For this study, we concentrated on the evaluation of the power spectral density (PSD) and Higuchi fractal dimension (HFD) measures. Relevant clinical impact of the specific features of neurodynamics identification stands primarily in the potential of classifying cortical parcels according to their neurodynamics as well as enhancing the effectiveness of neuromodulation interventions to cure symptoms secondary to neuronal activity unbalances.

Genetic Influence on Gyral Peaks

Genetic mechanisms have been hypothesized to be a major determinant in the formation of cortical folding. Although there is an increasing number of studies examining the heritability of cortical folding, most of them focus on sulcal pits rather than gyral peaks. Gyral peaks, which reflect the highest local foci on gyri and are consistent across individuals, remain unstudied in terms of heritability. To address this knowledge gap, we used high-resolution data from the Human Connectome Project (HCP) to perform classical twin analysis and estimate the heritability of gyral peaks across various brain regions. Our results showed that the heritability of gyral peaks was heterogeneous across different cortical regions, but relatively symmetric between hemispheres. We also found that pits and peaks are different in a variety of anatomic and functional measures. Further, we explored the relationship between the levels of heritability and the formation of cortical folding by utilizing the evolutionary timeline of gyrification. Our findings indicate that the heritability estimates of both gyral peaks and sulcal pits decrease linearly with the evolution timeline of gyrification. This suggests that the cortical folds which formed earlier during gyrification are subject to stronger genetic influences than the later ones. Moreover, the pits and peaks coupled by their time of appearance are also positively correlated in respect of their heritability estimates. These results fill the knowledge gap regarding genetic influences on gyral peaks and significantly advance our understanding of how genetic factors shape the formation of cortical folding. The comparison between peaks and pits suggests that peaks are not a simple morphological mirror of pits but could help complete the understanding of folding patterns.

Temporal fingerprints of cortical gyrification in marmosets and humans

Recent neuroimaging studies in humans have reported distinct temporal dynamics of gyri and sulci, which may be associated with putative functions of cortical gyrification. However, the complex folding patterns of the human cortex make it difficult to explain temporal patterns of gyrification. In this study, we used the common marmoset as a simplified model to examine the temporal characteristics and compare them with the complex gyrification of humans. Using a brain-inspired deep neural network, we obtained reliable temporal-frequency fingerprints of gyri and sulci from the awake rs-fMRI data of marmosets and humans. Notably, the temporal fingerprints of one region successfully classified the gyrus/sulcus of another region in both marmosets and humans. Additionally, the temporal-frequency fingerprints were remarkably similar in both species. We then analyzed the resulting fingerprints in several domains and adopted the Wavelet Transform Coherence approach to characterize the gyro-sulcal coupling patterns. In both humans and marmosets, sulci exhibited higher frequency bands than gyri, and the two were temporally coupled within the same range of phase angles. This study supports the notion that gyri and sulci possess unique and evolutionarily conserved features that are consistent across functional areas, and advances our understanding of the functional role of cortical gyrification.

Dynamic Visualization of Gyral and Sulcal Stereoelectroencephalographic Contacts in Humans

Stereoelectroencephalography (SEEG) is a neurosurgical method to survey electrophysiological activity within the brain to treat disorders such as Epilepsy. In this stereotactic approach, leads are implanted through straight trajectories to survey both cortical and sub-cortical activity.Visualizing the recorded locations covering sulcal and gyral activity while staying true to the cortical architecture is challenging due to the folded, three-dimensional nature of the human cortex.To overcome this challenge, we developed a novel visualization concept, allowing investigators to dynamically morph between the subjects' cortical reconstruction and an inflated cortex representation. This inflated view, in which gyri and sulci are viewed on a smooth surface, allows better visualization of electrodes buried within the sulcus while staying true to the underlying cortical architecture.Clinical relevance- These visualization techniques might also help guide clinical decision-making when defining seizure onset zones or resections for patients undergoing SEEG monitoring for intractable epilepsy.

Dynamic Visualization of Gyral and Sulcal Stereoelectroencephalographic contacts in Humans

Stereoelectroencephalography (SEEG) is a neurosurgical method to survey electrophysiological activity within the brain to treat disorders such as Epilepsy. In this stereotactic approach, leads are implanted through straight trajectories to survey both cortical and sub-cortical activity. Visualizing the recorded locations covering sulcal and gyral activity while staying true to the cortical architecture is challenging due to the folded, three-dimensional nature of the human cortex. To overcome this challenge, we developed a novel visualization concept, allowing investigators to dynamically morph between the subjects' cortical reconstruction and an inflated cortex representation. This inflated view, in which gyri and sulci are viewed on a smooth surface, allows better visualization of electrodes buried within the sulcus while staying true to the underlying cortical architecture.

Clinical relevance—

These visualization techniques might also help guide clinical decision-making when defining seizure onset zones or resections for patients undergoing SEEG monitoring for intractable epilepsy.

Modeling functional difference between gyri and sulci within intrinsic connectivity networks

Recently, the functional roles of the human cortical folding patterns have attracted increasing interest in the neuroimaging community. However, most existing studies have focused on the gyro-sulcal functional relationship on a whole-brain scale but possibly overlooked the localized and subtle functional differences of brain networks. Actually, accumulating evidences suggest that functional brain networks are the basic unit to realize the brain function; thus, the functional relationships between gyri and sulci still need to be further explored within different functional brain networks. Inspired by these evidences, we proposed a novel intrinsic connectivity network (ICN)-guided pooling-trimmed convolutional neural network (I-ptFCN) to revisit the functional difference between gyri and sulci. By testing the proposed model on the task functional magnetic resonance imaging (fMRI) datasets of the Human Connectome Project, we found that the classification accuracy of gyral and sulcal fMRI signals varied significantly for different ICNs, indicating functional heterogeneity of cortical folding patterns in different brain networks. The heterogeneity may be contributed by sulci, as only sulcal signals show heterogeneous frequency features across different ICNs, whereas the frequency features of gyri are homogeneous. These results offer novel insights into the functional difference between gyri and sulci and enlighten the functional roles of cortical folding patterns.

Gyral and sulcal connectivity in the human cerebral cortex

The rapid evolution of image acquisition and data analytic methods has established in vivo whole-brain tractography as a routine technology over the last 20 years. Imaging-based methods provide an additional approach to classic neuroanatomical studies focusing on biomechanical principles of anatomical organization and can in turn overcome the complexity of inter-individual variability associated with histological and tractography studies. In this work we propose a novel, reliable framework for determining brain tracts resolving the anatomical variance of brain regions. We distinguished 4 region types based on anatomical considerations: (i) gyral regions at borders between cortical communities; (ii) gyral regions within communities; (iii) sulcal regions at invariant locations across subjects; and (iv) other sulcal regions. Region types showed strikingly different anatomical and connection properties. Results allowed complementing the current understanding of the brain’s communication structure with a model of its anatomical underpinnings.

Normal development of the brain: a survey of joint structural-functional brain studies

Joint structural-functional (S-F) developmental studies present a novel approach to address the complex neuroscience questions on how the human brain works and how it matures. Joint S-F biomarkers have the inherent potential to model effectively the brain's maturation, fill the information gap in temporal brain atlases, and demonstrate how the brain's performance matures during the lifespan. This review presents the current state of knowledge on heterochronous and heterogeneous development of S-F links during the maturation period. The S-F relationship has been investigated in early-matured unimodal and prolonged-matured transmodal regions of the brain using a variety of structural and functional biomarkers and data acquisition modalities. Joint S-F unimodal studies have employed auditory and visual stimuli, while the main focus of joint S-F transmodal studies has been resting-state and cognitive experiments. However, nonsignificant associations between some structural and functional biomarkers and their maturation show that designing and developing effective S-F biomarkers is still a challenge in the field. Maturational characteristics of brain asymmetries have been poorly investigated by the joint S-F studies, and the results were partially inconsistent with previous nonjoint ones. The inherent complexity of the brain performance can be modeled using multifactorial and nonlinear techniques as promising methods to simulate the impact of age on S-F relations considering their analysis challenges.

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.

Joint Analysis of Functional and Structural Connectomes Between Preterm and Term Infant Brains via Canonical Correlation Analysis With Locality Preserving Projection

Preterm is a worldwide problem that affects infants' lives significantly. Moreover, the early impairment is more than limited to isolated brain regions but also to global and profound negative outcomes later, such as cognitive disorder. Therefore, seeking the differences of brain connectome between preterm and term infant brains is a vital step for understanding the developmental impairment caused by preterm. Existing studies revealed that studying the relationship between brain function and structure, and further investigating their differentiable connectomes between preterm and term infant brains is a way to comprehend and unveil the differences that occur in the preterm infant brains. Therefore, in this article, we proposed a novel canonical correlation analysis (CCA) with locality preserving projection (LPP) approach to investigate the relationship between brain functional and structural connectomes and how such a relationship differs between preterm and term infant brains. CCA is proposed to study the relationship between functional and structural connections, while LPP is adopted to identify the distinguishing features from the connections which can differentiate the preterm and term brains. After investigating the whole brain connections on a fine-scale connectome approach, we successfully identified 89 functional and 97 structural connections, which mostly contributed to differentiate preterm and term infant brains from the functional MRI (fMRI) and diffusion MRI (dMRI) of the public developing Human Connectome Project (dHCP) dataset. By further exploring those identified connections, the results innovatively revealed that the identified functional connections are short-range and within the functional network. On the contrary, the identified structural connections are usually remote connections across different functional networks. In addition, these connectome-level results show the new insights that longitudinal functional changes could deviate from longitudinal structural changes in the preterm infant brains, which help us better understand the brain-behavior changes in preterm infant brains.

Gyral-sulcal contrast in intrinsic functional brain networks across task performances

Functional mechanism of the brain and its relationship with the brain structural substrate have been an interest for multiple disciplines for centuries. Recently, gyri and sulci, two basic cortical folding patterns, have been demonstrated to act different functional roles. Specifically, a variety of functional MRI (fMRI) studies have consistently suggested that gyri represent a global functional center while sulci serve as a local functional unit under either resting state or task stimulus, which are further supported by brain structural analysis reporting that gyri have thicker cortex and denser long-distance axonal fibers. However, the consistency of such gyral-sulcal functional difference across different task stimuli, as well as its association with task conditions, remains to be explored. To this end, we used intrinsic networks as the testbed for cross-task comparison, and adopted a computational framework of dictionary learning and sparse representation of whole-brain fMRI signals to systematically examine the potential gyral-sulcal difference in signal representation residual (SRR) which reflected the degree of global functional communication. Using all seven task-based fMRI datasets in Human Connectome Project Q1 release, we found that within the intrinsic functional networks, the fMRI SRR was significantly smaller on gyral regions than on sulcal regions across different task stimuli, indicating that gyral regions were more involved in global functions of the brain and interregional communications. Moreover, the magnitudes of such gyral-sulcal difference varied across task conditions and intrinsic networks. Our work adds novel explanation and insight to the existing knowledge of functional differences between gyri and sulci.

Functional Parcellation of Human Brain Precuneus Using Density-Based Clustering

The human precuneus is involved in many high-level cognitive functions, which strongly suggests the existence of biologically meaningful subdivisions. However, the functional parcellation of the precuneus needs much to be investigated. In this study, we developed an eigen clustering (EIC) approach for the parcellation using precuneus-cortical functional connectivity from fMRI data of the Human Connectome Project. The EIC approach is robust to noise and can automatically determine the cluster number. It is consistently demonstrated that the human precuneus can be subdivided into six symmetrical and connected parcels. The anterior and posterior precuneus participate in sensorimotor and visual functions, respectively. The central precuneus with four subregions indicates a media role in the interaction of the default mode, dorsal attention, and frontoparietal control networks. The EIC-based functional parcellation is free of the spatial distance constraint and is more functionally coherent than parcellation using typical clustering algorithms. The precuneus subregions had high accordance with cortical morphology and revealed good functional segregation and integration characteristics in functional task-evoked activations. This study may shed new light on the human precuneus function at a delicate level and offer an alternative scheme for human brain parcellation.

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.

Differential Patterns of Gyral and Sulcal Morphological Changes During Normal Aging Process

The cerebral cortex is a highly convoluted structure with distinct morphologic features, namely the gyri and sulci, which are associated with the functional segregation or integration in the human brain. During the lifespan, the brain atrophy that is accompanied by cognitive decline is a well-accepted aging phenotype. However, the detailed patterns of cortical folding change during aging, especially the changing age-dependencies of gyri and sulci, which is essential to brain functioning, remain unclear. In this study, we investigated the morphology of the gyral and sulcal regions from pial and white matter surfaces using MR imaging data of 417 healthy participants across adulthood to old age (21–92 years). To elucidate the age-related changes in the cortical pattern, we fitted cortical thickness and intrinsic curvature of gyri and sulci using the quadratic model to evaluate their age-dependencies during normal aging. Our findings show that comparing to gyri, the sulcal thinning is the most prominent pattern during the aging process, and the gyrification of pial and white matter surfaces were also affected differently, which implies the vulnerability of functional segregation during aging. Taken together, we propose a morphological model of aging that may provide a framework for understanding the mechanisms underlying gray matter degeneration.

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.

The Cerebral Cortex is Bisectionally Segregated into Two Fundamentally Different Functional Units of Gyri and Sulci

The human cerebral cortex is highly folded into diverse gyri and sulci. Accumulating evidences suggest that gyri and sulci exhibit anatomical, morphological, and connectional differences. Inspired by these evidences, we performed a series of experiments to explore the frequency-specific differences between gyral and sulcal neural activities from resting-state and task-based functional magnetic resonance imaging (fMRI) data. Specifically, we designed a convolutional neural network (CNN) based classifier, which can differentiate gyral and sulcal fMRI signals with reasonable accuracies. Further investigations of learned CNN models imply that sulcal fMRI signals are more diverse and more high frequency than gyral signals, suggesting that gyri and sulci truly play different functional roles. These differences are significantly associated with axonal fiber wiring and cortical thickness patterns, suggesting that these differences might be deeply rooted in their structural and cellular underpinnings. Further wavelet entropy analyses demonstrated the validity of CNN-based findings. In general, our collective observations support a new concept that the cerebral cortex is bisectionally segregated into 2 functionally different units of gyri and sulci.

Temporal Variability of Cortical Gyral-Sulcal Resting State Functional Activity Correlates With Fluid Intelligence

The human cerebral cortex is highly convoluted as convex gyri and concave sulci. In the past decades, extensive studies have consistently revealed substantial differences between gyri and sulci in terms of genetics, anatomy, morphology, axonal fiber connections, and function. Although interesting findings have been reported to date to elucidate the functional difference between gyri and sulci, the temporal variability of functional activity, which could explain individual differences in learning and higher-order cognitive functions, and as well as differences in gyri and sulci, remains to be explored. The present study explored the temporal variability of cortical gyral-sulcal resting state functional activity and its association with fluid intelligence measures on the Human Connectome Project dataset. We found that the temporal variance of resting state fMRI BOLD signal was significantly larger in gyri than in sulci. We also found that the temporal variability of certain regions including middle frontal cortex, inferior parietal lobe and visual cortex was positively associated with fluid intelligence. Moreover, those regions were predominately located in gyri rather than in sulci. This study reports initial evidence for temporal variability difference of functional activity between gyri and sulci, and its association with fluid intelligence measures, and thus provides novel insights to understand the mechanism and functional relevance of gyri and sulci.

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.

Deep Learning Models Unveiled Functional Difference Between Cortical Gyri and Sulci

It is largely unknown whether there is functional role difference between cortical gyral and sulcal regions. Recent advancements in neuroimaging studies demonstrate clear difference of structural connection profiles in gyral and sulcal areas, suggesting possible functional role difference in these convex and concave cortical regions. To explore and confirm such possible functional difference, we design and apply a powerful deep learning model of convolutional neural networks (CNN) that has been proven to be superior in learning discriminative and meaningful patterns on fMRI. By using the CNN model, gyral and sulcal fMRI signals are learned and predicted, and the prediction performance is adopted to demonstrate the functional difference between gyri and sulci. By using the Human Connectome Project (HCP) fMRI data and macaque brain fMRI data, an average of 83% and 90% classification accuracy has been achieved to separate gyral/sulcal HCP task fMRI signals at the population and individual subject level, respectively; 81% and 86% classification accuracy for resting state fMRI signals at the group and individual subject level, respectively. In addition, 78% classification accuracy has been achieved to separate gyral/sulcal resting state fMRI signals in macaque brains. Importantly, further analysis reveals that the discriminative features that are learned by CNNs to differentiate gyral/sulcal fMRI signals can be meaningfully interpreted, thus unveiling the fundamental functional difference between cortical gyri and sulci. That is, gyri are more global functional integration centers with simpler lower frequency signal components, while sulci are more local processing units with more complex higher frequency signal components.

A Cortical Folding Pattern-Guided Model of Intrinsic Functional Brain Networks in Emotion Processing

There have been increasing studies demonstrating that emotion processing in humans is realized by the interaction within or among the large-scale intrinsic functional brain networks. Identifying those meaningful intrinsic functional networks based on task-based functional magnetic resonance imaging (task fMRI) with specific emotional stimuli and responses, and exploring the underlying functional working mechanisms of interregional neural communication within the intrinsic functional networks are thus of great importance to understand the neural basis of emotion processing. In this paper, we propose a novel cortical folding pattern-guided model of intrinsic networks in emotion processing: gyri serve as global functional connection centers that perform interregional neural communication among distinct regions via long distance dense axonal fibers, and sulci serve as local functional units that directly communicate with neighboring gyri via short distance fibers and indirectly communicate with other distinct regions via the neighboring gyri. We test the proposed model by adopting a computational framework of dictionary learning and sparse representation of emotion task fMRI data of 68 subjects in the publicly released Human Connectome Project. The proposed model provides novel insights of functional mechanisms in emotion processing.

Migraine with Aura: Surface-Based Analysis of the Cerebral Cortex with Magnetic Resonance Imaging

Objective

Previous migraine studies have reported gray matter alterations in various cortical regions with conflicting results. This study aimed to explore a cortical morphometric difference in migraineurs with aura (MA) compared to healthy subjects (HS) and to delineate a possible difference between the cortical morphological features and different aura phenotypes.

Materials and Methods

Forty-eight MA and 30 HS that were balanced by sex, age, and educational level were selected for this study. T2-weighted and three-dimensional T1-weighted magnetic resonance imaging (MRI) of the brain were acquired using a 1.5T MRI scanner. Surface-based morphometry from the MRI data was used to identify differences between the MA and HS group, and then between MA subgroups. The MA group was subdivided into migraineurs who experienced only visual aura (MVA) and migraineurs who had visual, somatosensory and dysphasic symptoms (MVA+).

Results

The MVA+ group had significantly reduced cortical surface area of the left rostral middle frontal cortex compared with the MVA group (p < 0.001). Migraine patients had significantly reduced volume of the left fusiform gyrus relative to HS (p < 0.001). Also, the sulcal depth increased at the level of the left temporal pole in the MVA+ group relative to the MVA group (p < 0.001). The vertex-by-vertex analysis did not exhibit any significant difference in cortical thickness between MA and HS, and between MVA+ and MVA, when corrected for multiple comparisons.

Conclusion

Migraineurs with aura demonstrates different morphometric features from HS in multiple cortical regions. MVA+ have different morphometric features in the left frontal and temporal lobe relative to MVA, which could be a source of distinct symptoms and serve as potential biomarkers of different MA subtypes.

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.

A Survey of White Matter Neurons at the Gyral Crowns and Sulcal Depths in the Rhesus Monkey

Gyrencephalic brains exhibit deformations of the six neocortical laminae at gyral crowns and sulcal depths, where the deeper layers are, respectively, expanded and compressed. The present study addresses: (1) the degree to which the underlying white matter neurons (WMNs) observe the same changes at gyral crowns and sulcal depths; and (2) whether these changes are consistent or variable across different cortical regions. WMNs were visualized by immunohistochemistry using the pan-neuronal label NeuN, and their density was quantified in eight rhesus monkey brains for four regions; namely, frontal (FR), superior frontal gyrus (SFG), parietal (Par) and temporal (TE). In all four regions, there were about 50% fewer WMNs in the sulcal depth, but there was also distinct variability from region to region. For the gyral crown, we observed an average density per 0.21 mm2 of 82 WMNs for the FR, 51 WMNs for SFG, 80 WMNs for Par and 93 WMNs for TE regions. By contrast, for the sulcal depth, the average number of WMNs per 0.21 mm2 was 41 for FR, 31 for cingulate sulcus (underlying the SFG), 54 for Par and 63 for TE cortical regions. Since at least some WMNs participate in cortical circuitry, these results raise the possibility of their differential influence on cortical circuitry in the overlying gyral and sulcal locations. The results also point to a possible role of WMNs in the differential vulnerability of gyral vs. sulcal regions in disease processes, and reinforce the increasing awareness of the WMNs as part of a complex, heterogeneous and structured microenvironment.

Elucidating functional differences between cortical gyri and sulci via sparse representation HCP grayordinate fMRI data

The highly convoluted cerebral cortex is characterized by two different topographic structures: convex gyri and concave sulci. Increasing studies have demonstrated that cortical gyri and sulci exhibit different structural connectivity patterns. Inspired by the intrinsic structural differences between gyri and sulci, in this paper, we present a data-driven framework based on sparse representation of fMRI data for functional network inferences, then examine the interactions within and across gyral and sulcal functional networks and finally elucidate possible functional differences using graph theory based properties. We apply the proposed framework to the high-resolution Human Connectome Project (HCP) grayordinate fMRI data. Extensive experimental results on both resting state fMRI data and task-based fMRI data consistently suggested that gyri are more functionally integrated, while sulci are more functionally segregated in the organizational architecture of cerebral cortex, offering novel understanding of the byzantine cerebral cortex.

Temporal Dynamics Assessment of Spatial Overlap Pattern of Functional Brain Networks Reveals Novel Functional Architecture of Cerebral Cortex

OBJECTIVE

Various studies in the brain mapping field have demonstrated that there exist multiple concurrent functional networks that are spatially overlapped and interacting with each other during specific task performance to jointly realize the total brain function. Assessing such spatial overlap patterns of functional networks (SOPFNs) based on functional magnetic resonance imaging (fMRI) has thus received increasing interest for brain function studies. However, there are still two crucial issues to be addressed. First, the SOPFNs are assessed over the entire fMRI scan assuming the temporal stationarity, while possibly time-dependent dynamics of the SOPFNs is not sufficiently explored. Second, the SOPFNs are assessed within individual subjects, while group-wise consistency of the SOPFNs is largely unknown.

METHODS

To address the two issues, we propose a novel computational framework of group-wise sparse representation of whole-brain fMRI temporal segments to assess the temporal dynamic spatial patterns of SOPFNs that are consistent across different subjects.

RESULTS

Experimental results based on the recently publicly released Human Connectome Project grayordinate task fMRI data demonstrate that meaningful SOPFNs exhibiting dynamic spatial patterns across different time periods are effectively and robustly identified based on the reconstructed concurrent functional networks via the proposed framework. Specifically, those SOPFNs locate significantly more on gyral regions than on sulcal regions across different time periods.

CONCLUSION

These results reveal novel functional architecture of cortical gyri and sulci.

SIGNIFICANCE

Moreover, these results help better understand functional dynamics mechanisms of cerebral cortex in the future.

Functional organization of the fusiform gyrus revealed with connectivity profiles

Within the object recognition-related ventral visual stream, the human fusiform gyrus (FG), which topographically connects the striate cortex to the inferior temporal lobe, plays a pivotal role in high-level visual/cognitive functions. However, though there are many previous investigations of distinct functional modules within the FG, the functional organization of the whole FG in its full functional heterogeneity has not yet been established. In the current study, a replicable functional organization of the FG based on distinct anatomical connectivity patterns was identified. The FG was parcellated into medial (FGm), lateral (FGl), and anterior (FGa) regions using diffusion tensor imaging. We validated the reasonability of such an organizational scheme from the perspective of resting-state whole brain functional connectivity patterns and the involvement of functional subnetworks. We found corroborating support for these three distinct modules, and suggest that the FGm serves as a transition region that combines multiple stimuli, the FGl is responsible for categorical recognition, and the FGa is involved in semantic understanding. These findings support two organizational functional transitions of the ventral temporal gyrus, a posterior/anterior direction of visual/semantic processing, and a media/lateral direction of high-level visual processing. Our results may facilitate a more detailed study of the human FG in the future. Hum Brain Mapp 37:3003-3016, 2016. © 2016 Wiley Periodicals, Inc.

Structural asymmetry of the human cerebral cortex: Regional and between-subject variability of surface area, cortical thickness, and local gyrification

Structural asymmetry varies across individuals, brain regions, and metrics of cortical organization. The current study investigated regional differences in asymmetry of cortical surface area, thickness, and local gyrification, and the extent of between-subject variability in these metrics, in a sample of healthy young adults (N=200). Between-subject variability in cortical structure may provide a means to assess the extent of biological flexibility or constraint of brain regions, and we explored the potential influence of this variability on the phenotypic expression of structural asymmetry. The findings demonstrate that structural asymmetries are nearly ubiquitous across the cortex, with differing regional organization for the three cortical metrics. This implies that there are multiple, only partially overlapping, maps of structural asymmetry. The results further indicate that the degree of asymmetry of a brain region can be predicted by the extent of the region's between-subject variability. These findings provide evidence that reduced biological constraint promotes the expression of strong structural asymmetry.

Neurostructural correlates of consistent and weak handedness

Various cognitive differences have been reported between consistent and weak handers, but little is known about the neurobiological factors that may be associated with this distinction. The current study examined cortical structural lateralization and corpus callosum volume in a large, well-matched sample of young adults (N = 164) to explore potential neurostructural bases for this hand group difference. The groups did not differ in corpus callosum volume. However, at the global hemispheric level, weak handers had reduced or absent asymmetries for grey and white matter volume, cortical surface area, thickness, and local gyrification, relative to consistent handers. Group differences were also observed for some regional hemispheric asymmetries, the most prominent of which was reduced or absent gyrification asymmetry for weak handers in a large region surrounding the central sulcus and extending into parietal association cortex. The findings imply that variations in handedness strength are associated with differences in structural lateralization, not only in somatomotor regions, but also in areas associated with high level cognitive control of action.

Sparse representation of HCP grayordinate data reveals novel functional architecture of cerebral cortex

The recently publicly released Human Connectome Project (HCP) grayordinate-based fMRI data not only has high spatial and temporal resolution, but also offers group-corresponding fMRI signals across a large population for the first time in the brain imaging field, thus significantly facilitating mapping the functional brain architecture with much higher resolution and in a group-wise fashion. In this article, we adopt the HCP grayordinate task-based fMRI (tfMRI) data to systematically identify and characterize task-based heterogeneous functional regions (THFRs) on cortical surface, i.e., the regions that are activated during multiple tasks conditions and contribute to multiple task-evoked systems during a specific task performance, and to assess the spatial patterns of identified THFRs on cortical gyri and sulci by applying a computational framework of sparse representations of grayordinate brain tfMRI signals. Experimental results demonstrate that both consistent task-evoked networks and intrinsic connectivity networks across all subjects and tasks in HCP grayordinate data are effectively and robustly reconstructed via the proposed sparse representation framework. Moreover, it is found that there are relatively consistent THFRs locating at bilateral parietal lobe, frontal lobe, and visual association cortices across all subjects and tasks. Particularly, those identified THFRs locate significantly more on gyral regions than on sulcal regions. These results based on sparse representation of HCP grayordinate data reveal novel functional architecture of cortical gyri and sulci, and might provide a foundation to better understand functional mechanisms of the human cerebral cortex in the future.

Anatomy-guided Dense Individualized and Common Connectivity-based Cortical Landmarks (A-DICCCOL)

Establishment of structural and functional correspondences of human brain that can be quantitatively encoded and reproduced across different subjects and populations is one of the key issues in brain mapping. As an attempt to address this challenge, our recently developed Dense Individualized and Common Connectivity-based Cortical Landmarks (DICCCOL) system reported 358 connectional landmarks, each of which possesses consistent DTI-derived white matter fiber connection pattern that is reproducible in over 240 healthy brains. However, the DICCCOL system can be substantially improved by integrating anatomical and morphological information during landmark initialization and optimization procedures. In this paper, we present a novel anatomy-guided landmark discovery framework that defines and optimizes landmarks via integrating rich anatomical, morphological, and fiber connectional information for landmark initialization, group-wise optimization and prediction, which are formulated and solved as an energy minimization problem. The framework finally determined 555 consistent connectional landmarks. Validation studies demonstrated that the 555 landmarks are reproducible, predictable, and exhibited reasonably accurate anatomical, connectional, and functional correspondences across individuals and populations and thus are named anatomy-guided DICCCOL or A-DICCCOL. This A-DICCCOL system represents common cortical architectures with anatomical, connectional, and functional correspondences across different subjects and would potentially provide opportunities for various applications in brain science.