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

Individualized functional connectome abnormalities obtained using two normative model unveil neurophysiological subtypes of obsessive compulsive disorder

The high heterogeneity observed among patients with obsessive-compulsive disorder (OCD) underscores the need to identify neurophysiological OCD subtypes to facilitate personalized diagnosis and treatment. In this study, our aim was to identify potential OCD subtypes based on individualized functional connectome abnormalities. We recruited a total of 99 patients with OCD and 104 healthy controls (HCs) matched for demographic characteristics. Individualized functional connectome abnormalities were obtained using normative models of functional connectivity strength (FCS) and used as features to unveil OCD subtypes. Sensitivity analyses were conducted to assess the reproducibility and robustness of the clustering outcomes. Patients exhibited significant intersubject heterogeneity in individualized functional connectome abnormalities. Two subtypes with distinct patterns of FCS abnormalities relative to HCs were identified. Subtype 1 patients primarily exhibited significantly decreased FCS in regions including the frontal gyrus, insula, hippocampus, and precentral/postcentral gyrus, whereas subtype 2 patients demonstrated increased FCS in widespread brain regions. When all patients were combined, no significant differences were observed. Additionally, the identified subtypes showed significant differences in age of onset. Furthermore, sensitivity analyses confirmed the reproducibility of our subtyping results. In conclusion, the identified OCD subtypes shed new light on the taxonomy and neurophysiological heterogeneity of OCD.

Maternal psychopathology is differentially associated with adolescent offspring neural response to reward given offspring ADHD risk

Reinforcement sensitivity is a hypothesized attention-deficit/hyperactivity disorder (ADHD) intermediate phenotype but its role in transgenerational transmission of ADHD-linked psychopathology risk is largely unknown. We examined, in a carefully phenotyped, N = 123 sample of adolescents (Mage = 15.27 years, SD = 0.984; 61.78% boys), whether (1) parental psychopathology is differentially associated with fMRI-indexed neural response to reward receipt and (2) both maternal and paternal psychopathology are associated with neural response to reward; across adolescents at-risk for and not at-risk for ADHD. Indices of parental psychopathology were differentially associated with adolescent offspring neural response to reward such that across measures, parental psychopathology was negatively or not associated with offspring superior frontal gyrus (SFG) response to reward receipt in adolescents at-risk for ADHD, but parental psychopathology was positively associated with offspring SFG response in adolescents not at-risk. Further, across measures, greater maternal psychopathology was associated with blunted adolescent SFG response to reward in adolescents at-risk for ADHD whereas greater maternal externalizing problems were linked to enhanced adolescent SFG response in adolescents not at-risk. Across measures, paternal psychopathology was not associated with adolescent response to reward, in either group. ADHD risk confers differential reward-related susceptibility to the effects of parental psychopathology. Results also show this association is nonspecific in terms of parental psychopathology type but is specific to maternal psychopathology.

COVID-19 lockdown effects on adolescent brain structure suggest accelerated maturation that is more pronounced in females than in males

Adolescence is a period of substantial social-emotional development, accompanied by dramatic changes to brain structure and function. Social isolation due to lockdowns that were imposed because of the COVID-19 pandemic had a detrimental impact on adolescent mental health, with the mental health of females more affected than males. We assessed the impact of the COVID-19 pandemic lockdowns on adolescent brain structure with a focus on sex differences. We collected MRI structural data longitudinally from adolescents prior to and after the pandemic lockdowns. The pre-COVID data were used to create a normative model of cortical thickness change with age during typical adolescent development. Cortical thickness values in the post-COVID data were compared to this normative model. The analysis revealed accelerated cortical thinning in the post-COVID brain, which was more widespread throughout the brain and greater in magnitude in females than in males. When measured in terms of equivalent years of development, the mean acceleration was found to be 4.2 y in females and 1.4 y in males. Accelerated brain maturation as a result of chronic stress or adversity during development has been well documented. These findings suggest that the lifestyle disruptions associated with the COVID-19 pandemic lockdowns caused changes in brain biology and had a more severe impact on the female than the male brain.

Disentangling the role of gray matter volume and concentration in autism spectrum disorder: A meta-analytic investigation of 25 years of voxel-based morphometry research

Despite over two decades of neuroimaging research, a unanimous definition of the pattern of structural variation associated with autism spectrum disorder (ASD) has yet to be found. One potential impeding issue could be the sometimes ambiguous use of measurements of variations in gray matter volume (GMV) or gray matter concentration (GMC). In fact, while both can be calculated using voxel-based morphometry analysis, these may reflect different underlying pathological mechanisms. We conducted a coordinate-based meta-analysis, keeping apart GMV and GMC studies of subjects with ASD. Results showed distinct and non-overlapping patterns for the two measures. GMV decreases were evident in the cerebellum, while GMC decreases were mainly found in the temporal and frontal regions. GMV increases were found in the parietal, temporal, and frontal brain regions, while GMC increases were observed in the anterior cingulate cortex and middle frontal gyrus. Age-stratified analyses suggested that such variations are dynamic across the ASD lifespan. The present findings emphasize the importance of considering GMV and GMC as distinct yet synergistic indices in autism research.

Behavioral and brain morphological changes before and after hemispherotomy

Hemispherotomy is an effective surgery for treating refractory epilepsy from diffuse unihemispheric lesions. To date, postsurgery neuroplastic changes supporting behavioral recovery after left or right hemispherotomy remain unclear. In the present study, we systematically investigated changes in gray matter volume (GMV) before and after surgery and further analyzed their relationships with behavioral scores in two large groups of pediatric patients with left and right hemispherotomy (29 left and 28 right). To control for the dramatic developmental effect during this stage, age-adjusted GMV within unaffected brain regions was derived voxel by voxel using a normative modeling approach with an age-matched reference cohort of 2115 healthy children. Widespread GMV increases in the contralateral cerebrum and ipsilateral cerebellum and GMV decreases in the contralateral cerebellum were consistently observed in both patient groups, but only the left hemispherotomy patients showed GMV decreases in the contralateral cingulate gyrus. Intriguingly, the GMV decrease in the contralateral cerebellum was significantly correlated with improvement in behavioral scores in the right but not the left hemispherotomy patients. Importantly, the preoperative voxelwise GMV features can be used to significantly predict postoperative behavioral scores in both patient groups. These findings indicate an important role of the contralateral cerebellum in the behavioral recovery following right hemispherotomy and highlight the predictive potential of preoperative imaging features in postoperative behavioral performance.

Dataset factors associated with age-related changes in brain structure and function in neurodevelopmental conditions

With brain structure and function undergoing complex changes throughout childhood and adolescence, age is a critical consideration in neuroimaging studies, particularly for those of individuals with neurodevelopmental conditions. However, despite the increasing use of large, consortium-based datasets to examine brain structure and function in neurotypical and neurodivergent populations, it is unclear whether age-related changes are consistent between datasets and whether inconsistencies related to differences in sample characteristics, such as demographics and phenotypic features, exist. To address this, we built models of age-related changes of brain structure (regional cortical thickness and regional surface area; N = 1218) and function (resting-state functional connectivity strength; N = 1254) in two neurodiverse datasets: the Province of Ontario Neurodevelopmental Network and the Healthy Brain Network. We examined whether deviations from these models differed between the datasets, and explored whether these deviations were associated with demographic and clinical variables. We found significant differences between the two datasets for measures of cortical surface area and functional connectivity strength throughout the brain. For regional measures of cortical surface area, the patterns of differences were associated with race/ethnicity, while for functional connectivity strength, positive associations were observed with head motion. Our findings highlight that patterns of age-related changes in the brain may be influenced by demographic and phenotypic characteristics, and thus future studies should consider these when examining or controlling for age effects in analyses.

Schizophrenia and Neurodevelopment: Insights From Connectome Perspective

BACKGROUND

Schizophrenia is conceptualized as a brain connectome disorder that can emerge as early as late childhood and adolescence. However, the underlying neurodevelopmental basis remains unclear. Recent interest has grown in children and adolescent patients who experience symptom onset during critical brain development periods. Inspired by advanced methodological theories and large patient cohorts, Chinese researchers have made significant original contributions to understanding altered brain connectome development in early-onset schizophrenia (EOS).

STUDY DESIGN

We conducted a search of PubMed and Web of Science for studies on brain connectomes in schizophrenia and neurodevelopment. In this selective review, we first address the latest theories of brain structural and functional development. Subsequently, we synthesize Chinese findings regarding mechanisms of brain structural and functional abnormalities in EOS. Finally, we highlight several pivotal challenges and issues in this field.

STUDY RESULTS

Typical neurodevelopment follows a trajectory characterized by gray matter volume pruning, enhanced structural and functional connectivity, improved structural connectome efficiency, and differentiated modules in the functional connectome during late childhood and adolescence. Conversely, EOS deviates with excessive gray matter volume decline, cortical thinning, reduced information processing efficiency in the structural brain network, and dysregulated maturation of the functional brain network. Additionally, common functional connectome disruptions of default mode regions were found in early- and adult-onset patients.

CONCLUSIONS

Chinese research on brain connectomes of EOS provides crucial evidence for understanding pathological mechanisms. Further studies, utilizing standardized analyses based on large-sample multicenter datasets, have the potential to offer objective markers for early intervention and disease treatment.

Dissecting Heterogeneity in Functional Network Connectivity Aberrations in Antipsychotic Medication-Naïve First Episode Psychosis Patients - A Normative Modeling Study

Importance

While there is a general consensus that functional connectome pathology is a key mechanism underlying psychosis spectrum disorders, the literature is plagued with inconsistencies and translation into clinical practice is non-existent. This is perhaps because group-level findings may not be accurate reflections of pathology at the individual patient level.

Objective

To characterize inter-individual heterogeneity in functional networks and investigate if normative values can be leveraged to identify biologically less heterogeneous subgroups of patients.

Design Setting and Participants

We used data collected in a case-control study conducted at the University of Alabama at Birmingham (UAB). We recruited antipsychotic medication-naïve first-episode psychosis patients from UAB outpatient, inpatient, and emergency room settings.

Main Outcomes and Measures

Individual-level patterns of deviations from a normative reference range in resting-state functional networks using the Yeo-17 atlas for parcellations.

Results

Statistical analyses included 108 medication-naïve first-episode psychosis patients. We found that there is a high level of inter-individual heterogeneity in resting-state network connectivity deviations from the normative reference range. Interestingly 48% of patients did not have any functional connectivity deviations, and no more than 11.1% of patients shared functional deviations between the same regions of interest. In a post hoc analysis, we grouped patients based on deviations into four theoretically possible groups. We discovered that all four groups do exist in our experimental data and showed that subgroups based on deviation profiles were significantly less heterogeneous compared to the overall group (positive deviation group: z= -2.88, p = 0.002; negative deviation group: z= -3.36, p<0.001).

Conclusions and Relevance

Our findings experimentally demonstrate that there is a high level of inter-individual heterogeneity in resting-state network pathology in first-episode psychosis patients which support the idea that group-level findings are not accurate reflections of pathology at the individual level. We also demonstrated that normative functional connectivity deviations may have utility for identifying biologically less heterogeneous subgroups of patients, even though they are not distinguishable clinically. Our findings constitute a significant step towards making precision psychiatry a reality, where patients are selected for treatments based on their individual biological characteristics.

KEY POINTS

Question: How heterogeneous is individual-level resting-state functional network pathology in patients suffering from a first psychotic episode? Can normative reference values in functional network connectivity be leveraged to identify biologically more homogenous subgroups of patients?Findings: We report that functional network pathology is highly heterogeneous, with no more than 11% of patients sharing functional deviations between the same regions of interest.Meaning: Normative modeling is a tool that can map individual neurobiological differences and enables the classification of a clinically heterogenous patient group into subgroups that are neurobiologically less heterogenous.

The future of disentangling the heterogeneity of autism with neuroimaging studies

Autism spectrum disorder (ASD) is a lifelong neurodevelopmental condition. Over the past decade, a considerable number of approaches have been developed to identify potential neuroimaging-based biomarkers of ASD and have uncovered specific neural mechanisms underlying behaviors associated with ASD. However, the substantial heterogeneity among those diagnosed with ASD hinders the development of biomarkers. Disentangling the heterogeneity of ASD is pivotal to improve quality of life for individuals with ASD by facilitating early diagnosis and individualized interventions for those who need support. In this Review, we discuss recent advances in neuroimaging that have facilitated the characterization of the heterogeneity of this condition from three frameworks: neurosubtyping, dimensional models, and normative models. In addition, we discuss the challenges, possible solutions, and clinical utility of these three frameworks. We argue that several factors need to be considered when parsing heterogeneity using neuroimaging, including co-occurring conditions, neurodevelopment, heredity and environment, and multi-site and multi-modality data. We close with a discussion of future directions for achieving a better understanding of the neural mechanisms underlying neurodevelopmental heterogeneity and the future of precision medicine in ASD.

Multiscale heterogeneity of white matter morphometry in psychiatric disorders

Background

Inter-individual variability in neurobiological and clinical characteristics in mental illness is often overlooked by classical group-mean case-control studies. Studies using normative modelling to infer person-specific deviations of grey matter volume have indicated that group means are not representative of most individuals. The extent to which this variability is present in white matter morphometry, which is integral to brain function, remains unclear.

Methods

We applied Warped Bayesian Linear Regression normative models to T1-weighted magnetic resonance imaging data and mapped inter-individual variability in person-specific white matter volume deviations in 1,294 cases (58% male) diagnosed with one of six disorders (attention-deficit/hyperactivity, autism, bipolar, major depressive, obsessive-compulsive and schizophrenia) and 1,465 matched controls (54% male) recruited across 25 scan sites. We developed a framework to characterize deviation heterogeneity at multiple spatial scales, from individual voxels, through inter-regional connections, specific brain regions, and spatially extended brain networks.

Results

The specific locations of white matter volume deviations were highly heterogeneous across participants, affecting the same voxel in fewer than 8% of individuals with the same diagnosis. For autism and schizophrenia, negative deviations (i.e., areas where volume is lower than normative expectations) aggregated into common tracts, regions and large-scale networks in up to 35% of individuals.

Conclusions

The prevalence of white matter volume deviations was lower than previously observed in grey matter, and the specific location of these deviations was highly heterogeneous when considering voxel-wise spatial resolution. Evidence of aggregation within common pathways and networks was apparent in schizophrenia and autism but not other disorders.

Brain health in diverse settings: How age, demographics and cognition shape brain function

Diversity in brain health is influenced by individual differences in demographics and cognition. However, most studies on brain health and diseases have typically controlled for these factors rather than explored their potential to predict brain signals. Here, we assessed the role of individual differences in demographics (age, sex, and education; n = 1298) and cognition (n = 725) as predictors of different metrics usually used in case-control studies. These included power spectrum and aperiodic (1/f slope, knee, offset) metrics, as well as complexity (fractal dimension estimation, permutation entropy, Wiener entropy, spectral structure variability) and connectivity (graph-theoretic mutual information, conditional mutual information, organizational information) from the source space resting-state EEG activity in a diverse sample from the global south and north populations. Brain-phenotype models were computed using EEG metrics reflecting local activity (power spectrum and aperiodic components) and brain dynamics and interactions (complexity and graph-theoretic measures). Electrophysiological brain dynamics were modulated by individual differences despite the varied methods of data acquisition and assessments across multiple centers, indicating that results were unlikely to be accounted for by methodological discrepancies. Variations in brain signals were mainly influenced by age and cognition, while education and sex exhibited less importance. Power spectrum activity and graph-theoretic measures were the most sensitive in capturing individual differences. Older age, poorer cognition, and being male were associated with reduced alpha power, whereas older age and less education were associated with reduced network integration and segregation. Findings suggest that basic individual differences impact core metrics of brain function that are used in standard case-control studies. Considering individual variability and diversity in global settings would contribute to a more tailored understanding of brain function.

A 10-Year Longitudinal Study of Brain Cortical Thickness in People with First-Episode Psychosis Using Normative Models

BACKGROUND

Clinical forecasting models have potential to optimize treatment and improve outcomes in psychosis, but predicting long-term outcomes is challenging and long-term follow-up data are scarce. In this 10-year longitudinal study, we aimed to characterize the temporal evolution of cortical correlates of psychosis and their associations with symptoms.

DESIGN

Structural magnetic resonance imaging (MRI) from people with first-episode psychosis and controls (n = 79 and 218) were obtained at enrollment, after 12 months (n = 67 and 197), and 10 years (n = 23 and 77), within the Thematically Organized Psychosis (TOP) study. Normative models for cortical thickness estimated on public MRI datasets (n = 42 983) were applied to TOP data to obtain deviation scores for each region and timepoint. Positive and Negative Syndrome Scale (PANSS) scores were acquired at each timepoint along with registry data. Linear mixed effects models assessed effects of diagnosis, time, and their interactions on cortical deviations plus associations with symptoms.

RESULTS

LMEs revealed conditional main effects of diagnosis and time × diagnosis interactions in a distributed cortical network, where negative deviations in patients attenuate over time. In patients, symptoms also attenuate over time. LMEs revealed effects of anterior cingulate on PANSS total, and insular and orbitofrontal regions on PANSS negative scores.

CONCLUSIONS

This long-term longitudinal study revealed a distributed pattern of cortical differences which attenuated over time together with a reduction in symptoms. These findings are not in line with a simple neurodegenerative account of schizophrenia, and deviations from normative models offer a promising avenue to develop biomarkers to track clinical trajectories over time.

Delineating the Heterogeneity of Alzheimer's Disease and Mild Cognitive Impairment Using Normative Models of Dynamic Brain Functional Networks

BACKGROUND

Alzheimer's disease (AD), which has been identified as the most common type of dementia, presents considerable heterogeneity in its clinical manifestations. Early intervention at the stage of mild cognitive impairment (MCI) holds potential in AD prevention. However, characterizing the heterogeneity of neurobiological abnormalities and identifying MCI subtypes pose significant challenges.

METHODS

We constructed sex-specific normative age models of dynamic brain functional networks and mapped the deviations of the brain characteristics for individuals from multiple datasets, including 295 patients with AD, 441 patients with MCI, and 1160 normal control participants. Then, based on these individual deviation patterns, subtypes for both AD and MCI were identified using the clustering method, and their similarities and differences were comprehensively assessed.

RESULTS

Individuals with AD and MCI were clustered into 2 subtypes, and these subtypes exhibited significant differences in their intrinsic brain functional phenotypes and spatial atrophy patterns, as well as in disease progression and cognitive decline trajectories. The subtypes with positive deviations in AD and MCI shared similar deviation patterns, as did those with negative deviations. There was a potential transformation of MCI with negative deviation patterns into AD, and participants with MCI had a more severe cognitive decline rate.

CONCLUSIONS

In this study, we quantified neurophysiological heterogeneity by analyzing deviation patterns from the dynamic functional connectome normative model and identified disease subtypes of AD and MCI using a comprehensive resting-state functional magnetic resonance imaging multicenter dataset. The findings provide new insights for developing early prevention and personalized treatment strategies for AD.

Normative modelling of molecular-based functional circuits captures clinical heterogeneity transdiagnostically in psychiatric patients

Advanced methods such as REACT have allowed the integration of fMRI with the brain's receptor landscape, providing novel insights transcending the multiscale organisation of the brain. Similarly, normative modelling has allowed translational neuroscience to move beyond group-average differences and characterise deviations from health at an individual level. Here, we bring these methods together for the first time. We used REACT to create functional networks enriched with the main modulatory, inhibitory, and excitatory neurotransmitter systems and generated normative models of these networks to capture functional connectivity deviations in patients with schizophrenia, bipolar disorder (BPD), and ADHD. Substantial overlap was seen in symptomatology and deviations from normality across groups, but these could be mapped into a common space linking constellations of symptoms through to underlying neurobiology transdiagnostically. This work provides impetus for developing novel biomarkers that characterise molecular- and systems-level dysfunction at the individual level, facilitating the transition towards mechanistically targeted treatments.

Regional patterns of human cortex development colocalize with underlying neurobiology

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

Disentangling the Individual-Shared and Individual-Specific Subspace of Altered Brain Functional Connectivity in Autism Spectrum Disorder

BACKGROUND

Despite considerable effort toward understanding the neural basis of autism spectrum disorder (ASD) using case-control analyses of resting-state functional magnetic resonance imaging data, findings are often not reproducible, largely due to biological and clinical heterogeneity among individuals with ASD. Thus, exploring the individual-shared and individual-specific altered functional connectivity (AFC) in ASD is important to understand this complex, heterogeneous disorder.

METHODS

We considered 254 individuals with ASD and 295 typically developing individuals from the Autism Brain Imaging Data Exchange to explore the individual-shared and individual-specific subspaces of AFC. First, we computed AFC matrices of individuals with ASD compared with typically developing individuals. Then, common orthogonal basis extraction was used to project AFC of ASD onto 2 subspaces: an individual-shared subspace, which represents altered connectivity patterns shared across ASD, and an individual-specific subspace, which represents the remaining individual characteristics after eliminating the individual-shared altered connectivity patterns.

RESULTS

Analysis yielded 3 common components spanning the individual-shared subspace. Common components were associated with differences of functional connectivity at the group level. AFC in the individual-specific subspace improved the prediction of clinical symptoms. The default mode network-related and cingulo-opercular network-related magnitudes of AFC in the individual-specific subspace were significantly correlated with symptom severity in social communication deficits and restricted, repetitive behaviors in ASD.

CONCLUSIONS

Our study decomposed AFC of ASD into individual-shared and individual-specific subspaces, highlighting the importance of capturing and capitalizing on individual-specific brain connectivity features for dissecting heterogeneity. Our analysis framework provides a blueprint for parsing heterogeneity in other prevalent neurodevelopmental conditions.

A Test of the Behavioral Model of Tic Disorders Using a Dynamical Systems Framework

Tic disorders are a class of neurodevelopmental disorders characterized by involuntary motor and/or vocal tics. It has been hypothesized that tics function to reduce aversive premonitory urges (i.e., negative reinforcement) and that suppression-based behavioral interventions such as habit reversal training (HRT) and exposure and response prevention (ERP) disrupt this process and facilitate urge reduction through habituation. However, previous findings regarding the negative reinforcement hypothesis and the effect of suppression on the urge-tic relationship have been inconsistent. The present study applied a dynamical systems framework and within-subject time-series autoregressive models to examine the temporal dynamics of urges and tics and assess whether their relationship changes over time. Eleven adults with tic disorders provided continuous urge ratings during separate conditions in which they were instructed to tic freely or to suppress tics. During the free-to-tic conditions, there was considerable heterogeneity across participants in whether and how the urge-tic relationship followed a pattern consistent with the automatic negative reinforcement hypothesis. Further, little evidence for within-session habituation was seen; tic suppression did not result in a reduction in premonitory urges for most participants. Analysis of broader urge change metrics did show significant disruption to the urge pattern during suppression, which has implications for the current biobehavioral model of tics.

Resting-State Changes in Aging and Parkinson's Disease Are Shaped by Underlying Neurotransmission: A Normative Modeling Study

BACKGROUND

Human healthy and pathological aging is linked to a steady decline in brain resting-state activity and connectivity measures. The neurophysiological mechanisms that underlie these changes remain poorly understood.

METHODS

Making use of recent developments in normative modeling and availability of in vivo maps for various neurochemical systems, we tested in the UK Biobank cohort (n = 25,917) whether and how age- and Parkinson's disease-related resting-state changes in commonly applied local and global activity and connectivity measures colocalize with underlying neurotransmitter systems.

RESULTS

We found that the distributions of several major neurotransmitter systems including serotonergic, dopaminergic, noradrenergic, and glutamatergic neurotransmission correlated with age-related changes across functional activity and connectivity measures. Colocalization patterns in Parkinson's disease deviated from normative aging trajectories for these, as well as for cholinergic and GABAergic (gamma-aminobutyric acidergic) neurotransmission. The deviation from normal colocalization of brain function and GABAA correlated with disease duration.

CONCLUSIONS

These findings provide new insights into molecular mechanisms underlying age- and Parkinson's-related brain functional changes by extending the existing evidence elucidating the vulnerability of specific neurochemical attributes to normal aging and Parkinson's disease. The results particularly indicate that alongside dopamine and serotonin, increased vulnerability of glutamatergic, cholinergic, and GABAergic systems may also contribute to Parkinson's disease-related functional alterations. Combining normative modeling and neurotransmitter mapping may aid future research and drug development through deeper understanding of neurophysiological mechanisms that underlie specific clinical conditions.

Cortical thickness modeling and variability in Alzheimer's disease and frontotemporal dementia

BACKGROUND AND OBJECTIVE

Alzheimer's disease (AD) and frontotemporal dementia (FTD) show different patterns of cortical thickness (CTh) loss compared with healthy controls (HC), even though there is relevant heterogeneity between individuals suffering from each of these diseases. Thus, we developed CTh models to study individual variability in AD, FTD, and HC.

METHODS

We used the baseline CTh measures of 379 participants obtained from the structural MRI processed with FreeSurfer. A total of 169 AD patients (63 ± 9 years, 65 men), 88 FTD patients (64 ± 9 years, 43 men), and 122 HC (62 ± 10 years, 47 men) were studied. We fitted region-wise temporal models of CTh using Support Vector Regression. Then, we studied associations of individual deviations from the model with cerebrospinal fluid levels of neurofilament light chain (NfL) and 14-3-3 protein and Mini-Mental State Examination (MMSE). Furthermore, we used real longitudinal data from 144 participants to test model predictivity.

RESULTS

We defined CTh spatiotemporal models for each group with a reliable fit. Individual deviation correlated with MMSE for AD and with NfL for FTD. AD patients with higher deviations from the trend presented higher MMSE values. In FTD, lower NfL levels were associated with higher deviations from the CTh prediction. For AD and HC, we could predict longitudinal visits with the presented model trained with baseline data. For FTD, the longitudinal visits had more variability.

CONCLUSION

We highlight the value of CTh models for studying AD and FTD longitudinal changes and variability and their relationships with cognitive features and biomarkers.

Normative modelling of brain morphometry across the lifespan with CentileBrain: algorithm benchmarking and model optimisation

The value of normative models in research and clinical practice relies on their robustness and a systematic comparison of different modelling algorithms and parameters; however, this has not been done to date. We aimed to identify the optimal approach for normative modelling of brain morphometric data through systematic empirical benchmarking, by quantifying the accuracy of different algorithms and identifying parameters that optimised model performance. We developed this framework with regional morphometric data from 37 407 healthy individuals (53% female and 47% male; aged 3-90 years) from 87 datasets from Europe, Australia, the USA, South Africa, and east Asia following a comparative evaluation of eight algorithms and multiple covariate combinations pertaining to image acquisition and quality, parcellation software versions, global neuroimaging measures, and longitudinal stability. The multivariate fractional polynomial regression (MFPR) emerged as the preferred algorithm, optimised with non-linear polynomials for age and linear effects of global measures as covariates. The MFPR models showed excellent accuracy across the lifespan and within distinct age-bins and longitudinal stability over a 2-year period. The performance of all MFPR models plateaued at sample sizes exceeding 3000 study participants. This model can inform about the biological and behavioural implications of deviations from typical age-related neuroanatomical changes and support future study designs. The model and scripts described here are freely available through CentileBrain.

Understanding divergence: Placing developmental neuroscience in its dynamic context

Neurodevelopment is not merely a process of brain maturation, but an adaptation to constraints unique to each individual and to the environments we co-create. However, our theoretical and methodological toolkits often ignore this reality. There is growing awareness that a shift is needed that allows us to study divergence of brain and behaviour across conventional categorical boundaries. However, we argue that in future our study of divergence must also incorporate the developmental dynamics that capture the emergence of those neurodevelopmental differences. This crucial step will require adjustments in study design and methodology. If our ultimate aim is to incorporate the developmental dynamics that capture how, and ultimately when, divergence takes place then we will need an analytic toolkit equal to these ambitions. We argue that the over reliance on group averages has been a conceptual dead-end with regard to the neurodevelopmental differences. This is in part because any individual differences and developmental dynamics are inevitably lost within the group average. Instead, analytic approaches which are themselves new, or simply newly applied within this context, may allow us to shift our theoretical and methodological frameworks from groups to individuals. Likewise, methods capable of modelling complex dynamic systems may allow us to understand the emergent dynamics only possible at the level of an interacting neural system.

Using normative modeling to assess pharmacological treatment effect on brain state in patients with schizophrenia

Quantifying individual differences in neuroimaging metrics is attracting interest in clinical studies with mental disorders. Schizophrenia is diagnosed exclusively based on symptoms, and the biological heterogeneity makes it difficult to accurately assess pharmacological treatment effects on the brain state. Using the Cambridge Centre for Ageing and Neuroscience data set, we built normative models of brain states and mapped the deviations of the brain characteristics of each patient, to test whether deviations were related to symptoms, and further investigated the pharmacological treatment effect on deviation distributions. Specifically, we found that the patients can be divided into 2 groups: the normalized group had a normalization trend and milder symptoms at baseline, and the other group showed a more severe deviation trend. The baseline severity of the depression as well as the overall symptoms could predict the deviation of the static characteristics for the dorsal and ventral attention networks after treatment. In contrast, the positive symptoms could predict the deviations of the dynamic fluctuations for the default mode and dorsal attention networks after treatment. This work evaluates the effect of pharmacological treatment on static and dynamic brain states using an individualized approach, which may assist in understanding the heterogeneity of the illness pathology as well as the treatment response.

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.

Mapping Neurophysiological Subtypes of Major Depressive Disorder Using Normative Models of the Functional Connectome

BACKGROUND

Major depressive disorder (MDD) is a highly heterogeneous disorder that typically emerges in adolescence and can occur throughout adulthood. Studies aimed at quantitatively uncovering the heterogeneity of individual functional connectome abnormalities in MDD and identifying reproducibly distinct neurophysiological MDD subtypes across the lifespan, which could provide promising insights for precise diagnosis and treatment prediction, are still lacking.

METHODS

Leveraging resting-state functional magnetic resonance imaging data from 1148 patients with MDD and 1079 healthy control participants (ages 11-93), we conducted the largest multisite analysis to date for neurophysiological MDD subtyping. First, we characterized typical lifespan trajectories of functional connectivity strength based on the normative model and quantitatively mapped the heterogeneous individual deviations among patients with MDD. Then, we identified neurobiological MDD subtypes using an unsupervised clustering algorithm and evaluated intersite reproducibility. Finally, we validated the subtype differences in baseline clinical variables and longitudinal treatment predictive capacity.

RESULTS

Our findings indicated great intersubject heterogeneity in the spatial distribution and severity of functional connectome deviations among patients with MDD, which inspired the identification of 2 reproducible neurophysiological subtypes. Subtype 1 showed severe deviations, with positive deviations in the default mode, limbic, and subcortical areas and negative deviations in the sensorimotor and attention areas. Subtype 2 showed a moderate but converse deviation pattern. More importantly, subtype differences were observed in depressive item scores and the predictive ability of baseline deviations for antidepressant treatment outcomes.

CONCLUSIONS

These findings shed light on our understanding of different neurobiological mechanisms underlying the clinical heterogeneity of MDD and are essential for developing personalized treatments for this disorder.

Brain-charting autism and attention deficit hyperactivity disorder reveals distinct and overlapping neurobiology

Background

Autism and attention deficit hyperactivity disorder (ADHD) are heterogeneous neurodevelopmental conditions with complex underlying neurobiology. Despite overlapping presentation and sex-biased prevalence, autism and ADHD are rarely studied together, and sex differences are often overlooked. Normative modelling provides a unified framework for studying age-specific and sex-specific divergences in neurodivergent brain development.

Methods

Here we use normative modelling and a large, multi-site neuroimaging dataset to characterise cortical anatomy associated with autism and ADHD, benchmarked against models of typical brain development based on a sample of over 75,000 individuals. We also examined sex and age differences, relationship with autistic traits, and explored the co-occurrence of autism and ADHD (autism+ADHD).

Results

We observed robust neuroanatomical signatures of both autism and ADHD. Overall, autistic individuals showed greater cortical thickness and volume localised to the superior temporal cortex, whereas individuals with ADHD showed more global effects of cortical thickness increases but lower cortical volume and surface area across much of the cortex. The autism+ADHD group displayed a unique pattern of widespread increases in cortical thickness, and certain decreases in surface area. We also found evidence that sex modulates the neuroanatomy of autism but not ADHD, and an age-by-diagnosis interaction for ADHD only.

Conclusions

These results indicate distinct cortical differences in autism and ADHD that are differentially impacted by age, sex, and potentially unique patterns related to their co-occurrence.

Normative Modeling of Brain Morphometry Across the Lifespan Using CentileBrain: Algorithm Benchmarking and Model Optimization

We present an empirically benchmarked framework for sex-specific normative modeling of brain morphometry that can inform about the biological and behavioral significance of deviations from typical age-related neuroanatomical changes and support future study designs. This framework was developed using regional morphometric data from 37,407 healthy individuals (53% female; aged 3-90 years) following a comparative evaluation of eight algorithms and multiple covariate combinations pertaining to image acquisition and quality, parcellation software versions, global neuroimaging measures, and longitudinal stability. The Multivariate Factorial Polynomial Regression (MFPR) emerged as the preferred algorithm optimized using nonlinear polynomials for age and linear effects of global measures as covariates. The MFPR models showed excellent accuracy across the lifespan and within distinct age-bins, and longitudinal stability over a 2-year period. The performance of all MFPR models plateaued at sample sizes exceeding 3,000 study participants. The model and scripts described here are freely available through CentileBrain (https://centilebrain.org/).

Predictive Clinical Neuroscience Portal (PCNportal): instant online access to research-grade normative models for clinical neuroscientists

Background

The neurobiology of mental disorders remains poorly understood despite substantial scientific efforts, due to large clinical heterogeneity and to a lack of tools suitable to map individual variability. Normative modeling is one recently successful framework that can address these problems by comparing individuals to a reference population. The methodological underpinnings of normative modelling are, however, relatively complex and computationally expensive. Our research group has developed the python-based normative modelling package Predictive Clinical Neuroscience toolkit (PCNtoolkit) which provides access to many validated algorithms for normative modelling. PCNtoolkit has since proven to be a strong foundation for large scale normative modelling, but still requires significant computation power, time and technical expertise to develop.

Methods

To address these problems, we introduce PCNportal. PCNportal is an online platform integrated with PCNtoolkit that offers access to pre-trained research-grade normative models estimated on tens of thousands of participants, without the need for computation power or programming abilities. PCNportal is an easy-to-use web interface that is highly scalable to large user bases as necessary. Finally, we demonstrate how the resulting normalized deviation scores can be used in a clinical application through a schizophrenia classification task applied to cortical thickness and volumetric data from the longitudinal Northwestern University Schizophrenia Data and Software Tool (NUSDAST) dataset.

Results

At each longitudinal timepoint, the transferred normative models achieved a mean[std. dev.] explained variance of 9.4[8.8]%, 9.2[9.2]%, 5.6[7.4]% respectively in the control group and 4.7[5.5]%, 6.0[6.2]%, 4.2[6.9]% in the schizophrenia group. Diagnostic classifiers achieved AUC of 0.78, 0.76 and 0.71 respectively.

Conclusions

This replicates the utility of normative models for diagnostic classification of schizophrenia and showcases the use of PCNportal for clinical neuroimaging. By facilitating and speeding up research with high-quality normative models, this work contributes to research in inter-individual variability, clinical heterogeneity and precision medicine.

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.

Assessing attention and impulsivity in the variable stimulus duration and variable intertrial interval rodent continuous performance test schedules using dopamine receptor antagonists in female C57BL/6JRj mice

RATIONALE

Dopaminergic dysfunction is implicated in disorders of impulsivity and inattention. The rodent continuous performance test (rCPT) has been used to quantify changes in attention and impulsivity.

OBJECTIVE

To examine the roles of dopamine receptors in attention and impulsivity behaviours measured in the rCPT variable stimulus duration (vSD) and the variable intertrial interval schedules (vITI) using DA receptor antagonists.

METHODS

Two cohorts of 35 and 36 female C57BL/6JRj mice were examined separately in the rCPT, vSD, and vITI schedules, respectively. Both cohorts received antagonists of the following receptors: D1/5 (SCH23390, SCH: 0.01, 0.02, 0.04 mg/kg) and D2/3 (raclopride, RAC 0.03, 0.10, 0.30 mg/kg) in consecutive balanced Latin square designs with flanking reference measurements. The antagonists were subsequently examined for effects on locomotor activity.

RESULTS

SCH showed similar effects in both schedules, and the effects were reference-dependent in the vITI schedule. SCH reduced responding, but improved response accuracy, impulsivity, discriminability, and locomotor activity. RAC showed mixed effects on responsivity, but improved accuracy and discriminability. The discriminability improvement was driven by an increase in hit rate in the vITI schedule and a reduction in false alarm rate in the vSD schedule. RAC also decreased locomotor activity.

CONCLUSION

Both D1/5 and D2/3 receptor antagonism reduced responding, but the outcome on discriminability differed, stemming from individual effects on hit and false alarm rate, and the weight of omissions within the calculation. The effects of SCH and RAC suggest that endogenous DA increases responding and impulsivity, but reduces accuracy and shows mixed effects on discriminability.

Assessing attention and impulsivity in the variable stimulus duration and variable intertrial interval rodent continuous performance test schedules using noradrenaline receptor antagonists in female C57BL/6JRj mice

RATIONALE

Noradrenergic dysfunction is associated with disorders of impulsivity and inattention. The rodent continuous performance test (rCPT) quantifies changes in attention and impulsivity.

OBJECTIVE

To use NA receptor antagonists to examine the roles of NA on attention and impulsivity behaviours measured in the rCPT variable stimulus duration (vSD) and the variable intertrial interval (vITI) schedules.

METHODS

Two cohorts of 36 female C57BL/6JRj mice were examined separately in the rCPT vSD and vITI schedules. Both cohorts received antagonists of the following adrenoceptors: α1 (doxazosin, DOX: 1.0, 3.0, 10.0 mg/kg), α2 (yohimbine, YOH: 0.1, 0.3, 1.0 mg/kg), and β1/2 (propranolol, PRO: 1.0, 3.0, 10.0 mg/kg) in consecutive balanced Latin square designs with flanking reference measurements. The antagonists were subsequently examined for effects on locomotor activity.

RESULTS

DOX showed similar effects in both schedules, improving discriminability and accuracy, and reducing responding and impulsivity, and DOX also reduced locomotor activity. YOH showed prominent effects in the vSD schedule to increase responding and impulsivity, while impairing discriminability and accuracy. YOH did not affect locomotor activity. PRO increased responding and impulsivity, decreased accuracy, but did not affect discriminability or locomotor activity.

CONCLUSION

Antagonism of α2 or β1/2 adrenoceptors caused similar increases in responding and impulsivity and worsened attentional performance, while α1 adrenoceptor antagonism showed the opposite effects. Our results suggest that endogenous NA exerts bidirectional control of most behaviours in the rCPT. The parallel vSD and vITI studies showed a substantial overlap in effects, but also some differences that indicate differing sensitivity towards noradrenergic manipulations.

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.

Parsing altered gray matter morphology of depression using a framework integrating the normative model and non-negative matrix factorization

The high inter-individual heterogeneity in individuals with depression limits neuroimaging studies with case-control approaches to identify promising biomarkers for individualized clinical decision-making. We put forward a framework integrating the normative model and non-negative matrix factorization (NMF) to quantitatively assess altered gray matter morphology in depression from a dimensional perspective. The proposed framework parses altered gray matter morphology into overlapping latent disease factors, and assigns patients distinct factor compositions, thus preserving inter-individual variability. We identified four robust disease factors with distinct clinical symptoms and cognitive processes in depression. In addition, we showed the quantitative relationship between the group-level gray matter morphological differences and disease factors. Furthermore, this framework significantly predicted factor compositions of patients in an independent dataset. The framework provides an approach to resolve neuroanatomical heterogeneity in depression.

Mapping the neuroanatomical heterogeneity of OCD using a framework integrating normative model and non-negative matrix factorization

Obsessive-compulsive disorder (OCD) is a spectrum disorder with high interindividual heterogeneity. We propose a comprehensible framework integrating normative model and non-negative matrix factorization (NMF) to quantitatively estimate the neuroanatomical heterogeneity of OCD from a dimensional perspective. T1-weighted magnetic resonance images of 98 first-episode untreated patients with OCD and matched healthy controls (HCs, n = 130) were acquired. We derived individualized differences in gray matter morphometry using normative model and parsed them into latent disease factors using NMF. Four robust disease factors were identified. Each patient expressed multiple factors and exhibited a unique factor composition. Factor compositions of patients were significantly correlated with severity of symptom, age of onset, illness duration, and exhibited sex differences, capturing sources of clinical heterogeneity. In addition, the group-level morphological differences obtained with two-sample t test could be quantitatively derived from the identified disease factors, reconciling the group-level and subject-level findings in neuroimaging studies. Finally, we uncovered two distinct subtypes with opposite morphological differences compared with HCs from factor compositions. Our findings suggest that morphological differences of individuals with OCD are the unique combination of distinct neuroanatomical patterns. The proposed framework quantitatively estimating neuroanatomical heterogeneity paves the way for precision medicine in OCD.

Relationship Between Prediction Accuracy and Feature Importance Reliability: an Empirical and Theoretical Study

There is significant interest in using neuroimaging data to predict behavior. The predictive models are often interpreted by the computation of feature importance, which quantifies the predictive relevance of an imaging feature. Tian and Zalesky (2021) suggest that feature importance estimates exhibit low split-half reliability, as well as a trade-off between prediction accuracy and feature importance reliability across parcellation resolutions. However, it is unclear whether the trade-off between prediction accuracy and feature importance reliability is universal. Here, we demonstrate that, with a sufficient sample size, feature importance (operationalized as Haufe-transformed weights) can achieve fair to excellent split-half reliability. With a sample size of 2600 participants, Haufe-transformed weights achieve average intra-class correlation coefficients of 0.75, 0.57 and 0.53 for cognitive, personality and mental health measures respectively. Haufe-transformed weights are much more reliable than original regression weights and univariate FC-behavior correlations. Original regression weights are not reliable even with 2600 participants. Intriguingly, feature importance reliability is strongly positively correlated with prediction accuracy across phenotypes. Within a particular behavioral domain, there is no clear relationship between prediction performance and feature importance reliability across regression models. Furthermore, we show mathematically that feature importance reliability is necessary, but not sufficient, for low feature importance error. In the case of linear models, lower feature importance error is mathematically related to lower prediction error. Therefore, higher feature importance reliability might yield lower feature importance error and higher prediction accuracy. Finally, we discuss how our theoretical results relate with the reliability of imaging features and behavioral measures. Overall, the current study provides empirical and theoretical insights into the relationship between prediction accuracy and feature importance reliability.

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.

Normative models for neuroimaging markers: Impact of model selection, sample size and evaluation criteria

Modelling population reference curves or normative modelling is increasingly used with the advent of large neuroimaging studies. In this paper we assess the performance of fitting methods from the perspective of clinical applications and investigate the influence of the sample size. Further, we evaluate linear and non-linear models for percentile curve estimation and highlight how the bias-variance trade-off manifests in typical neuroimaging data. We created plausible ground truth distributions of hippocampal volumes in the age range of 45 to 80 years, as an example application. Based on these distributions we repeatedly simulated samples for sizes between 50 and 50,000 data points, and for each simulated sample we fitted a range of normative models. We compared the fitted models and their variability across repetitions to the ground truth, with specific focus on the outer percentiles (1st, 5th, 10th) as these are the most clinically relevant. Our results quantify the expected decreasing trend in variance of the volume estimates with increasing sample size. However, bias in the volume estimates only decreases a modest amount, without much improvement at large sample sizes. The uncertainty of model performance is substantial for what would often be considered large samples in a neuroimaging context and rises dramatically at the ends of the age range, where fewer data points exist. Flexible models perform better across sample sizes, especially for non-linear ground truth. Surprisingly large samples of several thousand data points are needed to accurately capture outlying percentiles across the age range for applications in research and clinical settings. Performance evaluation methods should assess both bias and variance. Furthermore, caution is needed when attempting to go near the ends of the age range captured by the source data set and, as is a well known general principle, extrapolation beyond the age range should always be avoided. To help with such evaluations of normative models we have made our code available to guide researchers developing or utilising normative models.

Evidence from "big data" for the default-mode hypothesis of ADHD: a mega-analysis of multiple large samples

We sought to identify resting-state characteristics related to attention deficit/hyperactivity disorder, both as a categorical diagnosis and as a trait feature, using large-scale samples which were processed according to a standardized pipeline. In categorical analyses, we considered 1301 subjects with diagnosed ADHD, contrasted against 1301 unaffected controls (total N = 2602; 1710 males (65.72%); mean age = 10.86 years, sd = 2.05). Cases and controls were 1:1 nearest neighbor matched on in-scanner motion and key demographic variables and drawn from multiple large cohorts. Associations between ADHD-traits and resting-state connectivity were also assessed in a large multi-cohort sample (N = 10,113). ADHD diagnosis was associated with less anticorrelation between the default mode and salience/ventral attention (B = 0.009, t = 3.45, p-FDR = 0.004, d = 0.14, 95% CI = 0.004, 0.014), somatomotor (B = 0.008, t = 3.49, p-FDR = 0.004, d = 0.14, 95% CI = 0.004, 0.013), and dorsal attention networks (B = 0.01, t = 4.28, p-FDR < 0.001, d = 0.17, 95% CI = 0.006, 0.015). These results were robust to sensitivity analyses considering comorbid internalizing problems, externalizing problems and psychostimulant medication. Similar findings were observed when examining ADHD traits, with the largest effect size observed for connectivity between the default mode network and the dorsal attention network (B = 0.0006, t = 5.57, p-FDR < 0.001, partial-r = 0.06, 95% CI = 0.0004, 0.0008). We report significant ADHD-related differences in interactions between the default mode network and task-positive networks, in line with default mode interference models of ADHD. Effect sizes (Cohen's d and partial-r, estimated from the mega-analytic models) were small, indicating subtle group differences. The overlap between the affected brain networks in the clinical and general population samples supports the notion of brain phenotypes operating along an ADHD continuum.

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.

Different functional alteration in attention-deficit/hyperactivity disorder across developmental age groups: A meta-analysis and an independent validation of resting-state functional connectivity studies

BACKGROUND

Attention-deficit/hyperactivity disorder (ADHD) is a highly complex and heterogeneous disorder. Abnormal brain connectivity in ADHD might be influenced by developmental ages which might lead to the lacking of significant spatial convergence across studies. However, the developmental patterns and mechanisms of ADHD brain connectivity remain to be fully uncovered.

METHODS

In the present study, we searched PubMed, Scopus, Web of Science, and Embase for seed-based whole-brain resting-state functional connectivity studies of ADHD published through October 12th, 2020. The seeds meeting inclusion criteria were categorized into the cortex group and subcortex group, as previous studies suggested that the cortex and subcortex have different temporal patterns of development. Activation likelihood estimation meta-analysis was performed to investigate the abnormal connectivity in different age groups (all-age group, younger: <12 years, older: ≥12 years). Moreover, significant convergence of reported foci was used as seeds for validation with our independent dataset.

RESULTS

As with previous studies, scarce results were found in the all-age group. However, we found that the younger group consistently exhibited hyper-connectivity between different parts of the cortex and left middle frontal gyrus, and hypo-connectivity between different parts of the cortex and left putamen/pallidus/amygdala. Whereas, the older group (mainly for adults) showed hyper-connectivity between the cortex and right precuneus/sub-gyral/cingulate gyrus. Besides, the abnormal cortico-cortical and cortico-subcortical functional connectivity in children, and the abnormal cortico-cortical functional connectivity in adults were verified in our independent dataset.

CONCLUSION

Our study emphasizes the importance of developmental age effects on the study of brain networks in ADHD. Further, we proposed that cortico-cortical and cortico-subcortical connectivity might play an important role in the pathophysiology of children with ADHD, while abnormal cortico-cortical connections were more important for adults with ADHD. This work provided a potential new insight to understand the neurodevelopmental mechanisms and possible clinical application of ADHD.

Impact of neovascular age-related macular degeneration on eye-movement control during scene viewing: Viewing biases and guidance by visual salience

Human vision requires us to analyze the visual periphery to decide where to fixate next. In the present study, we investigated this process in people with age-related macular degeneration (AMD). In particular, we examined viewing biases and the extent to which visual salience guides fixation selection during free-viewing of naturalistic scenes. We used an approach combining generalized linear mixed modeling (GLMM) with a-priori scene parcellation. This method allows one to investigate group differences in terms of scene coverage and observers' well-known tendency to look at the center of scene images. Moreover, it allows for testing whether image salience influences fixation probability above and beyond what can be accounted for by the central bias. Compared with age-matched normally sighted control subjects (and young subjects), AMD patients' viewing behavior was less exploratory, with a stronger central fixation bias. All three subject groups showed a salience effect on fixation selection-higher-salience scene patches were more likely to be fixated. Importantly, the salience effect for the AMD group was of similar size as the salience effect for the control group, suggesting that guidance by visual salience was still intact. The variances for by-subject random effects in the GLMM indicated substantial individual differences. A separate model exclusively considered the AMD data and included fixation stability as a covariate, with the results suggesting that reduced fixation stability was associated with a reduced impact of visual salience on fixation selection.

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.

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.

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.

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.

Site effects how-to and when: An overview of retrospective techniques to accommodate site effects in multi-site neuroimaging analyses

Site differences, or systematic differences in feature distributions across multiple data-acquisition sites, are a known source of heterogeneity that may adversely affect large-scale meta- and mega-analyses of independently collected neuroimaging data. They influence nearly all multi-site imaging modalities and biomarkers, and methods to compensate for them can improve reliability and generalizability in the analysis of genetics, omics, and clinical data. The origins of statistical site effects are complex and involve both technical differences (scanner vendor, head coil, acquisition parameters, imaging processing) and differences in sample characteristics (inclusion/exclusion criteria, sample size, ancestry) between sites. In an age of expanding international consortium research, there is a growing need to disentangle technical site effects from sample characteristics of interest. Numerous statistical and machine learning methods have been developed to control for, model, or attenuate site effects - yet to date, no comprehensive review has discussed the benefits and drawbacks of each for different use cases. Here, we provide an overview of the different existing statistical and machine learning methods developed to remove unwanted site effects from independently collected neuroimaging samples. We focus on linear mixed effect models, the ComBat technique and its variants, adjustments based on image quality metrics, normative modeling, and deep learning approaches such as generative adversarial networks. For each method, we outline the statistical foundation and summarize strengths and weaknesses, including their assumptions and conditions of use. We provide information on software availability and comment on the ease of use and the applicability of these methods to different types of data. We discuss validation and comparative reports, mention caveats and provide guidance on when to use each method, depending on context and specific research questions.

Cofluctuation analysis reveals aberrant default mode network patterns in adolescents and youths with autism spectrum disorder

Resting-state functional connectivity (rsFC) approaches provide informative estimates of the functional architecture of the brain, and recently-proposed cofluctuation analysis temporally unwraps FC at every moment in time, providing refined information for quantifying brain dynamics. As a brain network disorder, autism spectrum disorder (ASD) was characterized by substantial alteration in FC, but the contribution of moment-to-moment-activity cofluctuations to the overall dysfunctional connectivity pattern in ASD remains poorly understood. Here, we used the cofluctuation approach to explore the underlying dynamic properties of FC in ASD, using a large multisite resting-state functional magnetic resonance imaging (rs-fMRI) dataset (ASD = 354, typically developing controls [TD] = 446). Our results verified that the networks estimated using high-amplitude frames were highly correlated with the traditional rsFC. Moreover, these frames showed higher average amplitudes in participants with ASD than those in the TD group. Principal component analysis was performed on the activity patterns in these frames and aggregated over all subjects. The first principal component (PC1) corresponds to the default mode network (DMN), and the PC1 coefficients were greater in participants with ASD than those in the TD group. Additionally, increased ASD symptom severity was associated with the increased coefficients, which may result in excessive internally oriented cognition and social cognition deficits in individuals with ASD. Our finding highlights the utility of cofluctuation approaches in prevalent neurodevelopmental disorders and verifies that the aberrant contribution of DMN to rsFC may underline the symptomatology in adolescents and youths with ASD.

Identification of shared and differentiating genetic architecture for autism spectrum disorder, attention-deficit hyperactivity disorder and case subgroups

Attention-deficit hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) are highly heritable neurodevelopmental conditions, with considerable overlap in their genetic etiology. We dissected their shared and distinct genetic etiology by cross-disorder analyses of large datasets. We identified seven loci shared by the disorders and five loci differentiating them. All five differentiating loci showed opposite allelic directions in the two disorders and significant associations with other traits, including educational attainment, neuroticism and regional brain volume. Integration with brain transcriptome data enabled us to identify and prioritize several significantly associated genes. The shared genomic fraction contributing to both disorders was strongly correlated with other psychiatric phenotypes, whereas the differentiating portion was correlated most strongly with cognitive traits. Additional analyses revealed that individuals diagnosed with both ASD and ADHD were double-loaded with genetic predispositions for both disorders and showed distinctive patterns of genetic association with other traits compared with the ASD-only and ADHD-only subgroups. These results provide insights into the biological foundation of the development of one or both conditions and of the factors driving psychopathology discriminatively toward either ADHD or ASD.

The expanding horizons of network neuroscience: From description to prediction and control

The field of network neuroscience has emerged as a natural framework for the study of the brain and has been increasingly applied across divergent problems in neuroscience. From a disciplinary perspective, network neuroscience originally emerged as a formal integration of graph theory (from mathematics) and neuroscience (from biology). This early integration afforded marked utility in describing the interconnected nature of neural units, both structurally and functionally, and underscored the relevance of that interconnection for cognition and behavior. But since its inception, the field has not remained static in its methodological composition. Instead, it has grown to use increasingly advanced graph-theoretic tools and to bring in several other disciplinary perspectives-including machine learning and systems engineering-that have proven complementary. In doing so, the problem space amenable to the discipline has expanded markedly. In this review, we discuss three distinct flavors of investigation in state-of-the-art network neuroscience: (i) descriptive network neuroscience, (ii) predictive network neuroscience, and (iii) a perturbative network neuroscience that draws on recent advances in network control theory. In considering each area, we provide a brief summary of the approaches, discuss the nature of the insights obtained, and highlight future directions.

Nomograms of human hippocampal volume shifted by polygenic scores

Nomograms are important clinical tools applied widely in both developing and aging populations. They are generally constructed as normative models identifying cases as outliers to a distribution of healthy controls. Currently used normative models do not account for genetic heterogeneity. Hippocampal volume (HV) is a key endophenotype for many brain disorders. Here, we examine the impact of genetic adjustment on HV nomograms and the translational ability to detect dementia patients. Using imaging data from 35,686 healthy subjects aged 44-82 from the UK Biobank (UKB), we built HV nomograms using Gaussian process regression (GPR), which - compared to a previous method - extended the application age by 20 years, including dementia critical age ranges. Using HV polygenic scores (HV-PGS), we built genetically adjusted nomograms from participants stratified into the top and bottom 30% of HV-PGS. This shifted the nomograms in the expected directions by ~100 mm3 (2.3% of the average HV), which equates to 3 years of normal aging for a person aged ~65. Clinical impact of genetically adjusted nomograms was investigated by comparing 818 subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database diagnosed as either cognitively normal (CN), having mild cognitive impairment (MCI) or Alzheimer's disease (AD) patients. While no significant change in the survival analysis was found for MCI-to-AD conversion, an average of 68% relative decrease was found in intra-diagnostic-group variance, highlighting the importance of genetic adjustment in untangling phenotypic heterogeneity.

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.