Fast, accurate, and interpretable decoding of electrocorticographic signals using dynamic mode decomposition.

Dynamic mode (DM) decomposition decomposes spatiotemporal signals into basic oscillatory components (DMs). DMs can improve the accuracy of neural decoding when used with the nonlinear Grassmann kernel, compared to conventional power features. However, such kernel-based machine learning algorithms have three limitations: large computational time preventing real-time application, incompatibility with non-kernel algorithms, and low interpretability. Here, we propose a mapping function corresponding to the Grassmann kernel that explicitly transforms DMs into spatial DM (sDM) features, which can be used in any machine learning algorithm. Using electrocorticographic signals recorded during various movement and visual perception tasks, the sDM features were shown to improve the decoding accuracy and computational time compared to conventional methods. Furthermore, the components of the sDM features informative for decoding showed similar characteristics to the high-γ power of the signals, but with higher trial-to-trial reproducibility. The proposed sDM features enable fast, accurate, and interpretable neural decoding.

End-to-end deep learning approach to mouse behavior classification from cortex-wide calcium imaging.

Deep learning is a powerful tool for neural decoding, broadly applied to systems neuroscience and clinical studies. Interpretable and transparent models that can explain neural decoding for intended behaviors are crucial to identifying essential features of deep learning decoders in brain activity. In this study, we examine the performance of deep learning to classify mouse behavioral states from mesoscopic cortex-wide calcium imaging data. Our convolutional neural network (CNN)-based end-to-end decoder combined with recurrent neural network (RNN) classifies the behavioral states with high accuracy and robustness to individual differences on temporal scales of sub-seconds. Using the CNN-RNN decoder, we identify that the forelimb and hindlimb areas in the somatosensory cortex significantly contribute to behavioral classification. Our findings imply that the end-to-end approach has the potential to be an interpretable deep learning method with unbiased visualization of critical brain regions.

Manifold alteration between major depressive disorder and healthy control subjects using dynamic mode decomposition in resting-state fMRI data.

Background: The World Health Organization has reported that approximately 300 million individuals suffer from the mood disorder known as MDD. Non-invasive measurement techniques have been utilized to reveal the mechanism of MDD, with rsfMRI being the predominant method. The previous functional connectivity and energy landscape studies have shown the difference in the coactivation patterns between MDD and HCs. However, these studies did not consider oscillatory temporal dynamics. Methods: In this study, the dynamic mode decomposition, a method to compute a set of coherent spatial patterns associated with the oscillation frequency and temporal decay rate, was employed to investigate the alteration of the occurrence of dynamic modes between MDD and HCs. Specifically, The BOLD signals of each subject were transformed into dynamic modes representing coherent spatial patterns and discrete-time eigenvalues to capture temporal variations using dynamic mode decomposition. All the dynamic modes were disentangled into a two-dimensional manifold using t-SNE. Density estimation and density ratio estimation were applied to the two-dimensional manifolds after the two-dimensional manifold was split based on HCs and MDD. Results: The dynamic modes that uniquely emerged in the MDD were not observed. Instead, we have found some dynamic modes that have shown increased or reduced occurrence in MDD compared with HCs. The reduced dynamic modes were associated with the visual and saliency networks while the increased dynamic modes were associated with the default mode and sensory-motor networks. Conclusion: To the best of our knowledge, this study showed initial evidence of the alteration of occurrence of the dynamic modes between MDD and HCs. To deepen understanding of how the alteration of the dynamic modes emerges from the structure, it is vital to investigate the relationship between the dynamic modes, cortical thickness, and surface areas.

Sensor array design of optically pumped magnetometers for accurately estimating source currents.

An optically pumped magnetometer (OPM) is a new generation of magnetoencephalography (MEG) devices that is small, light, and works at room temperature. Due to these characteristics, OPMs enable flexible and wearable MEG systems. On the other hand, if we have a limited number of OPM sensors, we need to carefully design their sensor arrays depending on our purposes and regions of interests (ROIs). In this study, we propose a method that designs OPM sensor arrays for accurately estimating the cortical currents at the ROIs. Based on the resolution matrix of minimum norm estimate (MNE), our method sequentially determines the position of each sensor to optimize its inverse filter pointing to the ROIs and suppressing the signal leakage from the other areas. We call this method the Sensor array Optimization based on Resolution Matrix (SORM). We conducted simple and realistic simulation tests to evaluate its characteristics and efficacy for real OPM-MEG data. SORM designed the sensor arrays so that their leadfield matrices had high effective ranks as well as high sensitivities to ROIs. Although SORM is based on MNE, the sensor arrays designed by SORM were effective not only when we estimated the cortical currents by MNE but also when we did so by other methods. With real OPM-MEG data we confirmed its validity for real data. These analyses suggest that SORM is especially useful when we want to accurately estimate ROIs' activities with a limited number of OPM sensors, such as brain-machine interfaces and diagnosing brain diseases.

Generalizable neuromarker for autism spectrum disorder across imaging sites and developmental stages: A multi-site study.

Autism spectrum disorder (ASD) is a lifelong condition, and its underlying biological mechanisms remain elusive. The complexity of various factors, including inter-site and development-related differences, makes it challenging to develop generalizable neuroimaging-based biomarkers for ASD. This study used a large-scale, multi-site dataset of 730 Japanese adults to develop a generalizable neuromarker for ASD across independent sites and different developmental stages. Our adult ASD neuromarker achieved successful generalization for the US and Belgium adults and Japanese adults. The neuromarker demonstrated significant generalization for children and adolescents. We identified 141 functional connections (FCs) important for discriminating individuals with ASD from TDCs. Finally, we mapped schizophrenia (SCZ) and major depressive disorder (MDD) onto the biological axis defined by the neuromarker and explored the biological continuity of ASD with SCZ and MDD. We observed that SCZ, but not MDD, was located proximate to ASD on the biological dimension defined by the ASD neuromarker. The successful generalization in multifarious datasets and the observed relations of ASD with SCZ on the biological dimensions provide new insights for a deeper understanding of ASD.

Generalizable neuromarker for autism spectrum disorder across imaging sites and developmental stages: A multi-site study.

Autism spectrum disorder (ASD) is a lifelong condition, and its underlying biological mechanisms remain elusive. The complexity of various factors, including inter-site and development-related differences, makes it challenging to develop generalizable neuroimaging-based biomarkers for ASD. This study used a large-scale, multi-site dataset of 730 Japanese adults to develop a generalizable neuromarker for ASD across independent sites (U.S., Belgium, and Japan) and different developmental stages (children and adolescents). Our adult ASD neuromarker achieved successful generalization for the US and Belgium adults (area under the curve [AUC] = 0.70) and Japanese adults (AUC = 0.81). The neuromarker demonstrated significant generalization for children (AUC = 0.66) and adolescents (AUC = 0.71; all P<0.05, family-wise-error corrected). We identified 141 functional connections (FCs) important for discriminating individuals with ASD from TDCs. These FCs largely centered on social brain regions such as the amygdala, hippocampus, dorsomedial and ventromedial prefrontal cortices, and temporal cortices. Finally, we mapped schizophrenia (SCZ) and major depressive disorder (MDD) onto the biological axis defined by the neuromarker and explored the biological continuity of ASD with SCZ and MDD. We observed that SCZ, but not MDD, was located proximate to ASD on the biological dimension defined by the ASD neuromarker. The successful generalization in multifarious datasets and the observed relations of ASD with SCZ on the biological dimensions provide new insights for a deeper understanding of ASD.

Learning task-agnostic and interpretable subsequence-based representation of time series and its applications in fMRI analysis.

The recent success of sequential learning models, such as deep recurrent neural networks, is largely due to their superior representation-learning capability for learning the informative representation of a targeted time series. The learning of these representations is generally goal-directed, resulting in their task-specific nature, giving rise to excellent performance in completing a single downstream task but hindering between-task generalisation. Meanwhile, with increasingly intricate sequential learning models, learned representation becomes abstract to human knowledge and comprehension. Hence, we propose a unified local predictive model based on the multi-task learning paradigm to learn the task-agnostic and interpretable subsequence-based time series representation, allowing versatile use of learned representations in temporal prediction, smoothing, and classification tasks. The targeted interpretable representation could convey the spectral information of the modelled time series to the level of human comprehension. Through a proof-of-concept evaluation study, we demonstrate the empirical superiority of learned task-agnostic and interpretable representation over task-specific and conventional subsequence-based representation, such as symbolic and recurrent learning-based representation, in solving temporal prediction, smoothing, and classification tasks. These learned task-agnostic representations can also reveal the ground-truth periodicity of the modelled time series. We further propose two applications of our unified local predictive model in functional magnetic resonance imaging (fMRI) analysis to reveal the spectral characterisation of cortical areas at rest and reconstruct more smoothed temporal dynamics of cortical activations in both resting-state and task-evoked fMRI data, giving rise to robust decoding.

Distinctive alterations in the mesocorticolimbic circuits in various psychiatric disorders.

AIM: Increasing evidence suggests that psychiatric disorders are linked to alterations in the mesocorticolimbic dopamine-related circuits. However, the common and disease-specific alterations remain to be examined in schizophrenia (SCZ), major depressive disorder (MDD), and autism spectrum disorder (ASD). Thus, this study aimed to examine common and disease-specific features related to mesocorticolimbic circuits. METHODS: This study included 555 participants from four institutes with five scanners: 140 individuals with SCZ (45.0% female), 127 individuals with MDD (44.9%), 119 individuals with ASD (15.1%), and 169 healthy controls (HC) (34.9%). All participants underwent resting-state functional magnetic resonance imaging. A parametric empirical Bayes approach was adopted to compare estimated effective connectivity among groups. Intrinsic effective connectivity focusing on the mesocorticolimbic dopamine-related circuits including the ventral tegmental area (VTA), shell and core parts of the nucleus accumbens (NAc), and medial prefrontal cortex (mPFC) were examined using a dynamic causal modeling analysis across these psychiatric disorders. RESULTS: The excitatory shell-to-core connectivity was greater in all patients than in the HC group. The inhibitory shell-to-VTA and shell-to-mPFC connectivities were greater in the ASD group than in the HC, MDD, and SCZ groups. Furthermore, the VTA-to-core and VTA-to-shell connectivities were excitatory in the ASD group, while those connections were inhibitory in the HC, MDD, and SCZ groups. CONCLUSION: Impaired signaling in the mesocorticolimbic dopamine-related circuits could be an underlying neuropathogenesis of various psychiatric disorders. These findings will improve the understanding of unique neural alternations of each disorder and will facilitate identification of effective therapeutic targets.

Aberrant Large-Scale Network Interactions Across Psychiatric Disorders Revealed by Large-Sample Multi-Site Resting-State Functional Magnetic Resonance Imaging Datasets.

BACKGROUND AND HYPOTHESIS: Dynamics of the distributed sets of functionally synchronized brain regions, known as large-scale networks, are essential for the emotional state and cognitive processes. However, few studies were performed to elucidate the aberrant dynamics across the large-scale networks across multiple psychiatric disorders. In this paper, we aimed to investigate dynamic aspects of the aberrancy of the causal connections among the large-scale networks of the multiple psychiatric disorders. STUDY DESIGN: We applied dynamic causal modeling (DCM) to the large-sample multi-site dataset with 739 participants from 4 imaging sites including 4 different groups, healthy controls, schizophrenia (SCZ), major depressive disorder (MDD), and bipolar disorder (BD), to compare the causal relationships among the large-scale networks, including visual network, somatomotor network (SMN), dorsal attention network (DAN), salience network (SAN), limbic network (LIN), frontoparietal network, and default mode network. STUDY RESULTS: DCM showed that the decreased self-inhibitory connection of LIN was the common aberrant connection pattern across psychiatry disorders. Furthermore, increased causal connections from LIN to multiple networks, aberrant self-inhibitory connections of DAN and SMN, and increased self-inhibitory connection of SAN were disorder-specific patterns for SCZ, MDD, and BD, respectively. CONCLUSIONS: DCM revealed that LIN was the core abnormal network common to psychiatric disorders. Furthermore, DCM showed disorder-specific abnormal patterns of causal connections across the 7 networks. Our findings suggested that aberrant dynamics among the large-scale networks could be a key biomarker for these transdiagnostic psychiatric disorders.

Virtual reality-based real-time imaging reveals abnormal cortical dynamics during behavioral transitions in a mouse model of autism.

Functional connectivity (FC) can provide insight into cortical circuit dysfunction in neuropsychiatric disorders. However, dynamic changes in FC related to locomotion with sensory feedback remain to be elucidated. To investigate FC dynamics in locomoting mice, we develop mesoscopic Ca2+ imaging with a virtual reality (VR) environment. We find rapid reorganization of cortical FC in response to changing behavioral states. By using machine learning classification, behavioral states are accurately decoded. We then use our VR-based imaging system to study cortical FC in a mouse model of autism and find that locomotion states are associated with altered FC dynamics. Furthermore, we identify FC patterns involving the motor area as the most distinguishing features of the autism mice from wild-type mice during behavioral transitions, which might correlate with motor clumsiness in individuals with autism. Our VR-based real-time imaging system provides crucial information to understand FC dynamics linked to behavioral abnormality of neuropsychiatric disorders.

Verification of the brain network marker of major depressive disorder: Test-retest reliability and anterograde generalization performance for newly acquired data.

BACKGROUND: Recently, we developed a generalizable brain network marker for the diagnosis of major depressive disorder (MDD) across multiple imaging sites using resting-state functional magnetic resonance imaging. Here, we applied this brain network marker to newly acquired data to verify its test-retest reliability and anterograde generalization performance for new patients. METHODS: We tested the sensitivity and specificity of our brain network marker of MDD using data acquired from 43 new patients with MDD as well as new data from 33 healthy controls (HCs) who participated in our previous study. To examine the test-retest reliability of our brain network marker, we evaluated the intraclass correlation coefficients (ICCs) between the brain network marker-based classifier's output (probability of MDD) in two sets of HC data obtained at an interval of approximately 1 year. RESULTS: Test-retest correlation between the two sets of the classifier's output (probability of MDD) from HCs exhibited moderate reliability with an ICC of 0.45 (95 % confidence interval,0.13-0.68). The classifier distinguished patients with MDD and HCs with an accuracy of 69.7 % (sensitivity, 72.1 %; specificity, 66.7 %). LIMITATIONS: The data of patients with MDD in this study were cross-sectional, and the clinical significance of the marker, such as whether it is a state or trait marker of MDD and its association with treatment responsiveness, remains unclear. CONCLUSIONS: The results of this study reaffirmed the test-retest reliability and generalization performance of our brain network marker for the diagnosis of MDD.

Neurofeedback Training without Explicit Phantom Hand Movements and Hand-Like Visual Feedback to Modulate Pain: A Randomized Crossover Feasibility Trial.

Phantom limb pain is attributed to abnormal sensorimotor cortical representations, although the causal relationship between phantom limb pain and sensorimotor cortical representations suffers from the potentially confounding effects of phantom hand movements. We developed neurofeedback training to change sensorimotor cortical representations without explicit phantom hand movements or hand-like visual feedback. We tested the feasibility of neurofeedback training in fourteen patients with phantom limb pain. Neurofeedback training was performed in a single-blind, randomized, crossover trial using two decoders constructed using motor cortical currents measured during phantom hand movements; the motor cortical currents contralateral or ipsilateral to the phantom hand (contralateral and ipsilateral training) were estimated from magnetoencephalograms. Patients were instructed to control the size of a disk, which was proportional to the decoding results, but to not move their phantom hands or other body parts. The pain assessed by the visual analogue scale was significantly greater after contralateral training than after ipsilateral training. Classification accuracy of phantom hand movements significantly increased only after contralateral training. These results suggested that the proposed neurofeedback training changed phantom hand representation and modulated pain without explicit phantom hand movements or hand-like visual feedback, thus showing the relation between the phantom hand representations and pain. PERSPECTIVE: Our work demonstrates the feasibility of using neurofeedback training to change phantom hand representation and modulate pain perception without explicit phantom hand movements and hand-like visual feedback. The results enhance the mechanistic understanding of certain treatments, such as mirror therapy, that change the sensorimotor cortical representation.

Diffuse Optical Tomography Using fNIRS Signals Measured from the Skull Surface of the Macaque Monkey.

Diffuse optical tomography (DOT), as a functional near-infrared spectroscopy (fNIRS) technique, can estimate three-dimensional (3D) images of the functional hemodynamic response in brain volume from measured optical signals. In this study, we applied DOT algorithms to the fNIRS data recorded from the surface of macaque monkeys' skulls when the animals performed food retrieval tasks using either the left- or right-hand under head-free conditions. The hemodynamic response images, reconstructed by DOT with a high sampling rate and fine voxel size, demonstrated significant activations at the upper limb regions of the primary motor area in the central sulcus and premotor, and parietal areas contralateral to the hands used in the tasks. The results were also reliable in terms of consistency across different recording dates. Time-series analyses of each brain area revealed preceding activity of premotor area to primary motor area consistent with previous physiological studies. Therefore, the fNIRS-DOT protocol demonstrated in this study provides reliable 3D functional brain images over a period of days under head-free conditions for region-of-interest-based time-series analysis.

Predicting behavior through dynamic modes in resting-state fMRI data.

Dynamic properties of resting-state functional connectivity (FC) provide rich information on brain-behavior relationships. Dynamic mode decomposition (DMD) has been used as a method to characterize FC dynamics. However, it remains unclear whether dynamic modes (DMs), spatial-temporal coherent patterns computed by DMD, provide information about individual behavioral differences. This study established a methodological approach to predict individual differences in behavior using DMs. Furthermore, we investigated the contribution of DMs within each of seven specific frequency bands (0-0.1,...,0.6-0.7 Hz) for prediction. To validate our approach, we tested whether each of 59 behavioral measures could be predicted by performing multivariate pattern analysis on a Gram matrix, which was created using subject-specific DMs computed from resting-state functional magnetic resonance imaging (rs-fMRI) data of individuals. DMD successfully predicted behavior and outperformed temporal and spatial independent component analysis, which is the conventional data decomposition method for extracting spatial activity patterns. Most of the behavioral measures that were predicted with significant accuracy in a permutation test were related to cognition. We found that DMs within frequency bands <0.2 Hz primarily contributed to prediction and had spatial structures similar to several common resting-state networks. Our results indicate that DMD is efficient in extracting spatiotemporal features from rs-fMRI data.

Whole-brain propagating patterns in human resting-state brain activities.

Repetitive propagating activities in resting-state brain activities have been widely observed in various species and regions. Because they resemble the preceding brain activities during tasks, they are assumed to reflect past experiences embedded in neuronal circuits. "Whole-brain" propagating activities may also reflect a process that integrates information distributed over the entire brain, such as visual and motor information. Here we reveal whole-brain propagating activities from human resting-state magnetoencephalography (MEG) and electroencephalography (EEG) data. We simultaneously recorded the MEGs and EEGs and estimated the source currents from both measurements. Then using our recently proposed algorithm, we extracted repetitive spatiotemporal patterns from the source currents. The estimated patterns consisted of multiple frequency components, each of which transiently exhibited the frequency-specific resting-state networks (RSNs) of functional MRIs (fMRIs), such as the default mode and sensorimotor networks. A simulation test suggested that the spatiotemporal patterns reflected the phase alignment of the multiple frequency oscillators induced by the propagating activities along the anatomical connectivity. These results argue that whole-brain propagating activities transiently exhibited multiple RSNs in their multiple frequency components, suggesting that they reflected a process to integrate the information distributed over the frequencies and networks.

Clustering of Multiple Psychiatric Disorders Using Functional Connectivity in the Data-Driven Brain Subnetwork.

Recently, the dimensional approach has attracted much attention, bringing a paradigm shift to a continuum of understanding of different psychiatric disorders. In line with this new paradigm, we examined whether there was common functional connectivity related to various psychiatric disorders in an unsupervised manner without explicitly using diagnostic label information. To this end, we uniquely applied a newly developed network-based multiple clustering method to resting-state functional connectivity data, which allowed us to identify pairs of relevant brain subnetworks and subject cluster solutions accordingly. Thus, we identified four subject clusters, which were characterized as major depressive disorder (MDD), young healthy control (young HC), schizophrenia (SCZ)/bipolar disorder (BD), and autism spectrum disorder (ASD), respectively, with the relevant brain subnetwork represented by the cerebellum-thalamus-pallidum-temporal circuit. The clustering results were validated using independent datasets. This study is the first cross-disorder analysis in the framework of unsupervised learning of functional connectivity based on a data-driven brain subnetwork.

A multi-site, multi-disorder resting-state magnetic resonance image database.

Machine learning classifiers for psychiatric disorders using resting-state functional magnetic resonance imaging (rs-fMRI) have recently attracted attention as a method for directly examining relationships between neural circuits and psychiatric disorders. To develop accurate and generalizable classifiers, we compiled a large-scale, multi-site, multi-disorder neuroimaging database. The database comprises resting-state fMRI and structural images of the brain from 993 patients and 1,421 healthy individuals, as well as demographic information such as age, sex, and clinical rating scales. To harmonize the multi-site data, nine healthy participants ("traveling subjects") visited the sites from which the above datasets were obtained and underwent neuroimaging with 12 scanners. All participants consented to having their data shared and analyzed at multiple medical and research institutions participating in the project, and 706 patients and 1,122 healthy individuals consented to having their data disclosed. Finally, we have published four datasets: 1) the SRPBS Multi-disorder Connectivity Dataset 2), the SRPBS Multi-disorder MRI Dataset (restricted), 3) the SRPBS Multi-disorder MRI Dataset (unrestricted), and 4) the SRPBS Traveling Subject MRI Dataset.

Comparison of traveling-subject and ComBat harmonization methods for assessing structural brain characteristics.

Multisite magnetic resonance imaging (MRI) is increasingly used in clinical research and development. Measurement biases-caused by site differences in scanner/image-acquisition protocols-negatively influence the reliability and reproducibility of image-analysis methods. Harmonization can reduce bias and improve the reproducibility of multisite datasets. Herein, a traveling-subject (TS) dataset including 56 T1-weighted MRI scans of 20 healthy participants in three different MRI procedures-20, 19, and 17 subjects in Procedures 1, 2, and 3, respectively-was considered to compare the reproducibility of TS-GLM, ComBat, and TS-ComBat harmonization methods. The minimum participant count required for harmonization was determined, and the Cohen's d between different MRI procedures was evaluated as a measurement-bias indicator. The measurement-bias reduction realized with different methods was evaluated by comparing test-retest scans for 20 healthy participants. Moreover, the minimum subject count for harmonization was determined by comparing test-retest datasets. The results revealed that TS-GLM and TS-ComBat reduced measurement bias by up to 85 and 81.3%, respectively. Meanwhile, ComBat showed a reduction of only 59.0%. At least 6 TSs were required to harmonize data obtained from different MRI scanners, complying with the imaging protocol predetermined for multisite investigations and operated with similar scan parameters. The results indicate that TS-based harmonization outperforms ComBat for measurement-bias reduction and is optimal for MRI data in well-prepared multisite investigations. One drawback is the small sample size used, potentially limiting the applicability of ComBat. Investigation on the number of subjects needed for a large-scale study is an interesting future problem.

Common Brain Networks Between Major Depressive-Disorder Diagnosis and Symptoms of Depression That Are Validated for Independent Cohorts.

Large-scale neuroimaging data acquired and shared by multiple institutions are essential to advance neuroscientific understanding of pathophysiological mechanisms in psychiatric disorders, such as major depressive disorder (MDD). About 75% of studies that have applied machine learning technique to neuroimaging have been based on diagnoses by clinicians. However, an increasing number of studies have highlighted the difficulty in finding a clear association between existing clinical diagnostic categories and neurobiological abnormalities. Here, using resting-state functional magnetic resonance imaging, we determined and validated resting-state functional connectivity related to depression symptoms that were thought to be directly related to neurobiological abnormalities. We then compared the resting-state functional connectivity related to depression symptoms with that related to depression diagnosis that we recently identified. In particular, for the discovery dataset with 477 participants from 4 imaging sites, we removed site differences using our recently developed harmonization method and developed a brain network prediction model of depression symptoms (Beck Depression Inventory-II [BDI] score). The prediction model significantly predicted BDI score for an independent validation dataset with 439 participants from 4 different imaging sites. Finally, we found 3 common functional connections between those related to depression symptoms and those related to MDD diagnosis. These findings contribute to a deeper understanding of the neural circuitry of depressive symptoms in MDD, a hetero-symptomatic population, revealing the neural basis of MDD.

Multiple clustering for identifying subject clusters and brain sub-networks using functional connectivity matrices without vectorization.

In neuroscience, the functional magnetic resonance imaging (fMRI) is a vital tool to non-invasively access brain activity. Using fMRI, the functional connectivity (FC) between brain regions can be inferred, which has contributed to a number of findings of the fundamental properties of the brain. As an important clinical application of FC, clustering of subjects based on FC recently draws much attention, which can potentially reveal important heterogeneity in subjects such as subtypes of psychiatric disorders. In particular, a multiple clustering method is a powerful analytical tool, which identifies clustering patterns of subjects depending on their FC in specific brain areas. However, when one applies an existing multiple clustering method to fMRI data, there is a need to simplify the data structure, independently dealing with elements in a FC matrix, i.e., vectorizing a correlation matrix. Such a simplification may distort the clustering results. To overcome this problem, we propose a novel multiple clustering method based on Wishart mixture models, which preserves the correlation matrix structure without vectorization. The uniqueness of this method is that the multiple clustering of subjects is based on particular networks of nodes (or regions of interest, ROIs), optimized in a data-driven manner. Hence, it can identify multiple underlying pairs of associations between a subject cluster solution and a ROI sub-network. The key assumption of the method is independence among sub-networks, which is effectively addressed by whitening correlation matrices. We applied the proposed method to synthetic and fMRI data, demonstrating the usefulness and power of the proposed method.