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

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.

Anatomical and behavioral correlates of auditory perception in developmental dyslexia

Developmental dyslexia is typically associated with difficulties in basic auditory processing and in manipulating speech sounds. However, the neuroanatomical correlates of auditory difficulties in developmental dyslexia (DD) and their contribution to individual clinical phenotypes are still unknown. Recent intracranial electrocorticography findings associated processing of sound amplitude rises and speech sounds with posterior and middle superior temporal gyrus (STG), respectively. We hypothesize that regional STG anatomy will relate to specific auditory abilities in DD, and that auditory processing abilities will relate to behavioral difficulties with speech and reading. One hundred and ten children (78 DD, 32 typically developing, age 7-15 years) completed amplitude rise time and speech in noise discrimination tasks. They also underwent a battery of cognitive tests. Anatomical MRI scans were used to identify regions in which local cortical gyrification complexity correlated with auditory behavior. Behaviorally, amplitude rise time but not speech in noise performance was impaired in DD. Neurally, amplitude rise time and speech in noise performance correlated with gyrification in posterior and middle STG, respectively. Furthermore, amplitude rise time significantly contributed to reading impairments in DD, while speech in noise only explained variance in phonological awareness. Finally, amplitude rise time and speech in noise performance were not correlated, and each task was correlated with distinct neuropsychological measures, emphasizing their unique contributions to DD. Overall, we provide a direct link between the neurodevelopment of the left STG and individual variability in auditory processing abilities in neurotypical and dyslexic populations.

Species -shared and -unique gyral peaks on human and macaque brains

Cortical folding is an important feature of primate brains that plays a crucial role in various cognitive and behavioral processes. Extensive research has revealed both similarities and differences in folding morphology and brain function among primates including macaque and human. The folding morphology is the basis of brain function, making cross-species studies on folding morphology important for understanding brain function and species evolution. However, prior studies on cross-species folding morphology mainly focused on partial regions of the cortex instead of the entire brain. Previously, our research defined a whole-brain landmark based on folding morphology: the gyral peak. It was found to exist stably across individuals and ages in both human and macaque brains. Shared and unique gyral peaks in human and macaque are identified in this study, and their similarities and differences in spatial distribution, anatomical morphology, and functional connectivity were also dicussed.

Species -Shared and -Unique Gyral Peaks on Human and Macaque Brains

Cortical folding is an important feature of primate brains that plays a crucial role in various cognitive and behavioral processes. Extensive research has revealed both similarities and differences in folding morphology and brain function among primates including macaque and human. The folding morphology is the basis of brain function, making cross-species studies on folding morphology important for understanding brain function and species evolution. However, prior studies on cross-species folding morphology mainly focused on partial regions of the cortex instead of the entire brain. Previously, we defined a whole-brain landmark based on folding morphology: the gyral peak. It was found to exist stably across individuals and ages in both human and macaque brains. In this study, we identified shared and unique gyral peaks in human and macaque, and investigated the similarities and differences in the spatial distribution, anatomical morphology, and functional connectivity of them.

Shaping the brain: The emergence of cortical structure and folding

The cerebral cortex-the brain's covering and largest region-has increased in size and complexity in humans and supports higher cognitive functions such as language and abstract thinking. There is a growing understanding of the human cerebral cortex, including the diversity and number of cell types that it contains, as well as of the developmental mechanisms that shape cortical structure and organization. In this review, we discuss recent progress in our understanding of molecular and cellular processes, as well as mechanical forces, that regulate the folding of the cerebral cortex. Advances in human genetics, coupled with experimental modeling in gyrencephalic species, have provided insights into the central role of cortical progenitors in the gyrification and evolutionary expansion of the cerebral cortex. These studies are essential for understanding the emergence of structural and functional organization during cortical development and the pathogenesis of neurodevelopmental disorders associated with cortical malformations.

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.

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.

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.

Multivariate Concavity Amplitude Index (MCAI) for characterizing Heschl's gyrus shape

Heschl's gyrus (HG), which includes primary auditory cortex, is highly variable in its shape (i.e. gyrification patterns), between hemispheres and across individuals. Differences in HG shape have been observed in the context of phonetic learning skill and expertise, and of professional musicianship, among others. Two of the most common configurations of HG include single HG, where a single transverse temporal gyrus is present, and common stem duplications (CSD), where a sulcus intermedius (SI) arises from the lateral aspect of HG. Here we describe a new toolbox, called 'Multivariate Concavity Amplitude Index' (MCAI), which automatically assesses the shape of HG. MCAI works on the output of TASH, our first toolbox which automatically segments HG, and computes continuous indices of concavity, which arise when sulci are present, along the outer perimeter of an inflated representation of HG, in a directional manner. Thus, MCAI provides a multivariate measure of shape, which is reproducible and sensitive to small variations in shape. We applied MCAI to structural magnetic resonance imaging (MRI) data of N=181 participants, including professional and amateur musicians and from non-musicians. Former studies have shown large variations in HG shape in the former groups. We validated MCAI by showing high correlations between the dominant (i.e. highest) lateral concavity values and continuous visual assessments of the degree of lateral gyrification of the first gyrus. As an application of MCAI, we also replicated previous visually obtained findings showing a higher likelihood of bilateral CSDs in musicians. MCAI opens a wide range of applications in evaluating HG shape in the context of individual differences, expertise, disorder and genetics.

Testing the Ecophenotype Model: Cortical Structure Alterations in Conduct Disorder With Versus Without Childhood Maltreatment

BACKGROUND

Childhood maltreatment is common in youths with conduct disorder (CD), and both CD and maltreatment have been linked to neuroanatomical alterations. Nonetheless, our understanding of the contribution of maltreatment to the neuroanatomical alterations observed in CD remains limited. We tested the applicability of the ecophenotype model to CD, which holds that maltreatment-related psychopathology is (neurobiologically) distinct from psychopathology without maltreatment.

METHODS

Surface-based morphometry was used to investigate cortical volume, thickness, surface area, and gyrification in a mixed-sex sample of participants with CD (n = 114) and healthy control subjects (HCs) (n = 146), ages 9 to 18 years. Using vertexwise general linear models adjusted for sex, age, total intracranial volume, and site, the control group was compared with the overall CD group and the CD subgroups with (n = 49) versus without (n = 65) maltreatment (assessed by the Children's Bad Experiences interview). These subgroups were also directly compared.

RESULTS

The overall CD group showed lower cortical thickness in the right inferior frontal gyrus. CD youths with a history of maltreatment showed more widespread structural alterations relative to HCs, comprising lower thickness, volume, and gyrification in inferior and middle frontal regions. Conversely, CD youths with no history of maltreatment only showed greater left superior temporal gyrus folding relative to HCs. When contrasting the CD subgroups, those with maltreatment displayed lower right superior temporal gyrus volume, right precentral gyrus surface area, and gyrification in frontal, temporal, and parietal regions.

CONCLUSIONS

Consistent with the ecophenotype model, findings indicated that CD youths with versus without maltreatment differ neurobiologically. This highlights the importance of considering maltreatment history in neuroimaging studies of CD and other disorders.

Genetic and Environmental Variation in Continuous Phenotypes in the ABCD Study®

Twin studies yield valuable insights into the sources of variation, covariation and causation in human traits. The ABCD Study® (abcdstudy.org) was designed to take advantage of four universities known for their twin research, neuroimaging, population-based sampling, and expertise in genetic epidemiology so that representative twin studies could be performed. In this paper we use the twin data to: (i) provide initial estimates of heritability for the wide range of phenotypes assessed in the ABCD Study using a consistent direct variance estimation approach, assuring that both data and methodology are sound; and (ii) provide an online resource for researchers that can serve as a reference point for future behavior genetic studies of this publicly available dataset. Data were analyzed from 772 pairs of twins aged 9-10 years at study inception, with zygosity determined using genotypic data, recruited and assessed at four twin hub sites. The online tool provides twin correlations and both standardized and unstandardized estimates of additive genetic, and environmental variation for 14,500 continuously distributed phenotypic features, including: structural and functional neuroimaging, neurocognition, personality, psychopathology, substance use propensity, physical, and environmental trait variables. The estimates were obtained using an unconstrained variance approach, so they can be incorporated directly into meta-analyses without upwardly biasing aggregate estimates. The results indicated broad consistency with prior literature where available and provided novel estimates for phenotypes without prior twin studies or those assessed at different ages. Effects of site, self-identified race/ethnicity, age and sex were statistically controlled. Results from genetic modeling of all 53,172 continuous variables, including 38,672 functional MRI variables, will be accessible via the user-friendly open-access web interface we have established, and will be updated as new data are released from the ABCD Study. This paper provides an overview of the initial results from the twin study embedded within the ABCD Study, an introduction to the primary research domains in the ABCD study and twin methodology, and an evaluation of the initial findings with a focus on data quality and suitability for future behavior genetic studies using the ABCD dataset. The broad introductory material is provided in recognition of the multidisciplinary appeal of the ABCD Study. While this paper focuses on univariate analyses, we emphasize the opportunities for multivariate, developmental and causal analyses, as well as those evaluating heterogeneity by key moderators such as sex, demographic factors and genetic background.

Assessment of Characteristics of Imaging Biomarkers for Quantifying Anterior Cingulate Cortex Changes: A Twin Study of Middle- to Advanced-Aged Populations in East Asia

Background and Objectives: Our aim was to assess genetic and environmental effects on surface morphological parameters for quantifying anterior cingulate cortex (ACC) changes in middle- to advanced-age East Asians using twin analysis. Materials and Methods: Normal twins over 39 years old comprising 37 monozygotic pairs and 17 dizygotic pairs underwent 3-dimensional (3D) T1-weighted imaging of the brain at 3T. Freesurfer-derived ACC parameters including thickness, standard deviation of thickness (STDthickness), volume, surface area, and sulcal morphological parameters (folding, mean, and Gaussian curvatures) were calculated from 3D T1-weighted volume images. Twin analysis with a model involving phenotype variance components of additive genetic effects (A), common environmental effects (C), and unique environmental effects (E) was performed to assess the magnitude of each genetic and environmental influence on parameters. Results: Most parameters fit best with an AE model. Both thickness (A: left 0.73/right 0.71) and surface area (A: left 0.63/right 0.71) were highly heritable. STDthickness was low to moderately heritable (A: left 0.48/right 0.29). Volume was moderately heritable (A: left 0.37). Folding was low to moderately heritable (A: left 0.44/right 0.28). Mean curvature (A: left 0.37/right 0.65) and Gaussian curvature (A: right 0.79) were moderately to highly heritable. Right volume and left Gaussian curvature fit best with a CE model, indicating a relatively weak contribution of genetic factors to these parameters. Conclusions: When assessing ACC changes in middle- to advanced-age East Asians, one must keep in mind that thickness and surface area appear to be strongly affected by genetic factors, whereas sulcal morphological parameters tend to involve environmental factors.

Mapping the genetic architecture of cortical morphology through neuroimaging: progress and perspectives

Cortical morphology is a key determinant of cognitive ability and mental health. Its development is a highly intricate process spanning decades, involving the coordinated, localized expression of thousands of genes. We are now beginning to unravel the genetic architecture of cortical morphology, thanks to the recent availability of large-scale neuroimaging and genomic data and the development of powerful biostatistical tools. Here, we review the progress made in this field, providing an overview of the lessons learned from genetic studies of cortical volume, thickness, surface area, and folding as captured by neuroimaging. It is now clear that morphology is shaped by thousands of genetic variants, with effects that are region- and time-dependent, thereby challenging conventional study approaches. The most recent genome-wide association studies have started discovering common genetic variants influencing cortical thickness and surface area, yet together these explain only a fraction of the high heritability of these measures. Further, the impact of rare variants and non-additive effects remains elusive. There are indications that the quickly increasing availability of data from whole-genome sequencing and large, deeply phenotyped population cohorts across the lifespan will enable us to uncover much of the missing heritability in the upcoming years. Novel approaches leveraging shared information across measures will accelerate this process by providing substantial increases in statistical power, together with more accurate mapping of genetic relationships. Important challenges remain, including better representation of understudied demographic groups, integration of other 'omics data, and mapping of effects from gene to brain to behavior across the lifespan.

Variability and concordance of sulcal patterns in the orbitofrontal cortex: A twin study

Sulcogyral patterns have been identified in the orbitofrontal cortex (OFC) based on the continuity of the medial and lateral orbital sulci. Pattern types are named according to their frequency in the population, with Type I present in ∼60%, Type II in ∼25%, Type III in ∼10%, and Type IV in ∼5%. Previous work has demonstrated that psychiatric conditions with high estimated heritability (e.g. schizophrenia, bipolar disorder) are associated with reduced frequency of Type I patterns, but the general heritability of the OFC sulcogyral patterns is unknown. We examined concordance of OFC patterns in 304 monozygotic (MZ) twins relative to 172 dizygotic (DZ) twins using structural magnetic resonance imaging data. We find that the frequency of pattern types within MZ and DZ twins are similar and bilateral concordance rates across all pattern types in DZ twins were 14% and 21% for MZ twins. Results from follow-up analyses confirm that continuity in the rostral-caudal direction is an important source of variability within the OFC, and subtype analyses indicate that variability is present in other sulci that are not represented by overall OFC pattern type. Overall, these results suggest that OFC sulcogyral patterns may reflect important variance that is not genetic in origin.

Machine learning in neuroimaging: from research to clinical practice

Neuroimaging is critical in clinical care and research, enabling us to investigate the brain in health and disease. There is a complex link between the brain's morphological structure, physiological architecture, and the corresponding imaging characteristics. The shape, function, and relationships between various brain areas change during development and throughout life, disease, and recovery. Like few other areas, neuroimaging benefits from advanced analysis techniques to fully exploit imaging data for studying the brain and its function. Recently, machine learning has started to contribute (a) to anatomical measurements, detection, segmentation, and quantification of lesions and disease patterns, (b) to the rapid identification of acute conditions such as stroke, or (c) to the tracking of imaging changes over time. As our ability to image and analyze the brain advances, so does our understanding of its intricate relationships and their role in therapeutic decision-making. Here, we review the current state of the art in using machine learning techniques to exploit neuroimaging data for clinical care and research, providing an overview of clinical applications and their contribution to fundamental computational neuroscience.

Computational models of cortical folding: A review of common approaches

The process of gyrification, by which the brain develops the intricate pattern of gyral hills and sulcal valleys, is the result of interactions between biological and mechanical processes during brain development. Researchers have developed a vast array of computational models in order to investigate cortical folding. This review aims to summarize these studies, focusing on five essential elements of the brain that affect development and gyrification and how they are represented in computational models: (i) the constraints of skull, meninges, and cerebrospinal fluid; (ii) heterogeneity of cortical layers and regions; (iii) anisotropic behavior of subcortical fiber tracts; (iv) material properties of brain tissue; and (v) the complex geometry of the brain. Finally, we highlight areas of need for future simulations of brain development.

A comparative study of pre-alpha islands in the entorhinal cortex from selected primates and in lissencephaly

The entorhinal cortex (EC) is the main interface between the sensory association areas of the neocortex and the hippocampus. It is crucial for the evaluation and processing of sensory data for long-term memory consolidation, and shows damage in many brain diseases, e.g., neurodegenerative diseases, such as Alzheimer's disease and developmental disorders. The pre-alpha layer of the EC in humans (layer II) displays a remarkable distribution of neurons in islands. These cellular islands give rise to a portion of the perforant path - the major reciprocal data stream for neocortical information into the hippocampal formation. However, the functional relevance of the morphological appearance of the pre-alpha layer in cellular islands and the precise timing of their initial appearance during primate evolution are largely unknown. Here, we conducted a comparative study of the EC from 38 non-human primates and Homo sapiens and found a strong relationship between gyrification index (GI) and the presence of the pre-alpha cellular islands. The formation of cellular islands also correlated wih brain and body weight as well as neopallial volume. In the two human lissencephalic cases, the cellular islands in the pre-alpha layer were lacking. These findings emphasize the relationship between cortical folding and island formation in the entorhinal cortex from an evolutionary perspective, and suggest a role in the pathomechanism of developmental brain disorders. This article is protected by copyright. All rights reserved.

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.