Automatic identification of functional clusters in FMRI data using spatial dependence

A survey of brain functional network extraction methods using fMRI data

Functional network (FN) analyses play a pivotal role in uncovering insights into brain function and understanding the pathophysiology of various brain disorders. This paper focuses on classical and advanced methods for deriving brain FNs from functional magnetic resonance imaging (fMRI) data. We systematically review their foundational principles, advantages, shortcomings, and interrelations, encompassing both static and dynamic FN extraction approaches. In the context of static FN extraction, we present hypothesis-driven methods such as region of interest (ROI)-based approaches as well as data-driven methods including matrix decomposition, clustering, and deep learning. For dynamic FN extraction, both window-based and windowless methods are surveyed with respect to the estimation of time-varying FN and the subsequent computation of FN states. We also discuss the scope of application of the various methods and avenues for future improvements.

Identifying canonical and replicable multi-scale intrinsic connectivity networks in 100k+ resting-state fMRI datasets

Despite the known benefits of data-driven approaches, the lack of approaches for identifying functional neuroimaging patterns that capture both individual variations and inter-subject correspondence limits the clinical utility of rsfMRI and its application to single-subject analyses. Here, using rsfMRI data from over 100k individuals across private and public datasets, we identify replicable multi-spatial-scale canonical intrinsic connectivity network (ICN) templates via the use of multi-model-order independent component analysis (ICA). We also study the feasibility of estimating subject-specific ICNs via spatially constrained ICA. The results show that the subject-level ICN estimations vary as a function of the ICN itself, the data length, and the spatial resolution. In general, large-scale ICNs require less data to achieve specific levels of (within- and between-subject) spatial similarity with their templates. Importantly, increasing data length can reduce an ICN's subject-level specificity, suggesting longer scans may not always be desirable. We also find a positive linear relationship between data length and spatial smoothness (possibly due to averaging over intrinsic dynamics), suggesting studies examining optimized data length should consider spatial smoothness. Finally, consistency in spatial similarity between ICNs estimated using the full data and subsets across different data lengths suggests lower within-subject spatial similarity in shorter data is not wholly defined by lower reliability in ICN estimates, but may be an indication of meaningful brain dynamics which average out as data length increases.

Aberrant dynamic functional network connectivity in type 2 diabetes mellitus individuals

An increasing number of recent brain imaging studies are dedicated to understanding the neuro mechanism of cognitive impairment in type 2 diabetes mellitus (T2DM) individuals. In contrast to efforts to date that are limited to static functional connectivity, here we investigate abnormal connectivity in T2DM individuals by characterizing the time-varying properties of brain functional networks. Using group independent component analysis (GICA), sliding-window analysis, and k-means clustering, we extracted thirty-one intrinsic connectivity networks (ICNs) and estimated four recurring brain states. We observed significant group differences in fraction time (FT) and mean dwell time (MDT), and significant negative correlation between the Montreal Cognitive Assessment (MoCA) scores and FT/MDT. We found that in the T2DM group the inter- and intra-network connectivity decreases and increases respectively for the default mode network (DMN) and task-positive network (TPN). We also found alteration in the precuneus network (PCUN) and enhanced connectivity between the salience network (SN) and the TPN. Our study provides evidence of alterations of large-scale resting networks in T2DM individuals and shed light on the fundamental mechanisms of neurocognitive deficits in T2DM.

Multimodal subspace independent vector analysis captures latent subspace structures in large multimodal neuroimaging studies

We present Multimodal Subspace Independent Vector Analysis (MSIVA), a methodology to capture both joint and unique vector sources across multiple data modalities by defining linked and modality-specific subspaces. In particular, MSIVA enables the estimation of independent subspaces of various sizes within modalities and their one-to-one linkage to corresponding subspaces across modalities. We compare MSIVA to a fully unimodal initialization baseline and a fully multimodal initialization baseline, and evaluate all three approaches with five distinct subspace structures on synthetic and neuroimaging datasets. We first demonstrate that MSIVA and the unimodal baseline can identify the correct ground-truth subspace structures from the incorrect ones in multiple synthetic datasets, while the multimodal baseline fails at detecting high-dimensional subspace structures. We then show that MSIVA can better capture the latent subspace structure with the minimum loss value in two large multimodal neuroimaging datasets compared to the unimodal baseline. Our results from subsequent per-subspace canonical correlation analysis (CCA) and brain-phenotype modeling demonstrate that the sources from the optimal subspace structure are strongly associated with phenotype measures, including age, sex and schizophrenia-related effects. Our proposed methodology MSIVA can be applied to capture linked and unique biomarkers from multimodal neuroimaging data.

Shared and Disorder-Specific Alterations of Brain Temporal Dynamics in Obsessive-Compulsive Disorder and Schizophrenia

BACKGROUND

Obsessive-compulsive disorder (OCD) and schizophrenia have distinct but also overlapping symptoms. Few studies have examined the shared and disorder-specific disturbances in dynamic brain function in the 2 disorders.

STUDY DESIGN

Resting-state functional magnetic resonance imaging data of 31 patients with OCD and 49 patients with schizophrenia, all untreated, and 45 healthy controls (HCs) were analyzed using spatial group independent component (IC) analysis. Time-varying degree centrality patterns across the whole brain were clustered into 3 reoccurring states, and state transition metrics were obtained. We further explored regional temporal variability of degree centrality for each IC across all time windows.

STUDY RESULTS

Patients with OCD and patients with schizophrenia both showed decreased occurrence of a state having the highest centrality in the sensorimotor and auditory networks. Additionally, patients with OCD and patients with schizophrenia both exhibited reduced dynamics of degree centrality in the superior frontal gyrus than controls, while dynamic degree centrality of the cerebellum was lower in patients with schizophrenia than with OCD and HCs. Altered dynamics of degree centrality nominally correlated with symptom severity in both patient groups.

CONCLUSIONS

Our study provides evidence of transdiagnostic and clinically relevant functional brain abnormalities across OCD and schizophrenia in neocortex, as well as functional dynamic alterations in the cerebellum specific to schizophrenia. These findings add to the recognition of overlap in neocortical alterations in the 2 disorders, and indicate that cerebellar alterations in schizophrenia may be specifically important in schizophrenia pathophysiology via impact on cerebellar thalamocortical circuitry.

Psychopathic traits and altered resting-state functional connectivity in incarcerated adolescent girls

Previous work in incarcerated boys and adult men and women suggest that individuals scoring high on psychopathic traits show altered resting-state limbic/paralimbic, and default mode functional network properties. However, it is unclear whether similar results extend to high-risk adolescent girls with elevated psychopathic traits. This study examined whether psychopathic traits [assessed via the Hare Psychopathy Checklist: Youth Version (PCL:YV)] were associated with altered inter-network connectivity, intra-network connectivity (i.e., functional coherence within a network), and amplitude of low-frequency fluctuations (ALFFs) across resting-state networks among high-risk incarcerated adolescent girls (n = 40). Resting-state networks were identified by applying group independent component analysis (ICA) to resting-state fMRI scans, and a priori regions of interest included limbic, paralimbic, and default mode network components. We tested the association of psychopathic traits (PCL:YV Factor 1 measuring affective/interpersonal traits and PCL:YV Factor 2 assessing antisocial/lifestyle traits) to these three resting-state measures. PCL:YV Factor 1 scores were associated with increased low-frequency and decreased high-frequency fluctuations in components corresponding to the default mode network, as well as increased intra-network FNC in components corresponding to cognitive control networks. PCL:YV Factor 2 scores were associated with increased low-frequency fluctuations in sensorimotor networks and decreased high-frequency fluctuations in default mode, sensorimotor, and visual networks. Consistent with previous analyses in incarcerated adult women, our results suggest that psychopathic traits among incarcerated adolescent girls are associated with altered intra-network ALFFs-primarily that of increased low-frequency and decreased high-frequency fluctuations-and connectivity across multiple networks including paralimbic regions. These results suggest stable neurobiological correlates of psychopathic traits among women across development.

Aberrant resting-state functional connectivity associated with childhood trauma among juvenile offenders

Individuals with history of childhood trauma are characterized by aberrant resting-state limbic and paralimbic functional network connectivity. However, it is unclear whether specific subtypes of trauma (i.e., experienced vs observed or community) showcase differential effects. This study examined whether subtypes of childhood trauma (assessed via the Trauma Checklist [TCL] 2.0) were associated with aberrant intra-network amplitude of fluctuations and connectivity (i.e., functional coherence within a network), and inter-network connectivity across resting-state networks among incarcerated juvenile males (n = 179). Subtypes of trauma were established via principal component analysis of the TCL 2.0 and resting-state networks were identified by applying group independent component analysis to resting-state fMRI scans. We tested the association of subtypes of childhood trauma (i.e., TCL Factor 1 measuring experienced trauma and TCL Factor 2 assessing community trauma), and TCL Total scores to the aforementioned functional connectivity measures. TCL Factor 2 scores were associated with increased high-frequency fluctuations and increased intra-network connectivity in cognitive control, auditory, and sensorimotor networks, occurring primarily in paralimbic regions. TCL Total scores exhibited similar neurobiological patterns to TCL Factor 2 scores (with the addition of aberrant intra-network connectivity in visual networks), and no significant associations were found for TCL Factor 1. Consistent with previous analyses of community samples, our results suggest that childhood trauma among incarcerated juvenile males is associated with aberrant intra-network amplitude of fluctuations and connectivity across multiple networks including predominately paralimbic regions. Our results highlight the importance of accounting for traumatic loss, observed trauma, and community trauma in assessing neurobiological aberrances associated with adverse experiences in childhood, as well as the value of trained-rater trauma assessments compared to self-report.

Cortical thickness analysis in temporal lobe epilepsy using fully Bayesian spectral method in magnetic resonance imaging

BACKGROUND

Temporal lobe epilepsy (TLE) is the most common type of epilepsy associated with changes in the cerebral cortex throughout the brain. Magnetic resonance imaging (MRI) is widely used for detecting such anomalies; nevertheless, it produces spatially correlated data that cannot be considered by the usual statistical models. This study aimed to compare cortical thicknesses between patients with TLE and healthy controls by considering the spatial dependencies across different regions of the cerebral cortex in MRI.

METHODS

In this study, T1-weighted MRI was performed on 20 healthy controls and 33 TLE patients. Nineteen patients had a left TLE and 14 had a right TLE. Cortical thickness was measured for all individuals in 68 regions of the cerebral cortex based on images. Fully Bayesian spectral method was utilized to compare the cortical thickness of different brain regions between groups. Neural networks model was used to classify the patients using the identified regions.

RESULTS

For the left TLE patients, cortical thinning was observed in bilateral caudal anterior cingulate, lateral orbitofrontal (ipsilateral), the bilateral rostral anterior cingulate, frontal pole and temporal pole (ipsilateral), caudal middle frontal and rostral middle frontal (contralateral side). For the right TLE patients, cortical thinning was only observed in the entorhinal area (ipsilateral). The AUCs of the neural networks for classification of left and right TLE patients versus healthy controls were 0.939 and 1.000, respectively.

CONCLUSION

Alteration of cortical gray matter thickness was evidenced as common effect of epileptogenicity, as manifested by the patients in this study using the fully Bayesian spectral method by taking into account the complex structure of the data.

Inter-individual differences in baseline dynamic functional connectivity are linked to cognitive aftereffects of tDCS

Transcranial direct current stimulation (tDCS) has the potential to modulate cognitive training in healthy aging; however, results from various studies have been inconsistent. We hypothesized that inter-individual differences in baseline brain state may contribute to the varied results. We aimed to explore whether baseline resting-state dynamic functional connectivity (rs-dFC) and/or conventional resting-state static functional connectivity (rs-sFC) may be related to the magnitude of cognitive aftereffects of tDCS. To achieve this aim, we used data from our double-blind randomized sham-controlled cross-over tDCS trial in 25 healthy seniors in which bifrontal tDCS combined with cognitive training had induced significant behavioral aftereffects. We performed a backward regression analysis including rs-sFC/rs-dFC measures to explain the variability in the magnitude of tDCS-induced improvements in visual object-matching task (VOMT) accuracy. Rs-dFC analysis revealed four rs-dFC states. The occurrence rate of a rs-dFC state 4, characterized by a high correlation between the left fronto-parietal control network and the language network, was significantly associated with tDCS-induced VOMT accuracy changes. The rs-sFC measure was not significantly associated with the cognitive outcome. We show that flexibility of the brain state representing readiness for top-down control of object identification implicated in the studied task is linked to the tDCS-enhanced task accuracy.

Aberrant resting-state functional connectivity in incarcerated women with elevated psychopathic traits

Previous work in incarcerated men suggests that individuals scoring high on psychopathy exhibit aberrant resting-state paralimbic functional network connectivity (FNC). However, it is unclear whether similar results extend to women scoring high on psychopathy. This study examined whether psychopathic traits [assessed via the Hare Psychopathy Checklist - Revised (PCL-R)] were associated with aberrant inter-network connectivity, intra-network connectivity (i.e., functional coherence within a network), and amplitude of fluctuations across limbic and surrounding paralimbic regions among incarcerated women (n = 297). Resting-state networks were identified by applying group Independent Component Analysis to resting-state fMRI scans. We tested the association of psychopathic traits (PCL-R Factor 1 measuring interpersonal/affective psychopathic traits and PCL-R Factor 2 assessing lifestyle/antisocial psychopathic traits) to the three FNC measures. PCL-R Factor 1 scores were associated with increased low-frequency fluctuations in executive control and attentional networks, decreased high-frequency fluctuations in executive control and visual networks, and decreased intra-network FNC in default mode network. PCL-R Factor 2 scores were associated with decreased high-frequency fluctuations and default mode networks, and both increased and decreased intra-network functional connectivity in visual networks. Similar to previous analyses in incarcerated men, our results suggest that psychopathic traits among incarcerated women are associated with aberrant intra-network amplitude fluctuations and connectivity across multiple networks including limbic and surrounding paralimbic regions.

Multimodel Order Independent Component Analysis: A Data-Driven Method for Evaluating Brain Functional Network Connectivity Within and Between Multiple Spatial Scales

Background: While functional connectivity is widely studied, there has been little work studying functional connectivity at different spatial scales. Likewise, the relationship of functional connectivity between spatial scales is unknown. Methods: We proposed an independent component analysis (ICA)-based approach to capture information at multiple-model orders (component numbers), and to evaluate functional network connectivity (FNC) both within and between model orders. We evaluated the approach by studying group differences in the context of a study of resting-state functional magnetic resonance imaging (rsfMRI) data collected from schizophrenia (SZ) individuals and healthy controls (HC). The predictive ability of FNC at multiple spatial scales was assessed using support vector machine-based classification. Results: In addition to consistent predictive patterns at both multiple-model orders and single-model orders, unique predictive information was seen at multiple-model orders and in the interaction between model orders. We observed that the FNC between model orders 25 and 50 maintained the highest predictive information between HC and SZ. Results highlighted the predictive ability of the somatomotor and visual domains both within and between model orders compared with other functional domains. Also, subcortical-somatomotor, temporal-somatomotor, and temporal-subcortical FNCs had relatively high weights in predicting SZ. Conclusions: In sum, multimodel order ICA provides a more comprehensive way to study FNC, produces meaningful and interesting results, which are applicable to future studies. We shared the spatial templates from this work at different model orders to provide a reference for the community, which can be leveraged in regression-based or fully automated (spatially constrained) ICA approaches. Impact statement Multimodel order independent component analysis (ICA) provides a comprehensive way to study brain functional network connectivity within and between multiple spatial scales, highlighting findings that would have been ignored in single-model order analysis. This work expands upon and adds to the relatively new literature on resting functional magnetic resonance imaging-based classification and prediction. Results highlighted the differentiating power of specific intrinsic connectivity networks on classifying brain disorders of schizophrenia patients and healthy participants, at different spatial scales. The spatial templates from this work provide a reference for the community, which can be leveraged in regression-based or fully automated ICA approaches.

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.

White matter substrates of functional connectivity dynamics in the human brain

The contribution of structural connectivity to functional connectivity dynamics is still far from being elucidated. Herein, we applied track-weighted dynamic functional connectivity (tw-dFC), a model integrating structural, functional, and dynamic connectivity, on high quality diffusion weighted imaging and resting-state fMRI data from two independent repositories. The tw-dFC maps were analyzed using independent component analysis, aiming at identifying spatially independent white matter components which support dynamic changes in functional connectivity. Each component consisted of a spatial map of white matter bundles that show consistent fluctuations in functional connectivity at their endpoints, and a time course representative of such functional activity. These components show high intra-subject, inter-subject, and inter-cohort reproducibility. We provided also converging evidence that functional information about white matter activity derived by this method can capture biologically meaningful features of brain connectivity organization, as well as predict higher-order cognitive performance.

Multi-spatial-scale dynamic interactions between functional sources reveal sex-specific changes in schizophrenia

We introduce an extension of independent component analysis (ICA), called multiscale ICA, and design an approach to capture dynamic functional source interactions within and between multiple spatial scales. Multiscale ICA estimates functional sources at multiple spatial scales without imposing direct constraints on the size of functional sources, overcomes the limitation of using fixed anatomical locations, and eliminates the need for model-order selection in ICA analysis. We leveraged this approach to study sex-specific and sex-common connectivity patterns in schizophrenia. Results show dynamic reconfiguration and interaction within and between multi-spatial scales. Sex-specific differences occur (a) within the subcortical domain, (b) between the somatomotor and cerebellum domains, and (c) between the temporal domain and several others, including the subcortical, visual, and default mode domains. Most of the sex-specific differences belong to between-spatial-scale functional interactions and are associated with a dynamic state with strong functional interactions between the visual, somatomotor, and temporal domains and their anticorrelation patterns with the rest of the brain. We observed significant correlations between multi-spatial-scale functional interactions and symptom scores, highlighting the importance of multiscale analyses to identify potential biomarkers for schizophrenia. As such, we recommend such analyses as an important option for future functional connectivity studies.

Biclustering fMRI time series: a comparative study

BACKGROUND

The effectiveness of biclustering, simultaneous clustering of rows and columns in a data matrix, was shown in gene expression data analysis. Several researchers recognize its potentialities in other research areas. Nevertheless, the last two decades have witnessed the development of a significant number of biclustering algorithms targeting gene expression data analysis and a lack of consistent studies exploring the capacities of biclustering outside this traditional application domain.

RESULTS

This work evaluates the potential use of biclustering in fMRI time series data, targeting the Region × Time dimensions by comparing seven state-in-the-art biclustering and three traditional clustering algorithms on artificial and real data. It further proposes a methodology for biclustering evaluation beyond gene expression data analysis. The results discuss the use of different search strategies in both artificial and real fMRI time series showed the superiority of exhaustive biclustering approaches, obtaining the most homogeneous biclusters. However, their high computational costs are a challenge, and further work is needed for the efficient use of biclustering in fMRI data analysis.

CONCLUSIONS

This work pinpoints avenues for the use of biclustering in spatio-temporal data analysis, in particular neurosciences applications. The proposed evaluation methodology showed evidence of the effectiveness of biclustering in finding local patterns in fMRI time series data. Further work is needed regarding scalability to promote the application in real scenarios.

Cohesive parcellation of the human brain using resting-state fMRI

BACKGROUND

The data burden for resting-state fMRI analysis rises with increasing resolutions available at ultrahigh fields. Therefore, a fundamental preprocessing step in brain network analysis is to reduce the data, usually by performing some kind of data parcellation. Most functional parcellations based on rsfMRI connectivity are synthesized from the dense connectome. In contrast, most network analyses begin by reducing each parcel to a single exemplar time series. This disconnect between parcel formation and usage assumes that parcel exemplars adequately represent their member voxels, which is not always the case for commonly used parcellations.

Abnormal Dynamic Functional Connectivity in Patients With End-Stage Renal Disease

Dynamic functional connectivity (FC) analysis can capture time-varying properties of connectivity; however, studies focusing on dynamic FC in patients with end-stage renal disease (ESRD) are very limited. This is the first study to explore the dynamic aspects of whole-brain FC and topological properties in ESRD patients. Resting-state functional magnetic resonance imaging data were acquired from 100 ESRD patients [50 hemodialysis (HD) patients and 50 non-dialysis patients] and 64 healthy controls (HCs). Independent component analysis, a sliding-window approach and graph-theory methods were used to study the dynamic FC properties. The intrinsic brain FC were clustered into four configuration states. Compared with HCs, both patient groups spent longer time in State 3, in which decreased FC between subnetworks of the default mode network (DMN) and between the dorsal DMN and language network was observed in these patients, and a further reduction in FC between the DMN subnetworks was found in HD patients compared to non-dialysis patients. The number of transitions and the variability of global and local efficiency progressively decreased from that in HCs to that of non-dialysis patients to that of HD patients. The completion time of Trail Making Test A and Trail Making Test B positively correlated with the mean dwell time of State 3 and negatively correlated with the number of transitions in ESRD patients. Our findings suggest impaired functional flexibility of network connections and state-specific FC disruptions in patients with ESRD, which may underlie their cognitive deficits. HD may have an adverse effect on time-varying FC.

Individualized spatial network predictions using Siamese convolutional neural networks: A resting-state fMRI study of over 11,000 unaffected individuals

Individuals can be characterized in a population according to their brain measurements and activity, given the inter-subject variability in brain anatomy, structure-function relationships, or life experience. Many neuroimaging studies have demonstrated the potential of functional network connectivity patterns estimated from resting functional magnetic resonance imaging (fMRI) to discriminate groups and predict information about individual subjects. However, the predictive signal present in the spatial heterogeneity of brain connectivity networks is yet to be extensively studied. In this study, we investigate, for the first time, the use of pairwise-relationships between resting-state independent spatial maps to characterize individuals. To do this, we develop a deep Siamese framework comprising three-dimensional convolution neural networks for contrastive learning based on individual-level spatial maps estimated via a fully automated fMRI independent component analysis approach. The proposed framework evaluates whether pairs of spatial networks (e.g., visual network and auditory network) are capable of subject identification and assesses the spatial variability in different network pairs' predictive power in an extensive whole-brain analysis. Our analysis on nearly 12,000 unaffected individuals from the UK Biobank study demonstrates that the proposed approach can discriminate subjects with an accuracy of up to 88% for a single network pair on the test set (best model, after several runs), and 82% average accuracy at the subcortical domain level, notably the highest average domain level accuracy attained. Further investigation of our network's learned features revealed a higher spatial variability in predictive accuracy among younger brains and significantly higher discriminative power among males. In sum, the relationship among spatial networks appears to be both informative and discriminative of individuals and should be studied further as putative brain-based biomarkers.

Altered brain functional network dynamics in classic trigeminal neuralgia: a resting-state functional magnetic resonance imaging study

Background

Accumulating studies have indicated a wide range of brain alterations with respect to the structure and function of classic trigeminal neuralgia (CTN). Given the dynamic nature of pain experience, the exploration of temporal fluctuations in interregional activity covariance may enhance the understanding of pain processes in the brain. The present study aimed to characterize the temporal features of functional connectivity (FC) states as well as topological alteration in CTN.

Methods

Resting-state functional magnetic resonance imaging and three-dimensional T1-weighted images were obtained from 41 CTN patients and 43 matched healthy controls (HCs). After group independent component analysis, sliding window based dynamic functional network connectivity (dFNC) analysis was applied to investigate specific FC states and related temporal properties. Then, the dynamics of the whole brain topological organization were estimated by calculating the coefficient of variation of graph-theoretical properties. Further correlation analyses were performed between all these measurements and clinical data.

Results

Two distinct states were identified. Of these, the state 2, characterized by complicated coupling between default mode network (DMN) and cognitive control network (CC) and tight connections within DMN, was expressed more in CTN patients and presented as increased fractional windows and dwell time. Moreover, patients switched less frequently between states than HCs. Regarding the dynamic topological analysis, disruptions in global graph-theoretical properties (including network efficiency and small-worldness) were observed in patients, coupled with decreased variability in nodal efficiency of anterior cingulate cortex (ACC) in the salience network (SN) and the thalamus and caudate nucleus in the subcortical network (SC). The variation of topological properties showed negative correlation with disease duration and attack frequency.

Conclusions

The present study indicated disrupted flexibility of brain topological organization under persistent noxious stimulation and further highlighted the important role of “dynamic pain connectome” regions (including DMN/CC/SN) in the pathophysiology of CTN from the temporal fluctuation aspect. Additionally, the findings provided supplementary evidence for current knowledge about the aberrant cortical-subcortical interaction in pain development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s10194-021-01354-z.

How recent learning shapes the brain: Memory-dependent functional reconfiguration of brain circuits

The process of storing recently encoded episodic mnestic traces so that they are available for subsequent retrieval is accompanied by specific brain functional connectivity (FC) changes. In this fMRI study, we examined the early processing of memories in twenty-eight healthy participants performing an episodic memory task interposed between two resting state sessions. Memory performance was assessed through a forced-choice recognition test after the scanning sessions. We investigated resting state system configuration changes via Independent Component Analysis by cross-modeling baseline resting state spatial maps onto the post-encoding resting state, and post-encoding resting state spatial maps onto baseline. We identified both persistent and plastic components of the overall brain functional configuration between baseline and post-encoding. While FC patterns within executive, default mode, and cerebellar circuits persisted from baseline to post-encoding, FC within the visual circuit changed. A significant session × performance interaction characterized medial temporal lobe and prefrontal cortex FC with the visual circuit, as well as thalamic FC within the executive control system. Findings reveal early-stage FC changes at the system-level subsequent to a learning experience and associated with inter-individual variation in memory performance.

Phase fMRI defines brain resting-state functional hubs within central and posterior regions

From a brain functional connectivity (FC) matrix, we can identify the hub nodes by a new method of eigencentrality mapping, which not only counts for one node's centrality but also all other nodes' centrality values through correlation connections in an eigenvector of the FC matrix. For the resting-state functional MRI (fMRI) data (complex-valued EPI images in nature), both magnitude and phase images are useful for brain FC analysis. We herein report on brain functional hubness analysis by constructing the FC matrix from phase fMRI data and identifying the hub nodes by eigencentrality mapping. In our study, we collected a cohort of 160 complex-valued fMRI dataset (consisting of magnitude and phase in pairs), and performed independent component analysis (ICA), FC matrix calculation (in size of 50 × 50) and FC matrix eigen decomposition; thereby obtained the 50-node eigencentrality values in the eigenvector associated with the largest eigenvalue. We also compared the hub structures inferred from FC matrices under different thresholding. Alternatively, we obtained the geometric hubs among p value the 50 nodes involved in the FC matrix through the use of harmonic centrality metric. Our results showed that phase fMRI data analysis defines the resting-state brain functional hubs primarily in the central region (subcortex) and the posterior region (parieto-occipital lobes and cerebella). The brain central hubness was supported by the geometric central hubness, which, however, is distinct from the magnitude-inferred hubness in brain superior region (frontal and parietal lobes). Our findings pose a new understanding of (or a debate over) brain functional connectivity architecture.

Modular and state-relevant functional network connectivity in high-frequency eyes open vs eyes closed resting fMRI data

BACKGROUND

Resting-state fMRI (rs-fMRI) is employed to assess "functional connections" of signal between brain regions. However, multiple rs-fMRI paradigms and data-filtering strategies have been used, highlighting the need to explore BOLD signal across the spectrum. Rs-fMRI data is typically filtered at frequencies ranging between 0.008∼0.2 Hz to mitigate nuisance signal (e.g. cardiac, respiratory) and maximize neuronal BOLD signal. However, some argue neuronal BOLD signal may be parsed at higher frequencies.

NEW METHOD

To assess the contributions of rs-fMRI along the BOLD spectra on functional network connectivity (FNC) matrices, a spatially constrained independent component analysis (ICA) was performed at seven different frequency "bins", after which FNC values and FNC measures of matrix-randomness were assessed using linear mixed models.

RESULTS

Results show FNCs at higher-frequency bins display similar randomness to those from the typical frequency bins (0.01-0.15), while the largest values are in the 0.31-0.46 Hz bin. Eyes open (EO) vs eyes closed (EC) comparison found EC was less random than EO across most frequency bins. Further, FNC was greater in EC across auditory and cognitive control pairings while EO values were greater in somatomotor, visual, and default mode FNC.

COMPARISON WITH EXISTING METHODS

Effect sizes of frequency and resting-state paradigm vary from small to large, but the most notable results are specific to frequency ranges and resting-state paradigm with artifacts like motion displaying negligible effect sizes.

CONCLUSIONS

These suggest unique information may be derived from FNC matrices across frequencies and paradigms, but additional data is necessary prior to any definitive conclusions.

Altered brain functional network dynamics in obsessive-compulsive disorder

Obsessive–compulsive disorder (OCD) is a debilitating and disabling neuropsychiatric disorder, whose neurobiological basis remains unclear. Although traditional static resting‐state magnetic resonance imaging (rfMRI) studies have found aberrant functional connectivity (FC) in OCD, alterations in whole‐brain FC and topological properties in the context of brain dynamics remain relatively unexplored. The rfMRI data of 29 patients with OCD and 40 healthy controls were analyzed using group independent component analysis to obtain independent components (ICs) and a sliding‐window approach to generate dynamic functional connectivity (dFC) matrices. dFC patterns were clustered into three reoccurring states, and state transition metrics were obtained. Then, graph‐theory methods were applied to dFC matrices to calculate the variability of network topological organization. The occurrence of a state (State 1) with the highest modularity index and lowest mean FC between networks was increased significantly in OCD, and the fractional time in brain State 1 was positively correlated with anxiety level in patients. State 1 was characterized by having positive connections within default mode (DMN) and salience networks (SAN), and negative coupling between the two networks. Additionally, ICs belonging to DMN and SAN showed lower temporal variability of nodal degree centrality and efficiency in patients, which was related to longer illness duration and higher current obsession ratings. Our results provide evidence of clinically relevant aberrant dynamic brain activity in OCD. Increased functional segregation among networks and impaired functional flexibility in connections among brain regions in DMN and SAN may play important roles in the neuropathology of OCD.

Dynamic Resting-State Connectivity Differences in Eyes Open Versus Eyes Closed Conditions

Introduction: Previous studies have shown significant conditional differences between eyes open, fixated at an image (EO) and eyes closed (EC) in the acquired resting-state functional magnetic resonance imaging (rs-fMRI) data. Aim: We recently showed significant functional network connectivity (FNC) differences between EO and EC across a variety of networks. In this study, we aim at further evaluating differences in dynamic FNC (dFNC) between EO and EC. Materials and Methods: Rs-fMRI were collected from adolescents aged 9-15 years old during both EO and EC conditions, and dFNC was calculated by using the independent component analysis framework. Results: We found that out of five states (clusters), state 1 was observed to be more dominant in the EO condition, whereas state 2 was observed to be more dominant in the EC condition. States 1 and 2 showed significant differences in the mean dwell time based on false discovery rate, and states 1, 2, 3, and 4 differed in the frequency of occurrences. These results are consistent with our previous study of static connectivity in suggesting that EO and EC differences not only are relatively strong but also importantly reveal that these differences vary over time, especially in one particularly transient connectivity pattern. Conclusion: Our results manifest as changes in the proportion of time spent in unique functional connectivity patterns, and they show unique transient functional connectivity patterns in a subset of identified states. Overall, our findings indicate that both static and dynamic rs-fMRI connectivity patterns are strongly impacted by basic conditional differences such as EO and EC. Impact statement Our findings not only suggest that eyes open, fixated at an image (EO) and eyes closed (EC) condition-related resting state functional magnetic resonance imaging differences are relatively strong, but they also reveal an important attribute of these conditions that these differences vary over time, especially in one particularly transient connectivity pattern. Our results manifest as changes in the proportion of time spent in unique functional connectivity patterns, and they show unique transient functional connectivity patterns in a subset of identified states. We believe there is benefit in having the EO/EC as a contrast of interest in future studies, if time allows.

NeuroMark: An automated and adaptive ICA based pipeline to identify reproducible fMRI markers of brain disorders

•

Propose a new pipeline to link brain changes among different datasets, studies, and disorders.

•

Identify reproducible biomarkers in schizophrenia using independent data.

•

Find both common and unique brain impairments in schizophrenia and autism.

•

Reveal gradual changes from healthy controls to mild cognitive impairment to Alzheimer’s disease.

•

Obtain high classification accuracy (~90%) between bipolar disorder and major depressive disorder.

Many mental illnesses share overlapping or similar clinical symptoms, confounding the diagnosis. It is important to systematically characterize the degree to which unique and similar changing patterns are reflective of brain disorders. Increasing sharing initiatives on neuroimaging data have provided unprecedented opportunities to study brain disorders. However, it is still an open question on replicating and translating findings across studies. Standardized approaches for capturing reproducible and comparable imaging markers are greatly needed. Here, we propose a pipeline based on the priori-driven independent component analysis, NeuroMark, which is capable of estimating brain functional network measures from functional magnetic resonance imaging (fMRI) data that can be used to link brain network abnormalities among different datasets, studies, and disorders. NeuroMark automatically estimates features adaptable to each individual subject and comparable across datasets/studies/disorders by taking advantage of the reliable brain network templates extracted from 1828 healthy controls as guidance. Four studies including 2442 subjects were conducted spanning six brain disorders (schizophrenia, autism spectrum disorder, mild cognitive impairment, Alzheimer’s disease, bipolar disorder, and major depressive disorder) to evaluate validity of the proposed pipeline from different perspectives (replication of brain abnormalities, cross-study comparison, identification of subtle brain changes, and multi-disorder classification using identified biomarkers). Our results highlight that NeuroMark effectively identified replicated brain network abnormalities of schizophrenia across different datasets; revealed interesting neural clues on the overlap and specificity between autism and schizophrenia; demonstrated brain functional impairments present to varying degrees in mild cognitive impairments and Alzheimer's disease; and captured biomarkers that achieved good performance in classifying bipolar disorder and major depressive disorder.

Predicting Biological Gender and Intelligence From fMRI via Dynamic Functional Connectivity

OBJECTIVE

This paper explores the predictive capability of dynamic functional connectivity extracted from functional magnetic resonance imaging (fMRI) of the human brain, in contrast to static connectivity used in past research.

METHODS

Several state-of-the-art features extracted from static functional connectivity of the brain are employed to predict biological gender and intelligence using publicly available Human Connectome Project (HCP) database. Next, a novel tensor parallel factor (PARAFAC) decomposition model is proposed to decompose sequence of dynamic connectivity matrices into common connectivity components that are orthonormal to each other, common time-courses, and corresponding distinct subject-wise weights. The subject-wise loading of the components are employed to predict biological gender and intelligence using a random forest classifier (respectively, regressor) using 5-fold cross-validation.

RESULTS

The results demonstrate that dynamic functional connectivity can indeed classify biological gender with a high accuracy (0.94, where male identification accuracy was 0.87 and female identification accuracy was 0.97). It can also predict intelligence with less normalized mean square error (0.139 for fluid intelligence and 0.031 for fluid ability metrics) compared with other functional connectivity measures (the nearest mean square error were 0.147 and 0.037 for fluid intelligence and fluid ability metrics, respectively, using static connectivity approaches).

CONCLUSION

Our work is an important milestone for the understanding of non-stationary behavior of hemodynamic blood-oxygen level dependent (BOLD) signal in brain and how they are associated with biological gender and intelligence.

SIGNIFICANCE

The paper demonstrates that dynamic behavior of brain can contribute substantially towards forming a fingerprint of biological gender and intelligence.

Neurotransmitter imbalance dysregulates brain dynamic fluidity in frontotemporal degeneration

Frontotemporal degeneration (FTD) is characterized by reduced global brain flexibility along with GABAergic/glutamatergic neurotransmitter deficits. We aimed to assess the relationship between dynamical properties of time-varying whole-brain network connectivity as well as static large-scale networks and neurotransmitter imbalance using resting-state functional MRI and transcranial magnetic stimulation (TMS) in sixty-six patients with FTD. We assessed GABAergic and glutamatergic neurotransmission by TMS, considering short- and long-interval intracortical inhibition and intracortical facilitation, and large-scale networks connectivity as well as four indexes of meta-state dynamic fluidity: (1) number of distinct meta-states, (2) number of switches from one meta-state to another, (3) span of the realized meta-states, and (4) total distance traveled in the state space. No significant correlations between TMS parameters and large-scale networks connectivity were observed. However, we observed a significant correlation between short-interval intracortical inhibition-intracortical facilitation and four meta-states (all indexes p < 0.02, false discovery rate-corrected). This study suggests that neurotransmitter imbalance dysregulates brain dynamic fluidity, linking microscopic and macroscopic changes in FTD.

Static and dynamic functional connectivity analysis of cerebrovascular reactivity: An fMRI study

BACKGROUND

Cerebrovascular reactivity (CVR) is an important aspect of brain function, and as such it is important to understand relationship between CVR and functional connectivity.

METHODS

This research studied the role of CVR, or the brain's ability to react to vasoactive stimuli on brain functional connectivity by scanning subjects with blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) while they periodically inhale room air and a CO 2-enriched gas mixture. We developed a new metric to measure the effect of CVR on each intrinsic connectivity network (ICN), which contrasts to voxel-wise CVR. We also studied the changes in whole-brain connectivity patterns using both static functional network connectivity (sFNC) and dynamic FNC (dFNC).

RESULTS

We found that network connectivity is generally weaker during vascular dilation, which is supported by previous research. The dFNC analysis revealed that participants did not return to the pre-CO 2 inhalation state, suggesting that one-minute periods of room-air inhalation is not enough for the CO 2 effect to fully dissipate.

CONCLUSIONS

Cerebrovascular reactivity is one tool that the cerebrovascular system uses to ensure the constant, finely-tuned flow of oxygen to function properly. Understanding the relationship between CVR and brain dynamism can provide unique information about cerebrovascular diseases and general brain function. We observed that CVR has a wide, but consistent relationship to connectivity patterns between functional networks.

Identifying commonality and specificity across psychosis sub-groups via classification based on features from dynamic connectivity analysis

It is difficult to distinguish schizophrenia (SZ), schizoaffective disorder (SAD), and bipolar disorder with psychosis (BPP) as their clinical diagnoses rely on symptoms that overlap. In this paper, we investigate if there is biological evidence to support the symptom-based clinical categories by looking across the three disorders using dynamic connectivity measures, and provide meaningful characteristics on which brain functional connectivity measures are commonly or uniquely impaired. Large-sample functional magnetic resonance image (fMRI) datasets from 623 subjects including 238 healthy controls (HCs), 113 SZ patients, 132 SAD patients, and 140 BPP patients were analyzed. First, we computed whole-brain dynamic functional connectivity (DFC) using a sliding-window technique, and then extracted the individual connectivity states by applying our previously proposed decomposition-based DFC analysis method. Next, with the features from the dominant connectivity state, we assessed the clinical categories by performing both four-group (SZ, SAD, BPP and healthy control groups) and pair-wise classification using a support vector machine within cross-validation. Furthermore, we comprehensively summarized the shared and unique connectivity alterations among the disorders. In terms of the classification performance, our method achieved 69% in the four-group classification and >80% in the between-group classifications for the mean overall accuracy; and yielded 66% in the four-group classification and >80% in the between-group classifications for the mean balanced accuracy. Through summarizing the features that were automatically selected in the classifications, we found that among the three symptom-related disorders, their disorder-common impairments primarily included the decreased connectivity strength between thalamus and cerebellum and the increased strength between postcentral gyrus and thalamus. The disorder-unique changes included more various brain regions, mainly in the temporal and frontal gyrus. Our work demonstrates that dynamic functional connectivity provides biological evidence that both common and unique impairments exist in psychosis sub-groups.

Determining the number of states in dynamic functional connectivity using cluster validity indexes

BACKGROUND

Clustering analysis is employed in brain dynamic functional connectivity (dFC) to cluster the data into a set of dynamic states. These states correspond to different patterns of functional connectivity that iterate through time. Although several cluster validity index (CVI) methods to determine the best clustering partition exists, the appropriateness of methods to apply in the case of dynamic connectivity analysis has not been determined.

NEW METHOD

Currently employed indexes do not provide a crisp answer on what is the best number of clusters. In addition, there is a lack of CVI testing in the context of dFC data. This work tests a comprehensive set of twenty four cluster validity indexes applied to addiction data and suggest the best ones for clustering dynamic functional connectivity.

RESULTS

Out of the twenty four considered CVIs, Davies-Bouldin and Ray-Turi were the most suitable methods to find the number of clusters in both simulation and real data. The solution for these two CVIs is to find a local minimum critical point, which can be automated using computational algorithms.

COMPARISON WITH EXISTING METHODS

Elbow-Criterion, Silhouette and GAP-Statistic methods have been widely used in dFC studies. These methods are included among the tested CVIs where the performances of all twenty four CVIs are compared.

CONCLUSIONS

Davies-Bouldin and Ray-Turi CVIs showed better performance among a group of twenty four CVIs in determining the number of clusters to use in dFC analysis.

A resting-state network comparison of combat-related PTSD with combat-exposed and civilian controls

Resting-state functional connectivity (rsFC) is an emerging means of understanding the neurobiology of combat-related post-traumatic stress disorder (PTSD). However, most rsFC studies to date have limited focus to cognitively related intrinsic connectivity networks (ICNs), have not applied data-driven methodologies or have disregarded the effect of combat exposure. In this study, we predicted that group independent component analysis (GICA) would reveal group-wise differences in rsFC across 50 active duty service members with PTSD, 28 combat-exposed controls (CEC), and 25 civilian controls without trauma exposure (CC). Intranetwork connectivity differences were identified across 11 ICNs, yet combat-exposed groups were indistinguishable in PTSD vs CEC contrasts. Both PTSD and CEC demonstrated anatomically diffuse differences in the Auditory Vigilance and Sensorimotor networks compared to CC. However, intranetwork connectivity in a subset of three regions was associated with PTSD symptom severity among executive (left insula; ventral anterior cingulate) and right Fronto-Parietal (perigenual cingulate) networks. Furthermore, we found that increased temporal synchronization among visuospatial and sensorimotor networks was associated with worse avoidance symptoms in PTSD. Longitudinal neuroimaging studies in combat-exposed cohorts can further parse PTSD-related, combat stress-related or adaptive rsFC changes ensuing from combat.

Decreased Cross-Domain Mutual Information in Schizophrenia From Dynamic Connectivity States

The study of dynamic functional network connectivity (dFNC) has been important to understand the healthy and diseased brain. Recent developments model groups of functionally related brain structures (defined as functional domains) as entities that can send and receive information. A domain analysis starts by detecting a finite set of connectivity patterns known as domain states within each functional domain. Dynamic functional domain connectivity (DFDC) is a novel information theoretic framework for studying the temporal sequence of the domain states and the amount of information shared among domains. In this setting, the information flow among functional domains can be compared to the flow of bits among entities in a digital network. Schizophrenia is a chronic psychiatric disorder which is associated with how the brain processes information. Here, we employed the DFDC framework to analyze a dataset containing resting-state fMRI scans from 163 healthy controls (HCs) and 151 schizophrenia patients (SZs). As in other information theory methods, this study measured domain state probabilities, entropy within each DFDC and the cross-domain mutual information (CDMI) between pairs of DFDC. Results indicate that SZs show significantly higher (transformed) entropy than HCs in subcortical (SC)-SC; default mode network (DMN)-visual (VIS) and frontoparietal (FRN)-VIS DFDCs. SZs also show lower (transformed) CDMI between SC-VIS vs. SC-sensorimotor (SM), attention (ATTN)-VIS vs. ATTN-SM and ATTN-SM vs. ATTN-ATTN DFDC pairs after correcting for multiple comparisons. These results imply that different DFDC pairs function in a more independent manner in SZs compared to HCs. Our findings present evidence of higher uncertainty and randomness in SZ brain function.

Dynamic functional connectivity in schizophrenia and autism spectrum disorder: Convergence, divergence and classification

BACKGROUND

Over the recent years there has been a growing debate regarding the extent and nature of the overlap in neuropathology between schizophrenia (SZ) and autism spectrum disorder (ASD). Dynamic functional network connectivity (dFNC) is a recent analysis method that explores temporal patterns of functional connectivity (FC). We compared resting-state dFNC in SZ, ASD and healthy controls (HC), characterized the associations between temporal patterns and symptoms, and performed a three-way classification analysis based on dFNC indices.

METHODS

Resting-state fMRI was collected from 100 young adults: 33 SZ, 33 ASD, 34 HC. Independent component analysis (ICA) was performed, followed by dFNC analysis (window = 33 s, step = 1TR, k-means clustering). Temporal patterns were compared between groups, correlated with symptoms, and classified via cross-validated three-way discriminant analysis.

RESULTS

Both clinical groups displayed an increased fraction of time (FT) spent in a state of weak, intra-network connectivity [p < .001] and decreased FT in a highly-connected state [p < .001]. SZ further showed decreased number of transitions between states [p < .001], decreased FT in a widely-connected state [p < .001], increased dwell time (DT) in the weakly-connected state [p < .001], and decreased DT in the highly-connected state [p = .001]. Social behavior scores correlated with DT in the widely-connected state in SZ [r = 0.416, p = .043], but not ASD. Classification correctly identified SZ at high rates (81.8%), while ASD and HC at lower rates.

CONCLUSIONS

Results indicate a severe and pervasive pattern of temporal aberrations in SZ (specifically, being "stuck" in a state of weak connectivity), that distinguishes SZ participants from both ASD and HC, and is associated with clinical symptoms.

In search of multimodal brain alterations in Alzheimer's and Binswanger's disease

Structural and functional brain abnormalities have been widely identified in dementia, but with variable replicability and significant overlap. Alzheimer's disease (AD) and Binswanger's disease (BD) share similar symptoms and common brain changes that can confound diagnosis. In this study, we aimed to investigate correlated structural and functional brain changes in AD and BD by combining resting-state functional magnetic resonance imaging (fMRI) and diffusion MRI. A group independent component analysis was first performed on the fMRI data to extract 49 intrinsic connectivity networks (ICNs). Then we conducted a multi-set canonical correlation analysis on three features, functional network connectivity (FNC) between ICNs, fractional anisotropy (FA) and mean diffusivity (MD). Two inter-correlated components show significant group differences. The first component demonstrates distinct brain changes between AD and BD. AD shows increased cerebellar FNC but decreased thalamic and hippocampal FNC. Such FNC alterations are linked to the decreased corpus callosum FA. AD also has increased MD in the frontal and temporal cortex, but BD shows opposite alterations. The second component demonstrates specific brain changes in BD. Increased FNC is mainly between default mode and sensory regions, while decreased FNC is mainly within the default mode domain and related to auditory regions. The FNC changes are associated with FA changes in posterior/middle cingulum cortex and visual cortex and increased MD in thalamus and hippocampus. Our findings provide evidence of linked functional and structural deficits in dementia and suggest that AD and BD have both common and distinct changes in white matter integrity and functional connectivity.

•

This is the first study to explore multi-modalities changes in different dementia.

•

A multimodal fusion method is applied to identify joint components.

•

Brain abnormalities in different modalities are highly correlated.

•

Alzheimer's and Binswanger's disease share similar brain changes.

•

Alzheimer's and Binswanger's disease also have distinct brain changes.

Association between the oral microbiome and brain resting state connectivity in smokers

Recent studies have shown a critical role of the gastrointestinal microbiome in brain and behavior via the complex gut-microbiome-brain axis. However, the influence of the oral microbiome in neurological processes is much less studied, especially in response to the stimuli, such as smoking, within the oral microenvironment. Additionally, given the complex structural and functional networks in brain, our knowledge about the relationship between microbiome and brain function through specific brain circuits is still very limited. In this pilot study, we leveraged next generation sequencing for microbiome and functional neuroimaging technique to enable the delineation of microbiome-brain network links as well as their relationship to cigarette smoking. Thirty smokers and 30 age- and sex-matched nonsmokers were recruited for 16S sequencing of their oral microbial community. Among them, 56 subjects were scanned by resting-state functional magnetic resonance imaging to derive brain functional networks. Statistical analyses were performed to demonstrate the influence of smoking on the oral microbial composition, functional network connectivity, and the associations between microbial shifts and functional network connectivity alternations. Compared to nonsmokers, we found a significant decrease of beta diversity (P = 6 × 10^-3) in smokers and identified several classes (Betaproteobacteria, Spirochaetia, Synergistia, and Mollicutes) with significant alterations in microbial abundance. Pathway analysis on the predicted KEGG pathways shows that the microbiota with altered abundance are mainly involved in pathways related to cell processes, DNA repair, immune system, and neurotransmitters signaling. One brain functional network connectivity component was identified to have a significant difference between smokers and nonsmokers (P = 0.032), mainly including connectivity between brain default network and other task-positive networks. This brain functional component was also significantly associated with smoking related microbiota, suggesting a correlated cross-individual pattern between smoking-induced oral microbiome dysbiosis and brain functional connectivity alternation, possibly involving immunological and neurotransmitter signaling pathways. This work is the first attempt to link oral microbiome and brain functional networks, and provides support for future work in characterizing the role of oral microbiome in mediating smoking effects on brain activity.

Adolescent sex differences in cortico-subcortical functional connectivity during response inhibition

Numerous lines of evidence have shown that cognitive processes engaged during response inhibition tasks are associated with structure and functional integration of regions within fronto-parietal networks. However, while prior studies have started to characterize how intrinsic connectivity during resting state differs between boys and girls, comparatively less is known about how functional connectivity differs between males and females when brain function is exogenously driven by the processing demands of typical Go/No-Go tasks that assess both response inhibition and error processing. The purpose of this study was to characterize adolescent sex differences and possible changes in sexually dimorphic regional functional connectivity across adolescent development in both cortical and subcortical brain connectivity elicited during a visual Go/No-Go task. A total of 130 healthy adolescents (ages 12-25 years) performed a Go/No-Go task during functional magnetic resonance imaging. High model-order group independent component analysis was used to characterize whole-brain network functional connectivity during response inhibition and then a univariate technique used to evaluate differences related to sex and age. As predicted and similar to previously described findings from non-task-driven resting state connectivity studies, functional connectivity sex differences were observed in several subcortical regions, including the amygdala, caudate, thalamus, and cortical regions, including inferior frontal gyrus engaged most strongly during successful response inhibition and/or error processing. Importantly, adolescent boys and girls exhibited different normative profiles of age-related changes in several default mode networks of regions and anterior cingulate cortex. These results suggest that cortical-subcortical functional networks supporting response inhibition operate differently between sexes during adolescence.

Parcellation of the human hippocampus based on gray matter volume covariance: Replicable results on healthy young adults

The hippocampus is a key brain region that participates in a range of cognitive and affective functions, and is involved in the etiopathogenesis of numerous neuropsychiatric disorders. The structural complexity and functional diversity of the hippocampus suggest the existence of structural and functional subdivisions within this structure. For the first time, we parcellated the human hippocampus with two independent data sets, each of which consisted of 198 T1-weighted structural magnetic resonance imaging (sMRI) images of healthy young subjects. The method was based on gray matter volume (GMV) covariance, which was quantified by a bivariate voxel-to-voxel linear correlation approach, as well as a multivariate masked independent component analysis approach. We subsequently interrogated the relationship between the GMV covariance patterns and the functional connectivity patterns of the hippocampal subregions using sMRI and resting-state functional MRI (fMRI) data from the same participants. Seven distinct GMV covariance-based subregions were identified for bilateral hippocampi, with robust reproducibility across the two data sets. We further demonstrated that the structural covariance patterns of the hippocampal subregions had a correspondence with the intrinsic functional connectivity patterns of these subregions. Together, our results provide a topographical configuration of the hippocampus with converging structural and functional support. The resulting subregions may improve our understanding of the hippocampal connectivity and functions at a subregional level, which provides useful parcellations and masks for future neuroscience and clinical research on the structural and/or functional connectivity of the hippocampus.

A method for building a genome-connectome bipartite graph model

It has been widely shown that genomic factors influence both risk for schizophrenia and variation in functional brain connectivity. Moreover, schizophrenia is characterized by disrupted brain connectivity. In this work, we proposed a genome-connectome bipartite graph model to perform imaging genomic analysis. Functional network connectivity (FNC) was estimated after decomposing resting state functional magnetic resonance imaging data from both healthy controls (HC) and patients with schizophrenia (SZ) into spatial brain components using group independent component analysis (G-ICA). Then 83 FNC connections showing a group difference (HC vs SZ) were selected as fMRI nodes, and eighty-one schizophrenia-related single nucleotide polymorphisms (SNPs) were selected as genetic nodes respectively in the bipartite graph. Edges connecting pairs of genetic and fMRI nodes were defined based on the SNP-FNC associations across subjects evaluated by a general linear model. Results show that some SNP nodes in the bipartite graph have a high degree implying they are influential in modulating brain connectivity and may be more strongly associated with the risk of schizophrenia than other SNPs. A bi-clustering analysis detected a cluster with 15 SNPs interacting with 38 FNC connections, most of which were within or between somato-motor and visual brain areas. This suggests that the activity of these brain regions may be related to common SNPs and provides insights into the pathology of schizophrenia. The findings suggest that the SNP-FNC bipartite graph approach is a novel model to investigate genetic influences on functional brain connectivity in mental illness.

Dynamic functional network connectivity in Huntington's disease and its associations with motor and cognitive measures

Dynamic functional network connectivity (dFNC) is an expansion of traditional, static FNC that measures connectivity variation among brain networks throughout scan duration. We used a large resting-state fMRI (rs-fMRI) sample from the PREDICT-HD study (N = 183 Huntington disease gene mutation carriers [HDgmc] and N = 78 healthy control [HC] participants) to examine whole-brain dFNC and its associations with CAG repeat length as well as the product of scaled CAG length and age, a variable representing disease burden. We also tested for relationships between functional connectivity and motor and cognitive measurements. Group independent component analysis was applied to rs-fMRI data to obtain whole-brain resting state networks. FNC was defined as the correlation between RSN time-courses. Dynamic FNC behavior was captured using a sliding time window approach, and FNC results from each window were assigned to four clusters representing FNC states, using a k-means clustering algorithm. HDgmc individuals spent significantly more time in State-1 (the state with the weakest FNC pattern) compared to HC. However, overall HC individuals showed more FNC dynamism than HDgmc. Significant associations between FNC states and genetic and clinical variables were also identified. In FNC State-4 (the one that most resembled static FNC), HDgmc exhibited significantly decreased connectivity between the putamen and medial prefrontal cortex compared to HC, and this was significantly associated with cognitive performance. In FNC State-1, disease burden in HDgmc participants was significantly associated with connectivity between the postcentral gyrus and posterior cingulate cortex, as well as between the inferior occipital gyrus and posterior parietal cortex.

The spatial chronnectome reveals a dynamic interplay between functional segregation and integration

The brain is highly dynamic, reorganizing its activity at different interacting spatial and temporal scales, including variation within and between brain networks. The chronnectome is a model of the brain in which nodal activity and connectivity patterns change in fundamental and recurring ways over time. Most literature assumes fixed spatial nodes/networks, ignoring the possibility that spatial nodes/networks may vary in time. Here, we introduce an approach to calculate a spatially fluid chronnectome (called the spatial chronnectome for clarity), which focuses on the variations of networks coupling at the voxel level, and identify a novel set of spatially dynamic features. Results reveal transient spatially fluid interactions between intra- and internetwork relationships in which brain networks transiently merge and separate, emphasizing dynamic segregation and integration. Brain networks also exhibit distinct spatial patterns with unique temporal characteristics, potentially explaining a broad spectrum of inconsistencies in previous studies that assumed static networks. Moreover, we show anticorrelative connections to brain networks are transient as opposed to constant across the entire scan. Preliminary assessments using a multi-site dataset reveal the ability of the approach to obtain new information and nuanced alterations that remain undetected during static analysis. Patients with schizophrenia (SZ) display transient decreases in voxel-wise network coupling within visual and auditory networks, and higher intradomain coupling variability. In summary, the spatial chronnectome represents a new direction of research enabling the study of functional networks which are transient at the voxel level, and the identification of mechanisms for within- and between-subject spatial variability.

Group ICA for identifying biomarkers in schizophrenia: 'Adaptive' networks via spatially constrained ICA show more sensitivity to group differences than spatio-temporal regression

Brain functional networks identified from fMRI data can provide potential biomarkers for brain disorders. Group independent component analysis (GICA) is popular for extracting brain functional networks from multiple subjects. In GICA, different strategies exist for reconstructing subject-specific networks from the group-level networks. However, it is unknown whether these strategies have different sensitivities to group differences and abilities in distinguishing patients. Among GICA, spatio-temporal regression (STR) and spatially constrained ICA approaches such as group information guided ICA (GIG-ICA) can be used to propagate components (indicating networks) to a new subject that is not included in the original subjects. In this study, based on the same a priori network maps, we reconstructed subject-specific networks using these two methods separately from resting-state fMRI data of 151 schizophrenia patients (SZs) and 163 healthy controls (HCs). We investigated group differences in the estimated functional networks and the functional network connectivity (FNC) obtained by each method. The networks were also used as features in a cross-validated support vector machine (SVM) for classifying SZs and HCs. We selected features using different strategies to provide a comprehensive comparison between the two methods. GIG-ICA generally showed greater sensitivity in statistical analysis and better classification performance (accuracy 76.45 ± 8.9%, sensitivity 0.74 ± 0.11, specificity 0.79 ± 0.11) than STR (accuracy 67.45 ± 8.13%, sensitivity 0.65 ± 0.11, specificity 0.71 ± 0.11). Importantly, results were also consistent when applied to an independent dataset including 82 HCs and 82 SZs. Our work suggests that the functional networks estimated by GIG-ICA are more sensitive to group differences, and GIG-ICA is promising for identifying image-derived biomarkers of brain disease.

•

We investigate which group ICA method is more sensitive to group differences in networks and FNC.

•

We evaluate which group ICA method can result in better classification performance.

•

We compare the spatio-temporal regression (STR) and GIG-ICA using different feature selection and model building strategies and independent samples.

•

We use different feature selection and model building strategies and classification using independent samples.

•

Using fMRI data including schizophrenia patients, GIG-ICA shows better performance than the traditional STR method.

Resting state connectivity differences in eyes open versus eyes closed conditions

Functional magnetic resonance imaging data are commonly collected during the resting state. Resting state functional magnetic resonance imaging (rs-fMRI) is very practical and applicable for a wide range of study populations. Rs-fMRI is usually collected in at least one of three different conditions/tasks, eyes closed (EC), eyes open (EO), or eyes fixated on an object (EO-F). Several studies have shown that there are significant condition-related differences in the acquired data. In this study, we compared the functional network connectivity (FNC) differences assessed via group independent component analysis on a large rs-fMRI dataset collected in both EC and EO-F conditions, and also investigated the effect of covariates (e.g., age, gender, and social status score). Our results indicated that task condition significantly affected a wide range of networks; connectivity of visual networks to themselves and other networks was increased during EO-F, while EC was associated with increased connectivity of auditory and sensorimotor networks to other networks. In addition, the association of FNC with age, gender, and social status was observed to be significant only in the EO-F condition (though limited as well). However, statistical analysis did not reveal any significant effect of interaction between eyes status and covariates. These results indicate that resting-state condition is an important variable that may limit the generalizability of clinical findings using rs-fMRI.

Concurrent Cross-Sectional and Longitudinal Analyses of Multivariate White Matter Profiles and Clinical Functioning in Pre-Diagnosis Huntington Disease

BACKGROUND

Gray matter (GM) atrophy in the striatum and across the brain is a consistently reported feature of the Huntington Disease (HD) prodrome. More recently, widespread prodromal white matter (WM) degradation has also been detected. However, longitudinal WM studies are limited and conflicting, and most analyses comparing WM and clinical functioning have also been cross-sectional.

OBJECTIVE

We simultaneously assessed changes in WM and cognitive and motor functioning at various prodromal HD stages.

METHODS

Data from 1,336 (1,047 prodromal, 289 control) PREDICT-HD participants were analyzed (3,700 sessions). MRI images were used to create GM, WM, and cerebrospinal fluid probability maps. Using source-based morphometry, independent component analysis was applied to WM probability maps to extract covarying spatial patterns and their subject profiles. WM profiles were analyzed in two sets of linear mixed model (LMM) analyses: one to compare WM profiles across groups cross-sectionally and longitudinally, and one to concurrently compare WM profiles and clinical variables cross-sectionally and longitudinally within each group.

RESULTS

Findings illustrate widespread prodromal changes in GM-adjacent-WM, with premotor, supplementary motor, middle frontal and striatal changes early in the prodrome that subsequently extend sub-gyrally with progression. Motor functioning agreed most with WM until the near-onset prodromal stage, when Stroop interference was the best WM indicator. Across groups, Trail-Making Test part A outperformed other cognitive variables in its similarity to WM, particularly cross-sectionally.

CONCLUSIONS

Results suggest that distinct regions coincide with cognitive compared to motor functioning. Furthermore, at different prodromal stages, distinct regions appear to align best with clinical functioning. Thus, the informativeness of clinical measures may vary according to the type of data available (cross-sectional or longitudinal) as well as age and CAG-number.

A framework for linking resting-state chronnectome/genome features in schizophrenia: A pilot study

Multimodal, imaging-genomics techniques offer a platform for understanding genetic influences on brain abnormalities in psychiatric disorders. Such approaches utilize the information available from both imaging and genomics data and identify their association. Particularly for complex disorders such as schizophrenia, the relationship between imaging and genomic features may be better understood by incorporating additional information provided by advanced multimodal modeling. In this study, we propose a novel framework to combine features corresponding to functional magnetic resonance imaging (functional) and single nucleotide polymorphism (SNP) data from 61 schizophrenia (SZ) patients and 87 healthy controls (HC). In particular, the features for the functional and genetic modalities include dynamic (i.e., time-varying) functional network connectivity (dFNC) features and the SNP data, respectively. The dFNC features are estimated from component time-courses, obtained using group independent component analysis (ICA), by computing sliding-window functional network connectivity, and then estimating subject specific states from this dFNC data using a k-means clustering approach. For each subject, both the functional (dFNC states) and SNP data are selected as features for a parallel ICA (pICA) based imaging-genomic framework. This analysis identified a significant association between a SNP component (defined by large clusters of functionally related SNPs statistically correlated with phenotype components) and time-varying or dFNC component (defined by clusters of related connectivity links among distant brain regions distributed across discrete dynamic states, and statistically correlated with genomic components) in schizophrenia. Importantly, the polygenetic risk score (PRS) for SZ (computed as a linearly weighted sum of the genotype profiles with weights derived from the odds ratios of the psychiatric genomics consortium (PGC)) was negatively correlated with the significant dFNC component, which were mostly present within a state that exhibited a lower occupancy rate in individuals with SZ compared with HC, hence identifying a potential dFNC imaging biomarker for schizophrenia. Taken together, the current findings provide preliminary evidence for a link between dFNC measures and genetic risk, suggesting the application of dFNC patterns as biomarkers in imaging genetic association study.

Spatial dynamics within and between brain functional domains: A hierarchical approach to study time-varying brain function

The analysis of time-varying activity and connectivity patterns (i.e., the chronnectome) using resting-state magnetic resonance imaging has become an important part of ongoing neuroscience discussions. The majority of previous work has focused on variations of temporal coupling among fixed spatial nodes or transition of the dominant activity/connectivity pattern over time. Here, we introduce an approach to capture spatial dynamics within functional domains (FDs), as well as temporal dynamics within and between FDs. The approach models the brain as a hierarchical functional architecture with different levels of granularity, where lower levels have higher functional homogeneity and less dynamic behavior and higher levels have less homogeneity and more dynamic behavior. First, a high-order spatial independent component analysis is used to approximate functional units. A functional unit is a pattern of regions with very similar functional activity over time. Next, functional units are used to construct FDs. Finally, functional modules (FMs) are calculated from FDs, providing an overall view of brain dynamics. Results highlight the spatial fluidity within FDs, including a broad spectrum of changes in regional associations, from strong coupling to complete decoupling. Moreover, FMs capture the dynamic interplay between FDs. Patients with schizophrenia show transient reductions in functional activity and state connectivity across several FDs, particularly the subcortical domain. Activity and connectivity differences convey unique information in many cases (e.g., the default mode) highlighting their complementarity information. The proposed hierarchical model to capture FD spatiotemporal variations provides new insight into the macroscale chronnectome and identifies changes hidden from existing approaches.

Clustering of fMRI data: the elusive optimal number of clusters

Model-free methods are widely used for the processing of brain fMRI data collected under natural stimulations, sleep, or rest. Among them is the popular fuzzy c-mean algorithm, commonly combined with cluster validity (CV) indices to identify the ‘true’ number of clusters (components), in an unsupervised way. CV indices may however reveal different optimal c-partitions for the same fMRI data, and their effectiveness can be hindered by the high data dimensionality, the limited signal-to-noise ratio, the small proportion of relevant voxels, and the presence of artefacts or outliers. Here, the author investigated the behaviour of seven robust CV indices. A new CV index that incorporates both compactness and separation measures is also introduced. Using both artificial and real fMRI data, the findings highlight the importance of looking at the behavior of different compactness and separation measures, defined here as building blocks of CV indices, to depict a full description of the data structure, in particular when no agreement is found between CV indices. Overall, for fMRI, it makes sense to relax the assumption that only one unique c-partition exists, and appreciate that different c-partitions (with different optimal numbers of clusters) can be useful explanations of the data, given the hierarchical organization of many brain networks.

Examining stability of independent component analysis based on coefficient and component matrices for voxel-based morphometry of structural magnetic resonance imaging

Independent component analysis (ICA) on group-level voxel-based morphometry (VBM) produces the coefficient matrix and the component matrix. The former contains variability among multiple subjects for further statistical analysis, and the latter reveals spatial maps common for all subjects. ICA algorithms converge to local optimization points in practice and the mostly applied stability investigation approach examines the stability of the extracted components. We found that the practically stable components do not guarantee to produce the practically stable coefficients of ICA decomposition for the further statistical analysis. Consequently, we proposed a novel approach including two steps: (1), the stability index for the coefficient matrix and the stability index for the component matrix were examined, respectively; (2) the two indices were multiplied to analyze the stability of ICA decomposition. The proposed approach was used to study the sMRI data of Type II diabetes mellitus group and the healthy control group (HC). Group differences in VBM were found in the superior temporal gyrus. Besides, it was revealed that the VBMs of the region of the HC group were significantly correlated with Montreal Cognitive Assessment (MoCA) describing the level of cognitive disorder. In contrast to the widely applied approach to investigating the stability of the extracted components for ICA decomposition, we proposed to examine the stability of ICA decomposition by fusion the stability of both coefficient matrix and the component matrix. Therefore, the proposed approach can examine the stability of ICA decomposition sufficiently.