Dissection of brain-wide resting-state and functional somatosensory circuits by fMRI with optogenetic silencing

To further advance functional MRI (fMRI)-based brain science, it is critical to dissect fMRI activity at the circuit level. To achieve this goal, we combined brain-wide fMRI with neuronal silencing in well-defined regions. Since focal inactivation suppresses excitatory output to downstream pathways, intact input and suppressed output circuits can be separated. Highly specific cerebral blood volume-weighted fMRI was performed with optogenetic stimulation of local GABAergic neurons in mouse somatosensory regions. Brain-wide spontaneous somatosensory networks were found mostly in ipsilateral cortical and subcortical areas, which differed from the bilateral homotopic connections commonly observed in resting-state fMRI data. The evoked fMRI responses to somatosensory stimulation in regions of the somatosensory network were successfully dissected, allowing the relative contributions of spinothalamic (ST), thalamocortical (TC), corticothalamic (CT), corticocortical (CC) inputs, and local intracortical circuits to be determined. The ventral posterior thalamic nucleus receives ST inputs, while the posterior medial thalamic nucleus receives CT inputs from the primary somatosensory cortex (S1) with TC inputs. The secondary somatosensory cortex (S2) receives mostly direct CC inputs from S1 and a few TC inputs from the ventral posterolateral nucleus. The TC and CC input layers in cortical regions were identified by laminar-specific fMRI responses with a full width at half maximum of <150 µm. Long-range synaptic inputs in cortical areas were amplified approximately twofold by local intracortical circuits, which is consistent with electrophysiological recordings. Overall, whole-brain fMRI with optogenetic inactivation revealed brain-wide, population-based, long-range circuits, which could complement data typically collected in conventional microscopic functional circuit studies.