(C) PLOS One This story was originally published by PLOS One and is unaltered. . . . . . . . . . . Object knowledge representation in the human visual cortex requires a connection with the language system [1] ['Bo Liu', 'State Key Laboratory Of Cognitive Neuroscience', 'Learning', 'Idg Mcgovern Institute For Brain Research', 'Beijing Normal University', 'Beijing', 'Department Of Radiology', 'First Hospital Of Shanxi Medical University', 'Taiyuan', 'Shanxi'] Date: 2025-05 Abstract How world knowledge is stored in the human brain is a central question in cognitive neuroscience. Object knowledge effects have been commonly observed in higher-order sensory association cortices, with the role of language being highly debated. Using object color as a test case, we investigated whether communication with the language system plays a necessary role in knowledge neural representation in the visual cortex and corresponding behaviors, combining diffusion imaging (measuring white-matter structural integrity), functional MRI (fMRI; measuring functional neural representation of knowledge), and neuropsychological assessments (measuring behavioral integrity) in a group of patients who suffered from stroke (N = 33; 18 with left-hemisphere lesions, 11 with right-hemisphere lesions, and 4 with bilateral lesions). The structural integrity loss of the white-matter connection between the anterior temporal language region and the ventral visual cortex had a significant effect on the neural representation strength of object color knowledge in the ventral visual cortex and on object color knowledge behavior across modalities. These contributions could not be explained by the potential effects of the early visual perception pathway or potential confounding brain or cognitive variables. Our experiments reveal the contribution of the vision-language connection in the ventral occipitotemporal cortex (VOTC) object knowledge neural representation and object knowledge behaviors, highlighting the significance of the language-sensory system interface. Citation: Liu B, Wang X, Wang X, Li Y, Han Y, Lu J, et al. (2025) Object knowledge representation in the human visual cortex requires a connection with the language system. PLoS Biol 23(5): e3003161. https://doi.org/10.1371/journal.pbio.3003161 Academic editor: Huan Luo, Peking University, CHINA Received: October 19, 2024; Accepted: April 12, 2025; Published: May 20, 2025 Copyright: © 2025 Liu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: All the data necessary to reproduce Figs 2–4, S3–S5 and Tables 1, S4–S6 are provided in S1 Data. The reconstructed WM tract masks, the VOTC masks, and the MATLAB codes for RSA and probabilistic tractography analysis are publicly available at Zenodo (https://zenodo.org/records/15070254; DOI: 10.5281/zenodo.15070254). The neuroimaging data supporting the current study are not shared openly due to ethical constraints, but are available from the Ethics Committee of the First Hospital of Shanxi Medical University (phone: +86 351 4639242) for researchers who meet the criteria for access to confidential data. Funding: This work was supported by the STI2030-Major Project (https://en.most.gov.cn/, 2021ZD0204100 (2021ZD0204104) to Y.B.), the National Natural Science Foundation of China (https://www.nsfc.gov.cn/english/site_1/index.html, 31925020 and 82021004 to Y.B., 32171052 to X.S.W., 32071050 to X.Y.W.), and the Fundamental Research Funds for the Central Universities (http://en.moe.gov.cn/, to Y.B.). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist. Abbreviations: ATL, anterior temporal lobe; CSF, , cerebrospinal fluid; ; EPI, , echo-planar imaging; ; FA, fractional anisotropy; FDR, false discovery rate; FLAIR, fluid-attenuated inversion recovery; FLIRT, FMRIB’s linear image registration tool; fMRI, functional MRI; FMRIB, Functional Magnetic Resonance Imaging of the Brain; FNIRT, FMRIB’s non-linear image registration tool; FOV, field of view; GLM, general linear model; HARDI, high angular resolution diffusion imaging; IFOF, , inferior fronto-occipital fasciculus; ; ILF, , inferior longitudinal fasciculus ; JHU, Johns Hopkins University; LAG, left angular gyrus; LdlATL, left dorsolateral anterior temporal lobe; LIFGorb, left inferior frontal gyrus, orbital part; LMFG, left middle frontal gyrus; LpMTG, left posterior middle temporal gyrus; MMSE, Mini-Mental State Examination; MNI, Montreal Neurological Institute; OP, occipital pole; RDM, representational dissimilarity matrix; RF, radio frequency; ROI, , region of interest; ; RSA, representation similarity analysis; SD, standard deviation; SPACE, Sampling Perfection with Application optimized Contrast using different flip angle Evolution; TE, echo time; TR, repetition time; VOTC, ventral occipitotemporal cortex; WM, white-matter Introduction The human brain stores rich object knowledge that supports our interactions with objects through various sensorimotor modalities. How such memory is stored in the brain is one of the core questions in cognitive neuroscience. Decades of brain imaging and neuropsychological studies have converged on revealing that knowledge about specific object attributes is grounded in higher-order sensory association cortices, derived from perceptual experiences (see reviews in [1–3]), which are presumably then bound together to form higher-order object concepts [4,5]. One line of the foundational evidence comes from object color knowledge. Retrieving knowledge about object color, when presented with a name or a grayscale object image (e.g., red for rose) or when processing objects with salient diagnostic colors relative to other objects (e.g., rose versus chair), activates regions in the ventral occipitotemporal cortex (VOTC) that also engage in perceiving colors [6–11]. The classical cases with object color knowledge deficits [12–14] had lesions spanning the ventral temporal cortex, including, although not confined to, the occipitotemporal cortex (i.e., extending to the ventral medial temporal pole regions). Thus, the current consensus is that the VOTC is the neural correlate supporting color knowledge, related to visual color perception. Whether object perception functionality in the VOTC, and by extension, the perception-derived knowledge representation, is driven only by visual/sensory properties or modulated by various high-order computations, especially language, remains vigorously debated [e.g., 15–19]. For example, while Wang and colleagues [19] showed that large visual computation models that incorporate language supervision yield additional explanatory power of the human VOTC functionality over pure vision models, other researchers have shown that such advantages might be due to differences other than language involvements across models ([18]; see also discussion in [19]), and pure vision models can well explain many seemingly high order visual functionalities in the VOTC [e.g., 15–17]. A key challenge is that in a typical system, information structures derived from various levels of perception and language are difficult to disentangle (cat and dog are related both in visual and language spaces, and in tactile and auditory spaces). In the case of object color, a line of recent studies with congenitally blind individuals (i.e., resulting in full deprivation of color sensory experiences) has identified object color knowledge representation in the dorsal anterior temporal lobe (ATL), in addition to the VOTC. This representation is thus likely derived from language experience ([20–23]; see [24] for review); the dorsal ATL cluster is intrinsically and functionally connected to the language system and with the VOTC in healthy controls [22]. This anterior non-sensory color knowledge representation thus constitutes a strong candidate that may support normal object color knowledge in individuals with congenital blindness, or patients with cortical blindness from bilateral V1 lesions or with cerebral achromatopsia [22,25,26]. These findings raise two questions: (1) Is the VOTC’s object knowledge representation affected by such connections with the anterior, language-derived, representation? (2) Do such connections with the language system contribute to supporting object color knowledge behavior? In other words, is the neural representation in the VOTC sufficient to “know” that a banana shape should be yellow in color without connection with anterior language-derived regions? Suggestive evidence that VOTC functionality alone may be insufficient to support object color knowledge and thus may depend on its communication with higher-order regions comes from studies of patients with semantic dementia. Rogers and colleagues [27] reported that semantic dementia patients, who had diffuse atrophy predominantly in the ATL, exhibited severe deficits in object color knowledge, even in non-verbal visual tasks, despite having relatively preserved sensory association cortices and primary perceptual processing (assumed to reflect intact VOTC functionality). A recent study further revealed that object color knowledge deficits in semantic dementia patients were related to decreased structural connections between the left fusiform gyrus (spanning from the anterior to posterior ventral temporal lobe) and the early visual cortex (calcarine) [28]. Given that brain atrophy in semantic dementia patients may affect the functionality of surrounding or remote areas [29–31], it is unclear what aspects of this connection in relation to the VOTC account for the object color knowledge deficits in these patients. Here, we combine structural, diffusion, and functional MRI (fMRI) and behavioral measures to examine the role of visual‒language connections in object knowledge representation in a group of stroke patients. Specifically, we tested whether lesions affecting the connection between the VOTC and language systems affect the local neuro-functional strength of object color representation in the VOTC (Analysis series 1) and whether the connection contributes to supporting object color knowledge behavior (Analysis series 2). We therefore deliberately included stroke patients suffering from lesions outside the VOTC to minimize the confounding effects of local damage to the VOTC on its functionality. Both the patient (N = 33) and healthy control groups (N = 35) underwent diffusion imaging for white-matter (WM) structural integrity measurement, task-based fMRI experiments for neural representation measurement, and object color knowledge neuropsychological tests for behavioral measurement. Whether and how VOTC-language WM tract integrity affects neural representation in the VOTC-color-knowledge mask and object color behaviors were examined in detail. The analysis scheme is shown in Fig 1. PPT PowerPoint slide PNG larger image TIFF original image Download: Fig 1. Flowchart of the imaging and behavioral data analysis and the lesion distribution map of the 33 patients. The data processing workflow was as follows: For task fMRI data: ① The neural RDM (a) was constructed by correlating the activity patterns of each pair of fruits and vegetables (estimated from the task fMRI data) within a sphere (radius = 6 mm) centered on each voxel, via Pearson’s correlation distance. The object color RDM of each subject (b) was obtained from pairwise object color similarity ratings. Partial Spearman’s rank correlation was then computed between the neural and object color RDMs, controlling for low-level visual control, shape, and general semantic RDMs, to construct the object color RSA maps of healthy controls (c) and patients (g). ② In the functionally defined VOTC-color-perception mask (d), the top 300 selected voxels (i.e., voxels with the highest Fisher-Z transformed r values) in individual participants’ object color RSA maps were binarized to construct the individual-level mask (e). ③ The individual masks were overlaid across those of all healthy controls and a group-level threshold of 0.25 was applied to obtain the group-level functional mask (i.e., VOTC-color-knowledge mask, f). ④ The mean Fisher-Z transformed r value of each patient’s object color RSA map (g) within the VOTC-color-knowledge mask (f) was calculated to quantify the neural representation of object color knowledge in the VOTC. For HARDI data: ⑤ Probabilistic tractography was implemented using the estimated HARDI data (h) between each pair of ROIs (e.g., VOTC-LdlATL, i) in the native space. The resulting tracking map was normalized and standardized using the maximum voxel intensity of each image and then binarized at 0.1 to construct the individual-level mask (j). ⑥ A group-level threshold of 0.5 across all healthy controls within the explicit WM mask was applied to obtain the group-level WM mask (k). ⑦ The mean value of each patient’s FA map (l) within the group-level WM mask (k) was calculated to quantify the WM integrity of this fiber bundle. Correlation analysis: ⑧ The correlations between the WM integrity (mean FA value) of each WM tract and the neural representation strength of object color knowledge in the VOTC-color-knowledge mask across patients were calculated to identify the tracts relevant to the VOTC knowledge representation. ⑨ Object color behaviors were measured as the composite score of the grayscale picture-color word matching task (verbal) and the object color true/false judgment task (non-verbal). ⑩ ⑪ The correlations were computed between the WM integrity, the neural representation of object color knowledge, and the object color composite score, respectively, to identify the neural correlates associated with object color behavior. The upper-right panel shows the lesion distribution map of the 33 patients. The n value of each voxel denotes the number of patients with a lesion. Brain imaging results were visualized using BrainNet Viewer (version 1.7; https://www.nitrc.org/projects/bnv/; RRID: SCR_009446), MRIcron (version 1.0.20190902; https://www.nitrc.org/projects/mricron; RRID: SCR_002403), or MRIcroGL (version 1.2.20210317; https://www.nitrc.org/projects/mricrogl). Abbreviations: RDM, representational dissimilarity matrix; VOTC, ventral occipitotemporal cortex; RSA, representation similarity analysis; ROI, region of interest; HARDI, high angular resolution diffusion imaging; L, left; dlATL, dorsolateral anterior temporal lobe; WM, white-matter; FA, fractional anisotropy. https://doi.org/10.1371/journal.pbio.3003161.g001 Discussion Combining HARDI data, task fMRI data on object color processing, and object color behavioral data from 33 stroke patients, we made two main discoveries. First, across patients, the structural integrity of the VOTC-LdlATL connection exerts a robust effect on the neural representation strength of object color knowledge in the VOTC. Second, VOTC-LdlATL tract integrity accounts for object color knowledge behavioral task performance beyond the potential effects of VOTC representations. This effect is above and beyond its contribution to related cognitive processes (i.e., object recognition and color patch recognition). These two lines of results remained robust when we controlled for a wide range of potential confounding variables, including the broader effects of lesions, patient post-onset duration, and patients with specific stroke etiologies, and were not attributable to potential associative impairment in the early visual perception pathway. Our main results nicely accommodate the previous imaging and patient lesion evidence and offer critical clarifications. As described in the Introduction, object color representation has been observed in healthy individuals in both the VOTC, presumably derived from sensory experience, and the dorsal ATL, presumably derived from non-sensory, language experience (see review in Bi [24]). Whether the commonly assumed sensory-derived representations function independently and sufficiently to support object knowledge remains unresolved. Object color knowledge deficits have been reported in stroke cases with ventral occipitotemporal lesions (spanning anteriorly) and in semantic dementia patients with ATL lesions, indicating that neither the posterior representation nor the anterior representations might be sufficient to support intact object color knowledge behavior [12,14,27,28,44]. The results of a case study showing impaired color name comprehension following a stroke that functionally disconnected color-biased visual regions from the dorsal and lateral parts of the left ATL, are suggestive of the potential role of the VOTC-LdlATL connection [45]. Inspection of lesion locations across these stroke cases and semantic dementia patients revealed that their lesions may affect the VOTC-LdlATL WM tract to varying degrees. However, conclusions about implicated neural structures are difficult to draw, given the intrinsic difficulty of case studies. By showing that, indeed, this connection disruption leads to VOTC functional representation reduction and has unique explanatory power for object color knowledge behavioral deterioration, our findings well explain the batch of patient findings and answer questions posed by the dual knowledge systems [24]. Note that while the VOTC functional integrity was significantly modulated by its connection with the LdlATL, its correlation with object color behavior was not significant. This null effect should be interpreted with caution. As we deliberately selected stroke patients whose lesions largely spared the VOTC, this patient cohort may not be optimal for revealing the effects of the VOTC. Another possibility is that behavior is the readout of many types of representations, and the performance for judging object color knowledge tested here may (also) pull from other types of representations stored elsewhere, i.e., the integration of such VOTC representations and those in the LdlATL. The key point is that, in explaining behavior, the effects of the VOTC-LdlATL connection are beyond the potential effects of VOTC functionality. Can our findings regarding the VOTC-LdlATL connection simply be because the fMRI task involves mapping to color words, i.e., mapping object color knowledge in the VOTC to color words represented in the ATL? We reason that this is unlikely. The fMRI task was designed to require minimal verbal retrieval—only the color name “red”. The instruction was to judge whether the object depicted in the grayscale image was usually red. All patients had no difficulty recognizing the color red, as shown in the color word matching with the color patch task (see S4 Fig). By the null hypothesis that VOTC representation is sufficient, the patients could accomplish the task by activating the object color knowledge representation in the VOTC, comparing it with the red color representation there, and making a correct response. Furthermore, in the object color behavioral assessments, we also included a non-verbal task (true/false color judgment), where mapping to words is not intrinsic to the task. Yet we observed that the disruption of the VOTC-LdlATL connection not only affected VOTC functionality but also contributed to color knowledge behavior independently. What is the functionality of the observed VOTC-LdlATL connection? The connection observed here was obtained by WM tracking between color-perception clusters in the VOTC and the subregion of the ATL parcellated from the language-activation mask [36] in healthy controls. This particular tracking to the language system was motivated by the findings about the language-derived representations here [22]. Our findings on the connection between the VOTC and the LdlATL, and not the other language parcels in the frontal and posterior temporal regions, converge with previous findings on identifying disembodied and language-experience-related knowledge representations in the dorsal ATL [21–23], highlighting the central role of this particular neural structure, more so than the other language regions, in interfacing with the sensory-derived representations to support knowledge. Rich evidence has shown that different subdivisions of the ATL serve different aspects of semantic representation, with the ventral-medial region being more related to the multimodal binding of object knowledge and the dorsolateral aspect being more related to language-derived abstract representations ([21,22,43]; see Bi [24], Lambon Ralph and colleagues [5] for review). Indeed, our ATL-subregion analyses confirmed that the connection between the VOTC and the dorsolateral ATL is indeed effective, although the negative effects of the ventral ATL should not be overinterpreted due to suboptimal lesion coverage and imaging signals (Figs 4B and S5). To better situate the target tract in the larger WM pathways, we overlaid the target connection with the major WM tracts in the Johns Hopkins University (JHU) WM template and found that it partly comprises the inferior longitudinal fasciculus (ILF; 50.6% overlapping voxels), the inferior fronto-occipital fasciculus (IFOF; 24.1% overlapping voxels), and, at minimum, the forceps major (2.9%). These results are in accordance with results of previous studies showing that these large WM bundles are involved in supporting object knowledge: the IFOF has been consistently shown to be involved in higher-order object semantic processing across multiple modalities [46–48], and the ILF tends to be implicated in non-verbal object processing tasks [49,50], although its role in more abstract semantic relationships or mapping with object words is less clear [51–53]. It has been shown that such long-range WM connections are one candidate that supports high-dimensional semantic space by incorporating cortical representations that process different specific dimensions of information [54]. Consistent with this broad semantic function of these WM connections, we also observed that the integrity of the VOTC-LdlATL connection contributes significantly to object recognition performance. Importantly, the VOTC-LdlATL connection integrity explained object color knowledge beyond such object recognition effects. Considering the previous findings about different types of object color knowledge representations (non-sensory language derived versus color sensory derived; see review in Bi [24]), one likely mechanism of this connection is to form an integral object color knowledge representation by incorporating different types of color representations (sensory-derived in the VOTC and language-derived in the LATL). It offers a solution to the question posed by the dual knowledge representation framework [24,55,56], i.e., how two different forms of representations are integrated. This observation extends the previous notion of representing concepts by connecting different attributes [3], to representing object knowledge by connecting different representation forms of the same attribute. Another, non-mutually exclusive possibility, is that language representations play an online modulatory role in categorizing object properties, especially along perceptual dimensions such as color, where other differences need to be abstracted away [57,58]. The results also, more broadly, shed light on a classical debate of the effects of language on cognition [e.g., 19,59,60]. By providing positive evidence for the effects of the vision-language connection on object (non-verbal) knowledge representations, the results speak against a strong modular view of the two systems [e.g., 60]. These findings, by highlighting the role of perception-language system connections in knowledge loss, motivate rehabilitation strategies that target such components, such as cognitive trainings that involve vision-language-alignment and/or brain stimulation targeting the underlying connectional functionality. As in all lesion studies, the power of our analyses is constrained by the lesion coverage of the patient group. Here, we specifically avoided lesions affecting the VOTC in the interest of examining the contribution of non-posterior connections. The resulting group tended to exhibit a typical lesion distribution of middle cerebral artery strokes, with comparable impairment in connections of interest, and making inferences about other connections should be refrained. To conclude, addressing the question of object knowledge representation and taking object color as a test case, the neural representation in the sensory-derived VOTC is modulated by the disruption of its connection with the LdlATL in the language network, which also explains object color knowledge behavior better than the VOTC neural representation does. These results underscore the critical role of the VOTC-LdlATL connection in representing object knowledge, implicating its specific involvement in bridging the language-derived representation and the sensory-derived knowledge representation, which in turn modulates VOTC functionality. These findings highlight the intricate interaction between vision and language systems in the human brain and the importance of studying brain connections in supporting knowledge representation [54]. Materials and methods Participants Thirty-three stroke patients and 35 healthy controls participated in the fMRI and behavioral experiments (see Table 1 for the demographic data of the two groups). All participants were native Chinese speakers who received payments and provided written informed consent. The study has been conducted according to the principles expressed in the Declaration of Helsinki, and was approved by the Ethics Committee of First Hospital of Shanxi Medical University (Approval No. 2021-K035). Thirty-three stroke patients (25 males), without major lesions affecting the VOTC, were recruited from the First Hospital of Shanxi Medical University. The inclusion criteria were as follows: age 20–65 years; right-handed before stroke onset [61]; normal vision or corrected vision; at least 3 months after stroke; first symptomatic stroke (ischemic or intraparenchymal hemorrhagic); no contraindications for MRI; lesion location involving the cortex and/or subcortical WM and not the VOTC; no other neurological or psychiatric diseases; ability to perform simple cognitive tasks and understand instructions; and intact color vision on the color vision examination plates test [62, similar to the Ishihara test]. The Chinese version of the Mini-Mental State Examination [MMSE, 63] was used as a measure of general cognitive state. The detailed demographic information and lesion locations of each patient are presented in S1 Table. Thirty-five healthy controls (21 males) matched to stroke patients by age, gender, and years of education participated in the study, with the following inclusion criteria: right-handed; normal vision or corrected vision; no contraindications for MRI; no history of neurological or psychiatric disorders; and intact color perception on the color vision examination plates test. Note that one patient (ID: 018) had a premorbid distance vision problem (not improved with correction) and thus participated in the auditory version of the task fMRI experiment, where the stimuli were auditory words (names of fruits/vegetables) with otherwise the same presentation parameters as the visual version. In addition, while our patient inclusion criteria required first-time strokes, two patients’ (ID: 003, 007) scans revealed potential old lesions. To ensure maximum lesion coverage, these three patients were included in the study and were further excluded in the control analyses. Procedures Task-fMRI experiment. An object color judgment task-fMRI experiment was conducted to obtain object color neural representations in stroke patients and healthy controls. To minimize the potential confounding effect of task difficulty on the neural representations and measure meaningful neural activities, we tailored the task so that the patients could perform this task in the scanner with ease. We deliberately selected items that were familiar to the patients, adopted a simple yes/no judgment task, and performed the following procedure before the fMRI scanning to ensure that the participants recognized the grayscale stimuli of fruits and vegetables. The participants were first asked to name each stimulus displayed at the center of the computer screen. If a patient had a naming problem, he/she was given a word-to-picture matching session, where an auditory fruit/vegetable word was presented, and the patient needed to select the correct image out of the whole stimuli set. For the objects that the participants failed to name or match, they were provided with correct answers to ensure object recognition. Grayscale images of 16 common fruits and vegetables (400 × 400 pixels; visual angle of 6.22 × 6.22 degrees) were used as stimuli. The images were matched for mean luminance separately for the foregrounds and backgrounds using the SHINE toolbox [64]. The number of pixels was comparable across the pictures (range: 29,889−31,866 pixels, mean = 30895.25, SD = 674.15). The visual stimuli were displayed using a Visual & Audio Stimulation System (Sinorad SA-9939) on a screen (89 cm width, 50 cm height; 171 cm from the participant’s eyes). The participants viewed the display through an angled mirror (45°) attached to the headcoil. During scanning, the participants viewed the stimuli and verified whether the typical color of the skin of the fruit/vegetable was red (Fig 2A). The participants pressed a button with their right index finger to give a “yes” response and pressed another button with their right middle finger for a “no” response; stroke patients with lesions affecting the left hemisphere responded with their left hand. All participants performed 4 runs (400 s per run) of slow-event-related task fMRI scanning, except for one patient (ID: 002) who completed three runs. Each run consisted of 32 trials, consisting of a 1-s-long stimulus period and an 11-s-long fixation period. Each stimulus was presented twice within each run. The order of the 32 trials was pseudorandomized, with the restriction that no two consecutive trials were identical. Each run began with a 12 s fixation period and ended with a 4 s fixation period. Independent pseudo-randomizations were created for each run for each participant. The experimental procedure was conducted using the E-prime 2 (Psychology Software Tools, Pittsburgh, PA, USA). Individual behavioral RDMs of color, shape, and semantics were obtained for the 16 stimuli used in the fMRI task. The participants were asked to rate the similarity of these objects in terms of color, shape, and semantic relatedness (1: most similar, 7: most dissimilar), in a pairwise fashion (i.e., 120 pairs), respectively. Pairs were presented in pseudorandom order, with no time limit, and independent pseudo-randomizations were created for each feature rating. Before each rating, the participants practiced with 10 pairs of novel but familiar fruits and vegetables with similar and dissimilar examples to understand the rating instructions. Individual RDMs were used for RSA (see below). Because we deliberately ensured that the items used in the fMRI experiment were familiar to the patients, the obtained group-averaged RDMs of the controls and patients were highly correlated (color: rho = 0.92, p < 0.001; shape: rho = 0.95, p < 0.001; semantics: rho = 0.94, p < 0.001). Two items (chili and watermelon) were excluded from the analysis, because even healthy control participants generated highly variable colors for chili (red, green, yellow; all correct), and were uncertain about the color of watermelon referred to the skin (green) or the inside (red). Fourteen items were included in all the analyses below. Neuropsychological tests. The participants were administered a series of neuropsychological tests outside the scanner to evaluate their object color behavioral task performances. Object color knowledge tasks: Object color behavior was assessed on two tasks that varied in the modalities of input: a grayscale picture-color word matching task (verbal) and an object color true/false judgment (non-verbal) task. In the grayscale picture-color word matching task, participants were asked to judge whether the color word they heard matched the typical color of the vegetable/fruit in the grayscale picture. In the object color true/false judgment, a picture of a vegetable/fruit with typical or atypical colors appeared in the center of the screen, and participants were asked to judge whether the vegetable/fruit was typically colored. In both tasks, the same 30 familiar vegetables and fruits that appeared in different exemplars were used as stimuli. Each item appeared twice, once with the typical color or color word and once with the atypical one. The participants were asked to press the left mouse button for “yes” and the right mouse button for “no”. The presentation order of the items was pseudorandom, and the same items did not appear consecutively, nor did the same responses appear consecutively three times. Each trial had a response deadline of 60 s. Responses were scored 1 if correct and 0 if wrong or if no response was given within the time deadline. Peripheral control tasks of object and color recognition: To assess the (relative) neural specificity related to object color behaviors, we also carried out a word-picture matching task and a color patch matching task, which involved the two cognitive processes most relevant to object color processing: object recognition and color recognition. For the word-picture matching task, the stimuli included 20 familiar vegetables and fruits. In each trial, a chromatic picture of the vegetable/fruit appeared in the center of the screen, and an auditory vegetable/fruit word was played simultaneously. The participants were asked to judge whether the words matched the picture. Each item was presented twice, once in a matching condition and once in a non-matching condition. The stimulus presentation settings of this task were consistent with those of the object color knowledge tasks. In the color patch matching task, participants were asked to match the color word they heard to the 11 visually presented color patches (orange, green, gray, black, yellow, pink, purple, white, blue, brown, and red). The order of the color patches in each trial was pseudorandom, and each trial was presented without a time limit. Tasks with grayscale pictures always preceded those with colored pictures to mitigate task-related priming. Each participant was tested in multiple sessions in a quiet room. Each session lasted for 2 h with pauses for rest upon request. Stimulus presentation was controlled by E-prime 1 (Psychology Software Tools, Pittsburgh, PA, USA). One healthy control from the task fMRI participant cohort did not participate in the object color knowledge and control tasks; three healthy controls and one stroke patient did not participate in the color patch matching task; the patient with the visual problem (ID: 018) completed the auditory version of pairwise similarity ratings and the visual versions of other tasks (with corrected vision), but did not participate in the color patch matching task. In addition to the task fMRI cohort, another group of 17 healthy controls took part in the object color knowledge tasks and the word-picture matching task; one of them further finished the color patch matching task. Image acquisition All functional and structural MRI data were collected using a 3T Siemens Magnetom Skyra scanner with a 32-channel radio frequency (RF) head coil, at the Department of Magnetic Resonance Imaging, First Hospital of Shanxi Medical University. Foam padding was placed inside the RF head coil to minimize head movement. The scans included task-state fMRI, HARDI, high-resolution 3D T1-weighted images, 3D T2-weighted images, and 3D fluid-attenuated inversion recovery (FLAIR) T2-weighted images. A multi-band echo-planar imaging (EPI) sequence was used for the fMRI and HARDI images. The acquisition parameters for fMRI were as follows: slice planes scanned along the rectal gyrus, phase encoding direction from posterior to anterior, repetition time (TR) = 2,000 ms, echo time (TE) = 30 ms, flip angle = 90°, axial scanning, field of view (FOV) = 180 mm × 180 mm, matrix size = 72 × 72, slice thickness = 2.5 mm, voxel size = 2.5 × 2.5 × 2.5 mm3, and multi-band factor = 2. The acquisition parameters for HARDI were as follows: diffusion gradient b values set to 0, 1,000, and 2,000 s/mm2, b0 repeated 10 times, 64 directions each for b1000 and b2000, TR = 3,000 ms, TE = 100 ms, flip angle = 90°, axial scanning, FOV = 224 mm × 224 mm, matrix size = 112 × 112, slice thickness = 2 mm, voxel size = 2 × 2 × 2 mm3, and multi-band factor = 2. To maximize image resolution and increase the accuracy of lesion delineation, the T1-weighted, T2-weighted, and FLAIR T2-weighted images of this project all underwent 3D thin-layer scanning. The 3D T1-weighted images were acquired using a magnetization-prepared rapid gradient-echo sequence, with the following acquisition parameters: TR = 2,530 ms, TE = 2.88 ms, flip angle = 7°, sagittal scanning, FOV = 224 mm × 256 mm, matrix size = 224 × 256, interpolated to 448 × 512, slice thickness = 1 mm, and voxel size = 0.5 × 0.5 × 1 mm3. For the 3D T2-weighted images, a Sampling Perfection with Application optimized Contrast using different flip angle Evolution (SPACE) sequence was utilized, and the acquisition parameters were as follows: TR = 3,200 ms, TE = 408 ms, flip angle = 120°, sagittal scanning, FOV = 230 mm × 230 mm, matrix size = 256 × 256, slice thickness = 0.9 mm, and voxel size = 0.9 × 0.9 × 0.9 mm3. A SPACE Inversion Recovery sequence was used for the 3D FLAIR T2-weighted images, and the following acquisition parameters were used: TR = 5,000 ms, TE = 394 ms, flip angle = 120°, sagittal scanning, FOV = 250 mm × 250 mm, matrix size = 256 × 256, interpolated to 512 × 512, slice thickness = 1 mm, and voxel size = 0.5 × 0.5 × 1 mm3. Imaging data preprocessing fMRI data. Analyses were conducted using SPM12 (RRID:SCR_007037; https://www.fil.ion.ucl.ac.uk/spm12/), based on MATLAB (Mathworks; RRID: SCR_001622). For each task run, the first six images were discarded to mitigate the effects of scanner instability at the onset of scanning. The remaining images underwent the following procedures: (1) temporal slice timing correction to adjust for the differences in acquisition time between slices within a single whole-brain scan; (2) realignment of each image frame to the first frame of the sequence to reduce the impact of head motion; (3) registration of EPI images to 3D T1 images; (4) normalization of 3D T1 images to MNI space using the unified segmentation algorithm [65]; (5) spatial normalization of EPI images using the normalization parameters obtained in the previous step, and resampling to a voxel size of 2 × 2 × 2 mm3. Two healthy controls were excluded from the fMRI analysis: one self-reported an uncomfortable feeling in the head during scanning, and the other exhibited excessive head motion during the scans (> 2.5 mm/2.5°). Four healthy controls and four stroke patients showed excessive head motion in 1 or 2 runs (healthy controls: > 2.5 mm/2.5°, stroke patients: > 3 mm/3°), and we analyzed the remaining runs of fMRI data for these participants. HARDI data. Preprocessing was performed using the Functional Magnetic Resonance Imaging of the Brain (FMRIB) Analysis Group Software Library (FSL, version 6.07; Oxford University, United Kingdom), which included the following procedures: (1) Eddycorrect, which corrected for eddy current distortions and head movement; (2) BET, which was used for skull removal; (3) DTIFIT, which was used for building diffusion tensor models and calculating the FA maps. We subsequently registered the FA images to the T1 images in the native space using the FMRIB’s linear image registration tool (FLIRT), and then used the T1-to-MNI transformation warp image to register the FA images to the Montreal Neurological Institute (MNI) space, with a target voxel size of 2 × 2 × 2 mm3. The warp image was obtained from normalizing 3D T1 images to MNI space with FLIRT for linear rigid affine transformations, and FMRIB’s non-linear image registration tool (FNIRT) for non-linear transformations. The derived transformation parameters were then inverted and used to warp the region of interest (ROI) from the MNI space to the native diffusion space using nearest-neighbor interpolation. Lesion mask drawing and normalization. Lesions were manually drawn by a well-trained radiology resident (Y.H.) and reviewed by an experienced radiologist (Y.L.), using the ITK-SNAP software (www.itksnap.org). The specific delineation and spatial normalization protocol included: (1) coregistering and resampling T2-weighted and FLAIR T2-weighted images to the native space of 3D T1-weighted images; (2) visually referring to the corresponding T1-weighted, and T2-weighted images, drawing on the axial plane of FLAIR T2-weighted images slice by slice, with gliosis included within the lesion outline; (3) employing the normalization parameters (obtained in step 4 of the fMRI data’s preprocessing) to register the lesion images from native to MNI space, with a voxel size of 1 × 1 × 1 mm3; and (4) applying a 3 mm smoothing kernel to the lesion images in standard space to obtain the patient-specific brain injury map (with 1 indicating injury and 0 indicating no injury). All processes and resulting normalized images were further inspected by an experienced clinical investigator (B.L.) for quality control. Data analysis Defining VOTC masks of interest in healthy controls. GLM: For each participant, general linear modeling was conducted using preprocessed fMRI data. The general linear model (GLM) for each run included 16 regressors, each corresponding to one of the fruits or vegetables. The button presses and the six head motion parameters from each run were incorporated as regressors of no interest. Each regressor was convolved with the canonical SPM hemodynamic response function and high-pass filtered with a cutoff set at 128 s. The resulting t value maps for each vegetable/fruit compared with the baseline were obtained for subsequent multivariate pattern analysis. Object color RSA: Searchlight RSA analyses were carried out with individual participants [66,67]. The neural RDM was constructed using Pearson’s correlation distance of the activation patterns (t value maps) within a sphere (radius = 6 mm) centered on each GM voxel across 14 fruits and vegetables pairwise. The 14 × 14 neural RDM was then correlated with individuals’ object color similarity RDM using partial Spearman’s rank correlation, while controlling for low-level visual RDM (i.e., gist RDM), shape RDM, and general semantic RDM, to identify regions specifically encoding object color knowledge. The gist RDM measures the dissimilarity of low-dimensional visual representations across images, on the basis of spectral and coarsely localized information [68]. The resulting individual r maps were Fisher-Z transformed. Probability maps: To account for inter-participant variability, the VOTC functional mask representing object color knowledge was localized using probability maps across healthy controls [45,69]. For each healthy control, the obtained RSA map specific to object color knowledge was masked with a functionally defined, VOTC-color-perception mask (defined in Wang and colleagues [22], by contrasting chromatic stimuli to grayscale stimuli in a color perceptual localizer in 14 healthy controls). In each masked image, the top 300 selected voxels (i.e., voxels with the highest Fisher-Z transformed r values) were binarized to construct the individual-level mask. Probability maps were created by summing the individual-level masks and dividing them by the number of healthy controls. Voxels that had a probability greater than 0.25 were retained as a group-level functional mask (i.e., VOTC-color-knowledge mask). To confirm the stability of the results produced by these arbitrary thresholds, we conducted validation analyses, using functional masks with various combinations of individual-level thresholds (top n voxels: 200, 250, 300, 350, 400, 450, and 500) and group-level thresholds (probability across healthy controls: 0.15, 0.2, 0.25, 0.3, and 0.35). Defining WM masks of interest in healthy controls. In this phase of the study, we mapped the WM connections between the left language and visual perceptual systems, as well as the control connections, in healthy controls through quantitative analysis utilizing probabilistic tractography. Probabilistic tractography, an algorithm that investigates the probability distributions of fiber orientations at each voxel in the brain, offers the advantages of explicitly representing uncertainty [70] and better delineating crossing fibers in the brain [71,72]. Particularly for datasets with multiple b-values (e.g., the HARDI data we collected), the application of probabilistic tractography with a gamma distribution may prove advantageous for combining multiple shells, thereby increasing the accuracy of fiber orientation estimates in WM and areas adjacent to the cortex [73]. ROI selection: The ROIs included a set of left language-specific regions, defined by Fedorenko and colleagues (contrasting intact sentences to nonword lists in 220 healthy subjects), the functionally defined VOTC-color-perception mask (defined in Wang and colleagues [22], by contrasting chromatic stimuli to grayscale stimuli in a color perceptual localizer in 14 healthy subjects), and the early visual cortex (the OP (48#) in the Harvard-Oxford Atlas, probability > 0.2). The right hemisphere homologs of the left-hemisphere language regions, as defined by Fedorenko and colleagues, were also included to investigate the lateralization of the observed effects. Given that different subdivisions of the ATL may serve different aspects of semantic representation, especially along the dorso-ventral axis ([21,22,43]; see Bi [24], Lambon Ralph and colleagues [5] for review), we examined the specificity of the LdlATL by further including different ATL subregions as ROIs. These subregions were parcellated either by functional activation connectivity profiling [43] or anatomical profiling (i.e., Harvard-Oxford Atlas (probability > 0.2)). For the functional activation connectivity profiling, the ATL was parcellated into four subregions (superior dorsal, inferior dorsal, lateral, and ventromedial) based on coactivation clustering, and the two dorsal subregions in that study were combined here into a single dorsal subregion. For the anatomical profiling, according to the Harvard-Oxford Atlas (probability > 0.2), the ATL subregions included the temporal pole (8#), the anterior superior temporal gyrus (9#), the anterior middle temporal gyrus (11#), the anterior inferior temporal gyrus (14#), the anterior parahippocampal gyrus (34#), and the anterior temporal fusiform cortex (37#). These ROIs were transformed to the diffusion native space for each participant via the inverse of the linear and nonlinear transformations as previously described. Probabilistic tractography: Probabilistic tractography was carried out using healthy controls’ HARDI data with FMRIB’s Diffusion Toolbox (FDT v5.0, http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FDT) in FSL. We employed the ROI-to-ROI approach, in which tractography was implemented between each pair of ROIs for the following tracts: (1) the main tracts of interest—the connection between the VOTC and the left language regions; (2) the tract involved in lower-level visual perception—the connection between the VOTC and the early visual cortex; (3) the right hemisphere homologs of the main tracts of interest—the connections between the VOTC and the right homologs of the left language regions; and (4) the connections between the VOTC and the ATL subregions. The detailed processes are as follows. Tractography was estimated using FSL’s BEDPOSTX with default parameters to model WM fiber orientations and crossing fibers, where the diffusion coefficient was modeled using a gamma distribution. Fiber tracking was initiated in both directions (from seed to target and vice versa) in the diffusion space, and 5,000 streamlines were drawn from each voxel in the ROI using Probtrackx 2.0. A binarized cerebrospinal fluid (CSF) mask (generated by setting the probability threshold > 0.75 for the CSF probability mask in SPM12) was set as an exclusion mask for the analyses. The resulting images containing the output connectivity distribution were normalized to MNI space and standardized using the maximum voxel intensity of each image, resulting in a standardized image with voxel values ranging from 0 to 1. We then applied an individual-level threshold for reducing false-positive fiber tracks combined with a group-level threshold for the fiber tracking success rate [74,75]. At the individual level, the standardized path images were thresholded at 0.1 to remove voxels with a low connectivity probability. The resulting maps were then binarized for each individual and summed across individuals. Next, at the group level, fiber projections that existed in more than 50% of the individuals within the explicit WM mask (probability > 0.4 in the WM probability map in SPM12) were retained for subsequent analysis. The threshold was selected arbitrarily, which is a common practice in tracking studies for reporting group results in general anatomical assessments. We further conducted validation tests, using WM connection masks with various combinations of individual-level thresholds (connectivity probability: 0.05, 0.1, 0.15, 0.2) and group-level thresholds (probability: 0.35, 0.4, 0.45, 0.5, 0.55, 0.6), with or without the explicit WM mask. To verify the tracts we reconstructed, we overlaid each tract onto the 20 major tracts in the JHU WM tractography atlas (thresholded at 25%; http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/Atlases) in FSL. The overlap percentage was computed by dividing the number of common voxels by the total number of voxels in the template tract. VOTC object color neural representation in patients: Effects of the WM tract integrity and specificity tests. To investigate whether the lesions to the tracts connecting the left language regions and the VOTC resulted in impaired VOTC object color neural representation, we conducted partial Spearman’s rank correlation analyses between the WM integrity of each tract and the strength of the VOTC object color neural representation across patients, with each patient’s total lesion volume as a covariate. For WM integrity, we computed the mean FA value within each tract. For the strength of the VOTC object color neural representation, we extracted the mean values from the object color RSA results within the VOTC color-knowledge mask (defined in healthy controls above). A false discovery rate (FDR) with a threshold of q < 0.05 was adopted to correct for comparisons of multiple tracts. Given the potential overlaps among the masks of the language-VOTC tracts, we carried out a specificity test by excluding voxels that overlapped with other tracts and considering only the mean FAs of the tract-specific voxels and then computed Spearman’s rank correlations. To investigate whether the observed effects of interest were influenced by potential alternative neuroanatomical and/or cognitive variables, in addition to controlling for total lesion volume, we conducted a series of control analyses: (1) GM lesions. To examine whether the WM effects are inherited from relevant GM lesions, we performed partial correlations, including the lesion percentage in the GM regions at both ends of this connection as covariates. The GM areas were identical to the seed regions used in probabilistic tractography, and lesion percentages in these GM regions were calculated as the number of voxels with lesions divided by the total number of voxels in the region. (2) Visual perception pathway. Considering that the VOTC functionality may be affected by lesions in the WM tracts connecting to the early visual cortex, we thus further added the mean FA values of these WM connections as covariates. (3) Patient etiology. Data from three patients (ID: 003, 007, 018) warrant scrutiny, as described in the subsection “Participants”. To eliminate the potential influence of these special patients, we performed correlation analyses after their exclusion. We further included the time post-onset (month) as a covariate along with total lesion volume given the substantial variation in this variable (ranging from 3 to 252 months). Object color behavior in patients: Underlying neural correlates and specificity tests. Behavioral data analysis: For the two object color knowledge tasks (grayscale picture-color word matching, object color true/false judgment) and peripheral control tasks (word-picture matching, color patch matching), accuracy scores were analyzed. Given the variations along the demographic variables (i.e., age, gender, and years of education), we employed a statistical method to adjust for these variables. The Bayesian test for a Deficit allowing for Covariates [41] was used to compare a patient’s performance with that of controls and to obtain the point estimates of effect sizes for the difference between a patient and controls (controlling for the influence of three demographic variables) as standardized scores. Cases with one-tailed p < 0.05 were considered to have (potential) deficits. Correlation analyses: In this section, we examined how the VOTC and related connections support object color behavior. Hence, we performed partial Spearman’s rank correlation analyses between the strength of the VOTC object color neural representation or the mean FA value of the observed tracts and the neuropsychological test performance across patients, controlling for total lesion volume. As object color behaviors were assessed in two tasks that varied in terms of input modalities, we computed a composite score by averaging the effect size estimate of accuracy from the two tasks (controlling for the three demographic variables and considering the variation in healthy controls). When correlated with neural measures, the composite score was considered and reported in the main text. The results of each individual task are shown in S4–S6 Tables. Differences between the correlations for the two tasks were compared using Hotelling’s t test (the FZT computator, http://psych.unl.edu/psycrs/statpage/regression.html), with results reported in S6 Table. Specificity tests: To investigate whether the observed results for WM tract integrity in object color behavior were influenced by various types of potentially confounding variables, we also conducted a series of control analyses by including lesions to GM regions, visual perception pathway integrity, and patient etiology as covariates. Importantly, to assess whether the observed tracts of interest were specifically associated with object color knowledge relative to other cognitive processes, we included the standardized scores of the word-picture matching task and the color patch matching task as covariates in addition to total lesion volume. Acknowledgments We thank Weiwei Men for help with imaging methodology, Yuxing Fang for help with tractography methodology, and all members of the Department of Radiology at the First Hospital of Shanxi Medical University for data collection. We are also grateful to all research participants. [END] --- [1] Url: https://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.3003161 Published and (C) by PLOS One Content appears here under this condition or license: Creative Commons - Attribution BY 4.0. via Magical.Fish Gopher News Feeds: gopher://magical.fish/1/feeds/news/plosone/