The research, led by Zijuan Yu, Haoyu Wang, Yan Zhuo, and Zhentao Zuo, combined several complementary approaches in dexmedetomidine-sedated rats: functional magnetic resonance imaging (fMRI) during electrical forepaw stimulation, optogenetic activation of specific brain regions, a tuned-burst stimulation paradigm, and dynamic causal modeling of effective connectivity. The goal was to determine where, along the pathway from thalamus to cortex, temporally dense sensory input becomes constrained.

In the first experiment, twelve rats received bilateral forepaw electrical stimulation at 1, 4, 8, 25, and 55 Hz, delivered as 2-milliampere, 300-microsecond pulses in a block design with 15 seconds of stimulation followed by 45 seconds of rest. The team measured blood-oxygen-level-dependent (BOLD) responses in the primary somatosensory forelimb cortex (S1FL), the secondary somatosensory cortex (S2), the ventral posterolateral thalamic nucleus (VPL), and the posterior thalamic nucleus (Po).

The results revealed a striking dissociation between cortical and thalamic frequency profiles. In S1FL, peak BOLD responses were maximal at intermediate frequencies, reaching 15.4 percent at 4 Hz and 17.5 percent at 8 Hz, but dropped to just 2.4 percent at 25 Hz before a secondary rise to 8.8 percent at 55 Hz. Crucially, the 55-Hz peak was not accompanied by the sustained response seen at 4 to 8 Hz; the response duration, quantified as full width at half maximum, was significantly broader at 4 and 8 Hz than at any other frequency tested. In contrast, VPL and Po showed high-frequency-dominant patterns, with VPL peaking at 55 Hz at 5.4 percent and Po at 7.4 percent, and response duration in VPL did not vary significantly across frequencies.

To quantify how much thalamic activity was expressed as sustained cortical activity, the researchers calculated a sustained cortical expression index, dividing the late S1FL response by the peak VPL response. This index varied significantly across frequencies, with much higher values at 4 Hz and 8 Hz than at 25 Hz or 55 Hz. Latency analysis reinforced the ascending sequence: VPL consistently responded earlier than S1FL, with significantly shorter 5-percent-to-peak latencies at 1, 4, 8, and 55 Hz, confirming that thalamic responses preceded cortical ones across most stimulation frequencies.

The key question was whether the attenuated high-frequency cortical response reflected a local cortical limitation or a transfer failure. To answer it, the team bypassed the thalamus entirely. In one cohort of six rats, they injected an adeno-associated virus carrying channelrhodopsin under the CaMKIIα promoter into S1FL and delivered light pulses directly to the cortex at 1, 8, and 55 Hz. Direct cortical activation produced robust responses at all frequencies, with peaks of 3.87, 4.47, and 3.52 percent respectively, and the overall frequency effect was not significant. Although the 55-Hz response was modestly smaller than the 8-Hz response, it remained temporally sustained throughout the stimulation period, demonstrating that cortex retains the capacity to generate sustained high-frequency BOLD responses when activated directly.

In a separate cohort of six rats, the researchers instead stimulated VPL optogenetically and measured both local thalamic and downstream cortical responses. Local VPL responses increased with frequency, peaking at 2.24 percent at 55 Hz, significantly larger than at 1 and 8 Hz. The downstream S1FL response showed the opposite pattern: it peaked at 9.75 percent during 8-Hz stimulation but fell to just 1.96 percent at 55 Hz, a level statistically indistinguishable from the 1-Hz response. At 55 Hz, the cortical response was weak and dominated by its onset, while the thalamic response was strong. This dissociation places the main constraint downstream of preserved local VPL activation and upstream of sustained S1FL expression.

A tuned-burst paradigm then probed how temporally dense thalamic activity is reformatted before reaching cortex. The researchers embedded five-pulse 50-Hz bursts within slower carrier cycles of 1, 4, or 8 Hz, holding the intraburst frequency constant while varying the envelope. Local VPL responses were present across all carrier conditions with no significant carrier-frequency effect on peak or mean amplitude, yet S1FL responses were consistently dominated by an early transient peak followed by substantial decay. The initial cortical peak significantly exceeded the late stimulation-period response at every carrier frequency after Holm correction. The team describes this pattern as temporal compression, noting that the preferential cortical enhancement at 8 Hz seen during peripheral stimulation was not reproduced by uniform optogenetic burst activation, suggesting the intermediate-frequency preference may depend on spatiotemporal input patterns not reproduced by uniform thalamic activation.

Dynamic causal modeling of the peripheral stimulation data provided convergent network-level evidence. The intrinsic model identified a prominent baseline VPL-to-S1FL connection, and frequency-specific modulation showed that 8-Hz stimulation selectively enhanced VPL-to-S1FL effective connectivity with a posterior probability above 0.99, while 55-Hz stimulation produced no significant positive modulation and 1-Hz stimulation negatively modulated Po-to-S1FL coupling. The authors caution that because DCM was estimated from BOLD time series, the null effect at 55 Hz should be interpreted cautiously, as reduced signal-to-noise at high frequency may limit sensitivity; they regard the modeling as convergent support rather than standalone proof of reduced synaptic efficacy.

The findings carry implications for interpreting sensory-evoked fMRI more broadly. A reduced cortical BOLD response at high stimulation frequency does not necessarily indicate that local cortical tissue cannot generate a hemodynamic response; it may instead reflect reduced efficacy of pathway-level transfer from thalamus to cortex. The authors note that this could matter for understanding sensory disturbances reported in neurological disorders such as Parkinson’s disease, amyotrophic lateral sclerosis, and stroke, where an abnormal cortical response might reflect altered thalamocortical transmission rather than a purely local cortical deficit. They also emphasize limitations: the study was conducted under dexmedetomidine sedation, included only adult male rats, involved no behavioral or perceptual measurements, and localizes the constraint at the macroscopic BOLD level without identifying the cellular mechanism, which may involve short-term thalamocortical synaptic depression, feedforward inhibition, state-dependent relay dynamics, or neurovascular nonlinearities. The relationship between temporal compression and tactile frequency coding should therefore be treated as a mechanistic hypothesis rather than direct evidence about perception.

Subject of Research: Frequency-dependent limits of thalamocortical sensory transfer between the rat ventral posterolateral thalamus and primary somatosensory cortex

Article Title: Frequency-dependent thalamocortical transfer limits sustained cortical expression of sensory drive

Article References: Yu, Z., Wang, H., Zhuo, Y., & Zuo, Z. (2026). Frequency-dependent thalamocortical transfer limits sustained cortical expression of sensory drive. iScience, 29(11), Article 118035. https://doi.org/10.1016/j.isci.2026.118035

Image Credits: AI Generated

DOI: 10.1016/j.isci.2026.118035

Keywords: thalamocortical transfer, somatosensory cortex, thalamus, fMRI, BOLD response, optogenetics, dynamic causal modeling, tactile stimulation, rats, neuroscience, sensory coding, VPL