
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39231965
51942
10.1038/s41467-024-51942-1
Article
Sex differences in functional cortical organization reflect differences in network topology rather than cortical morphometry
http://orcid.org/0009-0004-8142-1655
Serio Bianca b.serio@fz-juelich.de

1234
http://orcid.org/0000-0002-7973-6752
Hettwer Meike D. 1234
http://orcid.org/0000-0001-8006-8678
Wiersch Lisa 125
Bignardi Giacomo 36
Sacher Julia 3478
http://orcid.org/0000-0001-7726-6470
Weis Susanne 12
Eickhoff Simon B. 123
http://orcid.org/0000-0003-2998-6849
Valk Sofie L. valk@cbs.mpg.de

1234
1 https://ror.org/02nv7yv05 grid.8385.6 0000 0001 2297 375X Institute of Neuroscience and Medicine, Brain & Behavior (INM-7), Research Centre Jülich, Jülich, Germany
2 https://ror.org/024z2rq82 grid.411327.2 0000 0001 2176 9917 Institute of Systems Neuroscience, Medical Faculty, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
3 grid.4372.2 0000 0001 2105 1091 Max Planck School of Cognition, Leipzig, Germany
4 https://ror.org/0387jng26 grid.419524.f 0000 0001 0041 5028 Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
5 https://ror.org/05dnene97 grid.250903.d 0000 0000 9566 0634 Brain-Based Predictive Modeling Lab, Feinstein Institutes for Medical Research, Glen Oaks, New York, NY USA
6 https://ror.org/00671me87 grid.419550.c 0000 0004 0501 3839 Language and Genetics Department, Max Planck Institute for Psycholinguistics, Nijmegen, The Netherlands
7 grid.411339.d 0000 0000 8517 9062 Leipzig Center for Female Health & Gender Medicine, Medical Faculty, University Clinic Leipzig, Leipzig, Germany
8 grid.411339.d 0000 0000 8517 9062 Clinic for Cognitive Neurology, University Medical Center Leipzig, Leipzig, Germany
4 9 2024
4 9 2024
2024
15 77146 12 2023
21 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Differences in brain size between the sexes are consistently reported. However, the consequences of this anatomical difference on sex differences in intrinsic brain function remain unclear. In the current study, we investigate whether sex differences in intrinsic cortical functional organization may be associated with differences in cortical morphometry, namely different measures of brain size, microstructure, and the geodesic distance of connectivity profiles. For this, we compute a low dimensional representation of functional cortical organization, the sensory-association axis, and identify widespread sex differences. Contrary to our expectations, sex differences in functional organization do not appear to be systematically associated with differences in total surface area, microstructural organization, or geodesic distance, despite these morphometric properties being per se associated with functional organization and differing between sexes. Instead, functional sex differences in the sensory-association axis are associated with differences in functional connectivity profiles and network topology. Collectively, our findings suggest that sex differences in functional cortical organization extend beyond sex differences in cortical morphometry.

Here, the authors show sex differences in the functional sensory-association axis are associated with differences in functional connectivity profiles and network topology rather than reflecting differences in brain size, microstructural organization, or geodesic distance.

Subject terms

Computational neuroscience
Network topology
Computational neuroscience
https://doi.org/10.13039/501100002347 Bundesministerium für Bildung und Forschung (Federal Ministry of Education and Research) https://doi.org/10.13039/501100004189 Max-Planck-Gesellschaft (Max Planck Society) issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Sex differences in human global brain size are robust and widely acknowledged1–7, but the functional implications of this anatomical difference are not well understood. Indeed, sex differences in intrinsic brain function are sometimes deemed small or negligible beyond differences attributed to brain size2. Nevertheless, diverging patterns of functional connectivity between males and females have been reported even when controlling for differences in brain size and most consistently in sensory and association regions5,8,9. These regions in fact represent the two anchors of a key principle of hierarchical functional organization, the sensory-association (S-A) axis, differentiating localized primary sensory/motor areas from a more distributed set of transmodal association regions, including regions belonging to the frontoparietal and default mode networks (DMN)10,11. However, the extent to which sex differences in intrinsic functional cortical organization may be explained by neuroanatomical differences relating to brain size remains unclear.

In order to understand how sex differences in brain size may pertain to sex differences in brain function, it is first necessary to understand the relevance of brain size for overall functional cortical organization. Brain size and its variability may have important consequences for the spatial distribution of sensory and association areas across the cortical mantle, as illustrated by clear scaling patterns over evolution12,13 and development14. In fact, over the past 4 million years, hominin evolution has not only shown a general trend of increasing body mass, but also an even more important relative increase in brain size12. According to the tethering hypothesis, the brain’s sensory systems, acting as anchors, may have exerted constraining pressures during the growth of the developing ancestral mammalian cortex13. In this way, evolutionary cortical expansion may have led to the emergence of the S-A axis, with association cortices distributed across the cortical mantle and untethered from sensory hierarchies. Patterns of expansion across cortical regions along the S-A axis are also observed across human development, with a more markedly distributed areal expansion across frontoparietal association regions relative to limbic and sensorimotor areas14—see also15 for a comprehensive review of the S-A axis’ neurodevelopment. Through the increase of overall brain size, the differential expansion of sensory and association areas could thus be an important product of mammalian evolution and development, further reflecting the hierarchical functional differentiation of these regions. In fact, this different scaling and reorganization of regions along the S-A axis also appears to reflect patterns of sex differences in cortical morphometry, that is, cortical shape and size.

Morphometric differences between male and female brains have been extensively reported, with males showing a greater absolute brain volume6. Different measures of brain size are commonly used in the literature (such as intracranial volume, total brain volume, and total surface area). Although these measures highly covary and are often used interchangeably, they quantify different morphometric features of the brain, with sex differences in “brain size” ranging from 8% to 13% depending on the selected measure6. The size and direction of sex effects also vary by neuroanatomical property, such as different tissue types, brain regions, and features (including cortical thickness, gyrification, and surface area)16. It must be noted that within-group variability in cortical morphometry –which is typically greater in males– is larger than between-group mean effects, meaning that individual differences within sex are larger than group-differences between sexes17. Although individual differences in total brain size seem to account for most differences in relative regional volumes3, some sex differences still remain statistically significant when the variance explained by total brain size is taken into account7. Therefore, there may be sex differences in the scaling of regional brain volume that go beyond linear associations with overall brain and body size. In fact, although local sex differences in cortical morphometry are typically small in size and vary across studies18, they have been reported in both sensory and association regions: A meta-analysis identified volumetric sex differences in multiple cortical regions, including the anterior and posterior cingulate gyri, precuneus, right frontal pole, inferior and middle frontal gyri, insular cortex, Heschl’s gyrus, and lateral occipital cortex6. Another study found greater gray matter volume in females in prefrontal and superior parietal cortices, whilst males showed greater volumes in ventral occipitotemporal regions19. These findings further depict the apparent relevance of the S-A axis as an axis of morphometric variability between the sexes. Developmental trajectories of anatomical change also appear regionally heterogeneous along the S-A axis, with higher rates of global cortical thickness change found in fronto-temporal association regions and lower rates found in sensory regions20. Morphometric cortical properties therefore seem to not only follow patterns of variation along the S-A axis, but also differ between the sexes. Yet, how exactly sex-specific differences in cortical morphometry may be relevant to differences in intrinsic brain function has not been directly explored.

Consistent with patterns of morphometric variation and sex differences, evidence points to sex differences in intrinsic functional connectivity (FC) in regions represented at the poles of the S-A axis5,8,9. In fact, despite generally controversial findings on sex differences in brain function, findings of stronger FC in females within the DMN5,21–24 and stronger FC in males within sensorimotor areas5,22,25 are consistent. Overlapping morphometric and functional patterns of sex differences along the S-A axis thus suggest that differentiation in functional cortical organization may be somewhat orchestrated by the cortical mantle’s morphometric properties. Indeed, the structure, size, and shape of the cortex not only physically support functional connections, but also determine their length. Short- and long-range connections, as measured by geodesic distance (the distance separating two regions along the cortical mantle) have in fact been found in sensory and association regions respectively26, thus also displaying patterns of variation along the S-A axis. With increasing distance between regions, cortical function also appears to change more rapidly in association regions relative to sensorimotor regions27. These patterns further mirror patterns of microstructural cortical variability identified by post-mortem histology26 and myelin-sensitive in vivo magnetic resonance imaging (MRI)26,27. As such, intrinsic functional activity, showing variability between the sexes and along the S-A axis, seems to be embedded within the cortical mantle and its microstructural organization. Accumulating evidence further supports the important role played by cortical geometric properties, including size and shape, in sculpting functional architecture. Established findings from graph theory suggest that a cortical functional network’s properties are largely determined by its spatial embedding, namely by the length of its connections28. Peaks of DMN clusters on the S-A axis also appear to be equidistantly distributed relative to primary areas10, in line with the hypothesized untethering of association cortices from sensory hierarchies during evolutionary expansion13. Furthermore, recent findings suggest that the spatial organization of intrinsic cortical functional activity is dominated by long wave-lengths of geometric eigenmodes29. This research builds on notions from neural field theory positing that brain shape physically constrains brain-wide functional dynamics by imposing boundaries on emerging functional signals30,31. In the context of sex differences in functional cortical organization, brain size also explains some—although not all—sex-specific variance in FC32. Together, these findings point to possible morphometric properties that may not only underpin cortical functional architecture, but also be at the root of sex differences in functional cortical organization.

In the current work, we therefore investigated the extent to which sex differences in intrinsic functional cortical organization may be associated with differences in cortical morphometry, namely different measures of brain size, microstructure, and the geodesic distance of connectivity profiles. To this end, we used multimodal imaging data (including resting state functional MRI and structural T1 and T2 images) of the Human Connectome Project (HCP) S1200 release33, consisting of healthy young adults self-reporting their biological sex. We began by computing the S-A axis as our measure of functional organization, given that it reflects multimodal mechanisms bridging morphometric, structural, and functional features of cortical organization, and given that sex differences in functional connectivity are typically found in regions situated at the extremities of this axis, i.e., in sensory and association regions. We then tested for sex differences along this low dimensional axis of hierarchical organization. Next, we identified the cortical morphometric properties potentially constraining the S-A axis, including different measures of brain size (recognizing total surface area as the most pertinent to the S-A axis), microstructural organization (a low dimensional microstructure profile covariance (MPC) axis), and the mean geodesic distance of connectivity profiles. We then probed associations between patterns of sex differences in cortical morphometry and patterns of sex differences in the S-A axis. Contrary to our expectations, we do not find evidence supporting a morphometric explanation of sex differences in functional organization, despite identifying sex differences in the morphometric properties per se. As such, we further probed potential functional features that may intrinsically underpin sex differences on the S-A axis, and our findings suggest that differences in FC profiles and network topology may be a more plausible explanation of sex differences in functional organization.

Results

Sex differences in the S-A axis of functional cortical organization

In order to investigate whether sex differences in functional organization may be related to sex differences in cortical morphometry, we began by constructing our measure of functional organization and testing for related sex differences. We thus computed the S-A axis at the individual level as our measure of functional organization in subjects of the HCP S1200 release33. For this, we applied a non-linear dimensionality reduction algorithm on FC Fisher r-to-z transformed matrices. We only considered the top 10% of the row-wise z values, representing each seed region’s top 10% maximally functionally connected regions34,35. We used this 90% threshold for consistency with previous studies10,36,37 and for its high test-retest reliability and reproducibility38,39. As such, we computed the well-replicated low dimensional axis of functional brain organization explaining the most variance in the data (21.86%)—spanning from unimodal (sensory, here particularly visual) regions to heteromodal (association) regions10—and defined it as the S-A axis (Fig. 1A). Mean cortex-wide spatial patterns of S-A axis loadings were overall highly correlated between the sexes, rspin = 0.996, pspin < 0.001. Then, to test for regional effects of sex on S-A axis loadings, we fitted a linear mixed effects model (LMM) including fixed effects of sex, age, and total surface area, and random nested effects of family relatedness and sibling status (see Methods for more information on the nested structure of the HCP data and the statistical modeling). We identified sex differences in the S-A axis loadings in 23.3% of cortical regions (93 out of 400 Schaefer parcels) with small to medium effect sizes (minimum effect size of β = 0.210; maximum effect size of β = 0.435), which were distributed across the seven intrinsic functional Yeo networks40 (Fig. 1B, C). Positive standardized beta coefficients (β values), depicted in blue, represent higher loadings in males relative to females on the S-A axis, whereas negative β-values, depicted in red, represent higher loadings in females relative to males. In Supplementary Fig. 1, we show that patterns of within-sex variability in S-A axis loadings are similar between males and females, with only a few regions showing statistically significant sex differences in variance.Fig. 1 The sensory-association (S-A) axis of functional cortical organization and its sex differences.

A Mean S-A axis loadings (spanning from visual to default-mode network regions) across sexes; B Thresholded β-map of linear mixed effect model (LMM) results showing false discovery rate (FDR)-corrected (q < 0.05) statistically significant effects of sex on S-A axis loadings, where blue represents higher male loadings and red represents higher female loadings; C Functional network breakdown of cortical areas showing statistically significant FDR-corrected sex differences in S-A axis loadings. The outer ring displays absolute proportions of statistically significant cortical areas by functional Yeo network, the inner ring displays absolute nested proportions by directionality of effects, where blue represents higher male loadings and red represents higher female loadings. β standardized beta coefficient. Source data are provided as a Source Data file.

Morphometric correlates of the S-A axis

We then investigated potential morphometric constrains on functional organization by probing associations between the S-A axis and different measures of brain size, microstructural organization, and the mean geodesic distance of connectivity profiles.

First, we tested for associations between the S-A axis loadings and three measures of brain size commonly used in the literature, namely intracranial volume (ICV), total cortical volume (TCV), and total surface area. More specifically, ICV represents the entire volume encapsulated by the cranium (i.e., including cerebrospinal fluid), TCV represents the total volume of gray and white matter within the neocortex (excluding subcortical structures), and total surface area represents the entire surface area of the neocortical mantle (see Methods for the exact computation of these measures). Sex differences in these measures of brain size and other anthropometric measurements (i.e., height, weight, and body mass index) are further reported in Supplementary Table 1. For each measure of brain size, we fitted an LMM to test for regional effects of brain size on S-A axis loadings (Supplementary Fig. 2), and we found total surface area to have the most widespread effects amongst the three tested brain size measures (Fig. 2D; Supplementary Fig. 2C).Fig. 2 Morphometric correlates of the sensory-association (S-A) axis of functional cortical organization across sexes.

A Mean S-A axis loadings (spanning from visual to default-mode regions) across sexes; B Mean microstructural profile covariance (MPC) axis loadings (spanning from sensory to paralimbic regions) across sexes; C Mean geodesic distance of connectivity profiles across sexes; D Thresholded β-map of linear mixed effect model (LMM) results showing false discovery rate (FDR)-corrected statistically significant effects (q < 0.05) of total surface area on the S-A axis loadings; E Thresholded β-map of LMM results showing FDR-corrected statistically significant effects of MPC axis loadings on the S-A axis loadings; F Thresholded β-map of LMM results showing FDR-corrected statistically significant effects of mean geodesic distance on the S-A axis loadings; G Spatial correlation between mean patterns of S-A axis loadings and mean patterns of MPC axis loadings (color-coded by functional Yeo network), tested by a two-sided Spearman correlation and corrected for spatial autocorrelation, r = 0.20, pspin = 0.037. Error band displays 95% confidence interval; H Spatial correlation between mean patterns of S-A axis loadings and mean patterns of mean geodesic distance (color-coded by Yeo network), tested by a two-sided Spearman correlation and corrected for spatial autocorrelation, r = 0.76, pspin < 0.001. Error band displays 95% confidence interval. β standardized beta coefficient. Source data are provided as a Source Data file.

Second, we computed the microstructural profile covariance (MPC) axis of organization at the individual level, which is a low dimensional representation of the similarity of microstructural intensity profiles, defined as the T1-weighted (T1w) over T2-weighted (T2w) tissue intensity ratio, across cortical regions and layers37,41,42. We computed the MPC axis by conducting nonlinear dimensionality reduction on MPC matrices34,35, which were obtained by sampling and correlating the intracortical microstructural intensity of 12 equivolumetric depth profiles (see Methods). Following the same approach used for computing the S-A axis, we selected the axis explaining the most variance in the data (25.81%)—spanning from sensory to paralimbic regions– defining it as the MPC axis (Fig. 2B). We specifically selected this low-dimensional representation of microstructural organization as it has been previously shown to covary with the low-dimensional representation of functional organization (i.e., the S-A axis)42. To test for whole-brain associations between the S-A and MPC axes, we correlated the spatial maps of the axes’ mean loadings (Fig. 2A, B) across all subjects (Fig. 2G; r = 0.20, pspin = 0.037). We further fitted an LMM to test for regional effects of MPC axis loadings on S-A axis loadings at the parcel level (Fig. 2E), and found small and localized associations between the S-A and MPC axes.

Third, we computed the mean geodesic distance of connectivity profiles at the individual level. The mean geodesic distance of connectivity profiles is the mean distance along the cortical mantle between each region and its top 10% maximally functionally connected regions. Group-level patterns (i.e., averaged across all subjects; Fig. 2C) revealed shorter geodesic distances in visual and somatomotor (sensory) regions, and longer distances in frontoparietal and DMN (association) regions. We also tested for whole-brain associations between the S-A axis and patterns of mean geodesic distance of connectivity profiles by correlating their spatial maps (Fig. 2A, C) averaged across all subjects (Fig. 2H; r = 0.76, pspin < 0.001). We also fitted an LMM to test for regional effects of mean geodesic distance on S-A axis loadings at the parcel level (Fig. 2F) and found strong and widespread associations between the S-A axis and mean geodesic distance of functional connectivity profiles.

Sex differences in morphometric correlates are not associated with sex differences in the S-A axis

After establishing the cortical morphometric correlates of the S-A axis, we probed whether sex differences in cortical morphometry may reflect sex differences in the S-A axis. To this effect, we first computed sex differences in the morphometric correlates of the S-A axis using independent LLMs. Males displayed a larger total surface area (1947.35 ± 1542.24 cm2) relative to females (1715.07 ± 1469.42 cm2), β = 1.21, p < .001, as reported in Supplementary Table 1. We also identified bidirectional statistically significant regional sex differences in MPC axis loadings (Fig. 3B) as well as in the mean geodesic distance of functional connectivity profiles (Fig. 3C).Fig. 3 Sex differences in the morphometric correlates of the sensory-association (S-A) axis.

A Thresholded β-map of linear mixed effect model (LMM) results showing false discovery rate (FDR)-corrected statistically significant effects (q < 0.05) of sex on the S-A axis, where blue represents higher male loadings and red represents higher female loadings; B Thresholded β-map of LMM results showing FDR-corrected statistically significant effects of sex on the microstructure profile covariance (MPC) axis; C Thresholded β-map of LMM results showing FDR-corrected statistically significant effects of sex on the mean geodesic distance of connectivity profiles. β standardized beta coefficient. Source data are provided as a Source Data file.

We then tested whether the identified sex differences in cortical morphometry may be associated with sex differences in the S-A axis. First, we tested whether sex differences in S-A axis loadings were moderated by total surface area. For this, we modeled an interaction term of sex by total surface area on the S-A axis loadings (Fig. 4A) and found no statistically significant effects across cortical regions. In Supplementary Fig. 3A-C, we further show that this interaction effect yields virtually the same effects (β values) when including versus excluding height as a covariate from the LMM (r = 0.99, pspin < 0.001). This suggests that height—being an anthropometric feature that systematically differs between the sexes—does not explain variance in the moderation of sex effects by total surface area on the S-A axis loadings either. We also plotted within-sex effects of total surface area on S-A axis loadings, showing similar although slightly diverging patterns of effects between males and females (r = 0.66, pspin = 0.002; Supplementary Fig. 3D–F). However, the divergence of patterns between sexes may not be strong or systematic enough to be interpreted as meaningful, as underlined by the lack of statistically significant sex by total surface area interaction effects on the S-A axis.Fig. 4 Associations between sex differences in sensory-association (S-A) axis and sex differences in cortical morphometry.

A Unthresholded β-map of linear mixed effect model (LMM) testing for sex by total surface area interaction effects on S-A axis (there were no statistically significant sex by total surface area interaction effects after false discovery rate (FDR) correction (q < 0.05)); B Unthresholded β-map of LMM testing for sex by MPC axis interaction effects on S-A axis (there were no statistically significant sex by MPC axis interaction effects after FDR correction); C Unthresholded β-map of LMM testing for sex by mean geodesic distance interaction effects on S-A axis (35 cortical areas showed statistically significant sex by mean geodesic distance interaction effects after FDR correction and are delineated in black); D Scatterplot displaying the spatial correlation between patterns of sex differences (β-maps) in S-A axis loadings and in MPC axis loadings (color-coded by Yeo network), tested by a two-sided Spearman correlation and corrected for spatial autocorrelation, r = −0.05, pspin = 0.398. Error band displays 95% confidence interval; E Scatterplot displaying the spatial correlation between patterns of sex differences (β-maps) in S-A axis loadings and in the mean geodesic distance of connectivity profiles (color-coded by Yeo network), tested by a two-sided Spearman correlation and corrected for spatial autocorrelation, r = −0.07, pspin = 0.326. Error band displays 95% confidence interval. β standardized beta coefficient. Source data are provided as a Source Data file.

We then tested for spatial associations between regional sex effects on S-A axis loadings (Fig. 3A) and regional sex effects on MPC axis loadings (Fig. 3B) by correlating the two β-maps. Here, we found no statistically significant association between these two patterns of sex differences (Fig. 4D; r = −0.05, pspin = 0.398). Similarly, we tested for spatial associations between regional sex effects on S-A axis loadings (Fig. 3A) and regional sex effects on the mean geodesic distance of connectivity profiles (Fig. 3C). Again, we found no statistically significant association between these two patterns of sex differences (Fig. 4E; r = −0.07, pspin = 0.326). These results together suggest that sex differences in the S-A axis are not overall systematically associated with sex differences in cortical morphometry. In addition to probing spatial associations between sex differences in cortical morphometry and sex differences in the S-A axis, we also modeled interaction terms of sex by MPC axis loadings (Fig. 4B) and sex by mean geodesic distance (Fig. 4C) to test for these interaction effects on S-A axis loadings. Here, we found no statistically significant sex by MPC axis interaction effects, but find 35 cortical regions showing statistically significant sex by geodesic distance interaction effects, specifically in the bilateral medial prefrontal cortex, temporal regions, and left dorsolateral regions (delineated in black on Fig. 4C).

As an additional sensitivity analysis, we found that including the MPC axis and mean geodesic distance as covariates to an LMM testing for sex effects on the S-A axis yields highly similar regional sex effects to those reported in Fig. 1 (for which the original LMM only included total surface area as a morphometric covariate, to control for brain size), as shown by the strong correlation of β-maps (r = 0.95, pspin < 0.001, Supplementary Fig. 4A). Similarly, the association between sex effects when including all morphometric covariates versus not including any (i.e., also excluding total surface area) remains high despite a small decrease in correlation strength (r = 0.80, pspin < 0.001, Supplementary Fig. 4B). These findings further suggest that sex differences in total surface area—representing brain size—only explain a minor amount of variance in sex differences in the S-A axis.

Intrinsic functional underpinnings of sex differences in the S-A axis

Given that sex differences in the morphometric correlates of the S-A axis did not appear to reflect sex differences in the S-A axis, we probed potential intrinsic functional underpinnings of sex differences in the S-A axis. Concretely, we investigated which sex differences in features of FC may be conserved through data reduction and reflected in the sex differences that we observe in the S-A axis loadings. For this, we tested for associations between sex differences in the S-A axis loadings and sex differences in three intrinsic features of FC, namely mean FC strength, FC profiles, and network topology.

First, we computed mean FC strength at the individual level from FC matrices, representing—for each parcel– the mean row-wise z values of each given seed region’s top 10% maximally functionally connected regions. We then fitted an LMM to test for local effects of sex on mean FC strength (Fig. 5A, B), which revealed—amongst other sex differences—greater female intrinsic FC in DMN regions and greater male intrinsic FC in somatomotor regions. To test associations between patterns of sex differences in the S-A axis loadings (Fig. 1B) and in FC strength (Fig. 5A), we spatially correlated the β-maps of the respective sex effects and did not detect a statistically significant association between sex differences in the S-A axis and sex differences in FC strength (Fig. 5C; r = −0.02, pspin = 0.483).Fig. 5 Intrinsic sex differences in functional connectivity (FC) strength.

A Thresholded β-map of linear mixed effect model (LMM) results showing false discovery rate (FDR)-corrected statistically significant effects (q < 0.05) of sex on mean FC strength; B Functional network breakdown of connections showing statistically significant FDR-corrected sex differences in mean FC strength. The outer ring displays absolute proportions of statistically significant cortical areas by functional Yeo network, the inner ring displays absolute nested proportions by directionality of effects, where blue represents greater male FC strength and red represents greater female FC strength; C Spatial correlation between patterns of sex effects in S-A axis loadings and patterns of sex effects in mean FC strength (color-coded by Yeo network), tested by a two-sided Spearman correlation and corrected for spatial autocorrelation, r = −0.02, pspin = 0.483. Error band displays 95% confidence interval. β standardized beta coefficient. Source data are provided as a Source Data file.

Second, we investigated whether sex differences in the S-A axis may be related to sex differences in FC profiles. We defined FC profiles at the individual level, for which we identified the top 10% maximally functionally connected regions. Using the Chi-square (χ2) test of independence, we assessed—for each possible pairwise connection along the 400 × 400 matrix—sex differences in a given target region’s odds of belonging to the top 10% maximally functionally connected regions of a given seed region. In this way, we specifically considered the data—at the individual level—that our data reduction algorithm has been applied on in order to probe whether the top 10% connections made by females and males may have differentially influenced the computation of the S-A axis and consequently of sex differences in its loadings. We identified the direction of these sex effects with the odds ratio (OR), where OR < 1 indicates a given region’s greater female odds and OR > 1 indicates a given region’s greater male odds. Out of the 160,000 tested functional connections, 2004 connections (corresponding to 1.25% of all connections) showed statistically significant sex differences in their odds of constituting a seed’s top 10% connections after FDR correction, suggesting that sex differences in S-A axis loadings may in part stem from differences in FC profiles, namely differences in which functional connections are the strongest. For connections showing statistically significant sex differences, we found an OR ranging from 0.00–0.64 in the case of greater female odds (OR < 1), and an OR ranging from 1.56 to 25.36 in the case of greater male odds (OR > 1). For illustrative purposes, we summarized spatial patterns of sex differences in FC profiles as the sum of connections showing sex differences per seed region (Fig. 6A), as well as the overall networks involved in sex differences in FC profiles (Fig. 6B). Figure 6C displays the connections between seed and target regions showing statistically significant sex differences in their FC profiles.Fig. 6 Intrinsic sex differences in functional connectivity (FC) profiles and network topology.

A Number of connections (per seed region) showing statistically significant false discovery rate (FDR)-corrected sex differences in their odds of belonging to the given seed’s top 10% connections; B Functional network breakdown of connections showing statistically significant FDR-corrected sex effects in their odds of belonging to the given seed’s top 10% connections. The outer ring displays absolute proportions of statistically significant cortical areas by functional Yeo network, the inner ring displays absolute nested proportions by directionality of effects, where blue represents higher male odds and red represents higher female odds; C Connections between seed and target regions showing statistically significant FDR-corrected sex differences in FC profiles (odds ratio (OR) > 1 meaning that males have higher odds than females of having a target region belong to a seed region’s top 10% connections; OR < 1 meaning that females have higher odds than males of having a target region belong to a seed region’s top 10% connections; connections are color coded by Yeo network and weighed by number of connections between the network pairs; D β values of linear mixed effect model (LMM) results for the sex contrast in between-network (BN) dispersion for each pairwise Yeo network comparison, where blue represents higher male BN dispersion and red represents higher female BN dispersion (no statistically significant sex effects after spin permutation and multiple comparisons Bonferroni correction; two-sided pspin < 0.001); E β values of LMM results for the sex contrast in within-network (WN) dispersion for each Yeo network (displayed as white dots), plotted on null distributions of β values derived from 1000 spin permutations, where positive β values represent higher male WN dispersion and negative β-values represents higher female WN dispersion. * indicates multiple comparisons Bonferroni-corrected (pspin < 0.004) two-sided statistical significance of the sex contrast for dispersion in the DMN, β = 0.241, pspin < 0.001. β standardized beta coefficient. Source data are provided as a Source Data file.

Finally, we investigated sex differences in network topology, namely the organization of functional networks along the S-A axis. We computed measures of between-network dispersion, quantifying the pairwise distance between all networks along the S-A axis, where a higher value indicates greater segregation of the given pair of networks and a lower value indicates greater integration of the given pair of networks. Here, we computed between-network dispersion for each possible network pair (21 pairs of Yeo networks in total) at the individual level43. We also computed a measure of within-network dispersion40, quantifying the spread of regions within a network along the S-A axis, where a higher value indicates greater segregation of the given network’s regions and a lower value indicates greater integration of the given network’s regions. Here, we computed within-network dispersion for all seven Yeo networks at the individual level. LMMs did not show any statistically significant sex difference in between-network dispersion for any of the network pairs (Fig. 6D). However, we found greater male within-network dispersion in the DMN, β = 0.241, pspin < 0.001 (Fig. 6E), revealing a greater spread of regions belonging to the DMN along the S-A axis in males. The full statistical results for the analysis of sex differences in network dispersion are summarized in Supplementary Table 2.

Discussion

In the current work, we investigated the extent to which sex differences in functional cortical organization may be associated with differences in cortical morphometry, namely different measures brain size (focusing on total surface area), microstructure, and the geodesic distance of connectivity profiles. We identified widespread sex differences in young adult functional cortical organization as defined by the S-A axis, which however did not appear to be systematically associated with sex differences in total surface area, microstructural organization, nor the mean geodesic distance of connectivity profiles. This finding is particularly striking given that the morphometric properties under study were all per se associated with the S-A axis and differed between sexes. We observed that sex differences in the S-A axis were instead related to differences in FC profiles and network topology, namely greater male dispersion within the DMN. Collectively, our findings suggest that sex differences in functional cortical organization extend beyond neuroanatomical sex differences pertaining to cortical morphometry.

Considering the common use of different measures of brain size in the literature, which yield different magnitudes of sex differences depending on the selected measure6, we tested the effects of different measures of brain size on the S-A axis, namely ICV, TCV, and total surface area. Here, given that total surface area had the most widespread effects on functional organization, we deemed it the most appropriate measure of brain size for our study and further included it as a covariate in our models throughout our analyses. The relevance of total surface area for our study is also supported by the theoretical assumptions motivating our research question, namely the relevance of cortical shape and geometry in constraining brain wide functional dynamics29–31 and thus sex differences in these features potentially underpinning sex differences in the S-A axis. In showing the diverging statistical effects of different measures of brain size, our findings highlight the risk of introducing noise when including an inadequate measure of brain size, particularly when statistically controlling for brain size in the detection of sex effects on brain structure and function32,44–47. The complex heterogeneity of neuroanatomical properties constituting brain size should not be undermined, also considering that morphometric features vary differently as a function of age, whereby for example total brain volume but not ICV is affected by atrophy6. As such, future research on sex differences should also carefully select the measure of brain size that is most conceptually and empirically pertinent to the research question under study in order to avoid introducing noise in the analyses.

By establishing morphometric correlates of the S-A axis in addition to total surface area, namely a microstructural axis of cortical organization37,42,48 and the mean geodesic distance of connectivity profiles10,48, our findings align with previous work and argue for the rooting of functional cortical organization in cortical structure and shape. We show a particularly strong association between the mean geodesic distance of connectivity profiles and the S-A axis, supporting the relevance of the cortical mantle’s shape in sculpting functional organization. This may be a product of the cortical mantle’s evolutionary expansion, where association regions are untethered from sensory hierarchies13, and long-range connections preserve the overall connectedness of cortical networks by facilitating the communication between distant areas28. Furthermore, as indexed by the MPC axis, microstructural organization appears to mildly covary with the S-A axis, supporting to some degree the well-established idea of structural constraints on brain function37,42,49. In our study, we obtained intensity profiles via the ratio of T1w over T2w imaging sequences, and although it is commonly used to measure myelin37,42,50, the T1w/T2w ratio has been described as an acceptable qualitative proxy for myelin in gray but not white matter51. It is indeed thought to capture unique features of microstructural tissue that appear largely independent of diffusion-based metrics, thus portraying a mix of neuroanatomical features beyond pure myelin52. We therefore consider the T1w/T2w ratio—and the resulting MPC axis—as a general measure of tissue microstructure, which may serve as a scaffold for functional organization.

After establishing morphometric correlates of the S-A axis, we addressed our primary aim of probing the extent to which sex differences in functional cortical organization may be reflected by sex differences in cortical morphometry. Our findings overall suggest that morphometric differences between the sexes are altogether not substantial contributors to sex differences in the S-A axis of functional organization. We did not find any statistical spatial associations between patterns of sex differences in the S-A axis and patterns of sex differences in the MPC axis nor in the mean geodesic distance of connectivity profiles. Although we observed slightly diverging results when including—as opposed to excluding—total surface area as a covariate in our model testing for sex differences in S-A axis loadings, we did not find a statistically significant interaction between sex and total surface area in reflecting S-A axis loadings. We also found no statistically significant sex by MPC axis interactions but found few statistically significant sex by mean geodesic distance interactions. This can be understood as partly mirroring the sex differences in functional connectivity profiles that we further explain below, given that geodesic distance was averaged based on the connectivity profiles—more specifically the top 10% functional connections—which we find sex differences in. Altogether, the negligeable relevance of cortical morphometry to sex differences in the S-A axis is striking given that morphometric properties appear per se to be associated with the S-A axis and to differ between sexes. The biological mechanisms underpinning different patterns of morphometric and functional sex differences may thus be independent from one another, suggesting that sex differences in functional cortical organization may extend beyond the connectome’s supporting shape and structure.

Given that sex differences in morphometric correlates of the S-A axis did not seem to be associated with sex differences in the S-A axis, we probed and found potential intrinsic functional underpinnings of sex differences in the S-A axis. Firstly, the sex differences that we observed in the S-A axis loadings were distributed across functional networks, notably in the DMN, frontoparietal, and ventral attention networks. This is consistent with a previous study in youth reporting that these association networks show greater individual variability in their functional topography relative to lower-order sensory networks, whilst contributing the most to sex classification8. We also observed sex differences in intrinsic FC strength, replicating previous widely established findings of greater FC in females within DMN regions21–23 and in males within somatomotor regions22,25. However, these patterns did not spatially overlap with patterns of sex differences in the S-A axis, suggesting that FC strength is not a feature of intrinsic FC that is captured by sex differences in our low dimensional representation of functional organization. Instead, we found that sex differences in the S-A axis were related to modest sex differences in FC profiles, which also presented qualitative sex differences in the proportional breakdown of networks involved. The strongest functional connections (top 10% connections) of females seemed to involve the DMN more than in males, whereas males displayed more top connections involving the somatomotor networks relative to females. This is also consistent with the few regional statistically significant interaction effects of sex by mean geodesic distance of FC profiles on S-A axis loadings. Indeed, given that sensory and association regions are respectively known to have primarily short- versus a mix of short- and long-range connections26, sex differences in FC profiles involving regions in the somatomotor network and the DMN may be related to differences in the geodesic distances of functional connectivity profiles between sexes. As such, our findings suggest that sex differences in the S-A axis may be better represented by sex differences in FC profiles, i.e., the topology of functional connections, than sex differences in FC strength alone. Furthermore, sex differences in the configuration of functional connections may not only underly the recurrence of sex differences in these networks21–23,25, but may also reflect sex differences in network topology.

We observed greater male dispersion relative to females within the DMN. This finding suggests that areas belonging to the DMN are represented further apart on the S-A axis (i.e., showing less similarity in their FC profiles) in males relative to females, which is also consistent with previous findings of generally more segregated male networks53. These network-specific topological sex differences may be related to greater female odds of connections within the DMN, and greater male odds of somatomotor connections with other networks. As such, the apparent sex differences in network topology and functional connectivity profiles—albeit small—provide a more interpretable system-level description of the sex differences observed in the S-A axis loadings, representing a key principle of macroscale cortical organization. Concretely, network topology, which represents the organization of functional communities within and between functional networks43, may reflect brain states54—and in our case, possible differences thereof at the group-level. Network topology has also been associated with different cognitive features including arousal55, awareness and consciousness56, behavior and task performance57, and cognitive flexibility58. The balance between integration and segregation is complex, dynamic, and necessary to maintain the brain’s metastability59 by reaching a point of equilibrium between global organization and local specialization49. The brain is a highly interconnected and metabolically expensive organ, and its organization is required to dynamically balance topological efficiency and energy utilization in response to transient cognitive and physiological demands60. Our findings of sex differences in network topology may therefore pertain to intricate sex differences not only in group-averaged brain states at rest, but also in global energy expenditure, which would reflect physiological differences.

Despite the insights gained through our study, some limitations must be acknowledged. Firstly, by focusing on self-reported biological sex, we did not test for effects of gender or gender-related variables on functional organization and its morphometric correlates. This is relevant given that there are multiple factors contributing to sex/gender differences, including both biological and social factors. Findings may indeed appear more nuanced if we move beyond the unrealistic assumption of a clear-cut sexual dimorphism of brain structure and function61, further considering the relevance of gender, steroid hormones, and the role of X and Y chromosomes62. Nevertheless, we intentionally focused on the dichotomous variable of sex to identify group-level effects, as our study aimed to investigate the correspondence (i.e., shared variance) between sex differences in cortical morphometry and sex differences in cortical functional organization, regardless of variance explained by gender. We recognize the limitations of using a binary variable given that differences between groups may be attributable to both sex- and gender-based influences62. However, quantifying these influences was beyond the scope of this study. Secondly, we focused on neocortical functional organization, excluding subcortical structures and the cerebellum despite their substantial contributions to whole brain organization through their notable structural integration with the cortex63. For example, the amygdala and hippocampus have shown a variable degree of both structural6 and functional64 sex differences. Nevertheless, our exclusive focus on the neocortex was motivated by the relevance of using the S-A axis as our measure of functional organization, which is commonly obtained by reducing the dimensionality of FC matrices of cortical data10. By using the S-A axis, our work took a system-level approach to identify sex differences embedded in a key macroscale organizational principle, going beyond previous research on functional differences between the sexes solely focusing on intrinsic brain function. Thirdly, the morphometric properties considered in our study are not exhaustive, overlooking the contributions of other morphometric measures such as local volumes of gray matter. The inclusion of the MPC axis37,42 and the mean geodesic distance of connectivity profiles10,48 was however supported by their theoretical and empirical relevance to functional cortical organization, particularly its low dimensional embedding. Finally, our findings are limited to sex differences in a healthy sample and would benefit from being replicated in a more inclusive sample that is more representative of the overall population. This would be additionally informative considering the notable sex differences observed in populations that would typically be excluded from healthy samples, for example individuals with neurodevelopmental and psychiatric disorders65. Nevertheless, the structure-function associations that we investigated in this work are rather fundamental, and their essence should thus be fairly well captured in our large (healthy) sample.

All in all, our study gives rise to a set of questions pertaining to the mechanisms underpinning sex differences in functional cortical organization, given that they do not appear to be rooted in cortical morphometric differences. Our findings instead suggest that sex differences in the S-A axis are, to some extent, intrinsically related to differences in FC profiles and network topology. Although these sex differences appear to be small, they may be meaningful for broader sex differences in functional cortical organization, and future research should explore factors driving males and females to form these few distinct functional connections that are associated with sex differences in the system-level organization of functional networks, notably of the DMN. Recognizing the human body as a complex system of systems, future work should also investigate other biological factors that may contribute to functional sex differences such as genes located on sex chromosomes19 and steroid hormones66,67. Environmental factors should equally be considered, as they may not only differ on average as a function of sex and gender, but may also differently affect brain function across the sexes through divergent mechanisms68. An example of this is stress, whereby sex differences in the stress response have been found to contribute to sex differences in brain function and psychopathology via epigenetic mechanisms69–71. Generally, it is crucial for research in neuroscience to systematically test for sex differences in brain structure and function, as well as their biological and environmental underpinnings, in order to produce more rigorous and representative findings, ultimately leading to a more translational body of knowledge68.

Methods

The current study complies with ethical regulations set by The Independent Research Ethics Committee at the Medical Faculty of the Heinrich-Heine-University of Duesseldorf.

Participants and study design

Our analyses were conducted on the publicly available data of healthy young adults from the Human Connectome Project (HCP) S1200 release (http://www.humanconnectome.org/)33. We selected subjects with available functional, T1, and T2 neuroimaging data, resulting in a final sample of 1000 individuals (536 females) with a mean age of 28.73 ± 3.71 years, used for all analyses. The sample included 284 monozygotic twins (MZ), 184 dizygotic twins (DZ), 443 non-twin siblings, and 89 unrelated individuals, and the sociodemographic breakdown of the participants is additionally reported in Supplementary Table 3. Subjects were all born in Missouri but recruited in an attempt to broadly reflect the racial and ethnic composition of the United States population. Recruitment efforts aimed to yield a subject pool capturing a wide range of variability—in socioeconomic and behavioral terms—in order to be representative of the general healthy population. The term “healthy” was thus broadly defined. Individuals with documented neurodevelopmental and psychiatric disorders, or reporting physiological illnesses such as high blood pressure or diabetes were excluded from the HCP study recruitment protocol, but not individuals who reported smoking, being overweight, or a history of recreational drug use or heavy drinking (if they had not experienced severe symptoms). Informed consent was obtained for all study subjects and recruitment procedures were approved by the Washington University institutional review board. More detailed information about the HCP study design and recruitment procedure is available elsewhere33,72. Despite the use of the term gender in the HCP Data Dictionary, we use the term sex in this article given that the HCP study collected self-reported information on biological sex instead of gender identification, as reported elsewhere73. We have not used genetic information to verify the self-reported sex.

Structural MRI acquisition and preprocessing

The HCP’s MRI data was acquired on a customized 3T Siemens Skyra ConnectomeScanner with a 32-channel head coil at Washington University across four scanning sessions held over two days. Structural MRI images were acquired on the same day via high resolution T1-weighted (T1w) and T2-weighted (T2w) sequences. Two separate T1w images were acquired and averaged, with identical scanning parameters using a 3D MPRAGE sequence (0.7 mm isovoxels, FOV = 224 mm, matrix = 320 × 320 mm, 256 sagittal slices; TR = 2400 ms, TE = 2.14 ms, TI  =  1000 ms, flip angle = 8°, BW = 210 Hz per pixel, ES = 7.6 ms). Two separate T2w images were acquired and averaged, with identical scanning parameters using a variable flip angle turbo spin-echo (3D T2-SPACE) sequence, with the same isotropic resolution, matrix, FOV, and slices as for the T1w sequence (TR = 3200 ms, TE = 565 ms, BW = 744 Hz per pixel, total turbo factor = 314). The preprocessing steps included co-registering the T1w and T2w images, bias field (B1) correction, registration to MNI space, segmentation, and surface reconstruction. See refs. 33,72,74 for more detail on the HCP’s MRI protocols and the FreeSurfer segmentation pipeline.

Functional MRI (fMRI) acquisition and preprocessing

The HCP’s fMRI data was collected after the structural sequences and following the HCP’s minimal processing pipeline, as described above. A total of 1 h of resting-state functional data was collected across four identical 15 min scanning sessions, equally split over two days (LR1, RL1, LR2, RL2), with a gradient echo EPI sequence at a resolution of 2 mm isotropic (FOV = 208 × 180 mm, matrix = 104 × 90 mm, 72 slices covering the whole brain, TR = 720 ms, TE = 33 ms, multiband factor of 8, FA = 52°). The multimodal surface matching algorithm (MSMAll) was used to co-register the data to the HCP template 32k_LR surface space, consisting of 32,492 nodes per hemisphere (59,412 nodes excluding the medial wall). A more detailed description of the resting state fMRI data acquisition and analysis protocol is available elsewhere74,75.

Functional connectivity (FC) and the sensory-association (S-A) axis of functional organization

Throughout this work, we used the Schaefer 400 parcellation (clustered into 7 networks: visual, somatomotor, dorsal attention, ventral attention, limbic, frontoparietal, DMN40). This widely used functionally-derived parcellation scheme was originally obtained via a gradient-weighted Markov Random Field model integrating local gradient and global similarity approaches76. The vertex-wise functional timeseries were therefore averaged within the Schaefer 400 cortical regions. FC matrices (400 × 400) were then computed at the individual level—per scanning session—by correlating cortical timeseries in a pairwise manner using the Pearson product moment. We normalized the correlation coefficients using Fisher’s z-transformation. Final FC matrices were obtained by averaging each subject’s matrices across their four scanning sessions. From these FC matrices and for each subject, we computed the S-A axis of functional organization, as described below.

We conducted data reduction on the FC matrices to yield macroscale gradients of functional organization10. For this, we used diffusion map embedding, a nonlinear manifold learning algorithm that reduces complex, high-dimensional structures of data (in our case affinity matrices) to low-dimensional representations combining geometry with the probability distribution of data points34. Thus, cortical regions that are strongly interconnected (i.e., whose timeseries show high correlations) are represented closer together in the resulting low dimensional manifold of FC data, whereas parcels with low covariance are represented farther apart, as indexed by the cortical regions’ gradient loadings. To this end, we used the BrainSpace Python toolbox35 to compute ten gradients with the following parameters: 90% threshold (i.e., only considering the top 10% row-wise z values of FC matrices, representing each seed region’s top 10% maximally functionally connected regions), α = 0.5 (α controls whether the geometry of the set is reflected in the low-dimensional embedding— i.e., the influence of the sampling points density on the manifold, where α = 0 (maximal influence) and α = 1 (no influence)), and t = 0 (t controls the scale of eigenvalues). These parameters were selected for consistency with previous studies and represent choices that are recommended to retain global relations between datapoints in the embedded space whilst being relatively robust to noise10,42. To confirm the robustness of the S-A axis computed at the 90% threshold, we further show with a sensitivity analysis that the mean S-A axis computed at the 90% threshold shows high correlations with mean S-A axes computed at different thresholds (from 10% to 90%, in steps of 10%), with r values ranging from 0.84 to 0.93 (see Supplementary Fig. 5). In order to increase comparability for further between-subject analyses, Procrustes alignment was used to align individual gradients to mean gradients. Mean gradients were computed by applying diffusion map embedding—with the same parameters listed above—to the mean FC matrix (i.e., FC matrices averaged across all subjects). The computation of these FC gradients was carried out independently per hemisphere (i.e., considering the top 10% row-wise z values of only half of the FC matrices, shaped 200 × 200) and the gradient loadings resulting from each hemisphere were subsequently concatenated. This decision was made for consistency and comparability reasons within our study, so that the top 10% functional connections selected for data reduction corresponded to those considered in the calculation of the mean geodesic distance of connectivity profiles—which were only computed per hemisphere– as described further below). We verified and confirmed the stability FC gradients when computing them per hemisphere versus at the whole brain level, as shown by the spatial correlation of mean gradient loadings (r = 0.98, pspin < 0.001). Finally, we took the well-replicated principal gradient explaining the most variance in the data and spanning from visual to DMN regions10, which we labeled the S-A axis and used to represent functional organization for subsequent analyses.

We also computed, for each subject, mean FC strength at the parcel level in a seed-wise fashion, by averaging the row-wise z values of each seed region’s top 10% maximally functionally connected regions—again per hemisphere—and subsequently concatenated the hemispheric mean FC strength values to reconstruct whole brain data.

Cortical microstructure and microstructural profile covariance (MPC)

Microstructural properties—including myelin and cellular characteristics—show depth-dependent variation along cortical columns, as reported by histology42,77,78 as well as in vivo and post mortem neuroimaging37,41,42,78, which illustrate a cortical hierarchy11. Similar to previous work37, we quantified cortical microstructure, or “microstructural profile intensity” (MPI), using the myelin-sensitive MRI contrast obtained from the T1w/T2w ratio from the HCP minimal processing pipeline described above74 (a reliability check is reported in the Supplementary Fig. 6). The T1w/T2w ratio uses the T2w image to correct for inhomogeneities in the T1w image50. Then, we followed the previously described protocol37,41,42 to compute our measurement of MPC, which reflects the variation of MPI, across cortical depths. In short, we generated 14 equivolumetric surfaces within the inner and outer cortical surfaces, then excluded the inner- and outer-most surfaces, thus remaining with 12 surfaces representing cortical layers. Surface generation was based on a model compensating for cortical folding by altering the pairwise Euclidean distance (ρ) of intracortical surfaces throughout the cortex and thus preserving fractional volume between the surfaces. For each surface, ρ was calculated as defined in Eq. 1.1 ρ=1Aout−Ain⋅−Ain+αAout2+1−αAin2

for which α denotes a fraction of the total volume of the segment that the surface accounts for, while Aout and Ain respectively denote the surface areas of the outer and inner cortical surfaces.

Across the whole cortex and from the outer to the inner surfaces, we systematically sampled MPI values layer-wise for each of the 64,984 vertices of the HCP template 32k_LR surface space, which we then averaged within each of the 400 Schaefer cortical areas, per layer. Following a previously described protocol41, we constructed subject level 400 × 400 matrices using pairwise Pearson partial correlation on the MPI profiles of cortical parcels (i.e., correlating the MPI values across 12 layers between parcels), controlling for overall mean cortical MPI, followed by log transformation. We then used these matrices to compute MPC gradients by following the same procedure and using the same toolbox and parameters as for computing the FC gradients34,35, as described above and previously done37,41,42. For consistency and comparability with the computation of FC gradients and mean geodesic distance of connectivity profiles, we also computed MPC gradients independently per hemisphere and subsequently concatenated the gradient loadings resulting from each hemisphere. We verified and confirmed the stability MPC gradients when computing them per hemisphere versus at the whole brain level, as shown by the spatial correlation of mean gradient loadings (r = 0.99, pspin < 0.001). We also selected the principal gradient of MPC explaining the most variance in the data and spanning from sensory to paralimbic regions, which we labeled the MPC axis and used to represent microstructural organization in subsequent analyses.

Measures of brain size

In our analyses we included different measures of brain size typically used in the literature, including intracranial volume (ICV), total cortical volume (TCV), and total surface area. For ICV, we used the FreeSurfer output measure IntraCranialVol, which is an estimate of ICV based on the Talairach transform. We computed our own measure of TCV by summing the volumes of the TotCort_GM_Vol and Tot_WM_Vol FreeSurfer output measures, which we considered relevant to our study’s focus on cortical functional organization (thus excluding the volumes of subcortical structures). We computed total surface area by using the FreeSurfer mri_surf2surf tool to resample cortical gray matter surface for each subject.

Geodesic distance of connectivity profiles

Geodesic distances, representing the shortest distance between two vertices along the folded cortical mantle’s curvature, were computed using the Micapipe toolbox79, and following the previously described protocol80. In short, geodesic distance matrices were computed for each subject along their native cortical midsurface. The first step consisted in defining a centroid vertex for each cortical parcel, identified as the vertex having the shortest summed Euclidean distance from all other vertices within the parcel. Then, Dijkstra’s algorithm81 was used to compute geodesic distances between the centroid vertices and all other vertices on the on the native midusrface mesh. The vertex-wise geodesic distance values were then averaged within each parcel to form the geodesic distance matrices. From these individual matrices, we finally averaged—parcel-wise—the geodesic distance values of each seed cortical region’s top 10% maximally functionally connected regions per hemisphere, thus obtaining for each subject the mean geodesic distance of functional connectivity profiles by region.

Statistical analysis

Given that the HCP sample includes different levels of kinship, we used linear mixed effects models (LMMs) to account for sibling status (MZ, DZ, non-twin siblings) and family relatedness. In fact, all LMMs mentioned in this work consistently included sex, age, and total surface area as covariates (unless otherwise mentioned), and controlled for random nested effects of family relatedness and sibling status. In addition, effects on cortical data obtained via LMMs underwent false discovery rate (FDR) correction (q < 0.05), thus correcting for multiple comparisons across the 400 Schaefer cortical regions. Throughout this work, we also tested for associations in brain-wide patterns displayed in the form of cortical maps, for which we used Spearman-rank correlation followed by spin-permutation tests to control for spatial autocorrelation82. All statistical tests were two-sided.

After computing the S-A axis of functional cortical organization, we tested for sex differences in the S-A axis loadings with an LMM. Then, we investigated which measure of brain size (out of ICV, TCV, and total surface area) had the largest effects on the S-A axis loadings using separate LMMs (respectively only including ICV, TCV, or total surface area as a covariate, in addition to sex, age and the random nested effect of family relatedness and sibling status). The reason underlying our decision to systematically include total surface area as a covariate in all our LMMs (as the measure of brain size) is that it showed the most widespread effects on the S-A axis loadings out of the three tested measures. Then, we investigated associations between the S-A axis and cortical morphometry, namely the MPC axis and the mean geodesic distance of connectivity profiles, using both LMMs and Spearman-rank spatial correlations of cortical maps.

To probe whether sex differences in cortical morphometry may be associated with sex differences in the S-A axis, we tested whether sex differences in the S-A axis loadings were moderated by total surface area by modeling an additional interaction term of sex by total surface area on the S-A axis loadings within the original LMM. We also tested for sex differences in the MPC axis and in the mean geodesic distance of connectivity profiles, and conducted Spearman-rank correlations of cortical β-maps for the sex contrast in the S-A axis and in the morphometric measures. We also modeled two additional interaction terms within the original LMMs of sex by MPC axis loadings and sex by mean geodesic distance to show their effects on the S-A axis loadings. Finally, we conducted sensitivity analyses to test for sex effects on the S-A axis yielded by an LMM including all morphometric measures as covariates (i.e., including the MPC axis and the mean geodesic distance of connectivity profiles, in addition to total surface area), as well as an LMM not including any morphometric measures as covariates (i.e., also excluding total surface area). We then tested the similarity of both these sex effects with the original sex effects on the S-A axis with a Spearman-rank spatial correlation of the cortical β-maps.

In order to probe the potential intrinsic functional underpinnings of sex differences in the S-A axis, we tested for sex differences in FC strength (also with an LMM), as well as sex differences in FC profiles, i.e., the presence of sex differences in the top 10% maximally functionally connected regions used to compute the S-A axis. To this end, we built 400 × 400 binary matrices at the subject level—based on the subjects’ individual FC matrix z values—in which we marked in a seed-wise fashion (along the matrix rows) whether the given cortical region (along the matrix column) belongs to the given seed’s top 10% maximally functionally connected regions, where 1 indicated that the parcel belongs to the seed’s top 10% maximally functionally connected regions and 0 indicated that the parcel does not belong to the seed’s top 10% maximally functionally connected regions. We then summed the binary matrices separately within sexes in order to fill 160,000 contingency matrices—one for each cell (i.e., functional connection) of the 400 × 400 FC matrix—denoting the number of males and females for which a given cortical region belongs or does not belong to the seed region’s top 10% maximally functionally connected regions (see Table 1 for a visual representation of the contingency matrix structure).Table 1 Contingency matrix structure for the computation of sex differences in functional connectivity profiles

	Cortical region belongs to the seed region’s top 10% maximally functionally connected regions	Cortical region does not belong to the seed region’s top 10% maximally functionally connected regions	
Males	Cm	NCm	
Females	Cf	NCf	
Cm and Cf respectively denote the number of males and females for which the given cortical region (corresponding to the matrix column) constitutes the given seed region’s (corresponding to the matrix row) top 10% maximally functionally connected regions; NCm and NCf, respectively denote the number of males and females for which the given cortical region does not constitute the given seed region’s top 10% maximally functionally connected regions.

We then conducted the Chi-square (χ2) test of independence (degrees of freedom = 1) on each contingency table to test for sex differences in the odds of each parcel of belonging to the top 10% maximally functionally connected regions of each seed region. Given the large number of tests conducted here (400 × 400 = 160,000), we controlled for multiple comparisons using FDR correction. We quantified the size of these sex effects with the odds ratio (OR), calculated as defined in Eq. 2.2 OR=Cm/NCmCf/NCf

where OR > 1 indicates greater male odds—and OR < 1 indicates greater female odds—of a given region of belonging to a given seed’s top 10% maximally functionally connected regions.

We also tested for sex differences in network topology, i.e., how nodes are physically organized in networks and how networks are physically organized along the S-A axis. For this, we computed two measures of network dispersion: between- and within-network dispersion. Between-network dispersion is defined as the Euclidean distance between a pair of network centroids, where a higher value indicates that networks are more segregated from one another along the S-A axis. Within-network dispersion is defined as the sum squared Euclidean distance of network nodes (i.e., S-A axis regional loadings) to the network centroid, where a higher value indicates wider distribution and segregation of a given network’s nodes along the S-A axis. At the individual level, we thus computed between-network dispersion between all networks in a pairwise fashion (21 pairs), and within-network dispersion for all 7 networks, by defining network centroids as the median of the S-A axis loadings of all parcels belonging to a given network, following a previously described method43. Then, we computed sex differences in each of the 21 between-network dispersion metrics and 7 within-network dispersion metrics using LMMs. For each model, we computed a null distribution of β coefficients for sex differences using 1000 spherical rotations of the Schaefer parcellation scheme in order to shuffle the network labels82, against which we computed our p-value to determine statistical significance. We then assed pspinvalues against Bonferroni-corrected two-tailed α-levels of 0.001 (0.025/21) and 0.004 (0.025/7) for between-network and within-network dispersion sex contrasts, respectively.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Peer Review File

Reporting Summary

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-51942-1.

Acknowledgements

We want to thank the Human Connectome Project, Washington University, the University of Minnesota, and Oxford University Consortium (Principal Investigators: David Van Essen and Kamil Ugurbil; 1U54MH091657) originally funded by the 16 N.I.H. Institutes and Centers that support the N.I.H. Blueprint for Neuroscience Research; and by the McDonnell Center for Systems Neuroscience at Washington University. B.S., M.D.H., and G.B. were funded by the German Federal Ministry of Education and Research (BMBF) and the Max Planck Society. J.S. was funded by the Max Planck Society and University of Leipzig. L.W., S.W., and S.B.E. were funded by the European Union’s Horizon 2020 Research and Innovation Program (grant agreements 945539 [HBP SGA3], 826421 [VBC], and 101058516), the DFG (SFB 1451 and IRTG 2150), and the National Institute of Health (NIH; R01 MH074457). S.L.V. was supported by the Max Planck Society through the Otto Hahn Award.

Author contributions

Conceptualization: B.S. and S.L.V. Main analysis and visualization: B.S. Input on analysis: M.D.H., G.B., and S.L.V. Writing—original draft: B.S. Writing—review and editing: B.S., M.D.H., L.W., G.B., J.S., S.W., S.B.E., S.L.V. Supervision: S.L.V.

Peer review

Peer review information

Nature Communications thanks Elvisha Dhamala and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Funding

Open Access funding enabled and organized by Projekt DEAL.

Data availability

All data needed to evaluate the conclusions of the paper are present in the paper and in the Supplementary Materials, and are further available upon request. We obtained all data from the open-access Human Connectome Project (HCP) S1200 young adult sample. The HCP processed data are publicly available and can be directly downloaded at https://db.humanconnectome.org. Source data are provided with this paper.

Code availability

Analyses were conducted in Python and R: The code used in this manuscript is available at https://github.com/biancaserio/sex_diff_gradients (v1; https://zenodo.org/doi/10.5281/zenodo.12785462). The code and tutorials for functional gradient decomposition and to generate geodesic distances can further be found at https://brainspace.readthedocs.io/en/latest/index.html and https://micapipe.readthedocs.io/en/latest/ respectively.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Ankney CD The brain size/IQ debate Nature 1992 360 292 292 10.1038/360292b0 1448139
Ankney, C. D. The brain size/IQ debate. Nature 360, 292–292 (1992).1448139 10.1038/360292b0
2. Eliot L Dump the “dimorphism”: comprehensive synthesis of human brain studies reveals few male-female differences beyond size Neurosci. Biobehav. Rev. 2021 125 667 697 10.1016/j.neubiorev.2021.02.026 33621637
Eliot, L. et al. Dump the “dimorphism”: comprehensive synthesis of human brain studies reveals few male-female differences beyond size. Neurosci. Biobehav. Rev. 125, 667–697 (2021).33621637 10.1016/j.neubiorev.2021.02.026
3. Leonard CM Size matters: cerebral volume influences sex differences in neuroanatomy Cereb. Cortex 2008 18 2920 2931 10.1093/cercor/bhn052 18440950
Leonard, C. M. et al. Size matters: cerebral volume influences sex differences in neuroanatomy. Cereb. Cortex 18, 2920–2931 (2008).18440950 10.1093/cercor/bhn052
4. Peters M Sex differences in human brain size and the general meaning of differences in brain size Can. J. Psychol./Rev. 1991 45 507 10.1037/h0084307
Peters, M. Sex differences in human brain size and the general meaning of differences in brain size. Can. J. Psychol./Rev. 45, 507 (1991).10.1037/h0084307
5. Ritchie SJ Sex differences in the adult human brain: evidence from 5216 UK biobank participants Cereb. Cortex 2018 28 2959 2975 10.1093/cercor/bhy109 29771288
Ritchie, S. J. et al. Sex differences in the adult human brain: evidence from 5216 UK biobank participants. Cereb. Cortex 28, 2959–2975 (2018).29771288 10.1093/cercor/bhy109
6. Ruigrok AN A meta-analysis of sex differences in human brain structure Neurosci. Biobehav. Rev. 2014 39 34 50 10.1016/j.neubiorev.2013.12.004 24374381
Ruigrok, A. N. et al. A meta-analysis of sex differences in human brain structure. Neurosci. Biobehav. Rev. 39, 34–50 (2014).24374381 10.1016/j.neubiorev.2013.12.004
7. Williams CM Neuroanatomical norms in the UK Biobank: the impact of allometric scaling, sex, and age Hum. Brain Mapp. 2021 42 4623 4642 10.1002/hbm.25572 34268815
Williams, C. M. et al. Neuroanatomical norms in the UK Biobank: the impact of allometric scaling, sex, and age. Hum. Brain Mapp. 42, 4623–4642 (2021).34268815 10.1002/hbm.25572
8. Shanmugan S Sex differences in the functional topography of association networks in youth Proc. Natl Acad. Sci. 2022 119 e2110416119 10.1073/pnas.2110416119 35939696
Shanmugan, S. et al. Sex differences in the functional topography of association networks in youth. Proc. Natl Acad. Sci. 119, e2110416119 (2022).35939696 10.1073/pnas.2110416119
9. Weis S Sex classification by resting state brain connectivity Cereb. Cortex 2020 30 824 835 10.1093/cercor/bhz129 31251328
Weis, S. et al. Sex classification by resting state brain connectivity. Cereb. Cortex 30, 824–835 (2020).31251328 10.1093/cercor/bhz129
10. Margulies DS Situating the default-mode network along a principal gradient of macroscale cortical organization Proc. Natl Acad. Sci. 2016 113 12574 12579 10.1073/pnas.1608282113 27791099
Margulies, D. S. et al. Situating the default-mode network along a principal gradient of macroscale cortical organization. Proc. Natl Acad. Sci. 113, 12574–12579 (2016).27791099 10.1073/pnas.1608282113
11. Mesulam M-M From sensation to cognition Brain J. Neurol. 1998 121 1013 1052 10.1093/brain/121.6.1013
Mesulam, M.-M. From sensation to cognition. Brain J. Neurol. 121, 1013–1052 (1998).10.1093/brain/121.6.1013
12. Will M Different environmental variables predict body and brain size evolution in Homo Nat. Commun. 2021 12 4116 10.1038/s41467-021-24290-7 34238930
Will, M. et al. Different environmental variables predict body and brain size evolution in Homo. Nat. Commun. 12, 4116 (2021).34238930 10.1038/s41467-021-24290-7
13. Buckner RL Krienen FM The evolution of distributed association networks in the human brain Trends Cogn. Sci. 2013 17 648 665 10.1016/j.tics.2013.09.017 24210963
Buckner, R. L. & Krienen, F. M. The evolution of distributed association networks in the human brain. Trends Cogn. Sci. 17, 648–665 (2013).24210963 10.1016/j.tics.2013.09.017
14. Reardon P Normative brain size variation and brain shape diversity in humans Science 2018 360 1222 1227 10.1126/science.aar2578 29853553
Reardon, P. et al. Normative brain size variation and brain shape diversity in humans. Science 360, 1222–1227 (2018).29853553 10.1126/science.aar2578
15. Sydnor VJ Neurodevelopment of the association cortices: patterns, mechanisms, and implications for psychopathology Neuron 2021 109 2820 2846 10.1016/j.neuron.2021.06.016 34270921
Sydnor, V. J. et al. Neurodevelopment of the association cortices: patterns, mechanisms, and implications for psychopathology. Neuron 109, 2820–2846 (2021).34270921 10.1016/j.neuron.2021.06.016
16. Luders E Toga AW Sex differences in brain anatomy Prog. Brain Res. 2010 186 2 12 10.1016/B978-0-444-53630-3.00001-4
Luders, E. & Toga, A. W. Sex differences in brain anatomy. Prog. Brain Res. 186, 2–12 (2010).10.1016/B978-0-444-53630-3.00001-4
17. Wierenga LM Sex effects on development of brain structure and executive functions: greater variance than mean effects J. Cogn. Neurosci. 2019 31 730 753 10.1162/jocn_a_01375 30726177
Wierenga, L. M. et al. Sex effects on development of brain structure and executive functions: greater variance than mean effects. J. Cogn. Neurosci. 31, 730–753 (2019).30726177 10.1162/jocn_a_01375
18. Eliot L Remembering the null hypothesis when searching for brain sex differences Biol. sex. Differ. 2024 15 1 7 10.1186/s13293-024-00585-4 38178264
Eliot, L. Remembering the null hypothesis when searching for brain sex differences. Biol. sex. Differ. 15, 1–7 (2024).38178264 10.1186/s13293-024-00585-4
19. Liu S Integrative structural, functional, and transcriptomic analyses of sex-biased brain organization in humans Proc. Natl Acad. Sci. 2020 117 18788 18798 10.1073/pnas.1919091117 32690678
Liu, S. et al. Integrative structural, functional, and transcriptomic analyses of sex-biased brain organization in humans. Proc. Natl Acad. Sci. 117, 18788–18798 (2020).32690678 10.1073/pnas.1919091117
20. Raznahan A Patterns of coordinated anatomical change in human cortical development: a longitudinal neuroimaging study of maturational coupling Neuron 2011 72 873 884 10.1016/j.neuron.2011.09.028 22153381
Raznahan, A. et al. Patterns of coordinated anatomical change in human cortical development: a longitudinal neuroimaging study of maturational coupling. Neuron 72, 873–884 (2011).22153381 10.1016/j.neuron.2011.09.028
21. Allen EA A baseline for the multivariate comparison of resting-state networks Front. Syst. Neurosci. 2011 5 2 10.3389/fnsys.2011.00002 21442040
Allen, E. A. et al. A baseline for the multivariate comparison of resting-state networks. Front. Syst. Neurosci. 5, 2 (2011).21442040 10.3389/fnsys.2011.00002
22. Biswal BB Toward discovery science of human brain function Proc. Natl Acad. Sci. 2010 107 4734 4739 10.1073/pnas.0911855107 20176931
Biswal, B. B. et al. Toward discovery science of human brain function. Proc. Natl Acad. Sci. 107, 4734–4739 (2010).20176931 10.1073/pnas.0911855107
23. Bluhm RL Default mode network connectivity: effects of age, sex, and analytic approach Neuroreport 2008 19 887 891 10.1097/WNR.0b013e328300ebbf 18463507
Bluhm, R. L. et al. Default mode network connectivity: effects of age, sex, and analytic approach. Neuroreport 19, 887–891 (2008).18463507 10.1097/WNR.0b013e328300ebbf
24. Zhang X Gender differences are encoded differently in the structure and function of the human brain revealed by multimodal MRI Front. Hum. Neurosci. 2020 14 244 10.3389/fnhum.2020.00244 32792927
Zhang, X. et al. Gender differences are encoded differently in the structure and function of the human brain revealed by multimodal MRI. Front. Hum. Neurosci. 14, 244 (2020).32792927 10.3389/fnhum.2020.00244
25. Scheinost D Sex differences in normal age trajectories of functional brain networks Hum. Brain Mapp. 2015 36 1524 1535 10.1002/hbm.22720 25523617
Scheinost, D. et al. Sex differences in normal age trajectories of functional brain networks. Hum. Brain Mapp. 36, 1524–1535 (2015).25523617 10.1002/hbm.22720
26. Wang Y Long-range functional connections mirror and link microarchitectural and cognitive hierarchies in the human brain Cereb. Cortex 2023 33 1782 1798 10.1093/cercor/bhac172 35596951
Wang, Y. et al. Long-range functional connections mirror and link microarchitectural and cognitive hierarchies in the human brain. Cereb. Cortex 33, 1782–1798 (2023).35596951 10.1093/cercor/bhac172
27. Leech R Variation in spatial dependencies across the cortical mantle discriminates the functional behaviour of primary and association cortex Nat. Commun. 2023 14 5656 10.1038/s41467-023-41334-2 37704600
Leech, R. et al. Variation in spatial dependencies across the cortical mantle discriminates the functional behaviour of primary and association cortex. Nat. Commun. 14, 5656 (2023).37704600 10.1038/s41467-023-41334-2
28. Markov NT Cortical high-density counterstream architectures Science 2013 342 1238406 10.1126/science.1238406 24179228
Markov, N. T. et al. Cortical high-density counterstream architectures. Science 342, 1238406 (2013).24179228 10.1126/science.1238406
29. Pang, J. C., et al., Geometric constraints on human brain function. Nature 618, 566–574 (2023).
30. Robinson P Modal analysis of corticothalamic dynamics, electroencephalographic spectra, and evoked potentials Phys. Rev. E 2001 63 041909 10.1103/PhysRevE.63.041909
Robinson, P. et al. Modal analysis of corticothalamic dynamics, electroencephalographic spectra, and evoked potentials. Phys. Rev. E 63, 041909 (2001).10.1103/PhysRevE.63.041909
31. Wingeier BM Nunez PL Silberstein RB Spherical harmonic decomposition applied to spatial-temporal analysis of human high-density electroencephalogram Phys. Rev. E 2001 64 051916 10.1103/PhysRevE.64.051916
Wingeier, B. M., Nunez, P. L. & Silberstein, R. B. Spherical harmonic decomposition applied to spatial-temporal analysis of human high-density electroencephalogram. Phys. Rev. E 64, 051916 (2001).10.1103/PhysRevE.64.051916
32. Zhang C Functional connectivity predicts gender: evidence for gender differences in resting brain connectivity Hum. Brain Mapp. 2018 39 1765 1776 10.1002/hbm.23950 29322586
Zhang, C. et al. Functional connectivity predicts gender: evidence for gender differences in resting brain connectivity. Hum. Brain Mapp. 39, 1765–1776 (2018).29322586 10.1002/hbm.23950
33. Van Essen DC The WU-Minn human connectome project: an overview Neuroimage 2013 80 62 79 10.1016/j.neuroimage.2013.05.041 23684880
Van Essen, D. C. et al. The WU-Minn human connectome project: an overview. Neuroimage 80, 62–79 (2013).23684880 10.1016/j.neuroimage.2013.05.041
34. Coifman RR Lafon S Diffusion maps Appl. Comput. Harmonic Anal. 2006 21 5 30 10.1016/j.acha.2006.04.006
Coifman, R. R. & Lafon, S. Diffusion maps. Appl. Comput. Harmonic Anal. 21, 5–30 (2006).10.1016/j.acha.2006.04.006
35. Vos de Wael R BrainSpace: a toolbox for the analysis of macroscale gradients in neuroimaging and connectomics datasets Commun. Biol. 2020 3 103 10.1038/s42003-020-0794-7 32139786
Vos de Wael, R. et al. BrainSpace: a toolbox for the analysis of macroscale gradients in neuroimaging and connectomics datasets. Commun. Biol. 3, 103 (2020).32139786 10.1038/s42003-020-0794-7
36. Hong S-J Atypical functional connectome hierarchy in autism Nat. Commun. 2019 10 1022 10.1038/s41467-019-08944-1 30833582
Hong, S.-J. et al. Atypical functional connectome hierarchy in autism. Nat. Commun. 10, 1022 (2019).30833582 10.1038/s41467-019-08944-1
37. Valk SL Genetic and phylogenetic uncoupling of structure and function in human transmodal cortex Nat. Commun. 2022 13 2341 10.1038/s41467-022-29886-1 35534454
Valk, S. L. et al. Genetic and phylogenetic uncoupling of structure and function in human transmodal cortex. Nat. Commun. 13, 2341 (2022).35534454 10.1038/s41467-022-29886-1
38. Hong S-J Toward a connectivity gradient-based framework for reproducible biomarker discovery NeuroImage 2020 223 117322 10.1016/j.neuroimage.2020.117322 32882388
Hong, S.-J. et al. Toward a connectivity gradient-based framework for reproducible biomarker discovery. NeuroImage 223, 117322 (2020).32882388 10.1016/j.neuroimage.2020.117322
39. Knodt AR Test–retest reliability and predictive utility of a macroscale principal functional connectivity gradient Hum. Brain Mapp. 2023 44 6399 6417 10.1002/hbm.26517 37851700
Knodt, A. R. et al. Test–retest reliability and predictive utility of a macroscale principal functional connectivity gradient. Hum. Brain Mapp. 44, 6399–6417 (2023).37851700 10.1002/hbm.26517
40. Yeo, B. T. et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J. Neurophysiol. (2011).
41. Paquola C Shifts in myeloarchitecture characterise adolescent development of cortical gradients elife 2019 8 e50482 10.7554/eLife.50482 31724948
Paquola, C. et al. Shifts in myeloarchitecture characterise adolescent development of cortical gradients. elife 8, e50482 (2019).31724948 10.7554/eLife.50482
42. Paquola C Microstructural and functional gradients are increasingly dissociated in transmodal cortices PLoS Biol. 2019 17 e3000284 10.1371/journal.pbio.3000284 31107870
Paquola, C. et al. Microstructural and functional gradients are increasingly dissociated in transmodal cortices. PLoS Biol. 17, e3000284 (2019).31107870 10.1371/journal.pbio.3000284
43. Bethlehem RA Dispersion of functional gradients across the adult lifespan Neuroimage 2020 222 117299 10.1016/j.neuroimage.2020.117299 32828920
Bethlehem, R. A. et al. Dispersion of functional gradients across the adult lifespan. Neuroimage 222, 117299 (2020).32828920 10.1016/j.neuroimage.2020.117299
44. Dhamala E Proportional intracranial volume correction differentially biases behavioral predictions across neuroanatomical features, sexes, and development NeuroImage 2022 260 119485 10.1016/j.neuroimage.2022.119485 35843514
Dhamala, E. et al. Proportional intracranial volume correction differentially biases behavioral predictions across neuroanatomical features, sexes, and development. NeuroImage 260, 119485 (2022).35843514 10.1016/j.neuroimage.2022.119485
45. More, S., et al. Confound removal and normalization in practice: a neuroimaging based sex prediction case study. In Proc. Machine Learning and Knowledge Discovery in Databases. Applied Data Science and Demo Track: European Conference, ECML PKDD 2020, Ghent, Belgium, September 14–18, 2020, Proceedings, Part V (Springer, 2021).
46. Pintzka CW Marked effects of intracranial volume correction methods on sex differences in neuroanatomical structures: a HUNT MRI study Front. Neurosci. 2015 9 238 10.3389/fnins.2015.00238 26217172
Pintzka, C. W. et al. Marked effects of intracranial volume correction methods on sex differences in neuroanatomical structures: a HUNT MRI study. Front. Neurosci. 9, 238 (2015).26217172 10.3389/fnins.2015.00238
47. Sanchis-Segura C Sex differences in gray matter volume: how many and how large are they really? Biol. Sex. Differ. 2019 10 1 19 10.1186/s13293-019-0245-7 30606254
Sanchis-Segura, C. et al. Sex differences in gray matter volume: how many and how large are they really? Biol. Sex. Differ. 10, 1–19 (2019).30606254 10.1186/s13293-019-0245-7
48. Bignardi, G. et al. Pervasive inter-individual differences in the sensorimotor-association axis of cortical organization. bioRxiv, p. 2023.07. 13.548817 (2023).
49. Park H-J Friston K Structural and functional brain networks: from connections to cognition Science 2013 342 1238411 10.1126/science.1238411 24179229
Park, H.-J. & Friston, K. Structural and functional brain networks: from connections to cognition. Science 342, 1238411 (2013).24179229 10.1126/science.1238411
50. Glasser MF Van Essen DC Mapping human cortical areas in vivo based on myelin content as revealed by T1-and T2-weighted MRI J. Neurosci. 2011 31 11597 11616 10.1523/JNEUROSCI.2180-11.2011 21832190
Glasser, M. F. & Van Essen, D. C. Mapping human cortical areas in vivo based on myelin content as revealed by T1-and T2-weighted MRI. J. Neurosci. 31, 11597–11616 (2011).21832190 10.1523/JNEUROSCI.2180-11.2011
51. Sandrone S Mapping myelin in white matter with T1-weighted/T2-weighted maps: discrepancy with histology and other myelin MRI measures Brain Struct. Funct. 2023 228 525 535 10.1007/s00429-022-02600-z 36692695
Sandrone, S. et al. Mapping myelin in white matter with T1-weighted/T2-weighted maps: discrepancy with histology and other myelin MRI measures. Brain Struct. Funct. 228, 525–535 (2023).36692695 10.1007/s00429-022-02600-z
52. Uddin MN Comparisons between multi-component myelin water fraction, T1w/T2w ratio, and diffusion tensor imaging measures in healthy human brain structures Sci. Rep. 2019 9 2500 10.1038/s41598-019-39199-x 30792440
Uddin, M. N. et al. Comparisons between multi-component myelin water fraction, T1w/T2w ratio, and diffusion tensor imaging measures in healthy human brain structures. Sci. Rep. 9, 2500 (2019).30792440 10.1038/s41598-019-39199-x
53. Foo H Age-and sex-related topological organization of human brain functional networks and their relationship to cognition Front. Aging Neurosci. 2021 13 897 10.3389/fnagi.2021.758817
Foo, H. et al. Age-and sex-related topological organization of human brain functional networks and their relationship to cognition. Front. Aging Neurosci. 13, 897 (2021).10.3389/fnagi.2021.758817
54. Shine JM Poldrack RA Principles of dynamic network reconfiguration across diverse brain states NeuroImage 2018 180 396 405 10.1016/j.neuroimage.2017.08.010 28782684
Shine, J. M. & Poldrack, R. A. Principles of dynamic network reconfiguration across diverse brain states. NeuroImage 180, 396–405 (2018).28782684 10.1016/j.neuroimage.2017.08.010
55. Chang C Tracking brain arousal fluctuations with fMRI Proc. Natl Acad. Sci. 2016 113 4518 4523 10.1073/pnas.1520613113 27051064
Chang, C. et al. Tracking brain arousal fluctuations with fMRI. Proc. Natl Acad. Sci. 113, 4518–4523 (2016).27051064 10.1073/pnas.1520613113
56. Barttfeld P Signature of consciousness in the dynamics of resting-state brain activity Proc. Natl Acad. Sci. 2015 112 887 892 10.1073/pnas.1418031112 25561541
Barttfeld, P. et al. Signature of consciousness in the dynamics of resting-state brain activity. Proc. Natl Acad. Sci. 112, 887–892 (2015).25561541 10.1073/pnas.1418031112
57. Cole MW Multi-task connectivity reveals flexible hubs for adaptive task control Nat. Neurosci. 2013 16 1348 1355 10.1038/nn.3470 23892552
Cole, M. W. et al. Multi-task connectivity reveals flexible hubs for adaptive task control. Nat. Neurosci. 16, 1348–1355 (2013).23892552 10.1038/nn.3470
58. Douw L State-dependent variability of dynamic functional connectivity between frontoparietal and default networks relates to cognitive flexibility Neuroscience 2016 339 12 21 10.1016/j.neuroscience.2016.09.034 27687802
Douw, L. et al. State-dependent variability of dynamic functional connectivity between frontoparietal and default networks relates to cognitive flexibility. Neuroscience 339, 12–21 (2016).27687802 10.1016/j.neuroscience.2016.09.034
59. Tognoli E Kelso JS The metastable brain Neuron 2014 81 35 48 10.1016/j.neuron.2013.12.022 24411730
Tognoli, E. & Kelso, J. S. The metastable brain. Neuron 81, 35–48 (2014).24411730 10.1016/j.neuron.2013.12.022
60. Castrillon, G. et al. An energy costly architecture of neuromodulators for human brain evolution and cognition. Sci. Adv., 9, eadi7632 (2023).
61. Wiersch L Weis S Sex differences in the brain: More than just male or female Cogn. Neurosci. 2021 12 187 188 10.1080/17588928.2020.1867084 33406985
Wiersch, L. & Weis, S. Sex differences in the brain: More than just male or female. Cogn. Neurosci. 12, 187–188 (2021).33406985 10.1080/17588928.2020.1867084
62. Hyde JS The future of sex and gender in psychology: five challenges to the gender binary Am. Psychol. 2019 74 171 10.1037/amp0000307 30024214
Hyde, J. S. et al. The future of sex and gender in psychology: five challenges to the gender binary. Am. Psychol. 74, 171 (2019).30024214 10.1037/amp0000307
63. Schulte, J. et al. The global communication architecture of the human brain transcends the subcortical-cortical-cerebellar subdivisions. bioRxiv, 2023.07. 07.548139 (2023).
64. Stevens JS Hamann S Sex differences in brain activation to emotional stimuli: a meta-analysis of neuroimaging studies Neuropsychologia 2012 50 1578 1593 10.1016/j.neuropsychologia.2012.03.011 22450197
Stevens, J. S. & Hamann, S. Sex differences in brain activation to emotional stimuli: a meta-analysis of neuroimaging studies. Neuropsychologia 50, 1578–1593 (2012).22450197 10.1016/j.neuropsychologia.2012.03.011
65. Collaborators, G.M.D., Global, regional, and national burden of 12 mental disorders in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet Psychiatry 9:137–150 (2022).
66. Klein S Increased neural reactivity to emotional pictures in men with high hair testosterone concentrations Soc. Cogn. Affect. Neurosci. 2019 14 1009 1016 10.1093/scan/nsz067 31506704
Klein, S. et al. Increased neural reactivity to emotional pictures in men with high hair testosterone concentrations. Soc. Cogn. Affect. Neurosci. 14, 1009–1016 (2019).31506704 10.1093/scan/nsz067
67. Pritschet L Functional reorganization of brain networks across the human menstrual cycle Neuroimage 2020 220 117091 10.1016/j.neuroimage.2020.117091 32621974
Pritschet, L. et al. Functional reorganization of brain networks across the human menstrual cycle. Neuroimage 220, 117091 (2020).32621974 10.1016/j.neuroimage.2020.117091
68. Shansky RM Murphy AZ Considering sex as a biological variable will require a global shift in science culture Nat. Neurosci. 2021 24 457 464 10.1038/s41593-021-00806-8 33649507
Shansky, R. M. & Murphy, A. Z. Considering sex as a biological variable will require a global shift in science culture. Nat. Neurosci. 24, 457–464 (2021).33649507 10.1038/s41593-021-00806-8
69. Ratnu VS Emami MR Bredy TW Genetic and epigenetic factors underlying sex differences in the regulation of gene expression in the brain J. Neurosci. Res. 2017 95 301 310 10.1002/jnr.23886 27870402
Ratnu, V. S., Emami, M. R. & Bredy, T. W. Genetic and epigenetic factors underlying sex differences in the regulation of gene expression in the brain. J. Neurosci. Res. 95, 301–310 (2017).27870402 10.1002/jnr.23886
70. Kundakovic, M. Tickerhoof, M. M. Epigenetic mechanisms underlying sex differences in the brain and behavior. Trends Neurosci. 47, 19–35 (2023).
71. Slavich GM Sacher J Stress, sex hormones, inflammation, and major depressive disorder: extending social signal transduction theory of depression to account for sex differences in mood disorders Psychopharmacology 2019 236 3063 3079 10.1007/s00213-019-05326-9 31359117
Slavich, G. M. & Sacher, J. Stress, sex hormones, inflammation, and major depressive disorder: extending social signal transduction theory of depression to account for sex differences in mood disorders. Psychopharmacology 236, 3063–3079 (2019).31359117 10.1007/s00213-019-05326-9
72. Van Essen DC The Human Connectome Project: a data acquisition perspective Neuroimage 2012 62 2222 2231 10.1016/j.neuroimage.2012.02.018 22366334
Van Essen, D. C. et al. The Human Connectome Project: a data acquisition perspective. Neuroimage 62, 2222–2231 (2012).22366334 10.1016/j.neuroimage.2012.02.018
73. Dhamala E Shared functional connections within and between cortical networks predict cognitive abilities in adult males and females Hum. Brain Mapp. 2022 43 1087 1102 10.1002/hbm.25709 34811849
Dhamala, E. et al. Shared functional connections within and between cortical networks predict cognitive abilities in adult males and females. Hum. Brain Mapp. 43, 1087–1102 (2022).34811849 10.1002/hbm.25709
74. Glasser MF The minimal preprocessing pipelines for the Human Connectome Project Neuroimage 2013 80 105 124 10.1016/j.neuroimage.2013.04.127 23668970
Glasser, M. F. et al. The minimal preprocessing pipelines for the Human Connectome Project. Neuroimage 80, 105–124 (2013).23668970 10.1016/j.neuroimage.2013.04.127
75. Smith SM Resting-state fMRI in the human connectome project Neuroimage 2013 80 144 168 10.1016/j.neuroimage.2013.05.039 23702415
Smith, S. M. et al. Resting-state fMRI in the human connectome project. Neuroimage 80, 144–168 (2013).23702415 10.1016/j.neuroimage.2013.05.039
76. Schaefer A Local-global parcellation of the human cerebral cortex from intrinsic functional connectivity MRI Cereb. Cortex 2018 28 3095 3114 10.1093/cercor/bhx179 28981612
Schaefer, A. et al. Local-global parcellation of the human cerebral cortex from intrinsic functional connectivity MRI. Cereb. Cortex 28, 3095–3114 (2018).28981612 10.1093/cercor/bhx179
77. Schleicher A Observer-independent method for microstructural parcellation of cerebral cortex: a quantitative approach to cytoarchitectonics Neuroimage 1999 9 165 177 10.1006/nimg.1998.0385 9918738
Schleicher, A. et al. Observer-independent method for microstructural parcellation of cerebral cortex: a quantitative approach to cytoarchitectonics. Neuroimage 9, 165–177 (1999).9918738 10.1006/nimg.1998.0385
78. Zilles K Architectonics of the human cerebral cortex and transmitter receptor fingerprints: reconciling functional neuroanatomy and neurochemistry Eur. Neuropsychopharmacol. 2002 12 587 599 10.1016/S0924-977X(02)00108-6 12468022
Zilles, K. et al. Architectonics of the human cerebral cortex and transmitter receptor fingerprints: reconciling functional neuroanatomy and neurochemistry. Eur. Neuropsychopharmacol. 12, 587–599 (2002).12468022 10.1016/S0924-977X(02)00108-6
79. Cruces RR Micapipe: a pipeline for multimodal neuroimaging and connectome analysis Neuroimage 2022 263 119612 10.1016/j.neuroimage.2022.119612 36070839
Cruces, R. R. et al. Micapipe: a pipeline for multimodal neuroimaging and connectome analysis. Neuroimage 263, 119612 (2022).36070839 10.1016/j.neuroimage.2022.119612
80. Royer J An open MRI dataset for multiscale neuroscience Sci. Data 2022 9 569 10.1038/s41597-022-01682-y 36109562
Royer, J. et al. An open MRI dataset for multiscale neuroscience. Sci. Data 9, 569 (2022).36109562 10.1038/s41597-022-01682-y
81. Dijkstra, E. W., A note on two problems in connexion with graphs, In Edsger Wybe Dijkstra: His Life, Work, and Legacy. 287–290 (Association for Computing Machinery and Morgan & Claypool Publishers, 2022).
82. Alexander-Bloch AF On testing for spatial correspondence between maps of human brain structure and function Neuroimage 2018 178 540 551 10.1016/j.neuroimage.2018.05.070 29860082
Alexander-Bloch, A. F. et al. On testing for spatial correspondence between maps of human brain structure and function. Neuroimage 178, 540–551 (2018).29860082 10.1016/j.neuroimage.2018.05.070
