Shemesh Lab

Glioblastomas are aggressive brain tumors with dismal prognosis. One of the main bottlenecks for developing more effective therapies for glioblastoma stems from their histologic and molecular heterogeneity, leading to distinct tumor microenvironments and disease phenotypes. Effectively characterizing these features would improve the clinical management of glioblastoma. Glucose flux rates through glycolysis and mitochondrial oxidation have been recently shown to quantitatively depict glioblastoma proliferation in mouse models (GL261 and CT2A tumors) using dynamic glucose-enhanced (DGE) deuterium spectroscopy. However, the spatial features of tumor microenvironment phenotypes remain hitherto unresolved. Here, we develop a DGE Deuterium Metabolic Imaging (DMI) approach for profiling tumor microenvironments through glucose conversion kinetics. Using a multimodal combination of tumor mouse models, novel strategies for spectroscopic imaging and noise attenuation, and histopathological correlations, we show that tumor lactate turnover mirrors phenotype differences between GL261 and CT2A mouse glioblastoma, whereas recycling of the peritumoral glutamate-glutamine pool is a potential marker of invasion capacity in pooled cohorts, linked to secondary brain lesions. These findings were validated by histopathological characterization of each tumor, including cell density and proliferation, peritumoral invasion and distant migration, and immune cell infiltration. Our study bodes well for precision neuro-oncology, highlighting the importance of mapping glucose flux rates to better understand the metabolic heterogeneity of glioblastoma and its links to disease phenotypes.

The visual continuity illusion involves a shift in visual perception from static to dynamic vision modes when the stimuli arrive at high temporal frequency, and is critical for recognizing objects moving in the environment. However, how this illusion is encoded across the visual pathway remains poorly understood, with disparate frequency thresholds at retinal, cortical, and behavioural levels suggesting the involvement of other brain areas. Here, we employ a multimodal approach encompassing behaviour, whole-brain functional MRI, and electrophysiological measurements, for investigating the encoding of the continuity illusion in rats. Behavioural experiments report a frequency threshold of 18±2 Hz. Functional MRI reveal that superior colliculus signals transition from positive to negative at the behaviourally-driven threshold, unlike thalamic and cortical areas. Electrophysiological recordings indicate that these transitions are underpinned by neural activation/suppression. Lesions in the primary visual cortex reveal this effect to be intrinsic to the superior colliculus (under a cortical gain effect). Our findings highlight the superior colliculus crucial involvement in encoding temporal frequency shifts, especially the change from static to dynamic vision modes.

Understanding the dynamics of stability/plasticity balances during adulthood is pivotal for learning, disease, and recovery from injury. However, the brain-wide topography of sensory remapping remains unknown. Here, using a first-of-its-kind setup for delivering patterned visual stimuli in a rodent magnetic resonance imaging (MRI) scanner, coupled with biologically inspired computational models, we noninvasively mapped brain-wide propertiesreceptive fields (RFs) and spatial frequency (SF) tuning curvesthat were insofar only available from invasive electrophysiology or optical imaging. We then tracked the RF dynamics in the chronic visual deprivation model (VDM) of plasticity and found that light exposure progressively promoted a large-scale topographic remapping in adult rats. Upon light exposure, the initially unspecialized visual pathway progressively evidenced sharpened RFs (smaller and more spatially selective) and enhanced SF tuning curves. Our findings reveal that visual experience following VDM reshapes both structure and function of the visual system and shifts the stability/plasticity balance in adults.

Spontaneous fluctuations in functional magnetic resonance imaging (fMRI) signals correlate across distant brain areas, shaping functionally relevant intrinsic networks. However, the generative mechanism of fMRI signal correlations, and in particular the link with locally-detected ultra-slow oscillations, are not fully understood. To investigate this link, we record ultrafast ultrahigh field fMRI signals (9.4 Tesla, temporal resolution = 38 milliseconds) from female rats across three anesthesia conditions. Power at frequencies extending up to 0.3 Hz is detected consistently across rat brains and is modulated by anesthesia level. Principal component analysis reveals a repertoire of modes, in which transient oscillations organize with fixed phase relationships across distinct cortical and subcortical structures. Oscillatory modes are found to vary between conditions, resonating at faster frequencies under medetomidine sedation and reducing both in number, frequency, and duration with the addition of isoflurane. Peaking in power within clear anatomical boundaries, these oscillatory modes point to an emergent systemic property. This work provides additional insight into the origin of oscillations detected in fMRI and the organizing principles underpinning spontaneous long-range functional connectivity.

Noninvasively detecting and characterizing modulations in cellular scale micro-architecture remains a desideratum for contemporary neuroimaging. Diffusion MRI (dMRI) has become the mainstay methodology for probing microstructure, and, in ischemia, its contrasts have revolutionized stroke management. Diffusion kurtosis imaging (DKI) has been shown to significantly enhance the sensitivity of stroke detection compared to its diffusion tensor imaging (DTI) counterparts. However, the interpretation of DKI remains ambiguous as its contrast may arise from competing kurtosis sources related to the anisotropy of tissue components, diffusivity variance across components, and microscopic kurtosis (e.g., arising from cross-sectional variance, structural disorder, and restriction). Resolving these sources may be fundamental for developing more specific imaging techniques for stroke management, prognosis, and understanding its pathophysiology. In this study, we apply Correlation Tensor MRI (CTI) a double diffusion encoding (DDE) methodology recently introduced for deciphering kurtosis sources based on the unique information captured in DDE's diffusion correlation tensors to investigate the underpinnings of kurtosis measurements in acute ischemic lesions. Simulations for the different kurtosis sources revealed specific signatures for cross-sectional variance (representing neurite beading), edema, and cell swelling. Ex vivo CTI experiments at 16.4 T were then performed in an experimental photothrombotic stroke model 3 h post-stroke (N = 10), and successfully separated anisotropic, isotropic, and microscopic non-Gaussian diffusion sources in the ischemic lesions. Each of these kurtosis sources provided unique contrasts in the stroked area. Particularly, microscopic kurtosis was shown to be a primary \u201cdriver\u201d of total kurtosis upon ischemia; its large increases, coupled with decreases in anisotropic kurtosis, are consistent with the expected elevation in cross-sectional variance, likely linked to beading effects in small objects such as neurites. In vivo experiments at 9.4 T at the same time point (3 h post ischemia, N = 5) demonstrated the stability and relevance of the findings and showed that fixation is not a dominant confounder in our findings. In future studies, the different CTI contrasts may be useful to address current limitations of stroke imaging, e.g., penumbra characterization, distinguishing lesion progression form tissue recovery, and elucidating pathophysiological correlates.

Axon caliber plays a crucial role in determining conduction velocity and, consequently, in the timing and synchronization of neural activation. Noninvasive measurement of axon radii could have significant impact on the understanding of healthy and diseased neural processes. Until now, accurate axon radius mapping has eluded in vivo neuroimaging, mainly due to a lack of sensitivity of the MRI signal to micron-sized axons. Here, we show how - when confounding factors such as extra-axonal water and axonal orientation dispersion are eliminated - heavily diffusion-weighted MRI signals become sensitive to axon radii. However, diffusion MRI is only capable of estimating a single metric, the effective radius, representing the entire axon radius distribution within a voxel that emphasizes the larger axons. Our findings, both in rodents and humans, enable noninvasive mapping of critical information on axon radii, as well as resolve the long-standing debate on whether axon radii can be quantified.