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.

A neuroimaging measure to capture heterogeneous patterns of atrophy in Parkinson's disease and dementia with Lewy bodies

INTRODUCTION

Parkinson's disease (PD) and Dementia with Lewy bodies (DLB) show heterogeneous brain atrophy patterns which group-average analyses fail to capture. Neuroanatomical normative modelling overcomes this by comparing individuals to a large reference cohort. Patient-specific atrophy patterns are measured objectively and summarised to index overall neurodegeneration (the 'total outlier count'). We aimed to quantify patterns of neurodegenerative dissimilarity in participants with PD and DLB and evaluate the potential clinical relevance of total outlier count by testing its association with key clinical measures in PD and DLB.

MATERIALS AND METHODS

We included 108 participants with PD and 61 with DLB. PD participants were subclassified into high and low visual performers as this has previously been shown to stratify those at increased dementia risk. We generated z-scores from T1w-MRI scans for each participant relative to normative regional cortical thickness and subcortical volumes, modelled in a reference cohort (n = 58,836). Outliers (z < -1.96) were aggregated across 169 brain regions per participant. To measure dissimilarity, individuals' Hamming distance scores were calculated. We also examined total outlier counts between high versus low visual performance in PD; and PD versus DLB; and tested associations between these and cognition.

RESULTS

There was significantly greater inter-individual dissimilarity in brain-outlier patterns in PD poor compared to high visual performers (W = 522.5; p < 0.01) and in DLB compared to PD (W = 5649; p < 0.01). PD poor visual performers had significantly greater total outlier counts compared to high (β = -4.73 (SE = 1.30); t = -3.64; p < 0.01) whereas a conventional group-level GLM failed to identify differences. Higher total outlier counts were associated with poorer MoCA (β = -0.55 (SE = 0.27), t = -2.04, p = 0.05) and composite cognitive scores (β = -2.01 (SE = 0.79); t = -2.54; p = 0.02) in DLB, and visuoperception (β = -0.67 (SE = 0.19); t = -3.59; p < 0.01), in PD.

CONCLUSIONS

Neuroanatomical normative modelling shows promise as a clinically informative technique in PD and DLB, where patterns of atrophy are variable.

Population-wide cerebellar growth models of children and adolescents

In the past, the cerebellum has been best known for its crucial role in motor function. However, increasingly more findings highlight the importance of cerebellar contributions in cognitive functions and neurodevelopment. Using a total of 7240 neuroimaging scans from 4862 individuals, we describe and provide detailed, openly available models of cerebellar development in childhood and adolescence (age range: 6-17 years), an important time period for brain development and onset of neuropsychiatric disorders. Next to a traditionally used anatomical parcellation of the cerebellum, we generated growth models based on a recently proposed functional parcellation. In both, we find an anterior-posterior growth gradient mirroring the age-related improvements of underlying behavior and function, which is analogous to cerebral maturation patterns and offers evidence for directly related cerebello-cortical developmental trajectories. Finally, we illustrate how the current approach can be used to detect cerebellar abnormalities in clinical samples.

Smartphone keyboard dynamics predict affect in suicidal ideation

While digital phenotyping provides opportunities for unobtrusive, real-time mental health assessments, the integration of its modalities is not trivial due to high dimensionalities and discrepancies in sampling frequencies. We provide an integrated pipeline that solves these issues by transforming all modalities to the same time unit, applying temporal independent component analysis (ICA) to high-dimensional modalities, and fusing the modalities with linear mixed-effects models. We applied our approach to integrate high-quality, daily self-report data with BiAffect keyboard dynamics derived from a clinical suicidality sample of mental health outpatients. Applying the ICA to the self-report data (104 participants, 5712 days of data) revealed components related to well-being, anhedonia, and irritability and social dysfunction. Mixed-effects models (55 participants, 1794 days) showed that less phone movement while typing was associated with more anhedonia (β = -0.12, p = 0.00030). We consider this method to be widely applicable to dense, longitudinal digital phenotyping data.

From smartphone data to clinically relevant predictions: A systematic review of digital phenotyping methods in depression

BACKGROUND

Smartphone-based digital phenotyping enables potentially clinically relevant information to be collected as individuals go about their day. This could improve monitoring and interventions for people with Major Depressive Disorder (MDD). The aim of this systematic review was to investigate current digital phenotyping features and methods used in MDD.

METHODS

We searched PubMed, PsycINFO, Embase, Scopus and Web of Science (10/11/2023) for articles including: (1) MDD population, (2) smartphone-based features, (3) validated ratings. Risk of bias was assessed using several sources. Studies were compared within analysis goals (correlating features with depression, predicting symptom severity, diagnosis, mood state/episode, other). Twenty-four studies (9801 participants) were included.

RESULTS

Studies achieved moderate performance. Common themes included challenges from complex and missing data (leading to a risk of bias), and a lack of external validation.

DISCUSSION

Studies made progress towards relating digital phenotypes to clinical variables, often focusing on time-averaged features. Methods investigating temporal dynamics more directly may be beneficial for patient monitoring. European Research Council consolidator grant: 101001118, Prospero: CRD42022346264, Open Science Framework: https://osf.io/s7ay4.

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.

Alzheimer's disease heterogeneity revealed by neuroanatomical normative modeling

INTRODUCTION

Overlooking the heterogeneity in Alzheimer's disease (AD) may lead to diagnostic delays and failures. Neuroanatomical normative modeling captures individual brain variation and may inform our understanding of individual differences in AD-related atrophy.

METHODS

We applied neuroanatomical normative modeling to magnetic resonance imaging from a real-world clinical cohort with confirmed AD (n = 86). Regional cortical thickness was compared to a healthy reference cohort (n = 33,072) and the number of outlying regions was summed (total outlier count) and mapped at individual- and group-levels.

RESULTS

The superior temporal sulcus contained the highest proportion of outliers (60%). Elsewhere, overlap between patient atrophy patterns was low. Mean total outlier count was higher in patients who were non-amnestic, at more advanced disease stages, and without depressive symptoms. Amyloid burden was negatively associated with outlier count.

DISCUSSION

Brain atrophy in AD is highly heterogeneous and neuroanatomical normative modeling can be used to explore anatomo-clinical correlations in individual patients.

Extreme deviations from the normative model reveal cortical heterogeneity and associations with negative symptom severity in first-episode psychosis from the OPTiMiSE and GAP studies

There is currently no quantifiable method to predict long-term clinical outcomes in patients presenting with a first episode of psychosis. A major barrier to developing useful markers for this is biological heterogeneity, where many different pathological mechanisms may underly the same set of symptoms in different individuals. Normative modelling has been used to quantify this heterogeneity in established psychotic disorders by identifying regions of the cortex which are thinner than expected based on a normative healthy population range. These brain atypicalities are measured at the individual level and therefore potentially useful in a clinical setting. However, it is still unclear whether alterations in individual brain structure can be detected at the time of the first psychotic episode, and whether they are associated with subsequent clinical outcomes. We applied normative modelling of cortical thickness to a sample of first-episode psychosis patients, with the aim of quantifying heterogeneity and to use any pattern of cortical atypicality to predict symptoms and response to antipsychotic medication at timepoints from baseline up to 95 weeks (median follow-ups = 4). T1-weighted brain magnetic resonance images from the GAP and OPTiMiSE samples were processed with Freesurfer V6.0.0 yielding 148 cortical thickness features. An existing normative model of cortical thickness (n = 37,126) was adapted to integrate data from each clinical site and account for effects of gender and site. Our test sample consisted of control participants (n = 149, mean age = 26, SD = 6.7) and patient data (n = 295, mean age = 26, SD = 6.7), this sample was used for estimating deviations from the normative model and subsequent statistical analysis. For each individual, the 148 cortical thickness features were mapped to centiles of the normative distribution and converted to z-scores reflecting the distance from the population mean. Individual cortical thickness metrics of +/- 2.6 standard deviations from the mean were considered extreme deviations from the norm. We found that no more than 6.4% of psychosis patients had extreme deviations in a single brain region (regional overlap) demonstrating a high degree of heterogeneity. Mann-Whitney U tests were run on z-scores for each region and significantly lower z-scores were observed in FEP patients in the frontal, temporal, parietal and occipital lobes. Finally, linear mixed-effects modelling showed that negative deviations in cortical thickness in parietal and temporal regions at baseline are related to more severe negative symptoms over the medium-term. This study shows that even at the early stage of symptom onset normative modelling provides a framework to identify individualised cortical markers which can be used for early personalised intervention and stratification.

Smartphone keyboard dynamics predict affect in suicidal ideation

While digital phenotyping provides opportunities for unobtrusive, real-time mental health assessments, the integration of its modalities is not trivial due to high dimensionalities and discrepancies in sampling frequencies. We provide an integrated pipeline that solves these issues by transforming all modalities to the same time unit, applying temporal independent component analysis (ICA) to high-dimensional modalities, and fusing the modalities with linear mixed-effects models. We applied our approach to integrate high-quality, daily self-report data with BiAffect keyboard dynamics derived from a clinical suicidality sample of mental health outpatients. Applying the ICA to the self-report data (104 participants, 5712 days of data) revealed components related to well-being, anhedonia, and irritability and social dysfunction. Mixed-effects models (55 participants, 1794 days) showed that less phone movement while typing was associated with more anhedonia (β = -0.12, p = 0.00030). We consider this method to be widely applicable to dense, longitudinal digital phenotyping data.

Addressing Global Environmental Challenges to Mental Health Using Population Neuroscience: A Review

Importance

Climate change, pollution, urbanization, socioeconomic inequality, and psychosocial effects of the COVID-19 pandemic have caused massive changes in environmental conditions that affect brain health during the life span, both on a population level as well as on the level of the individual. How these environmental factors influence the brain, behavior, and mental illness is not well known.

Observations

A research strategy enabling population neuroscience to contribute to identify brain mechanisms underlying environment-related mental illness by leveraging innovative enrichment tools for data federation, geospatial observation, climate and pollution measures, digital health, and novel data integration techniques is described. This strategy can inform innovative treatments that target causal cognitive and molecular mechanisms of mental illness related to the environment. An example is presented of the environMENTAL Project that is leveraging federated cohort data of over 1.5 million European citizens and patients enriched with deep phenotyping data from large-scale behavioral neuroimaging cohorts to identify brain mechanisms related to environmental adversity underlying symptoms of depression, anxiety, stress, and substance misuse.

Conclusions and Relevance

This research will lead to the development of objective biomarkers and evidence-based interventions that will significantly improve outcomes of environment-related mental illness.

Unpacking the functional heterogeneity of the Emotional Face Matching Task: a normative modelling approach

Functional neuroimaging has contributed substantially to understanding brain function but is dominated by group analyses that index only a fraction of the variation in these data. It is increasingly clear that parsing the underlying heterogeneity is crucial to understand individual differences and the impact of different task manipulations. We estimate large-scale (N=7728) normative models of task-evoked activation during the Emotional Face Matching Task, which enables us to bind heterogeneous datasets to a common reference and dissect heterogeneity underlying group-level analyses. We apply this model to a heterogenous patient cohort, to map individual differences between patients with one or more mental health diagnoses relative to the reference cohort and determine multivariate associations with transdiagnostic symptom domains. For the face>shapes contrast, patients have a higher frequency of extreme deviations which are spatially heterogeneous. In contrast, normative models for faces>baseline have greater predictive value for individuals' transdiagnostic functioning.

A stable and replicable neural signature of lifespan adversity in the adult brain

Environmental adversities constitute potent risk factors for psychiatric disorders. Evidence suggests the brain adapts to adversity, possibly in an adversity-type and region-specific manner. However, the long-term effects of adversity on brain structure and the association of individual neurobiological heterogeneity with behavior have yet to be elucidated. Here we estimated normative models of structural brain development based on a lifespan adversity profile in a longitudinal at-risk cohort aged 25 years (n = 169). This revealed widespread morphometric changes in the brain, with partially adversity-specific features. This pattern was replicated at the age of 33 years (n = 114) and in an independent sample at 22 years (n = 115). At the individual level, greater volume contractions relative to the model were predictive of future anxiety. We show a stable neurobiological signature of adversity that persists into adulthood and emphasize the importance of considering individual-level rather than group-level predictions to explain emerging psychopathology.

Regional, circuit and network heterogeneity of brain abnormalities in psychiatric disorders

The substantial individual heterogeneity that characterizes people with mental illness is often ignored by classical case-control research, which relies on group mean comparisons. Here we present a comprehensive, multiscale characterization of the heterogeneity of gray matter volume (GMV) differences in 1,294 cases diagnosed with one of six conditions (attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, depression, obsessive-compulsive disorder and schizophrenia) and 1,465 matched controls. Normative models indicated that person-specific deviations from population expectations for regional GMV were highly heterogeneous, affecting the same area in <7% of people with the same diagnosis. However, these deviations were embedded within common functional circuits and networks in up to 56% of cases. The salience-ventral attention system was implicated transdiagnostically, with other systems selectively involved in depression, bipolar disorder, schizophrenia and attention-deficit/hyperactivity disorder. Phenotypic differences between cases assigned the same diagnosis may thus arise from the heterogeneous localization of specific regional deviations, whereas phenotypic similarities may be attributable to the dysfunction of common functional circuits and networks.

Fine-grained topographic organization within somatosensory cortex during resting-state and emotional face-matching task and its association with ASD traits

Sensory atypicalities are particularly common in autism spectrum disorders (ASD). Nevertheless, our knowledge about the divergent functioning of the underlying somatosensory region and its association with ASD phenotype features is limited. We applied a data-driven approach to map the fine-grained variations in functional connectivity of the primary somatosensory cortex (S1) to the rest of the brain in 240 autistic and 164 neurotypical individuals from the EU-AIMS LEAP dataset, aged between 7 and 30. We estimated the S1 connection topography ('connectopy') at rest and during the emotional face-matching (Hariri) task, an established measure of emotion reactivity, and accessed its association with a set of clinical and behavioral variables. We first demonstrated that the S1 connectopy is organized along a dorsoventral axis, mapping onto the S1 somatotopic organization. We then found that its spatial characteristics were linked to the individuals' adaptive functioning skills, as measured by the Vineland Adaptive Behavior Scales, across the whole sample. Higher functional differentiation characterized the S1 connectopies of individuals with higher daily life adaptive skills. Notably, we detected significant differences between rest and the Hariri task in the S1 connectopies, as well as their projection maps onto the rest of the brain suggesting a task-modulating effect on S1 due to emotion processing. All in all, variation of adaptive skills appears to be reflected in the brain's mesoscale neural circuitry, as shown by the S1 connectivity profile, which is also differentially modulated during rest and emotional processing.

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.

Personalising Alzheimer's Disease progression using brain atrophy markers

INTRODUCTION

Neuroanatomical normative modelling can capture individual variability in Alzheimer's Disease (AD). We used neuroanatomical normative modelling to track individuals' disease progression in people with mild cognitive impairment (MCI) and patients with AD.

METHODS

Cortical thickness and subcortical volume neuroanatomical normative models were generated using healthy controls (n~58k). These models were used to calculate regional Z-scores in 4361 T1-weighted MRI time-series scans. Regions with Z-scores <-1.96 were classified as outliers and mapped on the brain, and also summarised by total outlier count (tOC).

RESULTS

Rate of change in tOC increased in AD and in people with MCI who converted to AD and correlated with multiple non-imaging markers. Moreover, a higher annual rate of change in tOC increased the risk of MCI progression to AD. Brain Z-score maps showed that the hippocampus had the highest rate of atrophy change.

CONCLUSIONS

Individual-level atrophy rates can be tracked by using regional outlier maps and tOC.

Revealing Individual Neuroanatomical Heterogeneity in Alzheimer Disease Using Neuroanatomical Normative Modeling

BACKGROUND AND OBJECTIVES

Alzheimer disease (AD) is highly heterogeneous, with marked individual differences in clinical presentation and neurobiology. To explore this, we used neuroanatomical normative modeling to index regional patterns of variability in cortical thickness. We aimed to characterize individual differences and outliers in cortical thickness in patients with AD, people with mild cognitive impairment (MCI), and controls. Furthermore, we assessed the relationships between cortical thickness heterogeneity and cognitive function, β-amyloid, phosphorylated-tau, and ApoE genotype. Finally, we examined whether cortical thickness heterogeneity was predictive of conversion from MCI to AD.

METHODS

Cortical thickness measurements across 148 brain regions were obtained from T1-weighted MRI scans from 62 sites of the Alzheimer's Disease Neuroimaging Initiative. AD was determined by clinical and neuropsychological examination with no comorbidities present. Participants with MCI had reported memory complaints, and controls were cognitively normal. A neuroanatomical normative model indexed cortical thickness distributions using a separate healthy reference data set (n = 33,072), which used hierarchical Bayesian regression to predict cortical thickness per region using age and sex, while adjusting for site noise. Z-scores per region were calculated, resulting in a Z-score brain map per participant. Regions with Z-scores <-1.96 were classified as outliers.

RESULTS

Patients with AD (n = 206) had a median of 12 outlier regions (out of a possible 148), with the highest proportion of outliers (47%) in the parahippocampal gyrus. For 62 regions, over 90% of these patients had cortical thicknesses within the normal range. Patients with AD had more outlier regions than people with MCI (n = 662) or controls (n = 159) (F(2, 1,022) = 95.39, p = 2.0 × 10-16). They were also more dissimilar to each other than people with MCI or controls (F(2, 1,024) = 209.42, p = 2.2 × 10-16). A greater number of outlier regions were associated with worse cognitive function, CSF protein concentrations, and an increased risk of converting from MCI to AD within 3 years (hazard ratio 1.028, 95% CI 1.016-1.039, p = 1.8 × 10-16).

DISCUSSION

Individualized normative maps of cortical thickness highlight the heterogeneous effect of AD on the brain. Regional outlier estimates have the potential to be a marker of disease and could be used to track an individual's disease progression or treatment response in clinical trials.

Evidence for absence of links between striatal dopamine synthesis capacity and working memory capacity, spontaneous eye-blink rate, and trait impulsivity

Individual differences in striatal dopamine synthesis capacity have been associated with working memory capacity, trait impulsivity and spontaneous eye-blink rate (sEBR), as measured with readily available and easily administered, 'off-the-shelf' tests. Such findings have raised the suggestion that individual variation in dopamine synthesis capacity, estimated with expensive and invasive brain positron emission tomography (PET) scans, can be approximated with simple, more pragmatic tests. However, direct evidence for the relationship between these simple trait measures and striatal dopamine synthesis capacity has been limited and inconclusive. We measured striatal dopamine synthesis capacity using [¹⁸F]-FDOPA PET in a large sample of healthy volunteers (N=94) and assessed the correlation with simple, short tests of working memory capacity, trait impulsivity, and sEBR. We additionally explored the relationship with an index of subjective reward sensitivity. None of these trait measures correlated significantly with striatal dopamine synthesis capacity, nor did they have out-of-sample predictive power. Bayes Factor analyses indicated the evidence was in favour of absence of correlations for all but subjective reward sensitivity. These results warrant caution for using these off-the-shelf trait measures as proxies of striatal dopamine synthesis capacity.

Evidence for embracing normative modeling

In this work, we expand the normative model repository introduced in Rutherford et al., 2022a to include normative models charting lifespan trajectories of structural surface area and brain functional connectivity, measured using two unique resting-state network atlases (Yeo-17 and Smith-10), and an updated online platform for transferring these models to new data sources. We showcase the value of these models with a head-to-head comparison between the features output by normative modeling and raw data features in several benchmarking tasks: mass univariate group difference testing (schizophrenia versus control), classification (schizophrenia versus control), and regression (predicting general cognitive ability). Across all benchmarks, we show the advantage of using normative modeling features, with the strongest statistically significant results demonstrated in the group difference testing and classification tasks. We intend for these accessible resources to facilitate the wider adoption of normative modeling across the neuroimaging community.

The Link Between Autism and Sex-Related Neuroanatomy, and Associated Cognition and Gene Expression

OBJECTIVE

The male preponderance in prevalence of autism is among the most pronounced sex ratios across neurodevelopmental conditions. The authors sought to elucidate the relationship between autism and typical sex-differential neuroanatomy, cognition, and related gene expression.

METHODS

Using a novel deep learning framework trained to predict biological sex based on T₁-weighted structural brain images, the authors compared sex prediction model performance across neurotypical and autistic males and females. Multiple large-scale data sets comprising T₁-weighted MRI data were employed at four stages of the analysis pipeline: 1) pretraining, with the UK Biobank sample (>10,000 individuals); 2) transfer learning and validation, with the ABIDE data sets (1,412 individuals, 5-56 years of age); 3) test and discovery, with the EU-AIMS/AIMS-2-TRIALS LEAP data set (681 individuals, 6-30 years of age); and 4) specificity, with the NeuroIMAGE and ADHD200 data sets (887 individuals, 7-26 years of age).

RESULTS

Across both ABIDE and LEAP, features positively predictive of neurotypical males were on average significantly more predictive of autistic males (ABIDE: Cohen's d=0.48; LEAP: Cohen's d=1.34). Features positively predictive of neurotypical females were on average significantly less predictive of autistic females (ABIDE: Cohen's d=1.25; LEAP: Cohen's d=1.29). These differences in sex prediction accuracy in autism were not observed in individuals with ADHD. In autistic females, the male-shifted neurophenotype was further associated with poorer social sensitivity and emotional face processing while also associated with gene expression patterns of midgestational cell types.

CONCLUSIONS

The results demonstrate an increased resemblance in both autistic male and female individuals' neuroanatomy with male-characteristic patterns associated with typically sex-differential social cognitive features and related gene expression patterns. The findings hold promise for future research aimed at refining the quest for biological mechanisms underpinning the etiology of autism.

Patterns of connectome variability in autism across five functional activation tasks: findings from the LEAP project

BACKGROUND

Autism spectrum disorder (autism) is a complex neurodevelopmental condition with pronounced behavioral, cognitive, and neural heterogeneities across individuals. Here, our goal was to characterize heterogeneity in autism by identifying patterns of neural diversity as reflected in BOLD fMRI in the way individuals with autism engage with a varied array of cognitive tasks.

METHODS

All analyses were based on the EU-AIMS/AIMS-2-TRIALS multisite Longitudinal European Autism Project (LEAP) with participants with autism (n = 282) and typically developing (TD) controls (n = 221) between 6 and 30 years of age. We employed a novel task potency approach which combines the unique aspects of both resting state fMRI and task-fMRI to quantify task-induced variations in the functional connectome. Normative modelling was used to map atypicality of features on an individual basis with respect to their distribution in neurotypical control participants. We applied robust out-of-sample canonical correlation analysis (CCA) to relate connectome data to behavioral data.

RESULTS

Deviation from the normative ranges of global functional connectivity was greater for individuals with autism compared to TD in each fMRI task paradigm (all tasks p < 0.001). The similarity across individuals of the deviation pattern was significantly increased in autistic relative to TD individuals (p < 0.002). The CCA identified significant and robust brain-behavior covariation between functional connectivity atypicality and autism-related behavioral features.

CONCLUSIONS

Individuals with autism engage with tasks in a globally atypical way, but the particular spatial pattern of this atypicality is nevertheless similar across tasks. Atypicalities in the tasks originate mostly from prefrontal cortex and default mode network regions, but also speech and auditory networks. We show how sophisticated modeling methods such as task potency and normative modeling can be used toward unravelling complex heterogeneous conditions like autism.

Striatal connectopic maps link to functional domains across psychiatric disorders

Transdiagnostic approaches to psychiatry have significant potential in overcoming the limitations of conventional diagnostic paradigms. However, while frameworks such as the Research Domain Criteria have garnered significant enthusiasm among researchers and clinicians from a theoretical angle, examples of how such an approach might translate in practice to understand the biological mechanisms underlying complex patterns of behaviors in realistic and heterogeneous populations have been sparse. In a richly phenotyped clinical sample (n = 186) specifically designed to capture the complex nature of heterogeneity and comorbidity within- and between stress- and neurodevelopmental disorders, we use exploratory factor analysis on a wide range of clinical questionnaires to identify four stable functional domains that transcend diagnosis and relate to negative valence, cognition, social functioning and inhibition/arousal before replicating them in an independent dataset (n = 188). We then use connectopic mapping to map inter-individual variation in fine-grained topographical organization of functional connectivity in the striatum-a central hub in motor, cognitive, affective and reward-related brain circuits-and use multivariate machine learning (canonical correlation analysis) to show that these individualized topographic representations predict transdiagnostic functional domains out of sample (r = 0.20, p = 0.026). We propose that investigating psychiatric symptoms across disorders is a promising path to linking them to underlying biology, and can help bridge the gap between neuroscience and clinical psychiatry.

Closing the life-cycle of normative modeling using federated hierarchical Bayesian regression

Clinical neuroimaging data availability has grown substantially in the last decade, providing the potential for studying heterogeneity in clinical cohorts on a previously unprecedented scale. Normative modeling is an emerging statistical tool for dissecting heterogeneity in complex brain disorders. However, its application remains technically challenging due to medical data privacy issues and difficulties in dealing with nuisance variation, such as the variability in the image acquisition process. Here, we approach the problem of estimating a reference normative model across a massive population using a massive multi-center neuroimaging dataset. To this end, we introduce a federated probabilistic framework using hierarchical Bayesian regression (HBR) to complete the life-cycle of normative modeling. The proposed model provides the possibilities to learn, update, and adapt the model parameters on decentralized neuroimaging data. Our experimental results confirm the superiority of HBR in deriving more accurate normative ranges on large multi-site neuroimaging datasets compared to the current standard methods. In addition, our approach provides the possibility to recalibrate and reuse the learned model on local datasets and even on datasets with very small sample sizes. The proposed method will facilitate applications of normative modeling as a medical tool for screening the biological deviations in individuals affected by complex illnesses such as mental disorders.

Deviations from normative brain white and gray matter structure are associated with psychopathology in youth

Combining imaging modalities and metrics that are sensitive to various aspects of brain structure and maturation may help identify individuals that show deviations in relation to same-aged peers, and thus benefit early-risk-assessment for mental disorders. We used one timepoint multimodal brain imaging, cognitive, and questionnaire data from 1280 eight- to twenty-one-year-olds from the Philadelphia Neurodevelopmental Cohort. We estimated age-related gray and white matter properties and estimated individual deviation scores using normative modeling. Next, we tested for associations between the estimated deviation scores, and with psychopathology domain scores and cognition. More negative deviations in DTI-based fractional anisotropy (FA) and the first principal eigenvalue of the diffusion tensor (L1) were associated with higher scores on psychosis positive and prodromal symptoms and general psychopathology. A more negative deviation in cortical thickness (CT) was associated with a higher general psychopathology score. Negative deviations in global FA, surface area, L1 and CT were also associated with poorer cognitive performance. No robust associations were found between the deviation scores based on CT and DTI. The low correlations between the different multimodal magnetic resonance imaging-based deviation scores suggest that psychopathological burden in adolescence can be mapped onto partly distinct neurobiological features.

Contrasting Case-Control and Normative Reference Approaches to Capture Clinically Relevant Structural Brain Abnormalities in Patients With First-Episode Psychosis Who Are Antipsychotic Naive

Importance

To make progress toward precision psychiatry, it is crucial to move beyond case-control studies and instead capture individual variations and interpret them in the context of a normal range of biological systems.

Objective

To evaluate whether baseline deviations from a normative reference range in subcortical volumes are better predictors of antipsychotic treatment response than raw volumes in patients with first-episode psychosis (FEP) who were naive to antipsychotic medication.

Design, Setting, and Participants

In this prospective longitudinal study, patients with first-episode psychosis who were referred from different clinical settings (emergency department, inpatient units, and outpatient clinics) at the University of Alabama at Birmingham were included. A total of 286 patients were screened, 114 consented, 104 enrolled in the treatment trial, and 85 completed the trial. Patients were observed for 16 weeks. Controls were matched by age and sex. Data were collected between June 2016 and July 2021, and data were analyzed from August 2021 to June 2022.

Interventions

Risperidone on a flexible dosing scheme for 16 weeks. There was an option to switch to aripiprazole for excessive adverse effects.

Main Outcomes and Measures

The main outcome of this study was to evaluate, in patients with FEP who were naive to antipsychotic medication, the association of baseline raw volumes and volume deviations in subcortical brain regions with response to antipsychotic medication. Raw brain volumes or volume deviation changes after treatment were not examined.

Results

Of 190 included participants, 111 (58.4%) were male, and the mean (SD) age was 23.7 (5.5) years. Volumes and deviations were quantified in 98 patients with FEP, and data from 92 controls were used as comparison for case-control contrasts and reference curve calibration. In case-control contrasts, patients with FEP had lower raw thalamus (P = .002; F = 9.63; df = 1), hippocampus (P = .009; F = 17.23; df = 1), amygdala (P = .01; F = 6.55; df = 1), ventral diencephalon (P = .03; F = 4.84; df = 1), and brainstem volumes (P = .004; F = 8.39; df = 1). Of 98 patients, 36 patients with FEP (36%) displayed extreme deviations. Associations with treatment response significantly differed between raw volume and deviation measures in the caudate (z = -2.17; P = .03) and putamen (z = -2.15; P = .03).

Conclusions and Relevance

These data suggest that normative modeling allows capture of interindividual heterogeneity of regional brain volumes in patients with FEP and characterize structural pathology in a clinically relevant fashion. This holds promise for progress in precision medicine in psychiatry, where group-level studies have failed to derive reliable maps of structural pathology.

Cerebellar Atypicalities in Autism?

BACKGROUND

The cerebellum contains more than 50% of the brain's neurons and is involved in social cognition. Cerebellar anatomical atypicalities have repeatedly been reported in individuals with autism. However, studies have yielded inconsistent findings, likely because of a lack of statistical power, and did not capture the clinical and neuroanatomical diversity of autism. Our aim was to better understand cerebellar anatomy and its diversity in autism.

METHODS

We studied cerebellar gray matter morphology in 274 individuals with autism and 219 control subjects of a multicenter European cohort, EU-AIMS LEAP (European Autism Interventions-A Multicentre Study for Developing New Medications; Longitudinal European Autism Project). To ensure the robustness of our results, we conducted lobular parcellation of the cerebellum with 2 different pipelines in addition to voxel-based morphometry. We performed statistical analyses with linear, multivariate (including normative modeling), and meta-analytic approaches to capture the diversity of cerebellar anatomy in individuals with autism and control subjects. Finally, we performed a dimensional analysis of cerebellar anatomy in an independent cohort of 352 individuals with autism-related symptoms.

RESULTS

We did not find any significant difference in the cerebellum when comparing individuals with autism and control subjects using linear models. In addition, there were no significant deviations in our normative models in the cerebellum in individuals with autism. Finally, we found no evidence of cerebellar atypicalities related to age, IQ, sex, or social functioning in individuals with autism.

CONCLUSIONS

Despite positive results published in the last decade from relatively small samples, our results suggest that there is no striking difference in cerebellar anatomy of individuals with autism.

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.

The normative modeling framework for computational psychiatry

Normative modeling is an emerging and innovative framework for mapping individual differences at the level of a single subject or observation in relation to a reference model. It involves charting centiles of variation across a population in terms of mappings between biology and behavior, which can then be used to make statistical inferences at the level of the individual. The fields of computational psychiatry and clinical neuroscience have been slow to transition away from patient versus 'healthy' control analytic approaches, probably owing to a lack of tools designed to properly model biological heterogeneity of mental disorders. Normative modeling provides a solution to address this issue and moves analysis away from case-control comparisons that rely on potentially noisy clinical labels. Here we define a standardized protocol to guide users through, from start to finish, normative modeling analysis using the Predictive Clinical Neuroscience toolkit (PCNtoolkit). We describe the input data selection process, provide intuition behind the various modeling choices and conclude by demonstrating several examples of downstream analyses that the normative model may facilitate, such as stratification of high-risk individuals, subtyping and behavioral predictive modeling. The protocol takes ~1-3 h to complete.

Mapping dopaminergic projections in the human brain with resting-state fMRI

The striatum receives dense dopaminergic projections, making it a key region of the dopaminergic system. Its dysfunction has been implicated in various conditions including Parkinson’s disease (PD) and substance use disorder. However, the investigation of dopamine-specific functioning in humans is problematic as current MRI approaches are unable to differentiate between dopaminergic and other projections. Here, we demonstrate that ‘connectopic mapping’ – a novel approach for characterizing fine-grained, overlapping modes of functional connectivity – can be used to map dopaminergic projections in striatum. We applied connectopic mapping to resting-state functional MRI data of the Human Connectome Project (population cohort; N = 839) and selected the second-order striatal connectivity mode for further analyses. We first validated its specificity to dopaminergic projections by demonstrating a high spatial correlation (r = 0.884) with dopamine transporter availability – a marker of dopaminergic projections – derived from DaT SPECT scans of 209 healthy controls. Next, we obtained the subject-specific second-order modes from 20 controls and 39 PD patients scanned under placebo and under dopamine replacement therapy (L-DOPA), and show that our proposed dopaminergic marker tracks PD diagnosis, symptom severity, and sensitivity to L-DOPA. Finally, across 30 daily alcohol users and 38 daily smokers, we establish strong associations with self-reported alcohol and nicotine use. Our findings provide evidence that the second-order mode of functional connectivity in striatum maps onto dopaminergic projections, tracks inter-individual differences in PD symptom severity and L-DOPA sensitivity, and exhibits strong associations with levels of nicotine and alcohol use, thereby offering a new biomarker for dopamine-related (dys)function in the human brain.

Charting brain growth and aging at high spatial precision

Defining reference models for population variation, and the ability to study individual deviations is essential for understanding inter-individual variability and its relation to the onset and progression of medical conditions. In this work, we assembled a reference cohort of neuroimaging data from 82 sites (N=58,836; ages 2-100) and used normative modeling to characterize lifespan trajectories of cortical thickness and subcortical volume. Models are validated against a manually quality checked subset (N=24,354) and we provide an interface for transferring to new data sources. We showcase the clinical value by applying the models to a transdiagnostic psychiatric sample (N=1985), showing they can be used to quantify variability underlying multiple disorders whilst also refining case-control inferences. These models will be augmented with additional samples and imaging modalities as they become available. This provides a common reference platform to bind results from different studies and ultimately paves the way for personalized clinical decision-making.

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 T₁ -weighted MRI scans. Subsequently, we apply the selected rapid processing pipeline both to rapidly acquired multicontrast EPImix scans of 95 participants (which include T₁ -FLAIR, T₂ , T₂ *, T₂ -FLAIR, DWI and ADC contrasts, acquired in ~1 min), as well as to slower, more standard single-contrast T₁ -weighted scans of a subset of 66 participants. We quantify the correspondence between EPImix T₁ -FLAIR and single-contrast T₁ -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.

Brain aging in major depressive disorder: results from the ENIGMA major depressive disorder working group

Major depressive disorder (MDD) is associated with an increased risk of brain atrophy, aging-related diseases, and mortality. We examined potential advanced brain aging in adult MDD patients, and whether this process is associated with clinical characteristics in a large multicenter international dataset. We performed a mega-analysis by pooling brain measures derived from T1-weighted MRI scans from 19 samples worldwide. Healthy brain aging was estimated by predicting chronological age (18-75 years) from 7 subcortical volumes, 34 cortical thickness and 34 surface area, lateral ventricles and total intracranial volume measures separately in 952 male and 1236 female controls from the ENIGMA MDD working group. The learned model coefficients were applied to 927 male controls and 986 depressed males, and 1199 female controls and 1689 depressed females to obtain independent unbiased brain-based age predictions. The difference between predicted "brain age" and chronological age was calculated to indicate brain-predicted age difference (brain-PAD). On average, MDD patients showed a higher brain-PAD of +1.08 (SE 0.22) years (Cohen's d = 0.14, 95% CI: 0.08-0.20) compared with controls. However, this difference did not seem to be driven by specific clinical characteristics (recurrent status, remission status, antidepressant medication use, age of onset, or symptom severity). This highly powered collaborative effort showed subtle patterns of age-related structural brain abnormalities in MDD. Substantial within-group variance and overlap between groups were observed. Longitudinal studies of MDD and somatic health outcomes are needed to further assess the clinical value of these brain-PAD estimates.

The Normative Modeling Framework for Computational Psychiatry.. [PMID: 35650452 PMCID: PMC7613648 DOI: 10.1101/2021.08.08.455583]

Normative modeling is an emerging and innovative framework for mapping individual differences at the level of a single subject or observation in relation to a reference model. It involves charting centiles of variation across a population in terms of mappings between biology and behavior which can then be used to make statistical inferences at the level of the individual. The fields of computational psychiatry and clinical neuroscience have been slow to transition away from patient versus “healthy” control analytic approaches, likely due to a lack of tools designed to properly model biological heterogeneity of mental disorders. Normative modeling provides a solution to address this issue and moves analysis away from case-control comparisons that rely on potentially noisy clinical labels. In this article, we define a standardized protocol to guide users through, from start to finish, normative modeling analysis using the Predictive Clinical Neuroscience toolkit (PCNtoolkit). We describe the input data selection process, provide intuition behind the various modeling choices, and conclude by demonstrating several examples of down-stream analyses the normative model results may facilitate, such as stratification of high-risk individuals, subtyping, and behavioral predictive modeling. The protocol takes approximately 1-3 hours to complete.

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.

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.

Warped Bayesian Linear Regression for Normative Modelling of Big Data.. [PMID: 34798518 PMCID: PMC7613680 DOI: 10.1101/2021.04.05.438429]

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 scale normative modelling elegantly to big data cohorts and to correctly model non-Gaussian predictive distributions. In addition, this method provides also 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.

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.

Brain Connectome Mapping of Complex Human Traits and Their Polygenic Architecture Using Machine Learning

BACKGROUND

Mental disorders and individual characteristics such as intelligence and personality are complex traits sharing a largely unknown neuronal basis. Their genetic architectures are highly polygenic and overlapping, which is supported by heterogeneous phenotypic expression and substantial clinical overlap. Brain network analysis provides a noninvasive means of dissecting biological heterogeneity, yet its sensitivity, specificity, and validity in assessing individual characteristics relevant for brain function and mental health and their genetic underpinnings in clinical applications remain a challenge.

METHODS

In a machine learning approach, we predicted individual scores for educational attainment, fluid intelligence and dimensional measures of depression, anxiety, and neuroticism using functional magnetic resonance imaging-based static and dynamic temporal synchronization between large-scale brain network nodes in 10,343 healthy individuals from the UK Biobank. In addition to using age and sex to serve as our reference point, we also predicted individual polygenic scores for related phenotypes and 13 different neuroticism traits and schizophrenia.

RESULTS

Beyond high accuracy for age and sex, supporting the biological sensitivity of the connectome-based features, permutation tests revealed above chance-level prediction accuracy for trait-level educational attainment and fluid intelligence. Educational attainment and fluid intelligence were mainly negatively associated with static brain connectivity in frontal and default mode networks, whereas age showed positive correlations with a more widespread pattern. In contrast, prediction accuracy was at chance level for depression, anxiety, neuroticism, and polygenic scores across traits.

CONCLUSIONS

These novel findings provide a benchmark for future studies linking the genetic architecture of individual and mental health traits with functional magnetic resonance imaging-based brain connectomics.

Analysing brain networks in population neuroscience: a case for the Bayesian philosophy

Network connectivity fingerprints are among today's best choices to obtain a faithful sampling of an individual's brain and cognition. Widely available MRI scanners can provide rich information tapping into network recruitment and reconfiguration that now scales to hundreds and thousands of humans. Here, we contemplate the advantages of analysing such connectome profiles using Bayesian strategies. These analysis techniques afford full probability estimates of the studied network coupling phenomena, provide analytical machinery to separate epistemological uncertainty and biological variability in a coherent manner, usher us towards avenues to go beyond binary statements on existence versus non-existence of an effect, and afford credibility estimates around all model parameters at play which thus enable single-subject predictions with rigorous uncertainty intervals. We illustrate the brittle boundary between healthy and diseased brain circuits by autism spectrum disorder as a recurring theme where, we argue, network-based approaches in neuroscience will require careful probabilistic answers. This article is part of the theme issue 'Unifying the essential concepts of biological networks: biological insights and philosophical foundations'.

The genetic architecture of the human cerebral cortex

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder.

10,000 social brains: Sex differentiation in human brain anatomy

In human and nonhuman primates, sex differences typically explain much interindividual variability. Male and female behaviors may have played unique roles in the likely coevolution of increasing brain volume and more complex social dynamics. To explore possible divergence in social brain morphology between men and women living in different social environments, we applied probabilistic generative modeling to ~10,000 UK Biobank participants. We observed strong volume effects especially in the limbic system but also in regions of the sensory, intermediate, and higher association networks. Sex-specific brain volume effects in the limbic system were linked to the frequency and intensity of social contact, such as indexed by loneliness, household size, and social support. Across the processing hierarchy of neural networks, different conditions for social interplay may resonate in and be influenced by brain anatomy in sex-dependent ways.

Predicting individual clinical trajectories of depression with generative embedding

•

Patients with major depressive disorder (MDD) show variable clinical trajectories.

•

Generative embedding (GE) is used to predict clinical trajectories in MDD patients.

•

GE classifies patients with chronic depression vs. fast remission with 79% accuracy.

•

GE provides mechanistic interpretability and outperforms conventional measures.

•

Proof-of-concept that illustrates the potential of GE for clinical prediction.

Patients with major depressive disorder (MDD) show heterogeneous treatment response and highly variable clinical trajectories: while some patients experience swift recovery, others show relapsing-remitting or chronic courses. Predicting individual clinical trajectories at an early stage is a key challenge for psychiatry and might facilitate individually tailored interventions. So far, however, reliable predictors at the single-patient level are absent. Here, we evaluated the utility of a machine learning strategy – generative embedding (GE) – which combines interpretable generative models with discriminative classifiers. Specifically, we used functional magnetic resonance imaging (fMRI) data of emotional face perception in 85 MDD patients from the NEtherlands Study of Depression and Anxiety (NESDA) who had been followed up over two years and classified into three subgroups with distinct clinical trajectories. Combining a generative model of effective (directed) connectivity with support vector machines (SVMs), we could predict whether a given patient would experience chronic depression vs. fast remission with a balanced accuracy of 79%. Gradual improvement vs. fast remission could still be predicted above-chance, but less convincingly, with a balanced accuracy of 61%. Generative embedding outperformed classification based on conventional (descriptive) features, such as functional connectivity or local activation estimates, which were obtained from the same data and did not allow for above-chance classification accuracy. Furthermore, predictive performance of GE could be assigned to a specific network property: the trial-by-trial modulation of connections by emotional content. Given the limited sample size of our study, the present results are preliminary but may serve as proof-of-concept, illustrating the potential of GE for obtaining clinical predictions that are interpretable in terms of network mechanisms. Our findings suggest that abnormal dynamic changes of connections involved in emotional face processing might be associated with higher risk of developing a less favorable clinical course.

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.

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.