Improvement of diffusion tensor imaging-based tractography by free-water correction in nonedematous gliomas: assessment with brain mapping.

OBJECTIVE: The free-water correction algorithm (Freewater Estimator Using Interpolated Initialization [FERNET]) can be applied to standard diffusion tensor imaging (DTI) tractography to improve visualization of subcortical bundles in the peritumoral area of highly edematous brain tumors. Interest in its use for presurgical planning in purely infiltrative gliomas without peritumoral edema has never been evaluated. Using subcortical maps obtained with direct electrostimulation (DES) in awake surgery as a reference standard, the authors sought to 1) assess the accuracy of preoperative DTI-based tractography with FERNET in a series of nonedematous glioma patients, and 2) determine its potential usefulness in presurgical planning. METHODS: Based on DES-induced functional disturbances and tumor topography, the authors retrospectively reconstructed the putatively stimulated bundles and the peritumoral tracts of interest (various associative and projection pathways) of 12 patients. The tractography data obtained with and without FERNET were compared. RESULTS: The authors identified 21 putative tracts from 24 stimulation sites and reconstituted 49 tracts of interest. The number of streamlines of the putative tracts crossing the DES area was 26.8% higher (96.04 vs 75.75, p = 0.016) and their volume 20.4% higher (13.99 cm3 vs 11.62 cm3, p < 0.0001) with FERNET than with standard DTI. Additionally, the volume of the tracts of interest was 22.1% higher (9.69 cm3 vs 7.93 cm3, p < 0.0001). CONCLUSIONS: Free-water correction significantly increased the anatomical plausibility of the stimulated fascicles and the volume of tracts of interest in the peritumoral area of purely infiltrative nonedematous gliomas. Because of the functional importance of the peritumoral zone, applying FERNET to DTI could have potential implications on surgical planning and the safety of glioma resection.

Free-water correction DTI-based tractography in brain tumor surgery: assessment with functional and electrophysiological mapping of the white matter.

Peritumoral edema prevents fiber tracking from diffusion tensor imaging (DTI). A free-water correction may overcome this drawback, as illustrated in the case of a patient undergoing awake surgery for brain metastasis. The anatomical plausibility and accuracy of tractography with and without free-water correction were assessed with functional mapping and axono-cortical evoked-potentials (ACEPs) as reference methods. The results suggest a potential synergy between corrected DTI-based tractography and ACEPs to reliably identify and preserve white matter tracts during brain tumor surgery.

CrOssing fiber Modeling in the Peritumoral Area using dMRI (COMPARI).

Visualization of fiber tracts around the tumor is critical for neurosurgical planning and preservation of crucial structural connectivity during tumor resection. Biophysical modeling approaches estimate fiber tract orientations from differential water diffusivity information of diffusion MRI. However, the presence of edema and tumor infiltration presents a challenge to visualize crossing fiber tracts in the peritumoral region. Previous approaches proposed free water modeling to compensate for the effect of water diffusivity in edema, but those methods were limited in estimating complex crossing fiber tracts. We propose a new cascaded multi-compartment model to estimate tissue microstructure in the presence of edema and pathological contaminants in the area surrounding brain tumors. In our model (COMPARI), the isotropic components of diffusion signal, including free water and hindered water, were eliminated, and the fiber orientation distribution (FOD) of the remaining signal was estimated. In simulated data, COMPARI accurately recovered fiber orientations in the presence of extracellular water. In a dataset of 23 patients with highly edematous brain tumors, the amplitudes of FOD and anisotropic index distribution within the peritumoral region were higher with COMPARI than with a recently proposed multi-compartment constrained deconvolution model. In a selected patient with metastatic brain tumor, we demonstrated COMPARI's ability to effectively model and eliminate water from the peritumoral region. The white matter bundles reconstructed with our model were qualitatively improved compared to those of other models, and allowed the identification of crossing fibers. In conclusion, the removal of isotropic components as proposed with COMPARI improved the bio-physical modeling of dMRI in edema, thus providing information on crossing fibers, thereby enabling improved tractography in a highly edematous brain tumor. This model may improve surgical planning tools to help achieve maximal safe resection of brain tumors.

Artificial intelligence-based locoregional markers of brain peritumoral microenvironment.

In malignant primary brain tumors, cancer cells infiltrate into the peritumoral brain structures which results in inevitable recurrence. Quantitative assessment of infiltrative heterogeneity in the peritumoral region, the area where biopsy or resection can be hazardous, is important for clinical decision making. Here, we derive a novel set of Artificial intelligence (AI)-based markers capturing the heterogeneity of tumor infiltration, by characterizing free water movement restriction in the peritumoral region using Diffusion Tensor Imaging (DTI)-based free water volume fraction maps. We leverage the differences in the peritumoral region of metastasis and glioblastomas, the former consisting of vasogenic versus the latter containing infiltrative edema, to extract a voxel-wise deep learning-based peritumoral microenvironment index (PMI). Descriptive characteristics of locoregional hubs of uniformly high PMI values are then extracted as AI-based markers to capture distinct aspects of infiltrative heterogeneity. The proposed markers are utilized to stratify patients' survival and IDH1 mutation status on a population of 275 adult-type diffuse gliomas (CNS WHO grade 4). Our results show significant differences in the proposed markers between patients with different overall survival and IDH1 mutation status (t test, Wilcoxon rank sum test, linear regression; p < 0.01). Clustering of patients using the proposed markers reveals distinct survival groups (logrank; p < 10-5, Cox hazard ratio = 1.82; p < 0.005). Our findings provide a panel of markers as surrogates of infiltration that might capture novel insight about underlying biology of peritumoral microstructural heterogeneity, providing potential biomarkers of prognosis pertaining to survival and molecular stratification, with applicability in clinical decision making.

Longitudinal Abnormalities in White Matter Extracellular Free Water Volume Fraction and Neuropsychological Functioning in Patients with Traumatic Brain Injury.

Traumatic brain injury is a global public health problem associated with chronic neurological complications and long-term disability. Biomarkers that map onto the underlying brain pathology driving these complications are urgently needed to identify individuals at risk for poor recovery and to inform design of clinical trials of neuroprotective therapies. Neuroinflammation and neurodegeneration are two endophenotypes potentially associated with increases in brain extracellular water content, but the nature of extracellular free water abnormalities after neurotrauma and its relationship to measures typically thought to reflect traumatic axonal injury are not well characterized. The objective of this study was to describe the relationship between a neuroimaging biomarker of extracellular free water content and the clinical features of a cohort with primarily complicated mild traumatic brain injury. We analyzed a cohort of 59 adult patients requiring hospitalization for non-penetrating traumatic brain injury of all severities as well as 36 healthy controls. Patients underwent brain magnetic resonance imaging (MRI) at 2 weeks (n = 59) and 6 months (n = 29) post-injury, and controls underwent a single MRI. Of the participants with TBI, 50 underwent clinical neuropsychological assessment at 2 weeks and 28 at 6 months. For each subject, we derived a summary score representing deviations in whole brain white matter extracellular free water volume fraction (VF) and free water-corrected fractional anisotropy (fw-FA). The summary specific anomaly score (SAS) for VF was significantly higher in TBI patients at 2 weeks and 6 months post-injury relative to controls. SAS for VF exhibited moderate correlation with neuropsychological functioning, particularly on measures of executive function. These findings indicate abnormalities in whole brain white matter extracellular water fraction in patients with TBI and are an important step toward identifying and validating noninvasive biomarkers that map onto the pathology driving disability after TBI.

3D-QCNet - A pipeline for automated artifact detection in diffusion MRI images.

Artifacts are a common occurrence in Diffusion MRI (dMRI) scans. Identifying and removing them is essential to ensure the accuracy and viability of any post-processing carried out on these scans. This makes quality control (QC) a crucial first step prior to any analysis of dMRI data. Several QC methods for artifact detection exist, however they suffer from problems like requiring manual intervention and the inability to generalize across different artifacts and datasets. In this paper, we propose an automated deep learning (DL) pipeline that utilizes a 3D-Densenet architecture to train a model on diffusion volumes for automatic artifact detection. Our method is validated on 9000 volumes sourced from 7 large clinical datasets spanning different acquisition protocols (with different gradient directions, high and low b-values, single-shell and multi-shell acquisitions) from multiple scanners. Additionally, they represent diverse subject demographics including age, sex and the presence or absence of pathologies. Our QC method is found to accurately generalize across this heterogenous data by correctly detecting 92% artifacts on average across our test set. This consistent performance over diverse datasets underlines the generalizability of our method, which currently is a significant barrier hindering the widespread adoption of automated QC techniques. Thus, 3D-QCNet can be integrated into diffusion pipelines to effectively automate the arduous and time-intensive process of artifact detection.

Laterality and Sex Differences of Human Lateral Habenula Afferent and Efferent Fiber Tracts.

Introduction: The lateral habenula (LHb) is an epithalamic nucleus associated with negative valence and affective disorders. It receives input via the stria medullaris (SM) and sends output via the fasciculus retroflexus (FR). Here, we use tractography to reconstruct and characterize this pathway. Methods: Multi-shell human diffusion magnetic resonance imaging (dMRI) data was obtained from the human connectome project (HCP) (n = 20, 10 males) and from healthy controls (n = 10, 6 males) scanned at our institution. We generated LHb afferents and efferents using probabilistic tractography by selecting the pallidum as the seed region and the ventral tegmental area as the output target. Results: We were able to reconstruct the intended streamlines in all individuals from the HCP dataset and our dataset. Our technique also aided in identification of the LHb. In right-handed individuals, the streamlines were significantly more numerous in the left hemisphere (mean ratio 1.59 ± 0.09, p = 0.04). In left-handed individuals, there was no hemispheric asymmetry on average (mean ratio 1.00 ± 0.09, p = 1.0). Additionally, these streamlines were significantly more numerous in females than in males (619.9 ± 159.7 vs. 225.9 ± 66.03, p = 0.04). Conclusion: We developed a method to reconstruct the SM and FR without manual identification of the LHb. This technique enables targeting of these fiber tracts as well as the LHb. Furthermore, we have demonstrated that there are sex and hemispheric differences in streamline number. These findings may have therapeutic implications and warrant further investigation.

Connectomic assessment of injury burden and longitudinal structural network alterations in moderate-to-severe traumatic brain injury.

Traumatic brain injury (TBI) is a major public health problem. Caused by external mechanical forces, a major characteristic of TBI is the shearing of axons across the white matter, which causes structural connectivity disruptions between brain regions. This diffuse injury leads to cognitive deficits, frequently requiring rehabilitation. Heterogeneity is another characteristic of TBI as severity and cognitive sequelae of the disease have a wide variation across patients, posing a big challenge for treatment. Thus, measures assessing network-wide structural connectivity disruptions in TBI are necessary to quantify injury burden of individuals, which would help in achieving personalized treatment, patient monitoring, and rehabilitation planning. Despite TBI being a disconnectivity syndrome, connectomic assessment of structural disconnectivity has been relatively limited. In this study, we propose a novel connectomic measure that we call network normality score (NNS) to capture the integrity of structural connectivity in TBI patients by leveraging two major characteristics of the disease: diffuseness of axonal injury and heterogeneity of the disease. Over a longitudinal cohort of moderate-to-severe TBI patients, we demonstrate that structural network topology of patients is more heterogeneous and significantly different than that of healthy controls at 3 months postinjury, where dissimilarity further increases up to 12 months. We also show that NNS captures injury burden as quantified by posttraumatic amnesia and that alterations in the structural brain network is not related to cognitive recovery. Finally, we compare NNS to major graph theory measures used in TBI literature and demonstrate the superiority of NNS in characterizing the disease.

Tractography-Based Surgical Targeting for Thalamic Deep Brain Stimulation: A Comparison of Probabilistic vs Deterministic Fiber Tracking of the Dentato-Rubro-Thalamic Tract.

BACKGROUND: The ventral intermediate (VIM) thalamic nucleus is the main target for the surgical treatment of refractory tremor. Initial targeting traditionally relies on atlas-based stereotactic targeting formulas, which only minimally account for individual anatomy. Alternative approaches have been proposed, including direct targeting of the dentato-rubro-thalamic tract (DRTT), which, in clinical settings, is generally reconstructed with deterministic tracking. Whether more advanced probabilistic techniques are feasible on clinical-grade magnetic resonance acquisitions and lead to enhanced reconstructions is poorly understood. OBJECTIVE: To compare DRTT reconstructed with deterministic vs probabilistic tracking. METHODS: This is a retrospective study of 19 patients with essential tremor who underwent deep brain stimulation (DBS) with intraoperative neurophysiology and stimulation testing. We assessed the proximity of the DRTT to the DBS lead and to the active contact chosen based on clinical response. RESULTS: In the commissural plane, the deterministic DRTT was anterior (P < 10-4) and lateral (P < 10-4) to the DBS lead. By contrast, although the probabilistic DRTT was also anterior to the lead (P < 10-4), there was no difference in the mediolateral dimension (P = .5). Moreover, the 3-dimensional Euclidean distance from the active contact to the probabilistic DRTT was smaller vs the distance to the deterministic DRTT (3.32 ± 1.70 mm vs 5.01 ± 2.12 mm; P < 10-4). CONCLUSION: DRTT reconstructed with probabilistic fiber tracking was superior in spatial proximity to the physiology-guided DBS lead and to the empirically chosen active contact. These data inform strategies for surgical targeting of the VIM.

Tractography dissection variability: What happens when 42 groups dissect 14 white matter bundles on the same dataset?

White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human whole-brain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process.

Distinct tumor signatures using deep learning-based characterization of the peritumoral microenvironment in glioblastomas and brain metastases.

Tumor types are classically distinguished based on biopsies of the tumor itself, as well as a radiological interpretation using diverse MRI modalities. In the current study, the overarching goal is to demonstrate that primary (glioblastomas) and secondary (brain metastases) malignancies can be differentiated based on the microstructure of the peritumoral region. This is achieved by exploiting the extracellular water differences between vasogenic edema and infiltrative tissue and training a convolutional neural network (CNN) on the Diffusion Tensor Imaging (DTI)-derived free water volume fraction. We obtained 85% accuracy in discriminating extracellular water differences between local patches in the peritumoral area of 66 glioblastomas and 40 metastatic patients in a cross-validation setting. On an independent test cohort consisting of 20 glioblastomas and 10 metastases, we got 93% accuracy in discriminating metastases from glioblastomas using majority voting on patches. This level of accuracy surpasses CNNs trained on other conventional DTI-based measures such as fractional anisotropy (FA) and mean diffusivity (MD), that have been used in other studies. Additionally, the CNN captures the peritumoral heterogeneity better than conventional texture features, including Gabor and radiomic features. Our results demonstrate that the extracellular water content of the peritumoral tissue, as captured by the free water volume fraction, is best able to characterize the differences between infiltrative and vasogenic peritumoral regions, paving the way for its use in classifying and benchmarking peritumoral tissue with varying degrees of infiltration.

Free Water Volume Fraction: An Imaging Biomarker to Characterize Moderate-to-Severe Traumatic Brain Injury.

Traumatic brain injury (TBI) is a major clinical and public health problem with few therapeutic interventions successfully translated to the clinic. Identifying imaging-based biomarkers characterizing injury severity and predicting long-term functional and cognitive outcomes in TBI patients is crucial for treatment. TBI results in white matter (WM) injuries, which can be detected using diffusion tensor imaging (DTI). Trauma-induced pathologies lead to accumulation of free water (FW) in brain tissue, and standard DTI is susceptible to the confounding effects of FW. In this study, we applied FW DTI to estimate free water volume fraction (FW-VF) in patients with moderate-to-severe TBI and demonstrated its association with injury severity and long-term outcomes. DTI scans and neuropsychological assessments were obtained longitudinally at 3, 6, and 12 months post-injury for 34 patients and once in 35 matched healthy controls. We observed significantly elevated FW-VF in 85 of 90 WM regions in patients compared to healthy controls (p < 0.05). We then presented a patient-specific summary score of WM regions derived using Mahalanobis distance. We observed that MVF at 3 months significantly predicted functional outcome (p = 0.008), executive function (p = 0.005), and processing speed (p = 0.01) measured at 12 months and was significantly correlated with injury severity (p < 0.001). Our findings are an important step toward implementing MVF as a biomarker for personalized therapy and rehabilitation planning for TBI patients.

Enhanced Fiber Tractography Using Edema Correction: Application and Evaluation in High-Grade Gliomas.

BACKGROUND: A limitation of diffusion tensor imaging (DTI)-based tractography is peritumoral edema that confounds traditional diffusion-based magnetic resonance metrics. OBJECTIVE: To augment fiber-tracking through peritumoral regions by performing novel edema correction on clinically feasible DTI acquisitions and assess the accuracy of the fiber-tracks using intraoperative stimulation mapping (ISM), task-based functional magnetic resonance imaging (fMRI) activation maps, and postoperative follow-up as reference standards. METHODS: Edema correction, using our bi-compartment free water modeling algorithm (FERNET), was performed on clinically acquired DTI data from a cohort of 10 patients presenting with suspected high-grade glioma and peritumoral edema in proximity to and/or infiltrating language or motor pathways. Deterministic fiber-tracking was then performed on the corrected and uncorrected DTI to identify tracts pertaining to the eloquent region involved (language or motor). Tracking results were compared visually and quantitatively using mean fiber count, voxel count, and mean fiber length. The tracts through the edematous region were verified based on overlay with the corresponding motor or language task-based fMRI activation maps and intraoperative ISM points, as well as at time points after surgery when peritumoral edema had subsided. RESULTS: Volume and number of fibers increased with application of edema correction; concordantly, mean fractional anisotropy decreased. Overlay with functional activation maps and ISM-verified eloquence of the increased fibers. Comparison with postsurgical follow-up scans with lower edema further confirmed the accuracy of the tracts. CONCLUSION: This method of edema correction can be applied to standard clinical DTI to improve visualization of motor and language tracts in patients with glioma-associated peritumoral edema.

Analysis of Consistency in Structural and Functional Connectivity of Human Brain.

Analysis of structural and functional connectivity of brain has become a fundamental approach in neuroscientific research. Despite several studies reporting consistent similarities as well as differences for structural and resting state (rs) functional connectomes, a comparative investigation of connectomic consistency between the two modalities is still lacking. Nonetheless, connectomic analysis comprising both connectivity types necessitate extra attention as consistency of connectivity differs across modalities, possibly affecting the interpretation of the results. In this study, we present a comprehensive analysis of consistency in structural and rs-functional connectomes obtained from longitudinal diffusion MRI and rs-fMRI data of a single healthy subject. We contrast consistency of deterministic and probabilistic tracking with that of full, positive, and negative functional connectivities across various connectome generation schemes, using correlation as a measure of consistency.

Diffusion MRI tractography filtering techniques change the topology of structural connectomes.

Objective.The use of non-invasive techniques for the estimation of structural brain networks (i.e. connectomes) opened the door to large-scale investigations on the functioning and the architecture of the brain, unveiling the link between neurological disorders and topological changes of the brain network. This study aims at assessing if and how the topology of structural connectomes estimated non-invasively with diffusion MRI is affected by the employment of tractography filtering techniques in structural connectomic pipelines. Additionally, this work investigates the robustness of topological descriptors of filtered connectomes to the common practice of density-based thresholding.Approach.We investigate the changes in global efficiency, characteristic path length, modularity and clustering coefficient on filtered connectomes obtained with the spherical deconvolution informed filtering of tractograms and using the convex optimization modelling for microstructure informed tractography. The analysis is performed on both healthy subjects and patients affected by traumatic brain injury and with an assessment of the robustness of the computed graph-theoretical measures with respect to density-based thresholding of the connectome.Main results.Our results demonstrate that tractography filtering techniques change the topology of brain networks, and thus alter network metrics both in the pathological and the healthy cases. Moreover, the measures are shown to be robust to density-based thresholding.Significance.The present work highlights how the inclusion of tractography filtering techniques in connectomic pipelines requires extra caution as they systematically change the network topology both in healthy subjects and patients affected by traumatic brain injury. Finally, the practice of low-to-moderate density-based thresholding of the connectomes is confirmed to have negligible effects on the topological analysis.

Connectomic consistency: a systematic stability analysis of structural and functional connectivity.

OBJECTIVE: Connectomics, the study of brain connectivity, has become an indispensable tool in neuroscientific research as it provides insights into brain organization. Connectomes are generated using different modalities such as diffusion MRI to capture structural organization of the brain or functional MRI to elaborate brain's functional organization. Understanding links between structural and functional organizations is crucial in explaining how observed behavior emerges from the underlying neurobiological mechanisms. Many studies have investigated how these two organizations relate to each other; however, we still lack a comparative understanding on how much variation should be expected in the two modalities, both between people and within a single person across scans. APPROACH: In this study, we systematically analyzed the consistency of connectomes, that is the similarity between connectomes in terms of individual connections between brain regions and in terms of overall network topology. We present a comprehensive study of consistency in connectomes for a single subject examined longitudinally and across a large cohort of subjects cross-sectionally, in structure and function separately. Within structural connectomes, we compared connectomes generated by different tracking algorithms, parcellations, edge weighting schemes, and edge pruning techniques. In functional connectomes, we compared full, positive, and negative connectivity separately along with thresholding of weak edges. We evaluated consistency using correlation (incorporating information at the level of individual edges) and graph matching accuracy (evaluating connectivity at the level of network topology). We also examined the consistency of connectomes that are generated using different communication schemes. MAIN RESULTS: Our results demonstrate varying degrees of consistency for the two modalities, with structural connectomes showing higher consistency than functional connectomes. Moreover, we observed a wide variation in consistency depending on how connectomes are generated. SIGNIFICANCE: Our study sets a reference point for consistency of connectome types, which is especially important for structure-function coupling studies in evaluating mismatches between modalities.

Freewater estimatoR using iNtErpolated iniTialization (FERNET): Characterizing peritumoral edema using clinically feasible diffusion MRI data.

Characterization of healthy versus pathological tissue in the peritumoral area is confounded by the presence of edema, making free water estimation the key concern in modeling tissue microstructure. Most methods that model tissue microstructure are either based on advanced acquisition schemes not readily available in the clinic or are not designed to address the challenge of edema. This underscores the need for a robust free water elimination (FWE) method that estimates free water in pathological tissue but can be used with clinically prevalent single-shell diffusion tensor imaging data. FWE in single-shell data requires the fitting of a bi-compartment model, which is an ill-posed problem. Its solution requires optimization, which relies on an initialization step. We propose a novel initialization approach for FWE, FERNET, which improves the estimation of free water in edematous and infiltrated peritumoral regions, using single-shell diffusion MRI data. The method has been extensively investigated on simulated data and healthy dataset. Additionally, it has been applied to clinically acquired data from brain tumor patients to characterize the peritumoral region and improve tractography in it.

Graph Matching Based Connectomic Biomarker with Learning for Brain Disorders.

Advances in neuroimaging techniques such as diffusion MRI and functional MRI enabled evaluation of the brain as an information processing network that is called connectome. Connectomic analysis has led to numerous findings on the organization of the brain its pathological changes with diseases, providing imaging-based biomarkers that help in diagnosis and prognosis. A large majority of connectomic biomarkers benefit either from graph-theoretical measures that evaluate brain's network structure, or use standard metrics such as Euclidean distance or Pearson's correlation to show between-connectomes relations. However, such methods are limited in diagnostic evaluation of diseases, because they do not simultaneously measure the difference between individual connectomes, incorporate disease-specific patterns, and utilize network structure information. To address these limitations, we propose a graph matching based method to quantify connectomic similarity, which can be trained for diseases at functional systems level to provide a subject-specific biomarker assessing the disease. We validate our measure on a dataset of patients with traumatic brain injury and demonstrate that our measure achieves better separation between patients and controls compared to commonly used connectomic similarity measures. We further evaluate the vulnerability of the functional systems to the disease by utilizing the parameter tuning aspect of our method. We finally show that our similarity score correlates with clinical scores, highlighting its potential as a subject-specific biomarker for the disease.

Distance-based analysis of variance for brain connectivity.

The field of neuroimaging dedicated to mapping connections in the brain is increasingly being recognized as key for understanding neurodevelopment and pathology. Networks of these connections are quantitatively represented using complex structures, including matrices, functions, and graphs, which require specialized statistical techniques for estimation and inference about developmental and disorder-related changes. Unfortunately, classical statistical testing procedures are not well suited to high-dimensional testing problems. In the context of global or regional tests for differences in neuroimaging data, traditional analysis of variance (ANOVA) is not directly applicable without first summarizing the data into univariate or low-dimensional features, a process that might mask the salient features of high-dimensional distributions. In this work, we consider a general framework for two-sample testing of complex structures by studying generalized within-group and between-group variances based on distances between complex and potentially high-dimensional observations. We derive an asymptotic approximation to the null distribution of the ANOVA test statistic, and conduct simulation studies with scalar and graph outcomes to study finite sample properties of the test. Finally, we apply our test to our motivating study of structural connectivity in autism spectrum disorder.

Deviation from normative brain development is associated with symptom severity in autism spectrum disorder.

Background: Autism spectrum disorder (ASD) is a heterogeneous neurodevelopmental condition. The degree to which the brain development in ASD deviates from typical brain development, and how this deviation relates to observed behavioral outcomes at the individual level are not well-studied. We hypothesize that the degree of deviation from typical brain development of an individual with ASD would relate to observed symptom severity. Methods: The developmental changes in anatomical (cortical thickness, surface area, and volume) and diffusion metrics (fractional anisotropy and apparent diffusion coefficient) were compared between a sample of ASD (n = 247) and typically developing children (TDC) (n = 220) aged 6-25. Machine learning was used to predict age (brain age) from these metrics in the TDC sample, to define a normative model of brain development. This model was then used to compute brain age in the ASD sample. The difference between chronological age and brain age was considered a developmental deviation index (DDI), which was then correlated with ASD symptom severity. Results: Machine learning model trained on all five metrics accurately predicted age in the TDC (r = 0.88) and the ASD (r = 0.85) samples, with dominant contributions to the model from the diffusion metrics. Within the ASD group, the DDI derived from fractional anisotropy was correlated with ASD symptom severity (r = - 0.2), such that individuals with the most advanced brain age showing the lowest severity, and individuals with the most delayed brain age showing the highest severity. Limitations: This work investigated only linear relationships between five specific brain metrics and only one measure of ASD symptom severity in a limited age range. Reported effect sizes are moderate. Further work is needed to investigate developmental differences in other age ranges, other aspects of behavior, other neurobiological measures, and in an independent sample before results can be clinically applicable. Conclusions: Findings demonstrate that the degree of deviation from typical brain development relates to ASD symptom severity, partially accounting for the observed heterogeneity in ASD. Our approach enables characterization of each individual with reference to normative brain development and identification of distinct developmental subtypes, facilitating a better understanding of developmental heterogeneity in ASD.