The reliability and heritability of cortical folds and their genetic correlations across hemispheres

The heritability and structural correlates of resting-state fMRI complexity

The complexity of fMRI signals quantifies temporal dynamics of spontaneous neural activity, which has been increasingly recognized as providing important insights into cognitive functions and psychiatric disorders. However, its heritability and structural underpinnings are not well understood. Here, we utilize multi-scale sample entropy to extract resting-state fMRI complexity in a large healthy adult sample from the Human Connectome Project. We show that fMRI complexity at multiple time scales is heritable in broad brain regions. Heritability estimates are modest and regionally variable. We relate fMRI complexity to brain structure including surface area, cortical myelination, cortical thickness, subcortical volumes, and total brain volume. We find that surface area is negatively correlated with fine-scale complexity and positively correlated with coarse-scale complexity in most cortical regions, especially the association cortex. Most of these correlations are related to common genetic and environmental effects. We also find positive correlations between cortical myelination and fMRI complexity at fine scales and negative correlations at coarse scales in the prefrontal cortex, lateral temporal lobe, precuneus, lateral parietal cortex, and cingulate cortex, with these correlations mainly attributed to common environmental effects. We detect few significant associations between fMRI complexity and cortical thickness. Despite the non-significant association with total brain volume, fMRI complexity exhibits significant correlations with subcortical volumes in the hippocampus, cerebellum, putamen, and pallidum at certain scales. Collectively, our work establishes the genetic basis and structural correlates of resting-state fMRI complexity across multiple scales, supporting its potential application as an endophenotype for psychiatric disorders.

Heritability of functional gradients in the human subcortico-cortical connectivity

The human subcortex plays a pivotal role in cognition and is widely implicated in the pathophysiology of many psychiatric disorders. However, the heritability of functional gradients based on subcortico-cortical functional connectivity remains elusive. Here, leveraging twin functional MRI (fMRI) data from both the Human Connectome Project (n = 1023) and the Adolescent Brain Cognitive Development study (n = 936) datasets, we construct large-scale subcortical functional gradients and delineate an increased principal functional gradient pattern from unimodal sensory/motor networks to transmodal association networks. We observed that this principal functional gradient is heritable, and the strength of heritability exhibits a heterogeneous pattern along a hierarchical unimodal-transmodal axis in subcortex for both young adults and children. Furthermore, employing a machine learning framework, we show that this heterogeneous pattern of the principal functional gradient in subcortex can accurately discern the relationship between monozygotic twin pairs and dizygotic twin pairs with an accuracy of 76.2% (P < 0.001). The heritability of functional gradients is associated with the anatomical myelin proxied by MRI-derived T1-weighted/T2-weighted (T1w/T2w) ratio mapping in subcortex. This study provides new insights into the biological basis of subcortical functional hierarchy by revealing the structural and genetic properties of the subcortical functional gradients.

Phylogenetic differences in the morphology and shape of the central sulcus in great apes and humans: implications for the evolution of motor functions

The central sulcus divides the primary motor and somatosensory cortices in many anthropoid primate brains. Differences exist in the surface area and depth of the central sulcus along the dorso-ventral plane in great apes and humans compared to other primate species. Within hominid species, there are variations in the depth and aspect of their hand motor area, or knob, within the precentral gyrus. In this study, we used post-image analyses on magnetic resonance images to characterize the central sulcus shape of humans, chimpanzees (Pan troglodytes), gorillas (Gorilla gorilla), and orangutans (Pongo pygmaeus and Pongo abelii). Using these data, we examined the morphological variability of central sulcus in hominids, focusing on the hand region, a significant change in human evolution. We show that the central sulcus shape differs between great ape species, but all show similar variations in the location of their hand knob. However, the prevalence of the knob location along the dorso-ventral plane and lateralization differs between species and the presence of a second ventral motor knob seems to be unique to humans. Humans and orangutans exhibit the most similar and complex central sulcus shapes. However, their similarities may reflect divergent evolutionary processes related to selection for different positional and habitual locomotor functions.

Heritability of cerebellar subregion volumes in adolescent and young adult twins

Twin studies have found gross cerebellar volume to be highly heritable. However, whether fine-grained regional volumes within the cerebellum are similarly heritable is still being determined. Anatomical MRI scans from two independent datasets (QTIM: Queensland Twin IMaging, N = 798, mean age 22.1 years; QTAB: Queensland Twin Adolescent Brain, N = 396, mean age 11.3 years) were combined with an optimised and automated cerebellum parcellation algorithm to segment and measure 28 cerebellar regions. We show that the heritability of regional volumetric measures varies widely across the cerebellum (h 2 $$ {h}^2 $$ 47%-91%). Additionally, the good to excellent test-retest reliability for a subsample of QTIM participants suggests that non-genetic variance in cerebellar volumes is due primarily to unique environmental influences rather than measurement error. We also show a consistent pattern of strong associations between the volumes of homologous left and right hemisphere regions. Associations were predominantly driven by genetic effects shared between lobules, with only sparse contributions from environmental effects. These findings are consistent with similar studies of the cerebrum and provide a first approximation of the upper bound of heritability detectable by genome-wide association studies.

Sulcal Morphometry Predicts Mild Cognitive Impairment Conversion to Alzheimer's Disease

Background

Being able to differentiate mild cognitive impairment (MCI) patients who would eventually convert (MCIc) to Alzheimer's disease (AD) from those who would not (MCInc) is a key challenge for prognosis.

Objective

This study aimed to investigate the ability of sulcal morphometry to predict MCI progression to AD, dedicating special attention to an accurate identification of sulci.

Methods

Twenty-five AD patients, thirty-seven MCI and twenty-five healthy controls (HC) underwent a brain-MR protocol (1.5T scanner) including a high-resolution T1-weighted sequence. MCI patients underwent a neuropsychological assessment at baseline and were clinically re-evaluated after a mean of 2.3 years. At follow-up, 12 MCI were classified as MCInc and 25 as MCIc. Sulcal morphometry was investigated using the BrainVISA framework. Consistency of sulci across subjects was ensured by visual inspection and manual correction of the automatic labelling in each subject. Sulcal surface, depth, length, and width were retrieved from 106 sulci. Features were compared across groups and their classification accuracy in predicting MCI conversion was tested. Potential relationships between sulcal features and cognitive scores were explored using Spearman's correlation.

Results

The width of sulci in the temporo-occipital region strongly differentiated between each pair of groups. Comparing MCIc and MCInc, the width of several sulci in the bilateral temporo-occipital and left frontal areas was significantly altered. Higher width of frontal sulci was associated with worse performances in short-term verbal memory and phonemic fluency.

Conclusions

Sulcal morphometry emerged as a strong tool for differentiating HC, MCI, and AD, demonstrating its potential prognostic value for the MCI population.

A bipolar taxonomy of adult human brain sulcal morphology related to timing of fetal sulcation and trans-sulcal gene expression gradients

We developed a computational pipeline (now provided as a resource) for measuring morphological similarity between cortical surface sulci to construct a sulcal phenotype network (SPN) from each magnetic resonance imaging (MRI) scan in an adult cohort (N=34,725; 45-82 years). Networks estimated from pairwise similarities of 40 sulci on 5 morphological metrics comprised two clusters of sulci, represented also by the bipolar distribution of sulci on a linear-to-complex dimension. Linear sulci were more heritable and typically located in unimodal cortex; complex sulci were less heritable and typically located in heteromodal cortex. Aligning these results with an independent fetal brain MRI cohort (N=228; 21-36 gestational weeks), we found that linear sulci formed earlier, and the earliest and latest-forming sulci had the least between-adult variation. Using high-resolution maps of cortical gene expression, we found that linear sulcation is mechanistically underpinned by trans-sulcal gene expression gradients enriched for developmental processes.

Population-wise labeling of sulcal graphs using multi-graph matching

Population-wise matching of the cortical folds is necessary to compute statistics, a required step for e.g. identifying biomarkers of neurological or psychiatric disorders. The difficulty arises from the massive inter-individual variations in the morphology and spatial organization of the folds. The task is challenging both methodologically and conceptually. In the widely used registration-based techniques, these variations are considered as noise and the matching of folds is only implicit. Alternative approaches are based on the extraction and explicit identification of the cortical folds. In particular, representing cortical folding patterns as graphs of sulcal basins-termed sulcal graphs-enables to formalize the task as a graph-matching problem. In this paper, we propose to address the problem of sulcal graph matching directly at the population level using multi-graph matching techniques. First, we motivate the relevance of the multi-graph matching framework in this context. We then present a procedure for generating populations of artificial sulcal graphs, which allows us to benchmark several state-of-the-art multi-graph matching methods. Our results on both artificial and real data demonstrate the effectiveness of multi-graph matching techniques in obtaining a population-wise consistent labeling of cortical folds at the sulcal basin level.

A critical period plasticity framework for the sensorimotor-association axis of cortical neurodevelopment

To understand human brain development it is necessary to describe not only the spatiotemporal patterns of neurodevelopment but also the neurobiological mechanisms that underlie them. Human neuroimaging studies have provided evidence for a hierarchical sensorimotor-to-association (S-A) axis of cortical neurodevelopment. Understanding the biological mechanisms that underlie this program of development using traditional neuroimaging approaches has been challenging. Animal models have been used to identify periods of enhanced experience-dependent plasticity - 'critical periods' - that progress along cortical hierarchies and are governed by a conserved set of neurobiological mechanisms that promote and then restrict plasticity. In this review we hypothesize that the S-A axis of cortical development in humans is partly driven by the cascading maturation of critical period plasticity mechanisms. We then describe how recent advances in in vivo neuroimaging approaches provide a promising path toward testing this hypothesis by linking signals derived from non-invasive imaging to critical period mechanisms.

Temporo-basal sulcal connections: a manual annotation protocol and an investigation of sexual dimorphism and heritability

The temporo-basal region of the human brain is composed of the collateral, the occipito-temporal, and the rhinal sulci. We manually rated (using a novel protocol) the connections between rhinal/collateral (RS-CS), collateral/occipito-temporal (CS-OTS) and rhinal/occipito-temporal (RS-OTS) sulci, using the MRI of nearly 3400 individuals including around 1000 twins. We reported both the associations between sulcal polymorphisms as well with a wide range of demographics (e.g. age, sex, handedness). Finally, we also estimated the heritability, and the genetic correlation between sulcal connections. We reported the frequency of the sulcal connections in the general population, which were hemisphere dependent. We found a sexual dimorphism of the connections, especially marked in the right hemisphere, with a CS-OTS connection more frequent in females (approximately 35-40% versus 20-25% in males) and an RS-CS connection more common in males (approximately 40-45% versus 25-30% in females). We confirmed associations between sulcal connections and characteristics of incomplete hippocampal inversion (IHI). We estimated the broad sense heritability to be 0.28-0.45 for RS-CS and CS-OTS connections, with hints of dominant contribution for the RS-CS connection. The connections appeared to share some of their genetic causing factors as indicated by strong genetic correlations. Heritability appeared much smaller for the (rarer) RS-OTS connection.

Structural neuroimaging phenotypes and associated molecular and genomic underpinnings in autism: a review

Autism has been associated with differences in the developmental trajectories of multiple neuroanatomical features, including cortical thickness, surface area, cortical volume, measures of gyrification, and the gray-white matter tissue contrast. These neuroimaging features have been proposed as intermediate phenotypes on the gradient from genomic variation to behavioral symptoms. Hence, examining what these proxy markers represent, i.e., disentangling their associated molecular and genomic underpinnings, could provide crucial insights into the etiology and pathophysiology of autism. In line with this, an increasing number of studies are exploring the association between neuroanatomical, cellular/molecular, and (epi)genetic variation in autism, both indirectly and directly in vivo and across age. In this review, we aim to summarize the existing literature in autism (and neurotypicals) to chart a putative pathway from (i) imaging-derived neuroanatomical cortical phenotypes to (ii) underlying (neuropathological) biological processes, and (iii) associated genomic variation.

Heritability of cognitive abilities and regional brain structures in middle-aged to elderly East Asians

This study examined the single-nucleotide polymorphism heritability and genetic correlations of cognitive abilities and brain structural measures (regional subcortical volume and cortical thickness) in middle-aged and elderly East Asians (Korean) from the Gwangju Alzheimer's and Related Dementias cohort study. Significant heritability was found in memory function, caudate volume, thickness of the entorhinal cortices, pars opercularis, superior frontal gyri, and transverse temporal gyri. There were 3 significant genetic correlations between (i) the caudate volume and the thickness of the entorhinal cortices, (ii) the thickness of the superior frontal gyri and pars opercularis, and (iii) the thickness of the superior frontal and transverse temporal gyri. This is the first study to describe the heritability and genetic correlations of cognitive and neuroanatomical traits in middle-aged to elderly East Asians. Our results support the previous findings showing that genetic factors play a substantial role in the cognitive and neuroanatomical traits in middle to advanced age. Moreover, by demonstrating shared genetic effects on different brain regions, it gives us a genetic insight into understanding cognitive and brain changes with age, such as aging-related cognitive decline, cortical atrophy, and neural compensation.

Development of cortical folds in the human brain: An attempt to review biological hypotheses, early neuroimaging investigations and functional correlates

The folding of the human brain mostly takes place in utero, making it challenging to study. After a few pioneer studies looking into it in post-mortem foetal specimen, modern approaches based on neuroimaging have allowed the community to investigate the folding process in vivo, its normal progression, its early disturbances, and its relationship to later functional outcomes. In this review article, we aimed to first give an overview of the current hypotheses on the mechanisms governing cortical folding. After describing the methodological difficulties raised by its study in fetuses, neonates and infants with magnetic resonance imaging (MRI), we reported our current understanding of sulcal pattern emergence in the developing brain. We then highlighted the functional relevance of early sulcal development, through recent insights about hemispheric asymmetries and early factors influencing this dynamic such as prematurity. Finally, we outlined how longitudinal studies have started to relate early folding markers and the child's sensorimotor and cognitive outcome. Through this review, we hope to raise awareness on the potential of studying early sulcal patterns both from a fundamental and clinical perspective, as a window into early neurodevelopment and plasticity in relation to growth in utero and postnatal environment of the child.

Genetic determinants of individual variation in the superior temporal sulcus of chimpanzees (Pan troglodytes)

The superior temporal sulcus (STS) is a conserved fold that divides the middle and superior temporal gyri. In humans, there is considerable variation in the shape, folding pattern, lateralization, and depth of the STS that have been reported to be associated with social cognition and linguistic functions. We examined the role that genetic factors play on individual variation in STS morphology in chimpanzees. The surface area and depth of the STS were quantified in sample of 292 captive chimpanzees comprised of two genetically isolated population of individuals. The chimpanzees had been previously genotyped for AVPR1A and KIAA0319, two genes that play a role in social cognition and communication in humans. Single nucleotide polymorphisms in the KIAA0319 and AVPR1A genes were associated with average depth as well as asymmetries in the STS. By contrast, we found no significant effects of these KIA0319 and AVPR1A polymorphism on surface area and depth measures for the central sulcus. The overall findings indicate that genetic factors account for a small to moderate amount of variation in STS morphology in chimpanzees. These findings are discussed in the context of the role of the STS in social cognition and language in humans and their potential evolutionary origins.

Genetic variations within human gained enhancer elements affect human brain sulcal morphology

The expansion of the cerebral cortex is one of the most distinctive changes in the evolution of the human brain. Cortical expansion and related increases in cortical folding may have contributed to emergence of our capacities for high-order cognitive abilities. Molecular analysis of humans, archaic hominins, and non-human primates has allowed identification of chromosomal regions showing evolutionary changes at different points of our phylogenetic history. In this study, we assessed the contributions of genomic annotations spanning 30 million years to human sulcal morphology measured via MRI in more than 18,000 participants from the UK Biobank. We found that variation within brain-expressed human gained enhancers, regulatory genetic elements that emerged since our last common ancestor with Old World monkeys, explained more trait heritability than expected for the left and right calloso-marginal posterior fissures and the right central sulcus. Intriguingly, these are sulci that have been previously linked to the evolution of locomotion in primates and later on bipedalism in our hominin ancestors.

Genetic map of regional sulcal morphology in the human brain from UK biobank data

Genetic associations with macroscopic brain structure can provide insights into brain function and disease. However, specific associations with measures of local brain folding are largely under-explored. Here, we conducted large-scale genome- and exome-wide associations of regional cortical sulcal measures derived from magnetic resonance imaging scans of 40,169 individuals in UK Biobank. We discovered 388 regional brain folding associations across 77 genetic loci, with genes in associated loci enriched for expression in the cerebral cortex, neuronal development processes, and differential regulation during early brain development. We integrated brain eQTLs to refine genes for various loci, implicated several genes involved in neurodevelopmental disorders, and highlighted global genetic correlations with neuropsychiatric phenotypes. We provide an interactive 3D visualisation of our summary associations, emphasising added resolution of regional analyses. Our results offer new insights into the genetic architecture of brain folding and provide a resource for future studies of sulcal morphology in health and disease.

Unfolding the evolution of human cognition

Recent findings spanning fields, from braincases in paleoneurobiology to invivo measurements in cognitive neuroscience, provide insights into the evolution of cognition. Here, we integrate these findings and propose that studying small, evolutionarily new cortical structures has significant implications for identifying new links between neuroanatomical substrates and human-specific aspects of cognition.

Longitudinal Allometry of Sulcal Morphology in Health and Schizophrenia

Scaling between subcomponents of folding and total brain volume (TBV) in healthy individuals (HIs) is allometric. It is unclear whether this is true in schizophrenia (SZ) or first-episode psychosis (FEP). This study confirmed normative allometric scaling norms in HIs using discovery and replication samples. Cross-sectional and longitudinal diagnostic differences in folding subcomponents were then assessed using an allometric framework. Structural imaging from a longitudinal (Sample 1: HI and SZ, nHI Baseline = 298, nSZ Baseline = 169, nHI Follow-up = 293, nSZ Follow-up = 168, totaling 1087 images, all individuals ≥ 2 images, age 16-69 years) and a cross-sectional sample (Sample 2: nHI = 61 and nFEP = 89, age 10-30 years), all human males and females, is leveraged to calculate global folding and its nested subcomponents: sulcation index (SI, total sulcal/cortical hull area) and determinants of sulcal area: sulcal length and sulcal depth. Scaling of SI, sulcal area, and sulcal length with TBV in SZ and FEP was allometric and did not differ from HIs. Longitudinal age trajectories demonstrated steeper loss of SI and sulcal area through adulthood in SZ. Longitudinal allometric analysis revealed that both annual change in SI and sulcal area was significantly stronger related to change in TBV in SZ compared with HIs. Our results detail the first evidence of the disproportionate contribution of changes in SI and sulcal area to TBV changes in SZ. Longitudinal allometric analysis of sulcal morphology provides deeper insight into lifespan trajectories of cortical folding in SZ.SIGNIFICANCE STATEMENT Psychotic disorders are associated with deficits in cortical folding and brain size, but we lack knowledge of how these two morphometric features are related. We leverage cross-sectional and longitudinal samples in which we decompose folding into a set of nested subcomponents: sulcal and hull area, and sulcal depth and length. We reveal that, in both schizophrenia and first-episode psychosis, (1) scaling of subcomponents with brain size is different from expected scaling laws and (2) caution is warranted when interpreting results from traditional methods for brain size correction. Longitudinal allometric scaling points to loss of sulcal area as a principal contributor to loss of brain size in schizophrenia. These findings advance the understanding of cortical folding atypicalities in psychotic disorders.

Test-retest reliability of multi-parametric maps (MPM) of brain microstructure

Multiparameter mapping (MPM) is a quantitative MRI protocol that is promising for studying microstructural brain changes in vivo with high specificity. Reliability values are an important prior knowledge for efficient study design and facilitating replicable findings in development, aging and neuroplasticity research. To explore longitudinal reliability of MPM we acquired the protocol in 31 healthy young subjects twice over a rescan interval of 4 weeks. We assessed the within-subject coefficient of variation (WCV), the between-subject coefficient of variation (BCV), and the intraclass correlation coefficient (ICC). Using these metrics, we investigated the reliability of (semi-) quantitative magnetization transfer saturation (MT_sat), proton density (PD), transversal relaxation (R2*) and longitudinal relaxation (R1). To increase relevance for explorative studies in development and training-induced plasticity, we assess reliability both on local voxel- as well as ROI-level. Finally, we disentangle contributions and interplay of within- and between-subject variability to ICC and assess the optimal degree of spatial smoothing applied to the data. We reveal evidence that voxelwise ICC reliability of MPMs is moderate to good with median values in cortex (subcortical GM): MT: 0.789 (0.447) PD: 0.553 (0.264) R1: 0.555 (0.369) R2*: 0.624 (0.477). The Gaussian smoothing kernel of 2 to 4 mm FWHM resulted in optimal reproducibility. We discuss these findings in the context of longitudinal intervention studies and the application to research designs in neuroimaging field.

Browsing Multiple Subjects When the Atlas Adaptation Cannot Be Achieved via a Warping Strategy

Brain mapping studies often need to identify brain structures or functional circuits into a set of individual brains. To this end, multiple atlases have been published to represent such structures based on different modalities, subject sets, and techniques. The mainstream approach to exploit these atlases consists in spatially deforming each individual data onto a given atlas using dense deformation fields, which supposes the existence of a continuous mapping between atlases and individuals. However, this continuity is not always verified, and this "iconic" approach has limits. We present in this study an alternative, complementary, "structural" approach, which consists in extracting structures from the individual data, and comparing them without deformation. A "structural atlas" is thus a collection of annotated individual data with a common structure nomenclature. It may be used to characterize structure shape variability across individuals or species, or to train machine learning systems. This study exhibits Anatomist, a powerful structural 3D visualization software dedicated to building, exploring, and editing structural atlases involving a large number of subjects. It has been developed primarily to decipher the cortical folding variability; cortical sulci vary enormously in both size and shape, and some may be missing or have various topologies, which makes iconic approaches inefficient to study them. We, therefore, had to build structural atlases for cortical sulci, and use them to train sulci identification algorithms. Anatomist can display multiple subject data in multiple views, supports all kinds of neuroimaging data, including compound structural object graphs, handles arbitrary coordinate transformation chains between data, and has multiple display features. It is designed as a programming library in both C++ and Python languages, and may be extended or used to build dedicated custom applications. Its generic design makes all the display and structural aspects used to explore the variability of the cortical folding pattern work in other applications, for instance, to browse axonal fiber bundles, deep nuclei, functional activations, or other kinds of cortical parcellations. Multimodal, multi-individual, or inter-species display is supported, and adaptations to large scale screen walls have been developed. These very original features make it a unique viewer for structural atlas browsing.

The genetic architecture of human cortical folding

The folding of the human cerebral cortex is a highly genetically regulated process that allows for a much larger surface area to fit into the cranial vault and optimizes functional organization. Sulcal depth is a robust yet understudied measure of localized folding, previously associated with multiple neurodevelopmental disorders. Here, we report the first genome-wide association study of sulcal depth. Through the multivariate omnibus statistical test (MOSTest) applied to vertex-wise measures from 33,748 U.K. Biobank participants (mean age, 64.3 years; 52.0% female), we identified 856 genome-wide significant loci (P < 5 × 10−8). Comparisons with cortical thickness and surface area indicated that sulcal depth has higher locus yield, heritability, and effective sample size. There was a large amount of genetic overlap between these traits, with gene-based analyses indicating strong associations with neurodevelopmental processes. Our findings demonstrate sulcal depth is a promising neuroimaging phenotype that may enhance our understanding of cortical morphology.

Where the genome meets the connectome: Understanding how genes shape human brain connectivity

The integration of modern neuroimaging methods with genetically informative designs and data can shed light on the molecular mechanisms underlying the structural and functional organization of the human connectome. Here, we review studies that have investigated the genetic basis of human brain network structure and function through three complementary frameworks: (1) the quantification of phenotypic heritability through classical twin designs; (2) the identification of specific DNA variants linked to phenotypic variation through association and related studies; and (3) the analysis of correlations between spatial variations in imaging phenotypes and gene expression profiles through the integration of neuroimaging and transcriptional atlas data. We consider the basic foundations, strengths, limitations, and discoveries associated with each approach. We present converging evidence to indicate that anatomical connectivity is under stronger genetic influence than functional connectivity and that genetic influences are not uniformly distributed throughout the brain, with phenotypic variation in certain regions and connections being under stronger genetic control than others. We also consider how the combination of imaging and genetics can be used to understand the ways in which genes may drive brain dysfunction in different clinical disorders.

Towards Deciphering the Fetal Foundation of Normal Cognition and Cognitive Symptoms From Sulcation of the Cortex

Growing evidence supports that prenatal processes play an important role for cognitive ability in normal and clinical conditions. In this context, several neuroimaging studies searched for features in postnatal life that could serve as a proxy for earlier developmental events. A very interesting candidate is the sulcal, or sulco-gyral, patterns, macroscopic features of the cortex anatomy related to the fold topology-e.g., continuous vs. interrupted/broken fold, present vs. absent fold-or their spatial organization. Indeed, as opposed to quantitative features of the cortical sheet (e.g., thickness, surface area or curvature) taking decades to reach the levels measured in adult, the qualitative sulcal patterns are mainly determined before birth and stable across the lifespan. The sulcal patterns therefore offer a window on the fetal constraints on specific brain areas on cognitive abilities and clinical symptoms that manifest later in life. After a global review of the cerebral cortex sulcation, its mechanisms, its ontogenesis along with methodological issues on how to measure the sulcal patterns, we present a selection of studies illustrating that analysis of the sulcal patterns can provide information on prenatal dispositions to cognition (with a focus on cognitive control and academic abilities) and cognitive symptoms (with a focus on schizophrenia and bipolar disorders). Finally, perspectives of sulcal studies are discussed.

Imaging Transcriptomics of Brain Disorders

Noninvasive neuroimaging is a powerful tool for quantifying diverse aspects of brain structure and function in vivo, and it has been used extensively to map the neural changes associated with various brain disorders. However, most neuroimaging techniques offer only indirect measures of underlying pathological mechanisms. The recent development of anatomically comprehensive gene expression atlases has opened new opportunities for studying the transcriptional correlates of noninvasively measured neural phenotypes, offering a rich framework for evaluating pathophysiological hypotheses and putative mechanisms. Here, we provide an overview of some fundamental methods in imaging transcriptomics and outline their application to understanding brain disorders of neurodevelopment, adulthood, and neurodegeneration. Converging evidence indicates that spatial variations in gene expression are linked to normative changes in brain structure during age-related maturation and neurodegeneration that are in part associated with cell-specific gene expression markers of gene expression. Transcriptional correlates of disorder-related neuroimaging phenotypes are also linked to transcriptionally dysregulated genes identified in ex vivo analyses of patient brains. Modeling studies demonstrate that spatial patterns of gene expression are involved in regional vulnerability to neurodegeneration and the spread of disease across the brain. This growing body of work supports the utility of transcriptional atlases in testing hypotheses about the molecular mechanism driving disease-related changes in macroscopic neuroimaging phenotypes.

Sulcation of the intraparietal sulcus is related to symbolic but not non-symbolic number skills

•

The horizontal segment of intraparietal sulcus (HIPS) is one of the key functional regions for processing numbers.

•

Sulcal morphology is a qualitative feature of the brain determined in-utero and not affected by brain maturation and learning.

•

The HIPS sulcal pattern explains part of the variance in participant’s symbolic number comparison and math fluency abilities.

•

Participant’s non-symbolic number comparison abilities was not explained by HIPS sulcal pattern.

•

This association between HIPS sulcal pattern and symbolic number abilities was stable from childhood to young adulthood.

Understanding the constraints, including biological ones, that may influence mathematical development is of great importance because math ability is a key predictor of career success, income and even psychological well-being. While research in developmental cognitive neuroscience of mathematics has extensively studied the key functional regions for processing numbers, particularly the horizontal segment of intraparietal sulcus (HIPS), few studies have investigated the effects of early cerebral constraints on later mathematical abilities. In this pre-registered study, we investigated whether variability of the sulcal pattern of the HIPS, a qualitative feature of the brain determined in-utero and not affected by brain maturation and learning, accounts for individual difference in symbolic and non-symbolic number abilities. Seventy-seven typically developing school-aged children and 21 young adults participated in our study. We found that the HIPS sulcal pattern, (a) explains part of the variance in participant’s symbolic number comparison and math fluency abilities, and (b) that this association between HIPS sulcal pattern and symbolic number abilities was found to be stable from childhood to young adulthood. However, (c) we did not find an association between participant’s non-symbolic number abilities and HIPS sulcal morphology. Our findings suggest that early cerebral constraints may influence individual difference in math abilities, in addition to the well-established neuroplastic factors.

Sex Differences in Lifespan Trajectories and Variability of Human Sulcal and Gyral Morphology

Sex differences in the development and aging of human sulcal morphology have been understudied. We charted sex differences in trajectories and inter-individual variability of global sulcal depth, width, and length, pial surface area, exposed (hull) gyral surface area, unexposed sulcal surface area, cortical thickness, gyral span, and cortex volume across the lifespan in a longitudinal sample (700 scans, 194 participants 2 scans, 104 three scans, age range: 16-70 years) of neurotypical males and females. After adjusting for brain volume, females had thicker cortex and steeper thickness decline until age 40 years; trajectories converged thereafter. Across sexes, sulcal shortening was faster before age 40, while sulcal shallowing and widening were faster thereafter. Although hull area remained stable, sulcal surface area declined and was more strongly associated with sulcal shortening than with sulcal shallowing and widening. Males showed greater variability for cortex volume and lower variability for sulcal width. Our findings highlight the association between loss of sulcal area, notably through sulcal shortening, with cortex volume loss. Studying sex differences in lifespan trajectories may improve knowledge of individual differences in brain development and the pathophysiology of neuropsychiatric conditions.

Individual variation underlying brain age estimates in typical development

Typical brain development follows a protracted trajectory throughout childhood and adolescence. Deviations from typical growth trajectories have been implicated in neurodevelopmental and psychiatric disorders. Recently, the use of machine learning algorithms to model age as a function of structural or functional brain properties has been used to examine advanced or delayed brain maturation in healthy and clinical populations. Termed 'brain age', this approach often relies on complex, nonlinear models that can be difficult to interpret. In this study, we use model explanation methods to examine the cortical features that contribute to brain age modelling on an individual basis. In a large cohort of n = 768 typically-developing children (aged 3-21 years), we build models of brain development using three different machine learning approaches. We employ SHAP, a model-agnostic technique to identify sample-specific feature importance, to identify regional cortical metrics that explain errors in brain age prediction. We find that, on average, brain age prediction and the cortical features that explain model predictions are consistent across model types and reflect previously reported patterns of regions brain development. However, while several regions are found to contribute to brain age prediction error, we find little spatial correspondence between individual estimates of feature importance, even when matched for age, sex and brain age prediction error. We also find no association between brain age error and cognitive performance in this typically-developing sample. Overall, this study shows that, while brain age estimates based on cortical development are relatively robust and consistent across model types and preprocessing strategies, significant between-subject variation exists in the features that explain erroneous brain age predictions on an individual level.

Predicting Progression from Mild Cognitive Impairment to Alzheimer's Disease using MRI-based Cortical Features and a Two-State Markov Model

Magnetic resonance imaging (MRI) has a potential for early diagnosis of individuals at risk for developing Alzheimer's disease (AD). Cognitive performance in healthy elderly people and in those with mild cognitive impairment (MCI) has been associated with measures of cortical gyrification [1] and thickness (CT) [2], yet the extent to which sulcal measures can help to predict AD conversion above and beyond CT measures is not known. Here, we analyzed 721 participants with MCI from phases 1 and 2 of the Alzheimer's Disease Neuroimaging Initiative, applying a two-state Markov model to study the conversion from MCI to AD condition. Our preliminary results suggest that MRI-based cortical features, including sulcal morphometry, may help to predict conversion from MCI to AD.

Intergenerational Transmission of Cortical Sulcal Patterns from Mothers to their Children

Intergenerational effects are described as the genetic, epigenetic, as well as pre- and postnatal environmental influence parents have on their offspring's behavior, cognition, and brain. During fetal brain development, the primary cortical sulci emerge with a distinctive folding pattern that are under strong genetic influence and show little change of this pattern throughout postnatal brain development. We examined intergenerational transmission of cortical sulcal patterns by comparing primary sulcal patterns between children (N = 16, age 5.5 ± 0.81 years, 8 males) and their biological mothers (N = 15, age 39.72 ± 4.68 years) as well as between children and unrelated adult females. Our graph-based sulcal pattern comparison method detected stronger sulcal pattern similarity for child-mother pairs than child-unrelated pairs, where higher similarity between child-mother pairs was observed mostly for the right lobar regions. Our results also show that child-mother versus child-unrelated pairs differ for daughters and sons with a trend toward significance, particularly for the left hemisphere lobar regions. This is the first study to reveal significant intergenerational transmission of cortical sulcal patterns, and our results have important implications for the study of the heritability of complex behaviors, brain-based disorders, the identification of biomarkers, and targets for interventions.

A slower rate of sulcal widening in the brains of the nondemented oldest old

The relationships between aging and brain morphology have been reported in many previous structural brain studies. However, the trajectories of successful brain aging in the extremely old remain underexplored. In the limited research on the oldest old, covering individuals aged 85 years and older, there are very few studies that have focused on the cortical morphology, especially cortical sulcal features. In this paper, we measured sulcal width and depth as well as cortical thickness from T1-weighted scans of 290 nondemented community-dwelling participants aged between 76 and 103 years. We divided the participants into young old (between 76 and 84; mean = 80.35±2.44; male/female = 76/88) and oldest old (between 85 and 103; mean = 91.74±5.11; male/female = 60/66) groups. The results showed that most of the examined sulci significantly widened with increased age and that the rates of sulcal widening were lower in the oldest old. The spatial pattern of the cortical thinning partly corresponded with that of sulcal widening. Compared to females, males had significantly wider sulci, especially in the oldest old. This study builds a foundation for future investigations of neurocognitive disorders and neurodegenerative diseases in the oldest old, including centenarians.

Apolipoprotein E allele 4 effects on Single-Subject Gray Matter Networks in Mild Cognitive Impairment

There is evidence that gray matter networks are disrupted in Mild Cognitive Impairment (MCI) and associated with cognitive impairment and faster disease progression. However, it remains unknown how these alterations are related to the presence of Apolipoprotein E isoform E4 (ApoE4), the most prominent genetic risk factor for late-onset Alzheimer's disease (AD). To investigate this topic at the individual level, we explore the impact of ApoE4 and the disease progression on the Single-Subject Gray Matter Networks (SSGMNets) using the graph theory approach. Our data sample comprised 200 MCI patients selected from the ADNI database, classified as non-Converters and Converters (will progress into AD). Each group included 50 ApoE4-positive ('Carriers', ApoE4 + ) and 50 ApoE4-negative ('non-Carriers', ApoE4-). The SSGMNets were estimated from structural MRIs at two-time points: baseline and conversion. We investigated whether altered network topological measures at baseline and their rate of change (RoC) between baseline and conversion time points were associated with ApoE4 and disease progression. We also explored the correlation of SSGMNets attributes with general cognition score (MMSE), memory (ADNI-MEM), and CSF-derived biomarkers of AD (Aβ42, T-tau, and P-tau). Our results showed that ApoE4 and the disease progression modulated the global topological network properties independently but not in their RoC. MCI converters showed a lower clustering index in several regions associated with neurodegeneration in AD. The SSGMNets' topological organization was revealed to be able to predict cognitive and memory measures. The findings presented here suggest that SSGMNets could indeed be used to identify MCI ApoE4 Carriers with a high risk for AD progression.

Evidence in cortical folding patterns for prenatal predispositions to hallucinations in schizophrenia

All perception is a construction of the brain from sensory input. Our first perceptions begin during gestation, making fetal brain development fundamental to how we experience a diverse world. Hallucinations are percepts without origin in physical reality that occur in health and disease. Despite longstanding research on the brain structures supporting hallucinations and on perinatal contributions to the pathophysiology of schizophrenia, what links these two distinct lines of research remains unclear. Sulcal patterns derived from structural magnetic resonance (MR) images can provide a proxy in adulthood for early brain development. We studied two independent datasets of patients with schizophrenia who underwent clinical assessment and 3T MR imaging from the United Kingdom and Shanghai, China (n = 181 combined) and 63 healthy controls from Shanghai. Participants were stratified into those with (n = 79 UK; n = 22 Shanghai) and without (n = 43 UK; n = 37 Shanghai) hallucinations from the PANSS P3 scores for hallucinatory behaviour. We quantified the length, depth, and asymmetry indices of the paracingulate and superior temporal sulci (PCS, STS), which have previously been associated with hallucinations in schizophrenia, and constructed cortical folding covariance matrices organized by large-scale functional networks. In both ethnic groups, we demonstrated a significantly shorter left PCS in patients with hallucinations compared to those without, and to healthy controls. Reduced PCS length and STS depth corresponded to focal deviations in their geometry and to significantly increased covariance within and between areas of the salience and auditory networks. The discovery of neurodevelopmental alterations contributing to hallucinations establishes testable models for these enigmatic, sometimes highly distressing, perceptions and provides mechanistic insight into the pathological consequences of prenatal origins.

Dissimilarity in Sulcal Width Patterns in the Cortex can be Used to Identify Patients With Schizophrenia With Extreme Deficits in Cognitive Performance

Schizophrenia is a biologically complex disorder with multiple regional deficits in cortical brain morphology. In addition, interindividual heterogeneity of cortical morphological metrics is larger in patients with schizophrenia when compared to healthy controls. Exploiting interindividual differences in the severity of cortical morphological deficits in patients instead of focusing on group averages may aid in detecting biologically informed homogeneous subgroups. The person-based similarity index (PBSI) of brain morphology indexes an individual's morphometric similarity across numerous cortical regions amongst a sample of healthy subjects. We extended the PBSI such that it indexes the morphometric similarity of an independent individual (eg, a patient) with respect to healthy control subjects. By employing a normative modeling approach on longitudinal data, we determined an individual's degree of morphometric dissimilarity to the norm. We calculated the PBSI for sulcal width (PBSI-SW) in patients with schizophrenia and healthy control subjects (164 patients and 164 healthy controls; 656 magnetic resonance imaging scans) and associated it with cognitive performance and cortical sulcation index. A subgroup of patients with markedly deviant PBSI-SW showed extreme deficits in cognitive performance and cortical sulcation. Progressive reduction of PBSI-SW in the schizophrenia group relative to healthy controls was driven by these deviating individuals. By explicitly leveraging interindividual differences in the severity of PBSI-SW deficits, neuroimaging-driven subgrouping of patients is feasible. As such, our results pave the way for future applications of morphometric similarity indices for subtyping of clinical populations.

The Heritability of Cortical Folding: Evidence from the Human Connectome Project

The mechanisms underlying cortical folding are incompletely understood. Prior studies have suggested that individual differences in sulcal depth are genetically mediated, with deeper and ontologically older sulci more heritable than others. In this study, we examine FreeSurfer-derived estimates of average convexity and mean curvature as proxy measures of cortical folding patterns using a large (N = 1096) genetically informative young adult subsample of the Human Connectome Project. Both measures were significantly heritable near major sulci and primary fissures, where approximately half of individual differences could be attributed to genetic factors. Genetic influences near higher order gyri and sulci were substantially lower and largely nonsignificant. Spatial permutation analysis found that heritability patterns were significantly anticorrelated to maps of evolutionary and neurodevelopmental expansion. We also found strong phenotypic correlations between average convexity, curvature, and several common surface metrics (cortical thickness, surface area, and cortical myelination). However, quantitative genetic models suggest that correlations between these metrics are largely driven by nongenetic factors. These findings not only further our understanding of the neurobiology of gyrification, but have pragmatic implications for the interpretation of heritability maps based on automated surface-based measurements.