Characterizing the Contribution of Dependent Features in XAI Methods.

Explainable Artificial Intelligence (XAI) provides tools to help understanding how AI models work and reach a particular decision or outcome. It helps to increase the interpretability of models and makes them more trustworthy and transparent. In this context, many XAI methods have been proposed to make black-box and complex models more digestible from a human perspective. However, one of the main issues that XAI methods have to face especially when dealing with a high number of features is the presence of multicollinearity, which casts shadows on the robustness of the XAI outcomes, such as the ranking of informative features. Most of the current XAI methods either do not consider the collinearity or assume the features are independent which, in general, is not necessarily true. Here, we propose a simple, yet useful, proxy that modifies the outcome of any XAI feature ranking method allowing to account for the dependency among the features, and to reveal their impact on the outcome. The proposed method was applied to SHAP, as an example of XAI method which assume that the features are independent. For this purpose, several models were exploited for a well-known classification task (males versus females) using nine cardiac phenotypes extracted from cardiac magnetic resonance imaging as features. Principal component analysis and biological plausibility were employed to validate the proposed method. Our results showed that the proposed proxy could lead to a more robust list of informative features compared to the original SHAP in presence of collinearity.

Leukocyte Telomere Length and Cardiac Structure and Function: A Mendelian Randomization Study.

BACKGROUND: Existing research demonstrates the association of shorter leukocyte telomere length with increased risk of age-related health outcomes including cardiovascular diseases. However, the direct causality of these relationships has not been definitively established. Cardiovascular aging at an organ level may be captured using image-derived phenotypes of cardiac anatomy and function. METHODS AND RESULTS: In the current study, we use 2-sample Mendelian randomization to assess the causal link between leukocyte telomere length and 54 cardiac magnetic resonance imaging measures representing structure and function across the 4 cardiac chambers. Genetically predicted shorter leukocyte telomere length was causally linked to smaller ventricular cavity sizes including left ventricular end-systolic volume, left ventricular end-diastolic volume, lower left ventricular mass, and pulmonary artery. The association with left ventricular mass (ß =0.217, Pfalse discovery rate=0.016) remained significant after multiple testing adjustment, whereas other associations were attenuated. CONCLUSIONS: Our findings support a causal role for shorter leukocyte telomere length and faster cardiac aging, with the most prominent relationship with left ventricular mass.

Ischemic Heart Disease and Vascular Risk Factors Are Associated With Accelerated Brain Aging.

BACKGROUND: Ischemic heart disease (IHD) has been linked with poor brain outcomes. The brain magnetic resonance imaging-derived difference between predicted brain age and actual chronological age (brain-age delta in years, positive for accelerated brain aging) may serve as an effective means of communicating brain health to patients to promote healthier lifestyles. OBJECTIVES: The authors investigated the impact of prevalent IHD on brain aging, potential underlying mechanisms, and its relationship with dementia risk, vascular risk factors, cardiovascular structure, and function. METHODS: Brain age was estimated in subjects with prevalent IHD (n = 1,341) using a Bayesian ridge regression model with 25 structural (volumetric) brain magnetic resonance imaging features and built using UK Biobank participants with no prevalent IHD (n = 35,237). RESULTS: Prevalent IHD was linked to significantly accelerated brain aging (P < 0.001) that was not fully mediated by microvascular injury. Brain aging (positive brain-age delta) was associated with increased risk of dementia (OR: 1.13 [95% CI: 1.04-1.22]; P = 0.002), vascular risk factors (such as diabetes), and high adiposity. In the absence of IHD, brain aging was also associated with cardiovascular structural and functional changes typically observed in aging hearts. However, such alterations were not linked with risk of dementia. CONCLUSIONS: Prevalent IHD and coexisting vascular risk factors are associated with accelerated brain aging and risk of dementia. Positive brain-age delta representing accelerated brain aging may serve as an effective communication tool to show the impact of modifiable risk factors and disease supporting preventative strategies.

Telomere length is causally connected to brain MRI image derived phenotypes: A mendelian randomization study.

Recent evidence suggests that shorter telomere length (TL) is associated with neuro degenerative diseases and aging related outcomes. The causal association between TL and brain characteristics represented by image derived phenotypes (IDPs) from different magnetic resonance imaging (MRI) modalities remains unclear. Here, we use two-sample Mendelian randomization (MR) to systematically assess the causal relationships between TL and 3,935 brain IDPs. Overall, the MR results suggested that TL was causally associated with 193 IDPs with majority representing diffusion metrics in white matter tracts. 68 IDPs were negatively associated with TL indicating that longer TL causes decreasing in these IDPs, while the other 125 were associated positively (longer TL leads to increased IDPs measures). Among them, ten IDPs have been previously reported as informative biomarkers to estimate brain age. However, the effect direction between TL and IDPs did not reflect the observed direction between aging and IDPs: longer TL was associated with decreases in fractional anisotropy and increases in axial, radial and mean diffusivity. For instance, TL was positively associated with radial diffusivity in the left perihippocampal cingulum tract and with mean diffusivity in right perihippocampal cingulum tract. Our results revealed a causal role of TL on white matter integrity which makes it a valuable factor to be considered when brain age is estimated and investigated.

Microstructural MRI Correlates of Cognitive Impairment in Multiple Sclerosis: The Role of Deep Gray Matter.

Although cognitive impairment (CI) is frequently observed in people with multiple sclerosis (pwMS), its pathogenesis is still controversial. Conflicting results emerged concerning the role of microstructural gray matter (GM) damage especially when involving the deep GM structures. In this study, we aimed at evaluating whether differences in cortical and deep GM structures between apparently cognitively normal (ACN) and CI pwMS (36 subjects in total) are present, using an extensive set of diffusion MRI (dMRI) indices and conventional morphometry measures. The results revealed increased anisotropy and restriction over several deep GM structures in CI compared with ACN pwMS, while no changes in volume were present in the same areas. Conversely, reduced anisotropy/restriction values were detected in cortical regions, mostly the pericalcarine cortex and precuneus, combined with reduced thickness of the superior frontal gyrus and insula. Most of the dMRI metrics but none of the morphometric indices correlated with the Symbol Digit Modality Test. These results suggest that deep GM microstructural damage can be a strong anatomical substrate of CI in pwMS and might allow identifying pwMS at higher risk of developing CI.

Convergent and Discriminant Validity of Default Mode Network and Limbic Network Perfusion in Amnestic Mild Cognitive Impairment Patients.

BACKGROUND: Previous studies reported default mode network (DMN) and limbic network (LIN) brain perfusion deficits in patients with amnestic mild cognitive impairment (aMCI), frequently a prodromal stage of Alzheimer's disease (AD). However, the validity of these measures as AD markers has not yet been tested using MRI arterial spin labeling (ASL). OBJECTIVE: To investigate the convergent and discriminant validity of DMN and LIN perfusion in aMCI. METHODS: We collected core AD markers (amyloid-ß 42 [Aß42], phosphorylated tau 181 levels in cerebrospinal fluid [CSF]), neurodegenerative (hippocampal volumes and CSF total tau), vascular (white matter hyperintensities), genetic (apolipoprotein E [APOE] status), and cognitive features (memory functioning on Paired Associate Learning test [PAL]) in 14 aMCI patients. Cerebral blood flow (CBF) was extracted from DMN and LIN using ASL and correlated with AD features to assess convergent validity. Discriminant validity was assessed carrying out the same analysis with AD-unrelated features, i.e., somatomotor and visual networks' perfusion, cerebellar volume, and processing speed. RESULTS: Perfusion was reduced in the DMN (F = 5.486, p = 0.039) and LIN (F = 12.678, p = 0.004) in APOE ɛ4 carriers compared to non-carriers. LIN perfusion correlated with CSF Aß42 levels (r = 0.678, p = 0.022) and memory impairment (PAL, number of errors, r = -0.779, p = 0.002). No significant correlation was detected with tau, neurodegeneration, and vascular features, nor with AD-unrelated features. CONCLUSION: Our results support the validity of DMN and LIN ASL perfusion as AD markers in aMCI, indicating a significant correlation between CBF and amyloidosis, APOE ɛ4, and memory impairment.

Noise removal in resting-state and task fMRI: functional connectivity and activation maps.

OBJECTIVE: Blood-oxygenated-level dependent (BOLD)-based functional magnetic resonance imaging (fMRI) is a widely used non-invasive tool for mapping brain function and connectivity. However, the BOLD signal is highly affected by non-neuronal contributions arising from head motion, physiological noise and scanner artefacts. Therefore, it is necessary to recover the signal of interest from the other noise-related fluctuations to obtain reliable functional connectivity (FC) results. Several pre-processing pipelines have been developed, mainly based on nuisance regression and independent component analysis (ICA). The aim of this work was to investigate the impact of seven widely used denoising methods on both resting-state and task fMRI. APPROACH: Task fMRI can provide some ground truth given that the task administered has well established brain activations. The resulting cleaned data were compared using a wide range of measures: motion evaluation and data quality, resting-state networks and task activations, FC. MAIN RESULTS: Improved signal quality and reduced motion artefacts were obtained with all advanced pipelines, compared to the minimally pre-processed data. Larger variability was observed in the case of brain activation and FC estimates, with ICA-based pipelines generally achieving more reliable and accurate results. SIGNIFICANCE: This work provides an evidence-based reference for investigators to choose the most appropriate method for their study and data.

Usefulness of High Resolution T2-Weighted Images in the Evaluation and Surveillance of Vestibular Schwannomas? Is Gadolinium Needed?

OBJECTIVES: First aim of study was to compare the diagnostic accuracy of high resolution T2-WI (HRT2-WI) and gadolinium-enhanced T1-weighted image (Gd T1-WI) sequences in quantitative evaluation of vestibular schwannomas (VS). The second one was to determine through qualitative evaluation when Gd-injection should be recommended. METHODS: Two observers in consensus retrospectively reviewed 137 magnetic resonance imaging (MRI) scans of patients with histological diagnosis of VS (33 women, 26 men), 116 with both HRT2-WI and Gd T1-WI. The examinations were subdivided in: surveillance (10), pretreatment (30), posttreatment (43), and posttreatment follow-up (33) studies. Quantitative evaluation was based on size measurement of the lesion. Structural details of the lesion, facial nerve course, and involvement of the fundus of the internal auditory canal were assessed for the qualitative evaluation in both sequences. RESULTS: No statistically significant changes were demonstrated between size measurement in the HRT2-WI and Gd T1-WI (p = 0.329). Sensitivity, specificity, and accuracy of HRT2-WI in the detection of lesional size were 90.4, 98.9, 92.5%, respectively. HRT2-WI was worse for characterization of structural details in pretreatment, posttreatment, and posttreatment follow-up examinations. HRT2-WI showed better or equal capability for all the groups in the demonstration of the facial nerve course. In the evaluation of the involvement of the fundus of the internal auditory canal, HRT2-WI showed worse results or had a complementary role for the posttreatment and posttreatment follow-up groups, while in the other groups was considered mainly equal or complementary. Only in the surveillance group, HRT2-WI was never worse for all the criteria. CONCLUSION: Results for quantitative evaluation were similar in both the sequences. Gadolinium injection can be avoided only in surveillance studies.

Non-Aß-Dependent Factors Associated with Global Cognitive and Physical Function in Alzheimer's Disease: A Pilot Multivariate Analysis.

Recent literature highlights the importance of identifying factors associated with mild cognitive impairment (MCI) and Alzheimer's Disease (AD). Actual validated biomarkers include neuroimaging and cerebrospinal fluid assessments; however, we investigated non-Aß-dependent factors associated with dementia in 12 MCI and 30 AD patients. Patients were assessed for global cognitive function (Mini-Mental state examination-MMSE), physical function (Physical Performance Test-PPT), exercise capacity (6-min walking test-6MWT), maximal oxygen uptake (VO2max), brain volume, vascular function (flow-mediated dilation-FMD), inflammatory status (tumor necrosis factor-α ,TNF- α, interleukin-6, -10 and -15) and neurotrophin receptors (p75NTR and Tropomyosin receptor kinase A -TrkA). Baseline multifactorial information was submitted to two separate backward stepwise regression analyses to identify the variables associated with cognitive and physical decline in demented patients. A multivariate regression was then applied to verify the stepwise regression. The results indicated that the combination of 6MWT and VO2max was associated with both global cognitive and physical function (MMSE = 11.384 + (0.00599 × 6MWT) - (0.235 × VO2max)); (PPT = 1.848 + (0.0264 × 6MWT) + (19.693 × VO2max)). These results may offer important information that might help to identify specific targets for therapeutic strategies (NIH Clinical trial identification number NCT03034746).

Age-Specific 18F-FDG Image Processing Pipelines and Analysis Are Essential for Individual Mapping of Seizure Foci in Pediatric Patients with Intractable Epilepsy.

18F-FDG PET is an important tool for the presurgical assessment of children with drug-resistant epilepsy. Standard assessment is performed visually and is often subjective and highly user-dependent. Voxelwise statistics can be used to remove user-dependent biases by automatically identifying areas of significant hypo- or hypermetabolism associated with the epileptogenic area. In the clinical setting, this analysis is performed using commercially available software. These software packages suffer from two main limitations when applied to pediatric PET data: pediatric scans are spatially normalized to an adult standard template, and statistical comparisons use an adult control dataset. The aim of this work was to provide a reliable observer-independent pipeline for the analysis of pediatric 18F-FDG PET scans, as part of presurgical planning in epilepsy. Methods: A pseudocontrol dataset (19 subjects 6-9 y old, and 93 subjects 10-20 y old) was used to create two age-specific 18F-FDG PET pediatric templates in standard pediatric space. The 18F-FDG PET scans of 46 epilepsy patients (16 patients 6-9 y old, and 30 patients 10-17 y old) were retrospectively collated and analyzed using voxelwise statistics. This procedure was implemented with the standard pipeline available in the commercial software Scenium and an in-house Statistical Parametric Mapping, version 8 (SPM8), pipeline (including age-specific pediatric templates and reference database). A κ-test was used to assess the level of agreement between the findings of voxelwise analyses and the clinical diagnosis of each patient. The SPM8 pipeline was further validated using postsurgical seizure-free patients. Results: Improved agreement with the clinical diagnosis was reported using SPM8, in terms of focus localization, especially for the younger patient group: κ = 0.489 for Scenium versus 0.826 for SPM. The proposed pipeline also showed a sensitivity of about 70% in both age ranges for the localization of hypometabolic areas on pediatric 18F-FDG PET scans in postsurgical seizure-free patients. Conclusion: We showed that by creating age-specific templates and using pediatric control databases, our pipeline provides an accurate and sensitive semiquantitative method for assessing the 18F-FDG PET scans of patients under 18 y old.

Dual-echo ASL based assessment of motor networks: a feasibility study.

OBJECTIVE: Dual-echo arterial spin labeling (DE-ASL) technique has been recently proposed for the simultaneous acquisition of ASL and blood-oxygenation-level-dependent (BOLD)-functional magnetic resonance imaging (fMRI) data. The assessment of this technique in detecting functional connectivity at rest or during motor and motor imagery tasks is still unexplored both per-se and in comparison with conventional methods. The purpose is to quantify the sensitivity of the DE-ASL sequence with respect to the conventional fMRI sequence (cvBOLD) in detecting brain activations, and to assess and compare the relevance of node features in decoding the network structure. APPROACH: Thirteen volunteers were scanned acquiring a pseudo-continuous DE-ASL sequence from which the concomitant BOLD (ccBOLD) simultaneously to the ASL can be extracted. The approach consists of two steps: (i) model-based analyses for assessing brain activations at individual and group levels, followed by statistical analysis for comparing the activation elicited by the three sequences under two conditions (motor and motor imagery), respectively; (ii) brain connectivity graph-theoretical analysis for assessing and comparing the network models properties. MAIN RESULTS: Our results suggest that cvBOLD and ccBOLD have comparable sensitivity in detecting the regions involved in the active task, whereas ASL offers a higher degree of co-localization with smaller activation volumes. The connectivity results and the comparative analysis of node features across sequences revealed that there are no strong changes between rest and tasks and that the differences between the sequences are limited to few connections. SIGNIFICANCE: Considering the comparable sensitivity of the ccBOLD and cvBOLD sequences in detecting activated brain regions, the results demonstrate that DE-ASL can be successfully applied in functional studies allowing to obtain both ASL and BOLD information within a single sequence. Further, DE-ASL is a powerful technique for research and clinical applications allowing to perform quantitative comparisons as well as to characterize functional connectivity.

Exploring the Epileptic Brain Network Using Time-Variant Effective Connectivity and Graph Theory.

The application of time-varying measures of causality between source time series can be very informative to elucidate the direction of communication among the regions of an epileptic brain. The aim of the study was to identify the dynamic patterns of epileptic networks in focal epilepsy by applying multivariate adaptive directed transfer function (ADTF) analysis and graph theory to high-density electroencephalographic recordings. The cortical network was modeled after source reconstruction and topology modulations were detected during interictal spikes. First a distributed linear inverse solution, constrained to the individual grey matter, was applied to the averaged spikes and the mean source activity over 112 regions, as identified by the Harvard-Oxford Atlas, was calculated. Then, the ADTF, a dynamic measure of causality, was used to quantify the connectivity strength between pairs of regions acting as nodes in the graph, and the measure of node centrality was derived. The proposed analysis was effective in detecting the focal regions as well as in characterizing the dynamics of the spike propagation, providing evidence of the fact that the node centrality is a reliable feature for the identification of the epileptogenic zones. Validation was performed by multimodal analysis as well as from surgical outcomes. In conclusion, the time-variant connectivity analysis applied to the epileptic patients can distinguish the generator of the abnormal activity from the propagation spread and identify the connectivity pattern over time.

Ensemble average propagator-based detection of microstructural alterations after stroke.

PURPOSE: New analytical reconstruction techniques of diffusion weighted signal have been proposed. A previous work evidenced the exploitability of some indices derived from the simple harmonic oscillator-based reconstruction and estimation (3D-SHORE) model as numerical biomarkers of neural plasticity after stroke. Here, the analysis is extended to two additional indices: return to the plane/origin (RTPP/RTOP) probabilities. Moreover, several motor networks were introduced and the results were analyzed at different time scales. METHODS: Ten patients underwent three diffusion spectrum imaging (DSI) scans [1 week (tp1), 1 month (tp2) and 6 months (tp3) after stroke]. Ten matched controls underwent two DSI scans 1 month apart. 3D-SHORE was used for reconstructing the signal and the microstructural indices were derived. Tract-based analysis was performed along motor cortical, subcortical and transcallosal networks in the contralesional area. RESULTS: The optimal intra-class correlation coefficient (ICC) was obtained in the subcortical loop for propagator anisotropy (ICC [Formula: see text] 0.96), followed by generalized fractional anisotropy (ICC [Formula: see text] 0.94). The new indices reached the highest stability in the transcallosal network and performed well in the cortical and subcortical networks with the exception of RTOP in the cortical loop (ICC [Formula: see text] 0.59). They allowed discriminating patients from controls at the majority of the timescales. Finally, the regression model using indices calculated along the subcortical loop at tp1 resulted in the best prediction of clinical outcome. CONCLUSIONS: The whole set of microstructural indices provide measurements featuring high precision. The new indices allow discriminating patients from controls in all networks, except for RTPP in the cortical loop. Moreover, the 3D-SHORE indices in subcortical connections constitute a good regression model for predicting the clinical outcome at 6 months, supporting their suitability as numerical biomarkers for neuronal plasticity after stroke.

Reproducibility of EEG-fMRI results in a patient with fixation-off sensitivity.

Blood oxygenation level-dependent (BOLD) activation associated with interictal epileptiform discharges in a patient with fixation-off sensitivity (FOS) was studied using a combined electroencephalography-functional magnetic resonance imaging (EEG-fMRI) technique. An automatic approach for combined EEG-fMRI analysis and a subject-specific hemodynamic response function was used to improve general linear model analysis of the fMRI data. The EEG showed the typical features of FOS, with continuous epileptiform discharges during elimination of central vision by eye opening and closing and fixation; modification of this pattern was clearly visible and recognizable. During all 3 recording sessions EEG-fMRI activations indicated a BOLD signal decrease related to epileptiform activity in the parietal areas. This study can further our understanding of this EEG phenomenon and can provide some insight into the reliability of the EEG-fMRI technique in localizing the irritative zone.