Fig. 3: Brain-wide structural variation in gene expression in macaques.

a Pairwise differential expression matrix across the brain. Each matrix entry represents the upregulated (bottom left) and downregulated (top right) number of DEGs with a fold change ≥ 2 in expression level. Subcortical and occipital regions exhibit marked heterogeneous patterns (red) in contrast to the rest cortical regions (blue). DEGs, differentially expressed genes. b Heatmap of the numbers of DEGs between major brain structures (left) and cortical regions (right). Occipital regions exhibit marked inter-regional variations (red). GP globus pallidus, Tha thalamus, Pu putamen, Cd caudate, NA nucleus accumbens, Amy amygdala, HC hippocampus, CTX cortex, Pir piriform cortex. c Increased number of detected DEGs when increased the number of randomly-sampled areas within each lobe/region (color coded) (Permutation test, n = 1000). d Average number of DEGs in each lobe/region. In each box plot, the center line indicates the median, the white diamonds indicate the mean, the edges of the box indicate the 25th and 75th percentile (interquartile range, IQR) and the whiskers indicate last point within a 1.5× IQR (n = 561, frontal; 55, parietal; 351, temporal; 21, cingulate; 15, occipital; 36, insula). *p < 0.05, **p < 0.01, ***p < 0.001, pairwise two-sided Wilcoxon rank sum test (uncorrected for multiple comparisons), p values in d are provided as a Source Data. Source data are provided as a Source Data file.