La remodelación del glioblastoma de los circuitos neuronales humanos disminuye la supervivencia

Nov 23, 2023

Nature volume 617, pages 599–607 (2023)Cite this article

36k Accesses

1 Citations

767 Altmetric

Metrics details

Gliomas synaptically integrate into neural circuits1,2. Previous research has demonstrated bidirectional interactions between neurons and glioma cells, with neuronal activity driving glioma growth1,2,3,4 and gliomas increasing neuronal excitability2,5,6,7,8. Here we sought to determine how glioma-induced neuronal changes influence neural circuits underlying cognition and whether these interactions influence patient survival. Using intracranial brain recordings during lexical retrieval language tasks in awake humans together with site-specific tumour tissue biopsies and cell biology experiments, we find that gliomas remodel functional neural circuitry such that task-relevant neural responses activate tumour-infiltrated cortex well beyond the cortical regions that are normally recruited in the healthy brain. Site-directed biopsies from regions within the tumour that exhibit high functional connectivity between the tumour and the rest of the brain are enriched for a glioblastoma subpopulation that exhibits a distinct synaptogenic and neuronotrophic phenotype. Tumour cells from functionally connected regions secrete the synaptogenic factor thrombospondin-1, which contributes to the differential neuron–glioma interactions observed in functionally connected tumour regions compared with tumour regions with less functional connectivity. Pharmacological inhibition of thrombospondin-1 using the FDA-approved drug gabapentin decreases glioblastoma proliferation. The degree of functional connectivity between glioblastoma and the normal brain negatively affects both patient survival and performance in language tasks. These data demonstrate that high-grade gliomas functionally remodel neural circuits in the human brain, which both promotes tumour progression and impairs cognition.

Malignant brain tumours such as glioblastomas exist within the context of complex neural circuitry. Neuronal activity promotes glioma growth through both paracrine signalling (neuroligin-3 and brain-derived neurotrophic factor (BDNF)) and AMPAR (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor)-mediated excitatory electrochemical synapses1,2,3,4.

Likewise, glioblastomas influence neurons, inducing neuronal hyperexcitability through the secretion of non-synaptic glutamate and synaptogenic factors5,6 and reducing inhibitory interneurons7. Beyond preclinical models, we previously demonstrated in awake, resting patients that glioblastoma-infiltrated cortex exhibits increased neuronal excitability2. The mechanisms by which glioblastomas maintain the ability to engage with neuronal circuitry and alter cortical function remain incompletely understood9. Deciphering the processes by which gliomas remodel neural circuits may uncover therapeutic vulnerabilities for these lethal brain cancers. To address these gaps in knowledge, we performed intraoperative electrophysiology while patients engaged in language tasks: we analysed local field potentials in glioblastoma-infiltrated cortex during speech initiation, determined the decodability of neural responses and revealed biological drivers of synaptic enrichment in glioblastoma cells (Extended Data Fig. 1 and Supplementary Tables 1 and 2).

Glioblastomas and other high-grade gliomas interact with neural elements, resulting in cellular- and network-level changes10,11,12,13. While neurons within glioblastoma-infiltrated brain are hyperexcitable at rest, the extent of task-specific neuronal hyperexcitability and the ability to extract neural features from glioma-infiltrated cortex remain unclear. To examine cognitive task-specific neuronal activity from glioblastoma-infiltrated cortex, we selected a cohort of adult patients with cortically projecting tumours in the lateral prefrontal cortex (LPFC; Extended Data Fig. 2a). Electrocorticography (ECoG) electrodes were placed over tumour-infiltrated and normal-appearing cortex. ECoG signals filtered between 70–110 Hz were used for analysis of high-gamma band range power (HGp), which is strongly related to local neuronal population spikes14,15 and is increased by cortical hyperexcitability16. Spectral data demonstrated clear separation of frequencies across tumour and non-tumour electrodes (Fig. 1a and Extended Data Fig. 3a).

a, In participants with dominant hemisphere glioblastomas, we applied subdural ECoG over the posterior lateral frontal cortex during an audiovisual speech initiation task to assess circuit dynamics. Spectral data show the expected pattern of HGp increasing above 50 Hz in addition to clear separation of frequencies across tumour and non-tumour electrodes. b, The posterior lateral frontal cortex (outlined area) time series of HGp within tumour-infiltrated cortex between −600 ms and speech onset (0 ms). c, High-gamma (HG) recordings from averaged electrodes within each patient, while averaging the effect across the sampled region of cortex for an individual showing greater HG power within electrodes overlying tumour-infiltrated cortex (n = 14 patients, F1,21 = 25.562, P = 0.00005). Data are median (centre dot), first to third quartiles (bars) and the minimum and maximum points (whiskers). d, Electrodes were compared between non-tumour and tumour-infiltrated regions; the false-discovery rate (FDR)-corrected HGp demonstrates task-relevant hyperexcitability (P = 0.016). Data are mean ± s.e.m. e, Event-related spectral perturbations (ERSPs) during a naming task for low-frequency words (low freq., left column) and high-frequency words (high freq., middle column) in normal-appearing non-tumour regions (top row) and glioma-infiltrated (bottom row) cortex. Signals from high-frequency word trials were able to be decoded above chance in normal-appearing cortex (mean accuracy = 0.56, P = 0.000089) but not in glioma-infiltrated cortex (mean classifier accuracy = 0.49, P = 0.72) using a regularized logistic regression classifier with leave-one-participant-out cross-validation (right column). Data are mean ± 95% confidence interval. AU, arbitrary units. For b–d, statistical analysis was performed using two-sided linear mixed-effects models (b–d), and corrections for multiple comparisons were performed using FDR adjustment (b and d). For e, P values were determined using two-tailed Student's t-tests with Bonferroni multiple-comparison correction for the number of timepoints (left, ERSP) and one-sided Z-tests (right, classification accuracy).

Source data

ECoG was recorded from the dominant hemisphere LPFC during auditory and visual picture naming as an illustrative example of a well-defined cognitive neuronal circuit17. While patients were awake and speaking, HGp was recorded for single-electrode (Extended Data Fig. 3b) and group-level analysis. HGp data from control and non-tumour conditions demonstrate the expected neural time course of speech motor planning within the LPFC (Extended Data Figs. 2b and 3c,d), consistent with previously established models of speech initiation demonstrated in non-human primates and humans18,19. We next performed the same analysis focused only on electrode arrays recording from tumour-infiltrated cortex. Countering the theory that glioblastoma–synaptic integration may result in physiologically disorganized neural responses, we found task-relevant neural activity within the entire region of the tumour-infiltrated cortex, including cortical regions that are not typically implicated in speech production (Fig. 1b)—a notable finding that indicates tumour-induced functional remodelling of language circuitry. Similarly, we found that, across WHO grade 2–4 glioma subtypes, task-specific neuronal responses for speech initiation are maintained within the LPFC (Extended Data Fig. 2c,d). These findings demonstrate that neuronal activity within tumour-affected cortex is physiologically organized, including neuronal activity elicited by speech tasks in regions that are outside of regions that are typically involved in speech production.

In light of this finding of preserved task-evoked neural responses from tumour-infiltrated cortex, we next examined whether the magnitude of neural responses may differ in tumour-affected cortical language areas. We therefore pair-matched tumour-infiltrated and normal-appearing cortex (Extended Data Fig. 3d,e), demonstrating increased HGp during speech production in glioblastoma-infiltrated cortex, consistent with hyperexcitability (Fig. 1c,d).

Neural computations for speech vary by condition. Vocalization of infrequently used (low frequency) words, for example, requires a more intricate coordination of articulatory elements than that of commonly used (high frequency) words20,21. We therefore determined the decodability of neuronal signals from normal-appearing and glioblastoma-infiltrated cortex using a logistic regression classifier to distinguish between low-frequency and high-frequency word trial conditions (Fig. 1e). We implemented identical training and leave-one-participant-out cross-validation paradigms for both conditions. Normal-appearing cortex produced above-chance decoding between low- and high-frequency word trials. By contrast, glioblastoma-infiltrated cortex did not decode word trials above chance. These data further demonstrate that glioblastoma infiltration into the human cortex maintains task-specific neuronal responses, including neuronal hyperexcitability, yet tumour-affected cortex loses the ability to decode complex word conditions.

Having demonstrated that gliomas remodel neuronal circuits, we next examined whether specific molecularly defined glioma cellular subpopulations influence functional integration of the tumour into neural circuitry. Glioblastoma cells are heterogeneous22,23,24 and previous findings indicate that oligodendrocyte-precursor-cell-like subpopulations are enriched for synaptic gene expression2, whereas astrocyte-like subpopulations secrete synaptogenic factors8,25. Thus, functionally connected regions may vary within tumours and differences in functional connectivity between tumour regions may be due at least in part to varying subpopulations of glioma cells. With the goal of sampling functionally connected regions within gliomas, we measured neuronal oscillations within glioma-infiltrated brain using magnetoencephalography (MEG) and sampled primary patient glioblastoma tissues with varying functional connectivity during surgical tumour resection26,27,28. The connectivity of an individual voxel was derived by the mean imaginary coherence between the index voxel and the rest of the brain29,30,31. Intratumoural functional connectivity correlated with neuronal activity within tumour-infiltrated cortex and high functional connectivity (HFC) voxels were identified both within tumour regions that were contrast-enhancing or T2/FLAIR hyperintense on magnetic resonance imaging (MRI; Extended Data Fig. 4).

To investigate the differences between functionally connected, HFC and non-functionally connected low functional connectivity (LFC) tumour regions, we performed bulk and single-cell RNA sequencing (RNA-seq) analyses. Bulk RNA-seq transcriptomic analysis revealed upregulation in HFC tumour regions of genes that are involved in the assembly of neural circuits, including axon pathfinding genes (NTNG1, also known as netrin G1), synapse-associated genes (for example, SYNPO, also known as synaptopodin) and synaptogenic factors including a sevenfold upregulation of thrombospondin-1 (THBS1, encoding TSP-1). THBS1, which encodes a known synaptogenic factor that is secreted in the healthy brain by astrocytes32, was particularly interesting in the context of the observed remodelling of functional language circuitry described above (Extended Data Fig. 5a–c and Supplementary Table 3).

To further assess cellular subpopulation contribution to THBS1 expression, we performed single-cell sequencing analysis of biopsy samples from HFC and LFC tumour regions (Supplementary Table 4). Malignant tumour cells were inferred on the basis of the expression programs and detection of tumour-specific genetic alterations, including copy-number variants (Extended Data Fig. 6a–e). We found that 2.44% of all tumour cells expressed THBS1, and that HFC tumour cells expressed higher levels of THBS1 compared with LFC tumour cells (Fig. 2a and Extended Data Fig. 6g). Within LFC-region samples, THBS1 expression primarily derives from a non-tumour astrocyte population (Extended Data Fig. 6e–g). This suggests that, within low-connectivity intratumoural regions, astrocytes chiefly express THBS1, whereas, within HFC regions, high-grade glioma cells express THBS1 in addition to astrocytes and myeloid cells, which may promote the observed neural circuit remodelling (Extended Data Fig. 6h–j). Notably, myeloid cells, which include bone-marrow-derived macrophages, microglia, dendritic cells and neutrophils, chiefly comprise the glioblastoma tumour immune microenvironment (Extended Data Fig. 6d,e), and the microglial cell surface molecules CD36 and CD47 can function as TSP-1 receptors33,34. Although the role of TSP-1 in the tumour immune microenvironment is not yet clear, myeloid cell expression of TSP-1 suggests that multiple cell types in the tumour microenvironment of HFC regions may contribute to altered synaptic connectivity. Elevated expression of THBS1 within HFC regions was confirmed by protein-level analysis using HFC and LFC patient-derived glioblastoma biopsy tissues. Concordant with transcriptomic profiles, immunohistochemistry analysis demonstrated increased TSP-1 expression within HFC tissues (Extended Data Fig. 7a). Immunofluorescence and confocal microscopy analysis confirmed that malignant tumour cells express TSP-1 in HFC tissue (Fig. 2b). The fact that a subpopulation of malignant tumour cells in HFC regions produce TSP-1 suggests a differential potential of tumour cells in the HFC regions to promote synaptogenesis and thereby connectivity, consistent with the cancer biology principal that cellular subpopulations assume distinct roles within the heterogenous cancer ecosystem, which may be defined at least in part by functional connectivity measures.

a, Single-cell RNA-seq feature plot analysis of THBS1 in HFC (n = 6,666 cells, 3 participants) tissues; within HFC samples, THBS1 is primarily in glioblastoma cells (circled). b, TSP-1 immunofluorescence analysis of nestin-positive tumour cells in HFC and LFC tissues. n = 13 (HFC) and n = 11 (LFC) sections, 3 per group. P = 0.000073. Scale bar, 50 µm. The box plot shows the median (centre line), interquartile range (box limits) and minimum and maximum values (whiskers). c, The synapsin-1 puncta count in HFC and LFC glioblastoma tissue samples. n = 25 regions, 4 per group. P = 0.000014. Red, synapsin-1 (presynaptic puncta); white, neurofilament heavy and medium (neurons). Scale bar, 10 µm. Inset: magnified view of synapsin-1 puncta on neurons. Scale bar, 3 µm. d, PSD95 puncta count. n = 7 (HFC) and n = 9 (LFC) sections, 3 per group. P = 0.04. Red, PSD95 (postsynaptic puncta); white, neurofilament heavy and medium chains (NFH/M) (neurons). Scale bar, 10 µm. Inset: magnified view of PSD95 puncta on neurons. Scale bar, 3 µm. e, Representative confocal images showing synaptic punctum colocalization (yellow arrows). Red, synapsin-1; green, homer-1 (postsynaptic puncta); white, MAP2 (neurons); blue, 4′,6-diamidino-2-phenylindole dihydrochloride (DAPI). Scale bar, 10 µm. Quantification of the number of colocalized pre- and postsynaptic puncta (n = 13 (HFC) and n = 10 (LFC) regions, 2 per group; P = 0.005) and homer-1 puncta size in neuron–glioma co-culture (P = 0.000024). f, TSP-1 rescue of induced neuron (iN) organoids in co-culture with HFC and LFC cells for 6 h. Scale bar, 300 µm. Quantification of glioblastoma (GBM) cell integration measured on the basis of the fluorescence intensity of RFP-positive glioblastoma cells in the organoids. Significant differences between HFC and LFC groups (asterisks) and LFC and LFC + TSP-1 (hash) are indicated. n = 2 (HFC and LFC groups) and n = 1 (LFC + TSP-1 group). Scale bar, 300 µm. g, Representative MEA raster plots showing individual spikes (tick mark), bursts (cluster of spikes in blue) and synchronized network bursts (pink) after 48 h co-culture of cortical neurons (CN) with HFC and LFC cells (outlined in red and blue, respectively). Quantification of network burst frequency (Hz) (n = 2 (CN only), n = 3 (CN + HFC) and n = 4 (CN + LFC); P = 0.05) and network synchrony (area under normalized cross-correlation; n = 2 (CN only), n = 3 (CN + HFC) and n = 4 (CN + LFC); P = 0.0129 (CN versus CN + HFC); P = 0.0308 (CN + HFC versus CN + LFC)). Data are mean ± s.e.m. (b–g). P values were determined using two-tailed Student's t-tests (b–f) and one-way analysis of variance (ANOVA) with Tukey's post hoc test (g). *P < 0.05, **P < 0.01, ****P < 0.0001; NS, not significant.

Source data

Hypothesizing that this subpopulation of HFC glioma cells may promote synaptogenesis and consequent remodelling of connectivity as observed in glioma-associated language networks above, we next examined whether HFC-associated glioma cells promote structural synapse formation, similar to normal astrocytes35,36,37 and certain astrocyte-like glioblastoma cells8,25. We first analysed primary patient glioblastoma biopsies from HFC and LFC regions using immunohistochemistry and confocal microscopy. We found increased presynaptic neuronal puncta (synapsin-1; Fig. 2c) together with increased postsynaptic puncta density and cluster size on neurons (PSD95+neurofilament+), and synapsin–PSD95 puncta colocalization (Fig. 2d and Extended Data Fig. 7b) within HFC regions compared with LFC regions. Together, these data indicate increased synapse stability and/or synapse formation in high-connectivity regions of glioblastoma, supporting a role for TSP-1 in glioma-associated neural-circuit remodelling.

Primary patient-derived glioma cultures from HFC and LFC tumour regions were generated to perform further mechanistic experiments. We co-cultured high-grade glioma cells from HFC and LFC tumour regions with mouse hippocampal neurons to test the effects of TSP-1high- and TSP-1low-expressing primary patient-derived glioma cells on synaptic connectivity of neurons (Extended Data Fig. 7c,d). We then quantified the size of postsynaptic puncta (marked by the postsynaptic marker homer-1) and the number of colocalized pre- and postsynaptic puncta in HFC and LFC co-cultures with neurons. This demonstrated an increased number of colocalization points and postsynaptic homer-1+ punctum size in glioma and neuronal processes in HFC–neuron co-cultures compared with LFC–neuron co-cultures (Fig. 2e), additionally indicating a role for HFC glioma cells in synaptogenesis.

To further investigate the functional distinctions between malignant subpopulations isolated from HFC and LFC regions, we tested neuron–glioma interactions in a neuronal organoid model. We co-cultured HFC and LFC glioma cells with GFP-labelled human neuron organoids generated from an induced pluripotent stem (iPS) cell line integrated with a doxycycline-inducible human NGN2 transgene to drive neuronal differentiation38. Quantification of postsynaptic homer-1 in induced-neuron organoids revealed a relative increase in postsynaptic puncta density when co-cultured with HFC glioma cells compared with LFC glioma cells (Extended Data Fig. 7e). Live-cell imaging of neuronal organoids co-cultured with HFC and LFC glioma cells revealed that HFC glioma cultures exhibit prominent neuronal tropism and integrate extensively in the organoids, whereas LFC glioma cells displayed minimal integration with neuron organoids (Fig. 2f and Supplementary Videos 1 and 2). Notably, exogenous administration of TSP-1 to induced-neuron–LFC co-culture reversed this phenotype and promoted robust LFC glioma integration into the neuronal organoid (Fig. 2f and Supplementary Video 3), further implicating TSP-1 in neuron–glioma interactions. The electrophysical properties of TSP-1high-expressing cells in co-culture with neurons were analysed using multi-electrode array (MEA) electrophysiology. After co-culture for 48 h, the total number of network bursts (a measure of neuronal activity) from cortical neuron co-culture with TSP-1high-expressing HFC cells was increased relative to cortical neurons alone or under LFC co-culture conditions. Neurons in co-culture with HFC glioma cells also demonstrated increased network synchrony as measured by the area under normalized cross-correlation (the area under interelectrode cross-correlation normalized to the autocorrelations; Fig. 2g and Extended Data Fig. 7f).

Gliomas exhibit intratumoural heterogeneity with subpopulations of cancer cells assuming particular roles23,24. The human data presented above demonstrate localizational heterogeneity of functional integration in glioblastoma with normal brain circuity and suggest that, within intratumoural regions of HFC, a tumour subpopulation with synaptogenic properties exists. We next examined the structural synapses in TSP-1high-expressing HFC glioma cell-infiltrated mouse brain. RFP-labelled HFC or LFC glioma cells were stereotactically xenografted into the CA1 region of the mouse hippocampus2 (Fig. 3a). After a period of engraftment and growth, immuno-electron microscopy analysis identified neuron-to-neuron and neuron-to-glioma synapses2 (Fig. 3b). The total number of synapses (neuron-to-neuron and neuron-to-glioma combined) was significantly higher in HFC glioma xenografts than in LFC glioma xenografts (Fig. 3b and Extended Data Fig. 7g), further demonstrating a greater synaptogenic potential of glioma cells isolated from HFC patient tumour regions.

a, Representative micrograph showing RFP-labelled glioblastoma xenografted into the mouse hippocampus. Scale bar, 500 µm. b, Immuno-electron microscopy analysis of HFC or LFC cell xenografts. The asterisk denotes immuno-gold particle labelling of RFP. Postsynaptic density in RFP+ tumour cells (pseudocoloured red), synaptic cleft and clustered synaptic vesicles in apposing presynaptic neuron (pseudocoloured yellow) identify both neuron–glioma synapses in HFC-PDX (left) and neuron–neuron synapses in LFC-PDX (right). Quantification of the total number (neuron–neuron combined with neuron-glioma) of synapses per field of view in HFC/LFC xenografts. n = 4 mice per group. P = 0.0019. Scale bar, 1,000 nm. Data are median (centre line), with first and third quartiles (box limits) and the minimum and maximum points (whiskers). c, Representative immunohistochemistry images in glioblastoma tissues demonstrate increased Ki-67 protein expression in HFC samples. n = 13 (HFC) and n = 14 (LFC) regions, 4 per group. P = 0.04. Scale bar, 50 µm. d, Glioblastoma cells from HFC tissues show a marked increase in the proliferative index when co-cultured with mouse hippocampal neurons. n = 27 (HFC), n = 14 (HFC + neurons), n = 32 (LFC) and n = 14 regions (LFC + neurons), 3 per group. e, SEM images of HFC and LFC cells cultured in the presence or absence of neuronal conditioned medium (NCM) shows tumour microtubes (TMTs) that connect neighbouring cells through cytoplasmic extensions. Quantification of TMTs per cell. n = 10 (HFC), n = 8 (HFC + NCM), n = 17 (LFC) and n = 13 (LFC + NCM) regions, 2 per group. P = 0.0455. Scale bars, 20 µm (full fields) and 10 µm (magnified view). f, Quantification of the mean microtube length per spheroid. n = 11 (HFC), n = 16 (HFC + NCM), n = 5 (LFC) and n = 6 (LFC + NCM) spheroids, 1 per group. P = 0.000011. g, Representative SEM images showing TMTs and quantification of TMTs per cell from HFC shCtrl and HFC shTSP-1 conditions. n = 5 regions, 2 per group. P = 0.0012. Scale bar, 20 µm. h, Kaplan–Meier survival curves of mice bearing HFC or LFC xenografts. n = 4 (HFC) and n = 5 (LFC). P = 0.03. Data are mean ± s.e.m. (b–h). P values were determined using two-tailed Student's t-tests (b–g) and two-tailed log-rank analysis (h).*P < 0.05, **P < 0.01, ****P < 0.0001; NS, not significant.

Source data

Neurons promote glioma cell proliferation1,2,3,4 and we hypothesized that HFC cells may represent a cellular subpopulation within glioblastomas that are differentially regulated by neuronal factors. We found that primary patient biopsies from HFC and LFC regions demonstrated increased Ki-67 proliferative marker staining within HFC regions (Fig. 3c). To test whether HFC cells differentially proliferate in response to neuronal factors compared with LFC primary patient cultures, patient-derived HFC and LFC cells cultured alone or in co-culture with mouse hippocampal neurons were treated with 5-ethynyl-2′-deoxyuridine (EdU) overnight. HFC glioma cells exhibit a fivefold increase in proliferation when cultured with neurons. By contrast, the LFC glioma in vitro cell proliferation index (determined as the fraction of DAPI cells co-expressing EdU) is similar with and without hippocampal neurons in vitro (Fig. 3d and Extended Data Fig. 8). These results indicate that the ability of HFC cells to proliferate is contingent on the presence of neuronally secreted factors and that, in the absence of neuronal signals, they tend to acquire a dormant tumour phenotype.

Given the neuronal tropism exhibited by HFC glioma cells together with the concept that neural network integration requires invasion of brain parenchyma to reach and colocalize with neuronal elements, we next tested the effects of neuronal conditioned medium on invasion of HFC and LFC glioma cells using a spheroid invasion assay. LFC glioma cells demonstrated no differences in spheroid volume in the presence or absence of neuronal conditioned medium; however, HFC glioma cells exhibited an increased spheroid invasion area in response to neuronal conditioned medium. In addition to increased invasion area, HFC glioma cells extended long processes representing tumour microtubes in response to neuronal conditioned medium (Extended Data Fig. 9a). Tumour microtubes connect glioma cells in a gap-junction-coupled network1,39,40,41 through which neuronal-activity-induced currents are amplified2. Scanning electron microscopy (SEM) was performed on TSP-1high-expressing HFC and LFC cells in the presence or absence of neuronal conditioned medium, demonstrating robust cytoplasmic extensions connecting HFC cells (Fig. 3e). We also quantified the change in mean spheroid volume. We found that neuronal conditioned medium increased both invasion and microtube length in HFC but not LFC cultures (Fig. 3f and Extended Data Fig. 9a). Concordantly, the invasive marker MET was increased within HFC samples compared with LFC samples (Extended Data Fig. 9b,c). Primary patient-derived HFC cells were then transduced with a short hairpin RNA (shRNA) control or shRNA against THBS1 to knockdown TSP-1. Cell viability was confirmed using a live/dead assay with robust knockdown of our target protein (Extended Data Fig. 10a,b). Knockdown of THBS1 in HFC cells decreased the number of tumour microtubes relative to the control conditions (Fig. 3g), consistent with the known role for TSP-1 in tumour microtube formation39.

Glioblastoma cell invasion bears negative prognostic value. We therefore performed survival studies of mice that were orthotopically xenografted with patient-derived HFC or LFC glioma cells. Mice bearing HFC tumours exhibited greater tumour burden and shorter survival compared with LFC-tumour-xenografted mice (Fig. 3h and Extended Data Fig. 9d). Taken together, these results suggest that functionally connected intratumoural regions are enriched for a tumour cell population that is differentially responsive to neuronal signals and exhibits a proliferative, invasive and integrative phenotype in the neuronal microenvironment that negatively influences survival in a preclinical model.

We next investigated the effects of tumour-intrinsic functional connectivity on patient survival and cognition. First, we tested the hypothesis that gliomas exhibiting increased functional connectivity may be more aggressive, given the robust influence of neuronal activity on tumour progression2,3,4. We performed a human survival analysis of patients with newly diagnosed glioblastoma. After controlling for known correlates of survival (age, tumour volume, completion of chemotherapy and radiation, and extent of tumour resection)42, neural oscillations and functional connectivity were measured within tumour-infiltrated brain using MEG (Supplementary Tables 2 and 5). A Kaplan–Meier survival analysis illustrates an overall survival of 71 weeks for patients with functional connectivity compared with an overall survival of 123 weeks for participants without HFC voxels, illustrating a striking inverse relationship between survival and functional connectivity of the tumour (mean follow-up time, 50.5 months) (Extended Data Fig. 11a). To identify clinically relevant survival risk groups, we next used recursive partitioning survival analysis using the partDSA algorithm42,43,44. Within this analysis we controlled for important prognostic variables such as MGMT promoter methylation status. Overall survival risk was based on the interactive effects of all known prognostic variables (for example, age at diagnosis, sex, tumour location, chemotherapy, radiotherapy, the presence of functional connectivity within the tumour, pre- and post-operative tumour volume, and the extent of resection). The first division was based on known risk factors such as age and extent of tumour resection. Within this hierarchical model of partitioning, the degree of connectivity was identified as the next most important variable, which divided risk groups 2 and 3. Risk group 1 (black) had the worst outcomes and is the combination of patients older than 72 years or any age with less than 97% extent of tumour resection (subtotal resection). Risk group 3 (grey) had the best survival, and these are patients are younger than 62 years with over 97% extent of tumour resection and absence of functional connectivity in the tumour. Intermediate risk group 2 (red) revealed an interesting interaction between age and HFC. This group had two subsets: patients with over 97% resection of tumour and age younger than 72 years with intratumoural connectivity; and those between 62 and 72 years without functional integration (Fig. 4a,b). These results demonstrate the notable prognostic value of connectivity on survival. We next examined whether TSP-1, a secreted synaptogenic protein32,35, can be identified in the patient serum and whether circulating TSP-1 is correlated with functional connectivity. Circulating TSP-1 levels in the patient serum exhibited a notable positive correlation with intratumoural functional connectivity (Fig. 4c).

a,b, Schematic (a) and partDSA model (b) of overall survival in patients, incorporating the effects of glioblastoma intrinsic functional connectivity, therapeutic and clinical factors by recursive partitioning results into three risk groups. Risk group 1 (black) patients have the shortest survival, including a combination of (1) patients older than 72 and (2) patients younger than 72 with an extent of tumour resection (EOR) of less than 97%. Risk group 3 (grey) patients have the best survival, including patients who are younger than 62 with an extent of tumour resection of greater than 97% and no intratumoural connectivity. Intermediate risk group 2 (red) comprises a combination of patients with greater than 97% extent of resection and (1) an age of younger than 72 with tumour intrinsic connectivity or (2) patients between 62 and 72 years without connectivity. c, Linear regression statistics illustrate that serum TSP-1 is correlated with the extent of intratumoural functional connectivity. n = 56. P = 0.01. d, Representative MEA raster plots showing neuronal spikes (black tick marks), bursts (cluster of spikes in blue) and synchronized network bursts (pink) of neuron–HFC co-cultures (outlined in red) and 24–48 h exposure of neuron–HFC co-culture to (50 µM) GBP (outlined in orange). Quantification of the weighted mean firing rate (Hz) and network synchrony (area under normalized cross-correlation) from HFC and HFC + GBP glioma–neuron co-culture (weighted mean firing rate: n = 4 well, 2 per group; P = 0.04; area under normalized cross-correlation: n = 3 (HFC) and n = 4 (HFC + GBP); P = 0.007). e, Representative confocal images from neuron–HFC glioma co-culture showing a decrease in HFC cell proliferation after THBS1 knockdown using shRNA (n = 10 (HFC shCtrl) and n = 9 (HFC shTSP-1); P = 0.0068). Red, HNA (human nuclei); white, Ki-67. Scale bar, 30 µm. f, Representative confocal images from neuron–HFC glioma co-culture showing a decrease in HFC cell proliferation after gabapentin (32 µM) treatment for TSP-1 inhibition. n = 16 (HFC) and n = 15 (HFC + GBP), 2 per group. P = 0.0007. Red, HNA (human nuclei); white, Ki-67. Scale bar, 30 µm. g, Schematic for gabapentin treatment of HFC xenografted mice. i.p., intraperitoneal. h, Representative confocal images, and quantification demonstrating a decrease in the proliferation index (Ki-67+HNA+/HNA+) after gabapentin treatment in mice bearing HFC xenografts. n = 9 mice per group. P = 0.046. Red, HNA (human nuclei); white, Ki-67. Scale bar, 70 µm. Data are mean ± s.e.m. (d–f and h). P values were determined using two-sided linear regression analysis (c), and two-tailed (d–f) and one-tailed (h) Student's t-tests. *P < 0.05, **P < 0.01, ***P < 0.001; NS, not significant.

Source data

We hypothesized that, beyond survival, intratumoural functional connectivity may also influence cognition. We therefore performed visual picture and auditory naming testing in our cohort of patients with dominant hemisphere glioblastoma, given their correlation with aphasia in clinical populations45,46. Linear regression of the number of HFC voxels within tumours with language task performance demonstrated an inverse relationship between language cognitive performance and tumour functional connectivity (Extended Data Fig. 11b,c). Together, these findings suggest that functional integration of glioblastoma into neural circuits negatively influences cognition and survival.

Given the premise that TSP-1 serves as a regulator of neuronal activity-driven glioma growth, we sought to target TSP-1 therapeutically using gabapentin (GBP), which blocks the thrombospondin receptor α2δ-1 (ref. 47). In neuron–glioma co-cultures, individual spikes, bursts (cluster of spikes) and synchronized network bursts were reduced after 24–48 h exposure to GBP (Fig. 4d). Primary patient-derived HFC cells were transduced with shRNA-control or shRNA against THBS1 or treated with GBP. Pharmacological TSP-1 inhibition using GBP did not influence the proliferation of HFC cells grown alone in culture, verifying that there were no tumour cell-intrinsic effects of GBP (Extended Data Fig. 12). By contrast, genetic or pharmacological targeting of TSP-1 resulted in a marked decrease in proliferation of HFC glioma cells co-cultured with neurons (Fig. 4e,f). GBP administration to mice bearing HFC patient-derived xenografts (PDX) resulted in a marked decrease in glioma proliferation (Ki-67+HNA+/HNA+) in gabapentin-treated mice bearing HFC xenografts relative to vehicle-treated controls (Fig. 4g,h).

Integration of high-grade glioma into neural networks is manifested by bidirectional interactions whereby neuronal activity increases glioma growth1,2,3,4,48 and gliomas increase neuronal excitability5,6,7,8. To understand whether glioma–neuronal interactions influence neural circuit dynamics, we used short-range electrocorticography analysis of tumour-infiltrated cortex in humans to demonstrate language-task-specific activation as well as functional remodelling of language circuits. We further demonstrated that distinct intratumoural regions maintain functional connectivity through a subpopulation of TSP-1-expressing malignant cells (HFC glioma cells). This molecularly distinct glioma subpopulation is differentially responsive to neuronal signals, exhibiting a synaptogenic, proliferative, invasive and integrative profile. Previous research has demonstrated that neuronal activity promotes glioma proliferation through paracrine and synaptic signalling1,2,3,4, and we have now shown that patients with glioblastoma exhibiting functional connectivity between the tumour and the rest of the brain experience a shorter overall survival compared with patients without HFC. Pharmacological inhibition of TSP-1 decreases glioblastoma cell proliferation and network synchrony within the tumour microenvironment, highlighting a potential therapeutic strategy to be assessed in future clinical studies.

The neuronal microenvironment has emerged as a crucial regulator of glioma growth. Both paracrine signalling and connectivity remodelling may contribute to network-level changes in patients, affecting both cognition and survival. In patients, the role of neural network dynamics on survival and cognition remains poorly understood and how glioma–network interactions influence cognition remains unanswered. In fact, some studies using a heterogenous population of patients with both IDH wild type (WT) and mutant WHO grade III and IV gliomas have suggested that functional connectivity improves overall survival49,50,51; however, such previous research has been confounded by functional connectivity methods that are heavily influenced by the presence of tumour vascularity, limited spatial resolution and a heterogenous patient cohort. Nonetheless, the evidence in this study that glioblastomas remodel functional circuits and that functional connectivity negatively influences survival does not address direction of causality. It remains possible that glioma originating in functionally connected cortical regions are more strongly connected and may therefore exhibit greater network distribution, thereby encouraging distinct glioblastoma subpopulations with the ability to migrate52. A better understanding of the cross-talk between neurons and gliomas as well as how functional integration affects clinical outcomes may open the door to a range of pharmacological and neuromodulation therapeutic strategies focused on improving cognitive outcomes and survival.

Each participant in the study was recruited from a prospective registry of adults aged 18–85 with newly diagnosed frontal, temporal and parietal IDH-WT high-grade gliomas with detailed language assessments and baseline MEG recordings. Inclusionary criteria included the following: native English speaking, aged 18–85 years, and no previous history of psychiatric illness, neurological illness, or drug or alcohol abuse. All human electrocorticography data were obtained during lexical retrieval language tasks from 14 adult awake patients undergoing intraoperative brain mapping for surgical resection. Tumours from eight patients were used for RNA-seq experiments. Site-directed tumour biopsies from 19 patients were used for immunofluorescence/immunohistochemistry analysis and 24 patients were used for immunocytochemistry and cell-based functional assays. Tumours from eight patients were used for mouse xenograft experiments. All of the participants provided written informed consent to participate in this study, which was approved by the University of California, San Francisco (UCSF) institutional review board (IRB) for human research (UCSF CC-171027, CHR 17-23215) and performed in accordance with the Declaration of Helsinki.

We began by studying short-range circuit dynamics in a subset of 14 patients with dominant hemisphere glioblastoma infiltrating speech production areas of the inferior frontal lobe using ECoG in the intraoperative setting (Extended Data Fig. 1a). We then focused molecular studies on patients with surgically treated IDH-WT glioblastoma, performed extraoperative language assessments and imaginary coherence as a long-range measure of functional connectivity using MEG (Extended Data Fig. 1a,b). This enabled us to import functional connectivity data into the operating room in which we performed site-specific tissue biopsies of human glioma from regions with differing measures of functional connectivity for in vivo and in vitro cell biology experiments including primary patient cultures (n = 19 patients) and multimodal tissue profiling, including microscopy, sequencing, proteomics and patient-derived tumour xenografting (Extended Data Fig. 1c). This layered approach—combining clinical variables, cognition assessments, human and animal models of network dynamics, in addition to cell biology—served as a platform to study the clinical implications of glioma–neuron interactions (Extended Data Fig. 1d and Supplementary Tables 1 and 2).

The hemisphere of language dominance was determined using baseline magnetic source imaging. In brief, the participants sat in a 275-channel whole-head CTF Omega 2000 system (CTF Systems) sampling at 1,200 Hz while they performed an auditory-verb generation task. The resulting time series were then reconstructed in source space with an adaptive spatial filter after registration with high-resolution MRI. Finally, changes in beta-band activity during verb generation were compared across hemispheres to generate an overall laterality index. All of the participants were left-dominant and underwent electrophysiological recording of the left hemisphere. We implemented an intraoperative testing paradigm that was previously established9. Noise in the operating room was minimized through rigorous enforcement of the following: (1) all personnel were requested to cease verbal communication; (2) telephones and alarms were muted; and (3) surgical suction and all other non-essential machinery were temporarily shut down. A 15 inch laptop computer (60 Hz refresh rate) running a custom MATLAB script integrated with PsychToolbox 3 (http://psychtoolbox.org/) was placed 30 cm away from each participant. The script initiated a picture-naming task that consisted of a single block of 48 unique stimuli, each depicting a common object or animal through coloured line drawings. Each stimulus was presented at the point of central fixation and occupied 75% of the display. After presentation of each stimulus, the participants were required to vocalize a single word that best described the item.

Intraoperative photographs with and without subdural electrodes present were used to localize each electrode contact combined with stereotactic techniques9,53. Images were registered using landmarks from gyral anatomy and vascular arrangement to preoperative T1- and T2-weighted MRI scans. Tumour boundaries were localized on MRI scans and electrodes within 10 mm of necrotic tumour core tissue were identified as ‘tumour’ contacts. Electrodes overlying the hypointense core of the tumour extending from the contrast enhancing rim to the edge of FLAIR were considered to be tumour electrodes, and electrodes completely outside of any T1 post gadolinium or FLAIR signal were considered to be non-tumour or normal appearing by a trained co-author blinded to the electrophysiologic data2. Glioma-infiltrated regions were defined on the basis of two criteria previously established in the literature9, including mass-like region of T2-weighted FLAIR sequences signal. Imaging was confirmed by gross inspection of the cortex confirming dilation and/or an abnormal vascular pattern. Previous research has shown that regions of non-enhancing disease consist of infiltrating tumour cells intermixed with neurons and normal glial cells2,54. These labels were reviewed by the study principal investigator (S.L.H.-J.) and compared to labels derived during intraoperative stereotactic neuronavigation to reach a consensus (Brainlab).

Each participant received a training session 2 days before participation to ensure familiarity with the task. ECoG signals were acquired during a period after stopping the administration of anaesthetics (minimum drug wash out period of 20 min) and the patient was judged to be alert and awake after an extensive post-emergence wakefulness assessment to ensure adequate arousal55. Intraoperative tasks consisted of naming pictorial representations of common objects and animals (picture naming) and naming common objects and animals through auditory descriptions (auditory naming)56. Post-operative videos were reanalysed to ensure that all data were collected and correct responses only were included for analysis. Audio was sampled at 44.1 kHz from a dual-channel microphone placed 5 cm from the participant and electrophysiological signals were amplified (g.tec). Recordings were acquired at 4,800 Hz and downsampled to 1,200 Hz during the initial stages of processing. During offline analyses, audio and electrophysiological recordings were manually aligned, resampled and segmented into epochs (speech-locked). These epochs set time = 0 ms as speech onset and included ±2,000 ms for a total of 4,000 ms of signal per trial. Trials were discarded if (1) an incorrect response was given (including fillers and interjections) or (2) there was a greater than 2 s delay between stimulus presentation and response so as to maintain consistent trial dynamics and ensure that the neural signal indeed reflected the experimental manipulations. Channels with excessive noise artifacts were visually identified and removed if their kurtosis exceeded 5.0. After the rejection of artifactual channels, data were referenced to a common average, high-pass filtered at 0.1 Hz to remove slow-drift artifacts, and bandpass filtered between 70–110 Hz using a 300-Order FIR filter to focus the analyses on the high-gamma band range, which is strongly related to local mean population spiking rates. To extract the ERSPs, electrophysiological signals were first downsampled to 600 Hz, then high-pass filtered at 0.1 Hz to remove DC-offset and low-frequency drift, notch-filtered at 60 Hz and its harmonics to remove line noise, and bandpass-filtered between 70 and 170 Hz (that is, the high-gamma range) using a Hamming windowed sinc FIR filter. These signals were finally smoothed using a 100 ms Gaussian kernel, downsampled to 100 Hz and z-scored across each trial. Electrodes were subsequently rereferenced to the common average for each participant to facilitate group comparisons, and regions of interest were defined according to the Automated Anatomical Labelling Atlas (https://www.gin.cnrs.fr/en/tools/aal/). The location of grid implantation was solely directed by clinical indications. The accuracy of the final registration for each participant was independently confirmed using gyral and sulcal anatomy to triangulate the location of each electrode registered to the template surface and was then compared to intraoperative photographs of the actual cortex with the overlying grid(s)57. The HGp was then calculated using the square of the Hilbert transform on the filtered data. The HGp was then averaged across the resting-state time series, yielding a single measure of neural responsivity for each electrode contact. The HGp was then averaged across patients during the task response period, yielding a single measure of neuronal responsivity for each channel. The HGp levels were then compared between tumour and normal appearing channels. Linear mixed-effects modelling was used to perform statistical comparisons with repeated measures using the nlme package in R (v.3.1-161; https://cran.r-project.org/web/packages/nlme/citation.html). The signal's origin (that is, normal-appearing/glioma-infiltrated cortex) was modelled as a fixed effect and the participants were modelled as random effects. For continuous variables without repeated measures, t-tests were used. A threshold of P < 0.05 was used to denote statistical significance and corrections for multiple comparisons were made using the Bonferroni method.

To decode between low-frequency words (for example, rooster) and high-frequency words (for example, car), signals from normal-appearing and glioma-infiltrated electrodes were extracted from the anterior temporal lobe after participant-level registration to a common MNI atlas. Responses were time-locked to speech onset and the signal envelope was extracted using a Hilbert transform after applying a bandpass filter in the high-gamma range (70–170 Hz). Subsequently, an l2-regularized logistic regression classifier was trained (cost of 1) to distinguish neural responses during vocalization of low-frequency words (for example, rooster) from high-frequency words (for example, car). Model performance was determined by taking the accuracy on a held-out participant and averaging it across all folds (that is, leave-one-participant-out cross-validation) and statistical significance was determined by testing this accuracy against a binomial distribution. This process was conducted separately for normal-appearing and glioma-infiltrated cortex using an identical preprocessing, training and testing paradigm.

MEG recordings were performed according to an established protocol30,31. In brief, the study participants had continuous resting state MEG recorded with a 275-channel whole-head CTF Omega 2000 system (CTF Systems) using a sampling rate of 1,200 Hz. During resting-state recordings, the participants were awake with their eyes closed. Surface landmarks were co-registered to structural magnetic resonance images to generate the head shape. Within the alpha frequency band, an artifact-free 1 min epoch was selected for further analysis if the patient's head movement did not exceed 0.5 cm. This artifact-free, 1 min epoch was then analysed using the NUTMEG software suite (v.4; UCSF Biomagnetic Imaging Laboratory) to reconstruct whole-brain oscillatory activity from MEG sensors so as to construct functional connectivity (imaginary coherence (IC)) metrics54,58,59. Spatially normalized structural magnetic resonance images were used to overlay a volume-of-interest projection (grid size = 8 mm; approximately 3,000 voxels per participant) such that each voxel contained the entire time series of activity for that location derived by all the MEG sensor recordings. The time series within each voxel was then bandpass-filtered for the alpha band (8–12 Hz) and reconstructed in source space using a minimum-variance adaptive spatial filtering technique54,60. The alpha frequency band was selected because it was the most consistently identified peak in the power spectra from this sampling window in our patient series. Functional connectivity estimates were calculated using IC, a technique known to reduce overestimation biases in MEG data generated from common references, cross-talk and volume conduction26,28.

Resting-state MEG was also used to measure intratumoural gamma activity. A spatial beamformer was applied to extract neural signals at the voxel level from manually defined regions of interest corresponding to FLAIR signal abnormality (that is, within the infiltrative margin of the tumour)61. These source-space signals were then downsampled to 300 Hz, notch filtered at 60 Hz to remove line noise and rereferenced to the common average. Spectral activity from 1 to 50 Hz was estimated at each voxel using Thomson's multitaper method (pmtm in MATLAB R2021b) with 29 Slepian tapers. Next, gamma power from 30 to 50 Hz was computed at each voxel after subtracting the aperiodic component from each spectrum by fitting a Lorentzian function in semi-log space62. A point estimate of intratumoural gamma activity was subsequently computed by averaging the activity across all voxels for each participant and regressed against the corresponding number of manually counted intratumoural HFC nodes.

The functional connectivity of an individual voxel was derived by the mean IC between the index voxel and the rest of the brain, referenced to its contralesional pair30. It is possible that there are regions within gliomas with varying amounts of functional connectivity. Moreover, there are individual patients with more or less functional connectivity. We have addressed these differences in our experimental model. Intratumoural differences in functional connectivity were addressed by the following: in comparison to contralesional voxels, we used a two-tailed t-test to test the null hypothesis that the Z-transformed connectivity IC between the index voxel and non-tumour voxel is equal to the mean of the Z-transformed connectivity between all contralateral voxels and the same set of voxels. The resultant functional connectivity values were separated into tertiles: upper tertile (HFC) and lower tertile (LFC). Functional connectivity maps were created by projecting connectivity data onto each individual patient's preoperative structural magnetic resonance images and imported into the operating room neuronavigation console. Stereotactic site-directed biopsies from HFC (upper tertile) and LFC (lower tertile) intratumoural regions were taken and x, y, z coordinates determined using Brainlab neuro-navigation. Thus, only the extremes of intratumoural connectivity (high and low connectivity, HFC and LFC, respectively) were analysed for these experiments. Rather than raw values, each functional connectivity measure represents a Z-transformed value and it therefore remains likely that the HFC distinction for one patient does not perfectly coincide with the HFC distinction in another patient's tumour (intertumoural heterogeneity).

Pre-operative and post-operative tumour volumes were quantified using BrainLab Smartbrush (v.2.6; Brainlab). Pre-operative MRI scans were obtained within 24 h before resection, and post-operative scans were all obtained within 72 h after resection. Total contrast-enhancing tumour volumes were measured at both pre-operative and post-operative timepoints. The total contrast-enhancing tumour volume was measured on T1-weighted post-contrast images, and the non- enhancing tumour volume was measured on T2 or FLAIR sequences. Manual segmentation was performed with region-of-interest analysis ‘painting’ inclusion regions based on fluid-attenuated inversion-recovery (FLAIR) sequences from pre- and post-operative MRI scans to quantify tumour volume. The extent of resection was calculated as follows: (pre-operative tumour volume − post-operative tumour volume)/pre-operative tumour volume × 100%. Manual segmentations were performed for which the tumour volumetric measurements were verified for accuracy after an initial training period. Volumetric measurements were performed blinded to patients’ clinical outcomes. All of the patients in the cohort had available preoperative and postoperative MRI scans for analysis. To ensure that post-operative FLAIR signal was not surgically induced oedema or ischaemia, FLAIR pre- and post-operative MRIs were carefully compared alongside DWI sequences before including each region in the volume segmentation42. HFC voxels with T1 post gadolinium contrast enhancing tumour were considered to be HFC-positive for survival analysis.

One to two days before tumour resection, patients underwent baseline language evaluation, which consisted of naming pictorial representations of common objects and animals (picture naming) and naming common objects and animals through auditory descriptions (auditory naming). Visual picture naming and auditory stimulus naming testing were used given their known significance and clinical correlation with outcomes in clinical patient population63,64. The correct answers for these tasks (delivered on a laptop with a 15 inch monitor (60 Hz refresh rate) positioned two feet away from the seated patient in a quiet clinical setting) were matched on word frequency (that is, commonality within the English language) using SUBTLEXWF scores provided by the Elixcon project and content category. Task stimuli were randomized and presented using PsychToolbox. The task order was randomly selected by the psychometrist for each participant. Slides were manually advanced by the psychometrist either immediately after the participant provided a response or after 6 s if no response was given. The tasks were scored on a scale from 0 to 4 by a trained clinical research coordinator who was initially blinded to all clinical data (including imaging studies). No participants had uncorrectable visual or hearing loss. Details of the administration and scoring of auditory and picture naming language tasks can be found in previous studies27,55,65.

Tumour tissues with high (HFC) and low (LFC) functional connectivity sampled during surgery based on preoperative MEG were processed for quality control by a certified neuropathologist and were subsequently used to generate primary patient-derived cultures. Patient-matched samples were acquired from site-directed HFC and LFC intratumoural regions from the same patient. Intratumoural HFC and LFC tissues were dissociated both mechanically and enzymatically and then passed through a 40 µm filter to remove debris. The filtered cell suspension was then treated with ACK lysis buffer (Invitrogen) to remove red blood cells and subsequently cultured as free-floating neurospheres in a defined, serum-free medium designated tumour sphere culture medium, consisting of Dulbecco's modified Eagle's medium (DMEM-F12; Invitrogen), B27 (Invitrogen), N2 (Invitrogen), human-EGF (20 ng ml−1; Peprotech), human-FGF (20 ng ml−1; Peprotech). Normocin (InvivoGen) was also added to the cell culture medium in combination with penicillin–streptomycin (Invitrogen) to prevent mycoplasma, bacterial and fungal contaminations. Cell cultures were routinely tested for mycoplasma (PCR Mycoplasma Test Kit I/C, PromoCell) and no positive results were obtained (Extended Data Fig. 7d).

RNA was isolated from HFC (n = 3) and LFC (n = 4) tumour samples using the RNeasy Plus Universal Mini Kit (QIAGEN) and RNA quality was confirmed using the Advanced Analytical Fragment Analyzer. RNA-seq libraries were generated using the TruSeq Stranded RNA Library Prep Kit v2 (RS-122- 2001, Illumina) and 100 bp paired-end reads were sequenced on the Illumina HiSeq 2500 system to at least 26 million reads per sample at the Functional Genomics Core Facility at UCSF. Quality control of FASTQ files was performed using FASTQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Reads were trimmed with Trimmomatic (v.0.32)66 to remove leading and trailing bases with quality scores of less than 20 as well as any bases that did not have an average quality score of 20 within a sliding window of 4 bases. Any reads shorter than 72 bases after trimming were removed. Reads were subsequently mapped to the human reference genome GRCh38 (https://www.ncbi.nlm.nih.gov/assembly/GCF_000001405.39/)67 using HISAT268 (v.2.1.0) with the default parameters. For differential expression analysis, we extracted exon-level count data from the mapped HISAT2 output using featureCounts69. Differentially expression analysis was performed using DESeq270 using the apeglm parameter71 to accurately calculate log-transformed fold changes and setting a false-discovery rate of 0.05. Differentially expressed genes were identified as those with log-transformed fold changes of greater than 1 and an adjusted P value of less than 0.05. Unsupervised gene expression principal component analysis and volcano plots of IDH-WT glioblastoma (Extended Data Fig. 5b,c) revealed 144 differentially expressed genes between HFC and LFC tumour regions, including 40 genes involved in nervous system development (Supplementary Table 3).

Fresh tumour samples were acquired from the operating room and transported to the laboratory space in PBS and on ice. Tumour tissue was minced with #10 scalpels (Integra LifeSciences) and then digested in papain (Worthington Biochemical, LK003178) for 45 min at 37 °C. Digested tumour tissue was then incubated in red blood cell lysis buffer (eBioscience, 00-4300-54) for 10 min at room temperature. Finally, the samples were sequentially filtered through 70 μm and 40 μm filters to generate a single-cell suspension.

Single-cell suspensions of three patient-matched HFC and LFC tumour tissues were generated as described above and processed for single-cell RNA-seq using the Chromium Next GEM Single Cell 3′ GEM, Library & Gel Bead Kit v3.1 on the 10x Chromium controller (10x Genomics) using the manufacturer's recommended default protocol and settings, at a target cell recovery of 5,000 cells per sample. Although single-cell sequencing does not capture all cell types within the central nervous system microenvironment, the sequencing pipeline used in this study has been demonstrated to identify neurons and was therefore chosen for use in physiologically annotated fresh glioblastoma samples, compared with single-nucleus RNA-seq, which is commonly applied for frozen archived tissues72,73. One hundred base pair paired-end reads were sequenced on the Illumina NovaSeq 6000 system at the Center for Advanced Technology at the University of California San Francisco, and the resulting FASTQ files were processed using the CellRanger analysis suite (v.3.0.2; https://github.com/10XGenomics/cellranger) for alignment to the hg38 reference genome, identification of empty droplets, and determination of the count threshold for further analysis. A cell quality filter of greater than 500 features but fewer than 10,000 features per cell, and less than 20% of read counts attributed to mitochondrial genes, was used. Single-cell UMI count data were preprocessed in Seurat (v.3.0.1)74,75 using the sctransform workflow76, with scaling based on the regression of UMI count and the percentage of reads attributed to mitochondrial genes per cell. Dimensionality reduction was performed using principal component analysis and then principal component loadings were corrected for batch effects using Harmony77. Uniform manifold approximation and projection was performed on the reduced data with a minimum distance metric of 0.4 and Louvain clustering was performed using a resolution of 0.2. Marker selection was performed in Seurat using a minimum difference in the fraction of detection of 0.5 and a minimum log-transformed fold change of 0.5. We assessed the single-cell transcriptome from 6,666 HFC-region cells and 7,065 LFC-region cells (Supplementary Table 4).

After rehydration, 5.0 μm paraffin-embedded sections were processed for antigen retrieval followed by blocking and primary antibody incubation overnight at 4 °C. The following primary antibodies were used: rabbit anti-synapsin 1 (1:1,000, EMD Millipore), mouse anti-PSD95 (1:100, UC Davis), mouse anti-nestin (1:500, Abcam), mouse anti-neurofilament (M+H; 1:1,000, Novus Biologicals), mouse anti-TSP-1 (1:20, Invitrogen), rabbit anti-TSP-1 (1:50, Abcam), rabbit anti-MET (1:100, Abcam) and rabbit anti-Ki-67 (1:100, Abcam). We used species-specific secondary antibodies: Alexa 488 goat anti-chicken IgG, Alexa 488 goat anti-rabbit IgG, Alexa 568 goat anti-rabbit IgG, Alexa 568 goat anti-mouse IgG, Alexa 647 goat anti-rabbit IgG, all used at 1:500 (Invitrogen). After DAPI nuclear counter staining (Vector Laboratories, 1:1,000), coverslips were mounted with Fluoromount-G mounting medium (SouthernBiotech) for immunofluorescence analysis. The number of synapsin-1 and PSD95 puncta was quantified using spots (with automatic intensity maximum spot detection thresholds and a spot diameter of 1.0 µm) detection function of Imaris. The ratio of pre- and postsynaptic puncta was calculated by dividing the total number of synapsin-1 or PSD95 puncta on neurofilament-positive neurons to the total number of cells stained with DAPI in 135 μm × 135 μm field areas for quantification. Alternatively, the sections were incubated in DAB horseradish peroxidase (Vector Laboratories) for chemical colorimetric detection after incubation in ImmPress anti-rabbit IgG (Novus Biologicals) and counterstained with Harris haematoxylin for immunohistochemistry analysis.

Glioma cells were plated on poly-d-lysine and laminin-coated coverslips (Neuvitro) at a density of 10,000 cells per well in 24-well plates. Approximately 24 h later, 40,000 embryonic mouse hippocampal neurons (Gibco) were seeded on top of the glioma cells and maintained with serum-free Neurobasal medium supplemented with B27, gentamicin and GlutaMAX (Gibco). After 2 weeks of co-culture, cells were fixed with 4% paraformaldehyde (PFA) for 30 min at 4 °C and incubated in blocking solution (5% normal donkey and goat serum, 0.25% Triton X-100 in PBS) at room temperature for 1 h. Next, they were treated with primary antibodies diluted in the blocking solutions overnight at 4 °C. The following antibodies were used: rabbit anti-homer-1 (1:250, Pierce), mouse anti-synapsin-1 (1:200) and chicken anti-MAP2 (1:500, Abcam). The coverslips were then rinsed three times in PBS and incubated in secondary antibody solution (Alexa 488 goat anti-chicken IgG; Alexa 568 goat anti-mouse IgG, and Alexa 647 goat anti-rabbit IgG, all used at 1:500 (Invitrogen) in antibody diluent solution for 1 h at room temperature. The coverslips were rinsed three times in PBS and then mounted with VECTA antifade mounting medium with DAPI (Vector Laboratories).

Images were captured at 1,024 × 1,024 resolution using a ×10 objective on the Nikon C2 confocal microscope. The confocal microscope settings for the homer-1 Alexa488 and synapsin-1 Alexa647 channels were held constant across all of the samples that were used for the experiment. Collected images were then imported into Imaris software (Imaris v.9.2.1, Bitplane) and the threshold value for each channel was manually adjusted and the colocalized voxels of the synapsin-1 puncta with the homer-1 marker was detected by creating a colocalization channel using the built-in the colocalization module of the Imaris software. Furthermore, the colocalization events were quantified by running the built-in spot detection algorithm of Imaris in conjunction with the colocalization channel. Next, dendrites labelled by MAP2 was visualized in TRITC channel and reconstructed using the Filament tool of Imaris software; the number of colocalized puncta representing synapses were counted and presented as the number of synapsin-1- and homer-1-positive puncta per 10 μm of dendrite length. Areas of homer-1 immunolabelled synaptic puncta were reconstructed using Imaris software Surface tool on maximal-intensity projections. Surfaces were built using a surface area detail level of 0.1 µm, thresholding by absolute intensity and taking all voxel >1.0 into account. The area sizes of individual anti-homer1-immunostained puncta were analysed and the mean values were calculated.

Induced neuron organoids were generated from a WTC11 iPS cell clone integrated by human NGN2 transgene induction as described previously38,78. In brief, iNeuron organoids were generated by the transgenic human iPS cell WTC11 line by NGN2 induction through addition of 2 μg ml−1 doxycycline in the 1:1 mixture of Neurobasal and BrainPhys neural medium containing 1% B-27 supplement, 0.5% GlutaMAX, 0.2 μM compound E, 10 ng ml−1 BDNF and 10 ng ml−1 NT-3 for 10 days to induce neuronal differentiation. Next, neuron maturation was triggered by feeding the organoids with approximately 8-month-old organoid conditioned medium derived from astrocytes. Astrocytes were differentiated from the human iPS cell WTC11 line and cultured in a medium consisting of DMEM/F12 containing GlutaMAX, sodium bicarbonate, sodium pyruvate, N-2 supplement, B-27 supplement (Gibco), 2 μg ml−1 heparin, 10  ng ml−1 EGF and 10 ng ml−1 FGF2. Neuron organoids were characterized as postmitotic and stained for MAP2 and βIII-tubulin to validate neuronal induction efficiency. After 14 days of neuronal differentiation, HFC/LFC glioma cells labelled with RFP were added to the neuron organoid culture at a ratio of 1:3. Before iNeuron induction, the transgenic human iPS cell line WTC11 was transduced with GFP lentivirus. A Zeiss Cell Observer spinning-disc confocal microscope (Carl Zeiss) fitted with a temperature- and carbon-dioxide-controlled chamber was used to record live interactions of glioma cells with neuron organoids. Organoids were imaged every 10 min for a 6 h period, starting at the time of co-culture initiation, using a 10× objective with 0.4 NA. To assess the effect of exogenous TSP-1 on the functional integration between glioma cells and neurons, human recombinant TSP-1 (R&D Systems) was applied at a dose of 5 µg ml−1 to the LFC-neuron organoid co-culture. Live-cell image analyses were performed using ImageJ. In brief, a region of interest was drawn around each GFP-positive neuron organoid and the fluorescence intensity (integrated density) of the RFP-positive glioblastoma cells was measured in the outlined regions of interest for each of the indicated timepoints. At the end of two weeks, organoids from HFC and LFC co-cultures were embedded in OCT and sectioned at 10 μm thickness for homer-1 immunofluorescence staining. Determination of homer-1 expression was performed by analysing homer-1 puncta density of neuron-organoid-HFC and LFC co-cultures.

We prepared 24-well CytoView multi-electrode plates (Axion Biosystems) before the addition of cells by coating with poly-d-lysine (Thermo Fisher Scientific), laminin (Fisher Scientific) and fibronectin (Corning). In brief, 1 day before the establishment of cultures, a solution of 0.1 mg ml−1 of poly-d-lysine was added to the MEA plates at a volume of 100 μl per well and incubated at room temperature for 2 h. After 2 h, poly-d-lysine was aspirated and the plates were washed three times with sterile water, and allowed to air dry in a biosafety cabinet and stored at 4 °C. The next day, the plates were coated with 100 μl of 5 µg ml−1 laminin and 1 μg ml−1 fibronectin and incubated for 2 h at 37 °C before cell seeding.

Primary cortical cultures were established from E18 CD1 mice (Charles River Laboratories). Timed-pregnant CD1 dams were killed by CO2 euthanasia in accordance with UCSF Institutional Animal Care and Use Committee (IACUC). Dissection of complete cortex from E18 embryos was performed in ice-cold HBSS (Gibco) under a dissecting microscope (Zeiss). Dissected cortices were minced to 1 mm2 pieces and enzymatically digested in 5 ml of 0.25% trypsin reconstituted from 2.5% trypsin (Corning) in calcium- and magnesium-free Hank's Balanced Salt Solution (Worthington Biochemical Corporation) for 30 min at 37 °C. Then, 0.5 ml of 10 mg ml−1 of DNase (Sigma-Aldrich) was added in the last 5 min of dissociation. Mechanical dissociation was then carried about by trituration using fire-polished glass Pasteur pipettes until tissue was homogeneously suspended with no visible sections/aggregates and subsequently filtered through a 40 μm cell strainer (Thermo Fisher Scientific). Cells were collected by centrifugation at 500g for 5 min and the resulting cell pellet was resuspended in fresh complete BrainPhys culture medium (1× BrainPhys culture medium (StemCell Technologies) supplemented with B27 (Invitrogen), N2 (Invitrogen) and penicillin–streptomycin antibiotics (Invitrogen). Then, 10 μl of cell suspension was mixed 1:1 with Trypan Blue, and the viable cell concentration was quantified using a haemocytometer. Further dilution was performed to bring viable cell concentration to 100,000 cells per 10 μl. Droplets of 10 μl were then added directly over the electrode field of each pretreated MEA well and stored in the cell culture incubator for 1 h to allow cell adhesion. The wells were then carefully flooded with 500 μl complete BrainPhys medium and the cultures were maintained and allowed to mature in a tissue culture incubator with semi-weekly half-volume medium changes.

Spontaneous extracellular neuronal recordings were carried out using the Maestro Edge system with an integrated heating system and temperature controller (Axion Biosystems) in combination with the Axion 24-well CytoView MEA plates (each well housing a 4 × 4 16-channel electrode array that are 350 µm away from each other) and Axion Integrated Studio (AxIS) Navigator (v.3.5.2; Axion Biosystems). In brief, to record spontaneous neuronal activity, the Neural Real-Time module was used. The neuronal firing events/action potentials (herein referred to as the spike) was defined by applying an adaptive threshold crossing method, that sets the threshold for spike detection for each channel/electrode to 5 s.d. of the noise level79; activity exceeding this threshold was counted as a spike. Unless otherwise stated, all analysis considers only active channels, defined as channels exhibiting ≥5 spikes per min. Raw data files were obtained by sampling the channels simultaneously with a gain of 1,000× and a sampling frequency of 12.5 kHz per channel using a band-pass filter (200–3,000 Hz). To detect single-electrode bursting activity, an interspike interval threshold was used, setting the minimum number of spikes at 5 and the maximum interspike interval at 100 ms. The network bursting activity (simultaneous bursts at multiple MEA electrodes) was analysed by Neural Metric Tool (v.1.2.3; Axion Biosystems). For this purpose, the Adaptive algorithm was selected using the following settings: minimum number of spikes = 50 and minimum electrodes = 35%. Quantification of network synchrony was computed through AxIS software by calculating the area under the normalized cross-correlogram (AUNCC) as described previously80,81,82,83. AUNCC represents the area under interelectrode cross-correlation normalized to the autocorrelations, with higher values indicating greater synchronicity of the network. For additional neural data analysis, including mean firing rate of each electrode (the ratio of the total number of spikes per second and the total duration of recording (1,800 s)) and weighted mean firing rate (defined as the spike rate per well multiplied by the number of active electrodes in the associated well), raw data files were processed offline using the Statistics Compiler function in AxIS. Statistics Compiler output files were processed in Microsoft Excel (Microsoft) and with custom Python scripts to organize and extract individual parameter data for each well of each MEA plate and for data normalization. Raster plots illustrating spike histogram and network bursts were generated using Neural Metric Tool (Axion Biosystems).

Spontaneous neuronal activity from cortical cultures grown on MEA plates was recorded in 30 min sessions on days in vitro 1 (DIV1), DIV7 and DIV15. Bright-field images were captured at each of the above timepoints to assess the neuronal cell density and electrode coverage. Primary cortical neurons showed a constant maturation trend from DIV7 to DIV15, and the co-culture experiments were initiated when neurons showed a synchronous activity pattern network at DIV15. Baseline data were therefore recorded on DIV15 immediately before addition of glioma cells in the presence or absence of gabapentin. For glioma cell co-culture, a single-cell suspension from cultured neurospheres of primary patient-derived HFC and LFC were prepared and diluted to a viable cell concentration of 20,000 cells per 5 μl. Droplets of 5 μl were then plated on top of differentiating neurons in the MEA plate. After plating, glioma cells were allowed to adhere for approximately 1 h after which HFC cultures were exposed for next 24–48 h to a working concentration of 50 µM gabapentin84,85 diluted in complete BrainPhys medium or an equivalent amount of vehicle (sterile water) as a control. Each condition was run on two wells (experimental replicates). Neurons from two different embryos were used as biological triplicates (n = 2). Presented data from MEA recordings reflects well-wide averages from active electrodes, with the number of wells per condition represented by n values.

All in vivo experiments were conducted in accordance with the protocols approved by the UCSF Institutional Animal Care and Use Committee (IACUC) and performed in accordance with institutional guidelines. Animals were maintained under pathogen-free conditions, in temperature- and humidity-controlled housing, with free access to food and water, under a 12 h–12 h light–dark cycle. For brain tumour xenograft experiments, the IACUC does not set a limit on maximal tumour volume but rather on indications of morbidity. These limits were not exceeded in any of the experiments as mice were euthanized if they exhibited signs of neurological morbidity or lost 15% or more of their body weight.

For all xenograft studies, NSG mice (NOD-SCID-IL2R gamma chain-deficient, The Jackson Laboratory) were used. Male and female mice were used equally. For immuno-electron microscopy experiments, a single-cell suspension from cultured neurospheres of HFC and LFC (n = 2 each) labelled with red fluorescent protein (RFP) were prepared in sterile DMEM immediately before the xenograft procedure. Mice (n = 8; 2 biological replicates per patient line) at postnatal day 28–30 were anaesthetized with 1–4% isoflurane and placed into a stereotactic apparatus. The cranium was exposed through a midline incision under aseptic conditions. Approximately 50,000 cells in 2 µl sterile PBS were stereotactically implanted into the CA1 region of the hippocampus through a 31-gauge burr hole, using a digital pump at infusion rate of 0.4 µl min−1 and 31-gauge Hamilton syringe. Stereotactic coordinates used were as follows: 1.5 mm lateral to midline, 1.8 mm posterior to bregma, −1.4 mm deep to cranial surface. At the completion of infusion, the syringe needle was allowed to remain in place for a minimum of 2 min, then manually withdrawn at a rate of 0.875 mm min−1 to minimize backflow of the injected cell suspension. For survival studies, morbidity criteria used were either: reduction of weight by 15% initial weight, or clinical signs such as hunched posture, lethargy or persistent decumbency. Kaplan–Meier survival analysis using log-rank testing was performed to determine statistical significance.

Cell quantification was performed by a blinded investigator at 10–20× magnification using the Zeiss LSM800 scanning confocal microscope and Zen 2011 imaging software (Carl Zeiss). The area for quantification was selected for a 1-in-6 series of 40 μm coronal sections (240 µm apart from one another). Immunohistochemistry was performed on brain sections from HFC and LFC xenografts to stain for human nuclear antigen (HNA)-positive tumour cells. The tumour burden was evaluated using blinded rank-order analysis as previously reported86. For each mouse, the section with the maximal amount of tumour burden was selected as defined by the number of HNA-positive cells. From this section, a tiled ×10 image of the centre of the section (3 × 3 tiles stitched together) was created, therefore generating a single image for each mouse that represents the maximal amount of tumour burden. Next, the images generated from each mouse were compared and ranked in order from the image with the least to maximum number of cells. Subsequently, scores were assigned and ranged from 1 to 24 with 1 and 24 representing images with the lowest and highest number of HNA-positive cells, respectively. After scoring, the experimenter was unblinded and the scores from each image were assigned to each of the respective HFC or LFC experimental conditions. Statistical differences between HFC and LFC scores were evaluated using two-tailed unpaired Mann–Whitney tests.

Twelve weeks after xenografting, mice were euthanized by transcardial perfusion with Karnovsky's fixative: 2% glutaraldehyde (EMS, 16000) and 4% PFA (EMS, 15700) in 0.1 M sodium cacodylate (EMS, 12300), pH 7.4. Transmission electron microscopy (TEM) was performed in the tumour mass within the CA1 region of the hippocampus for all xenograft analysis. The samples were then post-fixed in 1% osmium tetroxide (EMS, 19100) for 1 h at 4 °C, washed three times with ultrafiltered water, then en bloc stained overnight at 4 °C. The samples were dehydrated in graded ethanol (50%, 75% and 95%) for 15 min each at 4 °C; the samples were then allowed to equilibrate to room temperature and were rinsed in 100% ethanol twice, followed by acetonitrile for 15 min. The samples were infiltrated with EMbed-812 resin (EMS, 14120) mixed 1:1 with acetonitrile for 2 h followed by 2:1 EMbed-812:acetonitrile overnight. The samples were then placed into EMbed-812 for 2 h, then placed into TAAB capsules filled with fresh resin, which were then placed into a 65 °C oven overnight. Sections were taken between 40 nm and 60 nm on a Leica Ultracut S (Leica) and mounted on 100-mesh Ni grids (EMS FCF100-Ni). For immunohistochemistry, microetching was done with 10% periodic acid and eluting of osmium with 10% sodium metaperiodate for 15 min at room temperature on parafilm. Grids were rinsed with water three times, followed by 0.5 M glycine quench, and then incubated in blocking solution (0.5% BSA, 0.5% ovalbumin in PBST) at room temperature for 20 min. Primary goat anti-RFP (1: 300, ABIN6254205) was diluted in the same blocking solution and incubated overnight at 4 °C. The next day, grids were rinsed in PBS three times, and incubated in secondary antibodies (1:10 10 nm gold-conjugated IgG, TED Pella, 15796) for 1 h at room temperature and rinsed with PBST followed by water. For each staining set, samples that did not contain any RFP-expressing cells were stained simultaneously to control for any non-specific binding. Grids were contrast stained for 30 s in 3.5% uranyl acetate in 50% acetone followed by staining in 0.2% lead citrate for 90 s. The samples were imaged using a JEOL JEM-1400 TEM at 120 kV and images were collected using a Gatan Orius digital camera.

Sections from the xenografted hippocampi of mice were imaged as above using TEM imaging. Here 101 sections of HFC xenografts across 4 mice and 104 sections of LFC xenografts across 3 mice were analysed. Electron microscopy images were taken at ×6,000 with a field of view of 15.75 μm2. Glioma cells were counted and analysed after unequivocal identification of a cluster of immunogold particle labelling with 15 or more particles. The total number of synapses, including neuron-to-neuron and neuron-to-glioma synapses (identified by: (1) the presence of synaptic vesicle clusters; (2) visually apparent synaptic cleft; and (3) identification of clear postsynaptic density in the glioma cell) were counted.

The EdU-incorporation assay was performed using an EdU assay kit (Invitrogen) according to the manufacturer's instructions. Patient-derived HFC/LFC glioma cells were seeded on poly-d-lysine- and laminin-coated coverslips at 10,000 cells per well of a 24-well plate. After 24 h of seeding, embryonic mouse hippocampal neurons were added to the glioma cells at 40,000 cells per well for the glioma-neuron co-culture group. After 72 h, glioma cells alone or in coculture with neurons were treated with 20 μM EdU overnight at 37 °C. Subsequently, the cells were fixed with 4% PFA and stained using the Click-iT EdU kit protocol. The proliferation index was then determined by quantifying the fraction of EdU-labelled cells/DAPI-labelled cells using confocal microscopy at ×10 magnification.

Glioma cell invasion was evaluated by performing an invasion assay using the Cultrex 3D Spheroid Cell Invasion Assay Kit (Trevigen) according to the manufacturer's protocol. In brief, 3,000 cells were resuspended in 50 μl of spheroid formation matrix solution (prepared in culture medium) in a round-bottom 96-well plate. Spheroids were allowed to form for 72 h and images were taken at ×10 magnification before addition of invasion matrix (0 h). Working on ice, 50 μl of invasion matrix was then added to each well and the plate was incubated at 37 °C. After 1 h of gel formation, for the glioma cell + mouse conditioned medium (mCM) group, 100 μl of mouse hippocampal neuron supernatant was added to the wells to assess the effect of neuronal secreted factors on the invasiveness of the glioma cells. After 24 h of incubation at 37 °C, invasions were observed under a microscope and images were taken at ×10. Microtube lengths as well as the area of each spheroid measured at 0 h (pre-invasion) and 24 h (post-invasion) were analysed using ImageJ and the difference was used to calculate the total area of cell invasion.

SEM was performed as described previously39 to investigate the tumour microtubes in greater detail. Primary patient-derived HFC/LFC glioma cells were seeded on poly-d-lysine- and laminin-coated coverslips (Neuvitro) at 10,000 cells per well of a 24-well plate. For the glioma cell + mCM and glioma cell + mCM + GBP groups, mouse neuronal conditioned medium was added to the wells in the absence or presence of gabapentin (32 µM), respectively, to assess the effect of neuronal secreted factors and gabapentin modulation on tumour microtube formation of glioma cells47,87,88. After culture for 1 week, the samples were fixed with 4% PFA for 30 min at 4 °C, washed with PBS and ultrafiltered water. After serial stepwise ethanol dehydration, coverslips were mounted on SEM stubs (TED Pella) using conductive adhesive tape (12 mm OD PELCO Tabs, TED Pella). The samples were then sputter-coated with a 2 nm layer of gold palladium. Tumour microtubes were observed under field emission SEM (Sigma 500; Carl Zeiss Microscopy) and micrographs were recorded at an accelerating voltage of 1.0 kV.

Primary patient-derived HFC glioma cells were seeded in 6-well plates at a density of 2 × 105 cells per well. After overnight incubation, cells were infected with lentiviral particles expressing a shRNA targeting THBS1 or a control scramble construct expressing green fluorescent protein (GFP) according to the manufacturer's protocol using polybrene transfection reagent and subsequently kept in a 5% CO2 incubator at 37 °C for 24 h. After 24 h, the medium was replaced with fresh complete culture medium. Cells were checked 72 h after transduction for GFP expression to evaluate the efficiency of the transduction. The lentiviral shRNA constructs targeting THBS1 (5′-AGACATCTTCCAAGCATATAA-3′) and control scrambled shRNA (5′-CCTAAGGTTAAGTCGCCCTCG-3′) were designed and constructed by VectorBuilder.

The cell viability of primary patient-derived HFC cells after THBS1 knockdown was evaluated using the Live/Dead Viability/Cytotoxicity kit (Molecular Probes). One-week after transduction, HFC cells expressing scramble shRNA or shRNA targeting THBS1 were seeded on poly-d-lysine- and laminin-coated coverslips (Neuvitro) at 10,000 cells per well of a 24-well plate. After 1 week of culture, a 500 μl cell-staining solution at a final concentration of 2 µM calcein AM and 4 µM ethidium bromide (EthD-1) in DPBS was added to each well, and the plates were incubated for 45 min at room temperature in dark. As an indicator of cell viability, the calcein-AM is metabolically converted by intracellular esterase activity resulting in the green, fluorescent product, calcein. EthD-1 is excluded from live cells but is readily taken up by dead cells and stains the DNA emitting red fluorescence. Live and dead cells were imaged from four random fields per well and were visualized under a fluorescence microscope. The percentage of live cells was calculated as the number of live cells (in green) divided by the number of total cells (green + red) per image field.

After 1 week of transduction, primary patient-derived high connectivity glioma cells expressing scramble shRNA or shRNA against THBS1 were seeded on poly-d-lysine- and laminin-coated coverslips (Neuvitro) at 10,000 cells per well of a 24-well plate. Approximately 24 h later, 40,000 embryonic mouse hippocampal neurons (Gibco) were seeded on top of the glioma cells and maintained with serum-free Neurobasal medium supplemented with B27, gentamicin and GlutaMAX (Gibco). After 1 week of co-culture, cells were fixed with 4% PFA for 30 min at 4 °C and incubated in blocking solution (5% normal donkey and goat serum, 0.25% Triton X-100 in PBS) at room temperature for 1 h. Next, cells were treated with primary antibodies diluted in the blocking solutions overnight at 4 °C. The following antibodies were used: rabbit anti-Ki-67 (1:500, Abcam) and human nuclear antigen (HNA; mouse anti-human nuclei, 235-1; 1:100, Millipore). The coverslips were then rinsed three times in PBS and incubated in secondary antibody solution (Alexa 488 goat anti-rabbit IgG, and Alexa 647 goat anti-mouse IgG) all used at 1:500 (Invitrogen) in antibody diluent solution for 1 h at room temperature. The coverslips were rinsed three times in PBS and then mounted with VECTA antifade mounting medium with DAPI (Vector Laboratories). To calculate the proliferation index, the total number of HNA-positive cells co-labelled with Ki-67 was divided by the total number of human nuclei-labelled cells per image field visualized using confocal microscopy at ×40 magnification.

The effect of TSP-1 inhibition by gabapentin on glioma cell proliferation was evaluated by Ki-67 immunofluorescence staining, as described above. In brief, primary patient-derived high connectivity glioma cells were seeded on poly-d-lysine- and laminin-coated coverslips (Neuvitro) at 10,000 cells per well of a 24-well plate. Approximately 24 h later, 40,000 embryonic mouse hippocampal neurons (Gibco) were seeded on top of the glioma cells and maintained with serum-free Neurobasal medium supplemented with B27, gentamicin and GlutaMAX (Gibco). The next day, cultures were treated either with vehicle (sterile water) or 32 µM gabapentin followed by daily half-medium switches of fresh 32 µM gabapentin until 1 week of coculture88. Subsequently, cells were fixed and immunostained for Ki-67 and HNA labelling for proliferation assessment as described above.

A single-cell suspension from cultured neurospheres of primary patient-derived HFC and LFC were prepared in sterile HBSS immediately before the xenograft procedure by dissociation with TrypLE (Thermo Fisher Scientific). Mice (4–6 biological replicates per patient line) at postnatal day 28–30 were anaesthetized with 1–4% isoflurane and placed into a stereotactic apparatus. The cranium was exposed through midline incision under aseptic conditions. Approximately 150,000 cells in 3 µl sterile HBSS were stereotactically implanted into the premotor cortex (M2) through a 26-gauge burr hole, using a digital pump at infusion rate of 1.0 μl min−1. Stereotactic coordinates used were as follows: 1.0 mm lateral to midline, 1.0 mm anterior to bregma, −1.0 mm deep to cortical surface. Four weeks after xenograft, HFC/LFC-bearing mice were treated with systemic administration of gabapentin (200 mg kg−1; Sigma-Aldrich; formulated in saline) through intraperitoneal injection for 28 consecutive days. Controls were treated with an identical volume of the relevant vehicle. At 8 weeks after xenograft, mice were euthanized and coronal brain sections at 40 µm were obtained for immunohistochemistry.

Cell quantification was performed by a blinded investigator at ×10–20 magnification using the Zeiss LSM800 scanning confocal microscope and Zen imaging software (Carl Zeiss). The area for quantification was selected for a 1-in-6 series of 40 μm coronal sections (240 µm apart from one another). Within 3 sections of maximal tumour burden, all HNA-positive (mouse anti-human nuclei, 235-1; 1:100, Millipore) cells were quantified to determine the tumour burden within the areas quantified. HNA-positive tumour cells were then assessed for double-labelling with Ki-67. To calculate the proliferation index (the percentage of proliferating tumour cells for each animal), the total number of HNA-positive cells co-labelled with Ki-67 across all areas quantified was divided by the total number of human nuclei-positive cells counted across all areas quantified. Differences in proliferation indices were calculated using unpaired, one-tailed Student's t-tests.

Peripheral blood samples from newly diagnosed patients with glioblastoma (n = 56) were collected and allowed to clot for 30 min at room temperature before centrifugation for 15 min at 1,000g. The serum was stored at −80 °C until analysis. The TSP-1 level was determined using the Quantikine immunosorbent assay kits according to the manufacturer's instructions (R&D Systems). To confirm the functional protein-level knockdown of THBS1, the TSP-1 level was also measured in cell culture supernatants collected from scramble and THBS1-shRNA transduced HFC cells (collected one-week after infection) using the Picokine ELISA kit (Boster Biological Technology).

Statistical tests were conducted using Prism (GraphPad) software unless otherwise indicated. The significance test of different groups was determined using Student's t-tests and one-way ANOVA with Tukey's post hoc tests. P ≤ 0.05 was considered to be statistically significant. Two-tailed unpaired Mann–Whitney tests were used to analyse the tumour cell burden of HFC and LFC hippocampal xenografts using the same significance denotation as above. Two-tailed log-rank analyses were used to analyse the statistical significance of Kaplan–Meier survival curves for human patients. All of the micrographs shown are representative of three independently conducted experiments, with similar results obtained. Statistical analyses for RNA-seq data are described above in the respective sections. The tumour volumes of patients with glioblastoma were calculated by manual segmentation with region-of-interest analysis ‘painting’ inclusion regions on the basis of FLAIR sequences. Volumetric measurements were made blinded to patients’ clinical outcomes. Manual segmentations were performed by co-authors A.E., A.A. and A.L. with tumour volumetrics verified for accuracy after an initial training period. Student's t-tests and χ2 tests were used to compare continuous and categorical variables between patient cohorts, respectively. Patient overall survival (OS) was defined as the time from the date of first surgery or original biopsy (if it occurred before surgery) until death or the last contact date. Alive patients were censored at the time of loss to follow-up or last follow-up date. Median follow-up was estimated using the reverse Kaplan–Meier method.

To identify clinically relevant survival risk groups in a multivariate setting, we used recursive partitioning analyses for survival data using the partDSA algorithm (v.0.9.14)42,43,44. Survival trees use recursive partitioning to divide patients into different risk groups. The Brier score was the chosen loss function for splitting and pruning. Such methods are nonparametric and, therefore, do not require the proportional hazards assumption. All known prognostic variables were included in the trees, including age at diagnosis, sex, MGMT promoter methylation status, tumour location, chemotherapy, radiotherapy, the presence of functional connectivity within the tumour, pre- and post-operative tumour volume and the extent of resection. MGMT methylation status was included in both univariate Cox proportional-hazard modelling in addition to multivariate recursive partitioning survival analysis. The tree that minimized the fivefold cross-validated error as well as the most parsimonious tree within one standard error of the overall minimum error were selected for review. Leaves of the resulting trees defined the final risk groups from which the corresponding Kaplan–Meier curves were generated. Median OS times and hazard ratios were generated and compared between risk groups using the Kaplan–Meier method and the Cox proportional hazards model, respectively. The proportional hazards assumption was verified.

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Bulk and single-cell RNA-seq data of primary patient-derived samples reported in this manuscript have been deposited at the NCBI Gene Expression Omnibus under the accession code GSE223065. The publicly available GRCh38 (hg38; https://www.ncbi.nlm.nih.gov/assembly/GCF_000001405.39/) was used in this study. Source data are provided with this paper.

Venkataramani, V. et al. Glutamatergic synaptic input to glioma cells drives brain tumour progression. Nature 573, 532–538 (2019).

Article ADS CAS PubMed Google Scholar

Venkatesh, H. S. et al. Electrical and synaptic integration of glioma into neural circuits. Nature 573, 539–545 (2019).

Article ADS CAS PubMed PubMed Central Google Scholar

Venkatesh, H. S. et al. Neuronal activity promotes glioma growth through neuroligin-3 secretion. Cell 161, 803–816 (2015).

Article CAS PubMed PubMed Central Google Scholar

Venkatesh, H. S. et al. Targeting neuronal activity-regulated neuroligin-3 dependency in high-grade glioma. Nature 549, 533–537 (2017).

Article ADS PubMed PubMed Central Google Scholar

Buckingham, S. C. et al. Glutamate release by primary brain tumors induces epileptic activity. Nat. Med. 17, 1269–1274 (2011).

Article CAS PubMed PubMed Central Google Scholar

Campbell, S. L., Buckingham, S. C. & Sontheimer, H. Human glioma cells induce hyperexcitability in cortical networks. Epilepsia 53, 1360–1370 (2012).

Article CAS PubMed PubMed Central Google Scholar

Campbell, S. L. et al. GABAergic disinhibition and impaired KCC2 cotransporter activity underlie tumor-associated epilepsy. Glia 63, 23–36 (2015).

Article PubMed Google Scholar

John Lin, C. C. et al. Identification of diverse astrocyte populations and their malignant analogs. Nat. Neurosci. 20, 396–405 (2017).

Article CAS PubMed Google Scholar

Aabedi, A. A. et al. Functional alterations in cortical processing of speech in glioma-infiltrated cortex. Proc. Natl Acad. Sci. USA 118, e2108959118 (2021).

Article CAS PubMed PubMed Central Google Scholar

Dalrymple, S. J. et al. Changes in proliferating cell nuclear antigen expression in glioblastoma multiforme cells along a stereotactic biopsy trajectory. Neurosurgery 35, 1036–1044 (1994).

Article CAS PubMed Google Scholar

Hadjiabadi, D. H. et al. Brain tumors disrupt the resting-state connectome. Neuroimage Clin. 18, 279–289 (2018).

Article PubMed PubMed Central Google Scholar

Krishna, S. & Hervey-Jumper, S. L. Neural regulation of cancer: cancer-induced remodeling of the central nervous system. Adv. Biol. 6, e2200047 (2022).

Stoecklein, V. M. et al. Resting-state fMRI detects alterations in whole brain connectivity related to tumor biology in glioma patients. Neuro Oncol. 22, 1388–1398 (2020).

Article PubMed PubMed Central Google Scholar

Crone, N. E., Boatman, D., Gordon, B. & Hao, L. Induced electrocorticographic gamma activity during auditory perception. Clin. Neurophysiol. 112, 565–582 (2001).

Article CAS PubMed Google Scholar

Steinschneider, M., Fishman, Y. I. & Arezzo, J. C. Spectrotemporal analysis of evoked and induced electroencephalographic responses in primary auditory cortex (A1) of the awake monkey. Cereb. Cortex 18, 610–625 (2008).

Article PubMed Google Scholar

Yizhar, O. et al. Neocortical excitation/inhibition balance in information processing and social dysfunction. Nature 477, 171–178 (2011).

Article ADS CAS PubMed PubMed Central Google Scholar

Aabedi, A. A. & Hervey-Jumper, S. Convergence of heteromodal lexical retrieval in the lateral prefrontal cortex. Sci. Rep. 11, 6305 (2020).

Chang, E. F. et al. Cortical spatio-temporal dynamics underlying phonological target detection in humans. J. Cogn. Neurosci. 23, 1437–1446 (2011).

Article PubMed Google Scholar

Romanski, L. M. et al. Dual streams of auditory afferents target multiple domains in the primate prefrontal cortex. Nat. Neurosci. 2, 1131–1136 (1999).

Article CAS PubMed PubMed Central Google Scholar

Chartier, J., Anumanchipalli, G. K., Johnson, K. & Chang, E. F. Encoding of articulatory kinematic trajectories in human speech sensorimotor cortex. Neuron 98, 1042–1054 (2018).

Article CAS PubMed PubMed Central Google Scholar

Criss, A. H., Aue, W. R. & Smith, L. The effects of word frequency and context variability in cued recall. J. Mem. Lang. 64, 119–132 (2011).

Article Google Scholar

Bhaduri, A. et al. Outer radial glia-like cancer stem cells contribute to heterogeneity of glioblastoma. Cell Stem Cell 26, 48–63 (2020).

Article CAS PubMed PubMed Central Google Scholar

Neftel, C. et al. An integrative model of cellular states, plasticity, and genetics for glioblastoma. Cell 178, 835–849 (2019).

Article CAS PubMed PubMed Central Google Scholar

Tirosh, I. & Suva, M. L. Tackling the many facets of glioblastoma heterogeneity. Cell Stem Cell 26, 303–304 (2020).

Article CAS PubMed Google Scholar

Yu, K. et al. PIK3CA variants selectively initiate brain hyperactivity during gliomagenesis. Nature 578, 166–171 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Guggisberg, A. G. et al. Mapping functional connectivity in patients with brain lesions. Ann. Neurol. 63, 193–203 (2008).

Article PubMed PubMed Central Google Scholar

Lee, A. T. et al. The impact of high functional connectivity network hub resection on language task performance in adult low- and high-grade glioma. J. Neurosurg. https://doi.org/10.3171/2020.1.JNS192267 (2020).

Nolte, G. et al. Identifying true brain interaction from EEG data using the imaginary part of coherency. Clin. Neurophysiol. 115, 2292–2307 (2004).

Article PubMed Google Scholar

Brookes, M. J. et al. Measuring functional connectivity using MEG: methodology and comparison with fcMRI. Neuroimage 56, 1082–1104 (2011).

Article PubMed Google Scholar

Tarapore, P. E. et al. Magnetoencephalographic imaging of resting-state functional connectivity predicts postsurgical neurological outcome in brain gliomas. Neurosurgery 71, 1012–1022 (2012).

Article PubMed Google Scholar

Tarapore, P. E. et al. Preoperative multimodal motor mapping: a comparison of magnetoencephalography imaging, navigated transcranial magnetic stimulation, and direct cortical stimulation. J. Neurosurg. 117, 354–362 (2012).

Article PubMed PubMed Central Google Scholar

Christopherson, K. S. et al. Thrombospondins are astrocyte-secreted proteins that promote CNS synaptogenesis. Cell 120, 421–433 (2005).

Article CAS PubMed Google Scholar

Daubon, T. et al. Deciphering the complex role of thrombospondin-1 in glioblastoma development. Nat. Commun. 10, 1146 (2019).

Article ADS PubMed PubMed Central Google Scholar

Jeanne, A., Schneider, C., Martiny, L. & Dedieu, S. Original insights on thrombospondin-1-related antireceptor strategies in cancer. Front. Pharmacol. 6, 252 (2015).

Article PubMed PubMed Central Google Scholar

Allen, N. J. & Eroglu, C. Cell biology of astrocyte-synapse interactions. Neuron 96, 697–708 (2017).

Article CAS PubMed PubMed Central Google Scholar

Bushong, E. A., Martone, M. E., Jones, Y. Z. & Ellisman, M. H. Protoplasmic astrocytes in CA1 stratum radiatum occupy separate anatomical domains. J. Neurosci. 22, 183–192 (2002).

Article CAS PubMed PubMed Central Google Scholar

Oberheim, N. A. et al. Uniquely hominid features of adult human astrocytes. J. Neurosci. 29, 3276–3287 (2009).

Article CAS PubMed PubMed Central Google Scholar

Fernandopulle, M. S. et al. Transcription factor-mediated differentiation of human iPSCs into neurons. Curr. Protoc. Cell Biol. 79, e51 (2018).

Article PubMed PubMed Central Google Scholar

Joseph, J. V. et al. TGF-β promotes microtube formation in glioblastoma through thrombospondin 1. Neuro Oncol. 24, 541–553 (2022).

Article CAS PubMed Google Scholar

Jung, E. et al. Tweety-homolog 1 drives brain colonization of gliomas. J. Neurosci. 37, 6837–6850 (2017).

Article CAS PubMed PubMed Central Google Scholar

Osswald, M. et al. Brain tumour cells interconnect to a functional and resistant network. Nature 528, 93–98 (2015).

Article ADS CAS PubMed Google Scholar

Molinaro, A. M. et al. Association of maximal extent of resection of contrast-enhanced and non-contrast-enhanced tumor with survival within molecular subgroups of patients with newly diagnosed glioblastoma. JAMA Oncol. 6, 495–503 (2020).

Article PubMed PubMed Central Google Scholar

Lostritto, K., Strawderman, R. L. & Molinaro, A. M. A partitioning deletion/substitution/addition algorithm for creating survival risk groups. Biometrics 68, 1146–1156 (2012).

Article MathSciNet PubMed MATH Google Scholar

Molinaro, A. M., Lostritto, K. & van der Laan, M. partDSA: deletion/substitution/addition algorithm for partitioning the covariate space in prediction. Bioinformatics 26, 1357–1363 (2010).

Article CAS PubMed PubMed Central Google Scholar

Bell, B. D., Seidenberg, M., Hermann, B. P. & Douville, K. Visual and auditory naming in patients with left or bilateral temporal lobe epilepsy. Epilepsy Res. 55, 29–37 (2003).

Article PubMed Google Scholar

Hamberger, M. J., McClelland, S. 3rd, McKhann, G. M. 2nd, Williams, A. C. & Goodman, R. R. Distribution of auditory and visual naming sites in nonlesional temporal lobe epilepsy patients and patients with space-occupying temporal lobe lesions. Epilepsia 48, 531–538 (2007).

Article PubMed Google Scholar

Eroglu, C. et al. Gabapentin receptor α2δ-1 is a neuronal thrombospondin receptor responsible for excitatory CNS synaptogenesis. Cell 139, 380–392 (2009).

Article CAS PubMed PubMed Central Google Scholar

Pan, Y. et al. NF1 mutation drives neuronal activity-dependent initiation of optic glioma. Nature 594, 277–282 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Belgers, V. et al. Postoperative oscillatory brain activity as an add-on prognostic marker in diffuse glioma. J. Neurooncol. 147, 49–58 (2020).

Article PubMed PubMed Central Google Scholar

Daniel, A. G. S. et al. Functional connectivity within glioblastoma impacts overall survival. Neuro. Oncol. 23, 412–421 (2020).

Derks, J. et al. Oscillatory brain activity associates with neuroligin-3 expression and predicts progression free survival in patients with diffuse glioma. J. Neurooncol. 140, 403–412 (2018).

Article CAS PubMed PubMed Central Google Scholar

Venkataramani, V. et al. Glioblastoma hijacks neuronal mechanisms for brain invasion. Cell 185, 2899–2917 (2022).

Article CAS PubMed Google Scholar

Delorme, A. & Makeig, S. EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 134, 9–21 (2004).

Article PubMed Google Scholar

Dalal, S. S. et al. Five-dimensional neuroimaging: localization of the time-frequency dynamics of cortical activity. Neuroimage 40, 1686–1700 (2008).

Article PubMed Google Scholar

Aabedi, A. A. et al. Assessment of wakefulness during awake craniotomy to predict intraoperative language performance. J. Neurosurg. https://doi.org/10.3171/2019.2.JNS183486 (2019).

Chang, E. F., Raygor, K. P. & Berger, M. S. Contemporary model of language organization: an overview for neurosurgeons. J. Neurosurg. 122, 250–261 (2015).

Article PubMed Google Scholar

Rolls, E. T., Huang, C. C., Lin, C. P., Feng, J. & Joliot, M. Automated Anatomical Labelling Atlas 3. Neuroimage 206, 116189 (2020).

Article PubMed Google Scholar

Hinkley, L. B. et al. Cognitive impairments in schizophrenia as assessed through activation and connectivity measures of magnetoencephalography (MEG) data. Front. Hum. Neurosci. 3, 73 (2010).

PubMed PubMed Central Google Scholar

Wipf, D. & Nagarajan, S. A unified Bayesian framework for MEG/EEG source imaging. Neuroimage 44, 947–966 (2009).

Article PubMed Google Scholar

Vrba, J. & Robinson, S. E. Signal processing in magnetoencephalography. Methods 25, 249–271 (2001).

Article CAS PubMed Google Scholar

Sekihara, K., Hild, K. E. 2nd & Nagarajan, S. S. A novel adaptive beamformer for MEG source reconstruction effective when large background brain activities exist. IEEE Trans. Biomed. Eng. 53, 1755–1764 (2006).

Article PubMed PubMed Central Google Scholar

Donoghue, T. et al. Parameterizing neural power spectra into periodic and aperiodic components. Nat. Neurosci. 23, 1655–1665 (2020).

Article CAS PubMed PubMed Central Google Scholar

Emerton, B. C., Gansler, D. A., Sandberg, E. H. & Jerram, M. Functional anatomic dissociation of description and picture naming in the left temporal lobe. Brain Imaging Behav. 8, 570–578 (2014).

Article PubMed Google Scholar

Hamberger, M. J., Seidel, W. T., Goodman, R. R., Perrine, K. & McKhann, G. M. Temporal lobe stimulation reveals anatomic distinction between auditory naming processes. Neurology 60, 1478–1483 (2003).

Article CAS PubMed Google Scholar

Wilson, S. M., Eriksson, D. K., Schneck, S. M. & Lucanie, J. M. A quick aphasia battery for efficient, reliable, and multidimensional assessment of language function. PLoS ONE 13, e0192773 (2018).

Article PubMed PubMed Central Google Scholar

Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).

Article CAS PubMed Google Scholar

Schneider, V. A. et al. Evaluation of GRCh38 and de novo haploid genome assemblies demonstrates the enduring quality of the reference assembly. Genome Res. 27, 849–864 (2017).

Article CAS PubMed PubMed Central Google Scholar

Kim, D., Paggi, J. M., Park, C., Bennett, C. & Salzberg, S. L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915 (2019).

Article CAS PubMed PubMed Central Google Scholar

Liao, Y., Smyth, G. K. & Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930 (2014).

Article CAS PubMed Google Scholar

Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).

Article PubMed PubMed Central Google Scholar

Zhu, A., Ibrahim, J. G. & Love, M. I. Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics 35, 2084–2092 (2019).

Article CAS PubMed Google Scholar

Bakken, T. E. et al. Single-nucleus and single-cell transcriptomes compared in matched cortical cell types. PLoS ONE 13, e0209648 (2018).

Article PubMed PubMed Central Google Scholar

Couturier, C. P. et al. Single-cell RNA-seq reveals that glioblastoma recapitulates a normal neurodevelopmental hierarchy. Nat. Commun. 11, 3406 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Butler, A., Hoffman, P., Smibert, P., Papalexi, E. & Satija, R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat. Biotechnol. 36, 411–420 (2018).

Article CAS PubMed PubMed Central Google Scholar

Stuart, T. et al. Comprehensive integration of single-cell data. Cell 177, 1888–1902 (2019).

Article CAS PubMed PubMed Central Google Scholar

Hafemeister, C. & Satija, R. Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. Genome Biol. 20, 296 (2019).

Article CAS PubMed PubMed Central Google Scholar

Korsunsky, I. et al. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat. Methods 16, 1289–1296 (2019).

Article CAS PubMed PubMed Central Google Scholar

Krencik, R. et al. Systematic three-dimensional coculture rapidly recapitulates interactions between human neurons and astrocytes. Stem Cell Rep. 9, 1745–1753 (2017).

Article CAS Google Scholar

Eisenman, L. N., Emnett, C. M., Mohan, J., Zorumski, C. F. & Mennerick, S. Quantification of bursting and synchrony in cultured hippocampal neurons. J. Neurophysiol. 114, 1059–1071 (2015).

Article CAS PubMed PubMed Central Google Scholar

Eggermont, J. J. Properties of correlated neural activity clusters in cat auditory cortex resemble those of neural assemblies. J. Neurophysiol. 96, 746–764 (2006).

Article PubMed Google Scholar

Passaro, A. P., Aydin, O., Saif, M. T. A. & Stice, S. L. Development of an objective index, neural activity score (NAS), reveals neural network ontogeny and treatment effects on microelectrode arrays. Sci. Rep. 11, 9110 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Pastore, V. P., Massobrio, P., Godjoski, A. & Martinoia, S. Identification of excitatory-inhibitory links and network topology in large-scale neuronal assemblies from multi-electrode recordings. PLoS Comput. Biol. 14, e1006381 (2018).

Article ADS PubMed PubMed Central Google Scholar

Ullo, S. et al. Functional connectivity estimation over large networks at cellular resolution based on electrophysiological recordings and structural prior. Front. Neuroanat. 8, 137 (2014).

Article PubMed PubMed Central Google Scholar

Cheng, J. K. & Chiou, L. C. Mechanisms of the antinociceptive action of gabapentin. J. Pharmacol. Sci. 100, 471–486 (2006).

Article ADS CAS PubMed Google Scholar

Roberto, M. et al. Cellular and behavioral interactions of gabapentin with alcohol dependence. J. Neurosci. 28, 5762–5771 (2008).

Article CAS PubMed PubMed Central Google Scholar

Keough, M. B. et al. An inhibitor of chondroitin sulfate proteoglycan synthesis promotes central nervous system remyelination. Nat. Commun. 7, 11312 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Cheng, C., Lau, S. K. & Doering, L. C. Astrocyte-secreted thrombospondin-1 modulates synapse and spine defects in the fragile X mouse model. Mol. Brain 9, 74 (2016).

Article PubMed PubMed Central Google Scholar

Crawford, D. C., Jiang, X., Taylor, A. & Mennerick, S. Astrocyte-derived thrombospondins mediate the development of hippocampal presynaptic plasticity in vitro. J. Neurosci. 32, 13100–13110 (2012).

Article CAS PubMed PubMed Central Google Scholar

Download references

This study was supported by NIH grants K08NS110919 and P50CA097257; Robert Wood Johnson Foundation grant 74259; the UCSF LoGlio Collective and Resonance Philanthropies; and U19 CA264339, Tom Paquin Brain Cancer Research Fund to S.L.H.-J.; the Sullivan Brain Cancer Fund to S. Krishna; NIH grants F30CA246808 and T32GM007618 to A.C.; NIH grant K99CA252001 to H.S.V.; NIH grant R01NS100440 to S.N.; NIH grant R00DC013828 to D.B.; NIH grants R01NS092597, DP1NS111132, P50CA165962, R01CA258384, R01CA263500 and U19CA264504, a Robert J. Kleberg Jr and Helen C. Kleberg Foundation grant and funding from Cancer Research UK to M.M.; American Brain Tumour Association grant MSSF1900021 and Lucien Rubinstein Award to N.A.; and the UCSF Physician Scientist Scholar Program, the UCSF Wolfe Meningioma Program Project, NIH grant K08CA212279 to D.R.R. We thank J. Phillips, A. Shai and the staff of the UCSF Brain Tumour Center Biorepository and Pathology Core; the staff within the Biological Imaging Development CoLab (BIDC) at UCSF for confocal microscopy and Imaris software support; T. Ozawa and R. Santos of the UCSF Brain Tumour Center Preclinical Therapeutics Core; N. Jumper for proofreading; and K. Probst for illustrations.

Department of Neurological Surgery, University of California, San Francisco, San Francisco, CA, USA

Saritha Krishna, Abrar Choudhury, Kyounghee Seo, Sofia Kakaizada, Anthony Lee, Alexander Aabedi, Galina Popova, Caroline Cao, Cesar Nava Gonzales, Rasika Sudharshan, Andrew Egladyous, Nyle Almeida, Yalan Zhang, Annette M. Molinaro, Andy G. S. Daniel, Jeanette Hyer, Edward F. Chang, Joanna J. Phillips, David R. Raleigh & Shawn L. Hervey-Jumper

Department of Neurology, Stanford University, Stanford, CA, USA

Michael B. Keough, Lijun Ni, Humsa S. Venkatesh, Kiarash Shamardani & Michelle Monje

Department of Psychology, University of Michigan, Ann Arbor, MI, USA

Benjamin Lipkin & David Brang

Weill Institute for Neurosciences, University of California San Francisco, San Francisco, CA, USA

Edward F. Chang & Shawn L. Hervey-Jumper

Department of Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, USA

Anne Findlay & Srikantan Nagarajan

Department of Pathology, University of California, San Francisco, San Francisco, CA, USA

Joanna J. Phillips & David R. Raleigh

Department of Radiation Oncology, University of California, San Francisco, San Francisco, USA

David R. Raleigh

Howard Hughes Medical Institute, Stanford, CA, USA

Michelle Monje

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

S. Krishna and S.L.H.-J. designed, conducted and analysed experiments. M.M. contributed to study design, data analysis and manuscript editing. A.C. and D.R. contributed to bulk and single cell transcriptomic analyses. M.B.K. and K. Shamardani contributed to in vivo PDX gabapentin experiment. K. Seo contributed to in vitro organoid experiments. L.N. contributed to electron microscopy data acquisition and analyses. S. Kakaizada and N.A. contributed to language and cognitive assessments. Y.Z., A.M.M., A.L. and A.A. contributed to human survival analyses. G.P. and J.H. contributed to the MEA experiment. J.J.P., C.C., C.G. and R.S. contributed to quantitative imaging analysis. A.E., A.F. and S.N. contributed to MEG data acquisition and quantification. H.S.V. contributed to orthotopic xenografting. E.F.C., B.L., A.G.S.D. and D.B. contributed to intraoperative electrocorticography analysis. S. Krishna and S.L.H.-J. wrote the manuscript. S.L.H.-J. conceived the project and supervised all aspects of the work.

Correspondence to Shawn L. Hervey-Jumper.

The authors declare no competing interests.

Nature thanks Michael Taylor, Frank Winkler and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

a, Schematic of study workflow. In human participants with dominant hemisphere gliomas, we applied subcortical high density electrocorticography during audiovisual speech initiation to assess tumour intrinsic neuronal circuit dynamics. Focusing on glioblastoma, we then assessed long-range functional connectivity using magnetoencephalography (MEG) imaginary coherence. b, Extra-operative language assessments were performed for correlation with biological assays. c, Long-range measure of tumour intrinsic functional connectivity identified regions of high and low connectivity for site specific biopsies which were used for in vivo and in vitro cell biology experiments. d, Multiple layered approach including clinical variables, cognition assessments, human and animal model network dynamics, and cell biology experiments serves as a platform for glioma influence on neural circuit dynamics.

a, Electrodes overlying normal appearing and glioma-infiltrated regions a cohort of 14 adult patients with cortically projecting glioma in the lateral prefrontal cortex. Electrodes over non-tumour regions are shown in white and those over tumour-infiltrative regions in black. b, Positive control conditions included speech initiation responses within non-infiltrated cortex of the left lateral prefrontal cortex (LPFC) for non-cortically projecting glioblastoma. (Left) Axial FLAIR MRI demonstrates tumour location within insular cortex. Hemisphere of language dominance on the left was performed according to study protocol. (Middle) Black outline illustrates LPFC with ECoG recordings obtained from electrode A24, denoted by the red dot. White star represents frontal lobe motor cortex. (Right) Identical to non-tumour comparisons for cortically projecting gliomas, speech responses demonstrate elevate high gamma power (HGp) prior to speech onset consistent with speech motor planning. Dark line and shaded region represent mean and 95% confidence interval, respectively. c, Selectivity of maintained tumour intrinsic task-specific cortical responses is identified across diffuse glioma subtypes (adapted from Aabedi et al. 2021). Spectral data show clear separation of frequencies across tumour (glioma-infiltrated) and non-tumour (normal-appearing) electrodes. Group level analysis of participants (n = 12) demonstrates speech initiation responses across WHO 2–4 diffuse glioma. d, Glioma subtype-specific speech initiation spectral responses for electrodes above normal-appearing and glioma-infiltrated cortex showing conserved phenotype with task-specific hyperexcitability observed only in participants with glioblastoma. Subtypes: grade 2 and 3 oligodendroglioma (n = 4), grade 2 and 3 astrocytoma (n = 4), and glioblastoma (n = 4). P value determined by two-tailed Student's t-test and corrections for multiple comparisons were made using the FDR method (a).

Source data

a, Spectral data of time vs frequency from 100–150 Hz for all tumour and non-tumour electrodes show the expected pattern of HgP increasing above 50 Hz. b, High gamma power (HGp) recording from single electrodes overlying tumour-infiltrated regions of brain. Dark line and shaded region represent mean and 95% confidence interval, respectively. c, Reconstructed time series of HGp from non-tumour electrodes demonstrating expected spatial and temporal pattern of neural activity within lateral prefrontal cortex. d, e, Electrodes matched to anatomical areas across tumour and non-tumour demonstrate hyperexcitability with glioma-infiltrated cortex. Plots are centred on median, with dark bars indicating the first and third quartiles, and whiskers the minimum and maximum points (n = 101 per group, P = 0.0163). P value determined by two-sided linear mixed-effects model with corrections for multiple comparisons made using the FDR method (b) and two-tailed Student's t-test (e). *P < 0.05.

Source data

a, Linear regression statistics illustrating that gamma power (a measure of neuronal activity) correlates with number of intratumoural high functional connectivity voxels in glioblastoma (n = 18 patients; P = 0.00002). Shaded area represents the 95% confidence interval predicted by the linear regression model. b, Sampling of functionally connected intratumoural regions using MEG was performed exclusively in participants with dominant hemisphere glioblastoma at the point of initial diagnosis. Site-directed tissue biopsies from HFC and LFC regions were taken as determined by MRI. Table illustrates site-specific sampling of each annotated specimen as it relates to contrast enhancing (CE) region and FLAIR tumour. Site-specific samples were acquired without regard for whether they originated from enhancing or FLAIR regions. c, While samples were not acquired based on whether they originated from contrast enhancing or FLAIR regions, the stereotactic coordinates of each sample were acquired. While 57.78% of HFC samples originated from contrast enhancing regions, this did not reach statistical significance. P = 0.1923 two-sided chi square, P = 0.235 two-sided Fisher's exact test. P value determined by two-sided linear regression analysis (a). NS, not significant.

Source data

a, Bulk RNA transcriptomic profile of HFC tissues showed a neurogenic signature including elevated (7-fold) expression of thrombospondin-1 (TSP-1) (n = 3-4 per group). b, Unsupervised principal component (PCA) analysis of bulk RNA sequencing data obtained from glioblastoma primary patient HFC (n = 3) and LFC (n = 4) samples. c, Volcano plots of IDH-WT glioblastoma samples revealed 144 differentially expressed genes between HFC and LFC tumour regions. The blue dots represent all differentially expressed genes, where differential expression is defined by the parameters: adjusted p-value < 0.05 and absolute log2fold change > 1. P value calculated using two-sided Wald test and adjusted for multiple comparisons with the Benjamini-Hochberg method.

Source data

a—c, Tumour cell validation using copy number variant assessment on three matched pairs of HFC and LFC samples from FC1 (SF#1), FC2 (SF#2), and FC3 (SF#3) glioblastoma patients. Trisomy 7 and monosomy 10 co-occur in most cells in FC1.Trisomy 7 is an early event, while monosomy 10 is a late event in FC2. FC3 patient sample contains no copy number variation but has high level amplification of NTRK2 gene. d, Single-cell RNA transcriptomic profile UMAP confirms distinct cell populations including non-tumour astrocytes and neurons. e, Gene enrichment profile used to identify each of the UMAP cell populations. f, g, Feature plot for TSP-1 in combined (HFC + LFC) and LFC (n = 7,065 cells, 3 participants) population; within LFC samples, TSP-1 expression is primarily from non-tumour astrocytes (suggesting that within low connectivity intratumoural regions, normal astrocytes secrete TSP-1 to generate connectivity mirroring normal physiology). h, Dot plots showing TSP-1 expression (grey to red scale) and percentage (number of cells expressing the gene) of TSP-1-positive cells in tumour cells and non-tumour astrocyte populations in HFC and LFC samples (n = 3 per group). Out of the total HFC tumour cells (n = 5325, 3 patients), 157 cells are TSP-1 positive accounting for a percentage of 2.95, while only 1.59% (51 cells out of a total of 3212 LFC tumour cells [n = 3 patients]) express TSP-1. However, in the non-tumour astrocyte population, the number of TSP-1-positive cells are higher in LFC (n = 34 out of a total of 41 astrocytes, accounting for 82.9%) compared to HFC (n = 15 out of a total of 20 astrocytes, accounting for 75%) samples. i, Violin plots illustrating significantly increased TSP-1 expression within HFC region glioblastoma cells relative to LFC regions (P = 1.4*10−7). j, Compared to HFC, slight trend of increased TSP-1 expression within non-tumour astrocytes of LFC population. However, this trend did not reach statistical significance, likely due to the small number of non-tumour astrocytes captured (P = 0.45). P values determined by two-tailed Student's t-test.

Source data

a, TSP-1 (immunohistochemistry) expression in arbitrary units (A.U.) in HFC and LFC tissues (n = 16 regions; 3 per group; P = 0.04). Scale bar, 50 µm. b, Representative confocal images of primary patient-derived HFC and LFC tissues showing regions of synaptic puncta colocalization (white arrows). Orange, synapsin-1 (presynaptic puncta); red, PSD95 (postsynaptic puncta); green, neurofilament (neurons). Scale bar, 15 µm. Quantification of the number of colocalized pre- and postsynaptic puncta (HFC: n = 7; LFC: n = 10 regions, 2 per group; P = 0.05). Quantification of postsynaptic PSD95 puncta size (HFC: n = 6; LFC: n = 9, 2 per group; P = 0.0441). Scale bar, 15 µm. c, Primary patient-derived cultures and mouse hippocampal neuron controls. Neurofilament (heavy and medium chains) and nestin antibodies used as specific markers to label mouse hippocampal neurons and glioblastoma cells, respectively in glioma-neuron co-culture. Left panel: mouse hippocampal neurons alone in culture for 14 days only express neurofilament (green) and not nestin (orange). Right panel: Nestin (orange) expression in GBM cells co-cultured with neurofilament (green) labelled mouse hippocampal neurons for 14 days. Scale bar, 100 µm. d, Cell cultures tested for mycoplasma using a commercially available kit (PCR Mycoplasma Test Kit I/C, PromoCell, Heidelberg, Germany) shows absence of a positive band at ~270 bp. Tested primary patient-derived lines shows internal control DNA at ~479 bp indicated a successfully performed PCR. e, Neuron organoids (GFP labelled) were generated from an iPSC cell line integrated with doxycycline inducible human NGN2 transgene and co-cultured with RFP labelled HFC and LFC cells (pseudo-coloured white) for two weeks. Quantification of postsynaptic Homer-1 puncta density (calculated by dividing the number of puncta measured with the area of the image field) in 2-week induced neuron (iN) organoid sections (n = 8 per group; P = 0.0009). Scale bar, 10 µm. f, Multi-electrode array of glioma-neuron co-culture and control conditions. Magnified view of multi-electrode array (MEA), showing RFP-labelled glioblastoma cells in co-culture with neurons (top row). Scale bar, 100 µm. Representative raster plot showing individual spikes/extracellular action potentials (tick mark), bursts (cluster of spikes in blue) and synchronized network bursts (pink) of mouse cortical neuron (DIV 18) only condition (bottom row). The cumulative trace above the raster plots depicts the population spike time histogram indicating the synchronized activity between the different electrodes (network burst). g, Structural synapses in primary patient-derived glioblastoma xenografts. Quantification of neuron-to-neuron (P = 0.1381) and neuron-glioma synapses (P = 0.0005) per high power field (hpf) in HFC and LFC xenografts (top row). Specificity negative controls for immuno-gold labelling (bottom row). (Left) HFC xenograft with secondary antibody only (no primary antibody) control and (right) non-glioma bearing negative control tissue demonstrating few randomly distributed immunogold particles across the tissue specimen. Scale bar, 1000 nm. Data presented as mean ± s.e.m. P values determined by two-tailed Student's t-test. *P < 0.05; **P < 0.01; ***P < 0.001; NS, not significant.

Source data

Primary patient glioblastoma cells from HFC regions illustrate marked increase in proliferation when co-cultured with mouse hippocampal neurons. a, Quantification of proliferation indices of HFC (n = 4) and LFC (n = 4) glioma cells alone in culture (in the absence of neurons) from individual patient lines, determined by quantifying the fraction of EdU labelled cells/DAPI labelled cells; P = 0.00005). b, Representative confocal images illustrating proliferating HFC and LFC glioma cells (EdU+, green) in the absence or presence of mouse hippocampal neurons (72h co-culture). Scale bar, 100 µm. Data presented as mean ± s.e.m. P value determined by two-tailed Student's t-test. ****P < 0.0001.

Source data

a, 3D spheroid invasion assay showing representative micrographs imaged 24 h after addition of invasion matrix. Analysis includes quantification of mean spheroid invasion area normalized for each sample to the initial (0 h) spheroid area (HFC: n = 11; HFC + NCM: n = 16; LFC: n = 5; LFC + NCM: n = 6 spheroids, 1 per group; P = 0.02). Scale bar, 200 µm. Data presented as mean ± s.e.m. P values determined by two-tailed Student's t-test. b, Representative confocal images of primary patient-derived HFC and LFC tissues showing MET-positive glioma cells (white arrows). Red, MET; orange, Nestin (HFC/LFC-GBM cells); blue, DAPI. Scale bar, 30 µm. Quantification of MET-positive glioma cells per high-power field (HFC: n = 8; LFC: n = 11, 3 per group; P = 0.0005). c, Representative immunohistochemistry images of MET staining in HFC and LFC tissues demonstrate increased tissue level protein expression (HFC: n = 26; LFC: n = 24 regions, 4 per group; P = 0.0329). Scale bar, 50 µm. Data presented as mean ± s.e.m. P values determined by two-tailed Student's t-test. d, Representative confocal images showing the diffuse infiltrative pattern of HFC cells in the hippocampus in comparison to the LFC cells. Quantification of tumour burden of HFC and LFC hippocampal xenografts using rank order analysis (HFC: n = 13; LFC: n = 11 regions; P = 0.002). Scale bar, 100 µm. Data presented as mean ± s.e.m. P value determined by two-tailed Mann-Whitney test. *P < 0.05; **P < 0.01; ***P < 0.001; NS, not significant.

Source data

a, Cell viability determined by live/dead cell assay. Representative images illustrating no significant cell death 2 weeks post-transduction of HFC cells with control (shCont) or TSP-1 shRNAs. Live (green) and dead (red) cells were imaged from four random fields per well and were visualized under a fluorescence microscope. The percentage of live cells were calculated as the number of live cells (in green) divided by the total number of cells (green + red) per image field (n = 4 regions per group; P = 0.0720). Scale bar, 100 µm. b, ELISA experiments performed on cell culture supernatants demonstrating strong reduction of TSP-1 expression after knockdown of TSP-1 in two different primary patient-derived HFC cell lines compared to the control scramble condition (n = 2 per group). Data presented as mean ± s.e.m. P values determined by two-tailed Student's t-test. N.S., not significant.

Source data

a, Kaplan-Meier human survival analysis illustrates 71-week overall survival for patients with HFC voxels as determined by contrast-enhanced T1-weighted images as compared to 123-weeks for participants without HFC voxels within their glioblastoma (mean follow-up months 50.5, range 4.9–155.9 months). b, Picture and auditory naming language task performance across the study population. c, Linear regression statistics illustrating a negative correlation between the number of intratumoural high functional connectivity voxels and baseline auditory and picture naming scores (n = 31, P = 0.0181). P value determined by two-tailed linear regression analysis.

Source data

Representative confocal images from HFC glioma monoculture showing no significant change in proliferation (as measured by the number of human nuclear antigen (HNA)-positive cells co-labelled with Ki67 divided by the total number of HNA-positive tumour cells counted across all areas quantified) upon pharmacological TSP-1 inhibition using (32 µM) gabapentin (HFC: n = 7; LFC: n = 9 regions, 2 per group; P = 0.50). Red, HNA (human nuclei); white, Ki67. Scale bar, 30 µm. Data presented as mean ± s.e.m. P values determined by two-tailed Student's t-test. NS, not significant.

Source data

Supplementary Tables 1–5.

3D co-culture of HFC glioma cells with induced neuron organoids.

3D co-culture of LFC glioma cells with induced neuron organoids.

3D co-culture of LFC glioma cells with induced neuron organoids in the presence of TSP-1.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Krishna, S., Choudhury, A., Keough, M.B. et al. Glioblastoma remodelling of human neural circuits decreases survival. Nature 617, 599–607 (2023). https://doi.org/10.1038/s41586-023-06036-1

Download citation

Received: 18 February 2021

Accepted: 31 March 2023

Published: 03 May 2023

Issue Date: 18 May 2023

DOI: https://doi.org/10.1038/s41586-023-06036-1

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.