Structural covariance networks are coupled to expression of genes enriched in supragranular layers of the human cortex.

Complex network topology is characteristic of many biological systems, including anatomical and functional brain networks (connectomes). Here, we first constructed a structural covariance network from MRI measures of cortical thickness on 296 healthy volunteers, aged 14-24 years. Next, we designed a new algorithm for matching sample locations from the Allen Brain Atlas to the nodes of the SCN. Subsequently we used this to define, transcriptomic brain networks by estimating gene co-expression between pairs of cortical regions. Finally, we explored the hypothesis that transcriptional networks and structural MRI connectomes are coupled. A transcriptional brain network (TBN) and a structural covariance network (SCN) were correlated across connection weights and showed qualitatively similar complex topological properties: assortativity, small-worldness, modularity, and a rich-club. In both networks, the weight of an edge was inversely related to the anatomical (Euclidean) distance between regions. There were differences between networks in degree and distance distributions: the transcriptional network had a less fat-tailed degree distribution and a less positively skewed distance distribution than the SCN. However, cortical areas connected to each other within modules of the SCN had significantly higher levels of whole genome co-expression than expected by chance. Nodes connected in the SCN had especially high levels of expression and co-expression of a human supragranular enriched (HSE) gene set that has been specifically located to supragranular layers of human cerebral cortex and is known to be important for large-scale, long-distance cortico-cortical connectivity. This coupling of brain transcriptome and connectome topologies was largely but not entirely accounted for by the common constraint of physical distance on both networks

Morphogenetic Principles of Brain Organisation in Health and Disease

Non-invasive neuroimaging methods, such as MRI, provide a window into the structure of the mammalian brain. However, despite the ubiquity of these methods, the biological interpretation of the information obtained using these tools remains elusive. In order to accurately link this macroscale data to microscale measurements, it is critical that the construct validity is high. This thesis provides novel analyses, pipelines and methods to: i) generate and validate maps of brain organisation obtained via MRI, and ii) demonstrate the utility of these methods in capturing elements of cognition and psychopathology. First, in Chapter 1, I review some of the neuroscientific context for the new methods presented, from cytoarchitecture to gene expression to connectomes. Chapters 2-4 introduce a new method, “Morphometric Similarity Mapping”, which captures the brain organisation of an individual by mapping the relationships of multiple features of the cerebral cortex. Chapter 2 focuses on the development of the analysis pipeline and the graph theoretical features of the resulting morphometric similarity networks (MSNs), with an emphasis on reproducibility. Chapter 3 highlights the generalisability of MSNs to the macaque monkey, linking MSNs to ex vivo tract tracing experiments and presenting new tools for processing non-human imaging data; as well as evidence that MSN topography is organised by cytoarchitectonic features. Chapter 4 is focused on determining the transcriptomic correlates of MSNs using publicly available gene expression maps, and on applying MSNs to examine the relationship between brain organisation and intelligence. Chapter 5 is dedicated to rigorous evaluation of the applicability of MSNs to measure specific disease-relevant phenotypes in 8 rare genetic disorder cohorts. This includes the validation of novel methods for utilising data from both single-cell sequencing technologies and differential gene expression experiments (in multiple tissue types) in analysing neuroimaging and bulk transcriptomic brain maps. Chapter 6 provides a brief summary and presents some ongoing and future projects expanding on this original work. It also importantly discusses a general framework of comparing brain maps, including MSNs and gene expression, as well as other canonical maps of brain structure and function. Altogether, this thesis presents and evaluates novel methods and applications for integrating multimodal neuroimaging data with genetic data derived from multiple tissue types and through various acquisition strategies. It also includes tools for performing these analyses in non-human primates, and pipelines for statistically comparing brain maps. These results not only provide insight into the manifestation of brain-related changes due to various components of human variation, but also provides a framework for evaluating this variation at multiple biological scales purely from non-invasive neuroimaging data

Gene transcription profiles associated with inter-modular hubs and connection distance in human functional magnetic resonance imaging networks.

Human functional magnetic resonance imaging (fMRI) brain networks have a complex topology comprising integrative components, e.g. long-distance inter-modular edges, that are theoretically associated with higher biological cost. Here, we estimated intra-modular degree, inter-modular degree and connection distance for each of 285 cortical nodes in multi-echo fMRI data from 38 healthy adults. We used the multivariate technique of partial least squares (PLS) to reduce the dimensionality of the relationships between these three nodal network parameters and prior microarray data on regional expression of 20 737 genes. The first PLS component defined a transcriptional profile associated with high intra-modular degree and short connection distance, whereas the second PLS component was associated with high inter-modular degree and long connection distance. Nodes in superior and lateral cortex with high inter-modular degree and long connection distance had local transcriptional profiles enriched for oxidative metabolism and mitochondria, and for genes specific to supragranular layers of human cortex. In contrast, primary and secondary sensory cortical nodes in posterior cortex with high intra-modular degree and short connection distance had transcriptional profiles enriched for RNA translation and nuclear components. We conclude that, as predicted, topologically integrative hubs, mediating long-distance connections between modules, are more costly in terms of mitochondrial glucose metabolism.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.PEV is supported by an MRC Bioinformatics Research Fellowship (MR/K020706/1). Functional MRI data acquisition was supported by a strategic award from the Wellcome Trust to the University of Cambridge (IMG, PBJ, ETB) and University College London (RJD, PF): the Neuroscience in Psychiatry Network (NSPN). Additional support was provided by the NIHR Cambridge Biomedical Research Centre. Access to gene expression data was provided by the Allen Institute for Brain Sciences Website: © 2015 Allen Institute for Brain Science. Allen Human Brain Atlas [Internet]. Available from: http://human.brain-map.org. FV is supported by a Gates Cambridge PhD studentship. KW is supported by the University of Cambridge MB/PhD Programme and the Wellcome Trust. We thank Gita Prabu, Roger Tait, Cinly Ooi, John Suckling and Becky Inkster for fMRI data collection and storage. ETB is employed half-time by the University of Cambridge and half-time by GlaxoSmithKline(GSK).This is the final version of the article. It first appeared from Royal Society Publishing via http://dx.doi.org/10.1098/rstb.2015.036

Synaptic and transcriptionally downregulated genes are associated with cortical thickness differences in autism.

Differences in cortical morphology-in particular, cortical volume, thickness and surface area-have been reported in individuals with autism. However, it is unclear what aspects of genetic and transcriptomic variation are associated with these differences. Here we investigate the genetic correlates of global cortical thickness differences (ΔCT) in children with autism. We used Partial Least Squares Regression (PLSR) on structural MRI data from 548 children (166 with autism, 295 neurotypical children and 87 children with ADHD) and cortical gene expression data from the Allen Institute for Brain Science to identify genetic correlates of ΔCT in autism. We identify that these genes are enriched for synaptic transmission pathways and explain significant variation in ΔCT. These genes are also significantly enriched for genes dysregulated in the autism post-mortem cortex (Odd Ratio (OR) = 1.11, Pcorrected  10-14), driven entirely by downregulated genes (OR = 1.87, Pcorrected  10-15). We validated the enrichment for downregulated genes in two independent data sets: Validation 1 (OR = 1.44, Pcorrected = 0.004) and Validation 2 (OR = 1.30; Pcorrected = 0.001). We conclude that transcriptionally downregulated genes implicated in autism are robustly associated with global changes in cortical thickness variability in children with autism

Morphometric Similarity Networks Detect Microscale Cortical Organization and Predict Inter-Individual Cognitive Variation.

Macroscopic cortical networks are important for cognitive function, but it remains challenging to construct anatomically plausible individual structural connectomes from human neuroimaging. We introduce a new technique for cortical network mapping based on inter-regional similarity of multiple morphometric parameters measured using multimodal MRI. In three cohorts (two human, one macaque), we find that the resulting morphometric similarity networks (MSNs) have a complex topological organization comprising modules and high-degree hubs. Human MSN modules recapitulate known cortical cytoarchitectonic divisions, and greater inter-regional morphometric similarity was associated with stronger inter-regional co-expression of genes enriched for neuronal terms. Comparing macaque MSNs with tract-tracing data confirmed that morphometric similarity was related to axonal connectivity. Finally, variation in the degree of human MSN nodes accounted for about 40% of between-subject variability in IQ. Morphometric similarity mapping provides a novel, robust, and biologically plausible approach to understanding how human cortical networks underpin individual differences in psychological functions

Gene transcription profiles associated with inter-modular hubs and connection distance in human fMRI networks

Graph theoretical methods have been widely used to investigate the topology of large-scale human brain networks constructed from resting state functional magnetic resonance imaging (fMRI). It has been demonstrated that such human functional connectomes have a complex topology comprising integrative components, such as hubs and inter-modular edges, that are associated with proxy markers of greater biological cost. In the absence of secure knowledge of the neurovascular mechanisms linking ensemble oscillations of neuronal populations to low frequency coupling or functional connectivity between regional fMRI time series, it has been challenging to validate fMRI network properties reductionistically. Supportive evidence to date has been mostly provided by analogous results on the relationships between integrative topology and biological cost in other nervous systems. Here, we use microarray data on brain regional expression of 20,737 genes to explore the relationships between fMRI network topology and transcription of genes annotated for biological processes and cellular components. We show that intra-modular degree and inter-modular degree are differently patterned in anatomical space, are differently associated with cytoarchitectonic classes of cortex, and are associated with distinct and statistically independent gene expression profiles. Genes strongly associated with nodes mediating many long-distance and inter-modular connections are significantly enriched for oxidative metabolism and mitochondria as well as for a subset of genes specifically enriched in human supragranular cortical layers. These results are directly supportive of the concept of high cost / high value network hubs in fMRI networks and point to the nascent opportunity to resolve the molecular and cellular substrates of human brain graphs

Cortical patterning of abnormal morphometric similarity in psychosis is associated with brain expression of schizophrenia-related genes.

Schizophrenia has been conceived as a disorder of brain connectivity, but it is unclear how this network phenotype is related to the underlying genetics. We used morphometric similarity analysis of MRI data as a marker of interareal cortical connectivity in three prior case-control studies of psychosis: in total, n = 185 cases and n = 227 controls. Psychosis was associated with globally reduced morphometric similarity in all three studies. There was also a replicable pattern of case-control differences in regional morphometric similarity, which was significantly reduced in patients in frontal and temporal cortical areas but increased in parietal cortex. Using prior brain-wide gene expression data, we found that the cortical map of case-control differences in morphometric similarity was spatially correlated with cortical expression of a weighted combination of genes enriched for neurobiologically relevant ontology terms and pathways. In addition, genes that were normally overexpressed in cortical areas with reduced morphometric similarity were significantly up-regulated in three prior post mortem studies of schizophrenia. We propose that this combined analysis of neuroimaging and transcriptional data provides insight into how previously implicated genes and proteins as well as a number of unreported genes in their topological vicinity on the protein interaction network may drive structural brain network changes mediating the genetic risk of schizophrenia.This study was supported by grants from the European Commission (PSYSCAN - Translating neuroimaging findings from research into clinical practice; ID: 603196) and the NIHR Cambridge Biomedical Research Centre (Mental Health). SEM holds a Henslow Fellowship at Lucy Cavendish College, University of Cambridge, funded by the Cambridge Philosophical Society. PEV was supported by the Medical Research Council (MR/K020706/1) and an MQ fellowship (MQF17_24) and is a Fellow of the Alan Turing Institute funded under the EPSRC grant EP/N510129/1. KJW was funded by an Alan Turing Institute Research Fellowship under EPSRC Research grant TU/A/000017. ETB is supported by a NIHR Senior Investigator Award

Shifts in myeloarchitecture characterise adolescent development of cortical gradients.

Funder: MQ: Transforming Mental Health; Grant(s): MQF17/24We studied an accelerated longitudinal cohort of adolescents and young adults (n = 234, two time points) to investigate dynamic reconfigurations in myeloarchitecture. Intracortical profiles were generated using magnetization transfer (MT) data, a myelin-sensitive magnetic resonance imaging contrast. Mixed-effect models of depth specific intracortical profiles demonstrated two separate processes i) overall increases in MT, and ii) flattening of the MT profile related to enhanced signal in mid-to-deeper layers, especially in heteromodal and unimodal association cortices. This development was independent of morphological changes. Enhanced MT in mid-to-deeper layers was found to spatially co-localise specifically with gene expression markers of oligodendrocytes. Interregional covariance analysis revealed that these intracortical changes contributed to a gradual differentiation of higher-order from lower-order systems. Depth-dependent trajectories of intracortical myeloarchitectural development contribute to the maturation of structural hierarchies in the human neocortex, providing a model for adolescent development that bridges microstructural and macroscopic scales of brain organisation

Schizotypy-Related Magnetization of Cortex in Healthy Adolescence Is Colocated With Expression of Schizophrenia-Related Genes.

BACKGROUND: Genetic risk is thought to drive clinical variation on a spectrum of schizophrenia-like traits, but the underlying changes in brain structure that mechanistically link genomic variation to schizotypal experience and behavior are unclear. METHODS: We assessed schizotypy using a self-reported questionnaire and measured magnetization transfer as a putative microstructural magnetic resonance imaging marker of intracortical myelination in 68 brain regions in 248 healthy young people (14-25 years of age). We used normative adult brain gene expression data and partial least squares analysis to find the weighted gene expression pattern that was most colocated with the cortical map of schizotypy-related magnetization. RESULTS: Magnetization was significantly correlated with schizotypy in the bilateral posterior cingulate cortex and precuneus (and for disorganized schizotypy, also in medial prefrontal cortex; all false discovery rate-corrected ps < .05), which are regions of the default mode network specialized for social and memory functions. The genes most positively weighted on the whole-genome expression map colocated with schizotypy-related magnetization were enriched for genes that were significantly downregulated in two prior case-control histological studies of brain gene expression in schizophrenia. Conversely, the most negatively weighted genes were enriched for genes that were transcriptionally upregulated in schizophrenia. Positively weighted (downregulated) genes were enriched for neuronal, specifically interneuronal, affiliations and coded a network of proteins comprising a few highly interactive "hubs" such as parvalbumin and calmodulin. CONCLUSIONS: Microstructural magnetic resonance imaging maps of intracortical magnetization can be linked to both the behavioral traits of schizotypy and prior histological data on dysregulated gene expression in schizophrenia