Moving beyond the mean: Subgroups and dimensions of brain activity and cognitive performance across domains

Neuroimaging-based variability in subtyping biomarkers for psychiatric heterogeneity

Neuroimaging-based subtyping is increasingly used to explain heterogeneity in psychiatric disorders. However, the clinical utility of these subtyping efforts remains unclear, and replication has been challenging. Here we examined how the choice of neuroimaging measures influences the derivation of neuro-subtypes and the consequences for clinical delineation. On a clinically heterogeneous dataset (total n = 566) that included controls (n = 268) and cases (n = 298) of psychiatric conditions, including individuals diagnosed with post-traumatic stress disorder (PTSD), traumatic brain injury (TBI), and comorbidity of both (PTSD&TBI), we identified neuro-subtypes among the cases using either structural, resting-state, or task-based measures. The neuro-subtypes for each modality had high internal validity but did not significantly differ in their clinical and cognitive profiles. We further show that the choice of neuroimaging measures for subtyping substantially impacts the identification of neuro-subtypes, leading to low concordance across subtyping solutions. Similar variability in neuro-subtyping was found in an independent dataset (n = 1642) comprised of major depression disorder (MDD, n = 848) and controls (n = 794). Our results suggest that the highly anticipated relationships between neuro-subtypes and clinical features may be difficult to discover.

Functional magnetic resonance imaging in schizophrenia: current evidence, methodological advances, limitations and future directions

Functional neuroimaging emerged with great promise and has provided fundamental insights into the neurobiology of schizophrenia. However, it has faced challenges and criticisms, most notably a lack of clinical translation. This paper provides a comprehensive review and critical summary of the literature on functional neuroimaging, in particular functional magnetic resonance imaging (fMRI), in schizophrenia. We begin by reviewing research on fMRI biomarkers in schizophrenia and the clinical high risk phase through a historical lens, moving from case-control regional brain activation to global connectivity and advanced analytical approaches, and more recent machine learning algorithms to identify predictive neuroimaging features. Findings from fMRI studies of negative symptoms as well as of neurocognitive and social cognitive deficits are then reviewed. Functional neural markers of these symptoms and deficits may represent promising treatment targets in schizophrenia. Next, we summarize fMRI research related to antipsychotic medication, psychotherapy and psychosocial interventions, and neurostimulation, including treatment response and resistance, therapeutic mechanisms, and treatment targeting. We also review the utility of fMRI and data-driven approaches to dissect the heterogeneity of schizophrenia, moving beyond case-control comparisons, as well as methodological considerations and advances, including consortia and precision fMRI. Lastly, limitations and future directions of research in the field are discussed. Our comprehensive review suggests that, in order for fMRI to be clinically useful in the care of patients with schizophrenia, research should address potentially actionable clinical decisions that are routine in schizophrenia treatment, such as which antipsychotic should be prescribed or whether a given patient is likely to have persistent functional impairment. The potential clinical utility of fMRI is influenced by and must be weighed against cost and accessibility factors. Future evaluations of the utility of fMRI in prognostic and treatment response studies may consider including a health economics analysis.

Sex differences in brain-behavior relationships in the first two years of life

Background

Evidence for sex differences in cognition in childhood is established, but less is known about the underlying neural mechanisms for these differences. Recent findings suggest the existence of brain-behavior relationship heterogeneities during infancy; however, it remains unclear whether sex underlies these heterogeneities during this critical period when sex-related behavioral differences arise.

Methods

A sample of 316 infants was included with resting-state functional magnetic resonance imaging scans at neonate (3 weeks), 1, and 2 years of age. We used multiple linear regression to test interactions between sex and resting-state functional connectivity on behavioral scores of working memory, inhibitory self-control, intelligence, and anxiety collected at 4 years of age.

Results

We found six age-specific, intra-hemispheric connections showing significant and robust sex differences in functional connectivity-behavior relationships. All connections are either with the prefrontal cortex or the temporal pole, which has direct anatomical pathways to the prefrontal cortex. Sex differences in functional connectivity only emerge when associated with behavior, and not in functional connectivity alone. Furthermore, at neonate and 2 years of age, these age-specific connections displayed greater connectivity in males and lower connectivity in females in association with better behavioral scores.

Conclusions

Taken together, we critically capture robust and conserved brain mechanisms that are distinct to sex and are defined by their relationship to behavioral outcomes. Our results establish brain-behavior mechanisms as an important feature in the search for sex differences during development.

Data-driven multivariate identification of gyrification patterns in a transdiagnostic patient cohort: A cluster analysis approach

BACKGROUND

Multivariate data-driven statistical approaches offer the opportunity to study multi-dimensional interdependences between a large set of biological parameters, such as high-dimensional brain imaging data. For gyrification, a putative marker of early neurodevelopment, direct comparisons of patterns among multiple psychiatric disorders and investigations of potential heterogeneity of gyrification within one disorder and a transdiagnostic characterization of neuroanatomical features are lacking.

METHODS

In this study we used a data-driven, multivariate statistical approach to analyze cortical gyrification in a large cohort of N = 1028 patients with major psychiatric disorders (Major depressive disorder: n = 783, bipolar disorder: n = 129, schizoaffective disorder: n = 44, schizophrenia: n = 72) to identify cluster patterns of gyrification beyond diagnostic categories.

RESULTS

Cluster analysis applied on gyrification data of 68 brain regions (DK-40 atlas) identified three clusters showing difference in overall (global) gyrification and minor regional variation (regions). Newly, data-driven subgroups are further discriminative in cognition and transdiagnostic disease risk factors.

CONCLUSIONS

Results indicate that gyrification is associated with transdiagnostic risk factors rather than diagnostic categories and further imply a more global role of gyrification related to mental health than a disorder specific one. Our findings support previous studies highlighting the importance of association cortices involved in psychopathology. Explorative, data-driven approaches like ours can help to elucidate if the brain imaging data on hand and its a priori applied grouping actually has the potential to find meaningful effects or if previous hypotheses about the phenotype as well as its grouping have to be revisited.

Spatial-extent inference for testing variance components in reliability and heritability studies

Clusterwise inference is a popular approach in neuroimaging to increase sensitivity, but most existing methods are currently restricted to the General Linear Model (GLM) for testing mean parameters. Statistical methods for testing variance components, which are critical in neuroimaging studies that involve estimation of narrow-sense heritability or test-retest reliability, are underdeveloped due to methodological and computational challenges, which would potentially lead to low power. We propose a fast and powerful test for variance components called CLEAN-V (CLEAN for testing Variance components). CLEAN-V models the global spatial dependence structure of imaging data and computes a locally powerful variance component test statistic by data-adaptively pooling neighborhood information. Correction for multiple comparisons is achieved by permutations to control family-wise error rate (FWER). Through analysis of task-fMRI data from the Human Connectome Project across five tasks and comprehensive data-driven simulations, we show that CLEAN-V outperforms existing methods in detecting test-retest reliability and narrow-sense heritability with significantly improved power, with the detected areas aligning with activation maps. The computational efficiency of CLEAN-V also speaks of its practical utility, and it is available as an R package.

Stratification of responses to tDCS intervention in a healthy pediatric population based on resting-state EEG profiles

Transcranial Direct Current Stimulation (tDCS) is a non-invasive neuromodulation technique with a wide variety of clinical and research applications. As increasingly acknowledged, its effectiveness is subject dependent, which may lead to time consuming and cost ineffective treatment development phases. We propose the combination of electroencephalography (EEG) and unsupervised learning for the stratification and prediction of individual responses to tDCS. A randomized, sham-controlled, double-blind crossover study design was conducted within a clinical trial for the development of pediatric treatments based on tDCS. The tDCS stimulation (sham and active) was applied either in the left dorsolateral prefrontal cortex or in the right inferior frontal gyrus. Following the stimulation session, participants performed 3 cognitive tasks to assess the response to the intervention: the Flanker Task, N-Back Task and Continuous Performance Test (CPT). We used data from 56 healthy children and adolescents to implement an unsupervised clustering approach that stratify participants based on their resting-state EEG spectral features before the tDCS intervention. We then applied a correlational analysis to characterize the clusters of EEG profiles in terms of participant's difference in the behavioral outcome (accuracy and response time) of the cognitive tasks when performed after a tDCS-sham or a tDCS-active session. Better behavioral performance following the active tDCS session compared to the sham tDCS session is considered a positive intervention response, whilst the reverse is considered a negative one. Optimal results in terms of validity measures was obtained for 4 clusters. These results show that specific EEG-based digital phenotypes can be associated to particular responses. While one cluster presents neurotypical EEG activity, the remaining clusters present non-typical EEG characteristics, which seem to be associated with a positive response. Findings suggest that unsupervised machine learning can be successfully used to stratify and eventually predict responses of individuals to a tDCS treatment.

Association of functional connectivity of the executive control network or default mode network with cognitive impairment in older adults with remitted major depressive disorder or mild cognitive impairment

Major depressive disorder (MDD) is associated with an increased risk of developing dementia. The present study aimed to better understand this risk by comparing resting state functional connectivity (rsFC) in the executive control network (ECN) and the default mode network (DMN) in older adults with MDD or mild cognitive impairment (MCI). Additionally, we examined the association between rsFC in the ECN or DMN and cognitive impairment transdiagnostically. We assessed rsFC alterations in ECN and DMN in 383 participants from five groups at-risk for dementia-remitted MDD with normal cognition (MDD-NC), non-amnestic mild cognitive impairment (naMCI), remitted MDD + naMCI, amnestic MCI (aMCI), and remitted MDD + aMCI-and from healthy controls (HC) or individuals with Alzheimer's dementia (AD). Subject-specific whole-brain functional connectivity maps were generated for each network and group differences in rsFC were calculated. We hypothesized that alteration of rsFC in the ECN and DMN would be progressively larger among our seven groups, ranked from low to high according to their risk for dementia as HC, MDD-NC, naMCI, MDD + naMCI, aMCI, MDD + aMCI, and AD. We also regressed scores of six cognitive domains (executive functioning, processing speed, language, visuospatial memory, verbal memory, and working memory) on the ECN and DMN connectivity maps. We found a significant alteration in the rsFC of the ECN, with post hoc testing showing differences between the AD group and the HC, MDD-NC, or naMCI groups, but no significant alterations in rsFC of the DMN. Alterations in rsFC of the ECN and DMN were significantly associated with several cognitive domain scores transdiagnostically. Our findings suggest that a diagnosis of remitted MDD may not confer functional brain risk for dementia. However, given the association of rs-FC with cognitive performance (i.e., transdiagnostically), rs-FC may help in stratifying this risk among people with MDD and varying degrees of cognitive impairment.

Within-subject reliability of concurrent TMS-fMRI during a single session

Concurrent transcranial magnetic stimulation with functional MRI (concurrent TMS-fMRI) allows real-time causative probing of brain connectivity. However, technical challenges, safety, and tolerability may limit the number of trials employed during a concurrent TMS-fMRI experiment. We leveraged an existing data set with 100 trials of active TMS compared to a sub-threshold control condition to assess the reliability of the evoked BOLD response during concurrent TMS-fMRI. This data will permit an analysis of the minimum number of trials that should be employed in a concurrent TMS-fMRI protocol in order to achieve reliable spatial changes in activity. Single-subject maps of brain activity were created by splitting the trials within the same experimental session into groups of 50, 40, 30, 25, 20, 15, or 10 trials, correlations (R) between t-maps derived from paired subsets of trials within the same individual were calculated as reliability. R was moderate-high for 50 trials (mean R = .695) and decreased as the number of trials decreased. Consistent with previous findings of high individual variability in the spatial patterns of evoked neuronal changes following a TMS pulse, the spatial pattern of Rs differed across participants, but regional R was correlated with the magnitude of TMS-evoked activity. These results demonstrate concurrent TMS-fMRI produces a reliable pattern of activity at the individual level at higher trial numbers, particularly within localized regions. The spatial pattern of reliability is individually idiosyncratic and related to the individual pattern of evoked changes.

Longer illness duration is associated with greater individual variability in functional brain activity in Schizophrenia, but not bipolar disorder

BACKGROUND

Individuals with schizophrenia exhibit greater inter-patient variability in functional brain activity during neurocognitive task performance. Some studies have shown associations of age and illness duration with brain function; however, the association of these variables with variability in brain function activity is not known. In order to better understand the progressive effects of age and illness duration across disorders, we examined the relationship with individual variability in brain activity.

METHODS

Neuroimaging and behavioural data were extracted from harmonized datasets collectively including 212 control participants, 107 individuals with bipolar disorder, and 232 individuals with schizophrenia (total n = 551). Functional activity in response to an N-back working memory task (2-back vs 1-back) was examined. Individual variability was quantified via the correlational distance of fMRI activity between participants; mean correlational distance of one participant in relation to all others was defined as a 'variability score'.

RESULTS

Greater individual variability was found in the schizophrenia group compared to the bipolar disorder and control groups (p = 1.52e-09). Individual variability was significantly associated with aging (p = 0.027), however, this relationship was not different across diagnostic groups. In contrast, in the schizophrenia sample only, a longer illness duration was associated with increased variability (p = 0.027).

CONCLUSION

An increase in variability was observed in the schizophrenia group related to illness duration, beyond the effects of normal aging, implying illness-related deterioration of cognitive networks. This has clinical implications for considering long-term trajectories in schizophrenia and progressive neural and cognitive decline which may be amiable to novel treatments.

A dynamic gradient architecture generates brain activity states

The human brain exhibits a diverse yet constrained range of activity states. While these states can be faithfully represented in a low-dimensional latent space, our understanding of the constitutive functional anatomy is still evolving. Here we applied dimensionality reduction to task-free and task fMRI data to address whether latent dimensions reflect intrinsic systems and if so, how these systems may interact to generate different activity states. We find that each dimension represents a dynamic activity gradient, including a primary unipolar sensory-association gradient underlying the global signal. The gradients appear stable across individuals and cognitive states, while recapitulating key functional connectivity properties including anticorrelation, modularity, and regional hubness. We then use dynamical systems modeling to show that gradients causally interact via state-specific coupling parameters to create distinct brain activity patterns. Together, these findings indicate that a set of dynamic, intrinsic spatial gradients interact to determine the repertoire of possible brain activity states.

Relationship of neurite architecture to brain activity during task-based fMRI

Functional MRI (fMRI) has been widely used to examine changes in neuronal activity during cognitive tasks. Commonly used measures of gray matter macrostructure (e.g., cortical thickness, surface area, volume) do not consistently appear to serve as structural correlates of brain function. In contrast, gray matter microstructure, measured using neurite orientation dispersion and density imaging (NODDI), enables the estimation of indices of neurite density (neurite density index; NDI) and organization (orientation dispersion index; ODI) in gray matter. Our study explored the relationship among neurite architecture, BOLD (blood-oxygen-level-dependent) fMRI, and cognition, using a large sample (n = 750) of young adults of the human connectome project (HCP) and two tasks that index more cortical (working memory) and more subcortical (emotion processing) targeting of brain functions. Using NODDI, fMRI, structural MRI and task performance data, hierarchical regression analyses revealed that higher working memory- and emotion processing-evoked BOLD activity was related to lower ODI in the right DLPFC, and lower ODI and NDI values in the right and left amygdala, respectively. Common measures of brain macrostructure (i.e., DLPFC thickness/surface area and amygdala volume) did not explain any additional variance (beyond neurite architecture) in BOLD activity. A moderating effect of neurite architecture on the relationship between emotion processing task-evoked BOLD response and performance was observed. Our findings provide evidence that neuro-/social-affective cognition-related BOLD activity is partially driven by the local neurite organization and density with direct impact on emotion processing. In vivo gray matter microstructure represents a new target of investigation providing strong potential for clinical translation.

When the time is right: Temporal dynamics of brain activity in healthy aging and dementia

Brain activity and communications are complex phenomena that dynamically unfold over time. However, in contrast with the large number of studies reporting neuroanatomical differences in activation relative to young adults, changes of temporal dynamics of neural activity during normal and pathological aging have been grossly understudied and are still poorly known. Here, we synthesize the current state of knowledge from MEG and EEG studies that aimed at specifying the effects of healthy and pathological aging on local and network dynamics, and discuss the clinical and theoretical implications of these findings. We argue that considering the temporal dynamics of brain activations and networks could provide a better understanding of changes associated with healthy aging, and the progression of neurodegenerative disease. Recent research has also begun to shed light on the association of these dynamics with other imaging modalities and with individual differences in cognitive performance. These insights hold great potential for driving new theoretical frameworks and development of biomarkers to aid in identifying and treating age-related cognitive changes.

Individualizing Representational Similarity Analysis

Representational similarity analysis (RSA) is a popular multivariate analysis technique in cognitive neuroscience that uses functional neuroimaging to investigate the informational content encoded in brain activity. As RSA is increasingly being used to investigate more clinically-geared questions, the focus of such translational studies turns toward the importance of individual differences and their optimization within the experimental design. In this perspective, we focus on two design aspects: applying individual vs. averaged behavioral dissimilarity matrices to multiple participants' neuroimaging data and ensuring the congruency between tasks when measuring behavioral and neural representational spaces. Incorporating these methods permits the detection of individual differences in representational spaces and yields a better-defined transfer of information from representational spaces onto multivoxel patterns. Such design adaptations are prerequisites for optimal translation of RSA to the field of precision psychiatry.