One challenge in contemporary neuroscience is to achieve an integrated understanding of the large-scale brain-wide interactions, particularly the spatiotemporal patterns of neural activity that give rise to functions and behavior. At present, little is known about the spatiotemporal properties of long-range neuronal networks. We examined brain-wide neural activity patterns elicited by stimulating ventral posteromedial (VPM) thalamo-cortical excitatory neurons through combined optogenetic stimulation and functional MRI (fMRI). We detected robust optogenetically evoked fMRI activation bilaterally in primary visual, somatosensory, and auditory cortices at low (1 Hz) but not high frequencies (5-40 Hz). Subsequent electrophysiological recordings indicated interactions over long temporal windows across thalamo-cortical, cortico-cortical, and interhemispheric callosal projections at low frequencies. We further observed enhanced visually evoked fMRI activation during and after VPM stimulation in the superior colliculus, indicating that visual processing was subcortically modulated by low-frequency activity originating from VPM. Stimulating posteromedial complex thalamo-cortical excitatory neurons also evoked brain-wide bloodoxygenation-level-dependent activation, although with a distinct spatiotemporal profile. Our results directly demonstrate that lowfrequency activity governs large-scale, brain-wide connectivity and interactions through long-range excitatory projections to coordinate the functional integration of remote brain regions. This low-frequency phenomenon contributes to the neural basis of long-range functional connectivity as measured by resting-state fMRI.fMRI | optogenetic | brain connectivity | low frequency | thalamus T he brain is a highly complex, interconnected structure with parallel and hierarchical networks distributed within and between neural systems (1, 2). This integrative architecture dictates the underlying principles of how brain-wide neural connectivity supports and organizes sensory, behavioral, and cognitive processes (3). Recent advances in structural (2, 4) and functional connectivity (5-12) mapping, as well as neural circuit modulatory tools, such as optogenetics (13, 14), permit detailed, high-resolution neural examinations at unprecedented scales. In particular, functional connectivity mapping using resting-state functional MRI (rsfMRI) produces noninvasive visualization of slow and spontaneous hemodynamic fluctuations in humans (5-9) and animals (10-12). Coherent low-frequency (<1 Hz) fluctuations across functionally coupled, large-scale networks is an intriguing property, although the exact underlying neural bases and functional significance remain unclear. Previous studies integrating large-scale electrical recordings, voltage-sensitive dye (VSD), and Ca 2+ -imaging techniques (15-18) support the hypothesis that slow, oscillating neural activity constrains and elicits these hemodynamic fluctuations. Furthermore, low-frequency activity can temporally synchronize remote brain region...