The impact of psychosocial adversity on brain and behaviour: an overview of existing knowledge and directions for future research.

Environmental experiences play a critical role in shaping the structure and function of the brain. Its plasticity in response to different external stimuli has been the focus of research efforts for decades. In this review, we explore the effects of adversity on brain's structure and function and its implications for brain development, adaptation, and the emergence of mental health disorders. We are focusing on adverse events that emerge from the immediate surroundings of an individual, i.e., microenvironment. They include childhood maltreatment, peer victimisation, social isolation, affective loss, domestic conflict, and poverty. We also take into consideration exposure to environmental toxins. Converging evidence suggests that different types of adversity may share common underlying mechanisms while also exhibiting unique pathways. However, they are often studied in isolation, limiting our understanding of their combined effects and the interconnected nature of their impact. The integration of large, deep-phenotyping datasets and collaborative efforts can provide sufficient power to analyse high dimensional environmental profiles and advance the systematic mapping of neuronal mechanisms. This review provides a background for future research, highlighting the importance of understanding the cumulative impact of various adversities, through data-driven approaches and integrative multimodal analysis techniques.

Estimating cortical thickness trajectories in children across different scanners using transfer learning from normative models.

This work illustrates the use of normative models in a longitudinal neuroimaging study of children aged 6-17 years and demonstrates how such models can be used to make meaningful comparisons in longitudinal studies, even when individuals are scanned with different scanners across successive study waves. More specifically, we first estimated a large-scale reference normative model using Hierarchical Bayesian Regression from N = 42,993 individuals across the lifespan and from dozens of sites. We then transfer these models to a longitudinal developmental cohort (N = 6285) with three measurement waves acquired on two different scanners that were unseen during estimation of the reference models. We show that the use of normative models provides individual deviation scores that are independent of scanner effects and efficiently accommodate inter-site variations. Moreover, we provide empirical evidence to guide the optimization of sample size for the transfer of prior knowledge about the distribution of regional cortical thicknesses. We show that a transfer set containing as few as 25 samples per site can lead to good performance metrics on the test set. Finally, we demonstrate the clinical utility of this approach by showing that deviation scores obtained from the transferred normative models are able to detect and chart morphological heterogeneity in individuals born preterm.

Personalizing progressive changes to brain structure in Alzheimer's disease using normative modeling.

INTRODUCTION: Neuroanatomical normative modeling captures individual variability in Alzheimer's disease (AD). Here we used normative modeling to track individuals' disease progression in people with mild cognitive impairment (MCI) and patients with AD. METHODS: Cortical and subcortical normative models were generated using healthy controls (n ≈ 58k). These models were used to calculate regional z scores in 3233 T1-weighted magnetic resonance imaging time-series scans from 1181 participants. Regions with z scores < -1.96 were classified as outliers mapped on the brain and summarized by total outlier count (tOC). RESULTS: tOC increased in AD and in people with MCI who converted to AD and also correlated with multiple non-imaging markers. Moreover, a higher annual rate of change in tOC increased the risk of progression from MCI to AD. Brain outlier maps identified the hippocampus as having the highest rate of change. DISCUSSION: Individual patients' atrophy rates can be tracked by using regional outlier maps and tOC. HIGHLIGHTS: Neuroanatomical normative modeling was applied to serial Alzheimer's disease (AD) magnetic resonance imaging (MRI) data for the first time. Deviation from the norm (outliers) of cortical thickness or brain volume was computed in 3233 scans. The number of brain-structure outliers increased over time in people with AD. Patterns of change in outliers varied markedly between individual patients with AD. People with mild cognitive impairment whose outliers increased over time had a higher risk of progression from AD.

Age-related brain deviations and aggression.

BACKGROUND: Disruptive behavior disorders (DBD) are heterogeneous at the clinical and the biological level. Therefore, the aims were to dissect the heterogeneous neurodevelopmental deviations of the affective brain circuitry and provide an integration of these differences across modalities. METHODS: We combined two novel approaches. First, normative modeling to map deviations from the typical age-related pattern at the level of the individual of (i) activity during emotion matching and (ii) of anatomical images derived from DBD cases (n = 77) and controls (n = 52) aged 8-18 years from the EU-funded Aggressotype and MATRICS consortia. Second, linked independent component analysis to integrate subject-specific deviations from both modalities. RESULTS: While cases exhibited on average a higher activity than would be expected for their age during face processing in regions such as the amygdala when compared to controls these positive deviations were widespread at the individual level. A multimodal integration of all functional and anatomical deviations explained 23% of the variance in the clinical DBD phenotype. Most notably, the top marker, encompassing the default mode network (DMN) and subcortical regions such as the amygdala and the striatum, was related to aggression across the whole sample. CONCLUSIONS: Overall increased age-related deviations in the amygdala in DBD suggest a maturational delay, which has to be further validated in future studies. Further, the integration of individual deviation patterns from multiple imaging modalities allowed to dissect some of the heterogeneity of DBD and identified the DMN, the striatum and the amygdala as neural signatures that were associated with aggression.

Patterns of sociocognitive stratification and perinatal risk in the child brain.

The expanding behavioral repertoire of the developing brain during childhood and adolescence is shaped by complex brain-environment interactions and flavored by unique life experiences. The transition into young adulthood offers opportunities for adaptation and growth but also increased susceptibility to environmental perturbations, such as the characteristics of social relationships, family environment, quality of schools and activities, financial security, urbanization and pollution, drugs, cultural practices, and values, that all act in concert with our genetic architecture and biology. Our multivariate brain-behavior mapping in 7,577 children aged 9 to 11 y across 585 brain imaging phenotypes and 617 cognitive, behavioral, psychosocial, and socioeconomic measures revealed three population modes of brain covariation, which were robust as assessed by cross-validation and permutation testing, taking into account siblings and twins, identified using genetic data. The first mode revealed traces of perinatal complications, including preterm and twin birth, eclampsia and toxemia, shorter period of breastfeeding, and lower cognitive scores, with higher cortical thickness and lower cortical areas and volumes. The second mode reflected a pattern of sociocognitive stratification, linking lower cognitive ability and socioeconomic status to lower cortical thickness, area, and volumes. The third mode captured a pattern related to urbanicity, with particulate matter pollution (PM25) inversely related to home value, walkability, and population density, associated with diffusion properties of white matter tracts. These results underscore the importance of a multidimensional and interdisciplinary understanding, integrating social, psychological, and biological sciences, to map the constituents of healthy development and to identify factors that may precede maladjustment and mental illness.

Accommodating site variation in neuroimaging data using normative and hierarchical Bayesian models.

The potential of normative modeling to make individualized predictions from neuroimaging data has enabled inferences that go beyond the case-control approach. However, site effects are often confounded with variables of interest in a complex manner and can bias estimates of normative models, which has impeded the application of normative models to large multi-site neuroimaging data sets. In this study, we suggest accommodating for these site effects by including them as random effects in a hierarchical Bayesian model. We compared the performance of a linear and a non-linear hierarchical Bayesian model in modeling the effect of age on cortical thickness. We used data of 570 healthy individuals from the ABIDE (autism brain imaging data exchange) data set in our experiments. In addition, we used data from individuals with autism to test whether our models are able to retain clinically useful information while removing site effects. We compared the proposed single stage hierarchical Bayesian method to several harmonization techniques commonly used to deal with additive and multiplicative site effects using a two stage regression, including regressing out site and harmonizing for site with ComBat, both with and without explicitly preserving variance caused by age and sex as biological variation of interest, and with a non-linear version of ComBat. In addition, we made predictions from raw data, in which site has not been accommodated for. The proposed hierarchical Bayesian method showed the best predictive performance according to multiple metrics. Beyond that, the resulting z-scores showed little to no residual site effects, yet still retained clinically useful information. In contrast, performance was particularly poor for the regression model and the ComBat model in which age and sex were not explicitly modeled. In all two stage harmonization models, predictions were poorly scaled, suffering from a loss of more than 90% of the original variance. Our results show the value of hierarchical Bayesian regression methods for accommodating site variation in neuroimaging data, which provides an alternative to harmonization techniques. While the approach we propose may have broad utility, our approach is particularly well suited to normative modeling where the primary interest is in accurate modeling of inter-subject variation and statistical quantification of deviations from a reference model.

Deep neural networks learn general and clinically relevant representations of the ageing brain.

The discrepancy between chronological age and the apparent age of the brain based on neuroimaging data - the brain age delta - has emerged as a reliable marker of brain health. With an increasing wealth of data, approaches to tackle heterogeneity in data acquisition are vital. To this end, we compiled raw structural magnetic resonance images into one of the largest and most diverse datasets assembled (n=53542), and trained convolutional neural networks (CNNs) to predict age. We achieved state-of-the-art performance on unseen data from unknown scanners (n=2553), and showed that higher brain age delta is associated with diabetes, alcohol intake and smoking. Using transfer learning, the intermediate representations learned by our model complemented and partly outperformed brain age delta in predicting common brain disorders. Our work shows we can achieve generalizable and biologically plausible brain age predictions using CNNs trained on heterogeneous datasets, and transfer them to clinical use cases.

Rapid processing and quantitative evaluation of structural brain scans for adaptive multimodal imaging.

Current neuroimaging acquisition and processing approaches tend to be optimised for quality rather than speed. However, rapid acquisition and processing of neuroimaging data can lead to novel neuroimaging paradigms, such as adaptive acquisition, where rapidly processed data is used to inform subsequent image acquisition steps. Here we first evaluate the impact of several processing steps on the processing time and quality of registration of manually labelled T1 -weighted MRI scans. Subsequently, we apply the selected rapid processing pipeline both to rapidly acquired multicontrast EPImix scans of 95 participants (which include T1 -FLAIR, T2 , T2 *, T2 -FLAIR, DWI and ADC contrasts, acquired in ~1 min), as well as to slower, more standard single-contrast T1 -weighted scans of a subset of 66 participants. We quantify the correspondence between EPImix T1 -FLAIR and single-contrast T1 -weighted scans, using correlations between voxels and regions of interest across participants, measures of within- and between-participant identifiability as well as regional structural covariance networks. Furthermore, we explore the use of EPImix for the rapid construction of morphometric similarity networks. Finally, we quantify the reliability of EPImix-derived data using test-retest scans of 10 participants. Our results demonstrate that quantitative information can be derived from a neuroimaging scan acquired and processed within minutes, which could further be used to implement adaptive multimodal imaging and tailor neuroimaging examinations to individual patients.

Beyond the average patient: how neuroimaging models can address heterogeneity in dementia.

Dementia is a highly heterogeneous condition, with pronounced individual differences in age of onset, clinical presentation, progression rates and neuropathological hallmarks, even within a specific diagnostic group. However, the most common statistical designs used in dementia research studies and clinical trials overlook this heterogeneity, instead relying on comparisons of group average differences (e.g. patient versus control or treatment versus placebo), implicitly assuming within-group homogeneity. This one-size-fits-all approach potentially limits our understanding of dementia aetiology, hindering the identification of effective treatments. Neuroimaging has enabled the characterization of the average neuroanatomical substrates of dementias; however, the increasing availability of large open neuroimaging datasets provides the opportunity to examine patterns of neuroanatomical variability in individual patients. In this update, we outline the causes and consequences of heterogeneity in dementia and discuss recent research that aims to tackle heterogeneity directly, rather than assuming that dementia affects everyone in the same way. We introduce spatial normative modelling as an emerging data-driven technique, which can be applied to dementia data to model neuroanatomical variation, capturing individualized neurobiological 'fingerprints'. Such methods have the potential to detect clinically relevant subtypes, track an individual's disease progression or evaluate treatment responses, with the goal of moving towards precision medicine for dementia.

Phenotyping the Preterm Brain: Characterizing Individual Deviations From Normative Volumetric Development in Two Large Infant Cohorts.

The diverse cerebral consequences of preterm birth create significant challenges for understanding pathogenesis or predicting later outcome. Instead of focusing on describing effects common to the group, comparing individual infants against robust normative data offers a powerful alternative to study brain maturation. Here we used Gaussian process regression to create normative curves characterizing brain volumetric development in 274 term-born infants, modeling for age at scan and sex. We then compared 89 preterm infants scanned at term-equivalent age with these normative charts, relating individual deviations from typical volumetric development to perinatal risk factors and later neurocognitive scores. To test generalizability, we used a second independent dataset comprising of 253 preterm infants scanned using different acquisition parameters and scanner. We describe rapid, nonuniform brain growth during the neonatal period. In both preterm cohorts, cerebral atypicalities were widespread, often multiple, and varied highly between individuals. Deviations from normative development were associated with respiratory support, nutrition, birth weight, and later neurocognition, demonstrating their clinical relevance. Group-level understanding of the preterm brain disguises a large degree of individual differences. We provide a method and normative dataset that offer a more precise characterization of the cerebral consequences of preterm birth by profiling the individual neonatal brain.

Sensory Modality-Independent Activation of the Brain Network for Language.

The meaning of a sentence can be understood, whether presented in written or spoken form. Therefore, it is highly probable that brain processes supporting language comprehension are at least partly independent of sensory modality. To identify where and when in the brain language processing is independent of sensory modality, we directly compared neuromagnetic brain signals of 200 human subjects (102 males) either reading or listening to sentences. We used multiset canonical correlation analysis to align individual subject data in a way that boosts those aspects of the signal that are common to all, allowing us to capture word-by-word signal variations, consistent across subjects and at a fine temporal scale. Quantifying this consistency in activation across both reading and listening tasks revealed a mostly left-hemispheric cortical network. Areas showing consistent activity patterns included not only areas previously implicated in higher-level language processing, such as left prefrontal, superior and middle temporal areas, and anterior temporal lobe, but also parts of the control network as well as subcentral and more posterior temporal-parietal areas. Activity in this supramodal sentence-processing network starts in temporal areas and rapidly spreads to the other regions involved. The findings indicate not only the involvement of a large network of brain areas in supramodal language processing but also that the linguistic information contained in the unfolding sentences modulates brain activity in a word-specific manner across subjects.SIGNIFICANCE STATEMENT The brain can extract meaning from written and spoken messages alike. This requires activity of both brain circuits capable of processing sensory modality-specific aspects of the input signals as well as coordinated brain activity to extract modality-independent meaning from the input. Using traditional methods, it is difficult to disentangle modality-specific activation from modality-independent activation. In this work, we developed and applied a multivariate methodology that allows for a direct quantification of sensory modality-independent brain activity, revealing fast activation of a wide network of brain areas, both including and extending beyond the core network for language.

Warped Bayesian linear regression for normative modelling of big data.

Normative modelling is becoming more popular in neuroimaging due to its ability to make predictions of deviation from a normal trajectory at the level of individual participants. It allows the user to model the distribution of several neuroimaging modalities, giving an estimation for the mean and centiles of variation. With the increase in the availability of big data in neuroimaging, there is a need to scale normative modelling to big data sets. However, the scaling of normative models has come with several challenges. So far, most normative modelling approaches used Gaussian process regression, and although suitable for smaller datasets (up to a few thousand participants) it does not scale well to the large cohorts currently available and being acquired. Furthermore, most neuroimaging modelling methods that are available assume the predictive distribution to be Gaussian in shape. However, deviations from Gaussianity can be frequently found, which may lead to incorrect inferences, particularly in the outer centiles of the distribution. In normative modelling, we use the centiles to give an estimation of the deviation of a particular participant from the 'normal' trend. Therefore, especially in normative modelling, the correct estimation of the outer centiles is of utmost importance, which is also where data are sparsest. Here, we present a novel framework based on Bayesian linear regression with likelihood warping that allows us to address these problems, that is, to correctly model non-Gaussian predictive distributions and scale normative modelling elegantly to big data cohorts. In addition, this method provides likelihood-based statistics, which are useful for model selection. To evaluate this framework, we use a range of neuroimaging-derived measures from the UK Biobank study, including image-derived phenotypes (IDPs) and whole-brain voxel-wise measures derived from diffusion tensor imaging. We show good computational scaling and improved accuracy of the warped BLR for certain IDPs and voxels if there was a deviation from normality of these parameters in their residuals. The present results indicate the advantage of a warped BLR in terms of; computational scalability and the flexibility to incorporate non-linearity and non-Gaussianity of the data, giving a wider range of neuroimaging datasets that can be correctly modelled.

Replicating extensive brain structural heterogeneity in individuals with schizophrenia and bipolar disorder.

Identifying brain processes involved in the risk and development of mental disorders is a major aim. We recently reported substantial interindividual heterogeneity in brain structural aberrations among patients with schizophrenia and bipolar disorder. Estimating the normative range of voxel-based morphometry (VBM) data among healthy individuals using a Gaussian process regression (GPR) enables us to map individual deviations from the healthy range in unseen datasets. Here, we aim to replicate our previous results in two independent samples of patients with schizophrenia (n1 = 94; n2 = 105), bipolar disorder (n1 = 116; n2 = 61), and healthy individuals (n1 = 400; n2 = 312). In line with previous findings with exception of the cerebellum our results revealed robust group level differences between patients and healthy individuals, yet only a small proportion of patients with schizophrenia or bipolar disorder exhibited extreme negative deviations from normality in the same brain regions. These direct replications support that group level-differences in brain structure disguise considerable individual differences in brain aberrations, with important implications for the interpretation and generalization of group-level brain imaging findings to the individual with a mental disorder.

Modelling brain development to detect white matter injury in term and preterm born neonates.

Premature birth occurs during a period of rapid brain growth. In this context, interpreting clinical neuroimaging can be complicated by the typical changes in brain contrast, size and gyrification occurring in the background to any pathology. To model and describe this evolving background in brain shape and contrast, we used a Bayesian regression technique, Gaussian process regression, adapted to multiple correlated outputs. Using MRI, we simultaneously estimated brain tissue intensity on T1- and T2-weighted scans as well as local tissue shape in a large cohort of 408 neonates scanned cross-sectionally across the perinatal period. The resulting model provided a continuous estimate of brain shape and intensity, appropriate to age at scan, degree of prematurity and sex. Next, we investigated the clinical utility of this model to detect focal white matter injury. In individual neonates, we calculated deviations of a neonate's observed MRI from that predicted by the model to detect punctate white matter lesions with very good accuracy (area under the curve > 0.95). To investigate longitudinal consistency of the model, we calculated model deviations in 46 neonates who were scanned on a second occasion. These infants' voxelwise deviations from the model could be used to identify them from the other 408 images in 83% (T2-weighted) and 76% (T1-weighted) of cases, indicating an anatomical fingerprint. Our approach provides accurate estimates of non-linear changes in brain tissue intensity and shape with clear potential for radiological use.

Heterogeneity in Brain Microstructural Development Following Preterm Birth.

Preterm-born children are at increased risk of lifelong neurodevelopmental difficulties. Group-wise analyses of magnetic resonance imaging show many differences between preterm- and term-born infants but do not reliably predict neurocognitive prognosis for individual infants. This might be due to the unrecognized heterogeneity of cerebral injury within the preterm group. This study aimed to determine whether atypical brain microstructural development following preterm birth is significantly variable between infants. Using Gaussian process regression, a technique that allows a single-individual inference, we characterized typical variation of brain microstructure using maps of fractional anisotropy and mean diffusivity in a sample of 270 term-born neonates. Then, we compared 82 preterm infants to these normative values to identify brain regions with atypical microstructure and relate observed deviations to degree of prematurity and neurocognition at 18 months. Preterm infants showed strikingly heterogeneous deviations from typical development, with little spatial overlap between infants. Greater and more extensive deviations, captured by a whole brain atypicality index, were associated with more extreme prematurity and predicted poorer cognitive and language abilities at 18 months. Brain microstructural development after preterm birth is highly variable between individual infants. This poorly understood heterogeneity likely relates to both the etiology and prognosis of brain injury.

Individual differences v. the average patient: mapping the heterogeneity in ADHD using normative models.

BACKGROUND: The present paper presents a fundamentally novel approach to model individual differences of persons with the same biologically heterogeneous mental disorder. Unlike prevalent case-control analyses, that assume a clear distinction between patient and control groups and thereby introducing the concept of an 'average patient', we describe each patient's biology individually, gaining insights into the different facets that characterize persistent attention-deficit/hyperactivity disorder (ADHD). METHODS: Using a normative modeling approach, we mapped inter-individual differences in reference to normative structural brain changes across the lifespan to examine the degree to which case-control analyses disguise differences between individuals. RESULTS: At the level of the individual, deviations from the normative model were frequent in persistent ADHD. However, the overlap of more than 2% between participants with ADHD was only observed in few brain loci. On average, participants with ADHD showed significantly reduced gray matter in the cerebellum and hippocampus compared to healthy individuals. While the case-control differences were in line with the literature on ADHD, individuals with ADHD only marginally reflected these group differences. CONCLUSIONS: Case-control comparisons, disguise inter-individual differences in brain biology in individuals with persistent ADHD. The present results show that the 'average ADHD patient' has limited informative value, providing the first evidence for the necessity to explore different biological facets of ADHD at the level of the individual and practical means to achieve this end.

Conceptualizing mental disorders as deviations from normative functioning.

Normative models are a class of emerging statistical techniques useful for understanding the heterogeneous biology underlying psychiatric disorders at the level of the individual participant. Analogous to normative growth charts used in paediatric medicine for plotting child development in terms of height or weight as a function of age, normative models chart variation in clinical cohorts in terms of mappings between quantitative biological measures and clinically relevant variables. An emerging body of literature has demonstrated that such techniques are excellent tools for parsing the heterogeneity in clinical cohorts by providing statistical inferences at the level of the individual participant with respect to the normative range. Here, we provide a unifying review of the theory and application of normative modelling for understanding the biological and clinical heterogeneity underlying mental disorders. We first provide a statistically grounded yet non-technical overview of the conceptual underpinnings of normative modelling and propose a conceptual framework to link the many different methodological approaches that have been proposed for this purpose. We survey the literature employing these techniques, focusing principally on applications of normative modelling to quantitative neuroimaging-based biomarkers in psychiatry and, finally, we provide methodological considerations and recommendations to guide future applications of these techniques. We show that normative modelling provides a means by which the importance of modelling individual differences can be brought from theory to concrete data analysis procedures for understanding heterogeneous mental disorders and ultimately a promising route towards precision medicine in psychiatry.

Correction: Conceptualizing mental disorders as deviations from normative functioning.

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

Frequency-specific directed interactions in the human brain network for language.

The brain's remarkable capacity for language requires bidirectional interactions between functionally specialized brain regions. We used magnetoencephalography to investigate interregional interactions in the brain network for language while 102 participants were reading sentences. Using Granger causality analysis, we identified inferior frontal cortex and anterior temporal regions to receive widespread input and middle temporal regions to send widespread output. This fits well with the notion that these regions play a central role in language processing. Characterization of the functional topology of this network, using data-driven matrix factorization, which allowed for partitioning into a set of subnetworks, revealed directed connections at distinct frequencies of interaction. Connections originating from temporal regions peaked at alpha frequency, whereas connections originating from frontal and parietal regions peaked at beta frequency. These findings indicate that the information flow between language-relevant brain areas, which is required for linguistic processing, may depend on the contributions of distinct brain rhythms.

Modeling longitudinal imaging biomarkers with parametric Bayesian multi-task learning.

Longitudinal imaging biomarkers are invaluable for understanding the course of neurodegeneration, promising the ability to track disease progression and to detect disease earlier than cross-sectional biomarkers. To properly realize their potential, biomarker trajectory models must be robust to both under-sampling and measurement errors and should be able to integrate multi-modal information to improve trajectory inference and prediction. Here we present a parametric Bayesian multi-task learning based approach to modeling univariate trajectories across subjects that addresses these criteria. Our approach learns multiple subjects' trajectories within a single model that allows for different types of information sharing, that is, coupling, across subjects. It optimizes a combination of uncoupled, fully coupled and kernel coupled models. Kernel-based coupling allows linking subjects' trajectories based on one or more biomarker measures. We demonstrate this using Alzheimer's Disease Neuroimaging Initiative (ADNI) data, where we model longitudinal trajectories of MRI-derived cortical volumes in neurodegeneration, with coupling based on APOE genotype, cerebrospinal fluid (CSF) and amyloid PET-based biomarkers. In addition to detecting established disease effects, we detect disease related changes within the insula that have not received much attention within the literature. Due to its sensitivity in detecting disease effects, its competitive predictive performance and its ability to learn the optimal parameter covariance from data rather than choosing a specific set of random and fixed effects a priori, we propose that our model can be used in place of or in addition to linear mixed effects models when modeling biomarker trajectories. A software implementation of the method is publicly available.