Cortical gene expression architecture links healthy neurodevelopment to the imaging, transcriptomics and genetics of autism and schizophrenia.

Human brain organization involves the coordinated expression of thousands of genes. For example, the first principal component (C1) of cortical transcription identifies a hierarchy from sensorimotor to association regions. In this study, optimized processing of the Allen Human Brain Atlas revealed two new components of cortical gene expression architecture, C2 and C3, which are distinctively enriched for neuronal, metabolic and immune processes, specific cell types and cytoarchitectonics, and genetic variants associated with intelligence. Using additional datasets (PsychENCODE, Allen Cell Atlas and BrainSpan), we found that C1-C3 represent generalizable transcriptional programs that are coordinated within cells and differentially phased during fetal and postnatal development. Autism spectrum disorder and schizophrenia were specifically associated with C1/C2 and C3, respectively, across neuroimaging, differential expression and genome-wide association studies. Evidence converged especially in support of C3 as a normative transcriptional program for adolescent brain development, which can lead to atypical supragranular cortical connectivity in people at high genetic risk for schizophrenia.

Childhood maltreatment influences adult brain structure through its effects on immune, metabolic, and psychosocial factors.

Childhood maltreatment (CM) leads to a lifelong susceptibility to mental ill-health which might be reflected by its effects on adult brain structure, perhaps indirectly mediated by its effects on adult metabolic, immune, and psychosocial systems. Indexing these systemic factors via body mass index (BMI), C-reactive protein (CRP), and rates of adult trauma (AT), respectively, we tested three hypotheses: (H1) CM has direct or indirect effects on adult trauma, BMI, and CRP; (H2) adult trauma, BMI, and CRP are all independently related to adult brain structure; and (H3) childhood maltreatment has indirect effects on adult brain structure mediated in parallel by BMI, CRP, and AT. Using path analysis and data from N = 116,887 participants in UK Biobank, we find that CM is related to greater BMI and AT levels, and that these two variables mediate CM's effects on CRP [H1]. Regression analyses on the UKB MRI subsample (N = 21,738) revealed that greater CRP and BMI were both independently related to a spatially convergent pattern of cortical effects (Spearman's ρ = 0.87) characterized by fronto-occipital increases and temporo-parietal reductions in thickness. Subcortically, BMI was associated with greater volume, AT with lower volume and CPR with effects in both directions [H2]. Finally, path models indicated that CM has indirect effects in a subset of brain regions mediated through its direct effects on BMI and AT and indirect effects on CRP [H3]. Results provide evidence that childhood maltreatment can influence brain structure decades after exposure by increasing individual risk toward adult trauma, obesity, and inflammation.

Transcriptional cartography integrates multiscale biology of the human cortex.

The cerebral cortex underlies many of our unique strengths and vulnerabilities, but efforts to understand human cortical organization are challenged by reliance on incompatible measurement methods at different spatial scales. Macroscale features such as cortical folding and functional activation are accessed through spatially dense neuroimaging maps, whereas microscale cellular and molecular features are typically measured with sparse postmortem sampling. Here, we integrate these distinct windows on brain organization by building upon existing postmortem data to impute, validate, and analyze a library of spatially dense neuroimaging-like maps of human cortical gene expression. These maps allow spatially unbiased discovery of cortical zones with extreme transcriptional profiles or unusually rapid transcriptional change which index distinct microstructure and predict neuroimaging measures of cortical folding and functional activation. Modules of spatially coexpressed genes define a family of canonical expression maps that integrate diverse spatial scales and temporal epochs of human brain organization - ranging from protein-protein interactions to large-scale systems for cognitive processing. These module maps also parse neuropsychiatric risk genes into subsets which tag distinct cyto-laminar features and differentially predict the location of altered cortical anatomy and gene expression in patients. Taken together, the methods, resources, and findings described here advance our understanding of human cortical organization and offer flexible bridges to connect scientific fields operating at different spatial scales of human brain research.

Connectome-based reservoir computing with the conn2res toolbox.

The connection patterns of neural circuits form a complex network. How signaling in these circuits manifests as complex cognition and adaptive behaviour remains the central question in neuroscience. Concomitant advances in connectomics and artificial intelligence open fundamentally new opportunities to understand how connection patterns shape computational capacity in biological brain networks. Reservoir computing is a versatile paradigm that uses high-dimensional, nonlinear dynamical systems to perform computations and approximate cognitive functions. Here we present conn2res: an open-source Python toolbox for implementing biological neural networks as artificial neural networks. conn2res is modular, allowing arbitrary network architecture and dynamics to be imposed. The toolbox allows researchers to input connectomes reconstructed using multiple techniques, from tract tracing to noninvasive diffusion imaging, and to impose multiple dynamical systems, from spiking neurons to memristive dynamics. The versatility of the conn2res toolbox allows us to ask new questions at the confluence of neuroscience and artificial intelligence. By reconceptualizing function as computation, conn2res sets the stage for a more mechanistic understanding of structure-function relationships in brain networks.

The neuropeptidergic connectome of C. elegans.

Efforts are ongoing to map synaptic wiring diagrams, or connectomes, to understand the neural basis of brain function. However, chemical synapses represent only one type of functionally important neuronal connection; in particular, extrasynaptic, "wireless" signaling by neuropeptides is widespread and plays essential roles in all nervous systems. By integrating single-cell anatomical and gene-expression datasets with biochemical analysis of receptor-ligand interactions, we have generated a draft connectome of neuropeptide signaling in the C. elegans nervous system. This network is characterized by high connection density, extended signaling cascades, autocrine foci, and a decentralized topology, with a large, highly interconnected core containing three constituent communities sharing similar patterns of input connectivity. Intriguingly, several key network hubs are little-studied neurons that appear specialized for peptidergic neuromodulation. We anticipate that the C. elegans neuropeptidergic connectome will serve as a prototype to understand how networks of neuromodulatory signaling are organized.

System-wide mapping of peptide-GPCR interactions in C. elegans.

Neuropeptides and peptide hormones are ancient, widespread signaling molecules that underpin almost all brain functions. They constitute a broad ligand-receptor network, mainly by binding to G protein-coupled receptors (GPCRs). However, the organization of the peptidergic network and roles of many peptides remain elusive, as our insight into peptide-receptor interactions is limited and many peptide GPCRs are still orphan receptors. Here we report a genome-wide peptide-GPCR interaction map in Caenorhabditis elegans. By reverse pharmacology screening of over 55,384 possible interactions, we identify 461 cognate peptide-GPCR couples that uncover a broad signaling network with specific and complex combinatorial interactions encoded across and within single peptidergic genes. These interactions provide insights into peptide functions and evolution. Combining our dataset with phylogenetic analysis supports peptide-receptor co-evolution and conservation of at least 14 bilaterian peptidergic systems in C. elegans. This resource lays a foundation for system-wide analysis of the peptidergic network.

Early adversity changes the economic conditions of mouse structural brain network organization.

Early adversity can change educational, cognitive, and mental health outcomes. However, the neural processes through which early adversity exerts these effects remain largely unknown. We used generative network modeling of the mouse connectome to test whether unpredictable postnatal stress shifts the constraints that govern the organization of the structural connectome. A model that trades off the wiring cost of long-distance connections with topological homophily (i.e., links between regions with shared neighbors) generated simulations that successfully replicate the rodent connectome. The imposition of early life adversity shifted the best-performing parameter combinations toward zero, heightening the stochastic nature of the generative process. Put simply, unpredictable postnatal stress changes the economic constraints that reproduce rodent connectome organization, introducing greater randomness into the development of the simulations. While this change may constrain the development of cognitive abilities, it could also reflect an adaptive mechanism that facilitates effective responses to future challenges.

Semantic Speech Networks Linked to Formal Thought Disorder in Early Psychosis.

BACKGROUND AND HYPOTHESIS: Mapping a patient's speech as a network has proved to be a useful way of understanding formal thought disorder in psychosis. However, to date, graph theory tools have not explicitly modelled the semantic content of speech, which is altered in psychosis. STUDY DESIGN: We developed an algorithm, "netts," to map the semantic content of speech as a network, then applied netts to construct semantic speech networks for a general population sample (N = 436), and a clinical sample comprising patients with first episode psychosis (FEP), people at clinical high risk of psychosis (CHR-P), and healthy controls (total N = 53). STUDY RESULTS: Semantic speech networks from the general population were more connected than size-matched randomized networks, with fewer and larger connected components, reflecting the nonrandom nature of speech. Networks from FEP patients were smaller than from healthy participants, for a picture description task but not a story recall task. For the former task, FEP networks were also more fragmented than those from controls; showing more connected components, which tended to include fewer nodes on average. CHR-P networks showed fragmentation values in-between FEP patients and controls. A clustering analysis suggested that semantic speech networks captured novel signals not already described by existing NLP measures. Network features were also related to negative symptom scores and scores on the Thought and Language Index, although these relationships did not survive correcting for multiple comparisons. CONCLUSIONS: Overall, these data suggest that semantic networks can enable deeper phenotyping of formal thought disorder in psychosis. Whilst here we focus on network fragmentation, the semantic speech networks created by Netts also contain other, rich information which could be extracted to shed further light on formal thought disorder. We are releasing Netts as an open Python package alongside this manuscript.

Molecular and network-level mechanisms explaining individual differences in autism spectrum disorder.

The mechanisms underlying phenotypic heterogeneity in autism spectrum disorder (ASD) are not well understood. Using a large neuroimaging dataset, we identified three latent dimensions of functional brain network connectivity that predicted individual differences in ASD behaviors and were stable in cross-validation. Clustering along these three dimensions revealed four reproducible ASD subgroups with distinct functional connectivity alterations in ASD-related networks and clinical symptom profiles that were reproducible in an independent sample. By integrating neuroimaging data with normative gene expression data from two independent transcriptomic atlases, we found that within each subgroup, ASD-related functional connectivity was explained by regional differences in the expression of distinct ASD-related gene sets. These gene sets were differentially associated with distinct molecular signaling pathways involving immune and synapse function, G-protein-coupled receptor signaling, protein synthesis and other processes. Collectively, our findings delineate atypical connectivity patterns underlying different forms of ASD that implicate distinct molecular signaling mechanisms.

Computational Models of Typical and Atypical Brain Network Development.

Over the last decade, the organization of brain networks at both micro- and macroscales has become a key focus of neuroscientific inquiry. This has revealed fundamental features of brain network organization-small-worldness, modularity, heavy-tailed degree distributions-and has highlighted how these structural features support brain function. However, the driving forces that shape brain networks over the course of development have begun to be explored only recently. Here, we review recent efforts to gain insights into the mechanisms of brain development through generative modeling of both macroscale human brain networks and microscale cellular connectomes in Caenorhabditis elegans and other organisms. We show how these mathematical models can begin to shed light on the biological processes that drive and constrain the development of brain networks. Finally, we show how generative network models can translate genetic and environmental differences into variability in developmental trajectories, leading to diverse cognitive and mental health outcomes in children and young people.

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.

Regional gene expression signatures are associated with sex-specific functional connectivity changes in depression.

The neural substrates of depression may differ in men and women, but the underlying mechanisms are incompletely understood. Here, we show that depression is associated with sex-specific patterns of abnormal functional connectivity in the default mode network and in five regions of interest with sexually dimorphic transcriptional effects. Regional differences in gene expression in two independent datasets explained the neuroanatomical distribution of abnormal connectivity. These gene sets varied by sex and were strongly enriched for genes implicated in depression, synapse function, immune signaling, and neurodevelopment. In an independent sample, we confirmed the prediction that individual differences in default mode network connectivity are explained by inferred brain expression levels for six depression-related genes, including PCDH8, a brain-specific protocadherin integral membrane protein implicated in activity-related synaptic reorganization. Together, our results delineate both shared and sex-specific changes in the organization of depression-related functional networks, with implications for biomarker development and fMRI-guided therapeutic neuromodulation.

Sexually divergent development of depression-related brain networks during healthy human adolescence.

Sexual differences in human brain development could be relevant to sex differences in the incidence of depression during adolescence. We tested for sex differences in parameters of normative brain network development using fMRI data on N = 298 healthy adolescents, aged 14 to 26 years, each scanned one to three times. Sexually divergent development of functional connectivity was located in the default mode network, limbic cortex, and subcortical nuclei. Females had a more "disruptive" pattern of development, where weak functional connectivity at age 14 became stronger during adolescence. This fMRI-derived map of sexually divergent brain network development was robustly colocated with i prior loci of reward-related brain activation ii a map of functional dysconnectivity in major depressive disorder (MDD), and iii an adult brain gene transcriptional pattern enriched for genes on the X chromosome, neurodevelopmental genes, and risk genes for MDD. We found normative sexual divergence in adolescent development of a cortico-subcortical brain functional network that is relevant to depression.

Atlas of lesion locations and postsurgical seizure freedom in focal cortical dysplasia: A MELD study.

OBJECTIVE: Drug-resistant focal epilepsy is often caused by focal cortical dysplasias (FCDs). The distribution of these lesions across the cerebral cortex and the impact of lesion location on clinical presentation and surgical outcome are largely unknown. We created a neuroimaging cohort of patients with individually mapped FCDs to determine factors associated with lesion location and predictors of postsurgical outcome. METHODS: The MELD (Multi-centre Epilepsy Lesion Detection) project collated a retrospective cohort of 580 patients with epilepsy attributed to FCD from 20 epilepsy centers worldwide. Magnetic resonance imaging-based maps of individual FCDs with accompanying demographic, clinical, and surgical information were collected. We mapped the distribution of FCDs, examined for associations between clinical factors and lesion location, and developed a predictive model of postsurgical seizure freedom. RESULTS: FCDs were nonuniformly distributed, concentrating in the superior frontal sulcus, frontal pole, and temporal pole. Epilepsy onset was typically before the age of 10 years. Earlier epilepsy onset was associated with lesions in primary sensory areas, whereas later epilepsy onset was associated with lesions in association cortices. Lesions in temporal and occipital lobes tended to be larger than frontal lobe lesions. Seizure freedom rates varied with FCD location, from around 30% in visual, motor, and premotor areas to 75% in superior temporal and frontal gyri. The predictive model of postsurgical seizure freedom had a positive predictive value of 70% and negative predictive value of 61%. SIGNIFICANCE: FCD location is an important determinant of its size, the age at epilepsy onset, and the likelihood of seizure freedom postsurgery. Our atlas of lesion locations can be used to guide the radiological search for subtle lesions in individual patients. Our atlas of regional seizure freedom rates and associated predictive model can be used to estimate individual likelihoods of postsurgical seizure freedom. Data-driven atlases and predictive models are essential for evidence-based, precision medicine and risk counseling in epilepsy.

Natural Language Processing markers in first episode psychosis and people at clinical high-risk.

Recent work has suggested that disorganised speech might be a powerful predictor of later psychotic illness in clinical high risk subjects. To that end, several automated measures to quantify disorganisation of transcribed speech have been proposed. However, it remains unclear which measures are most strongly associated with psychosis, how different measures are related to each other and what the best strategies are to collect speech data from participants. Here, we assessed whether twelve automated Natural Language Processing markers could differentiate transcribed speech excerpts from subjects at clinical high risk for psychosis, first episode psychosis patients and healthy control subjects (total N = 54). In-line with previous work, several measures showed significant differences between groups, including semantic coherence, speech graph connectivity and a measure of whether speech was on-topic, the latter of which outperformed the related measure of tangentiality. Most NLP measures examined were only weakly related to each other, suggesting they provide complementary information. Finally, we compared the ability of transcribed speech generated using different tasks to differentiate the groups. Speech generated from picture descriptions of the Thematic Apperception Test and a story re-telling task outperformed free speech, suggesting that choice of speech generation method may be an important consideration. Overall, quantitative speech markers represent a promising direction for future clinical applications.

A generative network model of neurodevelopmental diversity in structural brain organization.

The formation of large-scale brain networks, and their continual refinement, represent crucial developmental processes that can drive individual differences in cognition and which are associated with multiple neurodevelopmental conditions. But how does this organization arise, and what mechanisms drive diversity in organization? We use generative network modeling to provide a computational framework for understanding neurodevelopmental diversity. Within this framework macroscopic brain organization, complete with spatial embedding of its organization, is an emergent property of a generative wiring equation that optimizes its connectivity by renegotiating its biological costs and topological values continuously over time. The rules that govern these iterative wiring properties are controlled by a set of tightly framed parameters, with subtle differences in these parameters steering network growth towards different neurodiverse outcomes. Regional expression of genes associated with the simulations converge on biological processes and cellular components predominantly involved in synaptic signaling, neuronal projection, catabolic intracellular processes and protein transport. Together, this provides a unifying computational framework for conceptualizing the mechanisms and diversity in neurodevelopment, capable of integrating different levels of analysis-from genes to cognition.

Author Correction: Transcriptomic and cellular decoding of regional brain vulnerability to neurogenetic disorders.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Compulsivity is linked to reduced adolescent development of goal-directed control and frontostriatal functional connectivity.

A characteristic of adaptive behavior is its goal-directed nature. An ability to act in a goal-directed manner is progressively refined during development, but this refinement can be impacted by the emergence of psychiatric disorders. Disorders of compulsivity have been framed computationally as a deficit in model-based control, and have been linked also to abnormal frontostriatal connectivity. However, the developmental trajectory of model-based control, including an interplay between its maturation and an emergence of compulsivity, has not been characterized. Availing of a large sample of healthy adolescents (n = 569) aged 14 to 24 y, we show behaviorally that over the course of adolescence there is a within-person increase in model-based control, and this is more pronounced in younger participants. Using a bivariate latent change score model, we provide evidence that the presence of higher compulsivity traits is associated with an atypical profile of this developmental maturation in model-based control. Resting-state fMRI data from a subset of the behaviorally assessed subjects (n = 230) revealed that compulsivity is associated with a less pronounced change of within-subject developmental remodeling of functional connectivity, specifically between the striatum and a frontoparietal network. Thus, in an otherwise clinically healthy population sample, in early development, individual differences in compulsivity are linked to the developmental trajectory of model-based control and a remodeling of frontostriatal connectivity.

Transcriptomic and cellular decoding of regional brain vulnerability to neurogenetic disorders.

Neurodevelopmental disorders have a heritable component and are associated with region specific alterations in brain anatomy. However, it is unclear how genetic risks for neurodevelopmental disorders are translated into spatially patterned brain vulnerabilities. Here, we integrated cortical neuroimaging data from patients with neurodevelopmental disorders caused by genomic copy number variations (CNVs) and gene expression data from healthy subjects. For each of the six investigated disorders, we show that spatial patterns of cortical anatomy changes in youth are correlated with cortical spatial expression of CNV genes in neurotypical adults. By transforming normative bulk-tissue cortical expression data into cell-type expression maps, we link anatomical change maps in each analysed disorder to specific cell classes as well as the CNV-region genes they express. Our findings reveal organizing principles that regulate the mapping of genetic risks onto regional brain changes in neurogenetic disorders. Our findings will enable screening for candidate molecular mechanisms from readily available neuroimaging data.

Major Depressive Disorder Is Associated With Differential Expression of Innate Immune and Neutrophil-Related Gene Networks in Peripheral Blood: A Quantitative Review of Whole-Genome Transcriptional Data From Case-Control Studies.

BACKGROUND: Whole-genome transcription has been measured in peripheral blood samples as a candidate biomarker of inflammation associated with major depressive disorder. METHODS: We searched for all case-control studies on major depressive disorder that reported microarray or RNA sequencing measurements on whole blood or peripheral blood mononuclear cells. Primary datasets were reanalyzed, when openly accessible, to estimate case-control differences and to evaluate the functional roles of differentially expressed gene lists by technically harmonized methods. RESULTS: We found 10 eligible studies (N = 1754 depressed cases and N = 1145 healthy controls). Fifty-two genes were called significant by 2 of the primary studies (published overlap list). After harmonization of analysis across 8 accessible datasets (n = 1706 cases, n = 1098 controls), 272 genes were coincidentally listed in the top 3% most differentially expressed genes in 2 or more studies of whole blood or peripheral blood mononuclear cells with concordant direction of effect (harmonized overlap list). By meta-analysis of standardized mean difference across 4 studies of whole-blood samples (n = 1567 cases, n = 954 controls), 343 genes were found with false discovery rate <5% (standardized mean difference meta-analysis list). These 3 lists intersected significantly. Genes abnormally expressed in major depressive disorder were enriched for innate immune-related functions, coded for nonrandom protein-protein interaction networks, and coexpressed in the normative transcriptome module specialized for innate immune and neutrophil functions. CONCLUSIONS: Quantitative review of existing case-control data provided robust evidence for abnormal expression of gene networks important for the regulation and implementation of innate immune response. Further development of white blood cell transcriptional biomarkers for inflamed depression seems warranted.