Abstract
Background and Objective. Targeted speech therapy can lead to substantial naming improvement in some subjects with anomia following dominant-hemisphere stroke. We investigated whether treatment-induced improvement in naming is associated with poststroke preservation of structural neural network architecture. Methods. Twenty-four patients with poststroke chronic aphasia underwent 30 hours of speech therapy over a 2-week period and were assessed at baseline and after therapy. Whole brain maps of neural architecture were constructed from pretreatment diffusion tensor magnetic resonance imaging to derive measures of global brain network architecture (network small-worldness) and regional network influence (nodal betweenness centrality). Their relationship with naming recovery was evaluated with multiple linear regressions. Results. Treatment-induced improvement in correct naming was associated with poststroke preservation of global network small worldness and of betweenness centrality in temporal lobe cortical regions. Together with baseline aphasia severity, these measures explained 78% of the variability in treatment response. Conclusions. Preservation of global and left temporal structural connectivity broadly explains the variability in treatment-related naming improvement in aphasia. These findings corroborate and expand on previous classical lesion-symptom mapping studies by elucidating some of the mechanisms by which brain damage may relate to treated aphasia recovery. Favorable naming outcomes may result from the intact connections between spared cortical areas that are functionally responsive to treatment.
Keywords
Introduction
Stroke is one of the leading causes of disability worldwide and in the United States. 1 Aphasia is an impairment of language processing that often results from strokes affecting the dominant brain hemisphere. 2 A large number of subjects who suffered a left hemisphere stroke remain with some degree of long-term aphasia,2,3 and this persistent communication impairment is commonly recognized as one of the most debilitating sequelae from stroke. 4 One of the hallmark deficits of aphasia is the inability to name objects or people (anomia) 5 ; the severity of anomia is closely related to a poor quality of life.6,7
While many patients with aphasia exhibit some recovery of naming in the first days to weeks after the stroke, 30% to 40% persist with long-lasting naming impairments. 2 Speech therapy can be effective to improve language processing for some subjects with chronic aphasia.8-11 However, as is well recognized among clinicians, speech therapy can be ineffective for many others, and there are no well-defined means to predict which subjects may benefit from treatment. For example, in spite of decades of research, a relationship between linguistic factors and aphasia treatment success has not been revealed, suggesting that baseline behavioral testing may not provide clinicians the necessary prognostic indicators. 12 Accordingly, therapy management could be improved by a better understanding of the neurobiological mechanisms associated with aphasia and naming recovery.
Our group has demonstrated that patients with lesions encompassing areas commonly associated with lexical retrieval and phonological processing, such as Brodmann’s areas (BA) 37 and 39, are less likely to show treatment-related improvement in correct naming. 13 Furthermore, the functional modulation (ie, brain activation) of the perilesional left frontal and left temporal cortices are also associated with an increase in the number of items named correctly, 14 suggesting that the apparent preservation of cortical regions after the stroke is not, by itself, sufficient to enable recovery. Even though posttreatment brain activation among patients with aphasia differs from baseline, 15 some patients with apparently intact temporal or frontal cortical structures fail to achieve a good recovery. 16
Contemporary advances in structural brain mapping, particularly those employing structural magnetic resonance imaging (MRI), have greatly expanded our knowledge regarding the relationship between brain connectivity, neural network integrity, and language. Specifically, whole brain MRI diffusion tensor imaging (DTI) can be quantified to provide a map of neural architecture for each person, that is, the individual brain connectome. 17 Adapting connectome methods to the study of neural architecture in persons who suffered a stroke, our group has recently demonstrated that white matter destruction can lead to disconnection of remote and apparently intact cortical structures. 18 The destruction of white matter and subsequent disruption of brain connectivity are strong predictors of poststroke naming impairment. 19 Indeed, a compromised neural network can be as important as direct cortical necrosis in the ensuing severity of naming problems. 19
Thus, it is likely that changes occurring within white matter tracts leading to cortical disconnection may account for some of the variability in language recovery previously demonstrated by lesion symptom mapping studies. Specifically, preserved cortical connectivity may permit the functional recruitment of areas modulated by therapy and partially explain the structural framework for recovery.
This hypothesis is consistent with recent findings relating the importance of white matter tracts for language impairment and recovery. 20 Previous studies have shown that lesion load on language-relevant white matter tracts (namely, the arcuate and uncinate fasciculi) contribute to recovery, with changes between pre- and posttreatment of anomia.21-24 Abnormal patterns and posttreatment changes have also been shown to occur at the level of the functional connectome,25,26 which describes the synchronous activity of regions in an individual brain. All of the aforesaid results have encouraged the development of interventions incorporating noninvasive brain stimulation (eg, repetitive transcranial magnetic stimulation), which have preliminarily showed positive immediate-, short-, and long-term effects on naming performance.27-29
We hypothesized that the poststroke integrity of neural network architecture is a crucial determinant for naming recovery motivated by speech therapy. Here, we describe a cohort of patients with chronic aphasia and naming impairment after a stroke, who underwent a controlled period of intense speech therapy aimed at improving naming performance. Part of this cohort was previously examined by our group in a classical lesion-symptom mapping study. 30 Here, we complement that study by investigating whether direct cortical disconnection and neural architecture are the mechanisms associated with poststroke lesions that lead to impairments and permit recovery. All subjects were carefully assessed regarding their language abilities before and after the treatment. Instead of attempting to match treatment type to each individual’s cognitive-linguistic profile, all participants underwent anomia treatment using a protocol that has been amply described in the literature—semantic and phonological cueing hierarchies. 31 The rationale behind this approach is that the relationship between behavioral impairment, optimal treatment type, and treatment success is, at best, very poorly understood at this time. Much like traditional studies that employed lesion-symptom mapping to reveal the relationship between behavioral impairment and structural brain damage, the approach taken here assumes that residual neuroanatomy (ie, structural connectivity) relates to patients’ ability to improve with the aid of anomia treatment.
We conducted a detailed assessment of global and regional neural network architecture prior to treatment, with the intention to evaluate the relationship between the degree of poststroke network disruption and the ability to achieve improvement after speech therapy.
Methods
Subjects
We studied 24 individuals recruited from the local community with chronic aphasia as a result of a left hemisphere ischemic stroke that occurred at least 6 months prior to enrolling in this study. All subjects were right handed prior to the stroke and had no history of other neurological diseases. They also did not have a history or imaging evidence of other strokes. Patients with a history of medication-refractory seizures were not included in this study.
The behavioral data from 20 patients reported here were also included, in part, in 2 previous studies from our group.13,14 Those studies examined a subject population that overlapped with the subjects reported here. The question addressed here, the relationship between naming recovery and the integrity of structural brain networks, was not addressed in those studies.
All subjects signed an informed consent to participate in this study. They were recruited and tested at The University of South Carolina (USC). The USC Institutional Review Board approved this study.
Baseline Language Testing
A certified speech pathologist administered and scored language performance at baseline based on the Western Aphasia Battery (WAB) 32 and the Philadelphia Naming Test (PNT). 33 Together, these tests provide a measure of (a) aphasia type and severity (WAB) and (b) naming performance (PNT). Specifically, aphasia severity was defined in accordance with the WAB aphasia quotient (AQ), which is a measure of aphasia severity ranging from 0 to 100. The PNT is a computer-based assessment of naming, 33 which includes 175 pictures representing mid- and high-frequency nouns from a word frequency list compiled by Francis and colleagues. 34
Prior to therapy, all subjects underwent 2 independent complete testing sessions, no further than 2 to 3 days apart, in order to account for daily variability in performance. The baseline performance scores on the PNT were defined as the average of the 2 pretreatment testing sessions.
Baseline MRI Scanning
All subjects underwent MRI scanning at USC in a 3T Siemens Trio equipped with a 12-channel head coil in accordance with the following protocol: (a) T1-weighted images (3D MP-RAGE, TR = 2250 ms, TE = 4.15 ms, 256 × 256 matrix, 256 × 256 mm FOV, parallel imaging GRAPPA = 2, 80 reference lines, TA = 377 s); (b) Diffusion EPI scan (30-directions with b = 1000 s/mm2 and b = 2000 s/mm2, TR = 6100 ms, TE = 101 ms, 82 × 82 matrix, 222 × 222mm FOV, parallel imaging GRAPPA = 2, 80 45 contiguous 2.7 mm axial slices, TA = 390 s).
MRI scanning was performed at baseline before language therapy was applied.
Language Therapy
After baseline testing was concluded, all subjects underwent a period of speech therapy aimed at improving naming performance. Therapy was performed daily for a period of 2 weeks. During each 3-hour treatment session, the subject underwent clinician-administered anomia treatment using semantic and phonemic cueing hierarchies. 31 As indicated earlier, the treatment approach used here was not selected based on the cognitive-linguistic impairment profile of any single patient or a subgroup of patients. Rather, it was chosen given that studies that have employed similar methods have revealed its potential in improving naming across a large group of patients. 35 All subjects completed the treatment as planned without missing sessions.
The nouns targeted in treatment (n = 160) were not included on the PNT. Thus, naming improvement was qualified as an increase in naming of untrained items. During this time, none of the subjects participated in any other form of speech therapy, experienced any significant change in health status, or underwent any medication adjustments.
Posttherapy Language Testing
After the conclusion of treatment, all subjects were reevaluated with the PNT. This posttesting replicated the pretreatment evaluations, and was also scored by a speech language pathologist. The evaluation occurred no longer than 1 week after treatment completion. The 2 posttreatment evaluations occurred within 2 to 3 days of each other, scores again being the average of these posttherapy tests.
Preprocessing of MRI Data and Connectome Construction
The preprocessing steps involved in the construction of the individual connectome involved the following: (a) segmentation of the probabilistic gray matter map (from T1-weighted images), further divided into cortical anatomical regions of interest (ROIs) corresponding to BAs; (b) segmentation of the probabilistic white matter map from T1-weighted images; (c) registration of the individual white matter map and cortical ROIs into the individual DTI space; (d) probabilistic DTI fiber tracking; (e) iterative evaluation of the number of tractography streamlines connecting each possible pair of ROIs; (f) correction of each pair-wise connection strength (ie, number of streamlines) based on the volume of the connected ROIs and the distance travelled by the streamlines. A summary of these steps is demonstrated in Figures 1 and 2.

This figure summarizes the initial preprocessing steps employed to evaluate structural connectivity. Images from one representative subject are shown. The T1-weighted images are segmented into probabilistic maps of gray matter (panel A) and white matter (panel B). The probabilistic map of gray matter is then divided into anatomical regions of interest (ROIs) corresponding to Brodmann areas (panel C—each color representing a different ROI).

Following the initial preprocessing steps, special processing is then conducted to exclude the areas involved in the poststroke necrotic tissue, thus leaving only viable cortex (panel A—areas in yellow) and viable white matter (panel B—areas in green). Using this probabilistic white matter map as a guide, probabilistic white matter diffusion tensor imaging streamlines are reconstructed by seeding each region of interest (ROI) at a time, and computing the number of streamlines arriving at the other ROIs. Panel C illustrates the voxel-based counts of probabilistic streamlines generated by seeding Brodmann areas 6 (premotor cortex, with streamlines shown in “warm” color) and 48 (insula, streamlines in “cold” color) in the right hemisphere. The scale bars illustrate the voxel-based counts of streamlines.
We employed methods previously described by our group 19 to overcome the distortions caused by the stroke necrotic lesion, aiming to preserve the anatomical authenticity of gray and white matter tissues, while avoiding the computation of fiber tracking in areas corresponding to the necrotic tissue. The area of necrosis was manually drawn on T2 images. The 3-dimensional representation of this lesion was then used to calculate the lesion volume (number of voxels × voxel volume) for the purposes of our statistical analyses (see below). These methods are described in detail in the Supplemental Methods.
In summary, for each subject, the structural connectome was defined as a connectivity matrix
Measures of Neural Network Architecture
From the individual connectivity matrix, we derived a measure of the integrity of the global network, namely, normalized small worldness (NSW) which was calculated in accordance with the methods described by Brown and collaborators. 36 Specifically, for each subject’s structural network, 1000 random networks (of similar link weight distribution) were generated by means of a computerized algorithm. We then calculated the global clustering coefficient of each network, which is an average of how each node in the network is connected to its neighbors (ie, a proportion of how many links actually exist between a node and its neighbors and the maximum number of links that could exist). We then computed a ratio of the clustering coefficient from the real network to the average clustering coefficient from the random networks was calculated as the normalized clustering coefficient. We applied a similar procedure to calculate the ratio of characteristic path length (ie, the number of links counted in the shortest path to transverse all possible pairs of network nodes) between the real network and the 1000 random networks. NSW was calculated as the ratio between the normalized clustering and the normalized characteristic path length, thus representing one unique score per participant reflecting global network integrity.
Involvement of Temporal and Frontal Nodes in the Global Network
In order to assess the global network participation of temporal, parietal, and frontal nodes, we calculated the average nodal betweenness centrality (BC) of a series of ROIs. BC was calculated employing functions from the Brain Connectivity Toolbox. 37 BC is a measure of regional (ie, local, as opposed to global) network integration, and it is related to how many of the network’s shortest paths travel through the node, providing a measure of influence of each node into the network. 37 As stated above, a shortest path is the minimal network distance between each possible pair of nodes. Nodes with high BC are commonly traversed by shortest paths linking 2 other nodes, being thus very influential and central to the network. We selected left-sided ROIs based on recent models of speech processing anatomy, 38 and also based on conventional lobar neuroanatomy. Thus, ROIs were grouped as follows: left frontal ROIs (BA = [44, 45, 46]), left parietal ROIs (BA = [7, 39, 40, 41]), and left temporal ROIs (BA = [20, 21, 22, 37, 38]). These ROIs were selected to represent the majority of the dorsolateral aspects of the temporal and parietal lobes, as well as the posterior aspects of the inferior and middle temporal gyri. This approach permits an assessment of large representations of language areas within the left hemisphere, albeit with limited resolution to detect abnormalities within smaller or more localized subregions. This process was repeated to assess nodal BC of right hemisphere ROIs corresponding to the homologous BAs of the ROIs described above.
Statistical Analyses
We investigated the relationship between increase in correct naming and measures of neural architecture. We controlled for other confounding variables that could independently influence network architecture or naming improvement, such as age (at the time of the study enrollment), time since the stroke, overall aphasia severity (WAB-AQ), and the volume of the necrotic stroke lesion.
Treatment-related improvement in naming was defined based on the difference between pre- and posttreatment performance on the PNT, adjusted based on the potential for recovery, as follows: the number of new items correctly named on the PNT was divided by the total number of items in the PNT [175] minus the number of correctly named items at baseline. Thus, an improvement in naming 20 items was weighted differently if the baseline score had been 155 ([20/(175 − 155)] = 100% improvement), versus 75 ([20/(175 − 75)] = 20% improvement).
The relationship between naming improvement and neural architecture was investigated through multiple linear regressions with treatment-related improvement in naming set as the dependent variable, and measures of neural architecture along with controlling variables set as independent variables. The explanatory power of the resulting regression model was determined as R2 (proportional reduction in error) with the explanatory factors of interest entered into the regression analysis using a stepwise approach. Control for multiple comparisons was performed through Bonferroni correction: multiple linear regression analyses were performed to evaluate the relationship between the following variables: (a) baseline WAB and left hemisphere BCs; (b) baseline WAB and NSW; (c) PNT improvement and left hemisphere BCs; (d) PNT improvement and NSW; (e) PNT improvement, left BCs, and NSW; (f) PNT improvement and right and left BCs. Therefore, the level of statistical significance of a given model was set at P < .05/6 = .008.
Results
Demographics
The mean age of the study sample was 60.2 ± 11.9 years. Twelve subjects were male and the mean age at the time of the stroke was 55.5 ± 12.1 years. Twenty-one subjects were Caucasian and 3 were African American (Demographic details are further described in Supplemental Table 1).
Pretreatment Language Performance
As the purpose of this study was to determine the association between treatment-related improvement in correct naming and structural connectivity in aphasic patients, it was essential that the patient sample reflect a wide range of aphasia severities and types as well as variability in lesion location and extent of injury. The mean AQ for the group of subjects was 57.94 (SD = 28.27). Based on the WAB subscores, the aphasia types were: anomic, 41% (10/24); Broca’s, 33% (8/24); Wernicke’s, 8% (2/24); conduction, 8% (2/24); and global, 8% (2/24). At baseline, combining the 2 pretreatment testing sessions, among all subjects, the average number of correctly named items on the PNT was 59 ± 61 items, ranging from 0 to 154 items. The average number of attempts (ie, correct and error responses) was 118 ± 61 items, and the average number of no responses was 57 ± 61 items. Semantic paraphasias accounted for 9.6 ± 7.3 items, while phonemic paraphasias comprised 11.3 ± 14.7 items.
Posttreatment Language Performance
Subjects improved, on average, by correctly naming 6.8 ± 7.1 more items after treatment (taking into account the average between the 2 posttreatment sessions in comparison with the average from the 2 pretreatment sessions). The treatment-related improvement in naming, as defined above by ratio between number of correctly named items and the number of pretreatment incorrect responses, was 12.3 ± 15.6%, ranging from −4.5% (ie, worsening after treatment) to 20.5%. There was also an increase in the number of response attempts (12.7 ± 18.6 items) after treatment. Interestingly, on average, phonemic paraphasias increased by 1.9 ± 7.2 items, while semantic paraphasias remained unchanged (0 ± 5.5 items). No significant correlation was found between changes in correct naming and time (months) poststroke (r = .08, P = .71).
The pretreatment and posttreatment language performance data are presented in Supplemental Table 1.
Stroke Lesion Variability
All patients included in this study exhibited a chronic necrotic lesion corresponding to sequelae of liquefactive gliosis due to ischemic injury in the cortical or subcortical vascular territory perfused by the left middle cerebral artery. In all subjects, the lesion was clearly visible on T2-weighted images. As demonstrated by the overlap between lesion sites (Figure 3), the most common location of damage was the left insular cortex, followed by the left temporal and the left frontal opercular regions, and then by other cortical areas in the left dorsolateral frontal, parietal, and temporal regions. This overlap pattern is compatible with previously described lesion distributions in patients with aphasia. 39

The anatomical distribution of the necrotic poststroke brain lesions is demonstrated through a color-coded map (overlaid onto a cortical surface diagram) where each color represents the frequency, across all subjects, with which the corresponding anatomical region was lesioned.
Structural Connectomes and Left Hemisphere Connectivity Abnormalities
The structural connectome from each subject is demonstrated in Supplemental Figure 1. The average group connectome is demonstrated in Figure 4A. Even though the location of the necrotic lesion is heterogeneous among the study subjects, overall, the connections within the left hemisphere are less numerous or exhibit lower strength compared with connections within the right hemisphere. Figure 4A also demonstrates that the connectome can be abridged into a matrix where each entry corresponds to the connectivity strength between the row and column ROIs.

This diagram demonstrates the average connectome in panel A, with axis numbers corresponding to Brodmann’s areas (Bas). The scale bar demonstrates the link-wise strength, which corresponds to the log of the number of streamlines connecting the regions of interest (ROIs; corrected based on ROI volume and distance travelled by the streamlines). Panel B demonstrates the link-wise asymmetry in connectivity, illustrating how often links in the left hemisphere were present, in comparison with the homologous right hemisphere links that were present in 100% of subjects. A sphere in its center of mass illustrates the location of each BA ROI; straight lines connecting the nodes demonstrate links; the shade of gray of each link corresponds to its frequency of observation in comparison with the right hemisphere. The colors of the nodes illustrate their anatomical locations (light blue—parietal and inferior frontal; dark blue—frontal; yellow—temporal; dark red—pericingulate regions and insula; green—occipital).
Figure 4B demonstrates a quantification of the relative frequency of disconnection within the left hemisphere, shown as the relative frequency with which links in the left hemisphere were observed, including only links whose right hemisphere homologous was observed in all patients (ranging from 0—never to 100%—in all patients). As expected, less connectivity in the left hemisphere was observed involving areas with highest necrotic tissue overlay (insular, opercula, and dorsolateral regions).
Relationship Between Baseline Language Performance and Structural Connectivity
We observed a relationship trend between more severe aphasia (measured as lower WAB-AQ) and lower left temporal BC (F = 3.97, P = .011, model R = .76, model R2 = .58), controlling for age, time since the stroke, and stroke lesion size. However, this model was not statistically significant due to the control for multiple statistical comparisons. This trend suggests that loss of influence from left temporal nodes onto the general network was associated with more severe forms of aphasia at baseline (Figure 5 and Supplemental Table 2).

Scatter plots demonstrating the results from the linear regression models in relationship with baseline Western Aphasia Battery–aphasia quotient (WAB-AQ; upper panel) and treatment-related improvement in the naming (Philadelphia Naming Test [PNT] improvement—middle and lower panels) set as dependent variables. The adjusted models (controlling for age, time since the stroke, and lesion size) are shown on the left column, and the association between the independent variables related to global network measures (NSW—normalized small worldness) and regional left temporal network measures (BC—betweenness centrality) are shown on the right column. Supplemental Table 2 provides a comprehensive description of the statistical features of these models.
We did not observe a relationship between baseline aphasia severity and global network NSW (F = 2.7, P = .06; Supplemental Table 2). Moreover, NSW was not correlated with lesion volume (r = −.04, P = .9).
Relationship Between Naming Improvement and Regional and Global Structural Connectivity
We observed a robust relationship between treatment-related improvement in naming and a higher left temporal BC (controlling for age, time since the stroke, stroke lesion size, and WAB-AQ) (model F = 6.79, P = .00077, model R = .86, model R2 = .74; Figure 5, Supplemental Table 2). In this model, WAB-AQ was also associated with treatment-related improvement in naming (T = 2.45, P = .03); patients with less severe aphasia achieved a higher naming improvement.
We did not observe a significant relationship between treatment-related improvement in naming and regional network measures from the frontal or parietal areas (Supplemental Table 2). Lesion volume was not associated with treatment-induced improvements in naming (r = −.28, P = .18). Moreover, there was not a direct correlation between lesion burden in the temporal lobes and the network integration of the remaining temporal lobe cortex (BC; r = −.01, P = .9).
We also observed a robust relationship between PNT improvement and global network NSW. Patients with higher global NSW demonstrated greater treatment-related improvement in naming (controlling for age, time since the stroke, stroke lesion size, and WAB-AQ) (model F = 8.89, P = .0002, model R = .84, model R2 = .712; Figure 5, Supplemental Table 2). Similarly, WAB-AQ was also associated with PNT improvement.
Finally, if NSW and left temporal BC were both included as independent variables in the same model (again, controlling for age, time since the stroke, stroke lesion size, and WAB-AQ), this new model explained 78% of the variance observed in improvement in the PNT (F = 10.2, P = .00007, model R = .88, model R2 = .78). We did not evaluate the individual influence of each predictor in this model since left temporal BC and NSW were highly correlated (R = .62, P = .001), leading to inaccuracies in individual independent variable contribution due to multicolinearity.
The linear regression analysis combining nodal BC from left and right hemisphere BAs as independent variables demonstrated an association with improvement in the PNT (F = 5.62, P = .003, model R = .82, model R2 = .68); nonetheless, none of the right hemisphere regions (ie, right frontal, right temporal, and right parietal ROIs) were independently associated with improvement.
Discussion
In this study, we examined the relationship between regional and global brain neural network architecture in the chronic state after an ischemic stroke, and the ability to improve in naming performance after treatment in subjects with aphasia. Our findings can be summarized as described below.
First, preserved global network architecture (NSW) is associated with a more pronounced treatment-related improvement in naming. Second, if temporal lobe regions remain integrated with the remaining network (ie, temporal lobe BC), a higher treatment-related improvement is also more likely to occur. These 2 measures (NSW and left temporal BC) are highly associated with each other, but independently, each one of these measures could predict more than 70% of the variance in naming improvement. Third, subjects with more severe aphasia prior to treatment (lower WAB-AQ) are also less likely to improve with treatment. Combined with NSW and left temporal BC, a model composed by these 3 measures could predict 78% of the treatment-related improvement. Fourth, we did not observe a relationship between structural integrity of other brain regions (frontal or parietal) and treatment-related improvement.
Treatment-Related Improvement Is Still Possible in the Chronic Stages of Aphasia
Even though the prognosis for recovery for some patients with aphasia can be unfavorable, our study corroborates previous observations reporting that a subset of subjects may still experience recovery in the chronic stages of aphasia.8,40-42 As it has been suggested in treatment studies of aphasia,13,14 the identification of patients with the potential to recover could lead to a positive impact in the clinical management of aphasia, enabling therapy to a targeted population, thus providing access to speech therapy that is not routinely offered in the chronic poststroke stage, and the possibility of further recovery.
So far, the crucial limitation in this approach has been the identification of those patients who can actually improve. Besides aphasia severity, there are probably few behavioral variables that can discriminate patients with a good prognosis as accurately. For example, recovery has not been shown be related to age during therapy, age at the time of the stroke, necrotic lesion size, or even aphasia characteristics.43,44 Nonetheless, the patients who benefit from anomia therapy seem to have in common some basic neurobiological apparatus that enables improved naming as a result of treatment. Our findings suggest that the integrity of structural brain networks (eg, as reflected by NSW and BC measures) is a biomarker for the mechanisms that support treated aphasia recovery. As such, it can be quantified and potentially used in clinical settings to better understand and guide aphasia treatment in a more tailored fashion for each patient. Future studies should further assess cognitive variables such as motivation, attention, learning, and memory, as well as comprehension in order to better predict treatment-induced recovery. Moreover, potential research in this area should attempt to replicate the present framework by using alternative interventions. This includes not only other rehabilitation strategies but also different time frames. For example, language therapy used in this study was rather intense (3-hour training sessions, 5 days a week), albeit short in duration (2 weeks). It is possible that some patients require similarly intense training during longer periods of time in order for improvement to emerge on posttreatment testing. Importantly, we found no correlation between time poststroke and changes in naming before and after treatment, suggesting that the significant gains in naming that are reported here are probably not driven by spontaneous recovery. The systematic studying of treatment logistics (ie, treatment session length, duration, number/duration of breaks in between sessions, etc) is certainly warranted.
Structural Network Integrity as a Framework for Functional Recovery
Our group has recently demonstrated that cell loss and disruption of network connectivity is a pervasive phenomenon in patients with regional gliosis after an ischemic stroke. 18 White matter damage is well known to occur in experimental models of stroke, 45 but it has been underestimated in humans thus far due to methodological limitations in measuring structural connectivity. However, with improvement in spatial registration of brains with lesions 46 and in DTI fiber tracking techniques, it is now possible to measure radiographic correlates of Wallerian degeneration. 47 Furthermore, measurements of white matter pathways also allow for the comprehensive assessment of brain connectivity and construction of the brain connectome,48,49 which permits unprecedented analysis of whole-brain neural network architecture. The structural brain connectome is the scaffolding supporting brain function, and significant insight into this framework can be gained by investigating the relationship between behavioral measures and regional connectome variability, as well as graph theory network analyses of connectome conformation.17,50 Using brain connectome methodology, optimized to take into account anatomical deformations associated with poststroke gliosis, 46 our group recently demonstrated a direct association between disrupted connectivity and chronic stage aphasia severity. 19 To our knowledge, that was the first study assessing the relationship between connectome integrity and aphasia severity in a group study. Therefore, it is possible that structural connectivity also supports brain plasticity necessary for naming improvement. In fact, the extent of structural damage to gray matter components of the language network has been seen as the main limiting factor for treatment-induced naming improvement.13,30 Consequently, there is vast potential for correlating behavioral performance with the “strength” (ie, density of streamlines in the case of structural connectomes) between different brain areas. This connectivity approach complements conventional voxel-based lesion-symptom mapping (VLSM) 39 by elucidating some of the neurobiological mechanisms underlying the relationship between brain damage and loss of function. Specifically, this connectivity approach provides information about how white matter loss contributes to impairment and recovery. The connectivity approach should not be seen as antagonist to VLSM, but rather as its complement, because it permits the quantification of one of the elements of the lesion, that is, cortical disconnection. Since cortical disconnection can exist in regions that are located outside the lesion, the disconnection approach may hypothetically provide more information compared with VLSM in specific cases where disconnection may be more prominent than the necrotic lesion.
Moreover, the connectivity approach demonstrated here may also provide information about the underlying neural network responsible for orchestrating recovery. For instance, if a cluster of voxels in 2 different brain areas were shown to correlate with naming performance on VLSM, a correlation between behavior and the link between these 2 regions favors their concerted action in contributing to performance. On the contrary, if the link between areas relevant to behavior revealed by VLSM does not correlate with performance, their contribution is ought to be independent from one another. As a third possibility, if connectivity analyses demonstrate a significant contribution of areas not lesioned by the stroke, this would be suggestive of connectional diaschisis, the process by which structural and functional connectivity change between brain regions as a result of a distant lesion. 51 Our approach thus calls for follow-up studies directly comparing the ability of VLSM- and connectivity-derived measures to predict prognosis. Specifically, it prompts for the assessment as to whether there are subgroups of patients where one approach may be particularly suited at explaining recovery.
In this study, we demonstrated that preservation of global neural network architecture and, in particular, preservation of the left temporal neocortex influence onto the brain networks (ie, left temporal BC) are both important determinants of recovery. This was achieved by conducting analyses both at the whole-brain (NSW) and at the nodal level (BC). For the latter, we selected specific ROIs that would be relevant for naming. This was based on extensive research showing the involvement of left frontal, temporal, and parietal cortices as the core structures involved in naming. This is of course an arbitrary decision, as naming is a very complex function, relying on an extensive network that may extend beyond the structures selected for our analyses. In fact, recent evidence has proposed a potential role of dominant-homologous right hemisphere regions in aphasia recovery. 20 The 12 left hemisphere Brodmann areas we included, however, cover the vast majority of the language network and are in agreement with current neurobiological models of speech processing and naming.38,52,53 Taken together, therefore, our results show that preserved cortical connectivity is key for enabling anomia recovery, likely because it supports functional modulation during the reestablishment of language networks motivated by therapy.
Left Temporal Cortex Influence on the Global Network
We observed that the average BC from specific left temporal cortical ROIs was associated with naming recovery. Interestingly, nodal BC is a measure of participation, or influence, 50 of a given node on the global network, measuring how important that node is to the structural integrity of the network; BC denotes the ratio of the shortest paths within the network that involve a specific node. 37 That is, it quantifies how many routes between 2 other nodes require passing through a specific node. Therefore, BC is not necessarily high in highly connected nodes. For example, in a network with many random links, a node with a high number of connections may still remain isolated from most of the remaining connections. Conversely, if a node is central to the network and participates in many connections, it will have a high BC. In biological networks, a balance between efficiency and economy of connections likely determines the overall organization of the network. Areas with great influence (or high BC) are strategically placed to facilitate information transfer and reduce the randomness of connections.
To facilitate the visualization of network organization, Figure 6 provides examples of network configurations from 2 patients with relatively similar necrotic lesion burden, but with different network metrics and different levels of improvement in the PNT. The complete overview of network configurations from all patients is demonstrated in Supplemental Figure 2.

The network configurations of 2 representative patients are demonstrated as 2-dimensional circular diagrams, where each anatomical region of interest (ROI) is symbolized by a circle corresponding to the network node (with numbers indicating the Brodmann area [BA]), and the links between nodes are represented as curved lines. The nodes are organized in accordance with their anatomical distribution, that is, left ROIs are placed in the left side of the figure and ROIs are grouped based on lobes. The node size is proportional to the size of the anatomical ROI (log2[size of the ROI]). To improve visualization, the network was made sparser by preserving only the links above the 95% link-weight percentile. The color of the node represents the percentage of the ROI that was damaged by the stroke (in accordance with the color bar). Note that subject 8 has higher LT BC because there are more connections traversing the LT area (ie, to travel between 2 other ROIs, the path “bounces of” LT). FR, frontal regions; MT, medial and inferior temporal regions; LT, lateral temporal regions; PR, parietal regions; OC, occipital regions.
The Role of the Left Temporal Cortex in Language Recovery
Naming impairment may occur as a result of lesions affecting several aspects of the cortical language network. Hillis and colleagues demonstrated that several regions are essential for distinct processes underlying naming. 54 Specifically, anomia may arise from combined dysfunction involving the left anterior, inferior, and posterior middle/superior temporal cortex, posterior inferior frontal, and inferior parietal cortices.54,55 Dronkers and colleagues demonstrated that damage to the left midposterior middle temporal gyrus (MTG) prevents the retrieval of names associated with objects (ie, lexical-semantic retrieval). 56 These findings complement observations from Schwartz et al, demonstrating that semantic errors during naming are associated with damage to the left anterior and mid MTG,57,58 while phonological errors are associated with lack of integrity of the left supramarginal gyrus and inferior frontal cortex. 59 Furthermore, Hillis et al have demonstrated in several studies that BA 37 (ie, located in the posterior portion of the mediobasal and lateral surfaces of the temporal lobe, including the fusiform gyrus and inferior temporal gyrus) is crucial for amodal lexical retrieval.60-62 Our results are in accordance with these observations. They are also in accordance with a recent model of speech proposed by Hickok and Poeppel. 52
This model suggests that speech is associated with 2 distinct major processing pathways, a ventral stream and a dorsal stream. The white matter pathways connecting the posterior and anterior aspects of the temporal lobes likely correspond to the ventral stream, which processes sound-to-meaning relationships and supports lexical retrieval and semantic associations. Conversely, the dorsal stream, which links perisylvian parietotemporal, inferior frontal gyrus, anterior insula, and premotor cortex, supports articulation processing. 52
Our results demonstrate that temporal lobe regions directly involved in the ventral stream pathways (ie, within the middle inferior and superior temporal gyri) are associated with naming recovery. Moreover, recovery is not solely associated with the preservation of connectivity of these regions, but also with how much these regions influence the brain network. This is a novel observation that directly supports and complements the dual stream model, particularly as it relates to the role of temporal cortex in speech processing.
Interestingly, our findings regarding global network preservation suggest that connectivity integrity extending beyond the left temporal lobe is also important for aphasia recovery. This observation is in accordance with a recent study by Forkel et al, 20 which demonstrated that integrity of the left and right arcuate fasciculi play a role in spontaneous aphasia recovery. Together, their findings and ours support the notion that global network configurations, are determinants of aphasia recovery. These results prompt for future studies exploring the ability of structural connectivity measures in predicting treatment-induced recovery in an independent sample of patients, for example, undergoing different intervention regimes. It will also be important to assess structural connectivity in predicting score improvement among healthy controls. Usually, however, naming tasks such as the one used in this study tend to show a ceiling effect among control participants.
Global Network Integrity and Language Recovery
The relationship between global network integrity (NSW) and recovery can shed additional light on the importance of the preservation of nonrandom structural organization of neural networks for recovery. Interestingly, the mechanisms underlying preservation of global network architecture may be related to complex interaction between prestroke “brain health” and the disruptive effect of large stroke lesions on global neural network. Pre- and poststroke “brain-health” may be an important predictor of recovery and this may be a promising line for research in neurological rehabilitation. Clearly, further research is needed to substantiate this notion.
Limitations and Future Studies
The conclusions derived from the present study favor the use of graph theory measures derived from the structural connectivity in assessing prognosis of aphasia due to chronic stroke. Our findings corroborate the exciting and evolving field of brain connectomics, and we acknowledge a number of limitations that should be taken into account for future studies of this nature. First, we recruited a heterogeneous population of chronic stroke patients, pooling together different aphasic syndromes. Second, we focused on just one aspect of aphasia, namely, anomia. While these can be caveats from certain perspectives, they also constitute important advantages for the purposes of the present study. This is because by focusing on one particular behavioral measure across a diverse population of patients, we were able to obtain a broad range in language performance (ie, naming ability). The variability in such measure was used in network-based lesion-symptom mapping to identify which links were relevant in the recovery of anomia. Importantly, this is not an all-or-none, binomial categorization. Anomia and its recovery are graded phenomena, and so is the connectivity strength between 2 brain structures. Thus, by recruiting a heterogeneous sample, we demonstrate that brain connectivity can be useful in predicting changes in performance, despite the clinical label used to categorize each patient. This is especially important because performance variability is also observed across patients within the same clinical group. Nonetheless, future studies with larger sample sizes should attempt to determine whether connectivity measures are better suited for predicting language recovery in different subgroups of patients, as well as using various measures of language performance.
In addition, all patients underwent identical treatment, despite their individual baseline performance or type of aphasia. This constitutes a limitation inherent to rehabilitation studies seeking to provide a consistent and homogeneous intervention that reduces confounders. There is converging evidence demonstrating that the same patients reliably show improvements when undergoing different treatment modalities.63,64 Accordingly, future research could, for instance, employ more than one intervention and, by using a counterbalanced design, determine whether connectivity measures are able to superiorly predict recovery when patients undergo an intervention that better suits their individual needs. In the same context, we acknowledge that there are several ways to determine what exactly constitutes an improvement in naming. Here, we used a similar approach to Lasar and et al, 65 but alternative definitions can be considered.
Despite very recent advances offering an alternative approach, 66 structural connectivity relies on the segmentation of neural tissue into predefined regions of interest. This tethers analyses to a priori anatomical assumptions, depending on factors such as number of ROIs and volume. Such segmentation based on previously defined atlases can be convenient, however, as a large number of studies use these predetermined regions to describe structural and functional changes in the brain. In our study, we decided to employ Brodmann areas because of their widespread use both in research and clinical settings among patient populations with language deficits. We acknowledge, nonetheless, that the ROIs of this atlas are large, usually encompassing different brain regions that serve distinct functions (eg, anterior and posterior aspects of the temporal gyri). For this reason, future studies should attempt to replicate this approach using other atlases that, by segmenting the brain into a larger number of ROIs, ensure that each one of these areas is confined to a more restricted volume of neural tissue.
Another important issue to take into consideration is the time of assessment posttherapy. Here, we evaluated patients within a week after a 2-week intense intervention. Longitudinal follow-up of these patients is crucial, and more studies are needed to determine whether medium- and long-term outcomes can also be predicted using brain connectivity measures.
Conclusions
Preserved global neural network architecture, with maintenance of the temporal lobe influence in the configuration of the brain neural networks, is crucial to support naming recovery in individuals with chronic aphasia. These are neurobiological markers that explain more than three quarters of the variability in treatment-related naming improvement and suggest that preserved network integrity relates to potential brain plasticity the drives treatment induced improvements in naming in aphasia.
These results have the potential to encourage the translation of brain connectivity studies into real-life clinical contexts, by providing health practitioners with a tool to assess the potential for recovery. In doing so, physicians, speech pathologists, and caregivers can plan for relevant non-pharmacological interventions to subjects who are more likely to benefit from treatment, even in the chronic stages of aphasia. Better strategies can also be adopted for subjects with aphasia and their relatives to readjust their expectations about recovery and outcome, which has the potential to positively contribute to mental health and quality of life.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by The National Institute on Deafness and Other Communication Disorders (NIDCD) (Grant Numbers DC008355 and DC014021).
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
