White matter plasticity following cataract surgery in congenitally blind patients.

The visual system develops abnormally when visual input is absent or degraded during a critical period early in life. Restoration of the visual input later in life is generally thought to have limited benefit because the visual system will lack sufficient plasticity to adapt to and utilize the information from the eyes. Recent evidence, however, shows that congenitally blind adolescents can recover both low-level and higher-level visual function following surgery. In this study, we assessed behavioral performance in both a visual acuity and a face perception task alongside longitudinal structural white matter changes in terms of fractional anisotropy (FA) and mean diffusivity (MD). We studied congenitally blind patients with dense bilateral cataracts, who received cataract surgery at different stages of adolescence. Our goal was to differentiate between age- and surgery-related changes in both behavioral performance and structural measures to identify neural correlates which might contribute to recovery of visual function. We observed surgery-related long-term increases of structural integrity of late-visual pathways connecting the occipital regions with ipsilateral fronto-parieto-temporal regions or homotopic contralateral areas. Comparison to a group of age-matched healthy participants indicated that these improvements went beyond the expected changes in FA and MD based on maturation alone. Finally, we found that the extent of behavioral improvement in face perception was mediated by changes in structural integrity in late visual pathways. Our results suggest that sufficient plasticity remains in adolescence to partially overcome abnormal visual development and help localize the sites of neural change underlying sight recovery.

Early path dominance as a principle for neurodevelopment.

We perform targeted attack, a systematic computational unlinking of the network, to analyze its effects on global communication across the brain network through its giant cluster. Across diffusion magnetic resonance images from individuals in the UK Biobank, Adolescent Brain Cognitive Development Study and Developing Human Connectome Project, we find that targeted attack procedures on increasing white matter tract lengths and densities are remarkably invariant to aging and disease. Time-reversing the attack computation suggests a mechanism for how brains develop, for which we derive an analytical equation using percolation theory. Based on a close match between theory and experiment, our results demonstrate that tracts are limited to emanate from regions already in the giant cluster and tracts that appear earliest in neurodevelopment are those that become the longest and densest.

Development of the Alpha Rhythm Is Linked to Visual White Matter Pathways and Visual Detection Performance.

Alpha is the strongest electrophysiological rhythm in awake humans at rest. Despite its predominance in the EEG signal, large variations can be observed in alpha properties during development, with an increase in alpha frequency over childhood and adulthood. Here, we tested the hypothesis that these changes in alpha rhythm are related to the maturation of visual white matter pathways. We capitalized on a large diffusion MRI (dMRI)-EEG dataset (dMRI n = 2,747, EEG n = 2,561) of children and adolescents of either sex (age range, 5-21 years old) and showed that maturation of the optic radiation specifically accounts for developmental changes of alpha frequency. Behavioral analyses also confirmed that variations of alpha frequency are related to maturational changes in visual perception. The present findings demonstrate the close link between developmental variations in white matter tissue properties, electrophysiological responses, and behavior.

Temporal dynamics of white and gray matter plasticity during motor skill acquisition: a comparative diffusion tensor imaging and multiparametric mapping analysis.

Learning new motor skills relies on neural plasticity within motor and limbic systems. This study uniquely combined diffusion tensor imaging and multiparametric mapping MRI to detail these neuroplasticity processes. We recruited 18 healthy male participants who underwent 960 min of training on a computer-based motion game, while 14 were scanned without training. Diffusion tensor imaging, which quantifies tissue microstructure by measuring the capacity for, and directionality of, water diffusion, revealed mostly linear changes in white matter across the corticospinal-cerebellar-thalamo-hippocampal circuit. These changes related to performance and reflected different responses to upper- and lower-limb training in brain areas with known somatotopic representations. Conversely, quantitative MRI metrics, sensitive to myelination and iron content, demonstrated mostly quadratic changes in gray matter related to performance and reflecting somatotopic representations within the same brain areas. Furthermore, while myelin and iron-sensitive multiparametric mapping MRI was able to describe time lags between different cortical brain systems, diffusion tensor imaging detected time lags within the white matter of the motor systems. These findings suggest that motor skill learning involves distinct phases of white and gray matter plasticity across the sensorimotor network, with the unique combination of diffusion tensor imaging and multiparametric mapping MRI providing complementary insights into the underlying neuroplastic responses.

Clarifying Human White Matter.

Progress in magnetic resonance imaging (MRI) now makes it possible to identify the major white matter tracts in the living human brain. These tracts are important because they carry many of the signals communicated between different brain regions. MRI methods coupled with biophysical modeling can measure the tissue properties and structural features of the tracts that impact our ability to think, feel, and perceive. This review describes the fundamental ideas of the MRI methods used to identify the major white matter tracts in the living human brain.

DCP: A pipeline toolbox for diffusion connectome.

The brain structural network derived from diffusion magnetic resonance imaging (dMRI) reflects the white matter connections between brain regions, which can quantitatively describe the anatomical connection pattern of the entire brain. The development of structural brain connectome leads to the emergence of a large number of dMRI processing packages and network analysis toolboxes. However, the fully automated network analysis based on dMRI data remains challenging. In this study, we developed a cross-platform MATLAB toolbox named "Diffusion Connectome Pipeline" (DCP) for automatically constructing brain structural networks and calculating topological attributes of the networks. The toolbox integrates a few developed packages, including FSL, Diffusion Toolkit, SPM, Camino, MRtrix3, and MRIcron. It can process raw dMRI data collected from any number of participants, and it is also compatible with preprocessed files from public datasets such as HCP and UK Biobank. Moreover, a friendly graphical user interface allows users to configure their processing pipeline without any programming. To prove the capacity and validity of the DCP, two tests were conducted with using DCP. The results showed that DCP can reproduce the findings in our previous studies. However, there are some limitations of DCP, such as relying on MATLAB and being unable to fixel-based metrics weighted network. Despite these limitations, overall, the DCP software provides a standardized, fully automated computational workflow for white matter network construction and analysis, which is beneficial for advancing future human brain connectomics application research.

Time-division multiplexing (TDM) sequence removes bias in T2 estimation and relaxation-diffusion measurements.

PURPOSE: To compare the performance of multi-echo (ME) and time-division multiplexing (TDM) sequences for accelerated relaxation-diffusion MRI (rdMRI) acquisition and to examine their reliability in estimating accurate rdMRI microstructure measures. METHOD: The ME, TDM, and the reference single-echo (SE) sequences with six TEs were implemented using Pulseq with single-band (SB) and multi-band 2 (MB2) acceleration factors. On a diffusion phantom, the image intensities of the three sequences were compared, and the differences were quantified using the normalized RMS error (NRMSE). Shinnar-Le Roux (SLR) pulses were implemented for the SB-ME and SB-SE sequences to investigate the impact of slice profiles on ME sequences. For the in-vivo brain scan, besides the image intensity comparison and T2-estimates, different methods were used to assess sequence-related effects on microstructure estimation, including the relaxation diffusion imaging moment (REDIM) and the maximum-entropy relaxation diffusion distribution (MaxEnt-RDD). RESULTS: TDM performance was similar to the gold standard SE acquisition, whereas ME showed greater biases (3-4× larger NRMSEs for phantom, 2× for in-vivo). T2 values obtained from TDM closely matched SE, whereas ME sequences underestimated the T2 relaxation time. TDM provided similar diffusion and relaxation parameters as SE using REDIM, whereas SB-ME exhibited a 60% larger bias in the map and on average 3.5× larger bias in the covariance between relaxation-diffusion coefficients. CONCLUSION: Our analysis demonstrates that TDM provides a more accurate estimation of relaxation-diffusion measurements while accelerating the acquisitions by a factor of 2 to 3.

Reduced cross-scanner variability using vendor-agnostic sequences for single-shell diffusion MRI.

PURPOSE: To reduce the inter-scanner variability of diffusion MRI (dMRI) measures between scanners from different vendors by developing a vendor-neutral dMRI pulse sequence using the open-source vendor-agnostic Pulseq platform. METHODS: We implemented a standard EPI based dMRI sequence in Pulseq. We tested it on two clinical scanners from different vendors (Siemens Prisma and GE Premier), systematically evaluating and comparing the within- and inter-scanner variability across the vendors, using both the vendor-provided and Pulseq dMRI sequences. Assessments covered both a diffusion phantom and three human subjects, using standard error (SE) and Lin's concordance correlation to measure the repeatability and reproducibility of standard DTI metrics including fractional anisotropy (FA) and mean diffusivity (MD). RESULTS: Identical dMRI sequences were executed on both scanners using Pulseq. On the phantom, the Pulseq sequence showed more than a 2.5× reduction in SE (variability) across Siemens and GE scanners. Furthermore, Pulseq sequences exhibited markedly reduced SE in-vivo, maintaining scan-rescan repeatability while delivering lower variability in FA and MD (more than 50% reduction in cortical/subcortical regions) compared to vendor-provided sequences. CONCLUSION: The Pulseq diffusion sequence reduces the cross-scanner variability for both phantom and in-vivo data, which will benefit multi-center neuroimaging studies and improve the reproducibility of neuroimaging studies.

Reproducibility of Diffusion MRI-Based Tractography in the Fetal Brain.

BACKGROUND: Tractography based on diffusion MRI (dMRI) is a useful tool to study white matter of the developing brain. However, its application in fetal brains is limited due to motion artifacts and low resolution of in utero dMRI, leading to reduced reliability, which was scarcely investigated in previous studies. PURPOSE: To identify reliably traceable fibers in fetal brains and assess whether reproducibility varies with gestational age (GA) and varies between brain regions. STUDY TYPE: Prospective cohort study. SUBJECTS: A total of 44 healthy fetuses with GAs between 25 and 37 (31 ± 6). FIELD STRENGTH/SEQUENCE: 3-T, diffusion-weighted echo-planar imaging sequence (2-5 repeated dMRI scans within the same session per subject). ASSESSMENT: We fitted dMRI with constrained spherical deconvolution model and conducted tractography on eight fibers. We extracted volume, fractional anisotropy, and fiber count for each fiber and assessed the reproducibility of these metrics between repeated scans within each subject. Data were divided into two age-based subgroups (≤30 weeks, N = 28, and >30 weeks, N = 16) for further tests. STATISTICAL TESTS: The reproducibility were compared between fibers by analysis of variance and two-sample t tests. Multiple comparisons were corrected by the false discovery rate (5% was accepted). RESULTS: The reproducibility of the anterior thalamic radiation, inferior longitudinal fasciculus (ILF), genu of the corpus callosum (GCC), and body of the corpus callosum (BCC) significantly decreased with advancing GA (correlation coefficient = 0.525-0.823), as confirmed by group comparisons between fetuses in early GA (≤30 weeks) and late GA (>30 weeks) groups. Corticospinal tract, inferior fronto-occipital fasciculus, and GCC showed high reproducibility for fiber count (weighted dice average = 0.846 vs. 0.814), while BCC and ILF exhibited the lowest reproducibility in both age groups. DATA CONCLUSION: The study indicates that the reliability of fetal brain tractography depends on GA and varies among different fibers. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 2.

Brains of endurance athletes differ in the association areas but not in the primary areas.

Regular participation in sports results in a series of physiological adaptations. However, little is known about the brain adaptations to physical activity. Here we aimed to investigate whether young endurance athletes and non-athletes differ in the gray and white matter of the brain and whether cardiorespiratory fitness (CRF) is associated with these differences. We assessed the CRF, volumes of the gray and white matter of the brain using structural magnetic resonance imaging (sMRI), and brain white matter connections using diffusion magnetic resonance imaging (dMRI) in 20 young male endurance athletes and 21 healthy non-athletes. While total brain volume was similar in both groups, the white matter volume was larger and the gray matter volume was smaller in the athletes compared to non-athletes. The reduction of gray matter was located in the association areas of the brain that are specialized in processing of sensory stimuli. In the microstructure analysis, significant group differences were found only in the association tracts, for example, the inferior occipito-frontal fascicle (IOFF) showing higher fractional anisotropy and lower radial diffusivity, indicating stronger myelination in this tract. Additionally, gray and white matter brain volumes, as well as association tracts correlated with CRF. No changes were observed in other brain areas or tracts. In summary, the brain signature of the endurance athlete is characterized by changes in the integration of sensory and motor information in the association areas.

Translational Connectomics: overview of machine learning in macroscale Connectomics for clinical insights.

Connectomics is a neuroscience paradigm focused on noninvasively mapping highly intricate and organized networks of neurons. The advent of neuroimaging has led to extensive mapping of the brain functional and structural connectome on a macroscale level through modalities such as functional and diffusion MRI. In parallel, the healthcare field has witnessed a surge in the application of machine learning and artificial intelligence for diagnostics, especially in imaging. While reviews covering machine learn ing and macroscale connectomics exist for specific disorders, none provide an overview that captures their evolving role, especially through the lens of clinical application and translation. The applications include understanding disorders, classification, identifying neuroimaging biomarkers, assessing severity, predicting outcomes and intervention response, identifying potential targets for brain stimulation, and evaluating the effects of stimulation intervention on the brain and connectome mapping in patients before neurosurgery. The covered studies span neurodegenerative, neurodevelopmental, neuropsychiatric, and neurological disorders. Along with applications, the review provides a brief of common ML methods to set context. Conjointly, limitations in ML studies within connectomics and strategies to mitigate them have been covered.

Voxel-Wise Fusion of 3T and 7T Diffusion MRI Data to Extract more Accurate Fiber Orientations.

While 7T diffusion magnetic resonance imaging (dMRI) has high spatial resolution, its diffusion imaging quality is usually affected by signal loss due to B1 inhomogeneity, T2 decay, susceptibility, and chemical shift. In contrast, 3T dMRI has relative higher diffusion angular resolution, but lower spatial resolution. Combination of 3T and 7T dMRI, thus, may provide more detailed and accurate information about the voxel-wise fiber orientations to better understand the structural brain connectivity. However, this topic has not yet been thoroughly explored until now. In this study, we explored the feasibility of fusing 3T and 7T dMRI data to extract voxel-wise quantitative parameters at higher spatial resolution. After 3T and 7T dMRI data was preprocessed, respectively, 3T dMRI volumes were coregistered into 7T dMRI space. Then, 7T dMRI data was harmonized to the coregistered 3T dMRI B0 (b = 0) images. Last, harmonized 7T dMRI data was fused with 3T dMRI data according to four fusion rules proposed in this study. We employed high-quality 3T and 7T dMRI datasets (N = 24) from the Human Connectome Project to test our algorithms. The diffusion tensors (DTs) and orientation distribution functions (ODFs) estimated from the 3T-7T fused dMRI volumes were statistically analyzed. More voxels containing multiple fiber populations were found from the fused dMRI data than from 7T dMRI data set. Moreover, extra fiber directions were extracted in temporal brain regions from the fused dMRI data at Otsu's thresholds of quantitative anisotropy, but could not be extracted from 7T dMRI dataset. This study provides novel algorithms to fuse intra-subject 3T and 7T dMRI data for extracting more detailed information of voxel-wise quantitative parameters, and a new perspective to build more accurate structural brain networks.

Genetic risk for schizophrenia is associated with increased proportion of indirect connections in brain networks revealed by a semi-metric analysis: evidence from population sample stratified for polygenic risk.

Research studies based on tractography have revealed a prominent reduction of asymmetry in some key white-matter tracts in schizophrenia (SCZ). However, we know little about the influence of common genetic risk factors for SCZ on the efficiency of routing on structural brain networks (SBNs). Here, we use a novel recall-by-genotype approach, where we sample young adults from a population-based cohort (ALSPAC:N genotyped = 8,365) based on their burden of common SCZ risk alleles as defined by polygenic risk score (PRS). We compared 181 individuals at extremes of low (N = 91) or high (N = 90) SCZ-PRS under a robust diffusion MRI-based graph theoretical SBN framework. We applied a semi-metric analysis revealing higher SMR values for the high SCZ-PRS group compared with the low SCZ-PRS group in the left hemisphere. Furthermore, a hemispheric asymmetry index showed a higher leftward preponderance of indirect connections for the high SCZ-PRS group compared with the low SCZ-PRS group (PFDR < 0.05). These findings might indicate less efficient structural connectivity in the higher genetic risk group. This is the first study in a population-based sample that reveals differences in the efficiency of SBNs associated with common genetic risk variants for SCZ.

Impact of truncating diffusion MRI scans on diffusional kurtosis imaging.

OBJECTIVE: Diffusional kurtosis imaging (DKI) extends diffusion tensor imaging (DTI), characterizing non-Gaussian diffusion effects but requires longer acquisition times. To ensure the robustness of DKI parameters, data acquisition ordering should be optimized allowing for scan interruptions or shortening. Three methodologies were used to examine how reduced diffusion MRI scans impact DKI histogram-metrics: 1) the electrostatic repulsion model (OptEEM); 2) spherical codes (OptSC); 3) random (RandomTRUNC). MATERIALS AND METHODS: Pre-acquired diffusion multi-shell data from 14 female healthy volunteers (29±5 years) were used to generate reordered data. For each strategy, subsets containing different amounts of the full dataset were generated. The subsampling effects were assessed on histogram-based DKI metrics from tract-based spatial statistics (TBSS) skeletonized maps. To evaluate each subsampling method on simulated data at different SNRs and the influence of subsampling on in vivo data, we used a 3-way and 2-way repeated measures ANOVA, respectively. RESULTS: Simulations showed that subsampling had different effects depending on DKI parameter, with fractional anisotropy the most stable (up to 5% error) and radial kurtosis the least stable (up to 26% error). RandomTRUNC performed the worst while the others showed comparable results. Furthermore, the impact of subsampling varied across distinct histogram characteristics, the peak value the least affected (OptEEM: up to 5% error; OptSC: up to 7% error) and peak height (OptEEM: up to 8% error; OptSC: up to 11% error) the most affected. CONCLUSION: The impact of truncation depends on specific histogram-based DKI metrics. The use of a strategy for optimizing the acquisition order is advisable to improve DKI robustness to exam interruptions.

Sorting lung tumor volumes from 4D-MRI data using an automatic tumor-based signal reduces stitching artifacts.

PURPOSE: To investigate whether a novel signal derived from tumor motion allows more precise sorting of 4D-magnetic resonance (4D-MR) image data than do signals based on normal anatomy, reducing levels of stitching artifacts within sorted lung tumor volumes. METHODS: (4D-MRI) scans were collected for 10 lung cancer patients using a 2D T2-weighted single-shot turbo spin echo sequence, obtaining 25 repeat frames per image slice. For each slice, a tumor-motion signal was generated using the first principal component of movement in the tumor neighborhood (TumorPC1). Signals were also generated from displacements of the diaphragm (DIA) and upper and lower chest wall (UCW/LCW) and from slice body area changes (BA). Pearson r coefficients of correlations between observed tumor movement and respiratory signals were determined. TumorPC1, DIA, and UCW signals were used to compile image stacks showing each patient's tumor volume in a respiratory phase. Unsorted image stacks were also built for comparison. For each image stack, the presence of stitching artifacts was assessed by measuring the roughness of the compiled tumor surface according to a roughness metric (Rg). Statistical differences in weighted means of Rg between any two signals were determined using an exact permutation test. RESULTS: The TumorPC1 signal was most strongly correlated with superior-inferior tumor motion, and had significantly higher Pearson r values (median 0.86) than those determined for correlations of UCW, LCW, and BA with superior-inferior tumor motion (p < 0.05). Weighted means of ratios of Rg values in TumorPC1 image stacks to those in unsorted, UCW, and DIA stacks were 0.67, 0.69, and 0.71, all significantly favoring TumorPC1 (p = 0.02-0.05). For other pairs of signals, weighted mean ratios did not differ significantly from one. CONCLUSION: Tumor volumes were smoother in 3D image stacks compiled using the first principal component of tumor motion than in stacks compiled with signals based on normal anatomy.

Probing brain tissue microstructure with MRI: principles, challenges, and the role of multidimensional diffusion-relaxation encoding.

Diffusion MRI uses the random displacement of water molecules to sensitize the signal to brain microstructure and to properties such as the density and shape of cells. Microstructure modeling techniques aim to estimate these properties from acquired data by separating the signal between virtual tissue 'compartments' such as the intra-neurite and the extra-cellular space. A key challenge is that the diffusion MRI signal is relatively featureless compared with the complexity of brain tissue. Another challenge is that the tissue microstructure is wildly different within the gray and white matter of the brain. In this review, we use results from multidimensional diffusion encoding techniques to discuss these challenges and their tentative solutions. Multidimensional encoding increases the information content of the data by varying not only the b-value and the encoding direction but also additional experimental parameters such as the shape of the b-tensor and the echo time. Three main insights have emerged from such encoding. First, multidimensional data contradict common model assumptions on diffusion and T2 relaxation, and illustrates how the use of these assumptions cause erroneous interpretations in both healthy brain and pathology. Second, many model assumptions can be dispensed with if data are acquired with multidimensional encoding. The necessary data can be easily acquired in vivo using protocols optimized to minimize Cramér-Rao lower bounds. Third, microscopic diffusion anisotropy reflects the presence of axons but not dendrites. This insight stands in contrast to current 'neurite models' of brain tissue, which assume that axons in white matter and dendrites in gray matter feature highly similar diffusion. Nevertheless, as an axon-based contrast, microscopic anisotropy can differentiate gray and white matter when myelin alterations confound conventional MRI contrasts.

FIESTA: Autoencoders for accurate fiber segmentation in tractography.

White matter bundle segmentation is a cornerstone of modern tractography to study the brain's structural connectivity in domains such as neurological disorders, neurosurgery, and aging. In this study, we present FIESTA (FIbEr Segmentation in Tractography using Autoencoders), a reliable and robust, fully automated, and easily semi-automatically calibrated pipeline based on deep autoencoders that can dissect and fully populate white matter bundles. This pipeline is built upon previous works that demonstrated how autoencoders can be used successfully for streamline filtering, bundle segmentation, and streamline generation in tractography. Our proposed method improves bundle segmentation coverage by recovering hard-to-track bundles with generative sampling through the latent space seeding of the subject bundle and the atlas bundle. A latent space of streamlines is learned using autoencoder-based modeling combined with contrastive learning. Using an atlas of bundles in standard space (MNI), our proposed method segments new tractograms using the autoencoder latent distance between each tractogram streamline and its closest neighbor bundle in the atlas of bundles. Intra-subject bundle reliability is improved by recovering hard-to-track streamlines, using the autoencoder to generate new streamlines that increase the spatial coverage of each bundle while remaining anatomically correct. Results show that our method is more reliable than state-of-the-art automated virtual dissection methods such as RecoBundles, RecoBundlesX, TractSeg, White Matter Analysis and XTRACT. Our framework allows for the transition from one anatomical bundle definition to another with marginal calibration efforts. Overall, these results show that our framework improves the practicality and usability of current state-of-the-art bundle segmentation framework.

A face-off of MRI research sequences by their need for de-facing.

It is now widely known that research brain MRI, CT, and PET images may potentially be re-identified using face recognition, and this potential can be reduced by applying face-deidentification ("de-facing") software. However, for research MRI sequences beyond T1-weighted (T1-w) and T2-FLAIR structural images, the potential for re-identification and quantitative effects of de-facing are both unknown, and the effects of de-facing T2-FLAIR are also unknown. In this work we examine these questions (where applicable) for T1-w, T2-w, T2*-w, T2-FLAIR, diffusion MRI (dMRI), functional MRI (fMRI), and arterial spin labelling (ASL) sequences. Among current-generation, vendor-product research-grade sequences, we found that 3D T1-w, T2-w, and T2-FLAIR were highly re-identifiable (96-98%). 2D T2-FLAIR and 3D multi-echo GRE (ME-GRE) were also moderately re-identifiable (44-45%), and our derived T2* from ME-GRE (comparable to a typical 2D T2*) matched at only 10%. Finally, diffusion, functional and ASL images were each minimally re-identifiable (0-8%). Applying de-facing with mri_reface version 0.3 reduced successful re-identification to ≤8%, while differential effects on popular quantitative pipelines for cortical volumes and thickness, white matter hyperintensities (WMH), and quantitative susceptibility mapping (QSM) measurements were all either comparable with or smaller than scan-rescan estimates. Consequently, high-quality de-facing software can greatly reduce the risk of re-identification for identifiable MRI sequences with only negligible effects on automated intracranial measurements. The current-generation echo-planar and spiral sequences (dMRI, fMRI, and ASL) each had minimal match rates, suggesting that they have a low risk of re-identification and can be shared without de-facing, but this conclusion should be re-evaluated if they are acquired without fat suppression, with a full-face scan coverage, or if newer developments reduce the current levels of artifacts and distortion around the face.

Joint connectivity matrix independent component analysis: Auto-linking of structural and functional connectivities.

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.

Neural substrates of verbal memory impairment in schizophrenia: A multimodal connectomics study.

While verbal memory is among the most compromised cognitive domains in schizophrenia (SZ), its neural substrates remain elusive. Here, we explored the structural and functional brain network correlates of verbal memory impairment in SZ. We acquired diffusion and resting-state functional MRI data of 49 SZ patients, classified as having preserved (VMP, n = 22) or impaired (VMI, n = 26) verbal memory based on the List Learning task, and 55 healthy controls (HC). Structural and functional connectivity matrices were obtained and analyzed to assess associations with disease status (SZ vs. HC) and verbal memory impairment (VMI vs. VMP) using two complementary data-driven approaches: threshold-free network-based statistics (TFNBS) and hybrid connectivity independent component analysis (connICA). TFNBS showed altered connectivity in SZ patients compared with HC (p < .05, FWER-corrected), with distributed structural changes and functional reorganization centered around sensorimotor areas. Specifically, functional connectivity was reduced within the visual and somatomotor networks and increased between visual areas and associative and subcortical regions. Only a tiny cluster of increased functional connectivity between visual and bilateral parietal attention-related areas correlated with verbal memory dysfunction. Hybrid connICA identified four robust traits, representing fundamental patterns of joint structural-functional connectivity. One of these, mainly capturing the functional connectivity profile of the visual network, was significantly associated with SZ (HC vs. SZ: Cohen's d = .828, p < .0001) and verbal memory impairment (VMP vs. VMI: Cohen's d = -.805, p = .01). We suggest that aberrant connectivity of sensorimotor networks may be a key connectomic signature of SZ and a putative biomarker of SZ-related verbal memory impairment, in consistency with bottom-up models of cognitive disruption.