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

Obesity and the cerebral cortex: Underlying neurobiology in mice and humans

Obesity is a major modifiable risk factor for Alzheimer's disease (AD), characterized by progressive atrophy of the cerebral cortex. The neurobiology of obesity contributions to AD is poorly understood. Here we show with in vivo MRI that diet-induced obesity decreases cortical volume in mice, and that higher body adiposity associates with lower cortical volume in humans. Single-nuclei transcriptomics of the mouse cortex reveals that dietary obesity promotes an array of neuron-adverse transcriptional dysregulations, which are mediated by an interplay of excitatory neurons and glial cells, and which involve microglial activation and lowered neuronal capacity for neuritogenesis and maintenance of membrane potential. The transcriptional dysregulations of microglia, more than of other cell types, are like those in AD, as assessed with single-nuclei cortical transcriptomics in a mouse model of AD and two sets of human donors with the disease. Serial two-photon tomography of microglia demonstrates microgliosis throughout the mouse cortex. The spatial pattern of adiposity-cortical volume associations in human cohorts interrogated together with in silico bulk and single-nucleus transcriptomic data from the human cortex implicated microglia (along with other glial cells and subtypes of excitatory neurons), and it correlated positively with the spatial profile of cortical atrophy in patients with mild cognitive impairment and AD. Thus, multi-cell neuron-adverse dysregulations likely contribute to the loss of cortical tissue in obesity. The dysregulations of microglia may be pivotal to the obesity-related risk of AD.

Unraveling the molecular relevance of brain phenotypes: A comparative analysis of null models and test statistics

Correlating transcriptional profiles with imaging-derived phenotypes has the potential to reveal possible molecular architectures associated with cognitive functions, brain development and disorders. Competitive null models built by resampling genes and self-contained null models built by spinning brain regions, along with varying test statistics, have been used to determine the significance of transcriptional associations. However, there has been no systematic evaluation of their performance in imaging transcriptomics analyses. Here, we evaluated the performance of eight different test statistics (mean, mean absolute value, mean squared value, max mean, median, Kolmogorov-Smirnov (KS), Weighted KS and the number of significant correlations) in both competitive null models and self-contained null models. Simulated brain maps (n = 1,000) and gene sets (n = 500) were used to calculate the probability of significance (Psig) for each statistical test. Our results suggested that competitive null models may result in false positive results driven by co-expression within gene sets. Furthermore, we demonstrated that the self-contained null models may fail to account for distribution characteristics (e.g., bimodality) of correlations between all available genes and brain phenotypes, leading to false positives. These two confounding factors interacted differently with test statistics, resulting in varying outcomes. Specifically, the sign-sensitive test statistics (i.e., mean, median, KS, Weighted KS) were influenced by co-expression bias in the competitive null models, while median and sign-insensitive test statistics were sensitive to the bimodality bias in the self-contained null models. Additionally, KS-based statistics produced conservative results in the self-contained null models, which increased the risk of false negatives. Comprehensive supplementary analyses with various configurations, including realistic scenarios, supported the results. These findings suggest utilizing sign-insensitive test statistics such as mean absolute value, max mean in the competitive null models and the mean as the test statistic for the self-contained null models. Additionally, adopting the confounder-matched (e.g., coexpression-matched) null models as an alternative to standard null models can be a viable strategy. Overall, the present study offers insights into the selection of statistical tests for imaging transcriptomics studies, highlighting areas for further investigation and refinement in the evaluation of novel and commonly used tests.

Regional patterns of human cortex development colocalize with underlying neurobiology

Human brain morphology undergoes complex changes over the lifespan. Despite recent progress in tracking brain development via normative models, current knowledge of underlying biological mechanisms is highly limited. We demonstrate that human cerebral cortex development and aging trajectories unfold along patterns of molecular and cellular brain organization, traceable from population-level to individual developmental trajectories. During childhood and adolescence, cortex-wide spatial distributions of dopaminergic receptors, inhibitory neurons, glial cell populations, and brain-metabolic features explain up to 50% of variance associated with a lifespan model of regional cortical thickness trajectories. In contrast, modeled cortical change patterns during adulthood are best explained by cholinergic and glutamatergic neurotransmitter receptor and transporter distributions. These relationships are supported by developmental gene expression trajectories and translate to individual longitudinal data from over 8,000 adolescents, explaining up to 59% of developmental change at cohort- and 18% at single-subject level. Integrating neurobiological brain atlases with normative modeling and population neuroimaging provides a biologically meaningful path to understand brain development and aging in living humans.

Differences in the neural correlates of schizophrenia with positive and negative formal thought disorder in patients with schizophrenia in the ENIGMA dataset

Formal thought disorder (FTD) is a clinical key factor in schizophrenia, but the neurobiological underpinnings remain unclear. In particular, the relationship between FTD symptom dimensions and patterns of regional brain volume loss in schizophrenia remains to be established in large cohorts. Even less is known about the cellular basis of FTD. Our study addresses these major obstacles by enrolling a large multi-site cohort acquired by the ENIGMA Schizophrenia Working Group (752 schizophrenia patients and 1256 controls), to unravel the neuroanatomy of FTD in schizophrenia and using virtual histology tools on implicated brain regions to investigate the cellular basis. Based on the findings of previous clinical and neuroimaging studies, we decided to separately explore positive, negative and total formal thought disorder. We used virtual histology tools to relate brain structural changes associated with FTD to cellular distributions in cortical regions. We identified distinct neural networks positive and negative FTD. Both networks encompassed fronto-occipito-amygdalar brain regions, but positive and negative FTD demonstrated a dissociation: negative FTD showed a relative sparing of orbitofrontal cortical thickness, while positive FTD also affected lateral temporal cortices. Virtual histology identified distinct transcriptomic fingerprints associated for both symptom dimensions. Negative FTD was linked to neuronal and astrocyte fingerprints, while positive FTD also showed associations with microglial cell types. These results provide an important step towards linking FTD to brain structural changes and their cellular underpinnings, providing an avenue for a better mechanistic understanding of this syndrome.

Population Neuroscience: Principles and Advances

In population neuroscience, three disciplines come together to advance our knowledge of factors that shape the human brain: neuroscience, genetics, and epidemiology (Paus, Human Brain Mapping 31:891-903, 2010). Here, I will come back to some of the background material reviewed in more detail in our previous book (Paus, Population Neuroscience, 2013), followed by a brief overview of current advances and challenges faced by this integrative approach.

Cortical growth from infancy to adolescence in preterm and term-born children

Early life experiences can exert a significant influence on cortical and cognitive development. Very preterm birth exposes infants to several adverse environmental factors during hospital admission, which affect cortical architecture. However, the subsequent consequence of very preterm birth on cortical growth from infancy to adolescence has never been defined; despite knowledge of critical periods during childhood for establishment of cortical networks. Our aims were to: chart typical longitudinal cortical development and sex differences in cortical development from birth to adolescence in healthy term-born children; estimate differences in cortical development between children born at term and very preterm; and estimate differences in cortical development between children with normal and impaired cognition in adolescence. This longitudinal cohort study included children born at term (≥37 weeks' gestation) and very preterm (<30 weeks' gestation) with MRI scans at ages 0, 7 and 13 years (n = 66 term-born participants comprising 34 with one scan, 18 with two scans and 14 with three scans; n = 201 very preterm participants comprising 56 with one scan, 88 with two scans and 57 with three scans). Cognitive assessments were performed at age 13 years. Cortical surface reconstruction and parcellation were performed with state-of-the-art, equivalent MRI analysis pipelines for all time points, resulting in longitudinal cortical volume, surface area and thickness measurements for 62 cortical regions. Developmental trajectories for each region were modelled in term-born children, contrasted between children born at term and very preterm, and contrasted between all children with normal and impaired cognition. In typically developing term-born children, we documented anticipated patterns of rapidly increasing cortical volume, area and thickness in early childhood, followed by more subtle changes in later childhood, with smaller cortical size in females than males. In contrast, children born very preterm exhibited increasingly reduced cortical volumes, relative to term-born children, particularly during ages 0-7 years in temporal cortical regions. This reduction in cortical volume in children born very preterm was largely driven by increasingly reduced cortical thickness rather than area. This resulted in amplified cortical volume and thickness reductions by age 13 years in individuals born very preterm. Alterations in cortical thickness development were found in children with impaired language and memory. This study shows that the neurobiological impact of very preterm birth on cortical growth is amplified from infancy to adolescence. These data further inform the long-lasting impact on cortical development from very preterm birth, providing broader insights into neurodevelopmental consequences of early life experiences.

General and specific patterns of cortical gene expression as spatial correlates of complex cognitive functioning

Gene expression varies across the brain. This spatial patterning denotes specialised support for particular brain functions. However, the way that a given gene's expression fluctuates across the brain may be governed by general rules. Quantifying patterns of spatial covariation across genes would offer insights into the molecular characteristics of brain areas supporting, for example, complex cognitive functions. Here, we use principal component analysis to separate general and unique gene regulatory associations with cortical substrates of cognition. We find that the region-to-region variation in cortical expression profiles of 8235 genes covaries across two major principal components: gene ontology analysis suggests these dimensions are characterised by downregulation and upregulation of cell-signalling/modification and transcription factors. We validate these patterns out-of-sample and across different data processing choices. Brain regions more strongly implicated in general cognitive functioning (g; 3 cohorts, total meta-analytic N = 39,519) tend to be more balanced between downregulation and upregulation of both major components (indicated by regional component scores). We then identify a further 29 genes as candidate cortical spatial correlates of g, beyond the patterning of the two major components (|β| range = 0.18 to 0.53). Many of these genes have been previously associated with clinical neurodegenerative and psychiatric disorders, or with other health-related phenotypes. The results provide insights into the cortical organisation of gene expression and its association with individual differences in cognitive functioning.

Structural connectome architecture shapes the maturation of cortical morphology from childhood to adolescence

Cortical thinning is an important hallmark of the maturation of brain morphology during childhood and adolescence. However, the connectome-based wiring mechanism that underlies cortical maturation remains unclear. Here, we show cortical thinning patterns primarily located in the lateral frontal and parietal heteromodal nodes during childhood and adolescence, which are structurally constrained by white matter network architecture and are particularly represented using a network-based diffusion model. Furthermore, connectome-based constraints are regionally heterogeneous, with the largest constraints residing in frontoparietal nodes, and are associated with gene expression signatures of microstructural neurodevelopmental events. These results are highly reproducible in another independent dataset. These findings advance our understanding of network-level mechanisms and the associated genetic basis that underlies the maturational process of cortical morphology during childhood and adolescence.

Age-associated cortical similarity networks correlate with cell type-specific transcriptional signatures

Human brain structure shows heterogeneous patterns of change across adults aging and is associated with cognition. However, the relationship between cortical structural changes during aging and gene transcription signatures remains unclear. Here, using structural magnetic resonance imaging data of two separate cohorts of healthy participants from the Cambridge Centre for Aging and Neuroscience (n = 454, 18-87 years) and Dallas Lifespan Brain Study (n = 304, 20-89 years) and a transcriptome dataset, we investigated the link between cortical morphometric similarity network and brain-wide gene transcription. In two cohorts, we found reproducible morphometric similarity network change patterns of decreased morphological similarity with age in cognitive related areas (mainly located in superior frontal and temporal cortices), and increased morphological similarity in sensorimotor related areas (postcentral and lateral occipital cortices). Changes in morphometric similarity network showed significant spatial correlation with the expression of age-related genes that enriched to synaptic-related biological processes, synaptic abnormalities likely accounting for cognitive decline. Transcription changes in astrocytes, microglia, and neuronal cells interpreted most of the age-related morphometric similarity network changes, which suggest potential intervention and therapeutic targets for cognitive decline. Taken together, by linking gene transcription signatures to cortical morphometric similarity network, our findings might provide molecular and cellular substrates for cortical structural changes related to cognitive decline across adults aging.

Positive or negative environmental modulations on human brain development: the morpho-functional outcomes of music training or stress

In the last couple of decades, the study of human living brain has benefitted of neuroimaging and non-invasive electrophysiological techniques, which are particularly valuable during development. A number of studies allowed to trace the usual stages leading from pregnancy to adult age, and relate them to functional and behavioral measurements. It was also possible to explore the effects of some interventions, behavioral or not, showing that the commonly followed pathway to adulthood may be steered by external interventions. These events may result in behavioral modifications but also in structural changes, in some cases limiting plasticity or extending/modifying critical periods. In this review, we outline the healthy human brain development in the absence of major issues or diseases. Then, the effects of negative (different stressors) and positive (music training) environmental stimuli on brain and behavioral development is depicted. Hence, it may be concluded that the typical development follows a course strictly dependent from environmental inputs, and that external intervention can be designed to positively counteract negative influences, particularly at young ages. We also focus on the social aspect of development, which starts in utero and continues after birth by building social relationships. This poses a great responsibility in handling children education and healthcare politics, pointing to social accountability for the responsible development of each child.

Hemispheric asymmetry in cortical thinning reflects intrinsic organization of the neurotransmitter systems and homotopic functional connectivity

Hemispheric lateralization and its origins have been of great interest in neuroscience for over a century. The left-right asymmetry in cortical thickness may stem from differential maturation of the cerebral cortex in the two hemispheres. Here, we investigated the spatial pattern of hemispheric differences in cortical thinning during adolescence, and its relationship with the density of neurotransmitter receptors and homotopic functional connectivity. Using longitudinal data from IMAGEN study (N = 532), we found that many cortical regions in the frontal and temporal lobes thinned more in the right hemisphere than in the left. Conversely, several regions in the occipital and parietal lobes thinned less in the right (vs. left) hemisphere. We then revealed that regions thinning more in the right (vs. left) hemispheres had higher density of neurotransmitter receptors and transporters in the right (vs. left) side. Moreover, the hemispheric differences in cortical thinning were predicted by homotopic functional connectivity. Specifically, regions with stronger homotopic functional connectivity showed a more symmetrical rate of cortical thinning between the left and right hemispheres, compared with regions with weaker homotopic functional connectivity. Based on these findings, we suggest that the typical patterns of hemispheric differences in cortical thinning may reflect the intrinsic organization of the neurotransmitter systems and related patterns of homotopic functional connectivity.

Network connectivity and local transcriptomic vulnerability underpin cortical atrophy progression in Parkinson's disease

Parkinson's disease pathology is hypothesized to spread through the brain via axonal connections between regions and is further modulated by local vulnerabilities within those regions. The resulting changes to brain morphology have previously been demonstrated in both prodromal and de novo Parkinson's disease patients. However, it remains unclear whether the pattern of atrophy progression in Parkinson's disease over time is similarly explained by network-based spreading and local vulnerability. We address this gap by mapping the trajectory of cortical atrophy rates in a large, multi-centre cohort of Parkinson's disease patients and relate this atrophy progression pattern to network architecture and gene expression profiles. Across 4-year follow-up visits, increased atrophy rates were observed in posterior, temporal, and superior frontal cortices. We demonstrated that this progression pattern was shaped by network connectivity. Regional atrophy rates were strongly related to atrophy rates across structurally and functionally connected regions. We also found that atrophy progression was associated with specific gene expression profiles. The genes whose spatial distribution in the brain was most related to atrophy rate were those enriched for mitochondrial and metabolic function. Taken together, our findings demonstrate that both global and local brain features influence vulnerability to neurodegeneration in Parkinson's disease.

Neural Correlates of Positive and Negative Formal Thought Disorder in Individuals with Schizophrenia: An ENIGMA Schizophrenia Working Group Study

Formal thought disorder (FTD) is a key clinical factor in schizophrenia, but the neurobiological underpinnings remain unclear. In particular, relationship between FTD symptom dimensions and patterns of regional brain volume deficiencies in schizophrenia remain to be established in large cohorts. Even less is known about the cellular basis of FTD. Our study addresses these major obstacles based on a large multi-site cohort through the ENIGMA Schizophrenia Working Group (752 individuals with schizophrenia and 1256 controls), to unravel the neuroanatomy of positive, negative and total FTD in schizophrenia and their cellular bases. We used virtual histology tools to relate brain structural changes associated with FTD to cellular distributions in cortical regions. We identified distinct neural networks for positive and negative FTD. Both networks encompassed fronto-occipito-amygdalar brain regions, but negative FTD showed a relative sparing of orbitofrontal cortical thickness, while positive FTD also affected lateral temporal cortices. Virtual histology identified distinct transcriptomic fingerprints associated for both symptom dimensions. Negative FTD was linked to neuronal and astrocyte fingerprints, while positive FTD was also linked to microglial cell types. These findings relate different dimensions of FTD to distinct brain structural changes and their cellular underpinnings, improve our mechanistic understanding of these key psychotic symptoms.

WNT pathway in focal cortical dysplasia compared to perilesional nonlesional tissue in refractory epilepsies

BACKGROUND

Focal cortical dysplasia (FCD) is a malformation of cortical development that causes medical refractory seizures, and one of the main treatments may be surgical resection of the affected area of the brain. People affected by FCD may present with seizures of variable severity since childhood. Despite many medical treatments available, only surgery can offer cure. The pathophysiology of the disease is not yet understood; however, it is known that several gene alterations may play a role. The WNT/β-catenin pathway is closely related to the control and balance of cell proliferation and differentiation in the central nervous system. The aim of this study was to explore genes related to the WNT/β-catenin pathway in lesional and perilesional brain tissue in patients with FCD type II.

METHODS

Dysplastic and perilesional tissue from the primary dysplastic lesion of patients with FCD type IIa were obtained from two patients who underwent surgical treatment. The analysis of the relative expression of genes was performed by a qRT-PCR array (super array) containing 84 genes related to the WNT pathway.

RESULTS

Our results suggest the existence of molecular alteration in some genes of the WNT pathway in tissue with dysplastic lesions and of perilesional tissue. We call this tissue of normal-appearing adjacent cortex (NAAC). Of all genes analyzed, a large number of genes show similar behavior between injured, perilesional and control tissues. However, some genes have similar characteristics between the perilesional and lesional tissue and are different from the control brain tissue, presenting the perilesional tissue as a molecularly altered material.

CONCLUSION

Our results suggest that the perilesional area after surgical resection of tissue with cortical dysplasia presents molecular changes that may play a role in the recurrence of seizures in these patients. The perilesional tissue should receive expanded attention beyond the somatic mutations described and associated with FCD, such as mTOR, for example, to new signaling pathways that may play a crucial role in seizure recurrence.

General and specific patterns of cortical gene expression as spatial correlates of complex cognitive functioning

Gene expression varies across the brain. This spatial patterning denotes specialised support for particular brain functions. However, the way that a given gene's expression fluctuates across the brain may be governed by general rules. Quantifying patterns of spatial covariation across genes would offer insights into the molecular characteristics of brain areas supporting, for example, complex cognitive functions. Here, we use principal component analysis to separate general and unique gene regulatory associations with cortical substrates of cognition. We find that the region-to-region variation in cortical expression profiles of 8235 genes covaries across two major principal components : gene ontology analysis suggests these dimensions are characterised by downregulation and upregulation of cell-signalling/modification and transcription factors. We validate these patterns out-of-sample and across different data processing choices. Brain regions more strongly implicated in general cognitive functioning (g; 3 cohorts, total meta-analytic N = 39,519) tend to be more balanced between downregulation and upregulation of both major components (indicated by regional component scores). We then identify a further 41 genes as candidate cortical spatial correlates of g, beyond the patterning of the two major components (|β| range = 0.15 to 0.53). Many of these genes have been previously associated with clinical neurodegenerative and psychiatric disorders, or with other health-related phenotypes. The results provide insights into the cortical organisation of gene expression and its association with individual differences in cognitive functioning.

Case-control virtual histology elucidates cell types associated with cortical thickness differences in Alzheimer's disease

Many neuropsychiatric disorders are characterised by altered cortical thickness, but the cell types underlying these changes remain largely unknown. Virtual histology (VH) approaches map regional patterns of gene expression with regional patterns of MRI-derived phenotypes, such as cortical thickness, to identify cell types associated with case-control differences in those MRI measures. However, this method does not incorporate valuable information of case-control differences in cell type abundance. We developed a novel method, termed case-control virtual histology (CCVH), and applied it to Alzheimer's disease (AD) and dementia cohorts. Leveraging a multi-region gene expression dataset of AD cases (n = 40) and controls (n = 20), we quantified AD case-control differential expression of cell type-specific markers across 13 brain regions. We then correlated these expression effects with MRI-derived AD case-control cortical thickness differences across the same regions. Cell types with spatially concordant AD-related effects were identified through resampling marker correlation coefficients. Among regions thinner in AD, gene expression patterns identified by CCVH suggested fewer excitatory and inhibitory neurons, and greater proportions of astrocytes, microglia, oligodendrocytes, oligodendrocyte precursor cells, and endothelial cells in AD cases vs. controls. In contrast, original VH identified expression patterns suggesting that excitatory but not inhibitory neuron abundance was associated with thinner cortex in AD, despite the fact that both types of neurons are known to be lost in the disorder. Compared to original VH, cell types identified through CCVH are more likely to directly underlie cortical thickness differences in AD. Sensitivity analyses suggest our results are largely robust to specific analysis choices, like numbers of cell type-specific marker genes used and background gene sets used to construct null models. As more multi-region brain expression datasets become available, CCVH will be useful for identifying the cellular correlates of cortical thickness across neuropsychiatric illnesses.

Structural neuroimaging phenotypes and associated molecular and genomic underpinnings in autism: a review

Autism has been associated with differences in the developmental trajectories of multiple neuroanatomical features, including cortical thickness, surface area, cortical volume, measures of gyrification, and the gray-white matter tissue contrast. These neuroimaging features have been proposed as intermediate phenotypes on the gradient from genomic variation to behavioral symptoms. Hence, examining what these proxy markers represent, i.e., disentangling their associated molecular and genomic underpinnings, could provide crucial insights into the etiology and pathophysiology of autism. In line with this, an increasing number of studies are exploring the association between neuroanatomical, cellular/molecular, and (epi)genetic variation in autism, both indirectly and directly in vivo and across age. In this review, we aim to summarize the existing literature in autism (and neurotypicals) to chart a putative pathway from (i) imaging-derived neuroanatomical cortical phenotypes to (ii) underlying (neuropathological) biological processes, and (iii) associated genomic variation.

Neural Correlates of Positive and Negative Formal Thought Disorder in Individuals with Schizophrenia: An ENIGMA Schizophrenia Working Group Study

Formal thought disorder (FTD) is a key clinical factor in schizophrenia, but the neurobiological underpinnings remain unclear. In particular, relationship between FTD symptom dimensions and patterns of regional brain volume deficiencies in schizophrenia remain to be established in large cohorts. Even less is known about the cellular basis of FTD. Our study addresses these major obstacles based on a large multi-site cohort through the ENIGMA Schizophrenia Working Group (752 individuals with schizophrenia and 1256 controls), to unravel the neuroanatomy of positive, negative and total FTD in schizophrenia and their cellular bases. We used virtual histology tools to relate brain structural changes associated with FTD to cellular distributions in cortical regions. We identified distinct neural networks for positive and negative FTD. Both networks encompassed fronto-occipito-amygdalar brain regions, but negative FTD showed a relative sparing of orbitofrontal cortical thickness, while positive FTD also affected lateral temporal cortices. Virtual histology identified distinct transcriptomic fingerprints associated for both symptom dimensions. Negative FTD was linked to neuronal and astrocyte fingerprints, while positive FTD was also linked to microglial cell types. These findings relate different dimensions of FTD to distinct brain structural changes and their cellular underpinnings, improve our mechanistic understanding of these key psychotic symptoms.

Gene transcriptional expression of cortical thinning during childhood and adolescence

The cognitive and behavioral development of children and adolescents is closely related to the maturation of brain morphology. Although the trajectory of brain development has been depicted in detail, the underlying biological mechanism of normal cortical morphological development in childhood and adolescence remains unclear. By combining the Allen Human Brain Atlas dataset with two single-site magnetic resonance imaging data including 427 and 733 subjects from China and the United States, respectively, we performed partial least squares regression and enrichment analysis to explore the relationship between the gene transcriptional expression and the development of cortical thickness in childhood and adolescence. We found that the spatial model of normal cortical thinning during childhood and adolescence is associated with genes expressed predominantly in astrocytes, microglia, excitatory and inhibitory neurons. Top cortical development-related genes are enriched for energy-related and DNA-related terms and are associated with psychological and cognitive disorders. Interestingly, there is a great deal of similarity between the findings derived from the two single-site datasets. This fills the gap between early cortical development and transcriptomes, which promotes an integrative understanding of the potential biological neural mechanisms.

A biologically informed polygenic score of neuronal plasticity moderates the association between cognitive aptitudes and cortical thickness in adolescents

Although many studies of the adolescent brain identified positive associations between cognitive abilities and cortical thickness, little is known about mechanisms underlying such brain-behavior relationships. With experience-induced plasticity playing an important role in shaping the cerebral cortex throughout life, it is likely that some of the inter-individual variations in cortical thickness could be explained by genetic variations in relevant molecular processes, as indexed by a polygenic score of neuronal plasticity (PGS-NP). Here, we studied associations between PGS-NP, cognitive abilities, and thickness of the cerebral cortex, estimated from magnetic resonance images, in the Saguenay Youth Study (SYS, 533 females, 496 males: age=15.0 ± 1.8 years of age; cross-sectional), and the IMAGEN Study (566 females, 556 males; between 14 and 19 years; longitudinal). Using Gene Ontology, we first identified 199 genes implicated in neuronal plasticity, which mapped to 155,600 single nucleotide polymorphisms (SNPs). Second, we estimated their effect sizes from an educational attainment meta-GWAS to build a PGS-NP. Third, we examined a possible moderating role of PGS-NP in the relationship between performance intelligence quotient (PIQ), and its subtests, and the thickness of 34 cortical regions. In SYS, we observed a significant interaction between PGS-NP and object assembly vis-à-vis thickness in male adolescents (p = 0.026). A median-split analysis showed that, in males with a 'high' PGS-NP, stronger associations between object assembly and thickness were found in regions with larger age-related changes in thickness (r = 0.55, p = 0.00075). Although the interaction between PIQ and PGS-NP was non-significant (p = 0.064), we performed a similar median-split analysis. Again, in the high PGS-NP males, positive associations between PIQ and thickness were observed in regions with larger age-related changes in thickness (r = 0.40, p = 0.018). In the IMAGEN cohort, we did not replicate the first set of results (interaction between PGS-NP and cognitive abilities via-a-vis cortical thickness) while we did observe the same relationship between the brain-behaviour relationship and (longitudinal) changes in cortical thickness (Matrix reasoning: r = 0.63, p = 6.5e-05). No statistically significant results were observed in female adolescents in either cohort. Overall, these cross-sectional and longitudinal results suggest that molecular mechanisms involved in neuronal plasticity may contribute to inter-individual variations of cortical thickness related to cognitive abilities during adolescence in a sex-specific manner.

The neurobiology of autism spectrum disorder as it relates to twice exceptionality

There is a long-standing association between exceptional cognitive abilities of various sorts and neuropsychiatric illness, but it has historically largely been investigated in an exploratory and non-systematic way. One group in which this association has been investigated with more rigor is in subjects who have been identified as twice exceptional; an educational term describing subjects who are both gifted and diagnosed with a neuropsychiatric disorder. This term covers multiple conditions, but is of specific interest in particular in the study of autism spectrum disorder. Recent findings have led to the development of a hypothesis that a certain degree of the neurobiology associated with autism might even be advantageous for individuals and could lead to high giftedness, while becoming disadvantageous, once a certain threshold is surpassed. In this model, the same neurobiological mechanisms confer an increasing advantage up to a certain threshold, but become pathological past that point. Twice-exceptional individuals would be exactly at the inflection point, being highly gifted, but also symptomatic at the same time. Here, we review how existing neuroimaging literature on autism spectrum disorder can inform research on twice exceptionality specifically. We propose to study key neural networks with a robust implication in ASD to identify the neurobiology underlying twice-exceptionality. A better understanding of the neural mechanisms of twice exceptionality should help to better understand resilience and vulnerability to neurodevelopmental disorders and to. further support affected individuals.

Mapping Human Brain Pathways: Challenges and Opportunities in the Integration of Scales

The human brain is composed of a complex web of pathways. Diffusion magnetic resonance (MR) tractography is a neuroimaging technique that relies on the principle of diffusion to reconstruct brain pathways. Its tractography is broadly applicable to a range of problems as it is amenable for study in individuals of any age and from any species. However, it is well known that this technique can generate biologically implausible pathways, especially in regions of the brain where multiple fibers cross. This review highlights potential misconnections in two cortico-cortical association pathways with a focus on the aslant tract and inferior frontal occipital fasciculus. The lack of alternative methods to validate observations from diffusion MR tractography means there is a need to develop new integrative approaches to trace human brain pathways. This review discusses integrative approaches in neuroimaging, anatomical, and transcriptional variation as having much potential to trace the evolution of human brain pathways.

Brain-wide and cell-specific transcriptomic insights into MRI-derived cortical morphology in macaque monkeys

Integrative analyses of transcriptomic and neuroimaging data have generated a wealth of information about biological pathways underlying regional variability in imaging-derived brain phenotypes in humans, but rarely in nonhuman primates due to the lack of a comprehensive anatomically-defined atlas of brain transcriptomics. Here we generate complementary bulk RNA-sequencing dataset of 819 samples from 110 brain regions and single-nucleus RNA-sequencing dataset, and neuroimaging data from 162 cynomolgus macaques, to examine the link between brain-wide gene expression and regional variation in morphometry. We not only observe global/regional expression profiles of macaque brain comparable to human but unravel a dorsolateral-ventromedial gradient of gene assemblies within the primate frontal lobe. Furthermore, we identify a set of 971 protein-coding and 34 non-coding genes consistently associated with cortical thickness, specially enriched for neurons and oligodendrocytes. These data provide a unique resource to investigate nonhuman primate models of human diseases and probe cross-species evolutionary mechanisms.

Tracking Development of Connectivity in the Human Brain: Axons and Dendrites

The neuron doctrine laid the foundation for our current thinking about the structural and functional organization of the human brain. With the basic units of the nervous system-neurons-being physically separate, their connectivity relies on the conduction of action potentials in axons and their transmission across the synaptic cleft to the dendrites of other neurons. This study reviews available ex vivo data about the cellular composition of the human cerebral cortex, focusing on axons and dendrites, to conceptualize biological sources of signals detected in vivo with magnetic resonance imaging. To bridge the gap between ex vivo and in vivo observations, I then explain the basic principles of virtual histology, an approach that integrates spatially cell- or process-specific transcriptomic data with magnetic resonance signals to facilitate their neurobiological interpretation. Finally, I provide an overview of the initial insights gained in this manner in studies of brain development and maturation, in both health and disease.

Toward Best Practices for Imaging Transcriptomics of the Human Brain

Modern brainwide transcriptional atlases provide unprecedented opportunities for investigating the molecular correlates of brain organization, as quantified using noninvasive neuroimaging. However, integrating neuroimaging data with transcriptomic measures is not straightforward, and careful consideration is required to make valid inferences. In this article, we review recent work exploring how various methodological choices affect 3 main phases of imaging transcriptomic analyses, including 1) processing of transcriptional atlas data; 2) relating transcriptional measures to independently derived neuroimaging phenotypes; and 3) evaluating the functional implications of identified associations through gene enrichment analyses. Our aim is to facilitate the development of standardized and reproducible approaches for this rapidly growing field. We identify sources of methodological variability, key choices that can affect findings, and considerations for mitigating false positive and/or spurious results. Finally, we provide an overview of freely available open-source toolboxes implementing current best-practice procedures across all 3 analysis phases.

Virtual histology of morphometric similarity network after risperidone monotherapy and imaging-epigenetic biomarkers for treatment response in first-episode schizophrenia

BACKGROUND

Antipsychotic treatment has been conceived to alter brain connectivity, but it is unclear how the changes of network phenotypes relate to the underlying transcriptomics. Given DNA methylation (DNAm) may alter transcriptional levels, we further integrated an imaging-transcriptomic-epigenetic analysis to explore multi-omics treatment response biomarkers.

METHODS

Forty-two treatment-naive first-episode schizophrenia patients were scanned by TI weighted (T1W) imaging and DTI before and after 8-week risperidone monotherapy, and their peripheral blood genomic DNAm values were examined in parallel with MRI scanning. Morphometric similarity network (MSN) quantified with DTI and T1W data were used as a marker of treatment-related alterations in interareal cortical connectivity. We utilized partial least squares (PLS) to examine spatial associations between treatment-related MSN variations and cortical transcriptomic data obtained from the Allen Human Brain Atlas.

RESULTS

Longitudinal MSN alterations were related to treatment response on cognitive function and general psychopathology symptoms, while DNAm values of 59 PLS1 genes were on negative and positive symptoms. Virtual-histology transcriptomic analysis linked the MSN alterations with the neurobiological, cellular and metabolic pathways or processes, and assigned MSN-related genes to multiple cell types, specifying neurons and glial cells as contributing most to the transcriptomic associations of longitudinal changes in MSN.

CONCLUSIONS

We firstly reveal how brain-wide transcriptional levels and cell classes capture molecularly validated cortical connectivity alterations after antipsychotic treatment. Our findings represent a vital step towards the exploration of treatment response biomarkers on the basis of multiple omics rather than a single omics type as a strategy for advancing precise care.

Cortical profiles of numerous psychiatric disorders and normal development share a common pattern

The neurobiological bases of the association between development and psychopathology remain poorly understood. Here, we identify a shared spatial pattern of cortical thickness (CT) in normative development and several psychiatric and neurological disorders. Principal component analysis (PCA) was applied to CT of 68 regions in the Desikan-Killiany atlas derived from three large-scale datasets comprising a total of 41,075 neurotypical participants. PCA produced a spatially broad first principal component (PC1) that was reproducible across datasets. Then PC1 derived from healthy adult participants was compared to the pattern of CT differences associated with psychiatric and neurological disorders comprising a total of 14,886 cases and 20,962 controls from seven ENIGMA disease-related working groups, normative maturation and aging comprising a total of 17,697 scans from the ABCD Study® and the IMAGEN developmental study, and 17,075 participants from the ENIGMA Lifespan working group, as well as gene expression maps from the Allen Human Brain Atlas. Results revealed substantial spatial correspondences between PC1 and widespread lower CT observed in numerous psychiatric disorders. Moreover, the PC1 pattern was also correlated with the spatial pattern of normative maturation and aging. The transcriptional analysis identified a set of genes including KCNA2, KCNS1 and KCNS2 with expression patterns closely related to the spatial pattern of PC1. The gene category enrichment analysis indicated that the transcriptional correlations of PC1 were enriched to multiple gene ontology categories and were specifically over-represented starting at late childhood, coinciding with the onset of significant cortical maturation and emergence of psychopathology during the prepubertal-to-pubertal transition. Collectively, the present study reports a reproducible latent pattern of CT that captures interregional profiles of cortical changes in both normative brain maturation and a spectrum of psychiatric disorders. The pubertal timing of the expression of PC1-related genes implicates disrupted neurodevelopment in the pathogenesis of the spectrum of psychiatric diseases emerging during adolescence.

10-Year trajectories of depressive symptoms and subsequent brain health in middle-aged adults

Depressive symptoms differ in severity and stability over time. Trajectories depicting these changes, particularly those with high late-life depressive symptoms, have been associated with poor brain health at old age. To better understand these associations across the lifespan, we examined depressive symptoms trajectories in relation to brain health in middle age. We included 1676 participants from the ORACLE Study, all were expecting a child at baseline (mean age 32.8, 66.6% women). Depressive symptoms were assessed at baseline, 3 years and 10 years after baseline. Brain health (global brain volume, subcortical structures volume, white matter lesions, cerebral microbleeds, cortical thickness, cortical surface area) was assessed 15 years after baseline. Using k-means clustering, four depressive symptoms trajectories were identified: low, low increasing, decreasing, and high increasing symptoms. The high increasing trajectory was associated with smaller brain volume compared to low symptoms, not surviving multiple testing correction. The low increasing trajectory was associated with more cortical thickness in a small region encompassing the right lateral occipital cortex compared to low symptoms. These findings show that longitudinal depressive symptoms trajectories are only minimally associated with brain health in middle age, suggesting that associations may only emerge later in life.

Excitatory/inhibitory imbalance in autism: the role of glutamate and GABA gene-sets in symptoms and cortical brain structure

The excitatory/inhibitory (E/I) imbalance hypothesis posits that imbalance between excitatory (glutamatergic) and inhibitory (GABAergic) mechanisms underlies the behavioral characteristics of autism. However, how E/I imbalance arises and how it may differ across autism symptomatology and brain regions is not well understood. We used innovative analysis methods-combining competitive gene-set analysis and gene-expression profiles in relation to cortical thickness (CT) to investigate relationships between genetic variance, brain structure and autism symptomatology of participants from the AIMS-2-TRIALS LEAP cohort (autism = 359, male/female = 258/101; neurotypical control participants = 279, male/female = 178/101) aged 6-30 years. Using competitive gene-set analyses, we investigated whether aggregated genetic variation in glutamate and GABA gene-sets could be associated with behavioral measures of autism symptoms and brain structural variation. Further, using the same gene-sets, we corelated expression profiles throughout the cortex with differences in CT between autistic and neurotypical control participants, as well as in separate sensory subgroups. The glutamate gene-set was associated with all autism symptom severity scores on the Autism Diagnostic Observation Schedule-2 (ADOS-2) and the Autism Diagnostic Interview-Revised (ADI-R) within the autistic group. In adolescents and adults, brain regions with greater gene-expression of glutamate and GABA genes showed greater differences in CT between autistic and neurotypical control participants although in opposing directions. Additionally, the gene expression profiles were associated with CT profiles in separate sensory subgroups. Our results suggest complex relationships between E/I related genetics and autism symptom profiles as well as brain structure alterations, where there may be differential roles for glutamate and GABA.

Cortical thickness differences are associated with cellular component morphogenesis of astrocytes and excitatory neurons in nonsuicidal self-injuring youth

Nonsuicidal self-injury (NSSI) generally occurs in youth and probably progresses to suicide. An examination of cortical thickness differences (ΔCT) between NSSI individuals and controls is crucial to investigate potential neurobiological correlates. Notably, ΔCT are influenced by specific genetic factors, and a large proportion of cortical thinning is associated with the expression of genes that overlap in astrocytes and pyramidal cells. However, in NSSI youth, the mechanisms underlying the relations between the genetic and cell type-specific transcriptional signatures to ΔCT are unclear. Here, we studied the genetic association of ΔCT in NSSI youth by performing a partial least-squares regression (PLSR) analysis of gene expression data and 3D-T1 brain images of 45 NSSI youth and 75 controls. We extracted the top-10 Gene Ontology terms for the enrichment results of upregulated PLS component 1 genes related to ΔCT to conduct the cell-type classification and enrichment analysis. Enrichment of cell type-specific genes shows that cellular component morphogenesis of astrocytes and excitatory neurons accounts for the observed NSSI-specific ΔCT. We validated the main results in independent datasets to verify the robustness and specificity. We concluded that the brain ΔCT is associated with cellular component morphogenesis of astrocytes and excitatory neurons in NSSI youth.

Gene Expressions Preferentially Influence Cortical Thickness of Human Connectome Project Atlas Parcellated Regions in First-Episode Antipsychotic-Naïve Psychoses

Altered gene expressions may mechanistically link genetic factors with brain morphometric alterations. Existing gene expression studies have examined selected morphometric features using low-resolution atlases in medicated schizophrenia. We examined the relationship of gene expression with cortical thickness (CT), surface area (SA), and gray matter volume (GMV) of first-episode antipsychotic-naïve psychosis patients (FEAP = 85) and 81 controls, hypothesizing that gene expressions often associated with psychosis will differentially associate with different morphometric features. We explored such associations among schizophrenia and non-schizophrenia subgroups within FEAP group compared to controls. We mapped 360 Human Connectome Project atlas-based parcellations on brain MRI on to the publicly available brain gene expression data from the Allen Brain Institute collection. Significantly correlated genes were investigated using ingenuity pathway analysis to elucidate molecular pathways. CT but not SA or GMV correlated with expression of 1137 out of 15 633 genes examined controlling for age, sex, and average CT. Among these ≈19%, ≈39%, and 8% of genes were unique to FEAP, schizophrenia, and non-schizophrenia, respectively. Variants of 10 among these 1137 correlated genes previously showed genome-wide-association with schizophrenia. Molecular pathways associated with CT were axonal guidance and sphingosine pathways (common to FEAP and controls), selected inflammation pathways (unique to FEAP), synaptic modulation (unique to schizophrenia), and telomere extension (common to NSZ and healthy controls). We demonstrate that different sets of genes and molecular pathways may preferentially influence CT in different diagnostic groups. Genes with altered expressions correlating with CT and associated pathways may be targets for pathophysiological investigations and novel treatment designs.

Quantitative MRI maps of human neocortex explored using cell type-specific gene expression analysis

Quantitative magnetic resonance imaging (qMRI) allows extraction of reproducible and robust parameter maps. However, the connection to underlying biological substrates remains murky, especially in the complex, densely packed cortex. We investigated associations in human neocortex between qMRI parameters and neocortical cell types by comparing the spatial distribution of the qMRI parameters longitudinal relaxation rate (${R_{1}}$), effective transverse relaxation rate (${R_{2}}^{\ast }$), and magnetization transfer saturation (MTsat) to gene expression from the Allen Human Brain Atlas, then combining this with lists of genes enriched in specific cell types found in the human brain. As qMRI parameters are magnetic field strength-dependent, the analysis was performed on MRI data at 3T and 7T. All qMRI parameters significantly covaried with genes enriched in GABA- and glutamatergic neurons, i.e. they were associated with cytoarchitecture. The qMRI parameters also significantly covaried with the distribution of genes enriched in astrocytes (${R_{2}}^{\ast }$ at 3T, ${R_{1}}$ at 7T), endothelial cells (${R_{1}}$ and MTsat at 3T), microglia (${R_{1}}$ and MTsat at 3T, ${R_{1}}$ at 7T), and oligodendrocytes and oligodendrocyte precursor cells (${R_{1}}$ at 7T). These results advance the potential use of qMRI parameters as biomarkers for specific cell types.

Regional gene expression patterns are associated with task-specific brain activation during reward and emotion processing measured with functional MRI

The exploration of the spatial relationship between gene expression profiles and task-evoked response patterns known to be altered in neuropsychiatric disorders, for example depression, can guide the development of more targeted therapies. Here, we estimated the correlation between human transcriptome data and two different brain activation maps measured with functional magnetic resonance imaging (fMRI) in healthy subjects. Whole-brain activation patterns evoked during an emotional face recognition task were associated with topological mRNA expression of genes involved in cellular transport. In contrast, fMRI activation patterns related to the acceptance of monetary rewards were associated with genes implicated in cellular localization processes, metabolism, translation, and synapse regulation. An overlap of these genes with risk genes from major depressive disorder genome-wide association studies revealed the involvement of the master regulators TCF4 and PAX6 in emotion and reward processing. Overall, the identification of stable relationships between spatial gene expression profiles and fMRI data may reshape the prospects for imaging transcriptomics studies.

Neurofilament light-associated connectivity in young-adult Huntington's disease is related to neuronal genes

Upregulation of functional network connectivity in the presence of structural degeneration is seen in the premanifest stages of Huntington's disease (preHD) 10-15 years from clinical diagnosis. However, whether widespread network connectivity changes are seen in gene carriers much further from onset has yet to be explored. We characterized functional network connectivity throughout the brain and related it to a measure of disease pathology burden (CSF neurofilament light, NfL) and measures of structural connectivity in asymptomatic gene carriers, on average 24 years from onset. We related these measurements to estimates of cortical and subcortical gene expression. We found no overall differences in functional (or structural) connectivity anywhere in the brain comparing control and preHD participants. However, increased functional connectivity, particularly between posterior cortical areas, correlated with increasing CSF NfL level in preHD participants. Using the Allen Human Brain Atlas and expression-weighted cell-type enrichment analysis, we demonstrated that this functional connectivity upregulation occurred in cortical regions associated with regional expression of genes specific to neuronal cells. This relationship was validated using single-nucleus RNAseq data from post-mortem Huntington's disease and control brains showing enrichment of neuronal-specific genes that are differentially expressed in Huntington's disease. Functional brain networks in asymptomatic preHD gene carriers very far from disease onset show evidence of upregulated connectivity correlating with increased disease burden. These changes occur among brain areas that show regional expression of genes specific to neuronal GABAergic and glutamatergic cells.

Dynamic Structural Brain Changes in Anorexia Nervosa: A Replication Study, Mega-analysis, and Virtual Histology Approach

OBJECTIVE

Several, but not all, previous studies of brain structure in anorexia nervosa (AN) have reported reductions in gray matter volume and cortical thickness (CT) in acutely underweight patients, which seem to reverse upon weight gain. The biological mechanisms underlying these dynamic alterations remain unclear.

METHOD

In this structural magnetic resonance imaging study, we first replicated and extended previous results in (1) a larger independent sample of 75 acutely underweight adolescent and young adult female patients with AN (acAN; n = 54 rescanned longitudinally after partial weight restoration), 34 weight-recovered individuals with a history of AN (recAN), and 139 healthy controls (HC); and 2) a greater combined sample compiled of both our previous samples and the present replication sample (120 acAN [90 rescanned longitudinally], 68 recAN, and 207 HC). Next, we applied a "virtual histology" approach to the combined data, investigating relations between interregional profiles of differences in CT and profiles of cell-specific gene expression. Finally, we used the ENIGMA toolbox to relate aforementioned CT profiles to normative structural and functional connectomics.

RESULTS

We confirmed sizeable and widespread reductions of CT as well as volumes (and, to a lesser extent, surface area) in acAN and rapid increases related to partial weight restoration. No differences were detected between either short- or long-term weight-recovered patients and HC. The virtual histology analysis identified associations between gene expression profiles of S1 pyramidal cells and oligodendrocytes and brain regions with more marked differences in CT, whereas the remaining regions were those with a greater expression of genes specific to CA1 pyramidal, astrocytes, microglia, and ependymal cells. Furthermore, the most affected regions were also more functionally and structurally connected.

CONCLUSION

The overall data pattern deviates from findings in other psychiatric disorders. Both virtual histology and connectomics analyses indicated that brain regions most affected in AN are also the most energetically demanding.

Obesity-Related Neuroinflammation: Magnetic Resonance and Microscopy Imaging of the Brain

Obesity is a major risk factor of Alzheimer’s disease and related dementias. The principal feature of dementia is a loss of neurons and brain atrophy. The mechanistic links between obesity and the neurodegenerative processes of dementias are not fully understood, but recent research suggests that obesity-related systemic inflammation and subsequent neuroinflammation may be involved. Adipose tissues release multiple proinflammatory molecules (fatty acids and cytokines) that impact blood and vessel cells, inducing low-grade systemic inflammation that can transition to tissues, including the brain. Inflammation in the brain—neuroinflammation—is one of key elements of the pathobiology of neurodegenerative disorders; it is characterized by the activation of microglia, the resident immune cells in the brain, and by the structural and functional changes of other cells forming the brain parenchyma, including neurons. Such cellular changes have been shown in animal models with direct methods, such as confocal microscopy. In humans, cellular changes are less tangible, as only indirect methods such as magnetic resonance (MR) imaging are usually used. In these studies, obesity and low-grade systemic inflammation have been associated with lower volumes of the cerebral gray matter, cortex, and hippocampus, as well as altered tissue MR properties (suggesting microstructural variations in cellular and molecular composition). How these structural variations in the human brain observed using MR imaging relate to the cellular variations in the animal brain seen with microscopy is not well understood. This review describes the current understanding of neuroinflammation in the context of obesity-induced systemic inflammation, and it highlights need for the bridge between animal microscopy and human MR imaging studies.

Genetic mechanisms underlying brain functional homotopy: a combined transcriptome and resting-state functional MRI study

Functional homotopy, the high degree of spontaneous activity synchrony and functional coactivation between geometrically corresponding interhemispheric regions, is a fundamental characteristic of the intrinsic functional architecture of the brain. However, little is known about the genetic mechanisms underlying functional homotopy. Resting-state functional magnetic resonance imaging data from a discovery dataset (656 healthy subjects) and 2 independent cross-race, cross-scanner validation datasets (103 and 329 healthy subjects) were used to calculate voxel-mirrored homotopic connectivity (VMHC) indexing brain functional homotopy. In combination with the Allen Human Brain Atlas, transcriptome-neuroimaging spatial correlation analysis was conducted to identify genes linked to VMHC. We found 1,001 genes whose expression measures were spatially associated with VMHC. Functional enrichment analyses demonstrated that these VMHC-related genes were enriched for biological functions including protein kinase activity, ion channel regulation, and synaptic function as well as many neuropsychiatric disorders. Concurrently, specific expression analyses showed that these genes were specifically expressed in the brain tissue, in neurons and immune cells, and during nearly all developmental periods. In addition, the VMHC-associated genes were linked to multiple behavioral domains, including vision, execution, and attention. Our findings suggest that interhemispheric communication and coordination involve a complex interaction of polygenes with a rich range of functional features.

Auditory driven gamma synchrony is associated with cortical thickness in widespread cortical areas

OBJECTIVE

Gamma synchrony is a fundamental functional property of the cerebral cortex, impaired in multiple neuropsychiatric conditions (i.e. schizophrenia, Alzheimer's disease, stroke etc.). Auditory stimulation in the gamma range allows to drive gamma synchrony of the entire cortical mantle and to estimate the efficiency of the mechanisms sustaining it. As gamma synchrony depends strongly on the interplay between parvalbumin-positive interneurons and pyramidal neurons, we hypothesize an association between cortical thickness and gamma synchrony. To test this hypothesis, we employed a combined magnetoencephalography (MEG) - Magnetic Resonance Imaging (MRI) study.

METHODS

Cortical thickness was estimated from anatomical MRI scans. MEG measurements related to exposure of 40 Hz amplitude modulated tones were projected onto the cortical surface. Two measures of cortical synchrony were considered: (a) inter-trial phase consistency at 40 Hz, providing a vertex-wise estimation of gamma synchronization, and (b) phase-locking values between primary auditory cortices and whole cortical mantle, providing a measure of long-range cortical synchrony. A correlation between cortical thickness and synchronization measures was then calculated for 72 MRI-MEG scans.

RESULTS

Both inter-trial phase consistency and phase locking values showed a significant positive correlation with cortical thickness. For inter-trial phase consistency, clusters of strong associations were found in the temporal and frontal lobes, especially in the bilateral auditory and pre-motor cortices. Higher phase-locking values corresponded to higher cortical thickness in the frontal, temporal, occipital and parietal lobes.

DISCUSSION AND CONCLUSIONS

In healthy subjects, a thicker cortex corresponds to higher gamma synchrony and connectivity in the primary auditory cortex and beyond, likely reflecting underlying cell density involved in gamma circuitries. This result hints towards an involvement of gamma synchrony together with underlying brain structure in brain areas for higher order cognitive functions. This study contributes to the understanding of inherent cortical functional and structural brain properties, which might in turn constitute the basis for the definition of useful biomarkers in patients showing aberrant gamma synchronization.

Transcriptomic and cellular decoding of functional brain connectivity changes reveal regional brain vulnerability to pro- and anti-inflammatory therapies

BACKGROUND

Systemic inflammation induces acute changes in mood, motivation and cognition that closely resemble those observed in depressed individuals. However, the mechanistic pathways linking peripheral inflammation to depression-like psychopathology via intermediate effects on brain function remain incompletely understood.

METHODS

We combined data from 30 patients initiating interferon-α treatment for Hepatitis-C and 20 anti-tumour necrosis factor (TNF) therapy for inflammatory arthritis and used resting-state functional magnetic resonance imaging to investigate acute effects of each treatment on regional global brain connectivity (GBC). We leveraged transcriptomic data from the Allen Human Brain Atlas to uncover potential biological and cellular pathways underpinning regional vulnerability to GBC changes induced by each treatment.

RESULTS

Interferon-α and anti-TNF therapies both produced differential small-to-medium sized decreases in regional GBC. However, these were observed within distinct brain regions and the regional patterns of GBC changes induced by each treatment did not correlate suggesting independent underlying processes. Further, the spatial distribution of these differential GBC decreases could be captured by multivariate patterns of constitutive regional expression of genes respectively related to: i) neuroinflammation and glial cells; and ii) glutamatergic neurotransmission and neurons. The extent to which each participant expressed patterns of GBC changes aligning with these patterns of transcriptomic vulnerability also correlated with both acute treatment-induced changes in interleukin-6 (IL-6) and, for Interferon-α, longer-term treatment-associated changes in depressive symptoms.

CONCLUSIONS

Together, we present two transcriptomic models separately linking regional vulnerability to the acute effects of interferon-α and anti-TNF treatments on brain function to glial neuroinflammation and glutamatergic neurotransmission. These findings generate hypotheses about two potential brain mechanisms through which bidirectional changes in peripheral inflammation may contribute to the development/resolution of psychopathology.

Cell type-specific manifestations of cortical thickness heterogeneity in schizophrenia

Brain morphology differs markedly between individuals with schizophrenia, but the cellular and genetic basis of this heterogeneity is poorly understood. Here, we sought to determine whether cortical thickness (CTh) heterogeneity in schizophrenia relates to interregional variation in distinct neural cell types, as inferred from established gene expression data and person-specific genomic variation. This study comprised 1849 participants in total, including a discovery (140 cases and 1267 controls) and a validation cohort (335 cases and 185 controls). To characterize CTh heterogeneity, normative ranges were established for 34 cortical regions and the extent of deviation from these ranges was measured for each individual with schizophrenia. CTh deviations were explained by interregional gene expression levels of five out of seven neural cell types examined: (1) astrocytes; (2) endothelial cells; (3) oligodendrocyte progenitor cells (OPCs); (4) excitatory neurons; and (5) inhibitory neurons. Regional alignment between CTh alterations with cell type transcriptional maps distinguished broad patient subtypes, which were validated against genomic data drawn from the same individuals. In a predominantly neuronal/endothelial subtype (22% of patients), CTh deviations covaried with polygenic risk for schizophrenia (sczPRS) calculated specifically from genes marking neuronal and endothelial cells (r = -0.40, p = 0.010). Whereas, in a predominantly glia/OPC subtype (43% of patients), CTh deviations covaried with sczPRS calculated from glia and OPC-linked genes (r = -0.30, p = 0.028). This multi-scale analysis of genomic, transcriptomic, and brain phenotypic data may indicate that CTh heterogeneity in schizophrenia relates to inter-individual variation in cell-type specific functions. Decomposing heterogeneity in relation to cortical cell types enables prioritization of schizophrenia subsets for future disease modeling efforts.

What we learn about bipolar disorder from large-scale neuroimaging: Findings and future directions from the ENIGMA Bipolar Disorder Working Group

MRI-derived brain measures offer a link between genes, the environment and behavior and have been widely studied in bipolar disorder (BD). However, many neuroimaging studies of BD have been underpowered, leading to varied results and uncertainty regarding effects. The Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Bipolar Disorder Working Group was formed in 2012 to empower discoveries, generate consensus findings and inform future hypothesis-driven studies of BD. Through this effort, over 150 researchers from 20 countries and 55 institutions pool data and resources to produce the largest neuroimaging studies of BD ever conducted. The ENIGMA Bipolar Disorder Working Group applies standardized processing and analysis techniques to empower large-scale meta- and mega-analyses of multimodal brain MRI and improve the replicability of studies relating brain variation to clinical and genetic data. Initial BD Working Group studies reveal widespread patterns of lower cortical thickness, subcortical volume and disrupted white matter integrity associated with BD. Findings also include mapping brain alterations of common medications like lithium, symptom patterns and clinical risk profiles and have provided further insights into the pathophysiological mechanisms of BD. Here we discuss key findings from the BD working group, its ongoing projects and future directions for large-scale, collaborative studies of mental illness.

Effects of copy number variations on brain structure and risk for psychiatric illness: Large-scale studies from the ENIGMA working groups on CNVs

The Enhancing NeuroImaging Genetics through Meta-Analysis copy number variant (ENIGMA-CNV) and 22q11.2 Deletion Syndrome Working Groups (22q-ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA-CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q-ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest-ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi-site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene-dosage effects on distinct brain regions also emerged, providing further insight into genotype-phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This "genotype-first" approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

An overview of the first 5 years of the ENIGMA obsessive-compulsive disorder working group: The power of worldwide collaboration

Neuroimaging has played an important part in advancing our understanding of the neurobiology of obsessive-compulsive disorder (OCD). At the same time, neuroimaging studies of OCD have had notable limitations, including reliance on relatively small samples. International collaborative efforts to increase statistical power by combining samples from across sites have been bolstered by the ENIGMA consortium; this provides specific technical expertise for conducting multi-site analyses, as well as access to a collaborative community of neuroimaging scientists. In this article, we outline the background to, development of, and initial findings from ENIGMA's OCD working group, which currently consists of 47 samples from 34 institutes in 15 countries on 5 continents, with a total sample of 2,323 OCD patients and 2,325 healthy controls. Initial work has focused on studies of cortical thickness and subcortical volumes, structural connectivity, and brain lateralization in children, adolescents and adults with OCD, also including the study on the commonalities and distinctions across different neurodevelopment disorders. Additional work is ongoing, employing machine learning techniques. Findings to date have contributed to the development of neurobiological models of OCD, have provided an important model of global scientific collaboration, and have had a number of clinical implications. Importantly, our work has shed new light on questions about whether structural and functional alterations found in OCD reflect neurodevelopmental changes, effects of the disease process, or medication impacts. We conclude with a summary of ongoing work by ENIGMA-OCD, and a consideration of future directions for neuroimaging research on OCD within and beyond ENIGMA.

Linking Central Gene Expression Patterns and Mental States Using Transcriptomics and Large-Scale Meta-Analysis of fMRI Data: A Tutorial and Example Using the Oxytocin Signaling Pathway

The measurement of gene expression levels in the human brain can help accelerate our understanding of complex mental states and psychiatric illnesses. Mental states are typically associated with whole-brain networks; however, gene expression levels from postmortem brain samples have traditionally been measured in a limited number of brain regions due to resource limitations. The recent availability of whole-brain gene expression data from the Allen Human Brain Atlas (AHBA) provides the opportunity to generate gene expression patterns for over 20,000 genes. By linking these expression patterns with brain activity patterns that are associated with specific mental states, researchers can better understand which genes may support given mental states, via forward inference. Conversely, reverse inference can also be used to determine which mental state activation patterns are most strongly associated with a given gene expression map. This chapter provides a step-by-step guide on how to use the AHBA in conjunction with the NeuroSynth fMRI meta-analysis tool to identify the mental state correlates of specific gene expression patterns, using genes from oxytocin signaling pathway as an example. We also demonstrate how to perform an out-of-sample validation and assess the specificity of results for genes of interest.

Cortical thickness abnormalities in patients with post-traumatic stress disorder: A vertex-based meta-analysis

Neuroimaging studies report altered cortical thickness in patients with post-traumatic stress disorder (PTSD), but the results are inconsistent. Using anisotropic effect-size seed-based d mapping (AES-SDM) software with its recently-developed meta-analytic thickness mask, we conducted a meta-analysis of published studies which used whole-brain surface-based morphometry, in order to define consistent cortical thickness alterations in PTSD patients. Eleven studies with 438 patients and 396 controls were included. Compared with all controls, patients with PTSD showed increased cortical thickness in right superior temporal gyrus, and in left and right superior frontal gyrus; the former survived in subgroup analysis of adult patients, and in subgroup comparison with only non-PTSD trauma-exposed controls, the latter in subgroup comparison with only non-trauma-exposed healthy controls. Cortical thickness in right superior frontal gyrus was positively associated with percentage of female patients, and cortical thickness in left superior frontal gyrus was positively associated with symptom severity measured by the clinician-administered PTSD scale. These robust results may help to elucidate the pathophysiology of PTSD.

Consortium neuroscience of attention deficit/hyperactivity disorder and autism spectrum disorder: The ENIGMA adventure

Neuroimaging has been extensively used to study brain structure and function in individuals with attention deficit/hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) over the past decades. Two of the main shortcomings of the neuroimaging literature of these disorders are the small sample sizes employed and the heterogeneity of methods used. In 2013 and 2014, the ENIGMA-ADHD and ENIGMA-ASD working groups were respectively, founded with a common goal to address these limitations. Here, we provide a narrative review of the thus far completed and still ongoing projects of these working groups. Due to an implicitly hierarchical psychiatric diagnostic classification system, the fields of ADHD and ASD have developed largely in isolation, despite the considerable overlap in the occurrence of the disorders. The collaboration between the ENIGMA-ADHD and -ASD working groups seeks to bring the neuroimaging efforts of the two disorders closer together. The outcomes of case-control studies of subcortical and cortical structures showed that subcortical volumes are similarly affected in ASD and ADHD, albeit with small effect sizes. Cortical analyses identified unique differences in each disorder, but also considerable overlap between the two, specifically in cortical thickness. Ongoing work is examining alternative research questions, such as brain laterality, prediction of case-control status, and anatomical heterogeneity. In brief, great strides have been made toward fulfilling the aims of the ENIGMA collaborations, while new ideas and follow-up analyses continue that include more imaging modalities (diffusion MRI and resting-state functional MRI), collaborations with other large databases, and samples with dual diagnoses.

Genetic variation, brain, and intelligence differences

Individual differences in human intelligence, as assessed using cognitive test scores, have a well-replicated, hierarchical phenotypic covariance structure. They are substantially stable across the life course, and are predictive of educational, social, and health outcomes. From this solid phenotypic foundation and importance for life, comes an interest in the environmental, social, and genetic aetiologies of intelligence, and in the foundations of intelligence differences in brain structure and functioning. Here, we summarise and critique the last 10 years or so of molecular genetic (DNA-based) research on intelligence, including the discovery of genetic loci associated with intelligence, DNA-based heritability, and intelligence's genetic correlations with other traits. We summarise new brain imaging-intelligence findings, including whole-brain associations and grey and white matter associations. We summarise regional brain imaging associations with intelligence and interpret these with respect to theoretical accounts. We address research that combines genetics and brain imaging in studying intelligence differences. There are new, though modest, associations in all these areas, and mechanistic accounts are lacking. We attempt to identify growing points that might contribute toward a more integrated 'systems biology' account of some of the between-individual differences in intelligence.

Imaging transcriptomics: Convergent cellular, transcriptomic, and molecular neuroimaging signatures in the healthy adult human brain

The integration of transcriptomic and neuroimaging data, "imaging transcriptomics," has recently emerged to generate hypotheses about potential biological pathways underlying regional variability in neuroimaging features. However, the validity of this approach is yet to be examined in depth. Here, we sought to bridge this gap by performing transcriptomic decoding of the regional distribution of well-known molecular markers spanning different elements of the biology of the healthy human brain. Imaging transcriptomics identifies biological and cell pathways that are consistent with the known biology of a wide range of molecular neuroimaging markers. The extent to which it can capture patterns of gene expression that align well with elements of the biology of the neuroinflammatory axis, at least in healthy controls without a proinflammatory challenge, is inconclusive. Imaging transcriptomics might constitute an interesting approach to improve our understanding of the biological pathways underlying regional variability in a wide range of neuroimaging phenotypes.

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.