Simulating laminar neuroimaging data for a visual delayed match-to-sample task

Invasive electrophysiological and neuroanatomical studies in nonhuman mammalian experimental preparations have helped elucidate the lamina (layer) dependence of neural computations and interregional connections. Noninvasive functional neuroimaging can, in principle, resolve cortical laminae (layers), and thus provide insight into human neural computations and interregional connections. However human neuroimaging data are noisy and difficult to interpret; biologically realistic simulations can aid experimental interpretation by relating the neuroimaging data to simulated neural activity. We illustrate the potential of laminar neuroimaging by upgrading an existing large-scale, multiregion neural model that simulates a visual delayed match-to-sample (DMS) task. The new laminar-based neural unit incorporates spiny stellate, pyramidal, and inhibitory neural populations which are divided among supragranular, granular, and infragranular laminae (layers). We simulated neural activity which is translated into local field potential-like data used to simulate conventional and laminar fMRI activity. We implemented the laminar connectivity schemes proposed by Felleman and Van Essen (Cerebral Cortex, 1991) for interregional connections. The hemodynamic model that we employ is a modified version of one due to Heinzle et al. (Neuroimage, 2016) that incorporates the effects of draining veins. We show that the laminar version of the model replicates the findings of the existing model. The laminar model shows the finer structure in fMRI activity and functional connectivity. Laminar differences in the magnitude of neural activities are a prominent finding; these are also visible in the simulated fMRI. We illustrate differences between task and control conditions in the fMRI signal, and demonstrate differences in interregional laminar functional connectivity that reflect the underlying connectivity scheme. These results indicate that multi-layer computational models can aid in interpreting layer-specific fMRI, and suggest that increased use of laminar fMRI could provide unique and fundamental insights to human neuroscience.