Cell-Specific Gene-Expression Profiles and Cortical Thickness in the Human Brain

Brain atrophy progression in Parkinson's disease is shaped by connectivity and local vulnerability

Brain atrophy has been reported in the early stages of Parkinson's disease, but there have been few longitudinal studies. How intrinsic properties of the brain, such as anatomical connectivity, local cell-type distribution and gene expression combine to determine the pattern of disease progression also remains unknown. One hypothesis proposes that the disease stems from prion-like propagation of misfolded alpha-synuclein via the connectome that might cause varying degrees of tissue damage based on local properties. Here, we used MRI data from the Parkinson Progression Markers Initiative to map the progression of brain atrophy over 1, 2 and 4 years compared with baseline. We derived atrophy maps for four time points using deformation-based morphometry applied to T₁-weighted MRI from 120 de novo Parkinson's disease patients, 74 of whom had imaging at all four time points (50 Men: 24 Women) and 157 healthy control participants (115 Men: 42 Women). In order to determine factors that may influence neurodegeneration, we related atrophy progression to brain structural and functional connectivity, cell-type expression and gene ontology enrichment analyses. After regressing out the expected age and sex effects associated with normal ageing, we found that atrophy significantly progressed over 2 and 4 years in the caudate, nucleus accumbens, hippocampus and posterior cortical regions. This progression was shaped by both structural and functional brain connectivity. Also, the progression of atrophy was more pronounced in regions with a higher expression of genes related to synapses and was inversely related to the prevalence of oligodendrocytes and endothelial cells. In sum, we demonstrate that the progression of atrophy in Parkinson's disease is in line with the prion-like propagation hypothesis of alpha-synuclein and provide evidence that synapses may be especially vulnerable to synucleinopathy. In addition to identifying vulnerable brain regions, this study reveals different factors that may be implicated in the neurotoxic mechanisms leading to progression in Parkinson's disease. All brain maps generated here are available on request.

Standardizing workflows in imaging transcriptomics with the abagen toolbox

Gene expression fundamentally shapes the structural and functional architecture of the human brain. Open-access transcriptomic datasets like the Allen Human Brain Atlas provide an unprecedented ability to examine these mechanisms in vivo; however, a lack of standardization across research groups has given rise to myriad processing pipelines for using these data. Here, we develop the abagen toolbox, an open-access software package for working with transcriptomic data, and use it to examine how methodological variability influences the outcomes of research using the Allen Human Brain Atlas. Applying three prototypical analyses to the outputs of 750,000 unique processing pipelines, we find that choice of pipeline has a large impact on research findings, with parameters commonly varied in the literature influencing correlations between derived gene expression and other imaging phenotypes by as much as ρ ≥ 1.0. Our results further reveal an ordering of parameter importance, with processing steps that influence gene normalization yielding the greatest impact on downstream statistical inferences and conclusions. The presented work and the development of the abagen toolbox lay the foundation for more standardized and systematic research in imaging transcriptomics, and will help to advance future understanding of the influence of gene expression in the human brain.

Imaging Transcriptomics of Brain Disorders

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

Effects of severe obesity and sleeve gastrectomy on cortical thickness in adolescents

OBJECTIVE

Neurocognitive differences in pediatric obesity may be underpinned by cortical structural alterations. Differences in cortical thickness associated with severe obesity were examined, and preliminary evidence was sought for changes following vertical sleeve gastrectomy (VSG).

METHODS

A total of 18 adolescents with severe obesity (OB) and 17 without obesity (nOB), aged 14 to 21, underwent T1-weighted structural magnetic resonance imaging. A subset was scanned twice 5 months apart to compare cortical thickness following VSG (n = 6) with two control groups: wait-listed (n = 9) and nOB (n = 12).

RESULTS

At baseline, OB had a thinner cortex than nOB in motor and superior parietal cortices. At follow-up, VSG adolescents lost weight, the wait-listed group gained weight, and nOB did not change. Group × Time interactions indicated that VSG had cortical thinning in orbitofrontal, primary sensorimotor, superior, and middle temporal cortices and thickening in lingual, fusiform, and lateral occipital cortices. Wait-listed and nOB groups largely did not change.

CONCLUSIONS

Severe obesity is associated with a thinner cortex in motor and attentional function-associated regions. VSG resulted in cortical thinning in reward valuation, sensory, and perceptual regions and thickening in visual regions. Surgery-related changes in regions distinct from those associated with obesity suggest compensation, rather than normalization. These results provide preliminary evidence supporting structural neural alterations following sleeve gastrectomy.

CSF sTREM2 and Tau Work Together in Predicting Increased Temporal Lobe Atrophy in Older Adults

Neuroinflammation may be a key factor in brain atrophy in aging and age-related neurodegenerative disease. The objective of this study was to test the association between microglial expression of soluble Triggering Receptor Expressed on Myeloid Cells 2 (sTREM2), as a measure of neuroinflammation, and brain atrophy in cognitively unimpaired older adults. Brain magnetic resonance imagings (MRIs) and cerebrospinal fluid (CSF) sTREM2, total tau (t-tau), phosphorylated181 tau (p-tau), and Aβ42 were analyzed in 115 cognitively unimpaired older adults, classified according to the A/T/(N)-framework. MRIs were repeated after 2 (n = 95) and 4 (n = 62) years. High baseline sTREM2 was associated with accelerated cortical thinning in the temporal cortex of the left hemisphere, as well as bilateral hippocampal atrophy, independently of age, Aβ42, and tau. sTREM2-related atrophy only marginally increased with biomarker positivity across the AD continuum (A-T- #x2292; A+T- #x2292; A+T+) but was significantly stronger in participants with a high level of p-tau (T+). sTREM2-related cortical thinning correlated significantly with areas of high microglial-specific gene expression in the Allen Human Brain Atlas. In conclusion, increased CSF sTREM2 was associated with accelerated cortical and hippocampal atrophy in cognitively unimpaired older participants, particularly in individuals with tau pathology. This suggests a link between neuroinflammation, neurodegeneration, and amyloid-independent tauopathy.

Unveiling the neuroimaging-genetic intersections in the human brain

PURPOSE OF REVIEW

The prevalence of new public datasets of brain-wide and single-cell transcriptome data has created new opportunities to link neuroimaging findings with genetic data. The aim of this study is to present the different methodological approaches that have been used to combine this data.

RECENT FINDINGS

Drawing from various sources of open access data, several studies have been able to correlate neuroimaging maps with spatial distribution of brain expression. These efforts have enabled researchers to identify functional annotations of related genes, identify specific cell types related to brain phenotypes, study the expression of genes across life span and highlight the importance of selected brain genes in disease genetic networks.

SUMMARY

New transcriptome datasets and methodological approaches complement current neuroimaging work and will be crucial to improve our understanding of the biological mechanism that underlies many neurological conditions.

New insights into the dynamic development of the cerebral cortex in childhood and adolescence: Integrating macro- and microstructural MRI findings

Through dynamic transactional processes between genetic and environmental factors, childhood and adolescence involve reorganization and optimization of the cerebral cortex. The cortex and its development plays a crucial role for prototypical human cognitive abilities. At the same time, many common mental disorders appear during these critical phases of neurodevelopment. Magnetic resonance imaging (MRI) can indirectly capture several multifaceted changes of cortical macro- and microstructure, of high relevance to further our understanding of the neural foundation of cognition and mental health. Great progress has been made recently in mapping the typical development of cortical morphology. Moreover, newer less explored MRI signal intensity and specialized quantitative T2 measures have been applied to assess microstructural cortical development. We review recent findings of typical postnatal macro- and microstructural development of the cerebral cortex from early childhood to young adulthood. We cover studies of cortical volume, thickness, area, gyrification, T1-weighted (T1w) tissue contrasts such a grey/white matter contrast, T1w/T2w ratio, magnetization transfer and myelin water fraction. Finally, we integrate imaging studies with cortical gene expression findings to further our understanding of the underlying neurobiology of the developmental changes, bridging the gap between ex vivo histological- and in vivo MRI studies.

Pubertal Testosterone and the Structure of the Cerebral Cortex in Young Men

Adolescence is a period of brain maturation that may involve a second wave of organizational effects of sex steroids on the brain. Rodent studies suggest that, overall, organizational effects of gonadal steroid hormones decrease from the prenatal/perinatal period to adulthood. Here we used multimodal magnetic resonance imaging to investigate whether 1) testosterone exposure during adolescence (9-17 years) correlates with the structure of cerebral cortex in young men (n = 216, 19 years of age); 2) this relationship is modulated by the timing of testosterone surge during puberty. Our results showed that pubertal testosterone correlates with structural properties of the cerebral cortex, as captured by principal component analysis of T1 and T2 relaxation times, myelin water fraction, magnetization transfer ratio, fractional anisotropy and mean diffusivity. Many of the correlations between pubertal testosterone and the cortical structure were stronger in individuals with earlier (vs. later) testosterone surge. We also demonstrated that the strength of the relationship between pubertal testosterone and cortical structure across the cerebral cortex varies as a function of inter-regional profiles of gene expression specific to dendrites, axonal cytoskeleton, and myelin. This finding suggests that the cellular substrate underlying the relationships between pubertal testosterone and cerebral cortex involves both dendritic arbor and axon.

An expanding manifold in transmodal regions characterizes adolescent reconfiguration of structural connectome organization

Adolescence is a critical time for the continued maturation of brain networks. Here, we assessed structural connectome development in a large longitudinal sample ranging from childhood to young adulthood. By projecting high-dimensional connectomes into compact manifold spaces, we identified a marked expansion of structural connectomes, with strongest effects in transmodal regions during adolescence. Findings reflected increased within-module connectivity together with increased segregation, indicating increasing differentiation of higher-order association networks from the rest of the brain. Projection of subcortico-cortical connectivity patterns into these manifolds showed parallel alterations in pathways centered on the caudate and thalamus. Connectome findings were contextualized via spatial transcriptome association analysis, highlighting genes enriched in cortex, thalamus, and striatum. Statistical learning of cortical and subcortical manifold features at baseline and their maturational change predicted measures of intelligence at follow-up. Our findings demonstrate that connectome manifold learning can bridge the conceptual and empirical gaps between macroscale network reconfigurations, microscale processes, and cognitive outcomes in adolescent development.

Integration of human whole-brain transcriptome and neuroimaging data: Practical considerations of current available methods

The Allen Human Brain Atlas (AHBA) is the first example of human brain transcriptomic mappings and detailed anatomical annotations which, for the first time, has allowed the integration of human brain transcriptomics with neuroimaging. This has been made possible because the AHBA offered an original dataset that contains mRNA expression measures for >20,000 genes covering the whole brain and, critically, in a standard stereotaxic space. In recent years many different methods have been used to integrate this data set with brain imaging data, although this endeavour has lacked harmony in terms of the workflow of data processing and subsequent analyses. In this work we discuss five main issues that experience has highlighted as in need of thorough consideration when integrating the AHBA with neuroimaging. These concerns are corroborated by comparing the performance of three different publicly available methods in correlating the same measures of serotonin receptors density with the correspondent AHBA mRNA maps. In this representative case, we were able to show how these methods can lead to discrepant results, suggesting that processing options are not neutral. We believe that the field should take into serious consideration these issues as they could undermine reproducibility and, in the end, the intrinsic value of the AHBA. We also advise on possible strategies to overcome these discrepancies. Finally, we encourage authors towards practices that will improve reproducibility such as transparency in reporting, algorithm and data sharing, collaboration.

Virtual Histology of Cortical Thickness and Shared Neurobiology in 6 Psychiatric Disorders

IMPORTANCE

Large-scale neuroimaging studies have revealed group differences in cortical thickness across many psychiatric disorders. The underlying neurobiology behind these differences is not well understood.

OBJECTIVE

To determine neurobiologic correlates of group differences in cortical thickness between cases and controls in 6 disorders: attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder (BD), major depressive disorder (MDD), obsessive-compulsive disorder (OCD), and schizophrenia.

DESIGN, SETTING, AND PARTICIPANTS

Profiles of group differences in cortical thickness between cases and controls were generated using T1-weighted magnetic resonance images. Similarity between interregional profiles of cell-specific gene expression and those in the group differences in cortical thickness were investigated in each disorder. Next, principal component analysis was used to reveal a shared profile of group difference in thickness across the disorders. Analysis for gene coexpression, clustering, and enrichment for genes associated with these disorders were conducted. Data analysis was conducted between June and December 2019. The analysis included 145 cohorts across 6 psychiatric disorders drawn from the ENIGMA consortium. The numbers of cases and controls in each of the 6 disorders were as follows: ADHD: 1814 and 1602; ASD: 1748 and 1770; BD: 1547 and 3405; MDD: 2658 and 3572; OCD: 2266 and 2007; and schizophrenia: 2688 and 3244.

MAIN OUTCOMES AND MEASURES

Interregional profiles of group difference in cortical thickness between cases and controls.

RESULTS

A total of 12 721 cases and 15 600 controls, ranging from ages 2 to 89 years, were included in this study. Interregional profiles of group differences in cortical thickness for each of the 6 psychiatric disorders were associated with profiles of gene expression specific to pyramidal (CA1) cells, astrocytes (except for BD), and microglia (except for OCD); collectively, gene-expression profiles of the 3 cell types explain between 25% and 54% of variance in interregional profiles of group differences in cortical thickness. Principal component analysis revealed a shared profile of difference in cortical thickness across the 6 disorders (48% variance explained); interregional profile of this principal component 1 was associated with that of the pyramidal-cell gene expression (explaining 56% of interregional variation). Coexpression analyses of these genes revealed 2 clusters: (1) a prenatal cluster enriched with genes involved in neurodevelopmental (axon guidance) processes and (2) a postnatal cluster enriched with genes involved in synaptic activity and plasticity-related processes. These clusters were enriched with genes associated with all 6 psychiatric disorders.

CONCLUSIONS AND RELEVANCE

In this study, shared neurobiologic processes were associated with differences in cortical thickness across multiple psychiatric disorders. These processes implicate a common role of prenatal development and postnatal functioning of the cerebral cortex in these disorders.

Cellular correlates of cortical thinning throughout the lifespan

Cortical thinning occurs throughout the entire life and extends to late-life neurodegeneration, yet the neurobiological substrates are poorly understood. Here, we used a virtual-histology technique and gene expression data from the Allen Human Brain Atlas to compare the regional profiles of longitudinal cortical thinning through life (4004 magnetic resonance images [MRIs]) with those of gene expression for several neuronal and non-neuronal cell types. The results were replicated in three independent datasets. We found that inter-regional profiles of cortical thinning related to expression profiles for marker genes of CA1 pyramidal cells, astrocytes and, microglia during development and in aging. During the two stages of life, the relationships went in opposite directions: greater gene expression related to less thinning in development and vice versa in aging. The association between cortical thinning and cell-specific gene expression was also present in mild cognitive impairment and Alzheimer's Disease. These findings suggest a role of astrocytes and microglia in promoting and supporting neuronal growth and dendritic structures through life that affects cortical thickness during development, aging, and neurodegeneration. Overall, the findings contribute to our understanding of the neurobiology underlying variations in MRI-derived estimates of cortical thinning through life and late-life disease.

Adiposity-related insulin resistance and thickness of the cerebral cortex in middle-aged adults

The thickness of the cerebral cortex decreases with ageing. Recent research suggests that obesity and type 2 diabetes mellitus may accelerate this cortical thinning, and that obesity-related insulin resistance may be a shared mechanistic pathway. Ageing of the cerebral cortex demonstrates sex-specific trajectories, with a gradual shift towards accelerated thinning beginning in midlife. Here, we investigated whether adiposity-related insulin resistance is associated with lower thickness of the human cerebral cortex in a community-based sample of middle-aged adults. We studied 533 adult participants (36-65 years) from the Saguenay Youth Study. Adiposity was assessed with bioimpedance, and insulin resistance was evaluated from a fasting blood sample with the homeostatic model assessment of insulin resistance (HOMA-IR). Associations between adiposity-related insulin resistance (adiposity/IR) and cortical thickness were assessed with linear models, separately in males and females younger or older than 50 years. Potential biological underpinnings were investigated with virtual histology. Adiposity/IR was associated with lower cortical thickness in females older than 50 years but not in males or younger females. The strength of the association varied across the cerebral cortex, with regions of the lateral frontal and parietal cortices and the superior temporal cortex demonstrating most pronounced thinning. Based on virtual histology, adiposity/IR-related cortical thinning may involve neurones, astrocytes, oligodendrocytes and ependymal cells acting so that they lower the cortical potential for synaptogenesis, formation of dendritic spines, production of extracellular matrix and myelination. Adiposity-related insulin resistance is associated with lower cortical thickness in middle-aged women older than 50 years. This aspect of thinning may involve neuronal and glial cells in a way that lowers the capacity of the cerebral cortex for neuronal plasticity and maintenance of myelination.

Poor Self-Reported Sleep is Related to Regional Cortical Thinning in Aging but not Memory Decline-Results From the Lifebrain Consortium

We examined whether sleep quality and quantity are associated with cortical and memory changes in cognitively healthy participants across the adult lifespan. Associations between self-reported sleep parameters (Pittsburgh Sleep Quality Index, PSQI) and longitudinal cortical change were tested using five samples from the Lifebrain consortium (n = 2205, 4363 MRIs, 18–92 years). In additional analyses, we tested coherence with cell-specific gene expression maps from the Allen Human Brain Atlas, and relations to changes in memory performance. “PSQI # 1 Subjective sleep quality” and “PSQI #5 Sleep disturbances” were related to thinning of the right lateral temporal cortex, with lower quality and more disturbances being associated with faster thinning. The association with “PSQI #5 Sleep disturbances” emerged after 60 years, especially in regions with high expression of genes related to oligodendrocytes and S1 pyramidal neurons. None of the sleep scales were related to a longitudinal change in episodic memory function, suggesting that sleep-related cortical changes were independent of cognitive decline. The relationship to cortical brain change suggests that self-reported sleep parameters are relevant in lifespan studies, but small effect sizes indicate that self-reported sleep is not a good biomarker of general cortical degeneration in healthy older adults.

Cortical morphology at birth reflects spatiotemporal patterns of gene expression in the fetal human brain

Interruption to gestation through preterm birth can significantly impact cortical development and have long-lasting adverse effects on neurodevelopmental outcome. We compared cortical morphology captured by high-resolution, multimodal magnetic resonance imaging (MRI) in n = 292 healthy newborn infants (mean age at birth = 39.9 weeks) with regional patterns of gene expression in the fetal cortex across gestation (n = 156 samples from 16 brains, aged 12 to 37 postconceptional weeks [pcw]). We tested the hypothesis that noninvasive measures of cortical structure at birth mirror areal differences in cortical gene expression across gestation, and in a cohort of n = 64 preterm infants (mean age at birth = 32.0 weeks), we tested whether cortical alterations observed after preterm birth were associated with altered gene expression in specific developmental cell populations. Neonatal cortical structure was aligned to differential patterns of cell-specific gene expression in the fetal cortex. Principal component analysis (PCA) of 6 measures of cortical morphology and microstructure showed that cortical regions were ordered along a principal axis, with primary cortex clearly separated from heteromodal cortex. This axis was correlated with estimated tissue maturity, indexed by differential expression of genes expressed by progenitor cells and neurons, and engaged in stem cell differentiation, neuron migration, and forebrain development. Preterm birth was associated with altered regional MRI metrics and patterns of differential gene expression in glial cell populations. The spatial patterning of gene expression in the developing cortex was thus mirrored by regional variation in cortical morphology and microstructure at term, and this was disrupted by preterm birth. This work provides a framework to link molecular mechanisms to noninvasive measures of cortical development in early life and highlights novel pathways to injury in neonatal populations at increased risk of neurodevelopmental disorder.

Interruption to gestation through preterm birth can significantly impact cortical development and have long-lasting adverse effects on neurodevelopmental outcome. A large neuroimaging study of newborn infants reveals how their cortical structure at birth is associated with patterns of gene expression in the fetal cortex and how this relationship is affected by preterm birth.

Assessment of Neurobiological Mechanisms of Cortical Thinning During Childhood and Adolescence and Their Implications for Psychiatric Disorders

Importance

Many psychiatric disorders can be conceptualized as disorders of brain maturation during childhood and adolescence. Discovering the neurobiological underpinnings of brain maturation may elucidate molecular pathways of vulnerability and resilience to such disorders.

Objective

To investigate the underlying neurobiological mechanisms of age-associated cortical thinning during maturation and their implications for psychiatric disorders.

Design, Setting, and Participants

This multicohort analysis used data from 3 community-based studies. The Saguenay Youth Study provided data from 1024 adolescents who were recruited at a single site in Quebec, Canada. The IMAGEN cohort provided data from 1823 participants who were recruited in 8 European cities. The Brazil High Risk Cohort Study for the Development of Childhood Psychiatric Disorders provided data from 815 participants who were recruited in 2 Brazilian cities. Cortical thickness was estimated from the results of magnetic resonance imaging (MRI) scans, and age-associated cortical thinning was estimated in 34 cortical regions. Gene expression from the Allen Human Brain Atlas was aligned with the same regions. Similarities in the interregional profiles of gene expression and the profiles of age-associated cortical thinning were measured. The involvement of dendrites, dendritic spines, and myelin was tested using 3 gene panels. Enrichment for genes associated with psychiatric disorders was tested among the genes associated with thinning and their coexpression networks. Data analysis was conducted between March and October 2019.

Main Outcomes and Measures

MRI-derived estimates of age-associated cortical thinning and gene expression in 34 cortical regions.

Results

A total of 3596 individuals aged 9 to 21 years were included in this study. Of those, 1803 participants (50.1%) were female, and the mean (SD) age was 15.2 (2.6) years. Interregional profiles of age-associated cortical thinning were associated with interregional gradients in the expression of genes associated with dendrites, dendritic spines, and myelin; the variance in thinning explained by the gene panels across different points ranged from 0.45% to 10.55% for the dendrite panel, 0.00% to 9.98% for the spine panel, and 0.19% to 26.39% for the myelin panel. These genes and their coexpression networks were enriched for genes associated with several psychiatric disorders.

Conclusions and Relevance

In this study, genetic similarity between interregional variation in cortical thinning during maturation and multiple psychiatric disorders suggests overlapping molecular underpinnings. This finding adds to the understanding of the neurodevelopmental mechanisms of psychiatric disorders.

Associating Type 2 Diabetes Risk Factor Genes and FDG-PET Brain Metabolism in Normal Aging and Alzheimer's Disease

Background: Several studies have linked type 2 diabetes (T2D) to an increased risk of developing Alzheimer’s disease (AD). This has led to an interest in using antidiabetic treatments for the prevention of AD. However, the underlying mechanisms explaining the relationship between T2D and AD have not been completely elucidated.

Objective: Our objective was to examine cerebral ¹⁸F-fluorodeoxyglucose (FDG) uptake during normal aging and in AD patients in regions associated with diabetes genetic risk factor expression to highlight which genes may serve as potential targets for pharmaceutical intervention.

Methods: We calculated regional glucose metabolism differences in units of standardized uptake values (SUVR) for 386 cognitively healthy adults and 335 clinically probable AD patients. We then proceeded to extract gene-expression data from the publicly available Allen Human Brain Atlas (HBA) database. We used the nearest genes to 46 AD- and T2D-associated SNPs previously identified in the literature, and mapped their expression to the same 34 cortical regions in which we calculated SUVRs. SNPs with a donor consistency of 0.40 or greater were selected for further analysis. We evaluated the associations between SUVR and gene-expression across the brain.

Results: Of the 46 risk-factor genes, 15 were found to be significantly correlated with FDG-PET brain metabolism in healthy adults and probable AD patients after correction for multiple comparisons. Using multiple regression, we found that five genes explained a total of 72.5% of the SUVR variance across the healthy adult group regions, while four genes explained a total of 79.3% of the SUVR variance across the probable AD group regions. There were significant differences in whole-brain SUVR as a function of allele frequencies for two genes.

Conclusions: These results highlight the association between risk factor genes for T2D and regional glucose metabolism during both normal aging and in probable AD. Highlighted genes were associated with mitochondrial stability, vascular maintenance, and glucose intolerance. Pharmacological intervention of these pathways has the potential to improve glucose metabolism during normal again as well as in AD patients.

Individual Variation of Human Cortical Structure Is Established in the First Year of Life

BACKGROUND

Individual differences in cortical gray matter (GM) structure are associated with cognitive function and psychiatric disorders with developmental origins. Identifying when individual differences in cortical structure are established in childhood is critical for understanding the timing of abnormal cortical development associated with neuropsychiatric disorders.

METHODS

We studied the development of cortical GM and white matter volume, cortical thickness, and surface area using structural magnetic resonance imaging in two unique cohorts of singleton (121 male and 131 female) and twin (99 male and 83 female) children imaged longitudinally from birth to 6 years.

RESULTS

Cortical GM volume increases rapidly in the first year of life, with more gradual growth thereafter. Between ages 1 and 6 years, total surface area expands 29%, while average cortical thickness decreases about 3.5%. In both cohorts, a large portion of individual variation in cortical GM volume (81%-87%) and total surface area (73%-83%) at age 6 years is present by age 1 year. Regional heterogeneity of cortical thickness observed at age 6 is largely in place at age 1.

CONCLUSIONS

These findings indicate that individual differences in cortical GM structure are largely established by the end of the first year of life, following a period of rapid postnatal GM growth. This suggests that alterations in GM structure associated with psychiatric disorders with developmental origins may largely arise in the first year of life and that interventions to normalize or mitigate abnormal GM development may need to be targeted to very early childhood.

Thickness of the cerebral cortex shows positive association with blood levels of triacylglycerols carrying 18-carbon fatty acids

Perturbations in fatty acid (FA) metabolism as well as thinning of the cerebral cortex have been associated with cognitive decline in the elderly. Predominant FAs in the brain are docosahexaenoic acid (DHA) and arachidonic acid (ARA). Approximately 2-8% of esterified DHA and 3-5% of esterified ARA in the brain are replaced daily. DHA and ARA are derivatives of 18-carbon essential FAs, α-linolenic acid and linoleic acid, that must be imported into the brain from the circulation. In blood, FAs are primarily transported in triacylglycerols (TAGs) from which they can be released at the blood-brain-barrier and transported inside the brain. We show that circulating levels of TAGs carrying 18-carbon FAs are positively associated with cortical thickness in middle-aged adults. These associations are stronger in cortical regions with higher expression of genes regulating long-chain FA metabolism and cellular membranes, and cortical thickness in the same regions may be related to cognitive performance.

Maturation of the Human Cerebral Cortex During Adolescence: Myelin or Dendritic Arbor?

Previous in vivo studies revealed robust age-related variations in structural properties of the human cerebral cortex during adolescence. Neurobiology underlying these maturational phenomena is largely unknown. Here we employ a virtual-histology approach to gain insights into processes associated with inter-regional variations in cortical microstructure and its maturation, as indexed by magnetization transfer ratio (MTR). Inter-regional variations in MTR correlate with inter-regional variations in expression of genes specific to pyramidal cells (CA1) and ependymal cells; enrichment analyses indicate involvement of these genes in dendritic growth. On the other hand, inter-regional variations in the change of MTR during adolescence correlate with inter-regional profiles of oligodendrocyte-specific gene expression. Complemented by a quantitative hypothetical model of the contribution of surfaces associated with dendritic arbor (1631 m2) and myelin (48 m2), these findings suggest that MTR signals are driven mainly by macromolecules associated with dendritic arbor while maturational changes in the MTR signal are associated with myelination.

Cortical thickness and formal thought disorder in schizophrenia: An ultra high-field network-based morphometry study

BACKGROUND

Persistent formal thought disorder (FTD) is a core feature of schizophrenia. Recent cognitive and neuroimaging studies indicate a distinct mechanistic pathway underlying the persistent positive FTD (pFTD or disorganized thinking), though its structural determinants are still elusive. Using network-based cortical thickness estimates from ultra-high field 7-Tesla Magnetic Resonance Imaging (7T MRI), we investigated the structural correlates of pFTD.

METHODS

We obtained speech samples and 7T MRI anatomical scans from medicated clinically stable patients with schizophrenia (n = 19) and healthy controls (n = 20). Network-based morphometry was used to estimate the mean cortical thickness of 17 functional networks covering the entire cortical surface from each subject. We also quantified the vertexwise variability of thickness within each network to quantify the spatial coherence of the 17 networks, estimated patients vs. controls differences, and related the thickness of the affected networks to the severity of pFTD.

RESULTS

Patients had reduced thickness of the frontoparietal and default mode networks, and reduced spatial coherence affecting the salience and the frontoparietal control network. A higher burden of positive FTD related to reduced frontoparietal thickness and reduced spatial coherence of the salience network. The presence of positive FTD, but not its severity, related to the reduced thickness of the language network comprising of the superior temporal cortex.

CONCLUSIONS

These results suggest that cortical thickness of both cognitive control and language networks underlie the positive FTD in schizophrenia. The structural integrity of cognitive control networks is a critical determinant of the expressed severity of persistent FTD in schizophrenia.

Corticosteroids and Regional Variations in Thickness of the Human Cerebral Cortex across the Lifespan

Exposures to life stressors accumulate across the lifespan, with possible impact on brain health. Little is known, however, about the mechanisms mediating age-related changes in brain structure. We use a lifespan sample of participants (n = 21 251; 4-97 years) to investigate the relationship between the thickness of cerebral cortex and the expression of the glucocorticoid- and the mineralocorticoid-receptor genes (NR3C1 and NR3C2, respectively), obtained from the Allen Human Brain Atlas. In all participants, cortical thickness correlated negatively with the expression of both NR3C1 and NR3C2 across 34 cortical regions. The magnitude of this correlation varied across the lifespan. From childhood through early adulthood, the profile similarity (between NR3C1/NR3C2 expression and thickness) increased with age. Conversely, both profile similarities decreased with age in late life. These variations do not reflect age-related changes in NR3C1 and NR3C2 expression, as observed in 5 databases of gene expression in the human cerebral cortex (502 donors). Based on the co-expression of NR3C1 (and NR3C2) with genes specific to neural cell types, we determine the potential involvement of microglia, astrocytes, and CA1 pyramidal cells in mediating the relationship between corticosteroid exposure and cortical thickness. Therefore, corticosteroids may influence brain structure to a variable degree throughout life.

Individual-Level Identification of Gene Expression Associated with Volume Differences among Neocortical Areas

The human cerebral cortex is the source of many complex behaviors and is a vulnerable target of various neuropsychiatric disorders, but transcriptional profiles linked to cerebral cortical volume (CCV) differences across brain areas remain unknown. Here, we screened CCV-related genes using an across-sample spatial correlation analysis in 6 postmortem brains and then individually validated these correlations in 1091 subjects with different ages and ethnicities. We identified 62 genes whose transcriptional profiles were repeatedly associated with CCV in more than 90% of individuals. CCV-related genes were specifically expressed in neurons and in developmental periods from middle childhood to young adulthood, were enriched in ion channels and developmental processes, and showed significant overlap with genes linked to brain functional activity and mental disorders. The identified genes represent the conserved transcriptional architecture of the human cerebral cortex, suggesting a link between conserved gene transcription and neocortical structural properties.

Microglial regional heterogeneity and its role in the brain

Microglia have been recently shown to manifest a very interesting phenotypical heterogeneity across different regions in the mammalian central nervous system (CNS). However, the underlying mechanism and functional meaning of this phenomenon are currently unclear. Baseline diversities of adult microglia in their cell number, cellular and subcellular structures, molecular signature as well as relevant functions have been discovered. But recent transcriptomic studies using bulk RNAseq and single-cell RNAseq have produced conflicting results on region-specific signatures of microglia. It is highly speculative whether such spatial heterogeneity contributes to varying sensitivities of individual microglia to the same physiological and pathological signals in different CNS regions, and hence underlie their functional relevance for CNS disease development. This review aims to thoroughly summarize up-to-date knowledge on this specific topic and provide some insights on the potential underlying mechanisms, starting from microgliogenesis. Understanding regional heterogeneity of microglia in the context of their diverse neighboring neurons and other glia may provide an important clue for future development of innovative therapies for neuropsychiatric disorders.

Neural Transcription Correlates of Multimodal Cortical Phenotypes during Development

During development, cellular events such as cell proliferation, migration, and synaptogenesis determine the structural organization of the brain. These processes are driven in part by spatiotemporally regulated gene expression. We investigated how the genetic signatures of specific neural cell types shape cortical organization of the human brain throughout infancy and childhood. Using a transcriptional atlas and in vivo magnetic resonance imaging (MRI) data, we demonstrated time-dependent associations between the expression levels of neuronal and glial genes and cortical macro- and microstructure. Neonatal cortical phenotypes were associated with prenatal glial but not neuronal gene expression. These associations reflect cell migration and proliferation during fetal development. Childhood cortical phenotypes were associated with neuronal and astrocyte gene expression related to synaptic signaling processes, reflecting the refinement of cortical connections. These findings indicate that sequential developmental stages contribute to distinct MRI measures at different time points. This helps to bridge the gap between the genetic mechanisms driving cellular changes and widely used neuroimaging techniques.

Characterizing the Molecular Architecture of Cortical Regions Associated with High Educational Attainment in Older Individuals

Neuroimaging investigations have revealed interindividual variations in anatomy, metabolism, activity, and connectivity of specific cortical association areas through which years of education (YoE), as a common proxy of cognitive reserve, may operate in the face of age- or pathology-associated brain changes. However, the associated molecular properties of YoE-related brain regions and the biological pathways involved remain poorly understood. In the present study we first identified brain areas that showed an association between cortical thickness and YoE among 122 cognitively healthy older human individuals (87 female). We subsequently characterized molecular properties of these regions by studying brain-wide microarray measurements of regional gene expression. In accordance with previous studies, we observed that YoE were associated with higher cortical thickness in medial prefrontal, anterior cingulate, and orbitofrontal areas. Compared with the rest of the cortex, these regions exhibited a distinct gene expression profile characterized by relative upregulation of gene sets implicated in ionotropic and metabotropic neurotransmission as well as activation of immune response. Our genome-wide expression profile analysis of YoE-related brain regions points to distinct molecular pathways that may underlie a higher capacity for plastic changes in response to lifetime intellectual enrichment and potentially also a higher resilience to age-related pathologic brain changes.SIGNIFICANCE STATEMENT We combined a neuroimaging-based analysis with a transcriptome-wide gene expression approach to investigate the molecular-functional properties of cortical regions associated with educational attainment, as a commonly used proxy for cognitive reserve, in older individuals. The strongest association with education was observed in specific areas of the medial prefrontal cortex, and these areas exhibited a distinct gene expression profile characterized by relative upregulation of gene sets implicated in neurotransmission and immune responses. These findings complement previous neuroimaging studies in the field and point to novel biological pathways that may mediate the beneficial effects of high educational attainment on adaptability to cope with, or prevent, age-related brain changes. The identified genes and pathways now warrant further exploration in mechanistic studies.

A practical guide to linking brain-wide gene expression and neuroimaging data

The recent availability of comprehensive, brain-wide gene expression atlases such as the Allen Human Brain Atlas (AHBA) has opened new opportunities for understanding how spatial variations on molecular scale relate to the macroscopic neuroimaging phenotypes. A rapidly growing body of literature is demonstrating relationships between gene expression and diverse properties of brain structure and function, but approaches for combining expression atlas data with neuroimaging are highly inconsistent, with substantial variations in how the expression data are processed. The degree to which these methodological variations affect findings is unclear. Here, we outline a seven-step analysis pipeline for relating brain-wide transcriptomic and neuroimaging data and compare how different processing choices influence the resulting data. We suggest that studies using the AHBA should work towards a unified data processing pipeline to ensure consistent and reproducible results in this burgeoning field.

Transcriptomic characterization of MRI contrast with focus on the T1-w/T2-w ratio in the cerebral cortex

Magnetic resonance (MR) images of the brain are of immense clinical and research utility. At the atomic and subatomic levels, the sources of MR signals are well understood. However, we lack a comprehensive understanding of the macromolecular correlates of MR signal contrast. To address this gap, we used genome-wide measurements to correlate gene expression with MR signal intensity across the cerebral cortex in the Allen Human Brain Atlas (AHBA). We focused on the ratio of T1-weighted and T2-weighted intensities (T1-w/T2-w ratio image), which is considered to be a useful proxy for myelin content. As expected, we found enrichment of positive correlations between myelin-associated genes and the ratio image, supporting its use as a myelin marker. Genome-wide, there was an association with protein mass, with genes coding for heavier proteins expressed in regions with high T1-w/T2-w values. Oligodendrocyte gene markers were strongly correlated with the T1-w/T2-w ratio, but this was not driven by myelin-associated genes. Mitochondrial genes exhibit the strongest relationship, showing higher expression in regions with low T1-w/T2-w ratio. This may be due to the pH gradient in mitochondria as genes up-regulated by pH in the brain were also highly correlated with the ratio. While we corroborate associations with myelin and synaptic plasticity, differences in the T1-w/T2-w ratio across the cortex are more strongly linked to molecule size, oligodendrocyte markers, mitochondria, and pH. We evaluate correlations between AHBA transcriptomic measurements and a group averaged T1-w/T2-w ratio image, showing agreement with in-sample results. Expanding our analysis to the whole brain results in strong positive T1-w/T2-w correlations for immune system, inflammatory disease, and microglia marker genes. Genes with negative correlations were enriched for neuron markers and synaptic plasticity genes. Lastly, our findings are similar when performed on T1-w or inverted T2-w intensities alone. These results provide a molecular characterization of MR contrast that will aid interpretation of future MR studies of the brain.