Choosing Wavelet Methods, Filters, and Lengths for Functional Brain Network Construction

Wavelet methods are widely used to decompose fMRI, EEG, or MEG signals into time series representing neurophysiological activity in fixed frequency bands. Using these time series, one can estimate frequency-band specific functional connectivity between sensors or regions of interest, and thereby construct functional brain networks that can be examined from a graph theoretic perspective. Despite their common use, however, practical guidelines for the choice of wavelet method, filter, and length have remained largely undelineated. Here, we explicitly explore the effects of wavelet method (MODWT vs. DWT), wavelet filter (Daubechies Extremal Phase, Daubechies Least Asymmetric, and Coiflet families), and wavelet length (2 to 24) - each essential parameters in wavelet-based methods - on the estimated values of network diagnostics and in their sensitivity to alterations in psychiatric disease. We observe that the MODWT method produces less variable estimates than the DWT method. We also observe that the length of the wavelet filter chosen has a greater impact on the estimated values of network diagnostics than the type of wavelet chosen. Furthermore, wavelet length impacts the sensitivity of the method to detect differences between health and disease and tunes classification accuracy. Collectively, our results suggest that the choice of wavelet method and length significantly alters the reliability and sensitivity of these methods in estimating values of network diagnostics drawn from graph theory. They furthermore demonstrate the importance of reporting the choices utilized in neuroimaging studies and support the utility of exploring wavelet parameters to maximize classification accuracy in the development of biomarkers of psychiatric disease and neurological disorders.Comment: working pape

Disruption of transfer entropy and inter-hemispheric brain functional connectivity in patients with disorder of consciousness

Severe traumatic brain injury can lead to disorders of consciousness (DOC) characterized by deficit in conscious awareness and cognitive impairment including coma, vegetative state, minimally consciousness, and lock-in syndrome. Of crucial importance is to find objective markers that can account for the large-scale disturbances of brain function to help the diagnosis and prognosis of DOC patients and eventually the prediction of the coma outcome. Following recent studies suggesting that the functional organization of brain networks can be altered in comatose patients, this work analyzes brain functional connectivity (FC) networks obtained from resting-state functional magnetic resonance imaging (rs-fMRI). Two approaches are used to estimate the FC: the Partial Correlation (PC) and the Transfer Entropy (TE). Both the PC and the TE show significant statistical differences between the group of patients and control subjects; in brief, the inter-hemispheric PC and the intra-hemispheric TE account for such differences. Overall, these results suggest two possible rs-fMRI markers useful to design new strategies for the management and neuropsychological rehabilitation of DOC patients.Comment: 25 pages; 4 figures; 3 tables; 1 supplementary figure; 4 supplementary tables; accepted for publication in Frontiers in Neuroinformatic

Identifying functional network changing patterns in individuals at clinical high-risk for psychosis and patients with early illness schizophrenia: A group ICA study.

Although individuals at clinical high risk (CHR) for psychosis exhibit a psychosis-risk syndrome involving attenuated forms of the positive symptoms typical of schizophrenia (SZ), it remains unclear whether their resting-state brain intrinsic functional networks (INs) show attenuated or qualitatively distinct patterns of functional dysconnectivity relative to SZ patients. Based on resting-state functional magnetic imaging data from 70 healthy controls (HCs), 53 CHR individuals (among which 41 subjects were antipsychotic medication-naive), and 58 early illness SZ (ESZ) patients (among which 53 patients took antipsychotic medication) within five years of illness onset, we estimated subject-specific INs using a novel group information guided independent component analysis (GIG-ICA) and investigated group differences in INs. We found that when compared to HCs, both CHR and ESZ groups showed significant differences, primarily in default mode, salience, auditory-related, visuospatial, sensory-motor, and parietal INs. Our findings suggest that widespread INs were diversely impacted. More than 25% of voxels in the identified significant discriminative regions (obtained using all 19 possible changing patterns excepting the no-difference pattern) from six of the 15 interrogated INs exhibited monotonically decreasing Z-scores (in INs) from the HC to CHR to ESZ, and the related regions included the left lingual gyrus of two vision-related networks, the right postcentral cortex of the visuospatial network, the left thalamus region of the salience network, the left calcarine region of the fronto-occipital network and fronto-parieto-occipital network. Compared to HCs and CHR individuals, ESZ patients showed both increasing and decreasing connectivity, mainly hypo-connectivity involving 15% of the altered voxels from four INs. The left supplementary motor area from the sensory-motor network and the right inferior occipital gyrus in the vision-related network showed a common abnormality in CHR and ESZ groups. Some brain regions also showed a CHR-unique alteration (primarily the CHR-increasing connectivity). In summary, CHR individuals generally showed intermediate connectivity between HCs and ESZ patients across multiple INs, suggesting that some dysconnectivity patterns evident in ESZ predate psychosis in attenuated form during the psychosis risk stage. Hence, these connectivity measures may serve as possible biomarkers to predict schizophrenia progression

Learning and comparing functional connectomes across subjects

Functional connectomes capture brain interactions via synchronized fluctuations in the functional magnetic resonance imaging signal. If measured during rest, they map the intrinsic functional architecture of the brain. With task-driven experiments they represent integration mechanisms between specialized brain areas. Analyzing their variability across subjects and conditions can reveal markers of brain pathologies and mechanisms underlying cognition. Methods of estimating functional connectomes from the imaging signal have undergone rapid developments and the literature is full of diverse strategies for comparing them. This review aims to clarify links across functional-connectivity methods as well as to expose different steps to perform a group study of functional connectomes

Statistical Network Analysis for Functional MRI: Summary Networks and Group Comparisons

Comparing weighted networks in neuroscience is hard, because the topological properties of a given network are necessarily dependent on the number of edges of that network. This problem arises in the analysis of both weighted and unweighted networks. The term density is often used in this context, in order to refer to the mean edge weight of a weighted network, or to the number of edges in an unweighted one. Comparing families of networks is therefore statistically difficult because differences in topology are necessarily associated with differences in density. In this review paper, we consider this problem from two different perspectives, which include (i) the construction of summary networks, such as how to compute and visualize the mean network from a sample of network-valued data points; and (ii) how to test for topological differences, when two families of networks also exhibit significant differences in density. In the first instance, we show that the issue of summarizing a family of networks can be conducted by adopting a mass-univariate approach, which produces a statistical parametric network (SPN). In the second part of this review, we then highlight the inherent problems associated with the comparison of topological functions of families of networks that differ in density. In particular, we show that a wide range of topological summaries, such as global efficiency and network modularity are highly sensitive to differences in density. Moreover, these problems are not restricted to unweighted metrics, as we demonstrate that the same issues remain present when considering the weighted versions of these metrics. We conclude by encouraging caution, when reporting such statistical comparisons, and by emphasizing the importance of constructing summary networks.Comment: 16 pages, 5 figure