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 noninvasively accessible in humans using ultrahigh-resolution functional MRI (fMRI). Here, we applied multilayer graph theory to examine functional connectivity across different cortical depths in humans, using 7-T fMRI (1-mm3 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 interregional connections were limited to a single cortical depth only (“layer-by-layer matrices”) with 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 with 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.

With the advent of ultrahigh-resolution functional MRI (fMRI), increased noninvasive accessibility of the human cortical laminae has enabled more advanced study of the human cortex. One potential methodological approach to studying cortical laminae is through multilayer connectomics, whereby the human cortex is represented by several interconnected networks. We applied a multilayer graph theory to examine functional connectivity across different cortical depths in humans, using 7-T fMRI. We utilized global and local graph theory features that quantitatively characterize network attributes including network composition, nodal centrality, path-based measures, and hub segregation. We detected global and local functional differences between cortical depths that were augmented using multilayer connectomics compared with the layer-by-layer versions. Thus, multilayer connectomics could provide a methodological framework for studies on how the human cortex functions.

Investigating brain activity and function through network analyses has become an integral methodological foundation of neuroscience. Connectomics has yielded significant advances in understanding brain structure and function (Farahani et al., 2019; Milano et al., 2019). Modeling the brain as a system of nodes (brain regions) connected by edges (mathematical relationships)—often using graph theory—can be used to provide insight into brain characteristics and topological properties (Rubinov & Sporns, 2010). Brain networks can be derived from structural neuroimaging such as MRI or diffusion tensor imaging (DTI; structural connectomics; Griffa et al., 2013; Meoded et al., 2020; Yeh et al., 2021) or functional neuroimaging such as fMRI, EEG, or magnetoencephalography (MEG; functional connectomics; De Vico Fallani et al., 2014; Matthews & Hampshire, 2016; Sadaghiani et al., 2022; Smith et al., 2013; Xia & He, 2017). Both structural and functional connectomics have been used to understand disease models (Benito-León et al., 2019; Fleischer et al., 2019; Jacob et al., 2020; Kotlarz et al., 2022), aid in surgical mapping (Ahsan et al., 2020; Dadario et al., 2021; Gleichgerrcht et al., 2020; Hart et al., 2016), and characterize therapeutic effects of neuropsychiatric treatments (Caeyenberghs et al., 2017; Chen et al., 2020; Filippi et al., 2023; Lei et al., 2021; Tavakol et al., 2019; Yun & Kim, 2021).

An emerging field within connectomics, and more broadly graph theory, is the exploration of multilayer networks (Boccaletti et al., 2014; Kivelä et al., 2014). Multilayer networks are composed of individual layers of networks with interconnecting edges between different layers. Connections across layers can be solely between homologous nodes (multiplex) or between nodes regardless of layer or nodal position (multilayer). Multilayer connectomics enables the study of multifaceted and multimodal neuroimaging data, with the different groups of data divided into distinct layers of the connectivity matrix (Betzel & Bassett, 2017; De Domenico, 2017; Vaiana & Muldoon, 2020). For example, multilayer networks can be derived using correlations between different frequency bands of MEG recordings to identify the interplay between frequencies (Buldu & Porter, 2018). Additionally, different modalities such as MEG, fMRI, and diffusion MRI can be combined to identify patterns in brain processing (Breedt et al., 2023) or pathological dysfunction (Casas-Roma et al., 2022) that were not found in a traditional single-layer analysis. Thus, multilayer connectomics allows for the incorporation of multidimensional neuroimaging data and can identify relationships between distinct neuroimaging techniques and analyses.

One potential application of multilayer connectomics is understanding the hierarchical organization of the cerebral cortex. Neuroanatomical (Felleman & Van Essen, 1991; Rockland & Pandya, 1979; Zeki, 2018) and electrophysiological (Schroeder & Foxe, 2002; Schroeder et al., 2001) studies in animal models have identified that laminar input/output patterns can inform about bottom-up (feedforward) or top-down (feedback) processes between cortical regions. Despite its high resolution, a limitation of laminar electrophysiological recordings in comprehensive connectivity analyses is that the coverage area is typically very small. At the same time, mapping anatomical connections using fiber tracing has limited options for mapping of the postsynaptic targets (Rockland, 2019). Intracortical analyses of high-resolution fMRI data have the benefit that the coverage can be extended to the entire cerebral cortex. Advancements in ultra-high-resolution fMRI (≤1-mm3 voxel size) have enabled the sampling of functional signals from different depths of the cortical gray matter (Finn et al., 2019; Huber, Finn, et al., 2021; Norris & Polimeni, 2019; Polimeni et al., 2018). However, there are multiple challenges and unanswered questions for the feasibility of using cortical depth profiles of fMRI signals (Norris & Polimeni, 2019). Because deoxygenated blood also drains up to the cortical surface through the intracortical diving venules, fMRI voxels intersecting the superficial layers could also be affected by deeper neuronal activations (Markuerkiaga et al., 2016). Despite this limitation, studies using this emerging methodology have attempted to identify feedback and feedforward relationships noninvasively in the human brain (Chai et al., 2021; De Martino et al., 2015; Fracasso et al., 2018; Gau et al., 2020; Klein et al., 2018; Kok et al., 2016; Lankinen et al., 2023; Lawrence et al., 2019; Moerel et al., 2018, 2019; Muckli et al., 2015; Wu et al., 2018), akin to microscale recordings in animal models.

Even with the advent of high-resolution functional neuroimaging, multilayer connectomics have mostly focused on anatomical networks derived from structural MRI and DTI (Shamir & Assaf, 2021, 2023) due to their direct relationship to cortical architecture. For example, DTI and histological samples identified that cortical areas with similar laminar structure were more likely to be connected (Wei et al., 2018). Additionally, even in functional laminar studies, previous works have predominantly utilized task-based studies in predefined brain regions (Chai et al., 2021; De Martino et al., 2015; Finn et al., 2019; Fracasso et al., 2018; Gau et al., 2020; Klein et al., 2018; Kok et al., 2016; Lankinen et al., 2023; Lawrence et al., 2019; Moerel et al., 2018, 2019; Muckli et al., 2015; Polimeni et al., 2010; Wu et al., 2018). Consequently, the functional components of the whole-brain cortical depths continue to be underexplored.

In contrast to task-based studies, which primarily focus on specific cortical areas (Finn et al., 2021), a resting-state analysis enables whole-brain investigation of laminar organization (Huber, Finn, et al., 2021; Huber, Poser, et al., 2021). These resting-state connections have been shown to reflect anatomical connectivity (Adachi et al., 2012; Honey et al., 2009; Turk et al., 2016; van den Heuvel et al., 2016) and task-based networks (Di et al., 2013; Hermundstad et al., 2013). Thus, network differences within laminar resting-state fMRI networks represent functional differences between cortical depths. This work explores the laminar structure of the cortex using high-resolution resting-state fMRI and multilayer connectomics. We use a dual-pipeline approach in comparing the information extracted from layer-by-layer versus multilayer connectomics to test whether there are connectivity differences between cortical depths. We demonstrate the validity of multilayer functional laminar connectomics through showing that cortical depths have distinct graph theory characteristics that are more clearly identifiable through multilayer connectomics compared with the traditional single-layer methodology.

Participants

Thirty healthy adults (mean age ± standard deviation = 32.4 ± 10 years, 15 women, all right-handed) were recruited using an internal online recruiting platform. Participants were screened for vision problems, hearing problems, cognition-altering medications, and exclusions for MRI (metal in the body). Twenty-eight of the participants were native English speakers. Informed consent was obtained from all participants, and MRI safety screening forms were completed before each scan. The study design, protocol, and consent were approved by the Mass General Brigham Institutional Review Board.

Image Acquisition

Participants were measured in sets of 7.9-min resting-state fMRI scans occurring on different days (three to four sessions per participant). Twenty-three participants were measured in 12 resting-state scans. Seven participants had between 10 and 18 resting-state scans (10 scans: n = 1; 11 scans: n = 1; 13 scans: n = 2; 14 scans: n = 2, 18 scans: n = 1; Supporting Information Figure S1). The participants were instructed to avoid movement during the scans and to keep their eyes open and fixated on a fixation cross projected on a screen viewed through a mirror. The average duration of the sessions was around 2 hr. Breathing and heart rate were recorded using the built-in Siemens system at a sampling rate of 400 Hz. Inhalation and exhalation were measured with the Siemens respiratory-effort transducer attached to a respiratory belt. The heart rate was recorded using Siemens photoplethysmogram transducers on the participant's index finger.

The functional and structural neuroimaging data were acquired using a 7-T whole-body MRI scanner (MAGNETOM Terra, Siemens, Erlangen, Germany) with a home-built, custom-built, 64-channel array coil (Mareyam et al., 2020). To reduce participant head motion inside the scanner, MRI-compatible paddings were placed around the head and neck. In each imaging session, T1-weighted anatomical images were measured using a 0.75-mm isotropic multiecho magnetization-prepared rapid gradient echo (MPRAGE) pulse sequence (van der Kouwe et al., 2008; Zaretskaya et al., 2018) with TR = 2,530 ms; four echoes with TE of 1.72, 3.53, 5.34, and 7.15 ms; 7° flip angle, 240 × 240 mm2 FOV; and 224 sagittal slices. To help with pial surface placement by avoiding dura mater, T2-weighted anatomical images (voxel size = 0.83 × 0.83 × 0.80 mm, TR = 9,000 ms, TE = 269 ms, flip angle = 120°, FOV = 225 × 225 mm2, 270 sagittal slices) were acquired for 28 out of 30 participants in one of the imaging sessions. Resting-state functional imaging was collected using a T2*-weighted blipped-controlled aliasing in parallel imaging (CAIPI; Setsompop et al., 2012) simultaneous multislice (SMS) EPI sequence using multiband radiofrequency (RF) pulses (Setsompop et al., 2012) with 4× acceleration factor in phase-encoding direction, 3× acceleration factor in slice-encoding direction, TR = 2,800 ms, TE = 27.0 ms, isotropic 1-mm3 voxels, 78° flip angle, 192 × 192 mm2 FOV, 132 axial slices, anterior-to-posterior phase encoding direction, 1,446-Hz/pixel bandwidth, 0.82-ms nominal echo spacing, and fat suppression. In addition, to dewarp the functional data, an EPI scan was collected with identical parameters but with an opposite phase-encoding polarity (posterior-to-anterior [PA]-EPI) relative to the functional scans. For four participants with missing PA-EPI scans, the data were dewarped using a gradient-echo field map (TR = 1,040 ms, TE = 4.71 ms and 5.73 ms, isotropic 1.3-mm3 voxels; flip angle = 75°, FOV = 240 × 240 mm2, 120 slices, bandwidth = 303 Hz/pixel).

MRI Preprocessing

First, SPM12 (https://www.fil.ion.ucl.ac.uk/spm/ [SPM12-spm_preproc_run.m]; bias field correction, full-width at half-maximum [FWHM]: 18 mm, sampling distance: 2 mm, bias regularization: 1E−4) and customized MATLAB scripts were used to correct the bias field of the structural T1 and T2 images. Next, recon-all of FreeSurfer 6.0 (Fischl, 2012) with an extension for submillimeter 7-T data (Zaretskaya et al., 2018) was used to automatically create cortical reconstructions for each participant. An average of multiple T1-weighted anatomical volumes (three to four per participant) alongside a T2-weighted volume were used in the reconstruction to enhance the quality of the cortical surfaces. Nine intermediate surfaces were created between the white matter and pial surfaces with fixed relative distances, of which five were selected for the laminar analysis (described below). Lastly, the surfaces generated by recon-all were corrected manually for inaccuracies with Recon Edit of Freeview.

For the functional data, slice-timing and motion corrections were first implemented in FreeSurfer 7.1 (Fischl, 2012). Dewarping was then used to correct for geometric distortions caused by susceptibility-induced off-resonance fields. In dewarping, the off-resonance distortion field was estimated using the functional data and the PA-EPI scan collected with a reversed phase-encode blip; thus, the distortions are reversed in direction in respect to the scans (FreeSurfer: topup, applytopup; Andersson et al., 2003; Smith et al., 2004). For four participants that were missing the PA-EPI scan used above, the distortion field was estimated using the B0 field map scan in FreeSurfer 6.0 (FreeSurfer-epidewarp). The respiratory and heart rate artifacts were corrected using the RETROspective Image CORrection (RETROICOR) algorithm (third-order heart rate, respiratory, and multiplicative terms; Glover et al., 2000). Three participants were missing heart rate data, and, therefore, only respiratory recordings were used in RETROICOR. In addition, RETROICOR was not applied to five participants with missing respiratory and heart rate data. Functional data were then coregistered with the structural images using Boundary-Based Registration in FreeSurfer 6.0 (Greve & Fischl, 2009). By projecting each intersecting voxel onto the corresponding surface vertex using trilinear interpolation, the fMRI time series were then resampled onto the pial and white matter surfaces and the nine cortical depths between them.

From the nine intracortical surfaces, five alternating depths were selected starting closest to pial surface (depths 1–5, superficial to deep; Figure 1A). The outside surfaces (pial and white matter) were excluded to avoid partial volume effects from the cerebrospinal fluid and white matter, respectively. Additionally, depths included were alternated to minimize the potential for partial volume overlap between surfaces that would bias the correlation matrix generation. To measure this overlap, the distance from each voxel centroid (within cortical volume) to the cortical surface (white matter/pial surface) was calculated. The relative distance was defined so that the depth at pial surface was zero and one at the white matter border. Next, voxels intersecting each layer were picked and plotted with respect to their relative distances in a histogram. Supporting Information Figure S2 illustrates that taking every other layer limits the overlap between the layers, and thus leaking of information to adjacent layers. Additionally, to explore if temporal signal-to-noise ratio (tSNR) impacted connectivity matrix generation, the average tSNR per each cortical depth was calculated for each layer (Supporting Information Figure S3).

Figure 1.

Multilayer connectomic pipeline to analyze functional connectivity across different cortical depths. Here, “layers” refer to the dimensions of the connectivity matrices that represent fMRI signals gathered from different “cortical depths.” (A) The cortex was uniformly divided into five surfaces at different depths, as seen above projected on a 0.75-mm isotropic-resolution anatomical T1-weighted image. (B) The brain was parcellated into 148 regions of interest (ROIs) (74 per hemisphere) based on the Destrieux atlas in FreeSurfer (Destrieux et al., 2010; Fischl et al., 2004). The ROIs are shown on an inflated left-hemisphere cortical surface. (C) Schematic showing the difference between a layer-by-layer network and a multilayer network. In the layer-by-layer approach, each layer (network) is independent of other layers, while in the multilayer approach, the layers are interconnected. A sparser multilayer network is shown for visualization purposes. (D) Example matrix construction from both the layer-by-layer and multilayer approaches. While both approaches use matrices derived from Pearson correlations from the different layers, the multilayer approach generates a supra-adjacency matrix that also has correlations between different layers (shown in grayscale). (E) Example matrix construction for within-layer and between-layer matrices. For within-layer matrices, each submatrix is extracted individually for analysis. White areas represent connections excluded from the analysis.

Figure 1.

Multilayer connectomic pipeline to analyze functional connectivity across different cortical depths. Here, “layers” refer to the dimensions of the connectivity matrices that represent fMRI signals gathered from different “cortical depths.” (A) The cortex was uniformly divided into five surfaces at different depths, as seen above projected on a 0.75-mm isotropic-resolution anatomical T1-weighted image. (B) The brain was parcellated into 148 regions of interest (ROIs) (74 per hemisphere) based on the Destrieux atlas in FreeSurfer (Destrieux et al., 2010; Fischl et al., 2004). The ROIs are shown on an inflated left-hemisphere cortical surface. (C) Schematic showing the difference between a layer-by-layer network and a multilayer network. In the layer-by-layer approach, each layer (network) is independent of other layers, while in the multilayer approach, the layers are interconnected. A sparser multilayer network is shown for visualization purposes. (D) Example matrix construction from both the layer-by-layer and multilayer approaches. While both approaches use matrices derived from Pearson correlations from the different layers, the multilayer approach generates a supra-adjacency matrix that also has correlations between different layers (shown in grayscale). (E) Example matrix construction for within-layer and between-layer matrices. For within-layer matrices, each submatrix is extracted individually for analysis. White areas represent connections excluded from the analysis.

Close modal

Matrix Generation and Processing

Two parallel matrix processing pipelines were used to generate individual independent adjacency matrices for the layer-by-layer approach while creating one supra-adjacency matrix for the multilayer approach (Figures 1C and 1D). The layer-by-layer approach creates an independent network for each cortical depth while the multilayer approach results in five interconnected networks that combines all cortical depths.

One important distinction is between the terminology “depth” and “layer.” Here, depth refers to the anatomical depth in the cortex while layer refers to a specific network derived from a cortical depth. This distinction is critical to avoid equating a network layer with an associated cortical laminae.

The brain was parcellated into 148 regions of interest (ROIs; 74 per hemisphere) based on the Destrieux atlas in FreeSurfer (Destrieux et al., 2010; Fischl et al., 2004; Figure 1B). A detailed list of parcellations can be found in the Supporting Information Table S1.

Resting-state time series for each participant were concatenated across runs, leading to the following number of time points: 2,028 for n = 23 participants; 1,690 time points: n = 1; 1,859 time points: n = 1; 2,197 time points: n = 2; 2,366 time points: n = 2, and 3,042 time points: n = 1. Concatenated time series were detrended and filtered using a second-order Butterworth filter (high-pass: 0.01 Hz, low-pass: 0.1 Hz, MATLAB-filtfilt).

For the layer-by-layer approach, Pearson correlations were derived between ROIs within the same depth, resulting in 10,878 pairwise correlations from 148 nodes (ROIs) after removing 148 self-correlations (number of correlations = (nodes2 – diagonal nodes) / 2). Pearson correlation coefficients were normalized using Fisher's z-transformation, resulting in five 148 × 148 symmetric weighted connectivity matrices for each participant, that is, one matrix for each cortical depth (Figure 1D). Pearson correlation was used as opposed to partial correlations since partial correlations have been shown to perform poorly in networks with a large number of ROIs (Smith et al., 2011), the larger impact of noise and time series length on partial correlations (Liégeois et al., 2020; Matkovič et al., 2023), and the popularity of Pearson correlations in previous studies (Casas-Roma et al., 2022; Wang et al., 2014).

For the multilayer approach, Pearson correlations were derived between ROIs between and within all depths, resulting in 273,430 pairwise correlations from 740 nodes (140 ROIs × 5 depths) after removing 740 self-correlations. Pearson correlation coefficients were then normalized using Fisher's z-transformation with the final product being a 740 × 740 symmetric weighted connectivity matrix for each participant (Figure 1D).

For both approaches, individual matrices were normalized and thresholded at 2% intervals ranging from 2% to 40% graph density (ratio of edges present to total number of possible edges) to understand measure differences over a wide range of thresholds. Thresholding is required to minimize the effect spurious correlations and consider only positive correlations.

We also examined the account of within- and between-layer connections only in the context of the complete multilayer matrix. To this end, we draw two additional types of submatrices from the multilayer matrix, selectively concentrating on either their within-layer aspect (here, termed multilayer within-layer) or the between-layer aspects (termed multilayer between-layer) only. The multilayer within-layer matrices were derived by normalizing the supra-adjacency matrix, thresholding the matrix, and then extracting the nodes included in each individual layer (i.e., nodes 1 to 148 for layer 1), creating a 148 × 148 weighted connectivity matrix. The multilayer between-layer matrices were, in turn, derived by normalizing the supra-adjacency matrix, thresholding the matrix, and then zeroing the five diagonal matrices (from each cortical depth) composing within-layer connections, thus resulting in only between-layer connections. (It is worth noting that since the within-layer and between-layer connectivity matrices were extracted after thresholding, the analysis that requires normalization, that is, nonthresholded matrices, could not be conducted in the context of this analysis; Figure 1E.)

Edge Consistency and Variability

Connections (edges) within and between layers were explored to understand edge consistency and variability between participants. Edge consistency (Finn et al., 2015) was calculated by selecting the top 5% of edges with the lowest standard deviation in unthresholded multilayer networks. In contrast, edge variability (Menon & Krishnamurthy, 2019) was calculated by selecting the top 5% of edges with the highest standard deviation across participants. In both cases, edges in each layer were then summed and divided by the total number of significant edges (edges in the top 5%) to identify the percentage of significant edges in each layer.

Matrix Similarity

Matrix similarity was used to understand how matrices differed across layers. Thresholded (2%–40%) and normalized matrices were compared using cosine similarity,
where Xi and Xj are vectors of the upper triangular elements of two adjacency matrices (MATLAB-pdist2), with values ranging from −1 (maximal dissimilarity) to +1 (maximal similarity). Cosine similarity was shown to distinguish between matrices better than traditionally used Pearson correlation (Cabral et al., 2017; Menon & Krishnamurthy, 2019). Using both the layer-by-layer approach and within-layer matrices, each layer was compared with the other layers. Additionally, to understand how matrix generation differs between methods, the same layer was compared across layer-by-layer and within-layer approaches.

Connectomic Analysis

Global and nodal measures were calculated in MATLAB using the Brain Connectivity Toolbox (MATLAB Version R2022b; Rubinov & Sporns, 2010) on the Massachusetts Life Sciences Center Compute Cluster (DELL R440 servers with two Intel Xeon Silver 4214R 12-core CPUs). Global measures characterize the entire network while nodal measures characterize attributes of a specific node (ROI). Nodal measures can also be averaged to create a global measure. Measures can be grouped into four general categories to describe their overall network characterization: composition, centrality, integration, and segregation. Composition measures describe the topology of the network, while centrality measures detail-specific nodal importance for network function. Integration measures examine how information flows through the network, and segregation measures explore how the network is divided into functional components. Therefore, different measures can be used to understand different characteristics of the network. For example, decreased nodal and global average strength was found in maltreated children, indicating decreased overall brain connectivity (Puetz et al., 2017), while decreased clustering coefficient and global efficiency in patients with Parkinson's disease can signify deficits in brain network integration and segregation (Schill et al., 2023). Table 1 denotes the measures used in this work, and detailed explanation of each measure can be found in Rubinov and Sporns (2010). Additionally, small-worldness, a global quantifier that examines how “random” a network is organized, was also calculated on layer-by-layer networks since noisier data will appear more “random” (Humphries & Gurney, 2008). Small-worldness was not calculated on within-layer and between-layer connectivity matrices since the calculation requires normalization.

Table 1.

List of network measures used in this work organized by functional category

CompositionCentralityIntegrationSegregation
- Largest cluster size - Betweenness centrality (n- Characteristic path length - Maximized modularity 
- Graph density - Eigenvector centrality (n- Radius - Transitivity 
- Degree (n- Participation coefficient (n- Diameter - Clustering coefficient (n
- Strength (n - Global efficiency - Local efficiency (n
  - Assortativity  
CompositionCentralityIntegrationSegregation
- Largest cluster size - Betweenness centrality (n- Characteristic path length - Maximized modularity 
- Graph density - Eigenvector centrality (n- Radius - Transitivity 
- Degree (n- Participation coefficient (n- Diameter - Clustering coefficient (n
- Strength (n - Global efficiency - Local efficiency (n
  - Assortativity  

n: Denotes a nodal measure.

Statistical Analysis

To avoid the bias of selecting a single threshold, an area-under-the-curve (AUC) analysis was conducted to create a threshold-independent measure. The measures in Table 1 were calculated at each threshold (from 2% to 40% graph density in 2% intervals). The measure values at each threshold were then plotted against their threshold, and the area underneath the generated curve was calculated using a trapezoidal integration method (MATLAB-trapz).

For each global measure (and averaged nodal measure), the AUC value for each layer for all participants was compared using a one-way analysis of variance (ANOVA; MATLAB-anova1). A one-way ANOVA was also used to compare each nodal measure to find differences at each specific ROI (node). For the multilayer network and between-layer measures, only nodal values (both averaged and individual) were compared since global measures for the multilayer network (and thus between-layer measures) contain effects from every layer. For global values, a false discovery rate (FDR; Benjamini-Hochberg) correction (alpha = 0.05) was applied to account for multiple comparisons (Benjamini & Hochberg, 1995; Groppe, 2024). Additionally, for nodal values (nonaveraged), an FDR (Bonferroni-Holm method) correction (alpha = 0.01) was applied to account for multiple comparisons (Groppe, 2023; Holm, 1979).

Cortical Thickness Validation

One potential confounding factor using a whole-brain laminar analysis is that different brain regions have different cortical thicknesses (Barbas, 2015; Ding et al., 2009; Zachlod et al., 2020). Thus, comparing cortical thickness values of significant brain regions (defined above) can help evaluate whether our findings may be influenced by cortical thickness. Subsequently, cortical thickness values for each ROI for each participant were extracted using FreeSurfer and averaged across all subjects (Fischl, 2012). The distribution of significant nodes versus nonsignificant nodes for each nodal measure and pipeline with greater than 10 significant nodes were compared using a t test (MATLAB-ttest2).

Matrix Similarity and Edge Comparison

Edges from the multilayer matrix were compared to understand differences and similarities between participants and to see if the laminar connectomic methodology can distinguish different participants. High edge consistency indicates a similar connectivity pattern between participants, while a high edge variability increases the ability to distinguish been participants. Figure 2A and Table 2 show the percentage of consistent edge strengths between participants that each layer contains from the multilayer matrix. Layer 1 (derived from the depth closest to pial surface) has the largest number of consistent edges (29.52%), with connections between layer 1 and layer 5 (7.06%) being the most consistent between participants. In contrast, Figure 2B shows the edge variability between participants with the highest variability found in layer 5 (closest to white matter) overall (33.65%) and in within-layer connections (11.86%). It is important to clarify that edge consistency and variability are not mutually related, even though they provide complementary results above.

Figure 2.

(A) Edge consistency between each participant (multilayer matrix). The edges of each layer of the multilayer matrix were compared to find the edge strengths that had the lowest 5% standard deviation between participants. Higher values indicate a higher percentage of consistent edges, indicating consistent features between participants for those connections. (B) Edge variability between each participant (multilayer matrix). The edges of each layer of the multilayer matrix were compared to find the edge strengths that had the highest 5% standard deviation between participants. Higher values indicate a higher percentage of variable edges, indicating variable features between participants for those connections. (C) Cosine similarity between layer-by-layer and within-layer matrices (multilayer approach). Within participant, the matrix generation methods were compared using cosine similarity across a range of thresholds at each layer. Cosine similarity values range from −1 (maximal dissimilarity) to +1 (maximal similarity). The mean value at each threshold is plotted while the shaded region indicates the standard error. (D) Area-under-the-curve (AUC) measure in comparing layer-by-layer and within-layer matrix generation methods. Linear interpolation was used for visualization. The AUC from (C) is calculated using trapezoidal approximation. Higher values indicate higher similarity between methods while lower values indicate lower similarity between methods. The mean AUC value at each layer is plotted while the shaded region indicates the standard error.

Figure 2.

(A) Edge consistency between each participant (multilayer matrix). The edges of each layer of the multilayer matrix were compared to find the edge strengths that had the lowest 5% standard deviation between participants. Higher values indicate a higher percentage of consistent edges, indicating consistent features between participants for those connections. (B) Edge variability between each participant (multilayer matrix). The edges of each layer of the multilayer matrix were compared to find the edge strengths that had the highest 5% standard deviation between participants. Higher values indicate a higher percentage of variable edges, indicating variable features between participants for those connections. (C) Cosine similarity between layer-by-layer and within-layer matrices (multilayer approach). Within participant, the matrix generation methods were compared using cosine similarity across a range of thresholds at each layer. Cosine similarity values range from −1 (maximal dissimilarity) to +1 (maximal similarity). The mean value at each threshold is plotted while the shaded region indicates the standard error. (D) Area-under-the-curve (AUC) measure in comparing layer-by-layer and within-layer matrix generation methods. Linear interpolation was used for visualization. The AUC from (C) is calculated using trapezoidal approximation. Higher values indicate higher similarity between methods while lower values indicate lower similarity between methods. The mean AUC value at each layer is plotted while the shaded region indicates the standard error.

Close modal
Table 2.

Edge consistency and variability percentages for each layer derived from the multilayer matrix

LayerEdge consistency (%)Edge variability (%)
1 (Superficial) 29.52 8.25 
20.73 12.36 
15.82 19.11 
16.03 26.60 
5 (Deep) 17.89 33.65 
LayerEdge consistency (%)Edge variability (%)
1 (Superficial) 29.52 8.25 
20.73 12.36 
15.82 19.11 
16.03 26.60 
5 (Deep) 17.89 33.65 

Bolded values show the highest percentage for each measure.

After the edge analysis, connectivity matrices were compared within participants to understand the two competing matrix generation methodologies: layer-by-layer versus within-layer (multilayer) approach. Cosine similarity was used to examine similarities between different connectivity matrices. In comparing within-participant matrices across layers (Supporting Information Figure S4), layers were found to be similar, with the most distant layers (layer 1 to layer 5) having the lowest similarity in both layer-by-layer matrices and within-layer matrices. Matrices within the same layer and within participant were also compared across the matrix generation methods (Figures 2C and 2D). Layer 3 was the most consistent across the two methodologies while the peripheral layers (layer 1 and layer 5) differed the most between methods.

Single-Layer Results

Global.

Global network measures were calculated for layer-by-layer matrices. AUC values for each global measure can be found summarized in Supporting Information Table S2. Network diameter (p = 0.049) significantly increased from the most superficial layer 1 to the deepest layer 5 (Figure 3A). Largest cluster size (p = 0.0024), average betweenness centrality (p = 0.028), average local efficiency (p = 0.046), and eigenvector centrality (p = 0.028) significantly decreased from layer 1 to layer 5 (Figure 3B–E), with a peak in layer 2 for average betweenness centrality and average local efficiency. There were no significant differences in modularity, transitivity, characteristic path length, global efficiency, radius, assortativity, average degree centrality, average strength, average clustering coefficient, and average participation coefficient (Supporting Information Figures S5S6). Graph density and average degree centrality were constant across layers due to both measures being a direct function of thresholding (Supporting Information Figures S5S6). Additionally, small-worldness showed a general trend of decreasing with depth; however, there was no significant differences between layers (p = 0.2651), and small-worldness was greater than one (indicating a small-world network) for all thresholds except 40% graph density (Supporting Information Figure S7).

Figure 3.

AUC values across different layers for significant global measures (p ≤ 0.05) for the layer-by-layer analysis. Significance was calculated using a one-way ANOVA with an FDR correction (alpha = 0.05). Linear interpolation was used for visualization. The mean value across participants at each layer is plotted while the shaded region indicates the standard error. p-values shown are FDR corrected (Benjamini-Hochberg method, alpha = 0.05).

Figure 3.

AUC values across different layers for significant global measures (p ≤ 0.05) for the layer-by-layer analysis. Significance was calculated using a one-way ANOVA with an FDR correction (alpha = 0.05). Linear interpolation was used for visualization. The mean value across participants at each layer is plotted while the shaded region indicates the standard error. p-values shown are FDR corrected (Benjamini-Hochberg method, alpha = 0.05).

Close modal

Local.

Table 3 shows the number of nodes in brain regions in layer-by-layer matrices with significant differences between layers (FDR correction with alpha = 0.01; See Supporting Information Table S3 for specific values and regions). Degree centrality, strength, and eigenvector centrality had the greatest number of significant nodes (4/148; Table 3). The limbic region had more significant nodes than all other regions for each measure, except for clustering coefficient, which was tied with the temporal region (one significant node for each region). In all measures, the right hemisphere had more significant nodes than the left hemisphere (Table 3). In general, the most superficial layers (1 and 2) had the highest value for significant nodes (Table 3, Supporting Information Table S3, Supporting Information Figures S8S9). Significant nodes were distributed across node thickness levels (Supporting Information Figures S10S11).

Table 3.
Number of significant nodes within each brain region for a layer-by-layer analysis
Number of significant nodes within each brain region for a layer-by-layer analysis

Significance was calculated from the AUC values using a one-way ANOVA with an FDR correction (alpha = 0.01) to account for multiple comparisons (Groppe, 2023; Holm, 1979). Details of nodal mapping to each region can be found in Supporting Information Table S1. : hemisphere with the highest number of nodes; : measure with the highest number of nodes across both hemispheres; : region within each measure with the highest number of nodes.

Multilayer Results

Within-layer global.

Figure 4 shows global network measures calculated for within-layer matrices. AUC values for each global measure can be found summarized in Supporting Information Table S4. Characteristic path length (p = 0.014) and diameter (p < 0.001) all increased from layer 1 to layer 5 (Figure 4A–B). Largest cluster size (p < 0.001), graph density (p < 0.001), average degree centrality (p < 0.001), average strength (p < 0.001), average eigenvector centrality (p < 0.001), and average participation coefficient (p = 0.0011) significantly decreased with cortical depth (layer 1 to 5; Figure 4C–H). There were no significant differences for modularity, transitivity, global efficiency, radius, assortativity, average betweenness centrality, average clustering coefficient, and average local efficiency (Supporting Information Figures S12S13). In contrast to layer-by-layer results, graph density and average degree centrality were different across layers due to the within-layer matrix generation methodology, allowing each individual layer to have a different graph density.

Figure 4.

AUC values across different layers for significant global measures (p ≤ 0.05) for a within-layer analysis. Significance was calculated using a one-way ANOVA with an FDR correction (alpha = 0.05). Linear interpolation was used for visualization. The mean value across participants at each layer is plotted while the shaded region indicates the standard error. p-values shown are FDR corrected (Benjamini-Hochberg method, alpha = 0.05).

Figure 4.

AUC values across different layers for significant global measures (p ≤ 0.05) for a within-layer analysis. Significance was calculated using a one-way ANOVA with an FDR correction (alpha = 0.05). Linear interpolation was used for visualization. The mean value across participants at each layer is plotted while the shaded region indicates the standard error. p-values shown are FDR corrected (Benjamini-Hochberg method, alpha = 0.05).

Close modal

Within-layer local.

The number of nodes in brain regions with significant differences using within-layer matrices can be found in Table 4 (FDR correction with alpha = 0.01; see Supporting Information Table S5 for specific values and regions). Degree centrality had the largest number of significant nodes (15/148) followed by strength (9/148), participation coefficient (9/148), and local efficiency (8/148; Table 4). In all measures except participation coefficient, the limbic region had the most significant nodes; in participation coefficient, the temporal region had the most significant nodes (4/32). In all measures, the right hemisphere had more significant nodes than the left hemisphere (Table 4). For significant nodes, either layer 1 or 2 had the highest value (Table 4, Supporting Information Table S5, Supporting Information Figures S14S15). Significant nodes were spread across different node thicknesses (Supporting Information Figures S16S17).

Table 4.
Number of significant nodes within each brain region for a within-layer analysis
Number of significant nodes within each brain region for a within-layer analysis

Significance was calculated from the AUC values using a one-way ANOVA with an FDR correction (alpha = 0.01) to account for multiple comparisons (Groppe, 2023; Holm, 1979). Details of nodal mapping to each region can be found in Supporting Information Table S1. : hemisphere with the highest number of nodes; : measure with the highest number of nodes across both hemispheres; : region within each measure with the highest number of nodes.

Multilayer global.

Alongside measures for individual layers, global measures were calculated for the supra-adjacency matrix created using the multilayer approach. While only global values derived from nodal averages were statistically compared between layers, AUC values for all global measures are summarized in Supporting Information Table S6. Average degree centrality (p < 0.001), average strength (p = 0.0019), and average eigenvector centrality (p < 0.001) decreased from layer 1 to layer 5 with a slight peak at layer 2 (Figure 5A–C). Average betweenness centrality (p < 0.001) and average participation coefficient (p = 0.019) also decreased from layer 1 to layer 5 (Figure 5D–E). In contrast, average clustering coefficient tended to increase from layer 1 to layer 5 (p = 0.080; Supporting Information Figures S18S19, Supporting Information Table S6). Average local efficiency showed no significant difference between layers (Supporting Information Figures S18S19, Supporting Information Table S6). Additionally, layer-wise graph density, graph density derived from each individual layer within the multilayer network, decreased from layer 1 to layer 5 with a slight peak at layer 2 (Supporting Information Figure S20).

Figure 5.

AUC values across different layers for significant global measures (p ≤ 0.05) for a multilayer analysis. Significance was calculated using a one-way ANOVA with an FDR correction (alpha = 0.05). Linear interpolation was used for visualization. The mean value across participants at each layer is plotted while the shaded region indicates the standard error. p-values shown are FDR corrected (Benjamini-Hochberg method, alpha = 0.05).

Figure 5.

AUC values across different layers for significant global measures (p ≤ 0.05) for a multilayer analysis. Significance was calculated using a one-way ANOVA with an FDR correction (alpha = 0.05). Linear interpolation was used for visualization. The mean value across participants at each layer is plotted while the shaded region indicates the standard error. p-values shown are FDR corrected (Benjamini-Hochberg method, alpha = 0.05).

Close modal

Multilayer local.

Table 5 depicts the number of nodes in specific brain regions that were significantly different between layers for the multilayer analysis (FDR correction with alpha = 0.01; see Supporting Information Table S7 for specific values and regions). Betweenness centrality identified the most significant nodes (58/148) with 19/48 frontal nodes, 16/28 limbic nodes, 2/20 occipital nodes, 5/20 parietal nodes, and 16/32 temporal nodes (Table 5). The limbic region had highest percentage of nodes in all measures except clustering coefficient (4/32 in temporal) and local efficiency (no significant nodes found). Additionally, using betweenness centrality, the frontal region had the highest number of significant nodes (19/48). Again, in all measures, the right hemisphere had more nodes with significant differences between layers (Table 5). While most significant regions across measures were highest in the superficial layers (layers 1, 2), especially in betweenness centrality, the deepest layer (layer 5) had the highest values for nodes significant in clustering coefficient (Table 5, Supporting Information Table S7, Supporting Information Figure S21). The thickness of significant nodes was spread across the spectrum of thickness levels, with a preference toward thicker nodes for betweenness centrality (Supporting Information Figures S22S23).

Table 5.
Number of significant nodes within each brain region for a multilayer analysis
Number of significant nodes within each brain region for a multilayer analysis

Significance was calculated from the AUC values using a one-way ANOVA with an FDR correction (alpha = 0.01) to account for multiple comparisons (Groppe, 2023; Holm, 1979). Details of nodal mapping to each region can be found in Supporting Information Table S1. : hemisphere with the highest number of nodes; : measure with the highest number of nodes across both hemispheres; : region within each measure with the highest number of nodes.

Between-layer global.

Significant between-layer global measures can be found in Figure 6. AUC values for all global measures can be found in Supporting Information Table S8. Average degree centrality (p < 0.001), average strength (p = 0.0021), average eigenvector centrality (p < 0.001), average betweenness centrality (p < 0.001), and average participation coefficient (p = 0.032) decreased from superficial (layer 1) to deep (layer 5) with all except average participation coefficient showing a peak in layers 2 and 3 (Figures 6A–D and 6F). Average clustering coefficient (p = 0.0033) increased from layer 1 to layer 5 with a slight decrease from layer 1 to layer 2 (Figure 6E). Average local efficiency was the only measure that showed no significant difference between layers (Supporting Information Figures S24S25, Supporting Information Table S8).

Figure 6.

AUC values across different layers for significant global measures (p ≤ 0.05) for a between-layer analysis. Significance was calculated using a one-way ANOVA with an FDR correction (alpha = 0.05). Linear interpolation was used for visualization. The mean value across participants at each layer is plotted while the shaded region indicates the standard error. p-values shown are FDR corrected (Benjamini-Hochberg method, alpha = 0.05).

Figure 6.

AUC values across different layers for significant global measures (p ≤ 0.05) for a between-layer analysis. Significance was calculated using a one-way ANOVA with an FDR correction (alpha = 0.05). Linear interpolation was used for visualization. The mean value across participants at each layer is plotted while the shaded region indicates the standard error. p-values shown are FDR corrected (Benjamini-Hochberg method, alpha = 0.05).

Close modal

Between-layer local.

Brain regions with nodes that were significantly different using between-layer matrices can be found in Table 6 (FDR correction with alpha = 0.01; see Supporting Information Table S9 for specific values and regions). Betweenness centrality had the largest number of significant nodes (57/148) followed by clustering coefficient (22/148; Table 6). For betweenness centrality, most significant nodes had the highest values in the superficial layers (layers 1 and 2), while for clustering coefficient, all significant nodes were highest in layer 5 (Supporting Information Figure S26). In degree centrality, strength, eigenvector centrality, and betweenness centrality, the limbic region had the highest percentage of significant nodes within each region (Table 6). In betweenness centrality, frontal region nodes had the highest absolute number of significant nodes (20/48). In clustering coefficient (10/32) and participation coefficient (4/32), the temporal region had the most significant nodes (Table 6). Local efficiency had no significant nodes (Table 6, Supporting Information Table S9, Supporting Information Figure S26). The right hemisphere had more significant nodes than the left hemisphere for all measures (Table 6). Significant nodes were dispersed across different thickness levels, with betweenness centrality nodes leaning slightly more toward thicker regions (Supporting Information Figures S27S28).

Table 6.
Number of significant nodes within each brain region for a between-analysis
Number of significant nodes within each brain region for a between-analysis

Significance was calculated from the AUC values using a one-way ANOVA with an FDR correction (alpha = 0.01) to account for multiple comparisons (Groppe, 2023; Holm, 1979). Details of nodal mapping to each region can be found in Supporting Information Table S1. : hemisphere with the highest number of nodes; : measure with the highest number of nodes across both hemispheres; : region within each measure with the highest number of nodes.

Method Comparison

All four network measure methods (layer-by-layer, within-layer, multilayer, between-layer) identified global differences between layers (Figure 7). Layer-by-layer and within-layer methods showed an increase in an integration-based global measure (diameter) from layer 1 to layer 5; however, the within-layer approach identified an increase in characteristic path length as well. Similarly, layer-by-layer and within-layer approaches identified a decrease in largest cluster size from layer 1 to layer 5, with an additional decrease found in graph density using the within-layer method.

Figure 7.

Layer with the highest AUC value across different global network measures methods (p ≤ 0.05). Significance was calculated using a one-way ANOVA with an FDR correction (alpha = 0.05). White sections indicate no significant difference between layers. Gray sections indicate network measures that cannot separate effects of different layers.

Figure 7.

Layer with the highest AUC value across different global network measures methods (p ≤ 0.05). Significance was calculated using a one-way ANOVA with an FDR correction (alpha = 0.05). White sections indicate no significant difference between layers. Gray sections indicate network measures that cannot separate effects of different layers.

Close modal

For nodal-averaged global measures, which can be applied to all four methods, measures generally decreased from superficial layers (1 and 2) to deeper layers (Figure 7). The layer-by-layer methodology identified significant differences in three measures (average eigenvector centrality, average betweenness centrality, and average local efficiency) while the within-layer approach found significant differences in four measures (average degree centrality, average strength, average eigenvector centrality, and average participation coefficient). The multilayer approach also found differences in five measures (average degree centrality, average strength, average eigenvector centrality, average betweenness centrality, and average participation coefficient). Interestingly, the between-layer method identified the most measures with significant differences between layers using six measures, with two measures (average strength and average eigenvector centrality) peaking in the middle layer. Additionally, the between-layer approach was the only method to identify a difference in average clustering coefficient, which was the highest in layer 5.

Similar to global measures, the multilayer-based approaches (within-layer, multilayer, between-layer) identified more nodal differences between layers than the layer-by-layer approach (Tables 36). For example, the layer-by-layer method identified a maximum of four nodes with significant differences per measures compared with 15 for within-layer (degree centrality), 58 for multilayer (betweenness centrality), and 57 for between-layer (betweenness centrality). Despite this, in a majority of measures for all four methods, the limbic region had the greatest number of nodes with significant differences between layers. The right hemisphere also had more significant nodes across all methods and measures (Tables 36). The right (2.4127 ± 0.0306 mm; AVG ± SE) and left (2.4069 ± 0.0303 mm) hemispheres had comparable cortical thicknesses overall and across brain regions (Supporting Information Figure S29). Additionally, while the thickness of significant nodes versus nonsignificant nodes was significantly larger in betweenness centrality metrics (multilayer, between-layer), the absolute difference between significant and nonsignificant nodes was typically less than 1 mm (Supporting Information Figure S30).

Two measures that showed a considerable benefit from the multilayer-based approach were betweenness centrality and clustering coefficient (Figure 8). For example, the multilayer and between-layer methods showed substantial increase in the number of nodes that had significant differences between layers. Similarly, the number of nodes with significant differences between layers in clustering coefficient increased using multilayer and between-layer methods. More importantly, however, clustering coefficient in multilayer and between-layer approaches is the only measure to highlight the deepest layer as having the largest value. Likewise, the multilayer and between-layer methods are the only methods to include nodes that are the highest value in the middle layer (Supporting Information Figures S21 and S26).

Figure 8.

Nodes with significant differences between layers for each nodal measure pipeline (layer-by-layer, within-layer, multilayer, between-layer) for (A) betweenness centrality and (B) clustering coefficient. Significance was calculated from the AUC values using a one-way ANOVA with an FDR correction (alpha = 0.01) to account for multiple comparisons (Groppe, 2023; Holm, 1979). The colored section represents the layer with the highest value for the node. The nodes are based on the Destrieux atlas in FreeSurfer (Destrieux et al., 2010; Fischl et al., 2004). LH = left hemisphere; RH = right hemisphere.

Figure 8.

Nodes with significant differences between layers for each nodal measure pipeline (layer-by-layer, within-layer, multilayer, between-layer) for (A) betweenness centrality and (B) clustering coefficient. Significance was calculated from the AUC values using a one-way ANOVA with an FDR correction (alpha = 0.01) to account for multiple comparisons (Groppe, 2023; Holm, 1979). The colored section represents the layer with the highest value for the node. The nodes are based on the Destrieux atlas in FreeSurfer (Destrieux et al., 2010; Fischl et al., 2004). LH = left hemisphere; RH = right hemisphere.

Close modal

The graph theory analysis of high-resolution (7 T) resting-state fMRI revealed global and nodal network differences between cortical depths. Global integration measures (diameter, characteristic path length) were higher in deeper layers, while composition (largest cluster size, graph density, degree centrality, strength) and centrality (eigenvector centrality, betweenness centrality) measures were often higher in superficial layers (Figure 7).

Current literature exploring graph theory measures and overall laminar connectivity through networks in the human brain is very limited. Structural analysis of the human connectome using diffusion MRI combined with T1-weighted anatomical imaging found qualitative differences in degree, strength, and betweenness centrality nodal distributions across cortical depths (Shamir et al., 2022); however, network-wide global calculations were absent. In a functional approach, Deshpande et al. (2022) found no global differences between layers using mean blind deconvoluted Pearson correlations from resting-state fMRI. However, no threshold was used, enabling spurious correlations to impact the mean. Additionally, using the mean across the whole brain obfuscates any impact of a particular ROI. Our results, thus, significantly extend global and nodal network analysis of cortical architecture across the entire brain.

Our findings provide evidence of an advantage of applying multilayer graph theory to connectomic analysis. While differences between layers were seen across all methodologies, the multilayer approach provided a greater identification of these differences through identifying more measures with larger significant differences (Figure 7). Previous connectomic studies have shown a benefit of applying a multilayer network framework (Betzel & Bassett, 2017; De Domenico, 2017; Vaiana & Muldoon, 2020). For example, multilayer connectomics enables the integration of complex neuroimaging data (cortical lamina, frequency bands, multimodal neuroimaging; Betzel & Bassett, 2017; Breedt et al., 2023; Buldu & Porter, 2018; Casas-Roma et al., 2022; De Domenico, 2017; Shamir & Assaf, 2023; Vaiana & Muldoon, 2020) and the creation of new network features. New network features can be used to explain neuroscientific findings, as in this work, or even enhance machine learning workflows to better discriminate between disease states (Zhu et al., 2022). Future connectomics studies with available data should therefore consider using a multilayer framework to augment brain network modeling and analysis.

One particular benefit of multilayer analysis in laminar connectomics is the ability to discriminate between and incorporate the impact of within- and between-layer connections. While this comparison was limited to nodal averaged global measures and nodal values, there was still a stark comparison between within- and between-layer connections. When exploring layers individually (within-layer connections), the most superficial layer had the highest activation and is densely connected to itself while the deepest layer was relatively sparsely connected and took longer network paths to transmit information to different brain regions (Figures 4 and 7). However, in between layer connections, layer 1 becomes less important and the superficial middle layers (layers 2 and 3) become integral for cortical connectivity (Figures 67). The superficial middle layers (layers 2 and 3) had the highest between-layer degree, strength, and eigenvector centrality, indicating that both layers are densely connected to other layers. Additionally, layer 2's significantly larger betweenness centrality demonstrates that it is the most important layer for information flow between layers (Figures 67). Lastly, the deepest layer had the highest clustering coefficient, meaning that it had the highest likelihood of forming local hubs (triangles) with other layers. Thus, within- and between-layer analysis provides evidence of a highly active superficial layer that utilizes layers 2 and 3 to transmit information to other cortical layers.

One important note to contextualize the above findings is that while activity and hubs of information flow can be identified, the direction of information flow cannot be delineated. Cortical layer architecture can have diverse connectivity patterns across layers and hierarchal schemes (Felleman & Van Essen, 1991). Therefore, it is important to recognize the correlation nature of this work rather than infer causality. Furthermore, while our results primarily focus on layers with the highest measure/activity, this does not imply that other layers are inactive. This contextualization is notably important in the context of feedback/feedforward interpretations, a common framework for laminar analysis. In general, feedback is thought to target superficial and deep layers, and feedforward targets the middle layer (Barbas, 2015; Felleman & Van Essen, 1991; Rockland & Pandya, 1979). However, the fact that one area might also send information makes the interpretation less straightforward. Thus, the information transfer and clustering processes might not directly reflect the feedback/feedforward processes but complex interactions between them. This intricate feedback/feedforward interplay may also explain the connections between layer 1 and layer 5 (7.06%) being the most consistent between participants (Figure 2A), with both superficial and deep layers activated on opposite sides of the feedback/feedforward circuit (Miyashita, 2022). However, this feedback/feedforward interaction would need to be further validated and studied before drawing concrete conclusions.

The significance of connectivity patterns and characteristics for cortical depths differed across different brain regions. We identified extensive differences between layers in the frontal, limbic, and temporal brain regions (Figure 8, Tables 36). Interestingly, the limbic cortex, often with the most significant regions per measures, typically has less layers than other brain regions (Barbas, 2015). Thus, cellular architecture may play a role in the ability for laminar fMRI, and potentially multilayer connectomics, to detect differences between cortical layers. Cortical thickness may also play a role in detecting differences between layers (Fukutomi et al., 2018). The regions (frontal, limbic, and temporal) with the most differences were often the thickest regions (Supporting Information Figure S29), with our results overlapping with previous studies (Fukutomi et al., 2018). Additionally, other factors that may impact cortical function and detectability include the neurite density index, orientation dispersion index, and myelin (Fukutomi et al., 2018). However, Fukutomi's et al. (2018) findings show a varied distribution across regions for these measures. Despite this variation, hotspots in these metrics near the posterior-ventral part of the cingulate gyrus and transverse temporal sulcus often overlap with significant nodes in our work. Therefore, our results indicate the need to contextualize layer fMRI results within cortical metrics while providing a framework for potential regions (frontal, limbic, temporal) that may be suited for a whole-brain laminar analysis.

On a global network level, there were high levels of activation in superficial layers when compared with deeper layers, in line with previous resting-state fMRI analyses reporting increased activation patterns in superficial cortical depths (Guidi et al., 2020; Huber, Finn, et al., 2021; Huber, Poser, et al., 2021; Polimeni et al., 2010). This was also reflected in the higher graph density in superficial layers (Supporting Information Figure S20), despite marginally higher tSNR in deeper layers (Supporting Information Figure S3). Higher composition and centrality measures indicate a more robustly connected network (Figure 7). Additionally, at least within the same network (layer-by-layer and within-layer), deeper layers had significantly longer paths to transmit information, shown by higher global integration measures. It is, however, important to consider the pial vein bias (Polimeni et al., 2010), which increases gradient-echo BOLD signals from the deep to the more superficial parts of cortex. Further studies with alternative contrast mechanisms, which are less affected by the draining vein effect, are thus needed to determine whether the superficial cortical depths play a more critical role in the brain at rest, compared with the deeper aspects of the cortex or if this result is a function of signal strength and variance increasing toward the superficial depth.

Limitations

This study does have some limitations, both in terms of laminar analysis and connectomic analysis. Regarding our laminar analysis approach, the limitations of this study include biases associated with our fMRI pulse sequence, signal-to-noise ratio (SNR) bias, the number of cortical depths chosen, the depth-to-cortical layer incongruence, and the impact of resting-state versus task-based paradigms. This work uses gradient-echo BOLD pulse sequences, which may not be fully optimal for laminar analysis (Bandettini et al., 2021; Huber, Finn, et al., 2021). Gradient-echo BOLD can be influenced largely by large draining vessels (Markuerkiaga et al., 2016; Olman et al., 2007; Polimeni et al., 2010; Turner, 2002), while spin-echo BOLD (Duong et al., 2003; Uludaü et al., 2009; Yacoub et al., 2003; Zhao et al., 2004) and VAscular Space Occupancy (VASO; Chai et al., 2020, 2021) have been proposed as alternative fMRI contrasts for laminar analysis to address this large vein bias (Bandettini et al., 2021; Huber, Finn, et al., 2021). However, VASO and spin-echo BOLD have lower sensitivity and several practical challenges (Moerel et al., 2021). Similarly, SNR can vary at different cortical depths. For example, depths within the middle of the cortex will contain less tissue boundary effects compared with the depths near the pial and white matter surfaces (Blazejewska et al., 2019). This difference may be further exacerbated since the thickness and functionality of cortical layers can change based on the brain region (Barbas, 2015; Ding et al., 2009; Zachlod et al., 2020) and cortical curvature (Fatterpekar et al., 2003; Fischl & Dale, 2000; Hilgetag & Barbas, 2006; Van Essen & Maunsell, 1980). However, as shown above for nodal analysis, thickness varied across statistically significant nodes, suggesting that our results are not purely a function of cortical thickness since significant differences were identified in “thinner” nodes (Supporting Information Figures S10S11, S16S17, S22S23, S27S30). However, for multilayer (p < 0.001) and between-layer (p < 0.001) betweenness centrality (Supporting Information Figure S30), nodes with significant differences between layers had higher cortical thickness. This finding could either result from methodological constraints indicating an inability to detect differences at lower cortical thicknesses or a neurophysiological phenomenon of thicker nodes having larger functional differences between layers. Cortical curvature was not explored, and future laminar work should include the anatomical constraints of the cortex to address this. In addition to the location of the cortical depths chosen, the number of depths can affect the results. Other studies have used a smaller number of depths to ensure independence between depths (Sharoh et al., 2019), six depths to match the number of cortical layers (Pais-Roldán et al., 2023), or even a larger amount that showed an improved detection of cortical responses (Huber et al., 2017). The number of depths chosen should balance independence, cortical response detection, and computational demands from a higher depth count. The number of depths can also impact the role of partial volume effects due to voxel overlap. Furthermore, as mentioned above, the cortical depths do not directly equate to cytoarchitectural cortical layers. Lastly, this study used resting-state fMRI to study whole-brain connectivity; however, laminar resting-state fMRI activation patterns may be different than laminar task-based patterns (Pais-Roldán et al., 2023), limiting the broad applicability to task-based laminar paradigms. Despite potential activation pattern differences, the underlying anatomical basis of resting-state connections (Adachi et al., 2012; Honey et al., 2009; Turk et al., 2016; van den Heuvel et al., 2016) can still inform task-based paradigms. Ideally, a second dataset would be utilized to validate our results; however, few comparable datasets are available.

Regarding our connectomic analysis, limitations include the parcellation choice, network construction approach, thresholding methodology, and multilayer measure calculations. Parcellation choice can impact graph theory results (Albers et al., 2021; Arslan et al., 2018). This work used the Destrieux atlas in FreeSurfer (Destrieux et al., 2010; Fischl et al., 2004), which is based on anatomical nomenclature. However, an atlas derived from functional connectivity (Schaefer et al., 2018) or utilizing functional localizers for specific areas of interest (Nieto-Castañón & Fedorenko, 2012) may be more appropriate for a functional analysis study. Additionally, for laminar analysis, a custom atlas using laminar cytoarchitecture and cortical thickness may improve the accuracy of the results. Another impactful choice in connectomic methodology is how to construct the network from the fMRI time series. Pearson correlations perform better for network construction when using a large number of ROIs (Smith et al., 2011) or in noisier data (Liegeois et al., 2020; Matkovič et al., 2023). However, Pearson correlations also include indirect effects of ROIs, which can alter analysis. A popular alternative is partial correlation, which utilizes the inverse covariance matrix and, thus, excludes the indirect network effects. However, partial correlations can also include spurious connections (Berkson's paradox; Berkson, 1946) and tend to increase network construction complexity since partial correlations require regularization that has varying optimization parameters (Kim et al., 2015; Pervaiz et al., 2020). Thus, partial correlations may have future utility in laminar multilayer analysis; however, limited ROI number in relation to fMRI time series data points and lack of optimized laminar multilayer regularization parameters led to the Pearson correlation being used in this work. Also, as mentioned above, graph theory measures are directly impacted by thresholding the network (Osmanlıoğlu et al., 2020). The AUC analysis attempts to correct for this thresholding bias but still may be inadequate for eliminating the thresholding's effect on network characteristics. Additionally, network measures may be impacted as a result of SNR and network layer normalization (Mandke et al., 2018). For example, increased noise will transition network structure from small world to more random (Humphries & Gurney, 2008), which may occur as we measure deeper into the cortex. However, our results showed no significant differences in small-worldness between layers, indicating that this network structure change is not occurring in our work (Supporting Information Figure S7). Similarly, within multilayer approaches, normalization plays a key role since graph density can influence network properties. However, even with comparable tSNR (Supporting Information Figure S3), the graph density varied across cortical depths. While we believe this to be an intrinsic property of cortical connectivity having higher density in superficial cortical depths, as demonstrated by other studies (Logothetis et al., 2001), future work should explore different multilayer normalization schemes in laminar connectivity to more thoroughly parse through this effect (Mandke et al., 2018). Lastly, our statistical analysis of our measures may be limited to the use of ANOVA since it assumes normality and equal variance, which are sometimes violated by network measures. Additionally, when selecting peak values for each measure, the highest value was selected instead of using a planned contrast ANOVA.

Increased BOLD signal in superficial versus deeper layers may be due to vascular-related bias (Markuerkiaga et al., 2016; Olman et al., 2007; Pais-Roldán et al., 2023; Polimeni et al., 2010; Turner, 2002). One might conclude that the present results reflect vascular biases. The most superficial depth was excluded in this work to reduce this bias; however, the other layers will still have some effect of vascular draining. Additionally, even with the removal of the most superficial depth, the current most superficial depth may still be including superficial voxels that are sensitive to vascular-related bias. Despite this limitation, some composition and centrality measures peaked in layers 2 and 3, notably average strength, suggesting that some observed effects are not explainable by biases in superficial layers (Figure 7, Supporting Information Figures S10S11, S16S17, S22S23, S27S30). Lastly, even the utility of the multilayer approach to find more significant effects may be a result of the multilayer model being more sensitive to draining/signal confounds.

Conclusion

Our multilayer connectomics findings demonstrate global and nodal network differences between cortical depths that can be more aptly identified through the multilayer approach compared with traditional single-layer connectomics. These results demonstrate the validity of the multilayer connectomic framework on laminar fMRI and provide a methodological foundation for future multilayer laminar studies. Future work should further explore the intersection of connectomics and laminar studies and address current methodological constraints.

Supporting information for this article is available at https://doi.org/10.1162/netn_a_00436.

Parker Kotlarz: Conceptualization; Data curation; Formal analysis; Investigation; Methodology; Software; Validation; Visualization; Writing – original draft; Writing – review & editing. Kaisu Lankinen: Conceptualization; Investigation; Methodology; Software; Supervision; Visualization; Writing – review & editing. Maria Hakonen: Data curation; Investigation; Methodology; Software; Writing – review & editing. Tori Turpin: Investigation; Writing – review & editing. Jonathan R. Polimeni: Methodology; Software; Writing – review & editing. Jyrki Ahveninen: Conceptualization; Funding acquisition; Investigation; Methodology; Project administration; Supervision; Writing – review & editing.

Jyrki Ahveninen, National Institute on Deafness and Other Communication Disorders (https://dx.doi.org/10.13039/100000055), Award ID: R01DC017991. Jyrki Ahveninen, National Institute on Deafness and Other Communication Disorders (https://dx.doi.org/10.13039/100000055), Award ID: R01DC016765. Jyrki Ahveninen, National Institute on Deafness and Other Communication Disorders (https://dx.doi.org/10.13039/100000055), Award ID: R01DC016915. Not Applicable, National Institute of Biomedical Imaging and Bioengineering (https://dx.doi.org/10.13039/100000070), Award ID: P41-EB030006. Not Applicable, NIH Office of the Director (https://dx.doi.org/10.13039/100000052), Award ID: S10-OD023637. Award ID: K99DC022305, Maria Hakonen, National Institute on Deafness and Other Communication Disorders.

Multilayer connectomics:

A graph theory methodology that utilizes networks connected to other networks.

Cortical depth:

A measure of distance into the cortex.

Edge consistency:

Edges in the top 5% that are shared between subjects.

Edge variability:

Edges in the bottom 5% that are different between subjects.

Clustering coefficient:

A network measure that quantifies the extent of grouping of a node with other nodes in the network.

Betweenness centrality:

A network measure that quantifies how often a node is involved in a shortest path of the network.

Ultrahigh-resolution fMRI:

Functional neuroimaging that allows noninvasive submillimeter resolution of the brain.

Cortical laminae:

Individual cellular layers that comprise the human cortex.

Ultra-high-resolution fMRI:

Functional neuroimaging that allows noninvasive submillimeter resolution of the brain.

Cortical laminae:

Individual cellular layers that comprise the human cortex.

Adachi
,
Y.
,
Osada
,
T.
,
Sporns
,
O.
,
Watanabe
,
T.
,
Matsui
,
T.
,
Miyamoto
,
K.
, &
Miyashita
,
Y.
(
2012
).
Functional connectivity between anatomically unconnected areas is shaped by collective network-level effects in the macaque cortex
.
Cerebral Cortex
,
22
(
7
),
1586
1592
. ,
[PubMed]
Ahsan
,
S. A.
,
Chendeb
,
K.
,
Briggs
,
R. G.
,
Fletcher
,
L. R.
,
Jones
,
R. G.
,
Chakraborty
,
A. R.
, …
Sughrue
,
M. E.
(
2020
).
Beyond eloquence and onto centrality: A new paradigm in planning supratentorial neurosurgery
.
Journal of Neuro-oncology
,
146
(
2
),
229
238
. ,
[PubMed]
Albers
,
K. J.
,
Ambrosen
,
K. S.
,
Liptrot
,
M. G.
,
Dyrby
,
T. B.
,
Schmidt
,
M. N.
, &
Mørup
,
M.
(
2021
).
Using connectomics for predictive assessment of brain parcellations
.
NeuroImage
,
238
,
118170
. ,
[PubMed]
Andersson
,
J. L. R.
,
Skare
,
S.
, &
Ashburner
,
J.
(
2003
).
How to correct susceptibility distortions in spin-echo echo-planar images: Application to diffusion tensor imaging
.
NeuroImage
,
20
(
2
),
870
888
. ,
[PubMed]
Arslan
,
S.
,
Ktena
,
S. I.
,
Makropoulos
,
A.
,
Robinson
,
E. C.
,
Rueckert
,
D.
, &
Parisot
,
S.
(
2018
).
Human brain mapping: A systematic comparison of parcellation methods for the human cerebral cortex
.
NeuroImage
,
170
,
5
30
. ,
[PubMed]
Bandettini
,
P. A.
,
Huber
,
L.
, &
Finn
,
E. S.
(
2021
).
Challenges and opportunities of mesoscopic brain mapping with fMRI
.
Current Opinion in Behavioral Sciences
,
40
,
189
200
.
Barbas
,
H.
(
2015
).
General cortical and special prefrontal connections: Principles from structure to function
.
Annual Review of Neuroscience
,
38
,
269
289
. ,
[PubMed]
Benito-León
,
J.
,
Sanz-Morales
,
E.
,
Melero
,
H.
,
Louis
,
E. D.
,
Romero
,
J. P.
,
Rocon
,
E.
, &
Malpica
,
N.
(
2019
).
Graph theory analysis of resting-state functional magnetic resonance imaging in essential tremor
.
Human Brain Mapping
,
40
(
16
),
4686
4702
. ,
[PubMed]
Benjamini
,
Y.
, &
Hochberg
,
Y.
(
1995
).
Controlling the false discovery rate: A practical and powerful approach to multiple testing
.
Journal of the Royal Statistical Society: Series B (Methodological)
,
57
(
1
),
289
300
.
Berkson
,
J.
(
1946
).
Limitations of the application of fourfold table analysis to hospital data
.
Biometrics
,
2
(
3
),
47
53
. ,
[PubMed]
Betzel
,
R. F.
, &
Bassett
,
D. S.
(
2017
).
Multi-scale brain networks
.
NeuroImage
,
160
,
73
83
. ,
[PubMed]
Blazejewska
,
A. I.
,
Fischl
,
B.
,
Wald
,
L. L.
, &
Polimeni
,
J. R.
(
2019
).
Intracortical smoothing of small-voxel fMRI data can provide increased detection power without spatial resolution losses compared to conventional large-voxel fMRI data
.
NeuroImage
,
189
,
601
614
. ,
[PubMed]
Boccaletti
,
S.
,
Bianconi
,
G.
,
Criado
,
R.
,
Del Genio
,
C. I.
,
Gómez-Gardeñes
,
J.
,
Romance
,
M.
, …
Zanin
,
M.
(
2014
).
The structure and dynamics of multilayer networks
.
Physics Reports
,
544
(
1
),
1
122
. ,
[PubMed]
Breedt
,
L. C.
,
Santos
,
F. A. N.
,
Hillebrand
,
A.
,
Reneman
,
L.
,
van Rootselaar
,
A.-F.
,
Schoonheim
,
M. M.
, …
Douw
,
L.
(
2023
).
Multimodal multilayer network centrality relates to executive functioning
.
Network Neuroscience
,
7
(
1
),
299
321
. ,
[PubMed]
Buldú
,
J. M.
, &
Porter
,
M. A.
(
2018
).
Frequency-based brain networks: From a multiplex framework to a full multilayer description
.
Network Neuroscience
,
2
(
4
),
418
441
. ,
[PubMed]
Cabral
,
J.
,
Vidaurre
,
D.
,
Marques
,
P.
,
Magalhães
,
R.
,
Silva Moreira
,
P.
,
Miguel Soares
,
J.
, …
Kringelbach
,
M. L.
(
2017
).
Cognitive performance in healthy older adults relates to spontaneous switching between states of functional connectivity during rest
.
Scientific Reports
,
7
(
1
),
5135
. ,
[PubMed]
Caeyenberghs
,
K.
,
Verhelst
,
H.
,
Clemente
,
A.
, &
Wilson
,
P. H.
(
2017
).
Mapping the functional connectome in traumatic brain injury: What can graph metrics tell us?
NeuroImage
,
160
,
113
123
. ,
[PubMed]
Casas-Roma
,
J.
,
Martinez-Heras
,
E.
,
Solé-Ribalta
,
A.
,
Solana
,
E.
,
Lopez-Soley
,
E.
,
Vivó
,
F.
, …
Prados
,
F.
(
2022
).
Applying multilayer analysis to morphological, structural, and functional brain networks to identify relevant dysfunction patterns
.
Network Neuroscience
,
6
(
3
),
916
933
. ,
[PubMed]
Chai
,
Y.
,
Li
,
L.
,
Huber
,
L.
,
Poser
,
B. A.
, &
Bandettini
,
P. A.
(
2020
).
Integrated VASO and perfusion contrast: A new tool for laminar functional MRI
.
NeuroImage
,
207
,
116358
. ,
[PubMed]
Chai
,
Y.
,
Liu
,
T. T.
,
Marrett
,
S.
,
Li
,
L.
,
Khojandi
,
A.
,
Handwerker
,
D. A.
, …
Bandettini
,
P. A.
(
2021
).
Topographical and laminar distribution of audiovisual processing within human planum temporale
.
Progress in Neurobiology
,
205
,
102121
. ,
[PubMed]
Chen
,
V. C.-H.
,
Lin
,
K.-Y.
,
Tsai
,
Y.-H.
, &
Weng
,
J.-C.
(
2020
).
Connectome analysis of brain functional network alterations in breast cancer survivors with and without chemotherapy
.
PLOS ONE
,
15
(
5
),
e0232548
. ,
[PubMed]
Dadario
,
N. B.
,
Brahimaj
,
B.
,
Yeung
,
J.
, &
Sughrue
,
M. E.
(
2021
).
Reducing the cognitive footprint of brain tumor surgery
.
Frontiers in Neurology
,
12
,
711646
. ,
[PubMed]
De Domenico
,
M.
(
2017
).
Multilayer modeling and analysis of human brain networks
.
Gigascience
,
6
(
5
),
1
8
. ,
[PubMed]
De Martino
,
F.
,
Moerel
,
M.
,
Ugurbil
,
K.
,
Goebel
,
R.
,
Yacoub
,
E.
, &
Formisano
,
E.
(
2015
).
Frequency preference and attention effects across cortical depths in the human primary auditory cortex
.
Proceedings of the National Academy of Sciences of the United States of America
,
112
(
52
),
16036
16041
. ,
[PubMed]
De Vico Fallani
,
F.
,
Richiardi
,
J.
,
Chavez
,
M.
, &
Achard
,
S.
(
2014
).
Graph analysis of functional brain networks: Practical issues in translational neuroscience
.
Philosophical Transactions of the Royal Society B: Biological Sciences
,
369
(
1653
),
20130521
. ,
[PubMed]
Deshpande
,
G.
,
Wang
,
Y.
, &
Robinson
,
J.
(
2022
).
Resting state fMRI connectivity is sensitive to laminar connectional architecture in the human brain
.
Brain Informatics
,
9
(
1
),
2
. ,
[PubMed]
Destrieux
,
C.
,
Fischl
,
B.
,
Dale
,
A.
, &
Halgren
,
E.
(
2010
).
Automatic parcellation of human cortical gyri and sulci using standard anatomical nomenclature
.
NeuroImage
,
53
(
1
),
1
15
. ,
[PubMed]
Di
,
X.
,
Gohel
,
S.
,
Kim
,
E. H.
, &
Biswal
,
B. B.
(
2013
).
Task vs. rest-different network configurations between the coactivation and the resting-state brain networks
.
Frontiers in Human Neuroscience
,
7
,
493
. ,
[PubMed]
Ding
,
S.-L.
,
Van Hoesen
,
G. W.
,
Cassell
,
M. D.
, &
Poremba
,
A.
(
2009
).
Parcellation of human temporal polar cortex: A combined analysis of multiple cytoarchitectonic, chemoarchitectonic, and pathological markers
.
Journal of Comparative Neurology
,
514
(
6
),
595
623
. ,
[PubMed]
Duong
,
T. Q.
,
Yacoub
,
E.
,
Adriany
,
G.
,
Hu
,
X.
,
Ugurbil
,
K.
, &
Kim
,
S.-G.
(
2003
).
Microvascular BOLD contribution at 4 and 7 T in the human brain: Gradient-echo and spin-echo fMRI with suppression of blood effects
.
Magnetic Resonance in Medicine
,
49
(
6
),
1019
1027
. ,
[PubMed]
Farahani
,
F. V.
,
Karwowski
,
W.
, &
Lighthall
,
N. R.
(
2019
).
Application of graph theory for identifying connectivity patterns in human brain networks: A systematic review
.
Frontiers in Neuroscience
,
13
,
585
. ,
[PubMed]
Fatterpekar
,
G. M.
,
Delman
,
B. N.
,
Boonn
,
W. W.
,
Gultekin
,
S. H.
,
Fayad
,
Z. A.
,
Hoff
,
P. R.
, &
Naidich
,
T. P.
(
2003
).
MR microscopy of normal human brain
.
Magnetic Resonance Imaging Clinics of North America
,
11
(
4
),
641
653
. ,
[PubMed]
Felleman
,
D. J.
, &
Van Essen
,
D. C.
(
1991
).
Distributed hierarchical processing in the primate cerebral cortex
.
Cerebral Cortex
,
1
(
1
),
1
47
. ,
[PubMed]
Filippi
,
M.
,
Spinelli
,
E. G.
,
Cividini
,
C.
,
Ghirelli
,
A.
,
Basaia
,
S.
, &
Agosta
,
F.
(
2023
).
The human functional connectome in neurodegenerative diseases: Relationship to pathology and clinical progression
.
Expert Review of Neurotherapeutics
,
23
(
1
),
59
73
. ,
[PubMed]
Finn
,
E. S.
,
Huber
,
L.
, &
Bandettini
,
P. A.
(
2021
).
Higher and deeper: Bringing layer fMRI to association cortex
.
Progress in Neurobiology
,
207
,
101930
. ,
[PubMed]
Finn
,
E. S.
,
Huber
,
L.
,
Jangraw
,
D. C.
,
Molfese
,
P. J.
, &
Bandettini
,
P. A.
(
2019
).
Layer-dependent activity in human prefrontal cortex during working memory
.
Nature Neuroscience
,
22
(
10
),
1687
1695
. ,
[PubMed]
Finn
,
E. S.
,
Shen
,
X.
,
Scheinost
,
D.
,
Rosenberg
,
M. D.
,
Huang
,
J.
,
Chun
,
M. M.
, …
Constable
,
R. T.
(
2015
).
Functional connectome fingerprinting: Identifying individuals using patterns of brain connectivity
.
Nature Neuroscience
,
18
(
11
),
1664
1671
. ,
[PubMed]
Fischl
,
B.
(
2012
).
FreeSurfer
.
NeuroImage
,
62
(
2
),
774
781
. ,
[PubMed]
Fischl
,
B.
, &
Dale
,
A. M.
(
2000
).
Measuring the thickness of the human cerebral cortex from magnetic resonance images
.
Proceedings of the National Academy of Sciences of the United States of America
,
97
(
20
),
11050
11055
. ,
[PubMed]
Fischl
,
B.
,
van der Kouwe
,
A.
,
Destrieux
,
C.
,
Halgren
,
E.
,
Ségonne
,
F.
,
Salat
,
D. H.
, …
Dale
,
A. M.
(
2004
).
Automatically parcellating the human cerebral cortex
.
Cerebral Cortex
,
14
(
1
),
11
22
. ,
[PubMed]
Fleischer
,
V.
,
Radetz
,
A.
,
Ciolac
,
D.
,
Muthuraman
,
M.
,
Gonzalez-Escamilla
,
G.
,
Zipp
,
F.
, &
Groppa
,
S.
(
2019
).
Graph theoretical framework of brain networks in multiple sclerosis: A review of concepts
.
Neuroscience
,
403
,
35
53
. ,
[PubMed]
Fracasso
,
A.
,
Luijten
,
P. R.
,
Dumoulin
,
S. O.
, &
Petridou
,
N.
(
2018
).
Laminar imaging of positive and negative BOLD in human visual cortex at 7T
.
NeuroImage
,
164
,
100
111
. ,
[PubMed]
Fukutomi
,
H.
,
Glasser
,
M. F.
,
Zhang
,
H.
,
Autio
,
J. A.
,
Coalson
,
T. S.
,
Okada
,
T.
, …
Hayashi
,
T.
(
2018
).
Neurite imaging reveals microstructural variations in human cerebral cortical gray matter
.
NeuroImage
,
182
,
488
499
. ,
[PubMed]
Gau
,
R.
,
Bazin
,
P.-L.
,
Trampel
,
R.
,
Turner
,
R.
, &
Noppeney
,
U.
(
2020
).
Resolving multisensory and attentional influences across cortical depth in sensory cortices
.
eLife
,
9
,
e46856
. ,
[PubMed]
Gleichgerrcht
,
E.
,
Keller
,
S. S.
,
Drane
,
D. L.
,
Munsell
,
B. C.
,
Davis
,
K. A.
,
Kaestner
,
E.
, …
Bonilha
,
L.
(
2020
).
Temporal lobe epilepsy surgical outcomes can be inferred based on structural connectome hubs: A machine learning study
.
Annals of Neurology
,
88
(
5
),
970
983
. ,
[PubMed]
Glover
,
G. H.
,
Li
,
T. Q.
, &
Ress
,
D.
(
2000
).
Image-based method for retrospective correction of physiological motion effects in fMRI: RETROICOR
.
Magnetic Resonance in Medicine
,
44
(
1
),
162
167
. ,
[PubMed]
Greve
,
D. N.
, &
Fischl
,
B.
(
2009
).
Accurate and robust brain image alignment using boundary-based registration
.
NeuroImage
,
48
(
1
),
63
72
. ,
[PubMed]
Griffa
,
A.
,
Baumann
,
P. S.
,
Thiran
,
J.-P.
, &
Hagmann
,
P.
(
2013
).
Structural connectomics in brain diseases
.
NeuroImage
,
80
,
515
526
. ,
[PubMed]
Groppe
,
D.
(
2023
).
Bonferroni-Holm correction for multiple comparisons
. MATLAB Central File Exchange, Version 1.1.0.0. / 17 Sep 2012. https://www.mathworks.com/matlabcentral/fileexchange/28303-bonferroni-holm-correction-for-multiple-comparisons.
Groppe
,
D.
(
2024
).
fdr_bh
. MATLAB Central File Exchange, Version 2.3.0.0. / 19 December 2015.
Guidi
,
M.
,
Huber
,
L.
,
Lampe
,
L.
,
Merola
,
A.
,
Ihle
,
K.
, &
Möller
,
H. E.
(
2020
).
Cortical laminar resting-state signal fluctuations scale with the hypercapnic blood oxygenation level-dependent response
.
Human Brain Mapping
,
41
(
8
),
2014
2027
. ,
[PubMed]
Hart
,
M. G.
,
Price
,
S. J.
, &
Suckling
,
J.
(
2016
).
Connectome analysis for pre-operative brain mapping in neurosurgery
.
British Journal of Neurosurgery
,
30
(
5
),
506
517
. ,
[PubMed]
Hermundstad
,
A. M.
,
Bassett
,
D. S.
,
Brown
,
K. S.
,
Aminoff
,
E. M.
,
Clewett
,
D.
,
Freeman
,
S.
, …
Carlson
,
J. M.
(
2013
).
Structural foundations of resting-state and task-based functional connectivity in the human brain
.
Proceedings of the National Academy of Sciences of the United States of America
,
110
(
15
),
6169
6174
. ,
[PubMed]
Hilgetag
,
C. C.
, &
Barbas
,
H.
(
2006
).
Role of mechanical factors in the morphology of the primate cerebral cortex
.
PLoS Computational Biology
,
2
(
3
),
e22
. ,
[PubMed]
Holm
,
S.
(
1979
).
A simple sequentially rejective multiple test procedure
.
Scandinavian Journal of Statistics
,
6
(
2
),
65
70
. https://www.jstor.org/stable/4615733.
Honey
,
C. J.
,
Sporns
,
O.
,
Cammoun
,
L.
,
Gigandet
,
X.
,
Thiran
,
J. P.
,
Meuli
,
R.
, &
Hagmann
,
P.
(
2009
).
Predicting human resting-state functional connectivity from structural connectivity
.
Proceedings of the National Academy of Sciences of the United States of America
,
106
(
6
),
2035
2040
. ,
[PubMed]
Huber
,
L.
,
Finn
,
E. S.
,
Chai
,
Y.
,
Goebel
,
R.
,
Stirnberg
,
R.
,
Stöcker
,
T.
, …
Poser
,
B. A.
(
2021
).
Layer-dependent functional connectivity methods
.
Progress in Neurobiology
,
207
,
101835
. ,
[PubMed]
Huber
,
L.
,
Handwerker
,
D. A.
,
Jangraw
,
D. C.
,
Chen
,
G.
,
Hall
,
A.
,
Stüber
,
C.
, …
Bandettini
,
P. A.
(
2017
).
High-resolution CBV-fMRI allows mapping of laminar activity and connectivity of cortical input and output in human M1
.
Neuron
,
96
(
6
),
1253
1263
. ,
[PubMed]
Huber
,
L. R.
,
Poser
,
B. A.
,
Bandettini
,
P. A.
,
Arora
,
K.
,
Wagstyl
,
K.
,
Cho
,
S.
, …
Gulban
,
O. F.
(
2021
).
LayNii: A software suite for layer-fMRI
.
NeuroImage
,
237
,
118091
. ,
[PubMed]
Humphries
,
M. D.
, &
Gurney
,
K.
(
2008
).
Network 'small-world-ness': A quantitative method for determining canonical network equivalence
.
PLOS ONE
,
3
(
4
),
e0002051
. ,
[PubMed]
Jacob
,
Y.
,
Morris
,
L. S.
,
Huang
,
K.-H.
,
Schneider
,
M.
,
Rutter
,
S.
,
Verma
,
G.
, …
Balchandani
,
P.
(
2020
).
Neural correlates of rumination in major depressive disorder: A brain network analysis
.
NeuroImage: Clinical
,
25
,
102142
. ,
[PubMed]
Kim
,
J.
,
Wozniak
,
J. R.
,
Mueller
,
B. A.
, &
Pan
,
W.
(
2015
).
Testing group differences in brain functional connectivity: Using correlations or partial correlations?
Brain Connectivity
,
5
(
4
),
214
231
. ,
[PubMed]
Kivelä
,
M.
,
Arenas
,
A.
,
Barthelemy
,
M.
,
Gleeson
,
J. P.
,
Moreno
,
Y.
, &
Porter
,
M. A.
(
2014
).
Multilayer networks
.
Journal of Complex Networks
,
2
(
3
),
203
271
.
Klein
,
B. P.
,
Fracasso
,
A.
,
van Dijk
,
J. A.
,
Paffen
,
C. L. E.
,
Te Pas
,
S. F.
, &
Dumoulin
,
S. O.
(
2018
).
Cortical depth dependent population receptive field attraction by spatial attention in human V1
.
NeuroImage
,
176
,
301
312
. ,
[PubMed]
Kok
,
P.
,
Bains
,
L. J.
,
van Mourik
,
T.
,
Norris
,
D. G.
, &
de Lange
,
F. P.
(
2016
).
Selective activation of the deep layers of the human primary visual cortex by top-down feedback
.
Current Biology
,
26
(
3
),
371
376
. ,
[PubMed]
Kotlarz
,
P.
,
Nino
,
J. C.
, &
Febo
,
M.
(
2022
).
Connectomic analysis of Alzheimer's disease using percolation theory
.
Network Neuroscience
,
6
(
1
),
213
233
. ,
[PubMed]
Lankinen
,
K.
,
Ahlfors
,
S. P.
,
Mamashli
,
F.
,
Blazejewska
,
A. I.
,
Raij
,
T.
,
Turpin
,
T.
, …
Ahveninen
,
J.
(
2023
).
Cortical depth profiles of auditory and visual 7 T functional MRI responses in human superior temporal areas
.
Human Brain Mapping
,
44
(
2
),
362
372
. ,
[PubMed]
Lawrence
,
S. J.
,
Norris
,
D. G.
, &
de Lange
,
F. P.
(
2019
).
Dissociable laminar profiles of concurrent bottom-up and top-down modulation in the human visual cortex
.
eLife
,
8
,
e44422
. ,
[PubMed]
Lei
,
D.
,
Li
,
W.
,
Tallman
,
M. J.
,
Patino
,
L. R.
,
McNamara
,
R. K.
,
Strawn
,
J. R.
, …
DelBello
,
M. P.
(
2021
).
Changes in the brain structural connectome after a prospective randomized clinical trial of lithium and quetiapine treatment in youth with bipolar disorder
.
Neuropsychopharmacology
,
46
(
7
),
1315
1323
. ,
[PubMed]
Liégeois
,
R.
,
Santos
,
A.
,
Matta
,
V.
,
Van De Ville
,
D.
, &
Sayed
,
A. H.
(
2020
).
Revisiting correlation-based functional connectivity and its relationship with structural connectivity
.
Network Neuroscience
,
4
(
4
),
1235
1251
. ,
[PubMed]
Logothetis
,
N. K.
,
Pauls
,
J.
,
Augath
,
M.
,
Trinath
,
T.
, &
Oeltermann
,
A.
(
2001
).
Neurophysiological investigation of the basis of the fMRI signal
.
Nature
,
412
(
6843
),
150
157
. ,
[PubMed]
Mandke
,
K.
,
Meier
,
J.
,
Brookes
,
M. J.
,
O'Dea
,
R. D.
,
Van Mieghem
,
P.
,
Stam
,
C. J.
, …
Tewarie
,
P.
(
2018
).
Comparing multilayer brain networks between groups: Introducing graph metrics and recommendations
.
NeuroImage
,
166
,
371
384
. ,
[PubMed]
Mareyam
,
A.
,
Kirsch
,
J. E.
,
Chang
,
Y.
,
Madan
,
G.
, &
Wald
,
L. L.
(
2020
).
A 64-channel 7T array coil for accelerated brain MRI
.
Proceedings of the International Society for Magnetic Resonance in Medicine
,
28
(
0764
).
Markuerkiaga
,
I.
,
Barth
,
M.
, &
Norris
,
D. G.
(
2016
).
A cortical vascular model for examining the specificity of the laminar BOLD signal
.
NeuroImage
,
132
,
491
498
. ,
[PubMed]
Matkovič
,
A.
,
Anticevic
,
A.
,
Murray
,
J. D.
, &
Repovš
,
G.
(
2023
).
Static and dynamic fMRI-derived functional connectomes represent largely similar information
.
Network Neuroscience
,
7
(
4
),
1266
1301
. ,
[PubMed]
Matthews
,
P. M.
, &
Hampshire
,
A.
(
2016
).
Clinical concepts emerging from fMRI functional connectomics
.
Neuron
,
91
(
3
),
511
528
. ,
[PubMed]
Menon
,
S. S.
, &
Krishnamurthy
,
K.
(
2019
).
A comparison of static and dynamic functional connectivities for identifying subjects and biological sex using intrinsic individual brain connectivity
.
Scientific Reports
,
9
(
1
),
5729
. ,
[PubMed]
Meoded
,
A.
,
Goldenberg
,
N. A.
, &
Huisman
,
T. A. G. M.
(
2020
).
Structural connectomics: State of the art and applications in pediatric neurodevelopmental disorders, neuro-oncology, and arterial ischemic stroke
.
Journal of Pediatrics
,
221S
,
S37
S42
. ,
[PubMed]
Milano
,
M.
,
Guzzi
,
P. H.
, &
Cannataro
,
M.
(
2019
).
Network building and analysis in connectomics studies: A review of algorithms, databases and technologies
.
Network Modeling Analysis in Health Informatics and Bioinformatics
,
8
,
13
.
Miyashita
,
Y.
(
2022
).
Operating principles of the cerebral cortex as a six-layered network in primates: Beyond the classic canonical circuit model
.
Proceedings of the Japan Academy, Series B: Physical and Biological Sciences
,
98
(
3
),
93
111
. ,
[PubMed]
Moerel
,
M.
,
De Martino
,
F.
,
Uğurbil
,
K.
,
Formisano
,
E.
, &
Yacoub
,
E.
(
2018
).
Evaluating the columnar stability of acoustic processing in the human auditory cortex
.
Journal of Neuroscience
,
38
(
36
),
7822
7832
. ,
[PubMed]
Moerel
,
M.
,
De Martino
,
F.
,
Uğurbil
,
K.
,
Yacoub
,
E.
, &
Formisano
,
E.
(
2019
).
Processing complexity increases in superficial layers of human primary auditory cortex
.
Scientific Reports
,
9
(
1
),
5502
. ,
[PubMed]
Moerel
,
M.
,
Yacoub
,
E.
,
Gulban
,
O. F.
,
Lage-Castellanos
,
A.
, &
De Martino
,
F.
(
2021
).
Using high spatial resolution fMRI to understand representation in the auditory network
.
Progress in Neurobiology
,
207
,
101887
. ,
[PubMed]
Muckli
,
L.
,
De Martino
,
F.
,
Vizioli
,
L.
,
Petro
,
L. S.
,
Smith
,
F. W.
,
Ugurbil
,
K.
, …
Yacoub
,
E.
(
2015
).
Contextual feedback to superficial layers of V1
.
Current Biology
,
25
(
20
),
2690
2695
. ,
[PubMed]
Nieto-Castañón
,
A.
, &
Fedorenko
,
E.
(
2012
).
Subject-specific functional localizers increase sensitivity and functional resolution of multi-subject analyses
.
NeuroImage
,
63
(
3
),
1646
1669
. ,
[PubMed]
Norris
,
D. G.
, &
Polimeni
,
J. R.
(
2019
).
Laminar (f)MRI: A short history and future prospects
.
NeuroImage
,
197
,
643
649
. ,
[PubMed]
Olman
,
C. A.
,
Inati
,
S.
, &
Heeger
,
D. J.
(
2007
).
The effect of large veins on spatial localization with GE BOLD at 3 T: Displacement, not blurring
.
NeuroImage
,
34
(
3
),
1126
1135
. ,
[PubMed]
Osmanlıoğlu
,
Y.
,
Alappatt
,
J. A.
,
Parker
,
D.
, &
Verma
,
R.
(
2020
).
Connectomic consistency: A systematic stability analysis of structural and functional connectivity
.
Journal of Neural Engineering
,
17
(
4
),
045004
. ,
[PubMed]
Pais-Roldán
,
P.
,
Yun
,
S. D.
,
Palomero-Gallagher
,
N.
, &
Shah
,
N. J.
(
2023
).
Cortical depth-dependent human fMRI of resting-state networks using EPIK
.
Frontiers in Neuroscience
,
17
,
1151544
. ,
[PubMed]
Pervaiz
,
U.
,
Vidaurre
,
D.
,
Woolrich
,
M. W.
, &
Smith
,
S. M.
(
2020
).
Optimising network modelling methods for fMRI
.
NeuroImage
,
211
,
116604
. ,
[PubMed]
Polimeni
,
J. R.
,
Fischl
,
B.
,
Greve
,
D. N.
, &
Wald
,
L. L.
(
2010
).
Laminar analysis of 7T BOLD using an imposed spatial activation pattern in human V1
.
NeuroImage
,
52
(
4
),
1334
1346
. ,
[PubMed]
Polimeni
,
J. R.
,
Renvall
,
V.
,
Zaretskaya
,
N.
, &
Fischl
,
B.
(
2018
).
Analysis strategies for high-resolution UHF-fMRI data
.
NeuroImage
,
168
,
296
320
. ,
[PubMed]
Puetz
,
V. B.
,
Parker
,
D.
,
Kohn
,
N.
,
Dahmen
,
B.
,
Verma
,
R.
, &
Konrad
,
K.
(
2017
).
Altered brain network integrity after childhood maltreatment: A structural connectomic DTI-study
.
Human Brain Mapping
,
38
(
2
),
855
868
. ,
[PubMed]
Rockland
,
K. S.
(
2019
).
What do we know about laminar connectivity?
NeuroImage
,
197
,
772
784
. ,
[PubMed]
Rockland
,
K. S.
, &
Pandya
,
D. N.
(
1979
).
Laminar origins and terminations of cortical connections of the occipital lobe in the rhesus monkey
.
Brain Research
,
179
(
1
),
3
20
. ,
[PubMed]
Rubinov
,
M.
, &
Sporns
,
O.
(
2010
).
Complex network measures of brain connectivity: Uses and interpretations
.
NeuroImage
,
52
(
3
),
1059
1069
. ,
[PubMed]
Sadaghiani
,
S.
,
Brookes
,
M. J.
, &
Baillet
,
S.
(
2022
).
Connectomics of human electrophysiology
.
NeuroImage
,
247
,
118788
. ,
[PubMed]
Schaefer
,
A.
,
Kong
,
R.
,
Gordon
,
E. M.
,
Laumann
,
T. O.
,
Zuo
,
X.-N.
,
Holmes
,
A. J.
, …
Yeo
,
B. T. T.
(
2018
).
Local-global parcellation of the human cerebral cortex from intrinsic functional connectivity MRI
.
Cerebral Cortex
,
28
(
9
),
3095
3114
. ,
[PubMed]
Schill
,
J.
,
Simonyan
,
K.
,
Lang
,
S.
,
Mathys
,
C.
,
Thiel
,
C.
, &
Witt
,
K.
(
2023
).
Parkinson’s disease speech production network as determined by graph-theoretical network analysis
.
Network Neuroscience
,
7
(
2
),
712
730
. ,
[PubMed]
Schroeder
,
C. E.
, &
Foxe
,
J. J.
(
2002
).
The timing and laminar profile of converging inputs to multisensory areas of the macaque neocortex
.
Cognitive Brain Research
,
14
(
1
),
187
198
. ,
[PubMed]
Schroeder
,
C. E.
,
Lindsley
,
R. W.
,
Specht
,
C.
,
Marcovici
,
A.
,
Smiley
,
J. F.
, &
Javitt
,
D. C.
(
2001
).
Somatosensory input to auditory association cortex in the macaque monkey
.
Journal of Neurophysiology
,
85
(
3
),
1322
1327
. ,
[PubMed]
Setsompop
,
K.
,
Gagoski
,
B. A.
,
Polimeni
,
J. R.
,
Witzel
,
T.
,
Wedeen
,
V. J.
, &
Wald
,
L. L.
(
2012
).
Blipped-controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g-factor penalty
.
Magnetic Resonance in Medicine
,
67
(
5
),
1210
1224
. ,
[PubMed]
Shamir
,
I.
, &
Assaf
,
Y.
(
2021
).
An MRI-based, data-driven model of cortical laminar connectivity
.
Neuroinformatics
,
19
(
2
),
205
218
. ,
[PubMed]
Shamir
,
I.
, &
Assaf
,
Y.
(
2023
).
Expanding connectomics to the laminar level: A perspective
.
Network Neuroscience
,
7
(
2
),
377
388
. ,
[PubMed]
Shamir
,
I.
,
Tomer
,
O.
,
Krupnik
,
R.
, &
Assaf
,
Y.
(
2022
).
Modelling the laminar connectome of the human brain
.
Brain Structure & Function
,
227
(
6
),
2153
2165
. ,
[PubMed]
Sharoh
,
D.
,
van Mourik
,
T.
,
Bains
,
L. J.
,
Segaert
,
K.
,
Weber
,
K.
,
Hagoort
,
P.
, &
Norris
,
D. G.
(
2019
).
Laminar specific fMRI reveals directed interactions in distributed networks during language processing
.
Proceedings of the National Academy of Sciences of the United States of America
,
116
(
42
),
21185
21190
. ,
[PubMed]
Smith
,
S. M.
,
Jenkinson
,
M.
,
Woolrich
,
M. W.
,
Beckmann
,
C. F.
,
Behrens
,
T. E. J.
,
Johansen-Berg
,
H.
, …
Matthews
,
P. M.
(
2004
).
Advances in functional and structural MR image analysis and implementation as FSL
.
NeuroImage
,
23(Suppl 1)
,
S208
S219
. ,
[PubMed]
Smith
,
S. M.
,
Miller
,
K. L.
,
Salimi-Khorshidi
,
G.
,
Webster
,
M.
,
Beckmann
,
C. F.
,
Nichols
,
T. E.
, …
Woolrich
,
M. W.
(
2011
).
Network modelling methods for FMRI
.
NeuroImage
,
54
(
2
),
875
891
. ,
[PubMed]
Smith
,
S. M.
,
Vidaurre
,
D.
,
Beckmann
,
C. F.
,
Glasser
,
M. F.
,
Jenkinson
,
M.
,
Miller
,
K. L.
, …
Van Essen
,
D. C.
(
2013
).
Functional connectomics from resting-state fMRI
.
Trends in Cognitive Sciences
,
17
(
12
),
666
682
. ,
[PubMed]
Tavakol
,
S.
,
Royer
,
J.
,
Lowe
,
A. J.
,
Bonilha
,
L.
,
Tracy
,
J. I.
,
Jackson
,
G. D.
, …
Bernhardt
,
B. C.
(
2019
).
Neuroimaging and connectomics of drug-resistant epilepsy at multiple scales: From focal lesions to macroscale networks
.
Epilepsia
,
60
(
4
),
593
604
. ,
[PubMed]
Turk
,
E.
,
Scholtens
,
L. H.
, &
van den Heuvel
,
M. P.
(
2016
).
Cortical chemoarchitecture shapes macroscale effective functional connectivity patterns in macaque cerebral cortex
.
Human Brain Mapping
,
37
(
5
),
1856
1865
. ,
[PubMed]
Turner
,
R.
(
2002
).
How much cortex can a vein drain? Downstream dilution of activation-related cerebral blood oxygenation changes
.
NeuroImage
,
16
(
4
),
1062
1067
. ,
[PubMed]
Uludağ
,
K.
,
Müller-Bierl
,
B.
, &
Uğurbil
,
K.
(
2009
).
An integrative model for neuronal activity-induced signal changes for gradient and spin echo functional imaging
.
NeuroImage
,
48
(
1
),
150
165
. ,
[PubMed]
Vaiana
,
M.
, &
Muldoon
,
S. F.
(
2020
).
Multilayer brain networks
.
Journal of Nonlinear Science
,
30
(
5
),
2147
2169
.
van den Heuvel
,
M. P.
,
Scholtens
,
L. H.
,
Turk
,
E.
,
Mantini
,
D.
,
Vanduffel
,
W.
, &
Feldman Barrett
,
L.
(
2016
).
Multimodal analysis of cortical chemoarchitecture and macroscale fMRI resting-state functional connectivity
.
Human Brain Mapping
,
37
(
9
),
3103
3113
. ,
[PubMed]
van der Kouwe
,
A. J. W.
,
Benner
,
T.
,
Salat
,
D. H.
, &
Fischl
,
B.
(
2008
).
Brain morphometry with multiecho MPRAGE
.
NeuroImage
,
40
(
2
),
559
569
. ,
[PubMed]
Van Essen
,
D. C.
, &
Maunsell
,
J. H.
(
1980
).
Two-dimensional maps of the cerebral cortex
.
Journal of Comparative Neurology
,
191
(
2
),
255
281
. ,
[PubMed]
Wang
,
H. E.
,
Bénar
,
C. G.
,
Quilichini
,
P. P.
,
Friston
,
K. J.
,
Jirsa
,
V. K.
, &
Bernard
,
C.
(
2014
).
A systematic framework for functional connectivity measures
.
Frontiers in Neuroscience
,
8
,
405
. ,
[PubMed]
Wei
,
Y.
,
Scholtens
,
L. H.
,
Turk
,
E.
, &
van den Heuvel
,
M. P.
(
2018
).
Multiscale examination of cytoarchitectonic similarity and human brain connectivity
.
Network Neuroscience
,
3
(
1
),
124
137
. ,
[PubMed]
Wu
,
P.-Y.
,
Chu
,
Y.-H.
,
Lin
,
J.-F. L.
,
Kuo
,
W.-J.
, &
Lin
,
F.-H.
(
2018
).
Feature-dependent intrinsic functional connectivity across cortical depths in the human auditory cortex
.
Scientic Reports
,
8
(
1
),
13287
. ,
[PubMed]
Xia
,
M.
, &
He
,
Y.
(
2017
).
Functional connectomics from a “big data” perspective
.
NeuroImage
,
160
,
152
167
. ,
[PubMed]
Yacoub
,
E.
,
Duong
,
T. Q.
,
Van De Moortele
,
P.-F.
,
Lindquist
,
M.
,
Adriany
,
G.
,
Kim
,
S.-G.
, …
Hu
,
X.
(
2003
).
Spin-echo fMRI in humans using high spatial resolutions and high magnetic fields
.
Magnetic Resonance in Medicine
,
49
(
4
),
655
664
. ,
[PubMed]
Yeh
,
C.-H.
,
Jones
,
D. K.
,
Liang
,
X.
,
Descoteaux
,
M.
, &
Connelly
,
A.
(
2021
).
Mapping structural connectivity using diffusion MRI: Challenges and opportunities
.
Journal of Magnetic Resonance Imaging
,
53
(
6
),
1666
1682
. ,
[PubMed]
Yun
,
J.-Y.
, &
Kim
,
Y.-K.
(
2021
).
Graph theory approach for the structural-functional brain connectome of depression
.
Progress in Neuro-Psychopharmacology and Biological Psychiatry
,
111
,
110401
. ,
[PubMed]
Zachlod
,
D.
,
Rüttgers
,
B.
,
Bludau
,
S.
,
Mohlberg
,
H.
,
Langner
,
R.
,
Zilles
,
K.
, &
Amunts
,
K.
(
2020
).
Four new cytoarchitectonic areas surrounding the primary and early auditory cortex in human brains
.
Cortex
,
128
,
1
21
. ,
[PubMed]
Zaretskaya
,
N.
,
Fischl
,
B.
,
Reuter
,
M.
,
Renvall
,
V.
, &
Polimeni
,
J. R.
(
2018
).
Advantages of cortical surface reconstruction using submillimeter 7 T MEMPRAGE
.
NeuroImage
,
165
,
11
26
. ,
[PubMed]
Zeki
,
S.
(
2018
).
The rough seas of cortical cartography
.
Trends in Neurosciences
,
41
(
5
),
242
244
. ,
[PubMed]
Zhao
,
F.
,
Wang
,
P.
, &
Kim
,
S.-G.
(
2004
).
Cortical depth-dependent gradient-echo and spin-echo BOLD fMRI at 9.4T
.
Magnetic Resonance in Medicine
,
51
(
3
),
518
524
. ,
[PubMed]
Zhu
,
Y.
,
Zhu
,
G.
,
Li
,
B.
,
Yang
,
Y.
,
Zheng
,
X.
,
Xu
,
Q.
, &
Li
,
X.
(
2022
).
Abnormality of functional connections in the resting state brains of schizophrenics
.
Frontiers in Human Neuroscience
,
16
,
799881
. ,
[PubMed]

Author notes

Competing Interests: The authors have declared that no competing interests exist.

Handling Editor: Daniele Marinazzo

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. For a full description of the license, please visit https://creativecommons.org/licenses/by/4.0/legalcode.

Supplementary data