Reproducibility of graph measures at the subject level using resting-state fMRI

Topological Perturbations in the Functional Connectome Support the Deficit/Non-deficit Distinction in Antipsychotic Medication-Naïve First Episode Psychosis Patients

BACKGROUND

Heterogeneity in the etiology, pathophysiology, and clinical features of schizophrenia challenges clinicians and researchers. A helpful approach could be stratifying patients according to the presence or absence of clinical features of the deficit syndrome (DS). DS is characterized by enduring and primary negative symptoms, a clinically less heterogeneous subtype of the illness, and patients with features of DS are thought to present abnormal brain network characteristics, however, this idea has received limited attention. We investigated functional brain network topology in patients displaying deficit features and those who do not.

DESIGN

We applied graph theory analytics to resting-state functional magnetic resonance imaging data of 61 antipsychotic medication-naïve first episode psychosis patients, 18 DS and 43 non-deficit schizophrenia (NDS), and 72 healthy controls (HC). We quantified small-worldness, global and nodal efficiency measures, shortest path length, nodal local efficiency, and synchronization and contrasted them among the 3 groups.

RESULTS

DS presented decreased network integration and segregation compared to HC and NDS. DS showed lower global efficiency, longer global path lengths, and lower global local efficiency. Nodal efficiency was lower and the shortest path length was longer in DS in default mode, ventral attention, dorsal attention, frontoparietal, limbic, somatomotor, and visual networks compared to HC. Compared to NDS, DS showed lower efficiency and longer shortest path length in default mode, limbic, somatomotor, and visual networks.

CONCLUSIONS

Our data supports increasing evidence, based on topological perturbations of the functional connectome, that deficit syndrome may be a distinct form of the illness.

Behaviour-correlated profiles of cerebellar-cerebral functional connectivity observed in independent neurodevelopmental disorder cohorts

The cerebellum, through its connectivity with the cerebral cortex, plays an integral role in regulating cognitive and affective processes, and its dysregulation can result in neurodevelopmental disorder (NDD)-related behavioural deficits. Identifying cerebellar-cerebral functional connectivity (FC) profiles in children with NDDs can provide insight into common connectivity profiles and their correlation to NDD-related behaviours. 479 participants from the Province of Ontario Neurodevelopmental Disorders (POND) network (typically developing = 93, Autism Spectrum Disorder = 172, Attention Deficit/Hyperactivity Disorder = 161, Obsessive-Compulsive Disorder = 53, mean age = 12.2) underwent resting-state functional magnetic resonance imaging and behaviour testing (Social Communication Questionnaire, Toronto Obsessive-Compulsive Scale, and Child Behaviour Checklist - Attentional Problems Subscale). FC components maximally correlated to behaviour were identified using canonical correlation analysis. Results were then validated by repeating the investigation in 556 participants from an independent NDD cohort provided from a separate consortium (Healthy Brain Network (HBN)). Replication of canonical components was quantified by correlating the feature vectors between the two cohorts. The two cerebellar-cerebral FC components that replicated to the greatest extent were correlated to, respectively, obsessive-compulsive behaviour (behaviour feature vectors, rPOND-HBN = -0.97; FC feature vectors, rPOND-HBN = -0.68) and social communication deficit contrasted against attention deficit behaviour (behaviour feature vectors, rPOND-HBN = -0.99; FC feature vectors, rPOND-HBN = -0.78). The statistically stable (|z| > 1.96) features of the FC feature vectors, measured via bootstrap re-sampling, predominantly comprised of correlations between cerebellar attentional and control network regions and cerebral attentional, default mode, and control network regions. In both cohorts, spectral clustering on FC loading values resulted in subject clusters mixed across diagnostic categories, but no cluster was significantly enriched for any given diagnosis as measured via chi-squared test (p > 0.05). Overall, two behaviour-correlated components of cerebellar-cerebral functional connectivity were observed in two independent cohorts. This suggests the existence of generalizable cerebellar network differences that span across NDD diagnostic boundaries.

Modeling the neurocognitive dynamics of language across the lifespan

Healthy aging is associated with a heterogeneous decline across cognitive functions, typically observed between language comprehension and language production (LP). Examining resting-state fMRI and neuropsychological data from 628 healthy adults (age 18-88) from the CamCAN cohort, we performed state-of-the-art graph theoretical analysis to uncover the neural mechanisms underlying this variability. At the cognitive level, our findings suggest that LP is not an isolated function but is modulated throughout the lifespan by the extent of inter-cognitive synergy between semantic and domain-general processes. At the cerebral level, we show that default mode network (DMN) suppression coupled with fronto-parietal network (FPN) integration is the way for the brain to compensate for the effects of dedifferentiation at a minimal cost, efficiently mitigating the age-related decline in LP. Relatedly, reduced DMN suppression in midlife could compromise the ability to manage the cost of FPN integration. This may prompt older adults to adopt a more cost-efficient compensatory strategy that maintains global homeostasis at the expense of LP performances. Taken together, we propose that midlife represents a critical neurocognitive juncture that signifies the onset of LP decline, as older adults gradually lose control over semantic representations. We summarize our findings in a novel synergistic, economical, nonlinear, emergent, cognitive aging model, integrating connectomic and cognitive dimensions within a complex system perspective.

Creativity at rest: Exploring functional network connectivity of creative experts

The neuroscience of creativity seeks to disentangle the complex brain processes that underpin the generation of novel ideas. Neuroimaging studies of functional connectivity, particularly functional magnetic resonance imaging (fMRI), have revealed individual differences in brain network organization associated with creative ability; however, much of the extant research is limited to laboratory-based divergent thinking measures. To overcome these limitations, we compare functional brain connectivity in a cohort of creative experts (n = 27) and controls (n = 26) and examine links with creative behavior. First, we replicate prior findings showing reduced connectivity in visual cortex related to higher creative performance. Second, we examine whether this result is driven by integrated or segregated connectivity. Third, we examine associations between functional connectivity and vivid distal simulation separately in creative experts and controls. In accordance with past work, our results show reduced connectivity to the primary visual cortex in creative experts at rest. Additionally, we observe a negative association between distal simulation vividness and connectivity to the lateral visual cortex in creative experts. Taken together, these results highlight connectivity profiles of highly creative people and suggest that creative thinking may be related to, though not fully redundant with, the ability to vividly imagine the future.

Longitudinal Infant Functional Connectivity Prediction via Conditional Intensive Triplet Network

Longitudinal infant brain functional connectivity (FC) constructed from resting-state functional MRI (rs-fMRI) has increasingly become a pivotal tool in studying the dynamics of early brain development. However, due to various reasons including high acquisition cost, strong motion artifact, and subject dropout, there has been an extreme shortage of usable longitudinal infant rs-fMRI scans to construct longitudinal FCs, which hinders comprehensive understanding and modeling of brain functional development at early ages. To address this issue, in this paper, we propose a novel conditional intensive triplet network (CITN) for longitudinal prediction of the dynamic development of infant FC, which can traverse FCs within a long duration and predict the target FC at any specific age during infancy. Targeting at accurately modeling of the progression pattern of FC, while maintaining the individual functional uniqueness, our model effectively disentangles the intrinsically mixed age-related and identity-related information from the source FC and predicts the target FC by fusing well-disentangled identity-related information with the specific age-related information. Specifically, we introduce an intensive triplet auto-encoder for effective disentanglement of age-related and identity-related information and an identity conditional module to mix identity-related information with designated age-related information. We train the proposed model in a self-supervised way and design downstream tasks to help robustly disentangle age-related and identity-related features. Experiments on 464 longitudinal infant fMRI scans show the superior performance of the proposed method in longitudinal FC prediction in comparison with state-of-the-art approaches.

Astrocyte calcium dysfunction causes early network hyperactivity in Alzheimer's disease

Dysfunctions of network activity and functional connectivity (FC) represent early events in Alzheimer’s disease (AD), but the underlying mechanisms remain unclear. Astrocytes regulate local neuronal activity in the healthy brain, but their involvement in early network hyperactivity in AD is unknown. We show increased FC in the human cingulate cortex several years before amyloid deposition. We find the same early cingulate FC disruption and neuronal hyperactivity in App^NL-F mice. Crucially, these network disruptions are accompanied by decreased astrocyte calcium signaling. Recovery of astrocytic calcium activity normalizes neuronal hyperactivity and FC, as well as seizure susceptibility and day/night behavioral disruptions. In conclusion, we show that astrocytes mediate initial features of AD and drive clinically relevant phenotypes.

•

The cingulate cortex of humans and mice shows early functional deficits in AD

•

Astrocyte calcium signaling is decreased before the presence of amyloid plaques

•

Recovery of astrocyte calcium signals mitigates neuronal hyperactivity

•

Recovery of astrocytes normalizes cingulate connectivity and behavior disruptions

Reproducibility and replicability in neuroimaging data analysis

PURPOSE OF REVIEW

Machine learning solutions are being increasingly used in the analysis of neuroimaging (NI) data, and as a result, there is an increase in the emphasis of the reproducibility and replicability of these data-driven solutions. Although this is a very positive trend, related terminology is often not properly defined, and more importantly, (computational) reproducibility that refers to obtaining consistent results using the same data and the same code is often disregarded.

RECENT FINDINGS

We review the findings of a recent paper on the topic along with other relevant literature, and present two examples that demonstrate the importance of accounting for reproducibility in widely used software for NI data.

SUMMARY

We note that reproducibility should be a first step in all NI data analyses including those focusing on replicability, and introduce available solutions for assessing reproducibility. We add the cautionary remark that when not taken into account, lack of reproducibility can significantly bias all subsequent analysis stages.

Ranking of communities in multiplex spatiotemporal models of brain dynamics

As a relatively new field, network neuroscience has tended to focus on aggregate behaviours of the brain averaged over many successive experiments or over long recordings in order to construct robust brain models. These models are limited in their ability to explain dynamic state changes in the brain which occurs spontaneously as a result of normal brain function. Hidden Markov Models (HMMs) trained on neuroimaging time series data have since arisen as a method to produce dynamical models that are easy to train but can be difficult to fully parametrise or analyse. We propose an interpretation of these neural HMMs as multiplex brain state graph models we term Hidden Markov Graph Models. This interpretation allows for dynamic brain activity to be analysed using the full repertoire of network analysis techniques. Furthermore, we propose a general method for selecting HMM hyperparameters in the absence of external data, based on the principle of maximum entropy, and use this to select the number of layers in the multiplex model. We produce a new tool for determining important communities of brain regions using a spatiotemporal random walk-based procedure that takes advantage of the underlying Markov structure of the model. Our analysis of real multi-subject fMRI data provides new results that corroborate the modular processing hypothesis of the brain at rest as well as contributing new evidence of functional overlap between and within dynamic brain state communities. Our analysis pipeline provides a way to characterise dynamic network activity of the brain under novel behaviours or conditions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s41109-022-00454-2.

Spatial Stability of Functional Networks: A Measure to Assess the Robustness of Graph-Theoretical Metrics to Spatial Errors Related to Brain Parcellation

There is growing interest in studying human brain connectivity and in modelling the brain functional structure as a network. Brain network creation requires parcellation of the cerebral cortex to define nodes. Parcellation might be affected by possible errors due to inter- and intra-subject variability as a consequence of brain structural and physiological characteristics and shape variations related to ageing and diseases, acquisition noise, and misregistration. These errors could induce a knock-on effect on network measure variability. The aim of this study was to investigate spatial stability, a measure of functional connectivity variations induced by parcellation errors. We simulated parcellation variability with random small spatial changes and evaluated its effects on twenty-seven graph-theoretical measures. The study included subjects from three public online datasets. Two brain parcellations were performed using FreeSurfer with geometric atlases. Starting from these, 100 new parcellations were created by increasing the area of 30% of parcels, reducing the area of neighbour parcels, with a rearrangement of vertices. fMRI data were filtered with linear regression, CompCor, and motion correction. Adjacency matrices were constructed with 0.1, 0.2, 0.3, and 0.4 thresholds. Differences in spatial stability between datasets, atlases, and threshold were evaluated. The higher spatial stability resulted for Characteristic-path-length, Density, Transitivity, and Closeness-centrality, and the lower spatial stability resulted for Bonacich and Katz. Multivariate analysis showed a significant effect of atlas, datasets, and thresholds. Katz and Bonacich centrality, which was subject to larger variations, can be considered an unconventional graph measure, poorly implemented in the clinical field and not yet investigated for reliability assessment. Spatial stability (SS) is affected by threshold, and it decreases with increasing threshold for several measures. Moreover, SS seems to depend on atlas choice and scanning parameters. Our study highlights the importance of paying close attention to possible parcellation-related spatial errors, which may affect the reliability of functional connectivity measures.

The Effect of Aging on Brain Glucose Metabolic Connectivity Revealed by [18F]FDG PET-MR and Individual Brain Networks

Contrary to group-based brain connectivity analyses, the aim of this study was to construct individual brain metabolic networks to determine age-related effects on brain metabolic connectivity. Static 40–60 min [¹⁸F]FDG positron emission tomography (PET) images of 67 healthy subjects between 20 and 82 years were acquired with an integrated PET-MR system. Network nodes were defined by brain parcellation using the Schaefer atlas, while connectivity strength between two nodes was determined by comparing the distribution of PET uptake values within each node using a Kullback–Leibler divergence similarity estimation (KLSE). After constructing individual brain networks, a linear and quadratic regression analysis of metabolic connectivity strengths within- and between-networks was performed to model age-dependency. In addition, the age dependency of metrics for network integration (characteristic path length), segregation (clustering coefficient and local efficiency), and centrality (number of hubs) was assessed within the whole brain and within predefined functional subnetworks. Overall, a decrease of metabolic connectivity strength with healthy aging was found within the whole-brain network and several subnetworks except within the somatomotor, limbic, and visual network. The same decrease of metabolic connectivity was found between several networks across the whole-brain network and the functional subnetworks. In terms of network topology, a less integrated and less segregated network was observed with aging, while the distribution and the number of hubs did not change with aging, suggesting that brain metabolic networks are not reorganized during the adult lifespan. In conclusion, using an individual brain metabolic network approach, a decrease in metabolic connectivity strength was observed with healthy aging, both within the whole brain and within several predefined networks. These findings can be used in a diagnostic setting to differentiate between age-related changes in brain metabolic connectivity strength and changes caused by early development of neurodegeneration.

Dysregulation in Sphingolipid Signaling Pathways is Associated With Symptoms and Functional Connectivity of Pain Processing Brain Regions in Provoked Vestibulodynia

Provoked vestibulodynia (PVD) is a chronic pain disorder characterized by local hypersensitivity and severe pain with pressure localized to the vulvar vestibule. Despite decades of study, the lack of identified biomarkers has slowed the development of effective therapies. The primary aim of this study was to use metabolomics to identify novel biochemical mechanisms in vagina and blood underlying brain biomarkers and symptoms in PVD, thereby closing this knowledge gap. Using a cross-sectional case-control observational study design, untargeted and unbiased metabolomic profiling of vaginal fluid and plasma was performed in women with PVD compared to healthy controls. In women with PVD, we also obtained assessments of vulvar pain, vestibular and vaginal muscle tenderness, and 24-hour symptom intensity alongside resting-state brain functional connectivity of brain regions involved in pain processing and modulation. Compared to healthy controls, women with PVD demonstrated differences primarily in vaginal (but not plasma) concentrations of metabolites of the sphingolipid signaling pathways, suggesting localized effects in vagina and vulvar vestibule rather than systemic effects. Our findings reveal that dysregulation of sphingolipid metabolism in PVD is associated with increased vulvar pain and muscle tenderness, sexual dysfunction, and decreased functional connectivity strength in pain processing/modulatory brain regions. This data collectively suggests that alterations in sphingolipid signaling pathways are likely an important molecular biomarker in PVD that could lead to new targets for therapeutic intervention. PERSPECTIVE: This manuscript presents the results of a robust, unbiased molecular assessment of plasma and vaginal fluid samples in women with provoked vestibulodynia compared to healthy controls. The findings suggest that alterations in sphingolipid signaling pathways are associated with symptoms and brain biomarkers and may be an important molecular marker that could provide new targets for therapeutic intervention.

Reproducibility of graph measures at the subject level using resting-state fMRI

INTRODUCTION

Graph metrics have been proposed as potential biomarkers for diagnosis in clinical work. However, before it can be applied in a clinical setting, their reproducibility should be evaluated.

METHODS

This study systematically investigated the effect of two denoising pipelines and different whole-brain network constructions on reproducibility of subject-specific graph measures. We used the multi-session fMRI dataset from the Brain Genomics Superstruct Project consisting of 69 healthy young adults.

RESULTS

In binary networks, the test-retest variability for global measures was large at low density irrespective of the denoising strategy or the type of correlation. Weighted networks showed very low test-retest values (and thus a good reproducibility) for global graph measures irrespective of the strategy used. Comparing the test-retest values for different strategies, there were significant main effects of the type of correlation (Pearson correlation vs. partial correlation), the (partial) correlation value (absolute vs. positive vs. negative), and weight calculation (based on the raw (partial) correlation values vs. based on transformed Z-values). There was also a significant interaction effect between type of correlation and weight calculation. Similarly as for the binary networks, there was no main effect of the denoising pipeline.

CONCLUSION

Our results demonstrated that normalized global graph measures based on a weighted network using the absolute (partial) correlation as weight were reproducible. The denoising pipeline and the granularity of the whole-brain parcellation used to define the nodes were not critical for the reproducibility of normalized graph measures.