Cross-scanner reproducibility and harmonization of a diffusion MRI structural brain network: A traveling subject study of multi-b acquisition

Modeling the impact of MRI acquisition bias on structural connectomes: Harmonizing structural connectomes

One way to increase the statistical power and generalizability of neuroimaging studies is to collect data at multiple sites or merge multiple cohorts. However, this usually comes with site-related biases due to the heterogeneity of scanners and acquisition parameters, negatively impacting sensitivity. Brain structural connectomes are not an exception: Being derived from T1-weighted and diffusion-weighted magnetic resonance images, structural connectivity is impacted by differences in imaging protocol. Beyond minimizing acquisition parameter differences, removing bias with postprocessing is essential. In this work we create, from the exhaustive Human Connectome Project Young Adult dataset, a resampled dataset of different b-values and spatial resolutions, modeling a cohort scanned across multiple sites. After demonstrating the statistical impact of acquisition parameters on connectivity, we propose a linear regression with explicit modeling of b-value and spatial resolution, and validate its performance on separate datasets. We show that b-value and spatial resolution affect connectivity in different ways and that acquisition bias can be reduced using a linear regression informed by the acquisition parameters while retaining interindividual differences and hence boosting fingerprinting performance. We also demonstrate the generative potential of our model, and its generalization capability in an independent dataset reflective of typical acquisition practices in clinical settings.

Site effects in multisite fetal brain MRI: morphological insights into early brain development

OBJECTIVE

To evaluate multisite effects on fetal brain MRI. Specifically, to identify crucial acquisition factors affecting fetal brain structural measurements and developmental patterns, while assessing the effectiveness of existing harmonization methods in mitigating site effects.

MATERIALS AND METHODS

Between May 2017 and March 2022, T2-weighted fast spin-echo sequences in-utero MRI were performed on healthy fetuses from retrospectively recruited pregnant volunteers on four different scanners at four sites. A generalized additive model (GAM) was used to quantitatively assess site effects, including field strength (FS), manufacturer (M), in-plane resolution (R), and slice thickness (ST), on subcortical volume and cortical morphological measurements, including cortical thickness, curvature, and sulcal depth. Growth models were selected to elucidate the developmental trajectories of these morphological measurements. Welch's test was performed to evaluate the influence of site effects on developmental trajectories. The comBat-GAM harmonization method was applied to mitigate site-related biases.

RESULTS

The final analytic sample consisted of 340 MRI scans from 218 fetuses (mean GA, 30.1 weeks ± 4.4 [range, 21.7-40 weeks]). GAM results showed that lower FS and lower spatial resolution led to overestimations in selected brain regions of subcortical volumes and cortical morphological measurements. Only the peak cortical thickness in developmental trajectories was significantly influenced by the effects of FS and R. Notably, ComBat-GAM harmonization effectively removed site effects while preserving developmental patterns.

CONCLUSION

Our findings pinpointed the key acquisition factors in in-utero fetal brain MRI and underscored the necessity of data harmonization when pooling multisite data for fetal brain morphology investigations.

KEY POINTS

Question How do specific site MRI acquisition factors affect fetal brain imaging? Finding Lower FS and spatial resolution overestimated subcortical volumes and cortical measurements. Cortical thickness in developmental trajectories was influenced by FS and in-plane resolution. Clinical relevance This study provides important guidelines for the fetal MRI community when scanning fetal brains and underscores the necessity of data harmonization of cross-center fetal studies.

Harmonization of Structural Brain Connectivity through Distribution Matching

The increasing prevalence of multi-site diffusion-weighted magnetic resonance imaging (dMRI) studies potentially offers enhanced statistical power for investigating brain structure. However, these studies face challenges due to variations in scanner hardware and acquisition protocols. While several methods exist for dMRI data harmonization, few specifically address structural brain connectivity. We introduce a new distribution-matching approach to harmonizing structural brain connectivity across different sites and scanners. We evaluate our method using structural brain connectivity data from two distinct datasets of OASIS-3 and ADNI-2, comparing its performance to the widely used ComBat method. Our approach is meant to align the statistical properties of connectivity data from these two datasets. We examine the impact of harmonization on the correlation of brain connectivity with the Mini-Mental State Examination score and age. Our results demonstrate that our distribution-matching technique more effectively harmonizes structural brain connectivity, often producing stronger and more significant correlations compared to ComBat. Qualitative assessments illustrate the desired distributional alignment of ADNI-2 with OASIS-3, while quantitative evaluations confirm robust performance. This work contributes to the growing field of dMRI harmonization, potentially improving the reliability and comparability of structural connectivity studies that combine data from different sources in neuroscientific and clinical research.

Age Trajectories of the Structural Connectome in Child and Adolescent Offspring of Individuals With Bipolar Disorder or Schizophrenia

Background

Offspring of parents with severe mental illness (e.g., bipolar disorder or schizophrenia) are at elevated risk of developing psychiatric illness owing to both genetic predisposition and increased burden of environmental stress. Emerging evidence indicates a disruption of brain network connectivity in young offspring of patients with bipolar disorder and schizophrenia, but the age trajectories of these brain networks in this high-familial-risk population remain to be elucidated.

Methods

A total of 271 T1-weighted and diffusion-weighted scans were obtained from 174 offspring of at least 1 parent diagnosed with bipolar disorder (n = 74) or schizophrenia (n = 51) and offspring of parents without severe mental illness (n = 49). The age range was 8 to 23 years; 97 offspring underwent 2 scans. Anatomical brain networks were reconstructed into structural connectivity matrices. Network analysis was performed to investigate anatomical brain connectivity.

Results

Offspring of parents with schizophrenia had differential trajectories of connectivity strength and clustering compared with offspring of parents with bipolar disorder and parents without severe mental illness, of global efficiency compared with offspring of parents without severe mental illness, and of local connectivity compared with offspring of parents with bipolar disorder.

Conclusions

The findings of this study suggest that familial high risk of schizophrenia is related to deviations in age trajectories of global structural connectome properties and local connectivity strength.

MidRISH: Unbiased harmonization of rotationally invariant harmonics of the diffusion signal

Data harmonization is necessary for removing confounding effects in multi-site diffusion image analysis. One such harmonization method, LinearRISH, scales rotationally invariant spherical harmonic (RISH) features from one site ("target") to the second ("reference") to reduce confounding scanner effects. However, reference and target site designations are not arbitrary and resultant diffusion metrics (fractional anisotropy, mean diffusivity) are biased by this choice. In this work we propose MidRISH: rather than scaling reference RISH features to target RISH features, we project both sites to a mid-space. We validate MidRISH with the following experiments: harmonizing scanner differences from 37 matched patients free of cognitive impairment, and harmonizing acquisition and study differences on 117 matched patients free of cognitive impairment. We find that MidRISH reduces bias of reference selection while preserving harmonization efficacy of LinearRISH. Users should be cautious when performing LinearRISH harmonization. To select a reference site is to choose diffusion metric effect-size. Our proposed method eliminates the bias-inducing site selection step.

Tractometry of Human Visual White Matter Pathways in Health and Disease

Diffusion-weighted MRI (dMRI) provides a unique non-invasive view of human brain tissue properties. The present review article focuses on tractometry analysis methods that use dMRI to assess the properties of brain tissue within the long-range connections comprising brain networks. We focus specifically on the major white matter tracts that convey visual information. These connections are particularly important because vision provides rich information from the environment that supports a large range of daily life activities. Many of the diseases of the visual system are associated with advanced aging, and tractometry of the visual system is particularly important in the modern aging society. We provide an overview of the tractometry analysis pipeline, which includes a primer on dMRI data acquisition, voxelwise model fitting, tractography, recognition of white matter tracts, and calculation of tract tissue property profiles. We then review dMRI-based methods for analyzing visual white matter tracts: the optic nerve, optic tract, optic radiation, forceps major, and vertical occipital fasciculus. For each tract, we review background anatomical knowledge together with recent findings in tractometry studies on these tracts and their properties in relation to visual function and disease. Overall, we find that measurements of the brain's visual white matter are sensitive to a range of disorders and correlate with perceptual abilities. We highlight new and promising analysis methods, as well as some of the current barriers to progress toward integration of these methods into clinical practice. These barriers, such as variability in measurements between protocols and instruments, are targets for future development.

Implementing ABCD studyⓇ MRI sequences for multi-site cohort studies: Practical guide to necessary steps, preprocessing methods, and challenges

Large multi-site studies that combine magnetic resonance imaging (MRI) data across research sites present exceptional opportunities to advance neuroscience research. However, scanner or site variability and non-standardised image acquisition protocols, data processing and analysis pipelines can adversely affect the reliability and repeatability of MRI derived brain measures. We implemented a standardised MRI protocol based on that used in the Adolescent Brain Cognition Development (ABCD)Ⓡ study in two sites, and across four MRI scanners. Twice repeated measurements of a single healthy volunteer were obtained in two sites and in four 3T MRI scanners (vendors: Siemens, Philips, and GE). Imaging data included anatomical scans (T1 weighted, T2 weighted), diffusion weighted imaging (DWI) and resting state functional MRI (rs-fMRI). Standardised containerized pipelines were utilised to pre-process the data and different image quality metrics and test-retest variability of different brain metrics were evaluated. The implementation of the MRI protocols was possible with minor adjustments in acquisition (e.g. repetition time (TR), higher b-values) and exporting (DICOM formats) of images due to different technical performance of the scanners. This study provides practical insights into the implementation of standardised sequences and data processing for multisite studies, showcase the benefits of containerised preprocessing tools, and highlights the need for careful optimisation of multisite image acquisition.

Traveling Subject-Informed Harmonization Increases Reliability of Brain Diffusion Tensor and Neurite Mapping

Diffusion-weighted magnetic resonance imaging (dMRI) of brain has helped elucidate the microstructural changes of psychiatric and neurodegenerative disorders. Inconsistency between MRI models has hampered clinical application of dMRI-based metrics. Using harmonized dMRI data of 300 scans from 69 traveling subjects (TS) scanning the same individuals at multiple conditions with 13 MRI models and 2 protocols, the widely-used metrics such as diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI) were evaluated before and after harmonization with a combined association test (ComBat) or TS-based general linear model (TS-GLM). Results showed that both ComBat and TS-GLM significantly reduced the effects of the MRI site, model, and protocol for diffusion metrics while maintaining the intersubject biological effects. The harmonization power of TS-GLM based on TS data model is more powerful than that of ComBat. In conclusion, our research demonstrated that although ComBat and TS-GLM harmonization approaches were effective at reducing the scanner effects of the site, model, and protocol for DTI and NODDI metrics in WM, they exhibited high retainability of biological effects. Therefore, we suggest that, after harmonizing DTI and NODDI metrics, a multisite study with large cohorts can accurately detect small pathological changes by retaining pathological effects.

Structural Network Topology Reveals Higher Brain Resilience in Individuals with Preclinical Alzheimer's Disease

Introduction: The diagnosis of Alzheimer's disease (AD) requires the presence of amyloid and tau pathology, but it remains unclear how they affect the structural network in the pre-clinical stage. We aimed to assess differences in topological properties in cognitively normal (CN) individuals with varying levels of amyloid and tau pathology, as well as their association with AD pathology burden. Methods: A total of 68 CN individuals were included and stratified by normal/abnormal (-/+) amyloid (A) and tau (T) status based on positron emission tomography results, yielding three groups: A-T- (n = 19), A+T- (n = 28), and A+T+ (n = 21). Topological properties were measured from structural connectivity. Group differences and correlations with A and T were evaluated. Results: Compared with the A-T- group, the A+T+ group exhibited changes in the structural network topology. At the global level, higher assortativity was shown in the A+T+ group and was correlated with greater tau burden (r = 0.29, p = 0.02), while no difference in global efficiency was found across the three groups. At the local level, the A+T+ group showed disrupted topological properties in the left hippocampus compared with the A-T- group, characterized by lower local efficiency (p < 0.01) and a lower clustering coefficient (p = 0.014). Conclusions: The increased linkage in the higher level architecture of the white matter network reflected by assortativity may indicate increased brain resilience in the early pathological state. Our results encourage further investigation of the topological properties of the structural network in pre-clinical AD.

MidRISH: Unbiased harmonization of rotationally invariant harmonics of the diffusion signal

Objective

Data harmonization is necessary for removing confounding effects in multi-site diffusion image analysis. One such harmonization method, LinearRISH, scales rotationally invariant spherical harmonic (RISH) features from one site ("target") to the second ("reference") to reduce confounding scanner effects. However, reference and target site designations are not arbitrary and resultant diffusion metrics (fractional anisotropy, mean diffusivity) are biased by this choice. In this work we propose MidRISH: rather than scaling reference RISH features to target RISH features, we project both sites to a mid-space.

Methods

We validate MidRISH with the following experiments: harmonizing scanner differences from 37 matched patients free of cognitive impairment, and harmonizing acquisition and study differences on 117 matched patients free of cognitive impairment.

Conclusion

MidRISH reduces bias of reference selection while preserving harmonization efficacy of LinearRISH.

Significance

Users should be cautious when performing LinearRISH harmonization. To select a reference site is to choose diffusion metric effect-size. Our proposed method eliminates the bias-inducing site selection step.

Image harmonization: A review of statistical and deep learning methods for removing batch effects and evaluation metrics for effective harmonization

Magnetic resonance imaging and computed tomography from multiple batches (e.g. sites, scanners, datasets, etc.) are increasingly used alongside complex downstream analyses to obtain new insights into the human brain. However, significant confounding due to batch-related technical variation, called batch effects, is present in this data; direct application of downstream analyses to the data may lead to biased results. Image harmonization methods seek to remove these batch effects and enable increased generalizability and reproducibility of downstream results. In this review, we describe and categorize current approaches in statistical and deep learning harmonization methods. We also describe current evaluation metrics used to assess harmonization methods and provide a standardized framework to evaluate newly-proposed methods for effective harmonization and preservation of biological information. Finally, we provide recommendations to end-users to advocate for more effective use of current methods and to methodologists to direct future efforts and accelerate development of the field.

Different MRI structural processing methods do not impact functional connectivity computation

Resting-state functional magnetic resonance imaging (rs-fMRI) has become an increasingly popular technique. This technique can assess several features of brain connectivity, such as inter-regional temporal correlation (functional connectivity), from which graph measures of network organization can be derived. However, these measures are prone to a certain degree of variability depending on the analytical steps during preprocessing. Many studies have investigated the effect of different preprocessing steps on functional connectivity measures; however, no study investigated whether different structural reconstructions lead to different functional connectivity metrics. Here, we evaluated the impact of different structural segmentation strategies on functional connectivity outcomes. To this aim, we compared different metrics computed after two different registration strategies. The first strategy used structural information from the 3D T1-weighted image (unimodal), while the second strategy implemented a multimodal approach, where an additional registration step used the information from the T2-weighted image. The impact of these different approaches was evaluated on a sample of 58 healthy adults. As expected, different approaches led to significant differences in structural measures (i.e., cortical thickness, volume, and gyrification index), with the maximum impact on the insula cortex. However, these differences were only slightly translated to functional metrics. We reported no differences in graph measures and seed-based functional connectivity maps, but slight differences in the insula when we compared the mean functional strength for each parcel. Overall, these results suggested that functional metrics are only slightly different when using a unimodal compared to a multimodal approach, while the structural output can be significantly affected.

A comprehensive evaluation of multicentric reliability of single-subject cortical morphological networks on traveling subjects

Despite the prevalence of research on single-subject cerebral morphological networks in recent years, whether they can offer a reliable way for multicentric studies remains largely unknown. Using two multicentric datasets of traveling subjects, this work systematically examined the inter-site test-retest (TRT) reliabilities of single-subject cerebral morphological networks, and further evaluated the effects of several key factors. We found that most graph-based network measures exhibited fair to excellent reliabilities regardless of different analytical pipelines. Nevertheless, the reliabilities were affected by choices of morphological index (fractal dimension > sulcal depth > gyrification index > cortical thickness), brain parcellation (high-resolution > low-resolution), thresholding method (proportional > absolute), and network type (binarized > weighted). For the factor of similarity measure, its effects depended on the thresholding method used (absolute: Kullback-Leibler divergence > Jensen-Shannon divergence; proportional: Jensen-Shannon divergence > Kullback-Leibler divergence). Furthermore, longer data acquisition intervals and different scanner software versions significantly reduced the reliabilities. Finally, we showed that inter-site reliabilities were significantly lower than intra-site reliabilities for single-subject cerebral morphological networks. Altogether, our findings propose single-subject cerebral morphological networks as a promising approach for multicentric human connectome studies, and offer recommendations on how to determine analytical pipelines and scanning protocols for obtaining reliable results.

Comparing voxel- and feature-wise harmonization of complex graph measures from multiple sites for structural brain network investigation of aging

Complex graph theory measures of brain structural connectomes derived from diffusion weighted images (DWI) provide insight into the network structure of the brain. Further, as the number of available DWI datasets grows, so does the ability to investigate associations in these measures with major biological factors, like age. However, one key hurdle that remains is the presence of scanner effects that can arise from different DWI datasets and confound multisite analyses. Two common approaches to correct these effects are voxel-wise and feature-wise harmonization. However, it is still unclear how to best leverage them for graph-theory analysis of an aging population. Thus, there is a need to better characterize the impact of each harmonization method and their ability to preserve age related features. We investigate this by characterizing four complex graph theory measures (modularity, characteristic path length, global efficiency, and betweenness centrality) in 48 participants aged 55 to 86 from Baltimore Longitudinal Study of Aging (BLSA) and Vanderbilt Memory and Aging Project (VMAP) before and after voxel- and feature-wise harmonization with the Null Space Deep Network (NSDN) and ComBat, respectively. First, we characterize across dataset coefficients of variation (CoV) and find the combination of NSDN and ComBat causes the greatest reduction in CoV followed by ComBat alone then NSDN alone. Second, we reproduce published associations of modularity with age after correcting for other covariates with linear models. We find that harmonization with ComBat or ComBat and NSDN together improves the significance of existing age effects, reduces model residuals, and qualitatively reduces separation between datasets. These results reinforce the efficiency of statistical harmonization on the feature-level with ComBat and suggest that harmonization on the voxel-level is synergistic but may have reduced effect after running through the multiple layers of the connectomics pipeline. Thus, we conclude that feature-wise harmonization improves statistical results, but the addition of biologically informed voxel-based harmonization offers further improvement.

Removing the effects of the site in brain imaging machine-learning - Measurement and extendable benchmark

Multisite machine-learning neuroimaging studies, such as those conducted by the ENIGMA Consortium, need to remove the differences between sites to avoid effects of the site (EoS) that may prevent or fraudulently help the creation of prediction models, leading to impoverished or inflated prediction accuracy. Unfortunately, we have shown earlier that current Methods Aiming to Remove the EoS (MAREoS, e.g., ComBat) cannot remove complex EoS (e.g., including interactions between regions). And complex EoS may bias the accuracy. To overcome this hurdle, groups worldwide are developing novel MAREoS. However, we cannot assess their effectiveness because EoS may either inflate or shrink the accuracy, and MAREoS may both remove the EoS and degrade the data. In this work, we propose a strategy to measure the effectiveness of a MAREoS in removing different types of EoS. FOR MAREOS DEVELOPERS, we provide two multisite MRI datasets with only simple true effects (i.e., detectable by most machine-learning algorithms) and two with only simple EoS (i.e., removable by most MAREoS). First, they should use these datasets to fit machine-learning algorithms after applying the MAREoS. Second, they should use the formulas we provide to calculate the relative accuracy change associated with the MAREoS in each dataset and derive an EoS-removal effectiveness statistic. We also offer similar datasets and formulas for complex true effects and EoS that include first-order interactions. FOR MACHINE-LEARNING RESEARCHERS, we provide an extendable benchmark website to show: a) the types of EoS they should remove for each given machine-learning algorithm and b) the effectiveness of each MAREoS for removing each type of EoS. Relevantly, a MAREoS only able to remove the simple EoS may suffice for simple machine-learning algorithms, whereas more complex algorithms need a MAREoS that can remove more complex EoS. For instance, ComBat removes all simple EoS as needed for predictions based on simple lasso algorithms, but it leaves residual complex EoS that may bias the predictions based on standard support vector machine algorithms.

Exploring personalized structural connectomics for moderate to severe traumatic brain injury

Graph theoretical analysis of the structural connectome has been employed successfully to characterize brain network alterations in patients with traumatic brain injury (TBI). However, heterogeneity in neuropathology is a well-known issue in the TBI population, such that group comparisons of patients against controls are confounded by within-group variability. Recently, novel single-subject profiling approaches have been developed to capture inter-patient heterogeneity. We present a personalized connectomics approach that examines structural brain alterations in five chronic patients with moderate to severe TBI who underwent anatomical and diffusion magnetic resonance imaging. We generated individualized profiles of lesion characteristics and network measures (including personalized graph metric GraphMe plots, and nodal and edge-based brain network alterations) and compared them against healthy reference cases (N = 12) to assess brain damage qualitatively and quantitatively at the individual level. Our findings revealed alterations of brain networks with high variability between patients. With validation and comparison to stratified, normative healthy control comparison cohorts, this approach could be used by clinicians to formulate a neuroscience-guided integrative rehabilitation program for TBI patients, and for designing personalized rehabilitation protocols based on their unique lesion load and connectome.

Improved sensitivity and precision in multicentre diffusion MRI network analysis using thresholding and harmonization

PURPOSE

To investigate if network thresholding and raw data harmonization improve consistency of diffusion MRI (dMRI)-based brain networks while also increasing precision and sensitivity to detect disease effects in multicentre datasets.

METHODS

Brain networks were reconstructed from dMRI of five samples with cerebral small vessel disease (SVD; 629 patients, 166 controls), as a clinically relevant exemplar condition for studies on network integrity. We evaluated consistency of network architecture in age-matched controls, by calculating cross-site differences in connection probability and fractional anisotropy (FA). Subsequently we evaluated precision and sensitivity to disease effects by identifying connections with low FA in sporadic SVD patients relative to controls, using more severely affected patients with a pure form of genetically defined SVD as reference.

RESULTS

In controls, thresholding and harmonization improved consistency of network architecture, minimizing cross-site differences in connection probability and FA. In patients relative to controls, thresholding improved precision to detect disrupted connections by removing false positive connections (precision, before: 0.09-0.19; after: 0.38-0.70). Before harmonization, sensitivity was low within individual sites, with few connections surviving multiple testing correction (k = 0-25 connections). Harmonization and pooling improved sensitivity (k = 38), while also achieving higher precision when combined with thresholding (0.97).

CONCLUSION

We demonstrated that network consistency, precision and sensitivity to detect disease effects in SVD are improved by thresholding and harmonization. We recommend introducing these techniques to leverage large existing multicentre datasets to better understand the impact of disease on brain networks.

Multisite Harmonization of Structural DTI Networks in Children: An A-CAP Study

The analysis of large, multisite neuroimaging datasets provides a promising means for robust characterization of brain networks that can reduce false positives and improve reproducibility. However, the use of different MRI scanners introduces variability to the data. Managing those sources of variability is increasingly important for the generation of accurate group-level inferences. ComBat is one of the most promising tools for multisite (multiscanner) harmonization of structural neuroimaging data, but no study has examined its application to graph theory metrics derived from the structural brain connectome. The present work evaluates the use of ComBat for multisite harmonization in the context of structural network analysis of diffusion-weighted scans from the Advancing Concussion Assessment in Pediatrics (A-CAP) study. Scans were acquired on six different scanners from 484 children aged 8.00-16.99 years [Mean = 12.37 ± 2.34 years; 289 (59.7%) Male] ~10 days following mild traumatic brain injury (n = 313) or orthopedic injury (n = 171). Whole brain deterministic diffusion tensor tractography was conducted and used to construct a 90 x 90 weighted (average fractional anisotropy) adjacency matrix for each scan. ComBat harmonization was applied separately at one of two different stages during data processing, either on the (i) weighted adjacency matrices (matrix harmonization) or (ii) global network metrics derived using unharmonized weighted adjacency matrices (parameter harmonization). Global network metrics based on unharmonized adjacency matrices and each harmonization approach were derived. Robust scanner effects were found for unharmonized metrics. Some scanner effects remained significant for matrix harmonized metrics, but effect sizes were less robust. Parameter harmonized metrics did not differ by scanner. Intraclass correlations (ICC) indicated good to excellent within-scanner consistency between metrics calculated before and after both harmonization approaches. Age correlated with unharmonized network metrics, but was more strongly correlated with network metrics based on both harmonization approaches. Parameter harmonization successfully controlled for scanner variability while preserving network topology and connectivity weights, indicating that harmonization of global network parameters based on unharmonized adjacency matrices may provide optimal results. The current work supports the use of ComBat for removing multiscanner effects on global network topology.

Evaluating reproducibility and subject-specificity of microstructure-informed connectivity

Tractography enables identifying and evaluating the healthy and diseased brain's white matter pathways from diffusion-weighted magnetic resonance imaging data. As previous evaluation studies have reported significant false-positive estimation biases, recent microstructure-informed tractography algorithms have been introduced to improve the trade-off between specificity and sensitivity. However, a major limitation for characterizing the performance of these techniques is the lack of ground truth brain data. In this study, we compared the performance of two relevant microstructure-informed tractography methods, SIFT2 and COMMIT, by assessing the subject specificity and reproducibility of their derived white matter pathways. Specifically, twenty healthy young subjects were scanned at eight different time points at two different sites. Subject specificity and reproducibility were evaluated using the whole-brain connectomes and a subset of 29 white matter bundles. Our results indicate that although the raw tractograms are more vulnerable to the presence of false-positive connections, they are highly reproducible, suggesting that the estimation bias is subject-specific. This high reproducibility was preserved when microstructure-informed tractography algorithms were used to filter the raw tractograms. Moreover, the resulting track-density images depicted a more uniform coverage of streamlines throughout the white matter, suggesting that these techniques could increase the biological meaning of the estimated fascicles. Notably, we observed an increased subject specificity by employing connectivity pre-processing techniques to reduce the underlaying noise and the data dimensionality (using principal component analysis), highlighting the importance of these tools for future studies. Finally, no strong bias from the scanner site or time between measurements was found. The largest intraindividual variance originated from the sole repetition of data measurements (inter-run).