MRI provides a powerful tool for studying the functional and structural connections in the brain non-invasively. The technique of functional connectivity (FC) exploits the intrinsic temporal correlations of slow spontaneous signal fluctuations to characterise brain functional networks. In addition, diffusion MRI fibre-tracking can be used to study the white matter structural connections. In recent years, there has been considerable interest in combining these two techniques to provide an overall structural-functional description of the brain. In this work we applied the recently proposed super-resolution track-weighted imaging (TWI) methodology to demonstrate how whole-brain fibre-tracking data can be combined with FC data to generate a track-weighted (TW) FC map of FC networks. The method was applied to data from 8 healthy volunteers, and illustrated with (i) FC networks obtained using a seeded connectivity-based analysis (seeding in the precuneus/posterior cingulate cortex, PCC, known to be part of the default mode network), and (ii) with FC networks generated using independent component analysis (in particular, the default mode, attention, visual, and sensory-motor networks). TW-FC maps showed high intensity in white matter structures connecting the nodes of the FC networks. For example, the cingulum bundles show the strongest TW-FC values in the PCC seeded-based analysis, due to their major role in the connection between medial frontal cortex and precuneus/posterior cingulate cortex; similarly the superior longitudinal fasciculus was well represented in the attention network, the optic radiations in the visual network, and the corticospinal tract and corpus callosum in the sensory-motor network. The TW-FC maps highlight the white matter connections associated with a given FC network, and their intensity in a given voxel reflects the functional connectivity of the part of the nodes of the network linked by the structural connections traversing that voxel. They therefore contain a different (and novel) image contrast from that of the images used to generate them. The results shown in this study illustrate the potential of the TW-FC approach for the fusion of structural and functional data into a single quantitative image. This technique could therefore have important applications in neuroscience and neurology, such as for voxel-based comparison studies. <<<

The study of short association fibers is still an incomplete task due to their higher inter-subject variability and the smaller size of this kind of fibers in comparison to known long association bundles. However, their description is essential to understand human brain dysfunction and better characterize the human brain connectome. In this work, we present a multi-subject atlas of short association fibers, which was computed using a superficial white matter bundle identification method based on fiber clustering. To create the atlas, we used probabilistic tractography from one hundred subjects from the HCP database, aligned with non-linear registration. The method starts with an intra-subject clustering of short fibers (30-85 mm). Based on a cortical atlas, the intra-subject cluster centroids from all subjects are segmented to identify the centroids connecting each region of interest (ROI) of the atlas. To reduce computational load, the centroids from each ROI group are randomly separated into ten subgroups. Then, an inter-subject hierarchical clustering is applied to each centroid subgroup, followed by a second level of clustering to select the most-reproducible clusters across subjects for each ROI group. Finally, the clusters are labeled according to the regions that they connect, and clustered to create the final bundle atlas. The resulting atlas is composed of 525 bundles of superficial short association fibers along the whole brain, with 384 bundles connecting pairs of different ROIs and 141 bundles connecting portions of the same ROI. The reproducibility of the bundles was verified using automatic segmentation on three different tractogram databases. Results for deterministic and probabilistic tractography data show high reproducibility, especially for probabilistic tractography in HCP data. In comparison to previous work, our atlas features a higher number of bundles and greater cortical surface coverage. <<<

Alzheimer's disease (AD) is a chronic neurodegenerative disease characterized by progressive dementia, and amnestic mild cognitive impairment (aMCI) has been defined as a transitional stage between normal aging and AD. Accumulating evidence has shown that altered functional connectivity (FC) and structural connectivity (SC) in the default mode network (DMN) is the prominent hallmarks of AD. However, the relationship between the changes in SC and FC of the DMN is not yet clear. In the present study, we derived the FC and SC matrices of the DMN with functional magnetic resonance imaging (fMRI) and diffusion-weighted imaging (DWI) data and further assessed FC and SC abnormalities within a discovery dataset of 120 participants (39 normal controls, 34 patients with aMCI and 47 patients with AD), as well as a replication dataset of 122 participants (43 normal controls, 37 patients with aMCI and 42 patients with AD). Disrupted SC and FC were found among DMN components (e.g., the posterior cingulate cortex (PCC), medial prefrontal cortex (mPFC), and hippocampus) in patients in the aMCI and AD groups in the discovery dataset; most of the disrupted connections were also identified in the replication dataset. More importantly, some SC and FC elements were significantly correlated with the cognitive ability of patients with aMCI and AD. In addition, we found structural-functional decoupling between the PCC and the right hippocampus in patients in the aMCI and AD groups. These findings of the alteration of DMN connectivity in neurodegenerative cohorts deepen our understanding of the pathophysiological mechanisms of AD. <<<

OBJECTIVE: There is substantial variability in reported seizure outcome following pediatric epilepsy surgery, and lack of individualized predictive tools that could evaluate the probability of seizure freedom postsurgery. The aim of this study was to develop and validate a supervised machine learning (ML) model for predicting seizure freedom after pediatric epilepsy surgery. METHODS: This is a multicenter retrospective study of children who underwent epilepsy surgery at five pediatric epilepsy centers in North America. Clinical information, diagnostic investigations, and surgical characteristics were collected, and used as features to predict seizure-free outcome 1 year after surgery. The dataset was split randomly into 80% training and 20% testing data. Thirty-five combinations of five feature sets with seven ML classifiers were assessed on the training cohort using 10-fold cross-validation for model development. The performance of the optimal combination of ML classifier and feature set was evaluated in the testing cohort, and compared with logistic regression, a classical statistical approach. RESULTS: Of the 801 patients included, 61.3% were seizure-free 1 year postsurgery. During model development, the best combination was XGBoost ML algorithm with five features from the univariate feature set, including number of antiseizure medications, magnetic resonance imaging lesion, age at seizure onset, video-electroencephalography concordance, and surgery type, with a mean area under the curve (AUC) of .73 (95% confidence interval [CI] = .69-.77). The combination of XGBoost and univariate feature set was then evaluated on the testing cohort and achieved an AUC of .74 (95% CI = .66-.82; sensitivity = .87, 95% CI = .81-.94; specificity = .58, 95% CI = .47-.71). The XGBoost model outperformed the logistic regression model (AUC = .72, 95% CI = .63-.80; sensitivity = .72, 95% CI = .63-.82; specificity = .66, 95% CI = .53-.77) in the testing cohort (p = .005). SIGNIFICANCE: This study identified important features and validated an ML algorithm, XGBoost, for predicting the probability of seizure freedom after pediatric epilepsy surgery. Improved prognostication of epilepsy surgery is critical for presurgical counseling and will inform treatment decisions. <<<

Integration of the neural processes in the human brain is realized through interconnections that exist between different neural centers. These interconnections take place through white matter pathways. White matter tractography is currently the only available technique for the reconstruction of the anatomical connectivity in the human brain noninvasively and in vivo. The trajectory and terminations of white matter pathways are estimated from local orientations of nerve bundles. These orientations are obtained using measurements of water diffusion in the brain. In this article, the techniques for estimating fiber directions from diffusion measurements in the human brain are reviewed. Methods of white matter tractography are described, together with the current limitations of the technique, including sensitivity to image noise and partial voluming. The applications of white matter tractography to the topographical characterization of the white matter connections and the segmentation of specific white matter pathways, and corresponding functional units of gray matter, are discussed. In this context, the potential impact of white matter tractography in mapping the functional systems and subsystems in the human brain, and their interrelations, is described. Finally, the applications of white matter tractography to the study of brain disorders, including fiber tract localization in brains affected by tumors and the identification of impaired connectivity routes in neurologic and neuropsychiatric diseases, are discussed. <<<

Diffusion weighted (DW) MRI facilitates non-invasive quantification of tissue microstructure and, in combination with appropriate signal processing, three-dimensional estimates of fibrous orientation. In recent years, attention has shifted from the diffusion tensor model, which assumes a unimodal Gaussian diffusion displacement profile to recover fibre orientation (with various well-documented limitations), towards more complex high angular resolution diffusion imaging (HARDI) analysis techniques. Spherical deconvolution (SD) approaches assume that the fibre orientation density function (fODF) within a voxel can be obtained by deconvolving a 'common' single fibre response function from the observed set of DW signals. In practice, this common response function is not known a priori and thus an estimated fibre response must be used. Here the establishment of this single-fibre response function is referred to as 'calibration'. This work examines the vulnerability of two different SD approaches to inappropriate response function calibration: (1) constrained spherical harmonic deconvolution (CSHD)--a technique that exploits spherical harmonic basis sets and (2) damped Richardson-Lucy (dRL) deconvolution--a technique based on the standard Richardson-Lucy deconvolution. Through simulations, the impact of a discrepancy between the calibrated diffusion profiles and the observed ('Target') DW-signals in both single and crossing-fibre configurations was investigated. The results show that CSHD produces spurious fODF peaks (consistent with well known ringing artefacts) as the discrepancy between calibration and target response increases, while dRL demonstrates a lower over-all sensitivity to miscalibration (with a calibration response function for a highly anisotropic fibre being optimal). However, dRL demonstrates a reduced ability to resolve low anisotropy crossing-fibres compared to CSHD. It is concluded that the range and spatial-distribution of expected single-fibre anisotropies within an image must be carefully considered to ensure selection of the appropriate algorithm, parameters and calibration. Failure to choose the calibration response function carefully may severely impact the quality of any resultant tractography. <<<

White matter hyperintensities (WMH) of presumed vascular origin are frequently found in MRIs of healthy older adults. WMH are also associated with aging and cognitive decline. Here, we compared and validated three algorithms for WMH extraction: FreeSurfer (T1w), UBO Detector (T1w + FLAIR), and FSL's Brain Intensity AbNormality Classification Algorithm (BIANCA; T1w + FLAIR) using a longitudinal dataset comprising MRI data of cognitively healthy older adults (baseline N = 231, age range 64-87 years). As reference we manually segmented WMH in T1w, three-dimensional (3D) FLAIR, and two-dimensional (2D) FLAIR images which were used to assess the segmentation accuracy of the different automated algorithms. Further, we assessed the relationships of WMH volumes provided by the algorithms with Fazekas scores and age. FreeSurfer underestimated the WMH volumes and scored worst in Dice Similarity Coefficient (DSC = 0.434) but its WMH volumes strongly correlated with the Fazekas scores (r  = 0.73). BIANCA accomplished the highest DSC (0.602) in 3D FLAIR images. However, the relations with the Fazekas scores were only moderate, especially in the 2D FLAIR images (r  = 0.41), and many outlier WMH volumes were detected when exploring within-person trajectories (2D FLAIR: ~30%). UBO Detector performed similarly to BIANCA in DSC with both modalities and reached the best DSC in 2D FLAIR (0.531) without requiring a tailored training dataset. In addition, it achieved very high associations with the Fazekas scores (2D FLAIR: r  = 0.80). In summary, our results emphasize the importance of carefully contemplating the choice of the WMH segmentation algorithm and MR-modality. <<<

BACKGROUND: Diffusion MRI is the preferred non-invasive in vivo modality for the study of brain white matter connections. Tractography datasets contain 3D streamlines that can be analyzed to study the main brain white matter tracts. Fiber clustering methods have been used to automatically group similar fibers into clusters. However, due to inter-subject variability and artifacts, the resulting clusters are difficult to process for finding common connections across subjects, specially for superficial white matter. METHODS: We present an automatic method for labeling of short association bundles on a group of subjects. The method is based on an intra-subject fiber clustering that generates compact fiber clusters. Posteriorly, the clusters are labeled based on the cortical connectivity of the fibers, taking as reference the Desikan-Killiany atlas, and named according to their relative position along one axis. Finally, two different strategies were applied and compared for the labeling of inter-subject bundles: a matching with the Hungarian algorithm, and a well-known fiber clustering algorithm, called QuickBundles. RESULTS: Individual labeling was executed over four subjects, with an execution time of 3.6 min. An inspection of individual labeling based on a distance measure showed good correspondence among the four tested subjects. Two inter-subject labeling were successfully implemented and applied to 20 subjects and compared using a set of distance thresholds, ranging from a conservative value of 10 mm to a moderate value of 21 mm. Hungarian algorithm led to a high correspondence, but low reproducibility for all the thresholds, with 96 s of execution time. QuickBundles led to better correspondence, reproducibility and short execution time of 9 s. Hence, the whole processing for the inter-subject labeling over 20 subjects takes 1.17 h. CONCLUSION: We implemented a method for the automatic labeling of short bundles in individuals, based on an intra-subject clustering and the connectivity of the clusters with the cortex. The labels provide useful information for the visualization and analysis of individual connections, which is very difficult without any additional information. Furthermore, we provide two fast inter-subject bundle labeling methods. The obtained clusters could be used for performing manual or automatic connectivity analysis in individuals or across subjects. <<<

Phase-synchronization is a manifestation of interaction between neuronal groups measurable from LFP, EEG or MEG signals, however, volume conduction can cause the coherence and the phase locking value to spuriously increase. It has been shown that the imaginary component of the coherency (ImC) cannot be spuriously increased by volume-conduction of independent sources. Recently, it was proposed that the phase lag index (PLI), which estimates to what extent the phase leads and lags between signals from two sensors are nonequiprobable, improves on the ImC. Compared to ImC, PLI has the advantage of being less influenced by phase delays. However, sensitivity to volume-conduction and noise, and capacity to detect changes in phase-synchronization, is hindered by the discontinuity of the PLI, as small perturbations turn phase lags into leads and vice versa. To solve this problem, we introduce a related index, namely the weighted phase lag index (WPLI). Differently from PLI, in WPLI the contribution of the observed phase leads and lags is weighted by the magnitude of the imaginary component of the cross-spectrum. We demonstrate two advantages of the WPLI over the PLI, in terms of reduced sensitivity to additional, uncorrelated noise sources and increased statistical power to detect changes in phase-synchronization. Another factor that can affect phase-synchronization indices is sample-size bias. We show that, when directly estimated, both PLI and the magnitude of the ImC have typically positively biased estimators. To solve this problem, we develop an unbiased estimator of the squared PLI, and a debiased estimator of the squared WPLI. <<<