Imaging local genetic influences on cortical folding

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.

The longitudinal impact of type 2 diabetes on brain gyrification

Type 2 diabetes has an effect on brain structure, including cortical gyrification. The significance of these changes is better understood if assessed over time. However, there is a lack of studies assessing longitudinally the effect of this disease with complex aethology in gyrification. While changes in this feature have been associated mainly with genetic legacy, our study allowed to shed light on the effect of the variation of glycaemic profile over time in gyrification in this metabolic disease. In this longitudinal study, we analysed brain anatomical magnetic resonance images of 15 participants with type 2 diabetes and 13 healthy control participants to investigate the impact of this metabolic disease on the gyrification index over a 7-year period. We observed a significant interaction between time and group in six regions, four of which (left precentral gyrus, left gyrus rectus, left subcentral gyrus and sulci and right inferior temporal gyrus) showed an increase in gyrification in type 2 diabetes and a decrease in the control group and the two others (left pericallosal sulcus and right inferior frontal sulcus) the opposite pattern. The variation of the gyrification was correlated with the variation of the glycaemic profile. Following the interaction, the simple main effect of time in each group separately has shown that in the group with diabetes, there were more regions susceptible to alterations of gyrification. In sum, our results raise credit for the possibility that glycaemic control also might influence gyrification in type 2 diabetes.

Decreased cortical gyrification in major depressive disorder

BACKGROUND

Early neurodevelopmental deviations, such as abnormal cortical folding patterns, are candidate biomarkers of major depressive disorder (MDD). We aimed to investigate the association of MDD with the local gyrification index (LGI) in each cortical region at the whole-brain level, and the association of the LGI with clinical characteristics of MDD.

METHODS

We obtained T1-weighted images from 234 patients with MDD and 215 healthy controls (HCs). The LGI values from 66 cortical regions in the bilateral hemispheres were automatically calculated according to the Desikan-Killiany atlas. We compared the LGI values between the MDD and HC groups using analysis of covariance, including age, sex, and years of education as covariates. The association between the clinical characteristics and LGI values was investigated in the MDD group.

RESULTS

Compared with HCs, patients with MDD showed significantly decreased LGI values in the cortical regions, including the bilateral ventrolateral and dorsolateral prefrontal cortices, medial and lateral orbitofrontal cortices, insula, right rostral anterior cingulate cortex, and several temporal and parietal regions, with the largest effect size in the left pars triangularis (Cohen's f2 = 0.361; p = 1.78 × 10-13). Regarding the association of clinical characteristics with LGIs within the MDD group, recurrence and longer illness duration were associated with increased gyrification in several occipital and temporal regions, which showed no significant difference in LGIs between the MDD and HC groups.

CONCLUSIONS

These findings suggest that the LGI may be a relatively stable neuroimaging marker associated with MDD predisposition.

A comparison of structural morphometry in children and adults with persistent developmental stuttering

This cross-sectional study aimed to differentiate earlier occurring neuroanatomical differences that may reflect core deficits in stuttering versus changes associated with a longer duration of stuttering by analysing structural morphometry in a large sample of children and adults who stutter and age-matched controls. Whole-brain T1-weighted structural scans were obtained from 166 individuals who stutter (74 children, 92 adults; ages 3-58) and 191 controls (92 children, 99 adults; ages 3-53) from eight prior studies in our laboratories. Mean size and gyrification measures were extracted using FreeSurfer software for each cortical region of interest. FreeSurfer software was also used to generate subcortical volumes for regions in the automatic subcortical segmentation. For cortical analyses, separate ANOVA analyses of size (surface area, cortical thickness) and gyrification (local gyrification index) measures were conducted to test for a main effect of diagnosis (stuttering, control) and the interaction of diagnosis-group with age-group (children, adults) across cortical regions. Cortical analyses were first conducted across a set of regions that comprise the speech network and then in a second whole-brain analysis. Next, separate ANOVA analyses of volume were conducted across subcortical regions in each hemisphere. False discovery rate corrections were applied for all analyses. Additionally, we tested for correlations between structural morphometry and stuttering severity. Analyses revealed thinner cortex in children who stutter compared with controls in several key speech-planning regions, with significant correlations between cortical thickness and stuttering severity. These differences in cortical size were not present in adults who stutter, who instead showed reduced gyrification in the right inferior frontal gyrus. Findings suggest that early cortical anomalies in key speech planning regions may be associated with stuttering onset. Persistent stuttering into adulthood may result from network-level dysfunction instead of focal differences in cortical morphometry. Adults who stutter may also have a more heterogeneous neural presentation than children who stutter due to their unique lived experiences.

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.

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.

ABLE: Automated Brain Lines Extraction Based on Laplacian Surface Collapse

The archetypical folded shape of the human cortex has been a long-standing topic for neuroscientific research. Nevertheless, the accurate neuroanatomical segmentation of sulci remains a challenge. Part of the problem is the uncertainty of where a sulcus transitions into a gyrus and vice versa. This problem can be avoided by focusing on sulcal fundi and gyral crowns, which represent the topological opposites of cortical folding. We present Automated Brain Lines Extraction (ABLE), a method based on Laplacian surface collapse to reliably segment sulcal fundi and gyral crown lines. ABLE is built to work on standard FreeSurfer outputs and eludes the delineation of anastomotic sulci while maintaining sulcal fundi lines that traverse the regions with the highest depth and curvature. First, it segments the cortex into gyral and sulcal surfaces; then, each surface is spatially filtered. A Laplacian-collapse-based algorithm is applied to obtain a thinned representation of the surfaces. This surface is then used for careful detection of the endpoints of the lines. Finally, sulcal fundi and gyral crown lines are obtained by eroding the surfaces while preserving the connectivity between the endpoints. The method is validated by comparing ABLE with three other sulcal extraction methods using the Human Connectome Project (HCP) test-retest database to assess the reproducibility of the different tools. The results confirm ABLE as a reliable method for obtaining sulcal lines with an accurate representation of the sulcal topology while ignoring anastomotic branches and the overestimation of the sulcal fundi lines. ABLE is publicly available via https://github.com/HGGM-LIM/ABLE .

Alteration of a brain network with stable and strong functional connections in subjects with schizophrenia

It is widely accepted that there are some common network patterns in the human brain. However, the existence of stable and strong functional connections in the human brain and whether they change in schizophrenia is still a question. By setting 1% connections with the smallest coefficient of variation, we found a widespread brain functional network (frame network) in healthy people(n = 380, two datasets from public databases). We then explored the alterations in a medicated group (60 subjects with schizophrenia vs 71 matched controls) and a drug-naive first-episode group (68 subjects with schizophrenia vs 45 matched controls). A linear support vector classifier (SVC) was constructed to distinguish patients and controls using the medicated patients' frame network. We found most frame connections of healthy people had high strength, which were symmetrical and connected the left and right hemispheres. Conversely, significant differences in frame connections were observed in both patient groups, which were positively correlated with negative symptoms (mainly language dysfunction). Additionally, patients' frame network were more left-lateralized, concentrating on the left frontal lobe, and was quite accurate at distinguishing medicated patients from controls (classifier accuracy was 78.63%, sensitivity was 86.67%, specificity was 76.06%, and the area under the curve (AUC) was 0.83). Furthermore, the results were repeated in the drug-naive set (accuracy was 84.96%, sensitivity was 85.29%, specificity was 88.89%, and AUC was 0.93). These findings indicate that the abnormal pattern of frame network in subjects with schizophrenia might provide new insights into the dysconnectivity in schizophrenia.

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.

Altered Cerebral Curvature in Preterm Infants Is Associated with the Common Genetic Variation Related to Autism Spectrum Disorder and Lipid Metabolism

Preterm births are often associated with neurodevelopmental impairment. In the critical developmental period of the fetal brain, preterm birth disrupts cortical maturation. Notably, preterm birth leads to alterations in the fronto-striatal and temporal lobes and the limbic region. Recent advances in MRI acquisition and analysis methods have revealed an integrated approach to the genetic influence on brain structure. Based on imaging studies, we hypothesized that the altered cortical structure observed after preterm birth is associated with common genetic variations. We found that the presence of the minor allele at rs1042778 in OXTR was associated with reduced curvature in the right medial orbitofrontal gyrus (p < 0.001). The presence of the minor allele at rs174576 in FADS2 (p < 0.001) or rs740603 in COMT (p < 0.001) was related to reduced curvature in the left posterior cingulate gyrus. This study provides biological insight into altered cortical curvature at term-equivalent age, suggesting that the common genetic variations related to autism spectrum disorder (ASD) and lipid metabolism may mediate vulnerability to early cortical dysmaturation in preterm infants.

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.

Optimal Method for Fetal Brain Age Prediction Using Multiplanar Slices From Structural Magnetic Resonance Imaging

The accurate prediction of fetal brain age using magnetic resonance imaging (MRI) may contribute to the identification of brain abnormalities and the risk of adverse developmental outcomes. This study aimed to propose a method for predicting fetal brain age using MRIs from 220 healthy fetuses between 15.9 and 38.7 weeks of gestational age (GA). We built a 2D single-channel convolutional neural network (CNN) with multiplanar MRI slices in different orthogonal planes without correction for interslice motion. In each fetus, multiple age predictions from different slices were generated, and the brain age was obtained using the mode that determined the most frequent value among the multiple predictions from the 2D single-channel CNN. We obtained a mean absolute error (MAE) of 0.125 weeks (0.875 days) between the GA and brain age across the fetuses. The use of multiplanar slices achieved significantly lower prediction error and its variance than the use of a single slice and a single MRI stack. Our 2D single-channel CNN with multiplanar slices yielded a significantly lower stack-wise MAE (0.304 weeks) than the 2D multi-channel (MAE = 0.979, p < 0.001) and 3D (MAE = 1.114, p < 0.001) CNNs. The saliency maps from our method indicated that the anatomical information describing the cortex and ventricles was the primary contributor to brain age prediction. With the application of the proposed method to external MRIs from 21 healthy fetuses, we obtained an MAE of 0.508 weeks. Based on the external MRIs, we found that the stack-wise MAE of the 2D single-channel CNN (0.743 weeks) was significantly lower than those of the 2D multi-channel (1.466 weeks, p < 0.001) and 3D (1.241 weeks, p < 0.001) CNNs. These results demonstrate that our method with multiplanar slices accurately predicts fetal brain age without the need for increased dimensionality or complex MRI preprocessing steps.

Neurodevelopment of the association cortices: Patterns, mechanisms, and implications for psychopathology

The human brain undergoes a prolonged period of cortical development that spans multiple decades. During childhood and adolescence, cortical development progresses from lower-order, primary and unimodal cortices with sensory and motor functions to higher-order, transmodal association cortices subserving executive, socioemotional, and mentalizing functions. The spatiotemporal patterning of cortical maturation thus proceeds in a hierarchical manner, conforming to an evolutionarily rooted, sensorimotor-to-association axis of cortical organization. This developmental program has been characterized by data derived from multimodal human neuroimaging and is linked to the hierarchical unfolding of plasticity-related neurobiological events. Critically, this developmental program serves to enhance feature variation between lower-order and higher-order regions, thus endowing the brain's association cortices with unique functional properties. However, accumulating evidence suggests that protracted plasticity within late-maturing association cortices, which represents a defining feature of the human developmental program, also confers risk for diverse developmental psychopathologies.

Mapping gene transcription and neurocognition across human neocortex

Regulation of gene expression drives protein interactions that govern synaptic wiring and neuronal activity. The resulting coordinated activity among neuronal populations supports complex psychological processes, yet how gene expression shapes cognition and emotion remains unknown. Here, we directly bridge the microscale and macroscale by mapping gene expression patterns to functional activation patterns across the cortical sheet. Applying unsupervised learning to the Allen Human Brain Atlas and Neurosynth databases, we identify a ventromedial-dorsolateral gradient of gene assemblies that separate affective and perceptual domains. This topographic molecular-psychological signature reflects the hierarchical organization of the neocortex, including systematic variations in cell type, myeloarchitecture, laminar differentiation and intrinsic network affiliation. In addition, this molecular-psychological signature strengthens over neurodevelopment and can be replicated in two independent repositories. Collectively, our results reveal spatially covarying transcriptomic and cognitive architectures, highlighting the influence that molecular mechanisms exert on psychological processes.

A novel morphometric signature of brain alterations in type 2 diabetes: Patterns of changed cortical gyrification

Type 2 diabetes is a chronic disease that creates atrophic signatures in the brain, including decreases of total and regional volume of grey matter, white matter and cortical thickness. However, there is a lack of studies assessing cortical gyrification in type 2 diabetes. Changes in this emerging feature have been associated mainly with genetic legacy, but environmental factors may also play a role. Here, we investigated alterations of the gyrification index and classical morphometric measures in type 2 diabetes, a late acquired disease with complex aetiology with both underlying genetic and more preponderant environmental factors. In this cross‐sectional study, we analysed brain anatomical magnetic resonance images of 86 participants with type 2 diabetes and 40 healthy control participants, to investigate structural alterations in type 2 diabetes, including whole‐brain volumetric measures, local alterations of grey matter volume, cortical thickness and the gyrification index. We found concordant significant decrements in total and regional grey matter volume, and cortical thickness. Surprisingly, the cortical gyrification index was found to be mainly increased and mainly located in cortical sensory areas in type 2 diabetes. Moreover, alterations in gyrification correlated with clinical data, suggesting an influence of metabolic profile in alterations of gyrification in type 2 diabetes. Further studies should address causal influences of genetic and/or environmental factors in patterns of cortical gyrification in type 2 diabetes.

Best Practices in Structural Neuroimaging of Neurodevelopmental Disorders

Structural magnetic resonance imaging (sMRI) offers immense potential for increasing our understanding of how anatomical brain development relates to clinical symptoms and functioning in neurodevelopmental disorders. Clinical developmental sMRI may help identify neurobiological risk factors or markers that may ultimately assist in diagnosis and treatment. However, researchers and clinicians aiming to conduct sMRI studies of neurodevelopmental disorders face several methodological challenges. This review offers hands-on guidelines for clinical developmental sMRI. First, we present brain morphometry metrics and review evidence on typical developmental trajectories throughout adolescence, together with atypical trajectories in selected neurodevelopmental disorders. Next, we discuss challenges and good scientific practices in study design, image acquisition and analysis, and recent options to implement quality control. Finally, we discuss choices related to statistical analysis and interpretation of results. We call for greater completeness and transparency in the reporting of methods to advance understanding of structural brain alterations in neurodevelopmental disorders.

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.

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.

Shaping brain structure: Genetic and phylogenetic axes of macroscale organization of cortical thickness

The topology of the cerebral cortex has been proposed to provide an important source of constraint for the organization of cognition. In a sample of twins (n = 1113), we determined structural covariance of thickness to be organized along both a posterior-to-anterior and an inferior-to-superior axis. Both organizational axes were present when investigating the genetic correlation of cortical thickness, suggesting a strong genetic component in humans, and had a comparable organization in macaques, demonstrating they are phylogenetically conserved in primates. In both species, the inferior-superior dimension of cortical organization aligned with the predictions of dual-origin theory, and in humans, we found that the posterior-to-anterior axis related to a functional topography describing a continuum of functions from basic processes involved in perception and action to more abstract features of human cognition. Together, our study provides important insights into how functional and evolutionary patterns converge at the level of macroscale cortical structural organization.