Abstract
In graph theory, "multilayer networks" represent systems involving several interconnected topological levels. One example in neuroscience is the stratification of connections between different cortical depths or "laminae", which is becoming non-invasively accessible in humans using ultra-high-resolution functional MRI (fMRI). Here, we applied multilayer graph theory to examine functional connectivity across different cortical depths in humans, using 7T fMRI (1-mm(3) voxels; 30 participants). Blood oxygenation level dependent (BOLD) signals were derived from five depths between the white matter and pial surface. We compared networks where the inter-regional connections were limited to a single cortical depth only ("layer-by-layer matrices") to those considering all possible connections between areas and cortical depths ("multilayer matrix"). We utilized global and local graph theory features that quantitatively characterize network attributes including network composition, nodal centrality, path-based measures, and hub segregation. Detecting functional differences between cortical depths was improved using multilayer connectomics compared to the layer-by-layer versions. Superficial depths of the cortex dominated information transfer and deeper depths drove clustering. These differences were largest in frontotemporal and limbic regions. fMRI functional connectivity across different cortical depths may contain neurophysiologically relevant information; thus, multilayer connectomics could provide a methodological framework for studies on how information flows across this stratification.