Schematic overview of the gradient mapping protocol.

For each different imaging modality (DWI; BOLD; PET) a feature matrix was built by extracting voxel-wise features in the left and right pulvinar separately. Structural and functional connectivity to 400 cortical areas, as well as the expression profiles of 19 receptors and transporters, were extracted for each pulvinar voxel. Symmetric, affinity matrices comparing the feature profile similarity of each voxel against each other were estimated using cosine similarity and fed into the diffusion embedding algorithm. Additional details can be found in the main text.

Functional connectivity, structural connectivity, and receptor coexpression pulvinar gradients.

For each panel, the top row shows the scree plots of explained variance for each gradient. The gradients that have been considered for subsequent analyses are marked with red circles. The bottom row shows 3D reconstruction of the normalized gradient images overlaid on the thalamic outline (gray-shaded area). The thalamic volume has been obtained from the AAL3 atlas. Gradients are shown in posterior (left), superior (center), and anterior (right) views. A: anterior; P: posterior; L: left; R: right.

A) The relationship between pulvinar-cortical and cortico-cortical functional gradients.

The left panels show gradient-weighted connectivity values for the first three cortico-pulvinar functional connectivity gradients, each paired with the most correlated cortico-cortical functional connectivity gradient. Normalized values of each cortical ROI (Schaefer atlas, 400 parcels) are plotted on the cortical surface. The right panels show matrix plots of Pearson’s correlation values between cortico-cortical gradients (y-axis) and pulvinar-cortical gradient-weighted connectivity maps. Only values showing statistical significance (SA-corrected, FDR-adjusted p < 0.01) are shown. B) The relationship between pulvinar-cortical functional gradients and cortical connectivity networks. Violin plots of normalized gradient-weighted connectivity values grouped by seven intrinsic connectivity networks (as in Yeo et al., 2011); VN: visual network; SMN: sensorimotor network; DAN: dorsal attention network; VAN: ventral attention network; LN: limbic network; FPN: frontoparietal network, DMN: default-mode network.

A) Pulvinar correlates of the unimodal-transmodal cortical gradient.

Scatter plots illustrate the relationship between the first pulvinar-cortical connectivity gradient, corresponding to the unimodal-transmodal hierarchy of cortico-cortical connectivity, and its most correlated gradient values across other modalities. Due to the intrinsic sign indeterminacy of gradient values, absolute correlation values are considered. B) The relationship between correlated pulvinar gradients and discrete histological nuclei. Violin plots illustrating normalized gradient values grouped by histological nuclei (AAL atlas). MED: medial pulvinar; LAT: lateral pulvinar; ANT: anterior pulvinar; INF: inferior pulvinar. C) Structural connectivity and receptor coexpression patterns correlating with the unimodal-transmodal hierarchy. Left panels: gradient-weighted structural connectivity. Normalized values for each cortical ROI (Schaefer atlas, 400 parcels) are plotted on the cortical surface. Violin plots show values grouped by seven intrinsic connectivity networks (as in Yeo et al., 2011). VN: visual network; SMN: sensorimotor network; DAN: dorsal attention network; VAN: ventral attention network; LN: limbic network; FPN: frontoparietal network; DMN: default-mode network. Right panels: gradient-correlated receptor density values for the top 5 most correlated receptors. Details can be found in the main text.

A) Pulvinar correlates of the visual-to-sensorimotor cortical gradient.

Scatter plots illustrate the relationship between the third pulvinar-cortical connectivity gradient, corresponding to the visual-to-sensorimotor hierarchy of cortico-cortical connectivity, and its most correlated gradient values across other modalities. Due to the intrinsic sign indeterminacy of gradient values, absolute correlation values are considered. B) The relationship between correlated pulvinar gradients and discrete histological nuclei. Violin plots illustrating normalized gradient values grouped by histological nuclei (AAL atlas). MED: medial pulvinar; LAT: lateral pulvinar; ANT: anterior pulvinar; INF: inferior pulvinar. C) Structural connectivity and receptor coexpression patterns correlating with the visual-to-sensorimotor hierarchy. Left panels: gradient-weighted structural connectivity. Normalized values for each cortical ROI (Schaefer atlas, 400 parcels) are plotted on the cortical surface. Violin plots show values grouped by seven intrinsic connectivity networks (as in Yeo et al., 2011). VN: visual network; SMN: sensorimotor network; DAN: dorsal attention network; VAN: ventral attention network; LN: limbic network; FPN: frontoparietal network, DMN: default-mode network. Right panels: gradient-correlated receptor density values for the top 5 most correlated receptors. Details can be found in the main text.

K-means clustering of pulvinar gradients and their correspondence with histological pulvinar nuclei.

Panels to the left show the silhouette plots for left and right pulvinar clustering solutions; error bars are standard errors calculated across 50 resamples. Panels to the right show matrix plots of Dice similarity coefficients of pulvinar clusters against histological nuclei (AAL atlas). INF: inferior; ANT: anterior; LAT: lateral; MED: medial.

Cortico-cortical functional connectivity gradients.

The left panel shows the scree plot of explained variance for each gradient, where gradients that have been considered for subsequent analyses (1-5) are marked with red circles. The right panel shows normalized gradient values for gradients 4-5 plotted on the cortical surface; Gradients 1-3 are shown in Figure 3A.

Cortico-cortical structural connectivity gradients and their relationship to pulvinar-cortical structural gradients.

The left panel shows cortico-cortical structural gradients; normalized gradient values are plotted on the cortical surface (Schaefer atlas, 400 parcels). The right panel shows scree plots of explained variance, where gradients that have been considered for subsequent analyses (1-4) are marked with red circles. Bottom right: heatmaps of Pearson’s correlation values between pulvinar-cortical gradient-weighted connectivity (x-axis) and the corresponding four cortico-cortical structural connectivity gradients (y-axis).

Reliability and reproducibility of diffusion embedding of pulvinar-cortical connectivity data. A) Split-half stability analysis.

Violin plots show the distribution of correlation values across 100 split-half resamples (bars on the median and extreme values). B) Test-retest repeatability. C) Reproducibility on the supplementary dataset (LEMON).