Connectome reorganization associated with temporal lobe pathology and its surgical resection

Sara Larivière, Bo-yong Park, Jessica Royer, Jordan DeKraker, Alexander Ngo, Ella Sahlas, Judy Chen, Raúl Rodríguez-Cruces, Yifei Weng, Birgit Frauscher, Ruoting Liu, Zhengge Wang, Golia Shafiei, Bratislav Mišić, Andrea Bernasconi, Neda Bernasconi, Michael D. Fox, Zhiqiang Zhang, Boris C. Bernhardt
doi: https://doi.org/10.1101/2023.11.13.23298482

Abstract

Network neuroscience offers a unique framework to understand the organizational principles of the human brain. Despite recent progress, our understanding of how the brain is modulated by focal lesions remains incomplete. Resection of the temporal lobe is the most effective treatment to control seizures in pharmaco-resistant temporal lobe epilepsy (TLE), making this syndrome a powerful model to study lesional effects on network organization in young and middle-aged adults. Here, we assessed the downstream consequences of a focal lesion and its surgical resection on the brain’s structural connectome, and explored how this reorganization relates to clinical variables at the individual patient level. We included adults with pharmaco-resistant TLE (n = 37) who underwent anterior temporal lobectomy between two imaging time points, as well as age-and sex-matched healthy controls who underwent comparable imaging (n = 31). Core to our analysis was the projection of high-dimensional structural connectome data—derived from diffusion MRI tractography from each subject—into lower-dimensional gradients. We then compared connectome gradients in patients relative to controls before surgery, tracked surgically-induced connectome reconfiguration from pre-to postoperative time points, and examined associations to patient-specific clinical and imaging phenotypes. Before surgery, TLE presented with marked connectome changes in bilateral temporo-parietal regions, reflecting an increased segregation of the ipsilateral anterior temporal lobe from the rest of the brain. Surgery-induced connectome reorganization was localized to this temporo-parietal subnetwork, but primarily involved postoperative integration of contralateral regions with the rest of the brain. Using a partial least-squares analysis, we uncovered a latent clinical-imaging signature underlying this pre-to postoperative connectome reorganization, showing that patients who displayed postoperative integration in bilateral fronto-occipital cortices also had greater preoperative ipsilateral hippocampal atrophy, lower seizure frequency, and secondarily generalized seizures. Our results bridge the effects of focal brain lesions and their surgical resections with large-scale network reorganization and inter-individual clinical variability, thus offering new avenues to examine the fundamental malleability of the human brain.

Introduction

Human brain organization is increasingly conceptualized and analyzed from a network perspective, allowing great strides in understanding both health and disease. In recent years, these approaches have been guided by advances in neuroimaging acquisition and complex data analytics. In particular, the advent of diffusion MRI tractography has enabled the approximation of structural connections in vivo, and the systematic characterization of brain connectivity—so called connectomes.1 In parallel, a series of multivariate analytics has been developed to capture organizational principles of large-scale brain connectivity. These range from graph theoretical assessments of network topology2–4 to modular decompositions that identify interacting systems.5–7 More recently, dimensionality reduction methods that derive compact and continuous embeddings from high dimensional connectomes have begun to reveal salient spatial axes of connectivity.8,9 In the current work, we adopted the latter approach to establish how structural connectivity embeddings can be modulated by focal lesions.

We specifically study patients with pharmaco-resistant temporal lobe epilepsy (TLE), who suffer from seizures originating from pathological alterations in the temporal lobe despite anti-seizure medical treatment.10 Compared to patients with adequate seizure control, pharmaco-resistant patients present with elevated risk for comorbidity and mortality,11 impaired quality of life and wellbeing,12 cognitive dysfunction and psychiatric difficulties,13–16 as well as progressive brain damage.17–19 In these patients, randomized controlled trials have shown that resection of the affected temporal lobe is currently the most effective treatment to control seizures, to restore psychosocial functioning, and to improve quality of life.20,21 By applying dimensional decomposition approaches to structural connectivity data acquired in healthy individuals as well as patients before and after surgery, we can investigate the impact of focal lesions on connectome-level organization.

Our novel framework, thereby, complements previous investigations of structural connectivity alterations in TLE patients relative to controls before surgery, which have reported microstructural and architectural damage seldom limited to the temporal lobe or to the hemisphere ipsilateral to the seizure focus.22–25 Similarly, prior work has assessed topological changes based on graph theoretical parameterization, and has accumulated evidence for a modular rearrangement in networks that extends beyond the temporal lobe.26–30 Such findings are compatible with the notion that “focal” pathologies generally affect large-scale brain connectivity, and that large-scale compromise contributes to patient presentation and clinically relevant outcomes.31–35 In line with these findings, studies assessing brain network changes following epilepsy surgery suggest that a focal resection does not only impact the primary lesioned site, but rather has effects that cascade along direct and indirect connections into mutually interconnected networks.36–38 As such, TLE represents a powerful human lesion model to assess how focal pathology and its resection affect whole-brain organization.

Here, we leveraged pre-as well as postoperative diffusion MRI data acquired in a cohort of pharmaco-resistant TLE patients. Our analytical paradigm derived whole brain connectomes from diffusion tractographic modelling, followed by non-linear dimensionality reduction to derive a continuous coordinate system in which confluent spatial trends of structural connectivity variation can be investigated.8,39 We compared our patients to a group of healthy controls, to investigate preoperative alterations and we tracked connectome reorganization from pre-to postoperative time points in patients. Exploiting inter-patient heterogeneity, we furthermore examined associations of pre-to postoperative connectivity trends with clinical and other imaging phenotypes at the individual patient level.

Methods

Participants

We studied 37 consecutive people suffering from pharmaco-resistant TLE, who (i) underwent anterior temporal lobectomy as a treatment of their seizures at Jinling Hospital between 2009 and 2018, (ii) had postoperative histological confirmation of hippocampal sclerosis, (iii) underwent a research-dedicated, high-resolution 3T MRI before and after surgery that included T1-weighted (T1w) and diffusion MRI scans, and (iv) had at least 1 year of postoperative follow-up. Patients were diagnosed according to the classification of the International League Against Epilepsy based on a comprehensive examination that includes clinical history, seizure semiology, continuous video-electroencephalographic (EEG) telemetry recordings, neuroimaging, and neuropsychology. No patient had encephalitis, malformations of cortical development (e.g., tumors, vascular malformations), or a history of traumatic brain injury. We also studied 31 age-and sex-matched healthy individuals who underwent identical 3T MRI at one time point. Demographic and clinical information on both cohorts is provided in Table S1. Our study was approved by the research ethics board of Jinling Hospital, Nanjing University School of Medicine, and written informed consent was obtained from all participants.

MRI acquisition and processing

All MRI data were obtained on a Siemens Trio 3T scanner (Siemens, Erlangen, Germany) and included: (i) high-resolution 3D T1w MRI using a magnetization-prepared rapid gradient-echo sequence (repetition time [TR] = 2300 ms, echo time [TE] = 2.98 ms, flip angle = 9°, voxel size = 0.5 × 0.5 × 1 mm3, field of view [FOV] = 256 × 256 mm2, 176 slices) and (ii) diffusion MRI using a spin echo-based echo planar imaging sequence with 4 b0 images (dMRI, TR = 6100 ms, TE = 93 ms, flip angle = 90°, FOV = 240 × 240 mm2, voxel size = 0.94 × 0.94 × 3 mm3, b-value = 1000 s/mm2, diffusion directions = 120).

Multimodal MRI preprocessing was performed using micapipe, an open-access image processing and data fusion pipeline (v0.1.4; https://micapipe.readthedocs.io).40 In brief, native T1w structural images were deobliqued, reoriented to standard neuroscience orientation, corrected for intensity nonuniformity, intensity normalized, skull-stripped, and submitted to FreeSurfer (v6.0)41 to extract models of the cortical surfaces. Subcortical structures were segmented using FSL FIRST.42 Native DWI images were denoised, underwent b0 intensity normalization, and were corrected for susceptibility distortion, head motion, and eddy currents. Left and right hemispheric data of individuals with TLE were sorted into ipsilateral/contralateral to the seizure focus. To minimize confounds related to normal inter-hemispheric asymmetries, prior to sorting we z-scored data in patients relative to controls.

Surgical cavity mapping

Patient-specific surgical cavities were automatically segmented by registering both pre-and postoperative T1w images to the MNI152 standard template through linear transformations and subtracting the postoperative scan from the preoperative scan. Segmented cavities were visually inspected and manually edited to ensure that the extent of the resections was correctly identified.

Cavities were then mapped to the surface template, and a consensus label was generated as defined by the union of all segmentations.

Structural connectome generation

Structural connectomes were generated with MRtrix3 from pre-processed DWI data.43 We performed anatomically-constrained tractography using tissue types (cortical and subcortical grey matter, white matter, cerebrospinal fluid) segmented from each participant’s pre-processed T1w images registered to native DWI space. We estimated multi-tissue response functions, performed constrained spherical-deconvolution and intensity normalization, generated a tractogram with 40M streamlines (maximum tract length = 250; fractional anisotropy cut-off = 0.06), and applied spherical deconvolution informed filtering of tractograms (SIFT2) to reconstruct whole-brain streamlines weighted by cross-sectional multipliers. Reconstructed streamlines were mapped onto 400 similar-sized cortical parcels constrained within the boundaries of the Desikan-Killiany atlas44,45 and 14 regions in the subcortex and the hippocampus to produce subject-specific structural connectivity matrices. Group-average normative structural connectomes were defined using a distance-dependent thresholding procedure, which preserved the edge length distribution in individual patients,46 and were log transformed to reduce connectivity strength variance. As such, structural connectivity was defined by the number of streamlines between two regions (i.e., fiber density). Structural connectomes were thresholded to retain only the top 25% of connections.

Connectome gradient estimation

Cortex-and subcortex-wide structural connectome gradients were generated using BrainSpace (v0.1.10; https://brainspace.readthedocs.io).45 First, a group-level structural connectome (combining cortical and subcortical data) was obtained using an iterative, leave-one-out procedure (Figure 1). For both the group-level and the left-out subject-specific connectome, an affinity matrix was constructed with a normalized angle kernel, and eigenvectors were estimated via diffusion map embedding, an unsupervised nonlinear dimensionality reduction technique that projects connectome features into low-dimensional gradients.9,47 This technique is only controlled by only few parameters, computationally efficient, relatively robust to noise compared to other nonlinear techniques,48 and has been extensively used in the previous literature.8 Algorithm parameters were identical to those from prior applications, specifically α = 0.05 and diffusion time t = 0 to retain the global relations between data points in the embedded space and to emphasize direct connections.9 Interhemispheric connections were not included in the gradient computation; left and right gradients were generated separately and aligned with Procrustes.45 The eigenvectors estimated from this technique provide a connectivity coordinate system—a diffusion map—where Euclidean distances in the manifold correspond to diffusion times between the nodes of the network.47 In this gradient space, regions with similar connectivity profiles are closely located, whereas regions with different connectivity profiles are located farther apart. To analyze gradient changes within and across participants, each left-out subject-specific gradient was aligned to the group template gradient (estimated using data from all other participants) via Procrustes alignment,45 normalized, and sorted into ipsilateral/contralateral to the seizure focus. The scree plot describing the eigenvalue decay function revealed that four principal eigenvectors (i.e., gradients) corresponded to the cut-off point, after which eigenvalues leveled off at small values and only explained marginal information.

Figure 1.
Figure 1. Construction of structural connectome gradients.

Subject-specific structural connectome gradients were generated from diffusion MRI data and iteratively aligned to a group template (using data from all other participants). The first four gradients were retained.

Quantifying pre-and postoperative changes in structural gradients

We assessed cross-sectional connectome changes between individuals with TLE (preoperatively) and controls using univariate fixed effects models performed on each eigenvector independently.49 Eigenvector-specific score histograms were then generated across the entire cortex, and split into ipsilateral/contralateral as well as patients/controls to examine changes in the distribution of significant scores and to assess inter-hemispheric asymmetry. To signify an overall load of alterations, we also compared the aggregate of the first four eigenvectors in patients relative to controls using a multivariate fixed effect linear model. By statistically combining multidimensional gradients, the latter approach leverages their covariance to obtain a substantial gain in sensitivity, and thus unveil subthreshold properties not readily identified in a single eigenvector. Both uni-and multivariate linear models tested for differences between groups while controlling for effects of age and sex using BrainStat (v0.4.2; https://brainstat.readthedocs.io).49

To test for changes over time (i.e., pre-vs. postsurgery), we fitted univariate and multivariate linear mixed effects models, a flexible framework for the analysis of repeated measures.49,50 Our models controlled for effects of age and sex, and included a subject-specific random intercept to improve model fit in longitudinal designs. With this approach, we were able to test for pre-to postoperative within-subject individual and multivariate eigenvector changes.

Findings from all fixed and mixed effects models were corrected for multiple comparisons using the false discovery rate (FDR) procedure.51

4D gradient deformations

To simplify multidimensional changes into a single deformation scalar feature, we computed the Euclidean distance between the group template center and each subject-specific data point in the 4D gradient space (Figure 2B). We assessed patterns of eigenvector contraction (i.e., connectome integration) and expansion (i.e., connectome differentiation) relative to the centre of the gradient space. The template center was defined as the centroid of the first four eigenvectors. Shifts in structural connectivity patterns lead to a displacement of eigenvector scores in the 4D gradient space, which in turn affects their proximity to the center. Gradient contraction/expansion thus quantifies global brain reorganization in connectivity space.52,53 Gradient deformation was quantified only in regions of interest, as determined by areas that showed significant cross-sectional and longitudinal changes, respectively. The rationale was to determine whether preoperative and surgery-induced 4D gradient deformations were observed in areas that showed significant eigenvector changes.

Figure 2.
Figure 2. Preoperative structural connectome gradient changes in TLE.

(a) Univariate (left) and multivariate (right) fixed effect models compared gradient values in individuals with TLE (before surgery) and healthy controls. Significant gradient changes in patients relative to controls were observed in bilateral orbitofrontal and temporoparietal cortices, ipsilateral hippocampus, and contralateral amygdala (all pFDR < 0.05). (b) Multivariate gradient changes were simplified into a scalar feature of connectome deformations by computing the Euclidean distance (di,j) between the group template center (C) and each subject-specific data point in the 4D gradient space. Restricting this analysis to areas of significant gradient changes identified in (a) revealed increased segregation of temporo-parietal regions in TLE relative to controls.

Spatial associations between preoperative and surgery-induced reorganization

To assess whether preoperative alterations constrained surgery-induced changes, we performed spatial associations between cross-sectional (relative to controls) and longitudinal (pre-to postoperative) connectome findings. Statistical significance of spatial correlations was assessed using variogram-matching models that take into account the spatial dependencies in the data.45,54,55 This method generates surrogate brain maps with matched spatial autocorrelation to that of a target brain, and has previously been applied to cortical and subcortical structures. As recommended,54 surrogate maps were generated using surface-based geodesic distance between cortical regions and three-dimensional Euclidean distance between subcortical and cortical/subcortical regions. Variogram-matching null distributions were then generated from randomly shuffling surrogate maps while preserving the distance-dependent correlation between elements of the brain map. The empirical correlation coefficients were compared against the null distribution determined by the ensemble of spatially permuted correlation coefficients to estimate a p-value (termed pvar).

Clinical associations

We used partial least squares (PLS) analysis to unravel within-patient relationships between surgery-induced 4D gradient deformations and standard neuroimaging (i.e., hippocampal atrophy56) and clinical parameters (Figure 4A). PLS analysis is a multivariate associative technique that maximizes the covariance between two sets of variables.57 Briefly, these two variable sets are correlated with each other across patients. The resulting correlation matrix is submitted to singular value decomposition, which identifies linear combinations of the original variables to generate new latent variables that have maximum covariance.58 We evaluated the significance of our model using nonparametric methods: (i) permutation tests were used to assess statistical significance of each latent variable (10,000 permutations; termed pperm), (ii) bootstrap ratios were used to assess the reliability of singular vector weights (analogous to z-scores, such that 95% confidence interval corresponds to a bootstrap ratio of ±1.96 and ±2.58, respectively), and (iii) cross-validation was used to assess out-of-sample correlations between projected scores (using 100 randomized train-test splits of the original data, where 75% of the data was treated as a training set and 25% of the data was treated as an out-of-sample test set). Mathematical details of the analysis and inferential methods are described elsewhere.57,59

Results

Connectome gradient analysis at baseline

For every participant, we generated cortex-and subcortex-wide structural connectome gradients using non-linear dimensionality reduction.45 The first four gradients explained approximately 75.28 ± 1.89% (baseline/preoperatively) and 78.85 ± 1.92% (postoperatively) of variance across participants, with each gradient representing a different axis of connectivity variation (Figure 1). Across both time points and in both groups, gradients generally depicted smooth differentiation in structural connectivity profiles of frontal vs temporo-occipital regions (Gradient 1), orbitofrontal vs sensorimotor regions (Gradient 2), medial occipital vs anterior temporal lobe regions (Gradient 3), and sensory vs transmodal regions (Gradient 4).

Compared to healthy controls, presurgical TLE patients showed significant gradient changes in bilateral temporo-parietal (pFDR < 0.001) and orbitofrontal (pFDR < 0.005) cortices, together with the ipsilateral hippocampus (pFDR < 0.01) and the contralateral amygdala (pFDR < 0.005; Figure 2A). Considering each gradient separately, global histogram analysis revealed significant asymmetry between ipsilateral and contralateral gradient alterations in TLE. Most notably, ipsilateral mesiotemporal lobe regions were shifted toward the most extreme, segregated positions along Gradients 2, 3, and 4, whereas their contralateral homologues were shifted to a more central, integrated position (Figure S1). This ipsilateral mesiotemporal lobe segregation was further confirmed by quantifying their positional change in the 4D gradient space relative to controls. Constraining this analysis to areas of significant preoperative gradient changes, individuals with TLE showed highest deformation values (i.e., depicting increased expansion/segregation) in ipsilateral anterior temporal lobe (puncorr < 0.005) and parahippocampal gyrus (puncorr < 0.01; Figure 2B).

Tracking surgery-induced deformations

Directly comparing structural connectome changes before and after anterior temporal lobectomy (Figure 3A), patients showed significant multivariate gradient changes in posterior temporal regions near the resection site (pFDR < 0.05) and contralateral temporo-parietal cortices (pFDR < 0.05; Figure 3B). As evidenced in both univariate and 4D gradient deformation analyses, contralateral temporo-parietal regions were shifted towards more integrated positions in the 4D gradient space after surgery, particularly along the axis of Gradient 1 (pFDR < 0.05; Figure S2).

Figure 3.
Figure 3. Surgery-induced structural connectome gradient changes in TLE.

(a) Surgical cavities were automatically segmented from pre-and postoperative T1w MRIs. Cavity masks were excluded a posteriori in preoperative analysis and a priori in postoperative analyses. (b) Multivariate mixed effects models revealed a pattern of surgery-induced gradient changes encompassing ipsilateral posterior temporal and contralateral temporo-parietal regions (all pFDR < 0.05). (c) Analysis of 4D gradient deformations constrained to areas of multivariate gradient change revealed segregation of ipsilateral temporal regions but increased integration of contralateral temporo-parietal regions after surgery.

As all regions showing pre-to postoperative gradient changes fell within the boundaries of significant preoperative changes, we next assessed the extent to which preoperative connectome reorganization relates to surgery-induced changes. The degree of preoperative multivariate gradient changes was correlated with pre-to postoperative multivariate gradient changes across all regions (r = 0.74, pvar < 0.0001), as well as when considering only areas of significant preoperative (r = 0.72, pvar < 0.0001) and pre-to postoperative (r = 0.49, pvar = 0.05) gradient changes (Figure S3). In other words, regions with profound connectome alterations before surgery (relative to controls) corresponded to regions that also underwent the largest pre-to postoperative reorganization. We also found a significant, albeit weaker, negative spatial associations between whole-brain preoperative and surgery-induced 4D gradient deformations (r = –0.15, pvar < 0.005). In this case, a negative correlation indicates that regions that were initially segregated from the rest of the brain before surgery subsequently shifted towards a more integrated position in the connectome space after surgery.

Associations with clinical variables

Multivariate PLS analysis identified one statistically significant latent variable (LV-1 accounting for 44.65% of covariance, pperm < 0.05) that maximized the covariance between whole-brain 4D surgery-induced deformation patterns and standard neuroimaging parameters (e.g., hippocampal atrophy) and clinical measures (Figure 4A). Loadings (i.e., correlations of individual clinical measures with LV-1) revealed ipsilateral hippocampal atrophy in CA1-3 and dentate gyrus subfields (all r < –0.28), lower seizure frequency (r < –0.22), and presence of secondarily generalized seizures (r > 0.35) as the strongest contributors of LV-1. Patients presenting with these clinical features showed increased postoperative connectome integration in bilateral fronto-occipital cortices, ipsilateral postcentral gyrus, and contralateral mesiotemporal cortex (all bootstrap ratio < 95% confidence interval). Clinical measures and 4D surgery-induced gradient deformations were correlated at the individual patient level (r = 0.52), while out-of-sample correlations averaged r = 0.20 (and exceeded mean permuted out-of-sample correlations of r = – 0.028; Figure 4B).

Figure 4.
Figure 4. A clinical-imaging signature of surgery-induced connectome reorganization.

(a) Partial least-squares (PLS) analysis was used to relate patient-specific pre-to postoperative 4D gradient deformation with clinical variables. These two sets of variables are correlated across patients and subjected to singular value decomposition, yielding multiple latent variables. Bootstrap resampling evaluated the contribution of individual variables to the latent variable, whereas permutation tests and cross-validation assessed the pairing of the 4D surgery-induced deformation patterns and clinical variables. (b) The first latent variable (LV-1) accounted for 44.65% of the covariance between the imaging and clinical data (p < 0.05; top left). The contribution of individual clinical measures is shown using correlations between patient-specific clinical scores and scores on the multivariate pattern (loadings; top right). Error bars indicate bootstrap-estimated SEs. The contribution of individual regions is shown using bootstrap ratios (interpretable as z-scores; bottom left). Patients who displayed increased postoperative connectome integration had greater ipsilateral hippocampal atrophy in CA1-3 and DG subfields, lower seizure frequency, and secondarily generalized seizures. This association was confirmed by projecting individual patient data onto the weighted patterns (r = 0.52; bottom middle). Correlations between 4D surgery-induced deformations and clinical scores in cross-validated held-out data (mean ± SD r = 0.20 ± 0.29) were higher than those from the permuted null model (mean ± SD r = –0.028 ± 0.33; bottom right).

Discussion

Contemporary systems neuroscience is increasingly exploring advanced dimensionality reduction techniques to delineate the main organizational axes of the human connectome, and to derive potential clinical predictors.48 Here, we used connectome gradient mapping—an established dimensionality reduction method—to investigate local-global brain interactions in individuals with pharmaco-resistant TLE. We capitalized on this condition as a human lesion model to deepen our understanding of network-level consequences of focal lesions and their resection. We found that, before surgery, individuals with TLE showed an increased segregation of the ipsilateral temporal lobe (i.e., the disease epicenter) from the rest of the brain. By tracking pre-to postoperative connectome alterations, we showed that resection of the ipsilateral temporal lobe led to brain reorganization involving areas near the resection site, but also marked re-integration of contralateral temporo-parietal regions with the rest of the brain. We further provided proof-of-concept evidence of an association between inter-individual variations in connectome reorganization following surgery and clinically salient features. Taken together, our results demonstrate the potential of gradient mapping to bridge different scales of human brain organization and allow the study of the brain’s structural plasticity in response to focal lesions and their eventual resection in young adults.

Core to our analytical framework was the estimation of connectome gradients that index the spatial arrangement of macroscale organization along continuous cortico-subcortical axes. This layout serves as a meaningful coordinate system to contextualize different MRI-derived measures with other markers of neural organization. For instance, such gradient approaches have been used to investigate the associations between macroscale brain organization and transcriptomic,60 synaptic,61 molecular,62 microstructural,63,64 and structure-function9,65 factors. Gradient mapping studies aiming to better understand brain changes observed in psychiatric and neurological conditions—including epilepsy—are also starting to emerge. Indeed, recent stepwise functional connectivity studies and applications of dimensional decompositions have described an imbalance of functional integration and segregation between sensory and transmodal systems in TLE.66,67 Our study applied gradient mapping for the first time to structural connectivity data in epilepsy. Using diffusion MRI, we demonstrated shifts in macroscale connectome organization implicating strategic regions such as the disease epicenter itself, as well as distant, yet direct, temporo-parietal connections, reflecting an increased segregation of these regions from the rest of the brain. Our findings further extend previous diffusion MRI studies that have focused on specific characteristics of white matter tracts either locally or globally,25,68,69 and instead simultaneously evaluated structural connectome changes both at the lesion site and in distant regions. Our approach provides evidence that a condition with a focal lesion to the mesiotemporal lobe presents with changes that extend well beyond the site of primary pathology, affecting directly and indirectly structurally connected macroscale networks. Debate persists over the precise mechanisms that underlie changes in these structural networks, but such changes may reflect excitotoxic effects from seizure activity or deafferentation of mesiotemporal connections.70 In this context, a possible explanation for the observed segregation of the ipsilateral anterior temporal cortex may relate to hyperconnectivity within the epileptogenic mesiotemporal network, together with hypoconnectivity of the epileptogenic mesiotemporal network to the rest of the brain, or a combination of both. Under this account, topological “isolation” of the ipsilateral temporal lobe may represent a core connectional change that may incur locally sustained hyperexcitability and increased susceptibility to seizure activity.71 On the other hand, hypoconnectivity between ipsilateral temporal lobe regions and all other regions may be considered as a compensatory mechanism through which epileptogenic activity is contained, abrogated, and prevented from spreading.72,73

A common surgical procedure to alleviate seizures in patients with medically-intractable TLE is anterior temporal lobectomy; i.e., the resection of 4-4.5 cm of the temporal lobe (or up to 6.5 cm in case of non-dominant TLE).21 Here, we carried out longitudinal neuroimaging to track connectome reorganization before and after anterior temporal lobectomy. In our study, we showed that surgery-induced changes occur primarily in regions that were also affected preoperatively. The affected areas were bilateral, spanning regions near the resection site and contralateral temporo-parietal cortices. Albeit marginal, we showed evidence of postoperative connectome segregation of ipsilateral posterior temporal lobe regions, possibly pointing to ongoing Wallerian degeneration in nerve bundles disconnected after surgery. This process generally occurs in two phases, with acute dying back of axons within the first week postsurgery,74 followed by chronic myelin degradation lasting several months.75 With respect to the latter, there is ample evidence that ipsilateral temporal lobe fibers undergo irreversible Wallerian degeneration, with diffusion changes being apparent throughout the first year postsurgery.75–77 In line with these findings, a possible explanation for our observation of contralateral postoperative changes is the degeneration of commissural fibers connecting the two hippocampi and their afferent pathways.23 In contrast to preoperative alterations, changes in contralateral temporo-parietal cortices, which are expected to incur minimal downstream microstructural damage with surgery, were predominantly characterized by increased postoperative connectome integration. This shift in their connectivity profiles may be related to structural plasticity of the white matter in an active attempt to compensate for the resected grey matter and ongoing axonal/myelin degradation in homologous regions. Although such a relationship remains speculative, previous data have demonstrated that diffusion changes in contralateral fiber tracts fail to normalize following surgery.76 The absence of longitudinal data in controls in the current study may have hindered our ability to assess whether connectivity reorganization stabilizes, heals, or reorganizes over time relative to controls. Test-retest reliability studies of both diffusion MRI-derived connectomes78 and gradient-based frameworks48 have nevertheless shown high reliability in healthy adults.

Mining inter-individual variability in MRI data has been a cornerstone concept in neuroscience research.79,80 As a proof-of-concept, here we translated knowledge about connectome reorganization into individualized and clinically-meaningful taxonomies using latent clinical-imaging dimensions that link pre-to postsurgical shifts in connectome organization with heterogenous clinical measures. Our findings captured unique associations between connectome re-integration following surgery, primarily in bilateral fronto-occipital cortices, ipsilateral postcentral gyrus, and contralateral mesiotemporal regions, and greater ipsilateral hippocampal atrophy, lower seizure frequency, and presence of secondarily generalized seizures. As postsurgical seizure outcomes did not emerge as a significant contributor in our model, it is possible that a combination of both preoperative and postoperative features may better reflect seizure outcome, as suggested by a range of machine learning studies reporting associations between network measures and outcomes.29,81–86 Alternatively, outcome predictors may be identified from lower order latent variables, ideally with associated patterns of connectivity shifts centered on contralateral temporo-parietal and default-mode regions.73 Future studies assessing predictive features will ideally include multi-site data to assess predictive generalizability, enroll patients with longer-term seizure outcomes, and include both seizure and cognitive outcome measures to comprehensively assess prognostication of overall functioning and wellbeing. Nevertheless, the present results complement recent efforts to identify latent imaging-derived disease factors that not only variably co-express across TLE patients, but also have potential for individualized predictions.87

Taken together, these findings offer a new framework that can enrich our understanding of the effects of focal lesions and their eventual surgical resection on the macroscale organization of the human brain connectome. We dichotomized macroscale connectivity shifts before and after surgery in TLE, with the ipsilateral anterior temporal lobe—the target area for resection—being segregated from the rest of the brain prior to surgery. Secondary to surgery, contralateral temporo-parietal regions revealed connectivity shifts in a diametrically opposite direction, likely reflecting a compensatory response to restore or maintain structural connectivity. This data-driven connectome reorganization further accounted for a broad spectrum of clinical factors, potentially opening up new avenues for future neuroimaging studies to parse inter-individual differences in epilepsy and, consequently, may be translatable to individualized patient care. More broadly, a human lesion model, like the one described herein, offers a rare window into the building blocks of brain organization in both healthy and diseased states, and can ultimately allow causality to be inferred based on the effects of focal brain lesions in TLE.

Data Availability

All data produced in the present study are available upon reasonable request to the authors

Author contributions

S.L., Z.Z., and B.C.B contributed to the conception and design of the study; Y.W., Z.W, Z.Z, and R.L. contributed to the acquisition of data; S.L. and B.C.B contributed to data analysis and manuscript writing; revised and approved by all co-authors.

Potential conflicts of interest

None to disclose.

Data availability

Codes to conduct multimodal preprocessing (micapipe40), hippocampal segmentation (HippUnfold56), and statistical analysis (BrainStat49) are openly available.

Acknowledgments

S.L. acknowledges funding from Fonds de la Recherche du Québec–Santé (FRQ-S), the Canadian Institutes of Health Research (CIHR), and the Richard and Ann Sievers Neuroscience Award. B.Y.P. was supported by the National Research Foundation of Korea (NRF-2021R1F1A1052303; NRF-2022R1A5A7033499), Institute for Information and Communications Technology Planning and Evaluation (IITP) funded by the Korea Government (MSIT; No. 2022-0-00448, Deep Total Recall: Continual Learning for Human-Like Recall of Artificial Neural Networks; No. RS-2022-00155915, Artificial Intelligence Convergence Innovation Human Resources Development [Inha University]; No. 2021-0-02068, Artificial Intelligence Innovation Hub), and Institute for Basic Science (IBS-R015-D1). J.R. received support from the Canadian Open Neuroscience Platform (CONP) and CIHR. B.F. receives funding from FRQ-S (Chercheur-Boursier clinician Senior). A.N. acknowledges funding from FRQ-S. E.S. acknowledges funding from FRQ-S. G.S. acknowledges support from CIHR. B.M. acknowledges support from NSERC (Discovery Grant RGPIN 017-04265) and from the Canada Research Chairs Program (CRC). A.B. and N.B. were supported by FRQ-S and CIHR (MOP-57840, MOP-123520). Z.Z. was supported by the National Science Foundation of China (NSFC: 81422022; 863 project: 2014BAI04B05 and 2015AA020505) and the China Postdoctoral Science Foundation (2016M603064). B.C.B. acknowledges support from CIHR (FDN-154298, PJT-174995), SickKids Foundation (NI17-039), NSERC (Discovery-1304413), Azrieli Center for Autism Research of the Montreal Neurological Institute (ACAR), BrainCanada, FRQ-S, the Helmholtz International BigBrain Analytics and Learning Laboratory (Hiball), and the Canada Research Chairs program.

References

  1. 1.↵
    Sporns O, Tononi G, Kötter R. The Human Connectome: A Structural Description of the Human Brain. PLOS Computational Biology. 2005;1(4):e42. doi:10.1371/journal.pcbi.0010042
    OpenUrlCrossRefPubMed
  2. 2.↵
    Bassett DS, Bullmore ET. Small-World Brain Networks Revisited. Neuroscientist. 2017;23(5):499–516. doi:10.1177/1073858416667720
    OpenUrlCrossRefPubMed
  3. 3.↵
    Bullmore E, Sporns O. Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci. 2009;10(3):186–198. doi:10.1038/nrn2575
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    Bullmore ET, Bassett DS. Brain Graphs: Graphical Models of the Human Brain Connectome. Annual Review of Clinical Psychology. 2011;7(1):113–140. doi:10.1146/annurev-clinpsy-040510-143934
    OpenUrlCrossRefPubMed
  5. 5.↵
    Meunier D, Lambiotte R, Bullmore E. Modular and Hierarchically Modular Organization of Brain Networks. Frontiers in Neuroscience. 2010;4. Accessed October 10, 2023. https://www.frontiersin.org/articles/10.3389/fnins.2010.00200
  6. 6.
    Yeo BTT, Krienen FM, Sepulcre J, et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. Journal of Neurophysiology. 2011;106(3):1125–1165. doi:10.1152/jn.00338.2011
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    Uddin LQ, Yeo BTT, Spreng RN. Towards a Universal Taxonomy of Macro-scale Functional Human Brain Networks. Brain Topogr. 2019;32(6):926–942. doi:10.1007/s10548-019-00744-6
    OpenUrlCrossRefPubMed
  8. 8.↵
    Bernhardt BC, Smallwood J, Keilholz S, Margulies DS. Gradients in brain organization. Neuroimage. 2022;251:118987. doi:10.1016/j.neuroimage.2022.118987
    OpenUrlCrossRef
  9. 9.↵
    Margulies DS, Ghosh SS, Goulas A, et al. Situating the default-mode network along a principal gradient of macroscale cortical organization. Proc Natl Acad Sci U S A. 2016;113(44):12574–12579. doi:10.1073/pnas.1608282113
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    Chen Z, Brodie MJ, Liew D, Kwan P. Treatment Outcomes in Patients With Newly Diagnosed Epilepsy Treated With Established and New Antiepileptic Drugs: A 30-Year Longitudinal Cohort Study. JAMA Neurol. 2018;75(3):279–286. doi:10.1001/jamaneurol.2017.3949
    OpenUrlCrossRefPubMed
  11. 11.↵
    GBD 2016 Neurology Collaborators. Global, regional, and national burden of neurological disorders, 1990-2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet Neurol. 2019;18(5):459-480. doi:10.1016/S1474-4422(18)30499-X
    OpenUrlCrossRefPubMed
  12. 12.↵
    Langfitt JT, Westerveld M, Hamberger MJ, et al. Worsening of quality of life after epilepsy surgery: effect of seizures and memory decline. Neurology. 2007;68(23):1988–1994. doi:10.1212/01.wnl.0000264000.11511.30
    OpenUrlCrossRef
  13. 13.↵
    Helmstaedter C, Sonntag-Dillender M, Hoppe C, Elger CE. Depressed mood and memory impairment in temporal lobe epilepsy as a function of focus lateralization and localization. Epilepsy & Behavior. 2004;5(5):696–701. doi:10.1016/j.yebeh.2004.06.008
    OpenUrlCrossRefPubMedWeb of Science
  14. 14.
    Helmstaedter C, Kurthen M, Lux S, Reuber M, Elger CE. Chronic epilepsy and cognition: A longitudinal study in temporal lobe epilepsy. Annals of Neurology. 2003;54(4):425–432. doi:10.1002/ana.10692
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.
    Helmstaedter C, Kockelmann E. Cognitive Outcomes in Patients with Chronic Temporal Lobe Epilepsy. Epilepsia. 2006;47(s2):96–98. doi:10.1111/j.1528-1167.2006.00702.x
    OpenUrlCrossRefPubMed
  16. 16.↵
    Reyes A, Kaestner E, Ferguson L, et al. Cognitive phenotypes in temporal lobe epilepsy utilizing data-and clinically driven approaches: Moving toward a new taxonomy. Epilepsia. 2020;61(6):1211–1220. doi:10.1111/epi.16528
    OpenUrlCrossRef
  17. 17.↵
    Bernhardt BC, Worsley KJ, Kim H, Evans AC, Bernasconi A, Bernasconi N. Longitudinal and cross-sectional analysis of atrophy in pharmacoresistant temporal lobe epilepsy. Neurology. 2009;72(20):1747–1754. doi:10.1212/01.wnl.0000345969.57574.f5
    OpenUrlAbstract/FREE Full Text
  18. 18.
    Caciagli L, Bernasconi A, Wiebe S, Koepp MJ, Bernasconi N, Bernhardt BC. A meta-analysis on progressive atrophy in intractable temporal lobe epilepsy: Time is brain? Neurology. 2017;89(5):506–516. doi:10.1212/WNL.0000000000004176
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    Galovic M, de Tisi J, McEvoy AW, et al. Resective surgery prevents progressive cortical thinning in temporal lobe epilepsy. Brain. 2020;143(11):3262–3272. doi:10.1093/brain/awaa284
    OpenUrlCrossRef
  20. 20.↵
    Engel J, McDermott MP, Wiebe S, et al. Early surgical therapy for drug-resistant temporal lobe epilepsy: a randomized trial. JAMA. 2012;307(9):922–930. doi:10.1001/jama.2012.220
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    Wiebe S, Blume WT, Girvin JP, Eliasziw M, Effectiveness and Efficiency of Surgery for Temporal Lobe Epilepsy Study Group. A randomized, controlled trial of surgery for temporal-lobe epilepsy. N Engl J Med. 2001;345(5):311–318. doi:10.1056/NEJM200108023450501
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    Concha L, Kim H, Bernasconi A, Bernhardt BC, Bernasconi N. Spatial patterns of water diffusion along white matter tracts in temporal lobe epilepsy. Neurology. 2012;79(5):455–462. doi:10.1212/WNL.0b013e31826170b6
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    Concha L, Beaulieu C, Gross DW. Bilateral limbic diffusion abnormalities in unilateral temporal lobe epilepsy. Ann Neurol. 2005;57(2):188–196. doi:10.1002/ana.20334
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.
    Liu M, Bernhardt BC, Hong SJ, Caldairou B, Bernasconi A, Bernasconi N. The superficial white matter in temporal lobe epilepsy: a key link between structural and functional network disruptions. Brain. 2016;139(Pt 9):2431–2440. doi:10.1093/brain/aww167
    OpenUrlCrossRefPubMed
  25. 25.↵
    Yogarajah M, Powell HWR, Parker GJM, et al. Tractography of the parahippocampal gyrus and material specific memory impairment in unilateral temporal lobe epilepsy. Neuroimage. 2008;40(4):1755–1764. doi:10.1016/j.neuroimage.2007.12.046
    OpenUrlCrossRefPubMed
  26. 26.↵
    Bernhardt BC, Fadaie F, Liu M, et al. Temporal lobe epilepsy: Hippocampal pathology modulates connectome topology and controllability. Neurology. 2019;92(19):e2209–e2220. doi:10.1212/WNL.0000000000007447
    OpenUrlAbstract/FREE Full Text
  27. 27.
    Bonilha L, Nesland T, Martz GU, et al. Medial temporal lobe epilepsy is associated with neuronal fibre loss and paradoxical increase in structural connectivity of limbic structures. J Neurol Neurosurg Psychiatry. 2012;83(9):903–909. doi:10.1136/jnnp-2012-302476
    OpenUrlAbstract/FREE Full Text
  28. 28.
    DeSalvo MN, Douw L, Tanaka N, Reinsberger C, Stufflebeam SM. Altered structural connectome in temporal lobe epilepsy. Radiology. 2014;270(3):842–848. doi:10.1148/radiol.13131044
    OpenUrlCrossRefPubMed
  29. 29.↵
    Gleichgerrcht E, Keller SS, Drane DL, et al. Temporal Lobe Epilepsy Surgical Outcomes Can Be Inferred Based on Structural Connectome Hubs: A Machine Learning Study. Annals of Neurology. 2020;88(5):970–983. doi:10.1002/ana.25888
    OpenUrlCrossRefPubMed
  30. 30.↵
    He Y, Evans A. Graph theoretical modeling of brain connectivity. Current Opinion in Neurology. 2010;23(4):341. doi:10.1097/WCO.0b013e32833aa567
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    Fox MD. Mapping Symptoms to Brain Networks with the Human Connectome. N Engl J Med. 2018;379(23):2237–2245. doi:10.1056/NEJMra1706158
    OpenUrlCrossRefPubMed
  32. 32.
    Larivière S, Vos de Wael R, Paquola C, et al. Microstructure-Informed Connectomics: Enriching Large-Scale Descriptions of Healthy and Diseased Brains. Brain Connect. 2019;9(2):113–127. doi:10.1089/brain.2018.0587
    OpenUrlCrossRefPubMed
  33. 33.
    Schaper FLWVJ, Nordberg J, Cohen AL, et al. Mapping Lesion-Related Epilepsy to a Human Brain Network. JAMA Neurol. 2023;80(9):891. doi:10.1001/jamaneurol.2023.1988
    OpenUrlCrossRef
  34. 34.
    Tavakol S, Royer J, Lowe AJ, et al. Neuroimaging and connectomics of drug-resistant epilepsy at multiple scales: From focal lesions to macroscale networks. Epilepsia. 2019;60(4):593–604. doi:10.1111/epi.14688
    OpenUrlCrossRefPubMed
  35. 35.↵
    Thiebaut de Schotten M, Foulon C, Nachev P. Brain disconnections link structural connectivity with function and behaviour. Nat Commun. 2020;11(1):5094. doi:10.1038/s41467-020-18920-9
    OpenUrlCrossRef
  36. 36.↵
    Ji GJ, Zhang Z, Xu Q, et al. Connectome Reorganization Associated With Surgical Outcome in Temporal Lobe Epilepsy. Medicine (Baltimore). 2015;94(40):e1737. doi:10.1097/MD.0000000000001737
    OpenUrlCrossRef
  37. 37.
    Morgan VL, Rogers BP, González HFJ, Goodale SE, Englot DJ. Characterization of postsurgical functional connectivity changes in temporal lobe epilepsy. J Neurosurg. Published online June 14, 2019:1–11. doi:10.3171/2019.3.JNS19350
    OpenUrlCrossRef
  38. 38.↵
    Neal EG, Maciver S, Schoenberg MR, Vale FL. Surgical disconnection of epilepsy network correlates with improved outcomes. Seizure. 2020;76:56–63. doi:10.1016/j.seizure.2020.01.018
    OpenUrlCrossRef
  39. 39.↵
    Park BY, Paquola C, Bethlehem RAI, et al. Adolescent development of multiscale structural wiring and functional interactions in the human connectome. Proc Natl Acad Sci U S A. 2022;119(27):e2116673119. doi:10.1073/pnas.2116673119
    OpenUrlCrossRefPubMed
  40. 40.↵
    Cruces RR, Royer J, Herholz P, et al. Micapipe: A pipeline for multimodal neuroimaging and connectome analysis. NeuroImage. 2022;263:119612. doi:10.1016/j.neuroimage.2022.119612
    OpenUrlCrossRef
  41. 41.↵
    Dale AM, Fischl B, Sereno MI. Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage. 1999;9(2):179-194. doi:10.1006/nimg.1998.0395
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    Patenaude B, Smith SM, Kennedy DN, Jenkinson M. A Bayesian model of shape and appearance for subcortical brain segmentation. NeuroImage. 2011;56(3):907–922. doi:10.1016/j.neuroimage.2011.02.046
    OpenUrlCrossRefPubMedWeb of Science
  43. 43.↵
    Tournier JD, Smith R, Raffelt D, et al. MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation. Neuroimage. 2019;202:116137. doi:10.1016/j.neuroimage.2019.116137
    OpenUrlCrossRefPubMed
  44. 44.↵
    Desikan RS, Ségonne F, Fischl B, et al. An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage. 2006;31(3):968–980. doi:10.1016/j.neuroimage.2006.01.021
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    Vos de Wael R, Benkarim O, Paquola C, et al. BrainSpace: a toolbox for the analysis of macroscale gradients in neuroimaging and connectomics datasets. Commun Biol. 2020;3(1):1-10. doi:10.1038/s42003-020-0794-7
    OpenUrlCrossRefPubMed
  46. 46.↵
    Betzel RF, Griffa A, Hagmann P, Mišić B. Distance-dependent consensus thresholds for generating group-representative structural brain networks. Netw Neurosci. 2019;3(2):475–496. doi:10.1162/netn_a_00075
    OpenUrlCrossRef
  47. 47.↵
    Coifman RR, Lafon S. Diffusion maps. Applied and Computational Harmonic Analysis. 2006;21(1):5-30. doi:10.1016/j.acha.2006.04.006
    OpenUrlCrossRef
  48. 48.↵
    Hong SJ, Xu T, Nikolaidis A, et al. Toward a connectivity gradient-based framework for reproducible biomarker discovery. Neuroimage. 2020;223:117322. doi:10.1016/j.neuroimage.2020.117322
    OpenUrlCrossRef
  49. 49.↵
    Larivière S, Bayrak Ş, Vos de Wael R, et al. BrainStat: A toolbox for brain-wide statistics and multimodal feature associations. Neuroimage. 2023;266:119807. doi:10.1016/j.neuroimage.2022.119807
    OpenUrlCrossRef
  50. 50.↵
    Worsley KJ, Taylor J, Carbonell F, et al. A Matlab toolbox for the statistical analysis of univariate and multivariate surface and volumetric data using linear mixed effects models and random field theory. In: NeuroImage Organisation for Human Brain Mapping 2009 Annual Meeting. Vol 47. ; 2009:S102.
  51. 51.↵
    Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society: Series B (Methodological). 1995;57(1):289–300. doi:10.1111/j.2517-6161.1995.tb02031.x
    OpenUrlCrossRefPubMed
  52. 52.↵
    Park B yong, Bethlehem RA, Paquola C, et al. An expanding manifold in transmodal regions characterizes adolescent reconfiguration of structural connectome organization. Uddin LQ, Behrens TE, Fulcher BD, eds. eLife. 2021;10:e64694. doi:10.7554/eLife.64694
    OpenUrlCrossRef
  53. 53.↵
    Bethlehem RAI, Paquola C, Seidlitz J, et al. Dispersion of functional gradients across the adult lifespan. Neuroimage. 2020;222. doi:10.1016/j.neuroimage.2020.117299
    OpenUrlCrossRefPubMed
  54. 54.↵
    Burt JB, Helmer M, Shinn M, Anticevic A, Murray JD. Generative modeling of brain maps with spatial autocorrelation. Neuroimage. 2020;220:117038. doi:10.1016/j.neuroimage.2020.117038
    OpenUrlCrossRef
  55. 55.↵
    Leech R, Vos De Wael R, Váša F, et al. Variation in spatial dependencies across the cortical mantle discriminates the functional behaviour of primary and association cortex. Nat Commun. 2023;14(1):5656. doi:10.1038/s41467-023-41334-2
    OpenUrlCrossRef
  56. 56.↵
    DeKraker J, Haast RA, Yousif MD, et al. Automated hippocampal unfolding for morphometry and subfield segmentation with HippUnfold. eLife. 2022;11:e77945. doi:10.7554/eLife.77945
    OpenUrlCrossRef
  57. 57.↵
    McIntosh AR, Mišić B. Multivariate Statistical Analyses for Neuroimaging Data. Annu Rev Psychol. 2013;64(1):499–525. doi:10.1146/annurev-psych-113011-143804
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    Krishnan A, Williams LJ, McIntosh AR, Abdi H. Partial Least Squares (PLS) methods for neuroimaging: a tutorial and review. Neuroimage. 2011;56(2):455–475. doi:10.1016/j.neuroimage.2010.07.034
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    Kirschner M, Shafiei G, Markello RD, et al. Latent Clinical-Anatomical Dimensions of Schizophrenia. Schizophr Bull. 2020;46(6):1426–1438. doi:10.1093/schbul/sbaa097
    OpenUrlCrossRef
  60. 60.↵
    Valk SL, Xu T, Margulies DS, et al. Shaping brain structure: Genetic and phylogenetic axes of macroscale organization of cortical thickness. Sci Adv. 2020;6(39):eabb3417. doi:10.1126/sciadv.abb3417
    OpenUrlFREE Full Text
  61. 61.↵
    Wang XJ. Macroscopic gradients of synaptic excitation and inhibition in the neocortex. Nat Rev Neurosci. 2020;21(3):169–178. doi:10.1038/s41583-020-0262-x
    OpenUrlCrossRef
  62. 62.↵
    Vogel JW, La Joie R, Grothe MJ, et al. A molecular gradient along the longitudinal axis of the human hippocampus informs large-scale behavioral systems. Nat Commun. 2020;11(1):960. doi:10.1038/s41467-020-14518-3
    OpenUrlCrossRefPubMed
  63. 63.↵
    Huntenburg JM, Bazin PL, Margulies DS. Large-Scale Gradients in Human Cortical Organization. Trends in Cognitive Sciences. 2018;22(1):21–31. doi:10.1016/j.tics.2017.11.002
    OpenUrlCrossRefPubMed
  64. 64.↵
    Paquola C, Vos De Wael R, Wagstyl K, et al. Microstructural and functional gradients are increasingly dissociated in transmodal cortices. PLoS Biol. 2019;17(5):e3000284. doi:10.1371/journal.pbio.3000284
    OpenUrlCrossRefPubMed
  65. 65.↵
    Vázquez-Rodríguez B, Suárez LE, Markello RD, et al. Gradients of structure-function tethering across neocortex. Proc Natl Acad Sci U S A. 2019;116(42):21219–21227. doi:10.1073/pnas.1903403116
    OpenUrlAbstract/FREE Full Text
  66. 66.↵
    Fadaie F, Lee HM, Caldairou B, et al. Atypical functional connectome hierarchy impacts cognition in temporal lobe epilepsy. Epilepsia. 2021;62(11):2589–2603. doi:10.1111/epi.17032
    OpenUrlCrossRef
  67. 67.↵
    Xie K, Royer J, Larivière S, et al. Atypical connectome topography and signal flow in temporal lobe epilepsy. bioRxiv. Published online May 24, 2023:2023.05.23.541934. doi:10.1101/2023.05.23.541934
    OpenUrlAbstract/FREE Full Text
  68. 68.↵
    Ahmadi ME, Hagler DJ, McDonald CR, et al. Side matters: diffusion tensor imaging tractography in left and right temporal lobe epilepsy. AJNR Am J Neuroradiol. 2009;30(9):1740–1747. doi:10.3174/ajnr.A1650
    OpenUrlAbstract/FREE Full Text
  69. 69.↵
    Concha L, Beaulieu C, Collins DL, Gross DW. White-matter diffusion abnormalities in temporal-lobe epilepsy with and without mesial temporal sclerosis. J Neurol Neurosurg Psychiatry. 2009;80(3):312–319. doi:10.1136/jnnp.2007.139287
    OpenUrlAbstract/FREE Full Text
  70. 70.↵
    Thivard L, Lehéricy S, Krainik A, et al. Diffusion tensor imaging in medial temporal lobe epilepsy with hippocampal sclerosis. Neuroimage. 2005;28(3):682–690. doi:10.1016/j.neuroimage.2005.06.045
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    Bernhardt BC, Bernasconi N, Hong SJ, Dery S, Bernasconi A. Subregional Mesiotemporal Network Topology Is Altered in Temporal Lobe Epilepsy. Cerebral Cortex. 2016;26(7):3237–3248. doi:10.1093/cercor/bhv166
    OpenUrlCrossRefPubMed
  72. 72.↵
    Khambhati AN, Davis KA, Oommen BS, et al. Dynamic Network Drivers of Seizure Generation, Propagation and Termination in Human Neocortical Epilepsy. PLoS Comput Biol. 2015;11(12):e1004608. doi:10.1371/journal.pcbi.1004608
    OpenUrlCrossRefPubMed
  73. 73.↵
    Larivière S, Weng Y, Vos de Wael R, et al. Functional connectome contractions in temporal lobe epilepsy: Microstructural underpinnings and predictors of surgical outcome. Epilepsia. 2020;61(6):1221–1233. doi:10.1111/epi.16540
    OpenUrlCrossRef
  74. 74.↵
    Liu M, Gross DW, Wheatley BM, Concha L, Beaulieu C. The acute phase of Wallerian degeneration: longitudinal diffusion tensor imaging of the fornix following temporal lobe surgery. Neuroimage. 2013;74:128–139. doi:10.1016/j.neuroimage.2013.01.069
    OpenUrlCrossRefPubMed
  75. 75.↵
    Winston GP, Stretton J, Sidhu MK, Symms MR, Duncan JS. Progressive white matter changes following anterior temporal lobe resection for epilepsy. Neuroimage Clin. 2014;4:190–200. doi:10.1016/j.nicl.2013.12.004
    OpenUrlCrossRefPubMed
  76. 76.↵
    Concha L, Beaulieu C, Wheatley BM, Gross DW. Bilateral white matter diffusion changes persist after epilepsy surgery. Epilepsia. 2007;48(5):931–940. doi:10.1111/j.1528-1167.2007.01006.x
    OpenUrlCrossRefPubMed
  77. 77.↵
    McDonald CR, Hagler DJ, Girard HM, et al. Changes in fiber tract integrity and visual fields after anterior temporal lobectomy. Neurology. 2010;75(18):1631–1638. doi:10.1212/WNL.0b013e3181fb44db
    OpenUrlAbstract/FREE Full Text
  78. 78.↵
    Bonilha L, Gleichgerrcht E, Fridriksson J, et al. Reproducibility of the Structural Brain Connectome Derived from Diffusion Tensor Imaging. PLoS One. 2015;10(8):e0135247. doi:10.1371/journal.pone.0135247
    OpenUrlCrossRefPubMed
  79. 79.↵
    Llera A, Wolfers T, Mulders P, Beckmann CF. Inter-individual differences in human brain structure and morphology link to variation in demographics and behavior. eLife. 2019;8:e44443. doi:10.7554/eLife.44443
    OpenUrlCrossRef
  80. 80.↵
    Mueller S, Wang D, Fox MD, et al. Individual Variability in Functional Connectivity Architecture of the Human Brain. Neuron. 2013;77(3):586–595. doi:10.1016/j.neuron.2012.12.028
    OpenUrlCrossRefPubMedWeb of Science
  81. 81.↵
    Keller SS, Glenn GR, Weber B, et al. Preoperative automated fibre quantification predicts postoperative seizure outcome in temporal lobe epilepsy. Brain. 2017;140(1):68–82. doi:10.1093/brain/aww280
    OpenUrlCrossRefPubMed
  82. 82.
    Bonilha L, Jensen JH, Baker N, et al. The brain connectome as a personalized biomarker of seizure outcomes after temporal lobectomy. Neurology. 2015;84(18):1846–1853. doi:10.1212/WNL.0000000000001548
    OpenUrlAbstract/FREE Full Text
  83. 83.
    Gleichgerrcht E, Munsell B, Bhatia S, et al. Deep learning applied to whole-brain connectome to determine seizure control after epilepsy surgery. Epilepsia. 2018;59(9):1643–1654. doi:10.1111/epi.14528
    OpenUrlCrossRefPubMed
  84. 84.
    Morgan VL, Sainburg LE, Johnson GW, et al. Presurgical temporal lobe epilepsy connectome fingerprint for seizure outcome prediction. Brain Communications. 2022;4(3):fcac128. doi:10.1093/braincomms/fcac128
    OpenUrlCrossRef
  85. 85.
    Larivière S, Bernasconi A, Bernasconi N, Bernhardt BC. Connectome biomarkers of drug-resistant epilepsy. Epilepsia. 2021;62(1):6–24. doi:10.1111/epi.16753
    OpenUrlCrossRef
  86. 86.↵
    Bernhardt BC, Hong SJ, Bernasconi A, Bernasconi N. Magnetic resonance imaging pattern learning in temporal lobe epilepsy: Classification and prognostics. Annals of Neurology. 2015;77(3):436–446. doi:10.1002/ana.24341
    OpenUrlCrossRefPubMed
  87. 87.↵
    Lee HM, Fadaie F, Gill R, et al. Decomposing MRI phenotypic heterogeneity in epilepsy: a step towards personalized classification. Brain. 2022;145(3):897–908. doi:10.1093/brain/awab425
    OpenUrlCrossRef
Back to top
Posted November 13, 2023.
Download PDF

Supplementary Material

Data/Code
Email
Connectome reorganization associated with temporal lobe pathology and its surgical resection
Sara Larivière, Bo-yong Park, Jessica Royer, Jordan DeKraker, Alexander Ngo, Ella Sahlas, Judy Chen, Raúl Rodríguez-Cruces, Yifei Weng, Birgit Frauscher, Ruoting Liu, Zhengge Wang, Golia Shafiei, Bratislav Mišić, Andrea Bernasconi, Neda Bernasconi, Michael D. Fox, Zhiqiang Zhang, Boris C. Bernhardt
medRxiv 2023.11.13.23298482; doi: https://doi.org/10.1101/2023.11.13.23298482
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area