Cross-species functional alignment reveals evolutionary hierarchy within the connectome.

Evolution provides an important window into how cortical organization shapes function and vice versa. The complex mosaic of changes in brain morphology and functional organization that have shaped the mammalian cortex during evolution, complicates attempts to chart cortical differences across species. It limits our ability to fully appreciate how evolution has shaped our brain, especially in systems associated with unique human cognitive capabilities that lack anatomical homologues in other species. Here, we develop a function-based method for cross-species alignment that enables the quantification of homologous regions between humans and rhesus macaques, even when their location is decoupled from anatomical landmarks. Critically, we find cross-species similarity in functional organization reflects a gradient of evolutionary change that decreases from unimodal systems and culminates with the most pronounced changes in posterior regions of the default mode network (angular gyrus, posterior cingulate and middle temporal cortices). Our findings suggest that the establishment of the default mode network, as the apex of a cognitive hierarchy, has changed in a complex manner during human evolution - even within subnetworks.

Introduction

The human brain differs from other species in both scale and organization (Barton and Harvey, 2000; Krubitzer, 2009; Sousa et al., 2017; Van Essen et al., 2018). One way to understand the emergence of its unique functions is through comparative neuroimaging across different species (Mantini et al., 2012; Mars et al., 2014; Reid et al., 2016; Rilling, 2014; Van Essen et al., 2018). Evidence of similarities between neural systems in different species are assumed to reflect functions that may be relatively conserved across evolution (Kaas, 2012a; Krubitzer, 2007). In contrast, regions showing the greatest changes between humans and other species highlight neural changes that may account for features of cognition unique to humans (Ardesch et al., 2019; Buckner and Krienen, 2013; Patel et al., 2015). Traditionally, cross-species comparisons have depended on the identification of anatomical anchors (e.g. key cortical landmarks or common white matter tracts) and corresponding cortical features (e.g. myelination) in humans versus other species (Eichert et al., 2019; Goulas et al., 2014; Mars et al., 2018b; Van Essen et al., 2018; Van Essen and Dierker, 2007). These approaches have successfully identified putative homologous regions that are common to primates including macaques, marmosets and chimpanzees and humans (Chaplin et al., 2013; Donahue et al., 2018; Eichert et al., 2019; Van Essen and Dierker, 2007). These anatomy-based approaches have revolutionized our understanding of the organization of the mammalian cortex and suggest that the basic processes linked to moving and perceiving are relatively conserved across species (Hopkins et al., 2014; Krubitzer, 2007).

Although common anatomical landmarks highlight similarities in how the brain supports interactions with the immediate environment across primate species, evolution has also emphasized cognitive capabilities less closely tied to the processing of information in the ‘here and now’ (Donahue et al., 2018; Hill et al., 2010; Murphy et al., 2018; Smallwood et al., 2013; Sormaz et al., 2018). These include the ability to understand the hidden mental states of conspecifics (i.e. theory of mind), to solve problems in a creative manner, to use language to communicate intentions, and to explicitly imagine times and places not immediately available to perception (Berwick et al., 2013; Byrne, 1995; Hage and Nieder, 2016; Mansouri et al., 2017; Poerio et al., 2017). In humans, many of these processes are related to neural processing in transmodal regions of the so-called default mode and frontoparietal networks (Andrews-Hanna et al., 2014; Margulies et al., 2016). Understanding cross-species differences in transmodal regions is challenging because (i) these regions often lack clear anatomically-defined cross-species homologues (Buckner and Krienen, 2013; Mantini et al., 2013; Van Essen and Dierker, 2007) and (ii) evolution has not led to uniform cortical expansion, but has reorganized function in a heterogeneous manner that changes how regions communicate with one another (e.g. mosaic processes) (Barton and Harvey, 2000; Gómez-Robles et al., 2014; Smaers and Soligo, 2013).

A comprehensive account of how evolution has shaped cortical organization in humans, therefore, requires the mapping of how function has changed in regions of cortex where there are relatively few well-defined physical landmarks and where we may anticipate the largest cross-species differences. Contemporary accounts suggest that mammalian neural processing is organized along multiple hierarchies describing how information from distinct neural populations are integrated and segregated across the cortex (Buckner and Krienen, 2013). One important hierarchy reflects the process through which information from unimodal systems are bound together to form abstract, cross-modal, representational codes assumed to be important for multiple aspects of higher-order cognition (Mesulam, 1998). Recent observations suggest that this hierarchy is partly reflected in the intrinsic geometry of the cortex, such that regions of transmodal cortex occupy locations equidistant from unimodal systems (e.g. visual, auditory, and sensorimotor cortices) (Margulies et al., 2016). This hierarchy enables the transmodal cortex to integrate multiple sources of input and thus may underlie important aspects of higher-order cognitive functions.

To understand how evolution has shaped neural function in regions of the transmodal cortex, our study leverages recent advances in representing functional organization in a high-dimensional common space (Margulies et al., 2016). We build on prior work that shows it is possible to align across species by comparing neural activity during movie watching (Mantini et al., 2012) and that at a coarse level of analysis, similar neural hierarchies are observed in different primate species (Chaplin et al., 2013; Van Essen et al., 2018; Van Essen and Dierker, 2007). Since we were interested in charting functional differences in transmodal regions in different species we used resting-state functional connectivity data collected in two primate species (human and macaques) (Milham et al., 2018). Our connectome comparative method, which we refer to as joint-embedding, extracts the most similar dimensions of functional organization from both humans and macaques. Using this approach, we are able to empirically determine whether regions that in humans fall towards the apex of the unimodal-transmodal hierarchy have a different functional profile in macaques and thus shed important light how evolution has shaped neural function in regions that may be important for important aspects of human cognition.

Methods

Macaque data

All datasets in this study were from openly available sources. The macaque data stemmed from the recently established PRIME-DE (http://fcon_1000.projects.nitrc.org/indi/indiPRIME.html) (Milham et al., 2018). Three cohorts of macaque samples from PRIME-DE have been included in the present study. 1) Oxford data (anesthetized). The full dataset consisted of 20 rhesus macaque monkeys (macaca mulatta) scanned on a 3T with a 4-channel coil (Noonan et al., 2014). The resting-state fMRI (R-fMRI) data were collected while the animals were under anesthesia with 2 mm isotropic resolution, TR=2 s, 53.3 min (1600 volumes). No contrast-agent was used during the scans. Nineteen macaques with successful preprocessing and surface reconstruction were included in the present study (all males, age=4.01±0.98, weight=6.61 +/−2.04). 2) UC-Davis data (anesthetized). The full dataset consisted of 19 rhesus macaque monkeys (macaca mulatta) scanned on a Siemens Skyra 3T with a 4-channel clamshell coil. The resting-state fMRI data were collected with 1.4 × 1.4 × 1.4 mm resolution, TR=1.6s, 6.67 min (250 volumes) under anesthesia. No contrast-agent was used during the scans. Nineteen macaques were included in the present study (all female, age=20.38±0.93, weight=9.70±1.58). 3) Newcastle data (awake). The full data set consisted of 14 rhesus macaque monkeys (macaca mulatta) scanned on a Vertical Bruker 4.7T primate dedicated scanner (Baumann et al., 2015, 2011; Poirier et al., 2017; Rinne et al., 2017; Schönwiesner et al., 2015; Slater et al., 2016; Wilson et al., 2015). We included 10 animals (8 males, age=8.28±2.33, weight=11.76±3.38) who were scanned awake. The fMRI session was acquired with 1.2 × 1.2 × 1.2 mm resolution, TR=2s, 8.33-min per scan (250 volumes x 2 scan) for each animal. No contrast-agent was used during the scans.

Human data

The human dataset was selected from the unrelated participants of the HCP (Glasser et al., 2013). We selected the R-fMRI data from the unrelated participants (n = 178) in the HCP S500 release (Glasser et al., 2013). The first R-fMRI scan acquired on day one has been included in the current analysis, containing a 15-min run (phase encoding left-right) for each participant. The details of the acquisition and the preprocessing can be found at https://www.humanconnectome.org.

To ensure replicability, randomly split the human data into two subsets (subset HCP1, n=93, 46 females, age=29.23±3.49; subset HCP2, n=94, 36 females, age=28.99±3.43). These two subsets were grouped into two human and anesthetized macaque comparisons (HCP1-Oxford, HCP2-UCD) and two human and awake macaque comparisons (HCP1-Newcastle, HCP2-Newcastle). The following interspecies alignment analyses were replicated in all four comparison samples. We focus on the HCP1-Oxford sample in the main results and present the three other comparisons in Supplementary materials.

Preprocessing

The macaque monkey data were preprocessed using the customized HCP-like pipeline from DAF’s laboratory and the Computational Connectome System (Xu et al., 2015). The details of the data preprocessing were described previously (Xu et al., 2019, 2018). Briefly, the R-fMRI data were preprocessed including temporal compression, motion correction, 4D global scaling, nuisance regression using white matter (WM), and cerebrospinal fluid (CSF) signal and Friston-24 parameter models, bandpass filtering (0.01–0.1 Hz), detrending and co-registration to the native anatomical space. The data were then projected to the native mid-cortical surface and smoothed along the surface with FHWM=3mm. Finally, the preprocessed data were down-sampled to a 10k (10,242 vertices) resolution surface. Similar with the macaque preprocessing, the human data have been minimally preprocessed in the HCP pipeline in addition with the bandpass filtering (0.01–0.1 Hz), spatial smoothing along the surface (FWHM=6mm) and downsampling to the 10k (10,242 vertices) mid-cortical surface (Autio et al., 2019; Donahue et al., 2016).

Cross-species landmarks

The landmarks were selected based on the milestone study from Van Essen’ group (Van Essen and Dierker, 2007) and recent cross-species comparison based fMRI studies (Mars et al., 2011; Neubert et al., 2014; Sallet et al., 2013). Only potential landmarks that have been reported in at least two studies were included in the current work. The final set included 27 landmarks (Table S1). The area definition in humans was based upon the most recent multi-modal human parcellation (Glasser et al., 2016). For the landmark area in macaque, we first collected the area definitions from seven macaque atlases and used the vertices that at least overlapped within two atlases for the final macaque landmarks (Felleman and Van Essen, 1991; Ferry et al., 2000; Lewis and Van Essen, 2000; Markov et al., 2014; Paxinos and Franklin, 2012; Preuss and Goldman-Rakic, 1991; Van Essen et al., 2012). The details of the studies used to define the landmarks and the atlas references were listed in Table S1.

Joint-embedding

In previous work on manifold alignment, spectral embedding (e.g. diffusion maps) has demonstrated the ability to align the connectivity structure across individuals (Coifman and Lafon, 2006; Nenning et al., 2017). Recently, this approach has been used to characterize the connectivity topographies and capture the cortical gradients spanning along the unimodal (visual and somatomotor cortices) and transmodal regions (association cortex) within each species in human and macaque monkey (Haak et al., 2018; Margulies et al., 2016). Here, in order to align human and macaque monkey cortex, the challenge is to extract comparable cross-species components, rather than applying embeddings for each species individual and subsequently performing component matching. To address this challenge, we propose a joint-embedding approach to compute matched components (referred to as ‘gradients’) in human and macaque monkey.

First, we constructed a joint similarity matrix by concatenating within- and cross-species similarities of connectivity patterns (Fig. 1A), as defined in The diagonal within-species similarity matrices (W and W) are calculated using cosine similarity of row-wise thresholded functional connectivity at each vertex in each species (Margulies et al., 2016). The functional connectivity was calculated at the group-level by averaging the individual connectivity matrix first within each of the comparison samples. The off-diagonal cross-species similarity matrix W (and its transpose W) was calculated based on the landmark similarity profile of the functional connectivity pattern (Fig. S1). Specifically, similar to a previous study from Mars (Mars et al., 2018b), we first computed the thresholded vertex-to-vertex connectivity matrix (C and C) and averaged the vertex-wise connectivity to each landmark respectively to generate the vertex-to-landmark connectivity matrix (L and L) for each species. Based on these two connectivity matrices profile, we calculated the vertex-to-landmarks similarity matrix (S and S) within each species. That is, for each vertex within a species, the row i of matrix S is defined as the cosine similarity between row i of C and row i of L. Note that the 27 landmarks were matched homologous areas between human and macaque monkey, in other words, the columns of S and S are matched. Then we measured the cross-species similarity matrix W (and its transpose W) by comparing the similarity pattern to the homologous landmarks across species. To determine the threshold for the connectivity matrix within each species, we tested the sparsity thresholds at 1% to 10% and examined the distance of matched homologous landmarks between human and macaque in the resultant gradient space. The sparsity threshold 1% generated the most similar cross-species gradients and was employed in the final analysis.

Fig. 1.

Joint-embedding captures the common brain architecture between human and macaque monkey. A) Cross-species homologous landmarks defined by previous studies (Supplementary Table S1). B) Schematic diagram for constructing the cross-species functional common space. The joint-similarity matrix is concatenated by the vertex-wise within-species similarity matrix (diagonal) and between-species similarity matrix (off-diagonal). Spectral embedding was applied on the joint similarity matrix to extract N number of the matched components to construct the cross-species common space. C) Homologous landmark pairs are close in the joint-embedding space (Supplementary Fig. S2). D) Matched components (i.e. gradients on cortex) on human and macaque cortical surfaces highlight the homologous areas. E) The gradients were used as features in MSM for cortical surface alignment on sphere surfaces. F) The established alignment can predict the T1w/T2w map based on the other species. As visualized in a 2D density plot, the T1w/T2w in macaque (x -axis) shows significant spatial correlation with the human-to-macaque T1w/T2w prediction (y -axis).

Next, we applied the diffusion embedding algorithm on the concatenated matrix W, resulting in a set of components (Coifman and Lafon, 2006). Of note, the joint similarity matrix W is a symmetric matrix across species. The diagonal block matrices contain the within species connectivity profiles in human and macaque, encoding backbone connectivity structure (thresholded at top 1%), while the off-diagonal matrices provide a coupling across species via the comparable landmarks. Therefore, for each of the obtained components, the first half of entries correspond to the human vertices and the second half macaque vertices (Fig. 1A). Each component provides a set of matched cortical gradients covering the human and macaque cortices, which can be served as one of the dimensions of the common cross-species coordinate space. We first extracted the top 200 components and selected only the top k components to construct a gradient pool for the following surface matching procedure. Here, k is determined as the inflection point of eigenvalues (lambdas) on the scree plot (Fig. S2A). Twenty-five components (i.e. gradients) were selected in the HCP1-Oxford comparison sample (21 for HCP1-Newcastle, 18 for HCP2-UCDavis, 21 for HCP2-Newcastle).

Finally, we used the gradients from the above gradient pool as the surface features and aligned the human and macaque cortical surface with Multimodal Surface Matching (MSM) (Robinson et al., 2014). In order to avoid misalignment in the medial wall between human and macaque, we added the medial wall mask as an additional feature into MSM. The MSM configuration parameters ‘config_MSMsulc_pairwise’ was used in alignment. To optimize the number of gradients for the final alignment, we entered the top 5, 10, 15, and 20 gradients in MSM and determined the performance using 27 landmarks labels as the inspection standard. Top 15 components were finally selected for the alignment in comparison samples. Accordingly, these 15 components were used as gradient profiles to build the common coordinate space between human and macaque monkey (Fig. S3). It is worth noting that one of the macaque replication samples was scanned awake (Newcastle). To evaluate the state effect in the cross-species comparison, we calculated the similarity of the gradients’ profiles from HCP1-Newcastle cross-species analyses to those generated from the HCP1-Oxford cross-species analyses (Fig. S4). The Procrustes linear transform was applied to the raw gradients to match the order between two the two different analyses (awake human and macaque: HCP1-Newcastle, awake human and anesthetized macaque: HCP1-Oxford). The high similarity was observed in human gradients after the Procrustes transform (mean r=0.89). For macaque, the gradients were less similar between the anesthetized (Oxford) and awake (Newcastle) samples, though still relatively high (mean r=0.68).

We examined the alignment performance by applying the surface deformation to the myelin sensitive maps (i.e. T1w/T2w) and compared the aligned myelin prediction map with the actual T1w/T2w estimation in aligned species (Fig. 2D). In addition, several well-established human and macaque parcellations and networks can be registered well from human to macaque, vice versa (Fig. S5A). We also calculated the cross-species similarity matrix based on 15 gradient profiles at each vertex. We demonstrated the parcel-wise similarity matrix using the most recent multimodal parcellations for the human and its aligned human-to-macaque parcellation for the macaque (Fig. S6A). It can be seen that in general the cross-species similarity revealed greater similarity within networks than between networks (Fig. S6A).

Fig. 2.

The functional homology index (FCHI) reveals the cross-species similarity in network hierarchy. A) FCHI is calculated as the local maximum similarity of the functional gradients profile across species within corresponding searchlights (geodesic distance < 12 mm from the MSM matched vertex). B) FCHI exhibits high values in sensory cortices (e.g. visual, auditory, somatomotor), and lower in high-order association regions. C) FCHI reveals the network hierarchy (Yeo2011 networks). D) The spatial pattern of FCHI is significantly associated with the T1w/T2w map in human.

Functional Connectivity Homology Index (FCHI)

In order to quantify cross-species regional similarities of functional organization in the functional common space, we further developed the Functional Connectivity Homology Index (FCHI, Fig. 2A). Specifically, for each pair of coordinates identified as corresponding between species in MSM, we quantified the maximum cosine similarity of 15 gradients as FCHI across species within corresponding searchlights (radius = 12 mm on the midthickness surface). The searchlight approach mitigates the possibility of excessive topological constraints from MSM, while limiting the identification of matches that are unfeasible (e.g. MT and fusiform face area [FFA] in Fig. S5B). The maximum similarity within the corresponding searchlight quantified the highest likelihood that the functional gradients at each vertex in humans can be represented in macaque (Fig. 2B) and vice versa (Fig. S6B).

The activation possibility strength of BrainMap cognitive component

To quantify the relationship between functional homology and cognitive function, we employed a similar analysis as described in recent studies (Margulies et al., 2016; Wang et al., 2019). The human cognitive functions were represented using the activation possibility maps of 12 cognitive components from a previous large-scale meta-analysis based on the BrainMap database (Yeo et al., 2015). Specifically, we first grouped the macaque-to-human FCHI map into 10-percentile bins. For each of 12 cognitive components, the activation strength was normalized by dividing the sum of each component’ activation possibility and then sum within each of the 10 bins. The score in the heatmap represents the total activation possibility associated with a given cognitive component within each of the 10-percentile bin regions. The cognitive components were ordered based on the activation strength weighted by the log scale of percentile.

Evolutionary deformation and area expansion

The evolutionary surface expansion was calculated at each vertex based on the correspondence established in MSM. Specifically, we first estimated the vertex-wise surface area of the 32k standard surface mesh in native space for each of human and macaque individuals. We then resampled and smoothed (FWHM=6 mm in human and FWHM=3 mm in macaque) the area estimations to 10k surface using areal interpolation (Winkler et al., 2012). Next, the individual area maps were averaged across all the individuals to generate area maps for each of human (n=187) and macaque samples (n=48). After that, we estimated the macaque surface area at each of corresponding human vertices using the registration sphere in MSM (Winkler et al., 2012). The final relative area expansion was calculated by dividing the human surface area by the macaque surface area at each vertex on the human surface. Similarly, we calculated the relative area at each vertex on the macaque surface, suggesting the starting points of the expansion origin from macaque to human. To further demonstrate the evolutionary direction on the surface, we calculated macaque-to-human deformation vectors. To facilitate the visualization in highly folded regions (e.g. insular), we used the ‘very_inflated’ surface for both human and macaque monkey. Specifically, we identified the macaque-to-human coordinates for each of the vertices corresponding to the very_inflated macaque surface using MSM registration sphere. Next, we calculated the vector based on the MSM aligned coordinates from macaque to human. The length of the vector represents the strength of the evolutionary deformation along the direction from macaque to human surface.

Results

Joint-embedding approach captures the common brain architecture across species

To construct a common functional space for cross-species comparison, we extended the spectral embedding-based approach for mapping connectivity topographies (Haak et al., 2018; Langs et al., 2010; Margulies et al., 2016; Nenning et al., 2017). The framework of this cross-species method is visualized and introduced in the form of videos (http://fcon_1000.projects.nitrc.org/indi/PRIME/je_alignment.html).

In brief, we applied spectral embedding to a joint similarity matrix rather than computing embeddings for each species individually and subsequently performing alignment (Fig. 1). The joint similarity matrix was constructed by concatenating the following four submatrices: 1) two within-species similarity matrices (one for each species located along the diagonal), calculated using cosine similarity of thresholded functional connectivity at each vertex in each species, and 2) two off-species similarity matrices (macaque-to-human and its transpose human-to-macaque), calculated by cosine similarity of functional connectivity at each vertex with each of matched homologous landmarks, and treated as off diagonal matrices (Fig. 1A and B, Table S1) (Mars et al., 2011; Neubert et al., 2014; Sallet et al., 2013; Van Essen and Dierker, 2007). The joint-embedding results in a representation of functional connectivity shared between the two species in the form of components. Specifically, we extracted the matched components for both species, where each dimension represents the common feature along the respective embedding axis. We refer to this common space as the joint-embedding space and each dimension (component) as a ‘gradient’. The top 15 components were found to meet our minimum landmark alignment criteria (see Methods for details) and thus were retained for further analysis. We validated the matched inter-species common space by examining the similarity of the presumably anatomical landmarks in this space (Fig. 1C, Fig. S2). The Pearson correlations of the landmarks in joint-embedding space are all in the top 5% percentile of the pairwise correlations between human and macaque (ranging from r=0.877 [FEF] to r=0.999 [V1], p<0.001; Fig. S3).

The joint embedding highlights homologous regions and aligns myelin distribution across species

Our first set of analyses are aimed at establishing cross-species alignment using joint-embedding gradients. To demonstrate that the gradients can serve as common space for both species, we show that they capture the similarity in location of the well documented cross-species landmarks (Fig. 1C). After generating the joint-embedding space we projected the results (i.e. gradients) back onto human and macaque cortical surfaces. As shown in Fig. 1D, the homologous regions for both species are the apex of the same gradient (e.g. V1 as a negative nadir in gradient 1, motion-selective visual area (MT) as the positive apex in gradient 3) or occupy similar areas within the spectrum of a gradient (e.g. MT in gradient 1, V1 in gradient 3).

Next, we examined whether the joint embedding gradients provide a compact description of the distribution of myelin in both species. To establish the surface deformation between macaque and human cortex, the selected top 15 gradients as functional mesh features in Multimodal Surface Matching (MSM, Fig. 1E) for the surface registration between human and macaque (details of the model are in Methods) (Nenning et al., 2017; Robinson et al., 2014). We validated the alignment generated by joint embedding by testing its ability to generate a T1w/T2w (i.e. myelin sensitive) map for each species though predicting the T1w/T2w based on the other species following the established alignment (Glasser and Van Essen, 2011). We applied the surface human-to-macaque alignment to the myelin map that was calculated from HCP data and compared the aligned T1w/T2w predicated map to the actual macaque T1w/T2w map calculated based on Yerkes-19 template sample (Fig. 1F) (Donahue et al., 2016). The predicted map was similar to the actual T1w/T2w map (r=0.622, p<0.001 corrected). We also applied the macaque-to-human alignment to predict the human T1w/T2w using the macaque T1w/T2w map, yielding a similar, although weaker association (r=0.574, p<0.001 corrected).

Cross-species Functional Connectivity Homology Index (FCHI)

Having established that our cross species embedding adequately captures both the known homologous landmarks, and species-specific distributions of myelin in a reasonable manner, we next considered whether our approach can provide a description of how neural function differs between macaques and humans. To quantify regional similarities in functional organization across species, we developed the Functional Connectivity Homology Index (FCHI, Fig. 2A). For each pair of coordinates identified as corresponding between species using MSM, the FCHI quantifies the maximum similarity of functional gradient profiles across species within corresponding searchlights (radius = 12 mm along the surface). An advantage of using a searchlight approach is that it mitigates the possibility of excessive topological constraints from MSM, while limiting the identification of matches that are unfeasible. Fig. S5B demonstrates how the searchlight approach in identifying areas across species when their locations are anatomically decoupled (e.g. MT and fusiform face area) (Tsao et al., 2008; Yovel and Freiwald, 2013). The maximal similarity within the corresponding searchlight evaluated the highest likelihood that the functional gradients at each vertex in humans can be represented in macaque (Fig. 2B) and vice versa (Fig. S6B). These patterns were replicated in another three comparison samples and obtained the highly similar FCHI map using the awake macaque samples as well (Fig. S7).

Connectome wide differences in the organization of neural functions as demonstrated by the FCHI

The FCHI reveals the upper bounds of interspecies alignment that can be achieved based on functional organization when going from human-to-macaque, and from macaque-to-human. Here, the functional connectivity homology indices reflect the degree of common functional organization across both species. As shown in Fig. 2B, the cross-species homology index was highest in sensory areas including early visual cortices, the MT, auditory area, the fusiform face area (FFC) and the dorsal somatomotor areas. As expected, the prefrontal cortex, which is greatly expanded in human relative to the macaque, exhibited a relatively low FCHI. Importantly, a large set of high-order transmodal regions showed the lowest degree of homology index, in particular in dorsolateral prefrontal cortex (dlPFC), lateral temporal cortex (LTC), parietal gyrus regions (PG), and posterior medial cortex (PMC). This set of regions corresponds well with what is known as the default mode network (DMN) in humans. Of particular interest, the FCHI is at zero in PG regions, indicating that no regions can be found in macaques that have a functional organization similar to that of PG regions in humans. To quantify the FCHI at the level of networks, we averaged macaque-to-human similarity strength in each of seven human networks (Fig. 2C) (Yeo et al., 2011). The sensory networks, in particular the visual network, showed the highest FCHI, followed by the attention, limbic, frontoparietal networks with a moderate degree of homology index; the DMN showed the lowest homology index. Next, we examined whether cross species similarity in neural functions that are described by the FCHI are related to the cortical distribution of myelin estimated by T1w/T2w. We found that the T1w/T2w was significantly associated with the FCHI such that greater similarity between species was observed in regions that tended to be high in estimated levels of myelination (Fig. 2D, r=0.428, p<0.001, corrected). Together these analyses suggest that in functional terms regions that fall towards the apex of a unimodal-transmodal hierarchy in humans tend to have the greatest functional differences in macaques and that this difference may link to lower levels of myelin.

To evaluate the extent to which differences in FCHI along the cortical hierarchy were driven by the cortical geometry and distances to the landmarks, we measured the geodesic distance at each vertex to the nearest landmark and found a relatively low correlation (r=−0.331, Fig. S8A). It is worth noting that only the regions with the longest distance (dark red line in Fig. S8A, right panel) exhibited the skewed low FCHI scores. This suggested that the potential impact of the cortical geometry was not global but mainly driven by regions far away from the landmarks (i.e. PMC and lateral temporal lobe). We further performed a leave-one-out analysis and found that the FCHI map was highly stable across landmarks. All leave-one-out FCHI maps showed high similarity with the original FCHI map, with only validation yielding a value below 0.9 (A1: r=0.769, Fig. S8B).

In addition, we also evaluated how each of the gradients contribute to the cortical hierarchy of FCHI pattern. To this end, we calculated the FCHI based on the subset of the gradients (i.e. #1–5, #2–6, …, and #11–15) and compared to the FCHI map based on the full set of gradients (Fig. S9). We found that the first 5 gradients generated the highest similarity of FCHI pattern (r=0.864), followed by the second 5 gradients (r=0.822) and the third 5 gradients (r=0.625). The first 5 gradients also showed higher FCHI scores as compared to the last 5 gradients. These findings suggested that the early gradients captured more common functional modes between species and dominated spatial profile of the FCHI map.

The cross-species Functional Connectivity Homology Index characterizes the distributed local sensory hierarchies

Our analysis so far suggests that at a broad connectome wide level the functional differences between humans and macaques emerge in a hierarchical manner from unimodal to transmodal cortex. Next we examined whether the FCHI also can describe more local (i.e., within system) variations in cross-species functional organization. Here we focus on examples from the visual and somatomotor systems, as their hierarchical organizations are among the best understood (Buckner and Krienen, 2013; Felleman and Van Essen, 1991).

The distribution of the FCHI in the visual system can be seen in Fig. 3. It can be seen that in the early visual system, the primary visual area, V1, has the highest homology index, followed by the secondary visual area V2, the third visual area V3, and V3A; the lowest functional homologue values were observed in V4. This pattern suggests that FCHI reflects previously established visual processing order with increasing eccentricity (Felleman and Van Essen, 1991). Beyond early visual areas, in the ventral pathway, we found that the mean Functional Connectivity Homology Index progressively decreased with the complexity of information across ventral stream labeled by Brodmann Atlas areas: 17, 18, 19, 37, 20, 21 (Ungerleider and Haxby, 1994). Notably, consistent with the tethering hypothesis, which suggests that V1 and MT were both “molecular anchors ” for evolutionary processing, the FCHI for MT was comparable to V1, with lower scores being noted in surrounding areas (Fig. 3C). In somatomotor areas, the FCHI varied along the dorsal-ventral axis of somatotopic mapping (Glasser et al., 2016). The lower limb (i.e. foot) area has the greatest homology index, followed by trunk (i.e., body), and upper limb (i.e., hand); Such dorsal-to-ventral hierarchy from the lower limb to the upper limb may reflect brain adaptations for human bipedalism where the hands are more highly distinguished – for example with relatively long thumb for high precision grip in humans (Almécija et al., 2015; Oya et al., 2020). The ventral areas, eye and face/tongue areas showed the lowest FCHI (Fig. 3E). This might reflect the different evolutionary trajectory between human and NHP, as these areas evolved unique speech motor control for language in human evolution (Toda and Kudo, 2015).

Fig. 3.

Cross-species functional homology index characterizes distributed local sensory hierarchies. A) Parcel-wise FCHI maps in human and macaque monkey based on the recent parcellations (human: Glasser et al., 2016; macaque: Markov et al., 2014). B) FCHI reflects the hierarchy in the early visual processing system. Five areas defined by the Glasser parcellation are rendered on the surface and comprise the early visual areas ordered by the hierarchical streams of eccentricity mapping on x-axis (V1, V2, V3, V3A, V4). C) FCHI exhibits the highest scores in MT and MST with lower scores in MT neighboring areas (x-axis: ordered by the geodesic distance from MT area). D) FCHI decreases along the ventral visual processing hierarchy (BA17, BA18, BA19, BA37, BA20, BA 21). E) FCHI map varies along the dorsal-ventral axis of somatotopic mapping. The FCHI was averaged within labeled areas and visualized with 95% confidence intervals.

Having identified that the FCHI varies along local hierarchies, we evaluated whether this effect could be accounted for by cortical geometry rather than functional connectivity from fMRI data. To assess this, we repeated our analyses using the cortical distance matrices rather than functional connectivity data and calculated the homology index of intrinsic cortical geometry between human and macaque (Fig. S10). Notably, only in early visual areas did the geometry-based homology index demonstrate an association with hierarchy level, decreasing from V1 to higher-order areas. These results confirm that aside from early visual cortex, our primary findings above are driven by functional organization rather than cortical geometry.

The FCHI reveals the modular specialization of the subsystems in attention and frontoparietal networks

Our examination of transmodal cortex began with the frontoparietal and attention networks. In humans the frontoparietal and attention networks have been suggested to be dissociable into pairs of networks that are dissociable with respect to their functional proximity to unimodal versus transmodal systems (i.e. the default mode network) (Braga and Buckner, 2017; Dixon et al., 2018). In our analysis, the FCHI readily distinguished between the frontoparietal network-A (FN-A) and frontoparietal network-B (FN-B). Specifically, we found that FN-A which has stronger connections to the default mode network and exhibited lower FCHI scores (Fig. 4A); while, FN-B, which is more connected to the dorsal attention network in humans (Braga and Buckner, 2017; Dixon et al., 2018), exhibited higher scores. Similarly, a lower FCHI was observed in dorsal attention network-A (dATN-A) than dATN-B, which is more directly connected to retinotopic visual regions (Braga and Buckner, 2017). Finally, a similar pattern was observed in the two subnetworks in ventral attention network. Within both the frontoparietal and attention systems, therefore, we found a consistent pattern that the FCHI was lower in networks that are functionally less closely linked to unimodal systems (Caspari et al., 2015).

Fig. 4.

The functional homology index reveals the hierarchy of subsystems in the attention, frontoparietal and default mode networks. A) The FCHI shows lower scores in subsystem-A and higher scores in subsystem-B for dorsal attention, ventral attention and frontoparietal networks. B) FCHI differs among subsystems of the default mode network; the core DMN shows the lowest score, dorsal medial system an intermediate score, and medial temporal the highest score. Within each of the subsystems, the PG, PMC and LTC have the lowest FCHI scores. C) FCHI is associated with the level of ‘importance’ across subregions within the DMN (Spearman r=−0.648, p<0.001). The importance of the DMN subregions is established by a recent study based on the UK Biobank data.

The evolutionary hierarchy of subregions in default mode networks

Next, we examined the cross-species similarity of the DMN, which in humans is located at the apex of the principle cognitive hierarchy (Margulies et al., 2016). Similar to the frontoparietal network, the DMN exhibited overall differences in the FCHI among its subsystems (Andrews-Hanna et al., 2014, 2010; Braga and Buckner, 2017). In particular, the medial temporal system had the highest FCHI, the core DMN the lowest, and dorsal medial system intermediate (Fig. 4B). Across the systems, the angular gyrus (PG), posterior medial cortex (PMC) and the lateral temporal cortex (LTC) had the lowest FCHI scores (Fig. 4B).

In humans it has been recently established that regions within the DMN, primarily within the core DMN subregions, contain information regarding patterns of functional connectivity across the cortex as a whole (Kernbach et al., 2018). To understand whether cross species similarity in functional organization reflects this fine-grained distinction, we tested for associations between the FCHI and regions that reflect a high level of ‘importance’ within the DMN as identified by this prior analysis from the UK Biobank. We found that mean FCHI was negatively correlated with importance rank (Fig. 4C, Spearman r=−0.648, p<0.001), suggesting that the most important DMN subregions in the human had the least functional homology between species. Specifically, medial temporal system components, which were previously identified as having relatively low importance, were found to have high FCHI, while dorsal medial system components, which were identified as being of high importance, had low FCHI. Consistent with prior work suggesting divisions within the core subsystem with respect to high versus low importance, we found that vMPFC and DLPFC exhibited a higher FCHI, while temporal parietal junction (TPJ) and PMC had a lower FCHI (Patel et al., 2019). Together this analysis suggests that regions where the FCHI tends to be relatively lower, correspond to locations which in humans tend to be regions that are most important for reflecting global patterns of functional connectivity within the DMN (Kernbach et al., 2018).

While a number of recent studies have identified a “default-like ” transmodal network in nonhuman populations (e.g., macaque, marmoset, rodents), the extent to which this putative network functions in a manner akin to the human DMN remains an open question (Buckner and Margulies, 2019; Ghahremani et al., 2017; Hutchison and Everling, 2012; Mantini et al., 2011; Mantini and Vanduffel, 2013; Stafford et al., 2014). Here, we examined the similarity of functional organization of DMN subregions across species. This was accomplished by comparing gradient profiles in the common joint-embedding space using cosine similarity. Fig. 5A and B illustrates macaque-human similarities among DMN subregions that exceed our sparsity threshold (i.e., top 10% of pairwise human-macaque similarity over the entire cortex). Four DMN subregions in macaque (i.e., hippocampus [HC], vmPFC, dmPFC and vlPFC) were found to have functional organizations similar to those in humans. For each of the macaque DMN subregions, we identified its functional similarity with each vertex on the human cortex based on the degree of correspondence observed in their functional organization (i.e. gradient profile) (Fig. 5C). The macaque-to-human similarity maps seeded in hippocampus, vmPFC, dmPFC and vlPFC exhibited highly similar spatial patterns as human DMN. In contrast, the FCHI indicated low cross-species similarity in the macaque PMC (PCC), angular gyrus (i.e. PG), and retrosplenial cortex (RSC).

Fig. 5.

The cross-species similarity between humans and macaque monkeys. A) The human-macaque similarity of functional gradient profiles among DMN candidate subregions. The links in the diagram illustrate the similarity among DMN subregions that exceed the sparsity threshold (top 10% of pairwise human-macaque similarity across the entire cortex). B) The pairwise similarity matrix (cosine similarity) of DMN candidate subregions between humans and macaques. C) The cross-species similarity maps (cosine similarity) for each DMN candidate region seeded in macaque. The macaque-to-human similarity maps of HC, vmPFC, dmPFC and dlPFC regions show highly similar spatial patterns as human DMN (white border based on Yeo2011 networks). D) The macaque-to-human similarity maps are represented along the human principle connectivity gradient obtained based on the human HCP sample. Each line of the macaque DMN seed is smoothed by the locally weighted scatterplot smoothing (LOWESS) kernel. The positive distribution is visualized to explicitly demonstrate the extent to which DMN candidate regions in the macaque reached the human hierarchy apex.

Finally, given recent human studies suggesting the DMN is an apex transmodal network, situated at the furthest end of the macroscale sequence (Margulies et al., 2016), we examined the extent to which the DMN candidate regions in macaque may have evolved along the cortical processing hierarchy. To accomplish this, we used the human principal connectivity gradient as the reference hierarchy and computed the distribution of macaque-to-human similarity for each of macaque DMN seeds. As shown in Fig. 5D, in macaques hippocampus, vmPFC, dmPFC and vlPFC in macaque reached the human hierarchy apex, while the LTC, dPFC and dlPFC. PG, PMC and RSC, which in humans are close to the apex of the functional hierarchy, in macaques occupy lower positions in the unimodal to transmodal hierarchy.

Cross-species homology maps to cognitive functions

Together, our analysis highlights a pattern of increasing cross species differences in functional organization in regions that are thought to be the most transmodal in humans and which serve highly abstract functions. To quantify aspects of human cognition that are often associated with activity within regions of cross species difference, we compared the spatial distribution of the macaque-to-human FCHI to those provided by a large-scale meta-analysis of task fMRI experiments (Yeo et al., 2015). We first grouped the FCHI into 10-percentile bins and generated probability maps for brain activation under each that reflect the most likely cognitive functions. Probability strength for each cognitive component was normalized and averaged within each bin (see Methods). Fig. 6 shows these data in the form of a heat-map, with the components ordered in rows based on the possibility strengths (Margulies et al., 2016; Wang et al., 2019). The higher FCHI regions were associated with sensorimotor components (e.g., “visual ”, “auditory ”, “hand ”, and “face ”) whereas the lower FCHI regions were involved in high-order cognitive functions (e.g. “interoception ”, “emotion ”, “language ”, “reward ” and “dorsal attention ”). The activation for “working memory ”, “inhibition ” and “internal mentation ” were more likely to overlap with the extremes of the low FCHI regions. These findings establish that regions in which human neural organization is most different from recent common ancestors are related to a combination of executive functions and introspective processes, both of which reflect aspects of cognition that are assumed to be reasonably unique to our species.

Fig. 6.

Cross-species functional connectivity homology maps to the cognitive functions. Relationship between the FCHI map and twelve cognitive components based on the BrainMap meta-analysis database (Materials and Methods). In rows, the percentiles of the FCHI map are ordered from low to high. In columns, the cognitive components are ordered based on the normalized activation possibility strength weighted by the log scale of percentile. The higher FCHI regions were associated with sensorimotor components whereas the lower FCHI regions were involved in high-order cognitive functions.

The evolutionary surface area expansion and deformation reveals the network hierarchy

So far our analysis has established that regions of maximal cross species difference across humans and macaques, can be understood along a spectrum, from reasonably high levels of similarity in unimodal regions, to relatively low levels of similarity in regions of the transmodal cortex. This pattern was reflected as a functional shift in meta analytic data towards more abstract functions such as working memory or internal mentation that may be considered to be relatively unique to humans (Margulies et al., 2016). Theories such as the tethering hypothesis attempt to account for how the transmodal cortex gains the ability to support abstract functions, by assuming that this emerges from later development in the evolutionary process (Buckner and Krienen, 2013). To understand whether our index of cross-species similarity captures this hypothesized aspect of evolutionary change, our final analysis examined the correspondence between regions that are assumed to have expanded through evolution with the distribution of the FCHI. We first mapped the relative macaque-to-human changes in cortical area to the cortex (Fig. 7B). Sensory cortices expanded the least whereas the DMN and frontoparietal network expanded more than 20 times from macaque to human (Fig. 7B and C). In particular, the frontal cortex, temporoparietal junction, lateral temporal cortex, and the medial parietal cortex expanded the most from macaque to human. In addition, the spatial map describing surface area expansion is significantly correlated with the FCHI map (r=−.483, p<0.001 corrected) indicating that regions that have expanded the most in human are also those with the least cross species similarity as defined by the FCHI. To demonstrate the evolutionary expansion direction more explicitly, we visualized the macaque- to-human deformation vectors on the human inflated surface (Fig. 7A, center). In this figure the arrows describe the direction of change, and the color represents the degree of deformation. From macaque to human, the substantial expansion of the frontal cortex appears to push the parietal central regions in a posterior direction. The expansion of TPJ and LTC forced the temporal cortex to move posteriorly and occipital visual cortex is then squeezed into the medial side from macaque to human.

Fig. 7.

Evolutionary surface area expansion and deformation reveals the network hierarchy. A) Schematic diagram illustrating the evolutionary expansion from the macaque monkey to the human. Homologous anchors in somatomotor (area 3, area 4), primary visual areas (V1, V2) and MT areas were labeled on individual macaque and human surfaces, as well as the intermediate surfaces from macaque to human. Evolutionary expansion direction (i.e. macaque-to-human deformation vectors) is visualized in arrows on the human inflated surface (center). B) Surface areal expansion maps were calculated as the human area divided by macaque area at each of corresponding vertex on human and macaque surfaces. C) Sensory cortices expanded the least (10 times), whereas the high-order association cortex expanded the most from the macaque to the human.

Discussion

In this study, we used joint-embedding to create a common space that allowed us to assess the evolutionary changes in functional organization between species. In particular, we demonstrated differences in brain organization that have emerged through evolution can be understood in terms of variation along a hierarchy that reflects the transition of unimodal-transmodal systems. Our results reveal important clues as to how and why human cognition may differ from our close evolutionary ancestors.

We developed a Functional Connectivity Homology Index (FCHI) to quantify the likelihood that anatomically homologous regions share a common functional organization across species. At a global level, the topography of FCHI for humans and macaques had greater similarity in unimodal regions and lower similarity in systems linked to attention and more complex aspects of higher order cognition (Ardesch et al., 2019; Buckner and Krienen, 2013; Burt et al., 2018; Huntenburg et al., 2017). A more fine-grained analysis revealed important distinctions were present within canonical circuits in visual and sensory-motor territories. In these more specialized areas of cortex, we observed that the landmarks tend to fall in regions that were a local maximum for cross-species similarity, and that cross-species similarity declines in adjacent regions following well described local hierarchies (Buckner and Krienen, 2013; Kaas, 2012b; Krubitzer, 2009). This underscores the value of the landmark approach for identifying how systems directly concerned with input and output systems vary across species, while simultaneously highlighting the value of joint embedding as a means to describe functional similarity in other regions of cortex (Mars et al., 2018b; Van Essen and Dierker, 2007). Finally, we also found that the FCHI describes progressively greater differences across species in regions in which neural processing is thought to be important for more complex aspects of human cognition (e.g. attention, memory, internal mentation). Overall, our results establish that important features of the unimodal-transmodal hierarchy that have been observed in humans have emerged through evolution. This raises a number of important issues for our understanding how the organization of the cortex gives rise to uniquely human cognition.

First, using the FCHI as an index of evolutionary conservation, the present study provided clear support for mosaicism in the evolution of the human cortex, which suggests that evolutionary changes are not simultaneous across brain regions (Barton and Harvey, 2000; Gómez-Robles et al., 2014; Smaers and Soligo, 2013). Rather than simply differing across brain “modules ” or networks, the FCHI varied in a systematic fashion both between and within modules, notably in both cases indexing a hierarchical relationship. Locally within unimodal visual and sensorimotor cortical systems, the FCHI decreased away from points of well documented correspondence across species (Buckner and Krienen, 2013; García-Cabezas and Zikopoulos, 2019). Globally the FCHI showed a consistent spatial variation with cortical myelin that increased in the networks anchored in the unimodal cortex and decreased in networks that are important higher order functions (Margulies et al., 2016; Paquola et al., 2019). Within the DMN those regions with low FCHI are the same regions that in humans contain the most information about global brain dynamics. Together these results are consistent with a complex view of cross-species differentiation which impacts on cortical modules at both a local and global scale. They also highlight that at least parts of the mechanisms that lead to the apparent mosaic-like changes in cortical organization that emerged through evolution is also reflected in the neural hierarchy observed in humans (Burt et al., 2018; Goulas et al., 2014; Paquola et al., 2019; Sneve et al., 2018; Wang et al., 2019).

Second, the distribution of cross-species similarities identified by the FCHI is consistent with the tethering hypothesis (Buckner and Krienen, 2013; Kaas, 2012b; Krubitzer, 2009). We found that both the early visual cortex and the ventral visual pathway anchored from V1, and the hierarchy organization seen in MT and neighboring areas, were reflected as decreases in the FCHI. The tethering hypothesis suggests new functional capabilities have arisen through the gradual duplication, budding, and subdivision of brain areas (Buckner and Krienen, 2013). Specifically, through a process of cortical expansion, large parts of the cortical mantle are postulated to have progressively become untethered from direct roles in input and output system as they become increasingly distant from the constraints of molecular gradients that surrounded key “anchor ” regions (e.g., V1 and MT) (Buckner and Krienen, 2013; Rosa, 2002; Van Essen et al., 2018). We found that the FCHI reflects this hierarchy at both the local and global scale. For example, in canonical circuits in sensory and motor cortex, the observed radial pattern in the FCHI supports the view that evolutionary changes are increasingly important outside of the primary cortex. We also found the FCHI was lower in regions in which cortical expansion across species is thought to be largest (Mars et al., 2018b; Van Essen et al., 2018). Taken together, our findings support the tethering hypothesis in multiple functional aspects.

Third, our study provides important insight into the evolution of the commonly described DMN as is seen in humans. We found that nodes within the DMN show an important differential relationship in terms of their cross-species similarity. In particular, while some DMN regions are functionally similar in both species, two core regions of the DMN exhibited low similarity between humans and macaques: the PMC (PCC-PCU) and PG (parietal gyrus and angular gyrus). In contrast, the vmPFC shows a more similar organizational profile across species (Amiez et al.; Ghahremani et al., 2017). In humans it has recently been established that these regions with the lowest FCHI (e.g. PG and PCC) provide the greatest information about connectome wide patterns of information flow (Kernbach et al., 2018). Moreover, recent task-fMRI studies have determined that these regions can play a crucial role in the application of task-relevant information during working memory and cognitive flexibility tasks (Murphy et al., 2019, 2018; Vatansever et al., 2017). Together these results suggest that in functional terms, some aspects of what is called the default mode in humans are present in macaques, while others, particularly the posterior regions, have a much more unique functional profile in humans (Andrews-Hanna et al., 2014; Kernbach et al., 2018). One implication of the distribution of the FCHI across the DMN is that while the foundations of this system may be present in many species, the pattern seen in humans in which both anterior and posterior regions act in concert, may have emerged relatively later during human evolution. As a consequence, the functional pattern that represents the DMN as seen in humans, may be present in only a subset of our recent ancestors. It is also possible that these, as well as higher nonhuman primate species (e.g., chimpanzee) would exhibit transitional or intermediate DMN variants.

There are a number of issues that should be considered when interpreting our results. Our cross species joint-embedding is part of an emerging trend toward the use of common high-dimensional spaces as a tool to understand the evolutionary changes across species. Recent efforts have developed strategies for alignment based on white matter tracts, or myelin (T1w/T2w) maps (Eichert et al., 2019; Mars et al., 2018b), as well as task-based activations during movie viewing (Mantini et al., 2012). We anticipate that the joint-embedding approach used in this paper should be readily extensible to diffusion imaging data and encourage future work in this direction (Mars et al., 2018a, 2016). Regardless of which data modality is used, the use of a high-dimensional common space may help characterize functional similarities and differences across species, particularly in areas where clear anatomical landmarks are difficult or impossible to ascertain. In principle, this approach could be extended to multiple species (e.g., chimpanzee, baboon, marmoset etc.) to allow alignment across multiple non-human species including human primates, as well as other mammals, for example rodents. Future work using multimodal imaging in multiple species is required to disentangle function and anatomy in our understanding of brain evolution and development (Heuer et al., 2019; Heuer and Toro, 2019). It is also important to note that our analyses were carried out at the group-level. This decision was motivated by our desire to maximize signal to noise ratio, though proper implementation of this method in individual animals will require the additional optimization of methods (e.g., sufficient data collection for individual dataset, generation of individual-specific landmark masks) (Croxson et al., 2018; Xu et al., 2019, 2018). In addition, while the present work addressed issues regarding the reproducibility of findings by replicating analyses in independent samples, our understanding of their functional significance is based on a meta-analysis and so the precise functional meaning of the observed differences remains largely a matter of conjecture (Braga and Buckner, 2017; Laumann et al., 2015). In order to fully appreciate the significance of these cross-species differences in brain organization for human cognition, it will be necessary to understand how these patterns change across contexts, and, if possible, during active tasks states (Petit and Pouget, 2019; Sharma et al., 2019). Understanding how functional patterns change across a common embedded space during periods of task engagement, could provide invaluable insight into how evolution has shaped many important aspects of human cognition (Murphy et al., 2018).

In summary, our results provide novel insights into the mode in which evolutionary changes sculpt the cerebral cortex layout. Using the FCHI, a quantitative index of neural similarity across species, we established that important changes in how cortical regions communicate has emerged in a progressive fashion along the spectrum spanned by unimodal and transmodal regions, which reflect the structural and functional gradients that pertain to the human cerebral cortex. Importantly, our findings highlight the way that evolutionary changes might have contributed to the emergence of uniquely human higher-order cognition, and provide potential insights into the evolutionary roots of sub- systems within the attention, frontoparietal and default mode networks. Additionally, the posterior DMN, as the apex of a cognitive hierarchy, may have unique evolutionary adaptations and changed substantially from the most recent common ancestor of humans and macaques.