Evaluation of White Matter Integrity Utilizing the DELPHI (TMS-EEG) System

Abstract

Objective: Evaluation of brain white matter fibers connectivity damage in stroke and traumatic brain injury (TBI) subjects by direct electrophysiological imaging (DELPHI) which analyzes TMS (transcranial magnetic stimulation) evoked potentials. Methods: Study included 123 participants out of which: 53 subjects with white matter related pathologies of (39 stroke, 14 TBI) and 70 healthy age-related controls. All subjects underwent DELPHI brain network evaluations of TMS-EEG evoked potentials, and DTI scans for quantification of white matter microstructure fractional anisotropy. Results: DELPHI output measures show a significant difference between healthy and stroke\TBI groups, a multidimensional approach was able to classify healthy from unhealthy with a balanced accuracy of 0.81±0.02 and AUC of 0.88±0.01 . Moreover, a multivariant regression model of DELPHI output measures achieved prediction of white matter microstructure changes measured by fractional anisotropy (FA) with highest correlations observed for fibers proximal to stimulation area such as frontal corpus callosum (r=0.7±0.02) ,anterior internal capsule (r=0.7±0.02) and fronto-occipital fasciculus (r=0.65±0.03). Conclusions: These results indicate that features of TMS evoked response are correlated to white matter microstructure changes observed in pathological conditions such as stroke and TBI, and that a multi-dimensional approach combining these features in supervised learning methods serve as a strong indicator for abnormalities and changes in white matter integrity.

Introduction

The disruption of normal patterns of structural brain connectivity is believed to play a central role in the pathophysiology of many neurological and psychiatric disorders, such as dementia, movement disorders, stroke, traumatic brain injury (TBI), etc. White matter (WM) pathways consist of myelinated axonal structures that constitute the connectivity between different brain regions. Axonal injury and degeneration may occur even in the absence of tissue disruption. Therefore, patients may experience significant impairment despite the absence of abnormal findings on conventional CT or MRI. Moreover, focal imaging abnormalities that can be detected using CT and MRI are poor predictors of outcome (Niogi and Mukherjee, 2010; Adalbert and Coleman, 2013; Wang et al., 2016; Salvadores et al., 2017).

Axonal injury is a key determinant of clinical outcome in cases of brain injury and has been shown to be an important factor determining long-term motor, cognitive, and neuropsychiatric disability following brain injury (Mac Donald et al., 2007; Di Paola et al., 2010; Moura et al., 2019). Diagnostic tests that can discriminate significant axonal injury and degeneration are needed in order to accurately assess injury severity and effectively determine treatment and follow-up pathway.

WM pathways determined through diffusion tensor imaging (DTI) are traditionally considered to be the biophysical representation of axonal bundles and their myelin sheets. DTI streamlines between cortical and subcortical gray matter (GM) regions of interest (ROIs) can be used as a measure of the magnitude and strength of connection between ROIs (Basser et al., 1994; Le Bihan et al., 2001). Analyzing changes in brain connectivity using DTI tractography is widely used to evaluate axonal injury that is a hallmark of stroke and TBI (Arfanakis et al., 2002; Lo et al., 2010), in which WM tracks can be injured directly or indirectly through Wallerian degeneration (the anterograde distal degeneration of injured axons accompanied by demyelination) (Chen et al., 2017; Moura et al., 2019). DTI metrics are used to address specific or diffused WM damage that is frequently affected by stroke lesions in order to describe stroke focality and severity and predict rehabilitation potential (Puig et al., 2013; Puig et al., 2020). In TBI patients, a consistent reduction in fractional anisotropy (FA) has been typically found in areas affected by traumatic axonal injury (TAI). These regions include the subcortical WM of the frontal and temporal regions, the splenium of the corpus callosum, the posterior limb of the internal capsule, and the cerebral peduncles (Nakayama et al., 2006; Benson et al., 2007; Xu et al., 2007; Greenberg et al., 2008).

Combining transcranial magnetic stimulation (TMS) with electroencephalogram (EEG) is extensively used to study and assess cerebral reactivity and connectivity (Rossini et al., 2015; Ferreri et al., 2017; Tremblay et al., 2019). TMS is a non-invasive brain stimulation method that allows the study of human cortical function in vivo (Ilmoniemi et al., 1999; Hallett, 2007). TMS enables the modulation and exploration of brain functional status. Studies integrating TMS with EEG (TMS-EEG) have shown that TMS produces waves of activity that reverberate throughout the cortex and that are reproducible and reliable (Casarotto et al., 2010; Farzan et al., 2010; Kerwin et al., 2018), thus providing direct information about cortical excitability and connectivity with excellent temporal resolution (Ilmoniemi et al., 1997; Shafi et al., 2012; Chung et al., 2015; Pennisi et al., 2016; Bordet et al., 2017). TMS-EEG has been used to causally probe the dynamic effective connectivity of human brain networks (Rogasch and Fitzgerald, 2013; Ferreri et al., 2014; Kugiumtzis and Kimiskidis, 2015).

Direct electrophysiological imaging (DELPHI) is a new clinical methodology for evaluating fundamental physiological properties of brain network function by automatically combining TMS stimulation and EEG monitoring (Zifman et al., 2019). The DELPHI software algorithm extracts direct stimulation-related properties of brain network function, characterizing a profile of brain functional pathophysiology including properties of network integrity and plasticity.

Several studies have explored the correlation of TMS-evoked response to WM fibers and GM changes, indicating that this method may be useful for probing neurological disease and conditions that manifest changes in white or GM ROI’s (Kirton et al., 2010; Lanza et al., 2013; Pennisi et al., 2015).

The current study tests DELPHI technology and output parameters as a clinical tool for evaluation of specific neuronal tracks connectivity and displays its strong correlation to structural changes in WM fibers. The study population includes stroke and TBI patients that represent the population with connectivity disruption caused predominantly by WM changes and healthy controls.

Materials and Methods

Clinical Data Collection and Analysis

The study was carried out in accordance with the recommendation of the Shamir Medical center review board. The protocol was approved by the local institutional “Ethical Committee” as a retrospective study of data. All participants underwent the exact same MRI scan and DELPHI evaluation protocols.

Study Population

The study focused on two main groups: (a) healthy controls (HC) and (b) subjects diagnosed with WM injury-related conditions either (b1) ischemic stroke (cortical and subcortical) or (b2) TBI. All stroke/TBI subjects were considered to be at post-sub-acute stage of the injury, post-injury rehabilitation, and discharged from hospitalization. All subjects included in the study had updated documented medical history, demographics, clinical MRI evaluation, DTI scan, and DELPHI (TMS-EEG) scan performed within 2 weeks of each other. Healthy subject group inclusion criteria (a) were as follows: (1) age over 18 years, (2) no neurological or psychiatric disorder documented in medical history or self-report, (3) absence of any significant abnormal findings in clinical MRI evaluation, such as brain tumors, subdural hematoma, and other brain structural lesions related to diagnosed brain disease other than common age-related changes, and 4. no psychoactive or other brain directed medications. Stroke subject group (b1) inclusion criteria were as follows: (1) age over 18 years, (2) history of ischemic stroke over 6 months and less than 5 years prior to MRI, DTI, and TMS-EEG scans (post-sub-acute stage and rehabilitation), (3) injury was detectible by clinical MRI scan, and (4) no other psychiatric or neurological comorbidities. TBI subject group (b2) inclusion criteria were as follows: (1) age over 18 years, (2) history of TBI over 6 months and less than 5 years prior to MRI, DTI, and TMS-EEG scans (post-sub-acute stage and rehabilitation), (3) injury was detectible by clinical MRI scan, and 4. no other psychiatric or neurological comorbidities. Demographic information of the study groups is listed in Tables 1, 2. Due to the age differences of the TBI and stroke populations, the healthy population was divided into age-matched groups to serve as a valid control group for each of the WM conditions populations.

TABLE 1

Table 1. (A) Demographics of stroke and age-matched HC subjects. Age distribution was not significantly different (p > 0.05). (B) Overview of general injury location detected by clinical MRI.

TABLE 2

Table 2. (A) Demographics of TBI and HC subjects. Age distribution was not significantly different (p > 0.05). (B) Overview of general injury location detected by clinical MRI.

DTI Imaging

All subjects underwent brain MRI scan and DELPHI evaluation not more than 2 weeks apart. Imaging was performed with a 3 Tesla system (MAGNETOM Skyra; Siemens Healthineers, Erlangen, Germany) with 20 channels receiver head coil. Thirty diffusion-weighted images were scanned with different gradient directions (b = 1,000) and one volume without diffusion weighting, with the following parameters: TR = 10,300 ms, TE = 89 ms, voxel size = 1.8 × 1.8 mm, matrix = 128 × 128, no. of slices = 63, slice thickness = 2.2 mm.

DTI Analysis

DTI analysis was performed on the FA map calculated by Siemens post-processing software. For each subject, the WM atlas (ICBM-MORI white matter atlas (Mori et al., 2008) was registered to the FA map using SPM (version 12; The Wellcome Centre for Human Neuroimaging, UCL Queen Square Institute of Neurology, London, United Kingdom) and manually validated. Mean values for FA higher than 0.2 were calculated in WM regions according to the atlas (Wakana et al., 2004, 2007; Hua et al., 2008).

TMS-EEG

TMS was performed with a MagPro R30 stimulator (MagVenture, Denmark) and an MCF-B65-HO figure-8 Coil (MagVenture, Denmark). 32-Channel EEG data were obtained using TMS compatible BrainAmp DC amplifier (5 kHz sampling rate, ± 16.384 mV measurement range, analog low pass filter 1 kHz; Brain Products GmbH, Germany). These were connected to the waveguard^TM EEG cap (ANT Neuro, Netherland) with Ag-AgCl electrodes. The reference and ground electrodes were affixed to the ear lobes. EEG data were recorded using a BrainVision Recorder software (Brain Products GmbH, Germany).

Experimental Procedure

TMS coil was positioned over the left cortical motor (M1) region, at 45° toward the contralateral forehead according to the guidelines (Rossini et al., 2015). Single pulse (<0.3 Hz frequency and 1 Hz inhibitory frequency) stimulation was performed at 80% from rest motor threshold (RMT) intensity. Data acquisition, preprocessing, and cleaning of the TMS-evoked response include rejection of bad channels and epochs containing large artifacts, followed by bandpass FIR filter (0.5–45 Hz) as detailed by Zifman et al. (2019). A thin (0.5 mm) foam pad was attached to the TMS coil to minimize electrode movement and bone-conducted auditory artifact. Participants were instructed to keep their eyes closed throughout the examination to reduce ocular artifacts. The operator of the system conversed with subjects between the short stimulation protocol blocks in order to avoid drowsiness. Electrodes data were grouped to regional recording hotspots for analysis and statistical purposes: frontal: F3, F5—ipsilateral and F4, F6—contralateral to stimulation; parietal: C3, C5, CP1—ipsilateral and C4, C6, CP2—contralateral to stimulation; temporal CP5, CP3, CF5—ipsilateral and CP6, CP4, FC6—contralateral to stimulation; and occipital cortex: O1, PO3—ipsilateral and O2, PO4—contralateral to stimulation.

DELPHI Physiological Network Profile Analysis Parameters

DELPHI analyzes the regional and network TMS-evoked potential (TEP) measured as EEG patterns of single and history-dependent events as described by Zifman et al. (2019). TEP response refers to 20–300 ms after TMS stimulation and was described by four parameters: (1) the slope of early response, between 60 and 100 ms, referred to as early phase deflection (EPD); (2) the slope of late response, between 100 and 180 ms, referred to as late phase deflection (LPD); (3) waveform adherence (WFA) refers to the entire TEP response adherence to mean healthy age-related signal collected from previous studies; and (4) the normalized difference ratio between mean field potential of single pulse response (MFP_single) and mean field potential of inhibitory frequency of stimulation (MFP_i), indicating network short-term plasticity (STP) calculated as (MFP_single − MFP_i)/(MFP_single + MFP_i). All data processing and feature extraction were performed automatically by the DELPHI software algorithm.

Classification

Linear Support Vector Machine (SVM) classifier, with 50 permutations of stratified 5-fold cross-validation and class prior probabilities set to reflect realistic balanced proportions, was used for classification of population with WM-related conditions (TBI and stroke) and age-matched healthy population by DELPHI output measures as features vector. Linear SVMs are known as stable classifiers with low complexity that enables accommodation to outliers and are relatively insensitive to overtraining and curse-of-dimensionality (Lotte et al., 2007, 2018; Li et al., 2014). In order to achieve maximal signal-to-noise ratio (SNR) and decrease the feature vector length, only electrode hotspots that are closest to the stimulation site were used in the model construction (left/right temporal and parietal hotspots) as they are known to be the most reliable and are traditionally used in TMS-EEG studies (Lioumis et al., 2009; Ilmoniemi and Kicić, 2010). Feature vector included all four DELPHI output measures. After training the model on the train data set, the model’s classification performance was assessed on the test data set by extracting balanced accuracy, sensitivity (true positive rate), specificity (true negative rate), mean, and standard deviation (STDV) values over all permutations. Receiver Operating Characteristic (ROC) curve was constructed by plotting the mean true positive rate (sensitivity) against the mean false positive rate (1 - specificity) with the various cut-off thresholds.

Regression

WM ROIs FA was calculated for major fiber tracks, and a multivariate linear regression model was created to predict ROI’s FA high-dimensional data using DELPHI features vector (described above) in all study populations sets pooled together (healthy, stroke, and TBI populations) (Healy, 1980). The model was used with 50 permutations of 5-fold cross-validation, and prediction performance was assessed by Pearson’s correlation coefficients and root-mean-square error (RMSE), mean, and STDV over all permutations.

Statistical Analysis

Statistical data analysis to account for the differences between the groups was performed by two samples t-test, followed by the Bonferroni–Holm method for multiple comparisons correction. Corrected p-values are depicted in text, and the significance of the results following multiple comparisons correction is depicted by asterisks: ^∗p < 0.05, ^∗∗p < 0.01, ^∗∗∗p < 0.001, ns—non-significant. All data statistical analysis was made using MATLAB (R2020a; The Mathworks Inc., MA, United States).

Results

The current study focuses on TMS-EEG neurophysiological insight into brain network function and connectivity by correlating TMS-evoked EEG response to WM fibers integrity. For this reason, a heterogenic group of subjects who sustained structural connectivity injuries in potentially different areas were randomly selected. The study groups include post-rehabilitation, ischemic stroke subjects with lesions of different sizes, post-rehabilitation TBI subjects, and healthy subjects. Tables 1A,B, 2A,B summarize age, gender, hand dominance, and general brain region of injury.

Evaluation of WM Changes in Stroke/TBI

ROI-based comparison of WM tracks FA displayed, as expected, significantly reduced mean FA and elevated STDV for the majority of measured WM fibers both in the stroke group compared with HC and in the TBI group compared with age-matched HC (Table 3). The significant reduction of FA values in a wide range of ROIs is consistent with the heterogenicity of study population injuries and variability of fiber damage locations and sizes in both stroke and TBI.

TABLE 3

Table 3. ROI-based FA analysis of stroke and TBI compared with age-matched HC.

DELPHI Network Functional Analysis of Healthy vs. Stroke/TBI

Evaluating network function using DELPHI output parameters revealed a significant difference between the study groups. Significant differences were observed between stroke and age-matched HC for WFA in all regional hotspots (Figure 1A). EPD, which describes the early phase of response, did not display significant differences between the groups (Figure 1B); however, LPD displayed a significant difference for parietal hotspots (Figure 1C). STP index did not display significant differences (Figure 1D). DELPHI evaluation of TBI population compared with age-matched HC revealed a significant decrease in WFA in all recording hotspots (Figure 2A). EPD displayed significant differences in the left temporal and parietal areas, proximal to the stimulation area (Figure 2B). LPD and STP did not display any significant change (Figures 2C,D). A multidimensional approach combining DELPHI output measures using a linear SVM classifier is depicted in the model’s cross-validated classification performance. DELPHI TMS-EEG output parameters were able to differentiate HC from stroke subjects, at balanced accuracy, sensitivity, and specificity rates of 0.81 ± 0.03, 0.83 ± 0.05, and 0.8 ± 0.04, respectively. ROC area under the curve (AUC) was 0.87 ± 0.02 (Figure 3A). For TBI and HC, AUC was 0.87 ± 0.03, and classification balanced accuracy, sensitivity, and specificity rates were 0.82 ± 0.05, 0.83 ± 0.03, and 0.81 ± 0.08, respectively (Figure 3B). Pooling together all WM injury subjects (stroke and TBI) displayed a classification balanced accuracy of 0.81 ± 0.02, sensitivity rate of 0.82 ± 0.03, specificity of 0.8 ± 0.03, and AUC of 0.88 ± 0.01.

FIGURE 1

Figure 1. DELPHI output measures of HC vs. the stroke group at different recording sites. (A) WFA, (B) EPD, (C) LPD, (D) STP. Fl, frontal left; Fr, frontal right; Pl, parietal left; Pr, parietal right; Tl, temporal left; Tr, temporal right; Ol, occipital left; Or, occipital right. Mean values (thick lines) and SE (narrow lines) are depicted. Significant differences between populations are indicated by asterisks.

FIGURE 2

Figure 2. DELPHI output measures of HC vs. the TBI group at different recording sites. (A) WFA, (B) EPD, (C) LPD, (D) STP. Fl, frontal left; Fr, frontal right; Pl, parietal left; Pr, parietal right; Tl, temporal left; Tr, temporal right; Ol, occipital left; Or, occipital right. Mean values (thick lines) and SE (narrow lines) are depicted. Significant differences between populations are indicated by asterisks.

FIGURE 3

Figure 3. Classification performance of WM-related conditions using DELPHI output measures. Classification of stroke (A) and of TBI (B) and age-matched HC using DELPHI output measures as feature vector. Classification performance is depicted by the ROC curve and standard statistical measures, and balanced accuracy (accuracy), sensitivity, specificity, and AUC values are detailed in the figure.

Prediction of WM Abnormalities Using DELPHI

Next, we evaluated the capacity of DELPHI (TMS-EEG) output measures to predict WM-related connectivity abnormalities captured by DTI–FA values, using a multivariate linear regression model. All study groups were pooled together (HC, stroke, and TBI), and WM ROIs FA was calculated for major fiber tracks. Cross-validated results of the multivariate model exhibit high correlation between DELPHI output parameters and DTI–FA values of WM fibers. Highest correlations were observed for left and frontal tracts (Table 4), most significantly for frontal commissural fiber—corpus callosum (r = 0.7 ± 0.02, RMSE = 0.07) (Figure 4A), left projection fiber—anterior corona radiata (r = 0.66 ± 0.02, RMSE = 0.04) (Figure 4B), and left association fiber—superior fronto-occipital fasciculus (r = 0.65 ± 0.02, RMSE = 0.09) (Figure 4C).

TABLE 4

Table 4. WM ROIs FA values prediction by multivariate linear regression model of DELPHI output measures.

FIGURE 4

Figure 4. Prediction of WM abnormalities using DELPHI multivariate linear regression model on the pooled study groups collected data. Predicted FA values (Y-axis) were plotted against the true FA measured values (X-axis) for exemplary fibers. (A) Frontal commissural fiber—corpus callosum. (B) Left projection fiber—anterior corona radiata. (C) Left association fiber—superior fronto-occipital fasciculus.

Discussion

In the present study, we evaluated the performance of DELPHI output parameters as biomarkers for structural and functional brain connectivity abnormalities, captured by WM tracks integrity measures, in stroke and TBI patients. Study results demonstrate the changes between HC and TBI or stroke populations for DELPHI output measures. The most significant difference of both stroke and TBI patients from healthy age-matched controls was demonstrated by WFA that represents the general waveform of single pulse TEP response. Both EPD and LPD were most significantly changed between the HC and TBI/stroke groups in parietal and temporal hotspots, most probably due the higher SNR in the hotspots that are proximal to the stimulation area or contralateral to it. LPD was significantly changed in the stroke group alone, whereas EPD was significantly different in TBI alone. This may originate from the differences in age groups, as the TBI group and control are significantly younger, or might possibly arise from different natures of injuries that for stroke might be more focal and for TBI more diffused. These differences should be further explored and may potentially provide additional insight. The fact that WFA alone was significantly affected in both the TBI and stroke groups for all recording hotspots may indicate the notion that different injuries (location, sizes, etc.) are represented in different phases or latencies (late/early) of the TEP response and therefore were not as significantly noticeable in the EPD and LPD but rather in the WFA parameter that represents the entire response pattern capturing all possible changed parameters of the signal. DTI analysis demonstrated that, as expected, WM integrity was significantly compromised in both the TBI and stroke groups compared with HC in most ROIs. This is most probably due to the variability in injury locations and sizes. A multidimensional approach utilizing linear SVM model, which incorporates all DELPHI output measures, displays high classification performance of healthy vs. TBI or stroke populations with a balanced accuracy of 0.81 ± 0.02 and AUC of 0.88 ± 0.01. Moreover, a multivariate linear regression model revealed high correlation between DELPHI output measures and various WM fibers, which is consistent with the decrease of FA values in a wide range of ROIs in the stroke and TBI groups, compared with HC. Interestingly, higher sensitivity to WM integrity was observed in fibers that are in closest proximity to TMS stimulation site (M1L), namely, the highest values were recorded in the frontal and left fibers, which are proximal to the TMS stimulation site. These data support DELPHI (TMS-EEG) as a tool for neurophysiological evaluation of network connectivity and its potential as an affordable, safe and available tool for monitoring and evaluating brain network structural and functional abnormalities in specific neuronal circuits. These results are supported by previous studies (Manganotti et al., 2015; Lanza et al., 2017; Hordacre et al., 2019) and indicate that in order to get a more comprehensive and localized evaluation of functional and structural brain connectivity, integrity, and plasticity, at least two contralateral stimulation sites are required. Moreover, direct TMS stimulation of frontal, medial, and posterior areas may increase TMS-EEG sensitivity to both structural and functional ROIs. TMS-EEG technology has been shown to provide insight into the evaluation and monitoring of functional effective connectivity in various brain disorders, such as Alzheimer’s disease (AD), stroke, and TBI prognosis and rehabilitation prediction (Bortoletto et al., 2015; Manganotti et al., 2015; Tremblay et al., 2019). Our study supports the clinical utility of TMS-EEG in brain disorders of functional and structural effective connectivity, displays high sensitivity of DELPHI neurophysiological measures to WM-related structural connectivity changes as measured with DTI, and supports specific fiber localization correlated to the stimulation area. It is important to note that the current study classifiers did not consider age-related WM and GM changes, such as minor lacunar infracts, lesions, and brain atrophy, that may affect “brain health” status in the HC group; however, these measures were addressed in the regression model with FA analysis of WM fibers for all study groups. Further studies are required to address the functional sensitivity of changes in neuronal circuits and explore the potential predictive value of this tool for the early detection, predictive recovery, and prognosis of WM-related brain abnormalities and diseases. Future studies should include populations suffering from brain abnormalities that may or may not be detected by conventional clinical MRI/CT scans in order to establish detection of subtle WM changes in clinical conditions, such as mild TBI, cerebral small vessels disease (CSVD), etc., and record clinical measures of cognitive and physiological performance. The current study supports previous studies that display TEP as a biomarker for brain excitability and functional connectivity (Bagattini et al., 2019; Opie et al., 2019; Rossini et al., 2020). Here, we provide further evidence for the physiological relevance of probing specific networks (such as the motor network that was stimulated in the current study) with TMS-EEG technology in order to record structural or functional WM abnormalities in brain injuries or other diseases. These results demonstrate the potential clinical value of DELPHI output measures in differentiating brain-related diseases and abnormalities, providing a reliable, safe and easy-to-use tool for monitoring WM-related connectivity changes associated with brain-related abnormalities, such as brain injuries and degeneration.

Conclusion

The ability to facilitate between clinical observations, structural and functional brain connectivity, and network physiology is crucial in order to achieve an optimal clinical care. The data presented in this study support the concept that TMS-EEG physiological measures provide specific accessible insight into the human connectome that may help to narrow the gap between anatomical data, effective network function, and clinical observations, thus, providing a clinical tool for monitoring brain network function and brain health.

The DELPHI automated acquisition and analysis system can be used in order to monitor brain health throughout aging and may enable the early detection of abnormal pathophysiological changes leading to neurodegeneration.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary materials, further inquiries can be directed to the corresponding author/s.

Ethics Statement

The studies involving human participants were reviewed and approved by Shamir Medical center review board. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

NZ with OL-L, HF, and ID have determined the study design and criteria. NZ conducted most of the laboratory works. OL-L and HF designed most of the DELPHI software algorithm. HF and ID contributed in designing and writing most of the manuscript and submitting the manuscript for publication. SE was principal investigator in the study. ES was responsible for the MRI imaging analysis in the study. DH, DT, SE, and ES contributed in assuring the methods and the quality of the results and with reviewing the manuscript and approving for publication of the content. All authors contributed to the article and approved the submitted version.

Conflict of Interest

OL-L, NZ, DH, DT, HF, and ID are employed by QuantalX neuroscience.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Adalbert, R., and Coleman, M. P. (2013). Review: axon pathology in age-related neurodegenerative disorders. Neuropathol. Appl. Neurobiol. 39, 90–108. doi: 10.1111/j.1365-2990.2012.01308.x

Arfanakis, K., Haughton, V. M., Carew, J. D., Rogers, B. P., Dempsey, R. J., and Meyerand, M. E. (2002). Diffusion tensor MR imaging in diffuse axonal injury. AJNR Am. J. Neuroradiol. 23, 794–802.

Bagattini, C., Mutanen, T. P., Fracassi, C., Manenti, R., Cotelli, M., Ilmoniemi, R. J., et al. (2019). Predicting Alzheimer’s disease severity by means of TMS-EEG coregistration. Neurobiol. Aging 80, 38–45. doi: 10.1016/j.neurobiolaging.2019.04.008

Basser, P. J., Mattiello, J., and LeBihan, D. (1994). MR diffusion tensor spectroscopy and imaging. Biophys. J. 66, 259–267. doi: 10.1016/s0006-3495(94)80775-1

Benson, R. R., Meda, S. A., Vasudevan, S., Kou, Z., Govindarajan, K. A., Hanks, R. A., et al. (2007). Global white matter analysis of diffusion tensor images is predictive of injury severity in traumatic brain injury. J. Neurotrauma 24, 446–459. doi: 10.1089/neu.2006.0153

Bordet, R., Ihl, R., Korczyn, A. D., Lanza, G., Jansa, J., Hoerr, R., et al. (2017). Towards the concept of disease-modifier in post-stroke or vascular cognitive impairment: a consensus report. BMC Med. 15:107. doi: 10.1186/s12916-017-0869-6

Bortoletto, M., Veniero, D., Thut, G., and Miniussi, C. (2015). The contribution of TMS-EEG coregistration in the exploration of the human cortical connectome. Neurosci. Biobehav. Rev. 49, 114–124. doi: 10.1016/j.neubiorev.2014.12.014

Casarotto, S., Romero Lauro, L. J., Bellina, V., Casali, A. G., Rosanova, M., Pigorini, A., et al. (2010). EEG responses to TMS are sensitive to changes in the perturbation parameters and repeatable over time. PLoS One. 5:e10281. doi: 10.1371/journal.pone.0010281

Chen, Y. J., Nabavizadeh, S. A., Vossough, A., Kumar, S., Loevner, L. A., and Mohan, S. (2017). Wallerian degeneration beyond the corticospinal tracts: conventional and advanced MRI findings. J. Neuroimaging 27, 272–280. doi: 10.1111/jon.12404

Chung, S. W., Rogasch, N. C., Hoy, K. E., and Fitzgerald, P. B. (2015). Measuring brain stimulation induced changes in cortical properties using TMS-EEG. Brain Stimul. 8, 1010–1020. doi: 10.1016/j.brs.2015.07.029

Di Paola, M., Di Iulio, F., Cherubini, A., Blundo, C., Casini, A. R., Sancesario, G., et al. (2010). When, where, and how the corpus callosum changes in MCI and AD: a multimodal MRI study. Neurology 74, 1136–1142. doi: 10.1212/wnl.0b013e3181d7d8cb

Farzan, F., Barr, M. S., Levinson, A. J., Chen, R., Wong, W., Fitzgerald, P. B., et al. (2010). Reliability of long-interval cortical inhibition in healthy human subjects: a TMS-EEG study. J. Neurophysiol. 104, 1339–1346. doi: 10.1152/jn.00279.2010

Ferreri, F., Guerra, A., and Rossini, P. M. (2014). Neurophysiological markers of plastic brain reorganization following central and peripheral lesions. Arch. Ital. Biol. 152, 216–238.

Ferreri, F., Vecchio, F., Guerra, A., Miraglia, F., Ponzo, D., Vollero, L., et al. (2017). Age related differences in functional synchronization of EEG activity as evaluated by means of TMS-EEG coregistrations. Neurosci. Lett. 647, 141–146. doi: 10.1016/j.neulet.2017.03.021

Greenberg, G., Mikulis, D. J., Ng, K., DeSouza, D., and Green, R. E. (2008). Use of diffusion tensor imaging to examine subacute white matter injury progression in moderate to severe traumatic brain injury. Arch. Phys. Med. Rehabil. 89(Suppl. 12), S45–S50.

Hallett, M. (2007). Transcranial magnetic stimulation: a primer. Neuron 55, 187–199. doi: 10.1016/j.neuron.2007.06.026

Healy, J. D. (1980). Maximum likelihood estimation of a multivariate linear functional relationship. J. Multivar. Anal. 10, 243–251. doi: 10.1016/0047-259x(80)90016-0

Hordacre, B., Ghosh, R., Goldsworthy, M. R., and Ridding, M. C. (2019). Transcranial magnetic stimulation-EEG biomarkers of poststroke upper-limb motor function. J. Stroke Cerebrovasc. Dis. 28:104452. doi: 10.1016/j.jstrokecerebrovasdis.2019.104452

Hua, K., Zhang, J., Wakana, S., Jiang, H., Li, X., Reich, D. S., et al. (2008). Tract probability maps in stereotaxic spaces: analyses of white matter anatomy and tract-specific quantification. Neuroimage 39, 336–347. doi: 10.1016/j.neuroimage.2007.07.053

Ilmoniemi, R. J., and Kicić, D. (2010). Methodology for combined TMS and EEG. Brain Topogr. 22, 233–248. doi: 10.1007/s10548-009-0123-4

Ilmoniemi, R. J., Ruohonen, J., and Karhu, J. (1999). Transcranial magnetic stimulation–a new tool for functional imaging of the brain. Crit. Rev. Biomed. Eng. 27, 241–284.

Ilmoniemi, R. J., Virtanen, J., Ruohonen, J., Karhu, J., Aronen, H. J., Näätänen, R., et al. (1997). Neuronal responses to magnetic stimulation reveal cortical reactivity and connectivity. Neuroreport 8, 3537–3540. doi: 10.1097/00001756-199711100-00024

Kerwin, L. J., Keller, C. J., Wu, W., Narayan, M., and Etkin, A. (2018). Test-retest reliability of transcranial magnetic stimulation EEG evoked potentials. Brain Stimul. 11, 536–544. doi: 10.1016/j.brs.2017.12.010

Kirton, A., Deveber, G., Gunraj, C., and Chen, R. (2010). Cortical excitability and interhemispheric inhibition after subcortical pediatric stroke: plastic organization and effects of rTMS. Clin. Neurophysiol. 121, 1922–1929. doi: 10.1016/j.clinph.2010.04.021

Kugiumtzis, D., and Kimiskidis, V. K. (2015). Direct causal networks for the study of transcranial magnetic stimulation effects on focal epileptiform discharges. Int. J. Neural. Syst. 25:1550006. doi: 10.1142/s0129065715500069

Lanza, G., Bella, R., Giuffrida, S., Cantone, M., Pennisi, G., Spampinato, C., et al. (2013). Preserved transcallosal inhibition to transcranial magnetic stimulation in nondemented elderly patients with leukoaraiosis. Biomed. Res. Int. 2013:351680.

Lanza, G., Bramanti, P., Cantone, M., Pennisi, M., Pennisi, G., and Bella, R. (2017). Vascular cognitive impairment through the looking glass of transcranial magnetic stimulation. Behav. Neurol. 2017:1421326.

Le Bihan, D., Mangin, J. F., Poupon, C., Clark, C. A., Pappata, S., Molko, N., et al. (2001). Diffusion tensor imaging: concepts and applications. J. Magn. Reson. Imaging 13, 534–546.

Li, X., Chen, X., Yan, Y., Wei, W., and Wang, Z. J. (2014). Classification of EEG signals using a multiple kernel learning support vector machine. Sensors (Basel). 14, 12784–12802. doi: 10.3390/s140712784

Lioumis, P., Kicić, D., Savolainen, P., Mäkelä, J. P., and Kähkönen, S. (2009). Reproducibility of TMS-Evoked EEG responses. Hum. Brain Mapp. 30, 1387–1396. doi: 10.1002/hbm.20608

Lo, C. Y., Wang, P. N., Chou, K. H., Wang, J., He, Y., and Lin, C. P. (2010). Diffusion tensor tractography reveals abnormal topological organization in structural cortical networks in Alzheimer’s disease. J. Neurosci. 30, 16876–16885. doi: 10.1523/jneurosci.4136-10.2010

Lotte, F., Bougrain, L., Cichocki, A., Clerc, M., Congedo, M., Rakotomamonjy, A., et al. (2018). A review of classification algorithms for EEG-based brain-computer interfaces: a 10 year update. J. Neural. Eng. 15:031005. doi: 10.1088/1741-2552/aab2f2

Lotte, F., Congedo, M., Lécuyer, A., Lamarche, F., and Arnaldi, B. (2007). A review of classification algorithms for EEG-based brain-computer interfaces. J. Neural. Eng. 4, R1–R13.

Mac Donald, C. L., Dikranian, K., Bayly, P., Holtzman, D., and Brody, D. (2007). Diffusion tensor imaging reliably detects experimental traumatic axonal injury and indicates approximate time of injury. J. Neurosci. 27, 11869–11876. doi: 10.1523/jneurosci.3647-07.2007

Manganotti, P., Acler, M., Masiero, S., and Del Felice, A. (2015). TMS-evoked N100 responses as a prognostic factor in acute stroke. Funct. Neurol. 30, 125–130.

Mori, S., Oishi, K., Jiang, H., Jiang, L., Li, X., Akhter, K., et al. (2008). Stereotaxic white matter atlas based on diffusion tensor imaging in an ICBM template. Neuroimage 40, 570–582. doi: 10.1016/j.neuroimage.2007.12.035

Moura, L. M., Luccas, R., de Paiva, J. P. Q., Amaro, E., Leemans, A., Leite, C. D. C., et al. (2019). Diffusion tensor imaging biomarkers to predict motor outcomes in stroke: a narrative review. Front. Neurol. 10:445. doi: 10.3389/fneur.2019.00445

Nakayama, N., Okumura, A., Shinoda, J., Yasokawa, Y. T., Miwa, K., Yoshimura, S. I., et al. (2006). Evidence for white matter disruption in traumatic brain injury without macroscopic lesions. J. Neurol. Neurosurg. Psychiatry 77, 850–855. doi: 10.1136/jnnp.2005.077875

Niogi, S. N., and Mukherjee, P. (2010). Diffusion tensor imaging of mild traumatic brain injury. J. Head Trauma Rehabil. 25, 241–255.

Opie, G. M., Foo, N., Killington, M., Ridding, M. C., and Semmler, J. G. (2019). Transcranial magnetic stimulation-electroencephalography measures of cortical neuroplasticity are altered after mild traumatic brain injury. J. Neurotrauma 36, 2774–2784. doi: 10.1089/neu.2018.6353

Pennisi, G., Bella, R., and Lanza, G. (2015). Motor cortex plasticity in subcortical ischemic vascular dementia: what can TMS say? Clin. Neurophysiol. 126, 851–852. doi: 10.1016/j.clinph.2014.09.001

Pennisi, M., Lanza, G., Cantone, M., Ricceri, R., Spampinato, C., Pennisi, G., et al. (2016). Correlation between motor cortex excitability changes and cognitive impairment in vascular depression: pathophysiological insights from a longitudinal TMS study. Neural. Plast. 2016:8154969.

Puig, J., Blasco, G., Daunis-I-Estadella, J., Thomalla, G., Castellanos, M., Figueras, J., et al. (2013). Decreased corticospinal tract fractional anisotropy predicts long-term motor outcome after stroke. Stroke 44, 2016–2018. doi: 10.1161/strokeaha.111.000382

Puig, J., Ellis, M. J., Kornelsen, J., Figley, T. D., Figley, C. R., Daunis-I-Estadella, P., et al. (2020). Magnetic resonance imaging biomarkers of brain connectivity in predicting outcome after mild traumatic brain injury: a systematic review. J. Neurotrauma 37, 1761–1776. doi: 10.1089/neu.2019.6623

Rogasch, N. C., and Fitzgerald, P. B. (2013). Assessing cortical network properties using TMS-EEG. Hum. Brain Mapp. 34, 1652–1669. doi: 10.1002/hbm.22016

Rossini, P. M., Burke, D., Chen, R., Cohen, L. G., Daskalakis, Z., Di Iorio, R., et al. (2015). Non-invasive electrical and magnetic stimulation of the brain, spinal cord, roots and peripheral nerves: basic principles and procedures for routine clinical and research application. an updated report from an I.F.C.N. Committee. Clin. Neurophysiol. 126, 1071–1107. doi: 10.1016/j.clinph.2015.02.001

Rossini, P. M., Miraglia, F., Alù, F., Cotelli, M., Ferreri, F., Iorio, R. D., et al. (2020). Neurophysiological hallmarks of neurodegenerative cognitive decline: the study of brain connectivity as a biomarker of early dementia. J. Pers. Med. 10:34. doi: 10.3390/jpm10020034

Salvadores, N., Sanhueza, M., Manque, P., and Court, F. A. (2017). Axonal degeneration during aging and its functional role in neurodegenerative disorders. Front. Neurosci. 11:451. doi: 10.3389/fnins.2017.00451

Shafi, M. M., Westover, M. B., Fox, M. D., and Pascual-Leone, A. (2012). Exploration and modulation of brain network interactions with noninvasive brain stimulation in combination with neuroimaging. Eur. J. Neurosci. 35, 805–825. doi: 10.1111/j.1460-9568.2012.08035.x

Tremblay, S., Rogasch, N. C., Premoli, I., Blumberger, D. M., Casarotto, S., Chen, R., et al. (2019). Clinical utility and prospective of TMS-EEG. Clin. Neurophysiol. 130, 802–844.

Wakana, S., Caprihan, A., Panzenboeck, M. M., Fallon, J. H., Perry, M., Gollub, R. L., et al. (2007). Reproducibility of quantitative tractography methods applied to cerebral white matter. Neuroimage 36, 630–644. doi: 10.1016/j.neuroimage.2007.02.049

Wakana, S., Jiang, H., Nagae-Poetscher, L. M., van Zijl, P. C., and Mori, S. (2004). Fiber tract-based atlas of human white matter anatomy. Radiology 230, 77–87. doi: 10.1148/radiol.2301021640

Wang, Y., Liu, G., Hong, D., Chen, F., Ji, X., and Cao, G. (2016). White matter injury in ischemic stroke. Prog. Neurobiol. 141, 45–60.

Xu, J., Rasmussen, I. A., Lagopoulos, J., and Håberg, A. (2007). Diffuse axonal injury in severe traumatic brain injury visualized using high-resolution diffusion tensor imaging. J. Neurotrauma 24, 753–765. doi: 10.1089/neu.2006.0208

Zifman, N., Levy-Lamdan, O., Suzin, G., Efrati, S., Tanne, D., Fogel, H., et al. (2019). Introducing a novel approach for evaluation and monitoring of brain health across life span using direct non-invasive brain network electrophysiology. Front. Aging Neurosci. 11:248. doi: 10.3389/fnagi.2019.00248