Search

Identifying relationships among multiple datasets is an effective way to summarize information and has been growing in importance. In this paper, we propose a robust 3-step method for identifying the relationship structure among multiple datasets based on Independent Vector Analysis (IVA) and bootstrap-based hypothesis testing. Unlike previous approaches, our theory-backed method eliminates the need for user-defined thresholds and can effectively handle non-Gaussian data. It achieves this by incorporating higher-order statistics through IVA and employing an eigenvalue decomposition-based feature extraction approach without distributional assumptions. This way, our method estimates more interpretable components and effectively identifies the relationship structure using hierarchical clustering. Simulation results demonstrate the effectiveness of our method, as it achieves perfect Adjusted Mutual Information (AMI) for different values of the correlation between the components. When applied to multi-task fMRI data from patients with schizophrenia and healthy controls, our method successfully reveals activated brain regions associated with the disorder, and identifies the relationship structure of task datasets that matches our prior knowledge of the experiment. Moreover, our proposed method extends beyond task datasets, offering broad applicability in subgroup identification in neuroimaging and other domains.

Revealing associations among various structural and functional patterns of the brain can yield highly informative results about the healthy and disordered brain. Studies using neuroimaging data have more recently begun to utilize the information within as well as across various functional and anatomical domains (i.e., groups of brain networks). However, most whole-brain approaches assume similar complexity of interactions throughout the brain. Here we investigate the hypothesis that interactions between brain networks capture varying amounts of complexity, and that we can better capture this information by varying the complexity of the model subspace structure based on available training data. To do this, we employ a Bayesian optimization-based framework known as the Tree Parzen Estimator (TPE) to identify, exploit and analyze patterns of variation in the information encoded by temporal information extracted from functional magnetic resonance imaging (fMRI) subdomains of the brain. Using a repeated cross-validation procedure on a schizophrenia classification task, we demonstrate evidence that interactions between specific functional subdomains are better characterized by more sophisticated model architectures compared to less complicated ones required by the others for optimally contributing towards classification and understanding the brain's functional interactions. We show that functional subdomains known to be involved in schizophrenia require more complex architectures to optimally unravel discriminatory information about the disorder. Our study points to the need for adaptive, hierarchical learning frameworks that cater differently to the features from different subdomains, not only for a better prediction but also for enabling the identification of features predicting the outcome of interest.

BackgroundA number of studies in recent years have explored whole-brain dynamic connectivity using pairwise approaches. There has been less focus on trying to analyze brain dynamics in higher dimensions over time.MethodsWe introduce a new approach that analyzes time series trajectories to identify high traffic nodes in a high dimensional space. First, functional magnetic resonance imaging (fMRI) data are decomposed using spatial ICA to a set of maps and their associated time series. Next, density is calculated for each time point and high-density points are clustered to identify a small set of high traffic nodes. We validated our method using simulations and then implemented it on a real data set.ResultsWe present a novel approach that captures dynamics within a high dimensional space and also does not use any windowing in contrast to many existing approaches. The approach enables one to characterize and study the time series in a potentially high dimensional space, rather than looking at each component pair separately. Our results show that schizophrenia patients have a lower dynamism compared to healthy controls. In addition, we find patients spend more time in nodes associated with the default mode network and less time in components strongly correlated with auditory and sensorimotor regions. Interestingly, we also found that subjects oscillate between state pairs that show opposite spatial maps, suggesting an oscillatory pattern.ConclusionOur proposed method provides a novel approach to analyze the data in its native high dimensional space and can possibly provide new information that is undetectable using other methods.

ObjectiveMultimodal measurements of the same phenomena provide complementary information and highlight different perspectives, albeit each with their own limitations. A focus on a single modality may lead to incorrect inferences, which is especially important when a studied phenomenon is a disease. In this paper, we introduce a method that takes advantage of multimodal data in addressing the hypotheses of disconnectivity and dysfunction within schizophrenia (SZ).MethodsWe start with estimating and visualizing links within and among extracted multimodal data features using a Gaussian graphical model (GGM). We then propose a modularity-based method that can be applied to the GGM to identify links that are associated with mental illness across a multimodal data set. Through simulation and real data, we show our approach reveals important information about disease-related network disruptions that are missed with a focus on a single modality. We use functional MRI (fMRI), diffusion MRI (dMRI), and structural MRI (sMRI) to compute the fractional amplitude of low frequency fluctuations (fALFF), fractional anisotropy (FA), and gray matter (GM) concentration maps. These three modalities are analyzed using our modularity method.ResultsOur results show missing links that are only captured by the cross-modal information that may play an important role in disconnectivity between the components.ConclusionWe identified multimodal (fALFF, FA and GM) disconnectivity in the default mode network area in patients with SZ, which would not have been detectable in a single modality.SignificanceThe proposed approach provides an important new tool for capturing information that is distributed among multiple imaging modalities.

Network neuroscience, a multidisciplinary field merging insights from neuroscience and network theory, offers a profound understanding of neural network intricacies. However, the impact of varying node sizes on computed graph metrics in neuroimaging data remains underexplored. This study addresses this gap by adopting a data-driven methodology to delineate functional nodes and assess their influence on graph metrics. Using the Neuromark framework, automated independent component analysis is applied to resting state fMRI data, capturing functional network connectivity (FNC) matrices. Global and local graph metrics reveal intricate connectivity patterns, emphasizing the need for nuanced analysis. Notably, node sizes, computed based on voxel counts, contribute to a novel metric termed ‘node-metric coupling’ (NMC). Correlations between graph metrics and node dimensions are consistently observed. The study extends its analysis to a dataset comprising Alzheimer’s disease, mild cognitive impairment, and control subjects, showcasing the potential of NMC as a biomarker for brain disorders. The two key outcomes underscore the interplay between node sizes and resultant graph metrics within a given atlas, shedding light on an often-overlooked source of variability. Additionally, the study highlights the utility of NMC as a valuable biomarker, emphasizing the necessity of accounting for node sizes in future neuroimaging investigations. This work contributes to refining comparative studies employing diverse atlases and advocates for thoughtful consideration of intra-atlas node size in shaping graph metrics, paving the way for more robust neuroimaging research.

Most neuroimaging studies display results that represent only a tiny fraction of the collected data. While it is conventional to present ''only the significant results'' to the reader, here we suggest that this practice has several negative consequences for both reproducibility and understanding. This practice hides away most of the results of the dataset and leads to problems of selection bias and irreproducibility, both of which have been recognized as major issues in neuroimaging studies recently. Opaque, all-or-nothing thresholding, even if well-intentioned, places undue influence on arbitrary filter values, hinders clear communication of scientific results, wastes data, is antithetical to good scientific practice, and leads to conceptual inconsistencies. It is also inconsistent with the properties of the acquired data and the underlying biology being studied. Instead of presenting only a few statistically significant locations and hiding away the remaining results, studies should ''highlight'' the former while also showing as much as possible of the rest. This is distinct from but complementary to utilizing data sharing repositories: the initial presentation of results has an enormous impact on the interpretation of a study. We present practical examples and extensions of this approach for voxelwise, regionwise and cross-study analyses using publicly available data that was analyzed previously by 70 teams (NARPS; Botvinik-Nezer, et al., 2020), showing that it is possible to balance the goals of displaying a full set of results with providing the reader reasonably concise and ''digestible'' findings. In particular, the highlighting approach sheds useful light on the kind of variability present among the NARPS teams' results, which is primarily a varied strength of agreement rather than disagreement. Using a meta-analysis built on the informative ''highlighting'' approach shows this relative agreement, while one using the standard ''hiding'' approach does not. We describe how this simple but powerful change in practice—focusing on highlighting results, rather than hiding all but the strongest ones—can help address many large concerns within the field, or at least to provide more complete information about them. We include a list of practical suggestions for results reporting to improve reproducibility, cross-study comparisons and meta-analyses.

The brain is a complex, multiscale dynamical system composed of many interacting regions. Knowledge of the spatiotemporal organization of these interactions is critical for establishing a solid understanding of the brain's functional architecture and the relationship between neural dynamics and cognition in health and disease. The possibility of studying these dynamics through careful analysis of neuroimaging data has catalyzed substantial interest in methods that estimate time-resolved fluctuations in functional connectivity (often referred to as "dynamic" or time-varying functional connectivity; TVFC). At the same time, debates have emerged regarding the application of TVFC analyses to resting fMRI data, and about the statistical validity, physiological origins, and cognitive and behavioral relevance of resting TVFC. These and other unresolved issues complicate interpretation of resting TVFC findings and limit the insights that can be gained from this promising new research area. This article brings together scientists with a variety of perspectives on resting TVFC to review the current literature in light of these issues. We introduce core concepts, define key terms, summarize controversies and open questions, and present a forward-looking perspective on how resting TVFC analyses can be rigorously and productively applied to investigate a wide range of questions in cognitive and systems neuroscience.

BackgroundDynamic functional network connectivity (dFNC) of the brain has attracted considerable attention recently. Many approaches have been suggested to study dFNC with sliding window Pearson correlation (SWPC) being the most well-known. SWPC needs a relatively large sample size to reach a robust estimation but using large window sizes prevents us to detect rapid changes in dFNC.New methodHere we first calculate the gradients of each time series pair and use the magnitude of these gradients to calculate weighted average of shared trajectory (WAST) as a new estimator for dFNC.ResultsUsing WAST to compare healthy control and schizophrenia patients using a large dataset, we show disconnectivity between different regions associated with schizophrenia. In addition, WAST results reveals patients with schizophrenia stay longer in a connectivity state with negative connectivity between motor and sensory regions than do healthy controls.Comparison with existing methodsWe compare WAST with SWPC and multiplication of temporal derivatives (MTD) using different simulation scenarios. We show that WAST enables us to detect very rapid changes in dFNC (undetected by SWPC) while MTD performance is generally lower.ConclusionsAs large window sizes are unable to detect short states, using shorter window size is desirable if the estimator is robust enough. We provide evidence that WAST requires fewer samples (compared to SWPC) to reach a robust estimation. As a result, we were able to identify rapidly varying dFNC patterns undetected by SWPC while still being able to robustly estimate slower dFNC patterns.

The default mode network (DMN) plays a crucial role in internal self-processing, rumination, and social functions. Disruptions to DMN connectivity have been linked with early adversity and the emergence of psychopathology in adolescence and early adulthood. Herein, we investigate how subclinical psychiatric symptoms can impact DMN functional connectivity during the pubertal transition. Resting-state fMRI data were collected annually from 190 typically-developing youth (9–15 years-old) at three timepoints and within-network DMN connectivity was computed. We used latent growth curve modeling to determine how self-reported depressive and posttraumatic stress symptoms predicted rates of change in DMN connectivity over the three-year period. In the baseline model without predictors, we found no systematic changes in DMN connectivity over time. However, significant modulation emerged after adding psychopathology predictors; greater depressive symptomatology was associated with significant decreases in connectivity over time, whereas posttraumatic stress symptoms were associated with significant increases in connectivity over time. Follow-up analyses revealed that these effects were driven by connectivity changes involving the dorsal medial prefrontal cortex subnetwork. In conclusion, these data suggest that subclinical depressive and posttraumatic symptoms alter the trajectory of DMN connectivity, which may indicate that this network is a nexus of clinical significance in mental health disorders.

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.

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.

BackgroundThe N400 event-related potential (ERP) is triggered by meaningful stimuli that are incongruous, or unmatched, with their semantic context. Functional magnetic resonance imaging (fMRI) studies have identified brain regions activated by semantic incongruity, but their precise links to the N400 ERP are unclear. In schizophrenia (SZ), N400 amplitude reduction is thought to reflect overly broad associations in semantic networks, but the abnormalities in brain networks underlying deficient N400 remain unknown. We utilized joint independent component analysis (JICA) to link temporal patterns in ERPs to neuroanatomical patterns from fMRI and investigate relationships between N400 amplitude and neuroanatomical activation in SZ patients and healthy controls (HC).MethodsSZ patients (n = 24) and HC participants (n = 25) performed a picture-word matching task, in which words were either matched (APPLE→apple) by preceding pictures, or were unmatched by semantically related (in-category; IC, APPLE→lemon) or unrelated (out of category; OC, APPLE→cow) pictures, in separate ERP and fMRI sessions. A JICA "data fusion" analysis was conducted to identify the fMRI brain regions specifically associated with the ERP N400 component. SZ and HC loading weights were compared and correlations with clinical symptoms were assessed.ResultsJICA identified an ERP-fMRI "fused" component that captured the N400, with loading weights that were reduced in SZ. The JICA map for the IC condition showed peaks of activation in the cingulate, precuneus, bilateral temporal poles and cerebellum, whereas the JICA map from the OC condition was linked primarily to visual cortical activation and the left temporal pole. Among SZ patients, fMRI activity from the IC condition was inversely correlated with unusual thought content.ConclusionsThe neural networks associated with the N400 ERP response to semantic violations depends on conceptual relatedness. These findings are consistent with a distributed network underlying neural responses to semantic incongruity including unimodal visual areas as well as integrative, transmodal areas. Unusual thoughts in SZ may reflect impaired processing in transmodal hub regions such as the precuneus, leading to overly broad semantic associations.

Individuals at clinical high-risk (CHR) for psychosis are characterized by attenuated psychotic symptoms. Only a minority of CHR individuals convert to full-blown psychosis. Therefore, there is a strong interest in identifying neurobiological abnormalities underlying the psychosis risk syndrome. Dynamic functional connectivity (DFC) captures time-varying connectivity over short time scales, and has the potential to reveal complex brain functional organization. Based on resting-state functional magnetic resonance imaging (fMRI) data from 70 healthy controls (HCs), 53 CHR individuals, and 58 early illness schizophrenia (ESZ) patients, we applied a novel group information guided ICA (GIG-ICA) to estimate inherent connectivity states from DFC, and then investigated group differences. We found that ESZ patients showed more aberrant connectivities and greater alterations than CHR individuals. Results also suggested that disease-related connectivity states occurred in CHR and ESZ groups. Regarding the dominant state with the highest contribution to dynamic connectivity, ESZ patients exhibited greater impairments than CHR individuals primarily in the cerebellum, frontal cortex, thalamus and temporal cortex, while CHR and ESZ populations shared common aberrances mainly in the supplementary motor area, parahippocampal gyrus and postcentral cortex. CHR-specific changes were also found in the connections between the superior frontal gyrus and calcarine cortex in the dominant state. Our findings suggest that CHR individuals generally show an intermediate functional connectivity pattern between HCs and SZ patients but also have unique connectivity alterations.

Hallucinations characterize schizophrenia, with approximately 59% of patients reporting auditory hallucinations and 27% reporting visual hallucinations. Prior neuroimaging studies suggest that hallucinations are linked to disrupted communication across distributed (sensory, salience-monitoring and subcortical) networks. Yet, our understanding of the neurophysiological mechanisms that underlie auditory and visual hallucinations in schizophrenia remains limited. This study integrates two resting-state functional magnetic resonance imaging (fMRI) analysis methods - amplitudes of low-frequency fluctuations (ALFF) and functional network connectivity (FNC) - to explore the hypotheses that (1) abnormal FNC between salience and sensory (visual/auditory) networks underlies hallucinations in schizophrenia, and (2) disrupted hippocampal oscillations (as measured by hippocampal ALFF) beget changes in FNC linked to hallucinations. Our first hypothesis was supported by the finding that schizophrenia patients reporting hallucinations have higher FNC between the salience network and an associative auditory network relative to healthy controls. Hippocampal ALFF was negatively associated with FNC between primary auditory cortex and the salience network in healthy subjects, but was positively associated with FNC between these networks in patients reporting hallucinations. These findings provide indirect support favoring our second hypothesis. We suggest future studies integrate fMRI with electroencephalogram (EEG) and/or magnetoencephalogram (MEG) methods to directly probe the temporal relation between altered hippocampal oscillations and changes in cross-network functional communication.

We provided the dataset of brain connectome matrices, their similarities measures to self and others longitudinally, and Kessler's psychological distress scales (K10) including the response to each question. The dataset can be used to replicate the results of the manuscript titled “A longitudinal study of functional connectome uniqueness and its association with psychological distress in adolescence”. The functional connectome (whole-brain and 13 networks) matrices were calculated from the resting-state functional MRIs (rs-fMRIs). We collected rs-fMRI and Kessler's psychological distress scale (K10) in 77 adolescents longitudinally up to 9 times from 12 years of age every four months. After removal of data with excessive motion, 262 functional connectome matrices were provided with this paper. The 300 regions of interest (ROIs) were defined using the Greene lab brain atlas. The functional connectome matrices were calculated as correlations between time series from any pair of ROIs extracted from pre-processed fMRIs. This dataset could be potentially used to 1. Understand developmental changes in the functional brain connectivity, 2. As a normal control database of functional connectome matrices, 3. Develop and validate connectome and network-related analysing methods.

Alcohol use disorder (AUD) is a burden to society creating social and health problems. Detection of AUD and its effects on the brain are difficult to assess. This problem is enhanced by the comorbid use of other substances such as nicotine that has been present in previous studies. Recent machine learning algorithms have raised the attention of researchers as a useful tool in studying and detecting AUD. This work uses AUD and controls samples free of any other substance use to assess the performance of a set of commonly used machine learning classifiers detecting AUD from resting state functional network connectivity (rsFNC) derived from independent component analysis. The cohort used included 51 alcohol dependent subjects and 51 control subjects. Despite alcohol, none of the 102 subjects reported use of nicotine, cannabis or any other dependence or habit formation substance. Classification features consisted of whole brain rsFNC estimates undergoing a feature selection process using a random forest approach. Features were then fed to 10 different machine learning classifiers to be evaluated based on their classification performance. A neural network classifier showed the highest performance with an area under the curve (AUC) of 0.79. Other good performers with similar AUC scores were logistic regression, nearest neighbor, and support vector machine classifiers. The worst results were obtained with Gaussian process and quadratic discriminant analysis. The feature selection outcome pointed to functional connections between visual, sensorimotor, executive control, reward, and salience networks as the most relevant for classification. We conclude that AUD can be identified using machine learning classifiers in the absence of nicotine comorbidity.

Prior resting-state functional magnetic resonance imaging (fMRI) analyses have identified patterns of functional connectivity associated with hallucinations in schizophrenia (Sz). In this study, we performed an analysis of the mean amplitude of low-frequency fluctuations (ALFF) to compare resting state spontaneous low-frequency fluctuations in patients with Sz who report experiencing hallucinations impacting different sensory modalities. By exploring dynamics across 2 low-frequency passbands (slow-4 and slow-5), we assessed the impact of hallucination modality and frequency range on spatial ALFF variation. Drawing from a sample of Sz and healthy controls studied as part of the Functional Imaging Biomedical Informatics Research Network (FBIRN), we replicated prior findings showing that patients with Sz have decreased ALFF in the posterior brain in comparison to controls. Remarkably, we found that patients that endorsed visual hallucinations did not show this pattern of reduced ALFF in the back of the brain. These patients also had elevated ALFF in the left hippocampus in comparison to patients that endorsed auditory (but not visual) hallucinations. Moreover, left hippocampal ALFF across all the cases was related to reported hallucination severity in both the auditory and visual domains, and not overall positive symptoms. This supports the hypothesis that dynamic changes in the ALFF in the hippocampus underlie severity of hallucinations that impact different sensory modalities.

ObjectiveWe introduce a method for analyzing dynamically changing functional magnetic resonance imaging (fMRI) network connectivity estimates as they vary within and between broad functional domains. The method captures evidence of intertemporal directionality in cross joint functional-domain influence and extends standard whole-brain dynamic network connectivity approaches into additional functionally meaningful dimensions by evaluating transition probabilities between clustered intradomain and interdomain connectivity patterns.  Results: In applying this method to a large (N = 314) multisite resting-state fMRI dataset balanced between schizophrenia patients and healthy controls, we find evidence of joint functional domains that are global catalyzers, broadly shaping downstream functional relationships throughout the brain. Multiple interesting differences between patients and controls in both time-varying joint functional-domain connectivity patterns and in cross joint functional-domain intertemporal information flow were identified. Conclusion and Significance: Our proposed approach, thus, unifies the concepts of brain connectivity and interdomain connectivity and provides a powerful new way to evaluate functional connectivity data in the context of both the healthy and diseased brain.

Schizophrenia (SZ) is a debilitating mental illness characterized by adolescence or early adulthood onset of psychosis, positive and negative symptoms, as well as cognitive impairments. Despite a plethora of studies leveraging functional connectivity (FC) from functional magnetic resonance imaging (fMRI) to predict symptoms and cognitive impairments of SZ, the findings have exhibited great heterogeneity. We aimed to identify congruous and replicable connectivity patterns capable of predicting positive and negative symptoms as well as cognitive impairments in SZ. Predictable functional connections (FCs) were identified by employing an individualized prediction model, whose replicability was further evaluated across three independent cohorts (BSNIP, SZ = 174; COBRE, SZ = 100; FBIRN, SZ = 161). Across cohorts, we observed that altered FCs in frontal‐temporal‐cingulate‐thalamic network were replicable in prediction of positive symptoms, while sensorimotor network was predictive of negative symptoms. Temporal‐parahippocampal network was consistently identified to be associated with reduced cognitive function. These replicable 23 FCs effectively distinguished SZ from healthy controls (HC) across three cohorts (82.7%, 90.2%, and 86.1%). Furthermore, models built using these replicable FCs showed comparable accuracies to those built using the whole‐brain features in predicting symptoms/cognition of SZ across the three cohorts (r =.17–.33, p <.05). Overall, our findings provide new insights into the neural underpinnings of SZ symptoms/cognition and offer potential targets for further research and possible clinical interventions. [ABSTRACT FROM AUTHOR]

The longitudinal study of typical neurodevelopment is key for understanding deviations due to specific factors, such as psychopathology. However, research utilizing repeated measurements remains scarce. Resting-state functional magnetic resonance imaging (MRI) studies have traditionally examined connectivity as ‘static’ during the measurement period. In contrast, dynamic approaches offer a more comprehensive representation of functional connectivity by allowing for different connectivity configurations (time varying connectivity) throughout the scanning session. Our objective was to characterize the longitudinal developmental changes in dynamic functional connectivity in a population-based pediatric sample. Resting-state MRI data were acquired at the ages of 10 (range 8-to-12, n = 3,327) and 14 (range 13-to-15, n = 2,404) years old using a single, study-dedicated 3 Tesla scanner. A fully-automated spatially constrained group-independent component analysis (ICA) was applied to decompose multi-subject resting-state data into functionally homogeneous regions. Dynamic functional network connectivity (FNC) between all ICA time courses were computed using a tapered sliding window approach. We used a k-means algorithm to cluster the resulting dynamic FNC windows from each scan session into five dynamic states. We examined age and sex associations using linear mixed-effects models. First, independent from the dynamic states, we found a general increase in the temporal variability of the connections between intrinsic connectivity networks with increasing age. Second, when examining the clusters of dynamic FNC windows, we observed that the time spent in less modularized states, with low intra- and inter-network connectivity, decreased with age. Third, the number of transitions between states also decreased with age. Finally, compared to boys, girls showed a more mature pattern of dynamic brain connectivity, indicated by more time spent in a highly modularized state, less time spent in specific states that are frequently observed at a younger age, and a lower number of transitions between states. This longitudinal population-based study demonstrates age-related maturation in dynamic intrinsic neural activity from childhood into adolescence and offers a meaningful baseline for comparison with deviations from typical development. Given that several behavioral and cognitive processes also show marked changes through childhood and adolescence, dynamic functional connectivity should also be explored as a potential neurobiological determinant of such changes.

Clinical research employing functional magnetic resonance imaging (fMRI) is often conducted within the connectionist paradigm, focusing on patterns of connectivity between voxels, regions of interest (ROIs) or spatially distributed functional networks. Connectivity-based analyses are concerned with pairwise correlations of the temporal activation associated with restrictions of the whole-brain hemodynamic signal to locations of a priori interest. There is a more abstract question however that such spatially granular correlation-based approaches do not elucidate: Are the broad spatiotemporal organizing principles of brains in certain populations distinguishable from those of others? Global patterns (in space and time) of hemodynamic activation are rarely scrutinized for features that might characterize complex psychiatric conditions, aging effects or gender-among other variables of potential interest to researchers. We introduce a canonical, transparent technique for characterizing the role in overall brain activation of spatially scaled periodic patterns with given temporal recurrence rates. A core feature of our technique is the spatiotemporal spectral profile (STSP), a readily interpretable 2D reduction of the native four-dimensional brain × time frequency domain that is still "big enough" to capture important group differences in globally patterned brain activation. Its power to distinguish populations of interest is demonstrated on a large balanced multi-site resting fMRI dataset with nearly equal numbers of schizophrenia patients and healthy controls. Our analysis reveals striking differences in the spatiotemporal organization of brain activity that correlate with the presence of diagnosed schizophrenia, as well as with gender and age. To the best of our knowledge, this is the first demonstration that a 4D frequency domain analysis of full volume fMRI data exposes clinically or demographically relevant differences in resting-state brain function.

Network neuroscience, a multidisciplinary field merging insights from neuroscience and network theory, offers a profound understanding of neural network intricacies. However, the impact of varying node sizes on computed graph metrics in neuroimaging data remains underexplored. This study addresses this gap by adopting a data-driven methodology to delineate functional nodes and assess their influence on graph metrics. Using the Neuromark framework, automated independent component analysis is applied to resting state fMRI data, capturing functional network connectivity (FNC) matrices. Global and local graph metrics reveal intricate connectivity patterns, emphasizing the need for nuanced analysis. Notably, node sizes, computed based on voxel counts, contribute to a novel metric termed 'node-metric coupling' (NMC). Correlations between graph metrics and node dimensions are consistently observed. The study extends its analysis to a dataset comprising Alzheimer's disease, mild cognitive impairment, and control subjects, showcasing the potential of NMC as a biomarker for brain disorders. The two key outcomes underscore the interplay between node sizes and resultant graph metrics within a given atlas, shedding light on an often-overlooked source of variability. Additionally, the study highlights the utility of NMC as a valuable biomarker, emphasizing the necessity of accounting for node sizes in future neuroimaging investigations. This work contributes to refining comparative studies employing diverse atlases and advocates for thoughtful consideration of intra-atlas node size in shaping graph metrics, paving the way for more robust neuroimaging research. [ABSTRACT FROM AUTHOR]

Brain networks extracted by independent component analysis (ICA) from magnitude‐only fMRI data are usually denoised using various amplitude‐based thresholds. By contrast, spatial source phase (SSP) or the phase information of ICA brain networks extracted from complex‐valued fMRI data, has provided a simple yet effective way to perform the denoising using a fixed phase change. In this work, we extend the approach to magnitude‐only fMRI data to avoid testing various amplitude thresholds for denoising magnitude maps extracted by ICA, as most studies do not save the complex‐valued data. The main idea is to generate a mathematical SSP map for a magnitude map using a mapping framework, and the mapping framework is built using complex‐valued fMRI data with a known SSP map. Here we leverage the fact that the phase map derived from phase fMRI data has similar phase information to the SSP map. After verifying the use of the magnitude data of complex‐valued fMRI, this framework is generalized to work with magnitude‐only data, allowing use of our approach even without the availability of the corresponding phase fMRI datasets. We test the proposed method using both simulated and experimental fMRI data including complex‐valued data from University of New Mexico and magnitude‐only data from Human Connectome Project. The results provide evidence that the mathematical SSP denoising with a fixed phase change is effective for denoising spatial maps from magnitude‐only fMRI data in terms of retaining more BOLD‐related activity and fewer unwanted voxels, compared with amplitude‐based thresholding. The proposed method provides a unified and efficient SSP approach to denoise ICA brain networks in fMRI data. [ABSTRACT FROM AUTHOR]

It is becoming increasingly common to collect multiple related neuroimaging datasets either from different modalities or from different tasks and conditions. In addition, we have non-imaging data such as cognitive or behavioral variables, and it is through the association of these two sets of data—neuroimaging and non-neuroimaging—that we can understand and explain the evolution of neural and cognitive processes, and predict outcomes for intervention and treatment. Multiple methods for the joint analysis or fusion of multiple neuroimaging datasets or modalities exist; however, methods for the joint analysis of imaging and non-imaging data are still in their infancy. Current approaches for identifying brain networks related to cognitive assessments are still largely based on simple one-to-one correlation analyses and do not use the cross information available across multiple datasets. This work proposes two approaches based on independent vector analysis (IVA) to jointly analyze the imaging datasets and behavioral variables such that multivariate relationships across imaging data and behavioral features can be identified. The simulation results show that our proposed methods provide better accuracy in identifying associations across imaging and behavioral components than current approaches. With functional magnetic resonance imaging (fMRI) task data collected from 138 healthy controls and 109 patients with schizophrenia, results reveal that the central executive network (CEN) estimated in multiple datasets shows a strong correlation with the behavioral variable that measures working memory, a result that is not identified by traditional approaches. Most of the identified fMRI maps also show significant differences in activations across healthy controls and patients potentially providing a useful signature of mental disorders.

This report provides practical recommendations for the design and execution of multicenter functional MRI (MC-fMRI) studies based on the collective experience of the Function Biomedical Informatics Research Network (FBIRN). The study was inspired by many requests from the fMRI community to FBIRN group members for advice on how to conduct MC-fMRI studies. The introduction briefly discusses the advantages and complexities of MC-fMRI studies. Prerequisites for MC-fMRI studies are addressed before delving into the practical aspects of carefully and efficiently setting up a MC-fMRI study. Practical multisite aspects include: (i) establishing and verifying scan parameters including scanner types and magnetic fields, (ii) establishing and monitoring of a scanner quality program, (iii) developing task paradigms and scan session documentation, (iv) establishing clinical and scanner training to ensure consistency over time, (v) developing means for uploading, storing, and monitoring of imaging and other data, (vi) the use of a traveling fMRI expert, and (vii) collectively analyzing imaging data and disseminating results. We conclude that when MC-fMRI studies are organized well with careful attention to unification of hardware, software and procedural aspects, the process can be a highly effective means for accessing a desired participant demographics while accelerating scientific discovery.

This exploratory study set out to investigate dynamic functional connectivity (dFC) in patients with jerky and tremulous functional movement disorders (JT-FMD). The focus in this work is on dynamic brain states, which represent distinct dFC patterns that reoccur in time and across subjects. Resting-state fMRI data were collected from 17 patients with JT-FMD and 17 healthy controls (HC). Symptom severity was measured using the Clinical Global Impression-Severity scale. Depression and anxiety were measured using the Beck Depression Inventory (BDI) and Beck Anxiety Inventory (BAI), respectively. Independent component analysis was used to extract functional brain components. After computing dFC, dynamic brain states were determined for every subject using k-means clustering. Compared to HC, patients with JT-FMD spent more time in a state that was characterized predominantly by increasing medial prefrontal, and decreasing posterior midline connectivity over time. They also tended to visit this state more frequently. In addition, patients with JT-FMD transitioned significantly more often between different states compared to HC, and incorporated a state with decreasing medial prefrontal, and increasing posterior midline connectivity in their attractor, i.e., the cyclic patterns of state transitions. Altogether, this is the first study that demonstrates altered functional brain network dynamics in JT-FMD that may support concepts of increased self-reflective processes and impaired sense of agency as driving factors in FMD.

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.

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.

Fusing complementary information from different modalities can lead to the discovery of more accurate diagnostic biomarkers for psychiatric disorders. However, biomarker discovery through data fusion is challenging since it requires extracting interpretable and reproducible patterns from data sets, consisting of shared/unshared patterns and of different orders. For example, multi-channel electroencephalography (EEG) signals from multiple subjects can be represented as a third-order tensor with modes: subject, time, and channel, while functional magnetic resonance imaging (fMRI) data may be in the form of subject by voxel matrices. Traditional data fusion methods rearrange higher-order tensors, such as EEG, as matrices to use matrix factorization-based approaches. In contrast, fusion methods based on coupled matrix and tensor factorizations (CMTF) exploit the potential multi-way structure of higher-order tensors. The CMTF approach has been shown to capture underlying patterns more accurately without imposing strong constraints on the latent neural patterns, i.e., biomarkers. In this paper, EEG, fMRI, and structural MRI (sMRI) data collected during an auditory oddball task (AOD) from a group of subjects consisting of patients with schizophrenia and healthy controls, are arranged as matrices and higher-order tensors coupled along the subject mode, and jointly analyzed using structure-revealing CMTF methods [also known as advanced CMTF (ACMTF)] focusing on unique identification of underlying patterns in the presence of shared/unshared patterns. We demonstrate that joint analysis of the EEG tensor and fMRI matrix using ACMTF reveals significant and biologically meaningful components in terms of differentiating between patients with schizophrenia and healthy controls while also providing spatial patterns with high resolution and improving the clustering performance compared to the analysis of only the EEG tensor. We also show that these patterns are reproducible, and study reproducibility for different model parameters. In comparison to the joint independent component analysis (jICA) data fusion approach, ACMTF provides easier interpretation of EEG data by revealing a single summary map of the topography for each component. Furthermore, fusion of sMRI data with EEG and fMRI through an ACMTF model provides structural patterns; however, we also show that when fusing data sets from multiple modalities, hence of very different nature, preprocessing plays a crucial role.

Studies of resting state functional MRI (rs-fRMI) are increasingly focused on “dynamics”, or on those properties of brain activation that manifest and vary on timescales shorter than the scan's full duration. This shift in focus has led to a flurry of interest in developing hypothesis testing frameworks and null models applicable to the dynamical setting. Thus far however, these efforts have been weakened by a number of crucial shortcomings that are outlined and discussed in this article. We focus here on aspects of recently proposed null models that, we argue, are poorly formulated relative to the hypotheses they are designed to test, i.e., their potential role in separating functionally relevant BOLD signal dynamics from noise or intermittent background and maintenance type processes is limited by factors that are fundamental rather than merely quantitative or parametric. In this short position paper, we emphasize that (1) serious care must be exercised in building null models for rs-fMRI dynamics from distributionally stationary univariate or multivariate timeseries, i.e., timeseries whose values are each independently drawn from one pre-specified probability distribution; and (2) measures such as kurtosis that quantify over-concentration of observed values in the far tails of some reference distribution may not be particularly suitable for capturing signal features most plausibly contributing to functionally relevant brain dynamics. Other metrics targeted, for example, at capturing the type of epochal signal variation that is often viewed as a signature of brain responsiveness to stimuli or experimental tasks, could play a more scientifically clarifying role. As we learn more about the phenomenon of functionally relevant brain dynamics and its imaging correlates, scientifically meaningful null hypotheses and well-tuned null models will naturally emerge. We also revisit the important concept of distributional stationarity, discuss how it manifests within realizations vs. across multiple realizations, and provide guidance on the benefits and limitations of employing this type of stationarity in modeling the absence of functionally relevant temporal dynamics in resting state fMRI. We hope that the discussions herein are useful, and promote thoughtful consideration of these important issues.

We introduce a novel recurrent neural network (RNN) approach to account for temporal dynamics and dependencies in brain networks observed via functional magnetic resonance imaging (fMRI). Our approach directly parameterizes temporal dynamics through recurrent connections, which can be used to formulate blind source separation with a conditional (rather than marginal) independence assumption, which we call RNN-ICA. This formulation enables us to visualize the temporal dynamics of both first order (activity) and second order (directed connectivity) information in brain networks that are widely studied in a static sense, but not well-characterized dynamically. RNN-ICA predicts dynamics directly from the recurrent states of the RNN in both task and resting state fMRI. Our results show both task-related and group-differentiating directed connectivity.

Brain functional imaging data, especially functional magnetic resonance imaging (fMRI) data, have been employed to reflect functional integration of the brain. Alteration in brain functional connectivity (FC) is expected to provide potential biomarkers for classifying or predicting brain disorders. In this paper, we present a comprehensive review in order to provide guidance about the available brain FC measures and typical classification strategies. We survey the state-of-the-art FC analysis methods including widely used static functional connectivity (SFC) and more recently proposed dynamic functional connectivity (DFC). Temporal correlations among regions of interest (ROIs), data-driven spatial network and functional network connectivity (FNC) are often computed to reflect SFC from different angles. SFC can be extended to DFC using a sliding-window framework, and intrinsic connectivity states along the time-varying connectivity patterns are typically extracted using clustering or decomposition approaches. We also briefly summarize window-less DFC approaches. Subsequently, we highlight various strategies for feature selection including the filter, wrapper and embedded methods. In terms of model building, we include traditional classifiers as well as more recently applied deep learning methods. Moreover, we review representative applications with remarkable classification accuracy for psychosis and mood disorders, neurodevelopmental disorder, and neurological disorders using fMRI data. Schizophrenia, bipolar disorder, autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD), Alzheimer's disease and mild cognitive impairment (MCI) are discussed. Finally, challenges in the field are pointed out with respect to the inaccurate diagnosis labeling, the abundant number of possible features and the difficulty in validation. Some suggestions for future work are also provided.

The study of human brain connectivity, including structural connectivity (SC) and functional connectivity (FC), provides insights into the neurophysiological mechanism of brain function and its relationship to human behavior and cognition. Both types of connectivity measurements provide crucial yet complementary information. However, integrating these two modalities into a single framework remains a challenge, because of the differences in their quantitative interdependencies as well as their anatomical representations due to distinctive imaging mechanisms. In this study, we introduced a new method, joint connectivity matrix independent component analysis (cmICA), which provides a data‐driven parcellation and automated‐linking of SC and FC information simultaneously using a joint analysis of functional magnetic resonance imaging (MRI) and diffusion‐weighted MRI data. We showed that these two connectivity modalities produce common cortical segregation, though with various degrees of (dis)similarity. Moreover, we show conjoint FC networks and structural white matter tracts that directly link these cortical parcellations/sources, within one analysis. Overall, data‐driven joint cmICA provides a new approach for integrating or fusing structural connectivity and FC systematically and conveniently, and provides an effective tool for connectivity‐based multimodal data fusion in brain. [ABSTRACT FROM AUTHOR]

Mental disorders like schizophrenia are currently diagnosed by physicians/psychiatrists through clinical assessment and their evaluation of patient’s self-reported experiences as the illness emerges. There is great interest in identifying biological markers of prognosis at the onset of illness, rather than relying on the evolution of symptoms across time. Functional network connectivity, which indicates a subject's overall level of 'synchronicity' of activity between brain regions, demonstrates promise in providing individual subject predictive power. Many previous studies reported functional connectivity changes during resting-state using only functional magnetic resonance imaging (fMRI). Nevertheless, exclusive reliance on fMRI to generate such networks may limit the inference of the underlying dysfunctional connectivity, which is hypothesized to be a factor in patient symptoms, as fMRI measures connectivity via hemodynamics. Therefore, combination of connectivity assessments using fMRI and magnetoencephalography (MEG), which more directly measures neuronal activity, may provide improved classification of schizophrenia than either modality alone. Moreover, recent evidence indicates that metrics of dynamic connectivity may also be critical for understanding pathology in schizophrenia. In this work, we propose a new framework for extraction of important disease related features and classification of patients with schizophrenia based on using both fMRI and MEG to investigate functional network components in the resting state. Results of this study show that the integration of fMRI and MEG provides important information that captures fundamental characteristics of functional network connectivity in schizophrenia and is helpful for prediction of schizophrenia patient group membership. Combined fMRI/MEG methods, using static functional network connectivity analyses, improved classification accuracy relative to use of fMRI or MEG methods alone (by 15% and 12.45%, respectively), while combined fMRI/MEG methods using dynamic functional network connectivity analyses improved classification up to 5.12% relative to use of fMRI alone and up to 17.21% relative to use of MEG alone.