Abstract

To further advance functional magnetic resonance imaging (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 simulation 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 (VPL) receives ST inputs, while the posterior medial thalamic nucleus (POm) 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 VPL. 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 {micro}m. Long-range synaptic inputs in cortical areas were amplified approximately 2-fold 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.