Early deviation from normal structural connectivity: A novel intrinsic severity score for mild TBI

MVComp toolbox: MultiVariate Comparisons of brain MRI features accounting for common information across metrics

Multivariate approaches have recently gained in popularity to address the physiological unspecificity of neuroimaging metrics and to better characterize the complexity of biological processes underlying behavior. However, commonly used approaches are biased by the intrinsic associations between variables, or they are computationally expensive and may be more complicated to implement than standard univariate approaches. Here, we propose using the Mahalanobis distance (D2), an individual-level measure of deviation relative to a reference distribution that accounts for covariance between metrics. To facilitate its use, we introduce an open-source python-based tool for computing D2 relative to a reference group or within a single individual: the MultiVariate Comparison (MVComp) toolbox. The toolbox allows different levels of analysis (i.e., group- or subject-level), resolutions (e.g., voxel-wise, ROI-wise) and dimensions considered (e.g., combining MRI metrics or WM tracts). Several example cases are presented to showcase the wide range of possible applications of MVComp and to demonstrate the functionality of the toolbox. The D2 framework was applied to the assessment of white matter (WM) microstructure at 1) the group-level, where D2 can be computed between a subject and a reference group to yield an individualized measure of deviation. We observed that clustering applied to D2 in the corpus callosum yields parcellations that highly resemble known topography based on neuroanatomy, suggesting that D2 provides an integrative index that meaningfully reflects the underlying microstructure. 2) At the subject level, D2 was computed between voxels to obtain a measure of (dis)similarity. The loadings of each MRI metric (i.e., its relative contribution to D2) were then extracted in voxels of interest to showcase a useful option of the MVComp toolbox. These relative contributions can provide important insights into the physiological underpinnings of differences observed. Integrative multivariate models are crucial to expand our understanding of the complex brain-behavior relationships and the multiple factors underlying disease development and progression. Our toolbox facilitates the implementation of a useful multivariate method, making it more widely accessible.

A normative modeling approach to quantify white matter changes and predict functional outcomes in stroke patients

Objectives

The diverse nature of stroke necessitates individualized assessment, presenting challenges to case-control neuroimaging studies. The normative model, measuring deviations from a normal distribution, provides a solution. We aim to evaluate stroke-induced white matter microstructural abnormalities at group and individual levels and identify potential prognostic biomarkers.

Methods

Forty-six basal ganglia stroke patients and 46 healthy controls were recruited. Diffusion-weighted imaging and clinical assessment were performed within 7 days after stroke. We used automated fiber quantification to characterize intergroup alterations of segmental diffusion properties along 20 fiber tracts. Then each patient was compared to normative reference (46 healthy participants) by Mahalanobis distance tractometry for 7 significant fiber tracts. Mahalanobis distance-based deviation loads (MaDDLs) and fused MaDDLmulti were extracted to quantify individual deviations. We also conducted correlation and logistic regression analyses to explore relationships between MaDDL metrics and functional outcomes.

Results

Disrupted microstructural integrity was observed across the left corticospinal tract, bilateral inferior fronto-occipital fasciculus, left inferior longitudinal fasciculus, bilateral thalamic radiation, and right uncinate fasciculus. The correlation coefficients between MaDDL metrics and initial functional impairment ranged from 0.364 to 0.618 (p < 0.05), with the highest being MaDDLmulti. Furthermore, MaDDLmulti demonstrated a significant enhancement in predictive efficacy compared to MaDDL (integrated discrimination improvement [IDI] = 9.62%, p = 0.005) and FA (IDI = 34.04%, p < 0.001) of the left corticospinal tract.

Conclusion

MaDDLmulti allows for assessing behavioral disorders and predicting prognosis, offering significant implications for personalized clinical decision-making and stroke recovery. Importantly, our method demonstrates prospects for widespread application in heterogeneous neurological diseases.

MRI-based multivariate gray matter volumetric distance for predicting motor symptom progression in Parkinson's disease

While Parkinson's disease (PD)-related neurodegeneration is associated with structural changes in the brain, conventional magnetic resonance imaging (MRI) has proven less effective for clinical diagnosis due to its inability to reliably identify subtle changes early in the disease course. In this study, we aimed to develop a structural MRI-based biomarker to predict the rate of progression of motor symptoms in the early stages of PD. The study included 88 patients with PD and 120 healthy controls from the Parkinson's Progression Markers Initiative database; MRI at baseline and motor symptom scores assessed using the MDS-UPDRS-III at two time points (baseline and 48 months) were selected. Group-level volumetric analyses revealed that the volumetric reductions in the left striatum were associated with the decline in motor functioning. Then, we developed a patient-specific multivariate gray matter volumetric distance and demonstrated that it could significantly predict changes in motor symptom scores (P < 0.05). Further, we classified patients as relatively slower and faster progressors with 89% accuracy using a support vector machine classifier. Thus, we identified a promising structural MRI-based biomarker for predicting the rate of progression of motor symptoms and classifying patients based on motor symptom severity.

A novel MRI-based volumetric index for monitoring the motor symptoms in Parkinson's disease

BACKGROUND

Conventional MRI scans have limited usefulness in monitoring Parkinson's disease as they typically do not show any disease-specific brain abnormalities. This study aimed to identify an imaging biomarker for tracking motor symptom progression by using a multivariate statistical approach that can combine gray matter volume information from multiple brain regions into a single score specific to each PD patient.

METHODS

A cohort of 150 patients underwent MRI at baseline and had their motor symptoms tracked for up to 10 years using MDS-UPDRS-III, with motor symptoms focused on total and subscores, including rigidity, bradykinesia, postural instability, and gait disturbances, resting tremor, and postural-kinetic tremor. Gray matter volume extracted from MRI data was summarized into a patient-specific summary score using Mahalanobis distance, MGMV. MDS-UPDRS-III's progression and its association with MGMV were modeled via linear mixed-effects models over 5- and 10-year follow-up periods.

RESULTS

Over the 5-year follow-up, there was a significant increase (P < 0.05) in MDS-UPDRS-III total and subscores, except for postural-kinetic tremor. Over the 10-year follow-up, all MDS-UPDRS-III scores increased significantly (P < 0.05). A higher baseline MGMV was associated with a significant increase in MDS-UPDRS-III total, bradykinesia, postural instability and gait disturbances, and resting tremor (P < 0.05) over the 5-year follow-up, but only with total, bradykinesia, and postural instability and gait disturbances during the 10-year follow-up (P < 0.05).

CONCLUSIONS

Higher MGMV scores were linked to faster motor symptom progression, suggesting it could be a valuable marker for clinicians monitoring Parkinson's disease over time.

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.

Structural brain network deviations predict recovery after traumatic brain injury

OBJECTIVE

Traumatic brain injury results in diffuse axonal injury and the ensuing maladaptive alterations in network function are associated with incomplete recovery and persistent disability. Despite the importance of axonal injury as an endophenotype in TBI, there is no biomarker that can measure the aggregate and region-specific burden of axonal injury. Normative modeling is an emerging quantitative case-control technique that can capture region-specific and aggregate deviations in brain networks at the individual patient level. Our objective was to apply normative modeling in TBI to study deviations in brain networks after primarily complicated mild TBI and study its relationship with other validated measures of injury severity, burden of post-TBI symptoms, and functional impairment.

METHOD

We analyzed 70 T1-weighted and diffusion-weighted MRIs longitudinally collected from 35 individuals with primarily complicated mild TBI during the subacute and chronic post-injury periods. Each individual underwent longitudinal blood sampling to characterize blood protein biomarkers of axonal and glial injury and assessment of post-injury recovery in the subacute and chronic periods. By comparing the MRI data of individual TBI participants with 35 uninjured controls, we estimated the longitudinal change in structural brain network deviations. We compared network deviation with independent measures of acute intracranial injury estimated from head CT and blood protein biomarkers. Using elastic net regression models, we identified brain regions in which deviations present in the subacute period predict chronic post-TBI symptoms and functional status.

RESULTS

Post-injury structural network deviation was significantly higher than controls in both subacute and chronic periods, associated with an acute CT lesion and subacute blood levels of glial fibrillary acid protein (r = 0.5, p = 0.008) and neurofilament light (r = 0.41, p = 0.02). Longitudinal change in network deviation associated with change in functional outcome status (r = -0.51, p = 0.003) and post-concussive symptoms (BSI: r = 0.46, p = 0.03; RPQ: r = 0.46, p = 0.02). The brain regions where the node deviation index measured in the subacute period predicted chronic TBI symptoms and functional status corresponded to areas known to be susceptible to neurotrauma.

CONCLUSION

Normative modeling can capture structural network deviations, which may be useful in estimating the aggregate and region-specific burden of network changes induced by TAI. If validated in larger studies, structural network deviation scores could be useful for enrichment of clinical trials of targeted TAI-directed therapies.

Multi-tract multi-symptom relationships in pediatric concussion

Background

The heterogeneity of white matter damage and symptoms in concussion has been identified as a major obstacle to therapeutic innovation. In contrast, most diffusion MRI (dMRI) studies on concussion have traditionally relied on group-comparison approaches that average out heterogeneity. To leverage, rather than average out, concussion heterogeneity, we combined dMRI and multivariate statistics to characterize multi-tract multi-symptom relationships.

Methods

Using cross-sectional data from 306 previously concussed children aged 9-10 from the Adolescent Brain Cognitive Development Study, we built connectomes weighted by classical and emerging diffusion measures. These measures were combined into two informative indices, the first representing microstructural complexity, the second representing axonal density. We deployed pattern-learning algorithms to jointly decompose these connectivity features and 19 symptom measures.

Results

Early multi-tract multi-symptom pairs explained the most covariance and represented broad symptom categories, such as a general problems pair, or a pair representing all cognitive symptoms, and implicated more distributed networks of white matter tracts. Further pairs represented more specific symptom combinations, such as a pair representing attention problems exclusively, and were associated with more localized white matter abnormalities. Symptom representation was not systematically related to tract representation across pairs. Sleep problems were implicated across most pairs, but were related to different connections across these pairs. Expression of multi-tract features was not driven by sociodemographic and injury-related variables, as well as by clinical subgroups defined by the presence of ADHD. Analyses performed on a replication dataset showed consistent results.

Conclusions

Using a double-multivariate approach, we identified clinically-informative, cross-demographic multi-tract multi-symptom relationships. These results suggest that rather than clear one-to-one symptom-connectivity disturbances, concussions may be characterized by subtypes of symptom/connectivity relationships. The symptom/connectivity relationships identified in multi-tract multi-symptom pairs were not apparent in single-tract/single-symptom analyses. Future studies aiming to better understand connectivity/symptom relationships should take into account multi-tract multi-symptom heterogeneity.

Funding

Financial support for this work came from a Vanier Canada Graduate Scholarship from the Canadian Institutes of Health Research (G.I.G.), an Ontario Graduate Scholarship (S.S.), a Restracomp Research Fellowship provided by the Hospital for Sick Children (S.S.), an Institutional Research Chair in Neuroinformatics (M.D.), as well as a Natural Sciences and Engineering Research Council CREATE grant (M.D.).

Demonstration of functional rehabilitation treatment effects in children and young people after severe acquired brain injury

PURPOSE

To examine relationships between functional outcomes after pediatric acquired brain injury (ABI) and measures of rehabilitation dose.

METHODS

An observational study of children receiving residential neurorehabilitation after severe ABI.

RESULTS

Basic total rehabilitation dose shows a paradoxical inverse relationship to global outcome. This is due to confounding by both initial injury severity and length of stay, and variation in treatment content for a given total rehabilitation dose. Content-aware rehabilitation dose measures show robust positive correlations between fractions of rehabilitation treatment received and plausibly related aspects of outcome: specifically, between rates of recovery of gross motor function and the fraction of rehabilitation effort directed to active practice and motor learning. This relationship was robust to adjustment for therapists' expectations of recovery.

CONCLUSION

Content-aware measures of rehabilitation dose are robustly causally related to pertinent aspects of outcome. These findings are step toward a goal of comparative effectiveness research in pediatric neurorehabilitation.

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.

Improving the predictive potential of diffusion MRI in schizophrenia using normative models-Towards subject-level classification

Diffusion MRI studies consistently report group differences in white matter between individuals diagnosed with schizophrenia and healthy controls. Nevertheless, the abnormalities found at the group-level are often not observed at the individual level. Among the different approaches aiming to study white matter abnormalities at the subject level, normative modeling analysis takes a step towards subject-level predictions by identifying affected brain locations in individual subjects based on extreme deviations from a normative range. Here, we leveraged a large harmonized diffusion MRI dataset from 512 healthy controls and 601 individuals diagnosed with schizophrenia, to study whether normative modeling can improve subject-level predictions from a binary classifier. To this aim, individual deviations from a normative model of standard (fractional anisotropy) and advanced (free-water) dMRI measures, were calculated by means of age and sex-adjusted z-scores relative to control data, in 18 white matter regions. Even though larger effect sizes are found when testing for group differences in z-scores than are found with raw values (p < .001), predictions based on summary z-score measures achieved low predictive power (AUC < 0.63). Instead, we find that combining information from the different white matter tracts, while using multiple imaging measures simultaneously, improves prediction performance (the best predictor achieved AUC = 0.726). Our findings suggest that extreme deviations from a normative model are not optimal features for prediction. However, including the complete distribution of deviations across multiple imaging measures improves prediction, and could aid in subject-level classification.

Detecting microstructural deviations in individuals with deep diffusion MRI tractometry

Most diffusion magnetic resonance imaging studies of disease rely on statistical comparisons between large groups of patients and healthy participants to infer altered tissue states in the brain; however, clinical heterogeneity can greatly challenge their discriminative power. There is currently an unmet need to move away from the current approach of group-wise comparisons to methods with the sensitivity to detect altered tissue states at the individual level. This would ultimately enable the early detection and interpretation of microstructural abnormalities in individual patients, an important step towards personalized medicine in translational imaging. To this end, Detect was developed to advance diffusion magnetic resonance imaging tractometry towards single-patient analysis. By operating on the manifold of white-matter pathways and learning normative microstructural features, our framework captures idiosyncrasies in patterns along white-matter pathways. Our approach paves the way from traditional group-based comparisons to true personalized radiology, taking microstructural imaging from the bench to the bedside.

Quantitative assessment of structural and functional changes in temporal lobe epilepsy with hippocampal sclerosis

Background

Magnetic resonance imaging (MRI) changes in hippocampal sclerosis (HS) could be subtle in a significant proportion of mesial temporal lobe epilepsy (mTLE) patients. In this study, we aimed to document the structural and functional changes in the hippocampus and amygdala seen in HS patients.

Methods

Quantitative features of the hippocampus and amygdala were extracted from structural MRI data in 66 mTLE patients and 28 controls. Structural covariance analysis was undertaken using volumetric data from the amygdala and hippocampus. Functional connectivity (FC) measured using resting intracranial electroencephalography (EEG) was analyzed in 22 HS patients and 16 non-HS disease controls.

Results

Hippocampal atrophy was present in both MRI-positive and MRI-negative HS groups (Mann-Whitney U: 7.61, P<0.01; Mann-Whitney U: 6.51, P<0.01). Amygdala volumes were decreased in the patient group (Mann-Whitney U: 2.92, P<0.05), especially in MRI-negative HS patients (Mann-Whitney U: 2.75, P<0.05). The structural covariance analysis showed the normalized volumes of the amygdala and hippocampus were tightly coupled in both controls and HS patients (ρ_Spearman =0.72, P<0.01). FC analysis indicated that HS patients had significantly increased connectivity (Student's t: 2.58, P=0.03) within the hippocampus but decreased connectivity between the hippocampus and amygdala (Student's t: 3.33, P=0.01), particularly for MRI-negative HS patients.

Conclusions

Quantitative structural changes, including hippocampal atrophy and temporal pole blurring, are present in both MRI-positive and MRI-negative HS patients, suggesting the potential usefulness of incorporating quantitative analyses into clinical practice. HS is characterized by increased intra-hippocampal EEG synchronization and decreased coupling between the hippocampus and amygdala.

MRI network progression in mesial temporal lobe epilepsy related to healthy brain architecture

We measured MRI network progression in mesial temporal lobe epilepsy (mTLE) patients as a function of healthy brain architecture. Resting-state functional MRI and diffusion-weighted MRI were acquired in 40 unilateral mTLE patients and 70 healthy controls. Data were used to construct region-to-region functional connectivity, structural connectivity, and streamline length connectomes per subject. Three models of distance from the presumed seizure focus in the anterior hippocampus in the healthy brain were computed using the average connectome across controls. A fourth model was defined using regions of transmodal (higher cognitive function) to unimodal (perceptual) networks across a published functional gradient in the healthy brain. These models were used to test whether network progression in patients increased when distance from the anterior hippocampus or along a functional gradient in the healthy brain decreases. Results showed that alterations of structural and functional networks in mTLE occur in greater magnitude in regions of the brain closer to the seizure focus based on healthy brain topology, and decrease as distance from the focus increases over duration of disease. Overall, this work provides evidence that changes across the brain in focal epilepsy occur along healthy brain architecture.

Multivariate white matter alterations are associated with epilepsy duration

Previous studies investigating associations between white matter alterations and duration of temporal lobe epilepsy (TLE) have shown differing results, and were typically limited to univariate analyses of tracts in isolation. In this study, we apply a multivariate measure (the Mahalanobis distance), which captures the distinct ways white matter may differ in individual patients, and relate this to epilepsy duration. Diffusion MRI, from a cohort of 94 subjects (28 healthy controls, 33 left-TLE and 33 right-TLE), was used to assess the association between tract fractional anisotropy (FA) and epilepsy duration. Using ten white matter tracts, we analysed associations using the traditional univariate analysis (z-scores) and a complementary multivariate approach (Mahalanobis distance), incorporating multiple white matter tracts into a single unified analysis. For patients with right-TLE, FA was not significantly associated with epilepsy duration for any tract studied in isolation. For patients with left-TLE, the FA of two limbic tracts (ipsilateral fornix, contralateral cingulum gyrus) were significantly negatively associated with epilepsy duration (Bonferonni corrected p < .05). Using a multivariate approach we found significant ipsilateral positive associations with duration in both left, and right-TLE cohorts (left-TLE: Spearman's ρ = 0.487, right-TLE: Spearman's ρ = 0.422). Extrapolating our multivariate results to duration equals zero (i.e., at onset) we found no significant difference between patients and controls. Associations using the multivariate approach were more robust than univariate methods. The multivariate Mahalanobis distance measure provides non-overlapping and more robust results than traditional univariate analyses. Future studies should consider adopting both frameworks into their analysis in order to ascertain a more complete understanding of epilepsy progression, regardless of laterality.

Independent components of human brain morphology

Quantification of brain morphology has become an important cornerstone in understanding brain structure. Measures of cortical morphology such as thickness and surface area are frequently used to compare groups of subjects or characterise longitudinal changes. However, such measures are often treated as independent from each other. A recently described scaling law, derived from a statistical physics model of cortical folding, demonstrates that there is a tight covariance between three commonly used cortical morphology measures: cortical thickness, total surface area, and exposed surface area. We show that assuming the independence of cortical morphology measures can hide features and potentially lead to misinterpretations. Using the scaling law, we account for the covariance between cortical morphology measures and derive novel independent measures of cortical morphology. By applying these new measures, we show that new information can be gained; in our example we show that distinct morphological alterations underlie healthy ageing compared to temporal lobe epilepsy, even on the coarse level of a whole hemisphere. We thus provide a conceptual framework for characterising cortical morphology in a statistically valid and interpretable manner, based on theoretical reasoning about the shape of the cortex.

Investigating Brain Network Changes and Their Association With Cognitive Recovery After Traumatic Brain Injury: A Longitudinal Analysis

Traumatic brain injury (TBI) can result in acute cognitive deficits and diffuse axonal injury reflected in white matter brain network alterations, which may, or may not, later recover. Our objective is to first characterize the ways in which brain networks change after TBI and, second, investigate if those changes are associated with recovery of cognitive deficits. We aim to make initial progress in discerning the relationships between brain network changes, and their (dys)functional correlates. We analyze longitudinally acquired MRI from 23 TBI patients (two time points: 6 days, 12 months post-injury) and cross-sectional data from 28 controls to construct white matter brain networks. Cognitive assessment was also performed. Graph theory and regression analysis were applied to identify changed brain network metrics after injury that are associated with subsequent improvements in cognitive function. Sixteen brain network metrics were found to be discriminative of different post-injury phases. Eleven of those explain 90% (adjusted R 2) of the variability observed in cognitive recovery following TBI. Brain network metrics that had a high contribution to the explained variance were found in frontal and temporal cortex, additional to the anterior cingulate cortex. Our preliminary study suggests that network reorganization may be related to recovery of impaired cognitive function in the first year after a TBI.