Sex differences in brain homotopic co-activations: a meta-analytic study

An element of great interest in functional connectivity is ‘homotopic connectivity’ (HC), namely the connectivity between two mirrored areas of the two hemispheres, mainly mediated by the fibers of the corpus callosum. Despite a long tradition of studying sexual dimorphism in the human brain, to our knowledge only one study has addressed the influence of sex on HC. We investigated the issue of homotopic co-activations in women and men using a coordinate-based meta-analytic method and data from the BrainMap database. A first unexpected observation was that the database was affected by a sex bias: women-only groups are investigated less often than men-only ones, and they are more often studied in certain domains such as emotion compared to men, and less in cognition. Implementing a series of sampling procedures to equalize the size and proportion of the datasets, our results indicated that females exhibit stronger interhemispheric co-activation than males, suggesting that the female brain is less lateralized and more integrated than that of males. In addition, males appear to show less intense but more extensive co-activation than females. Some local differences also appeared. In particular, it appears that primary motor and perceptual areas are more co-activated in males, in contrast to the opposite trend in the rest of the brain. This argues for a multidimensional view of sex brain differences and suggests that the issue should be approached with more complex models than previously thought. Supplementary Information The online version contains supplementary material available at 10.1007/s00429-022-02572-0.


Introduction
Two properties were suggested to characterize brain connectivity: segregation and integration (Sporns, 2013). Segregation refers to a network organization characterized by clusters of nodes (neurons or areas) more strongly connected between them than with other clusters, promoting functional specialization (Jager and Postma 2003;Sporns 2013). Integration refers to the interaction between specialized brain regions that allows both long-range synchronization and flow of information (Friston 2002;Heuvel & Sporns 2013). The interaction between these two processes defines the pattern of correlation and anticorrelation between resting state networks (Fox et al. 2005a), as well as the functional lateralization of the two hemispheres (Gotts et al. 2013).
A relevant element of FC is homotopic connectivity (HC). HC refers to the (structural or functional) connectivity between two homologous areas of the two hemispheres, mainly mediated by the fibers of the corpus callosum (Goldstein et al. 2021;Mancuso et al. 2019a;Mollink et al. 2019). Crucially, HC can be taken as an index of local interhemispheric integration (Jin et al. 2020). Values of homotopic FC are higher than those of other forms of intra-and interhemispheric FC (Stark et al. 2008), and it correlates with the performance of cognitive tasks, such as visuospatial attention (Gracia-Tabuenca et al. 2018) and executive function (Zhao et al. 2020;Vallesi et al. 2021). Furthermore, altered patterns of homotopic FC were observed in a wide array of Chiara Bonelli and Lorenzo Mancuso contributed equally to the work. 1 3 pathologic conditions or abnormal states (Shan et al. 2021;Yu et al. 2018;Chen et al. 2016;Guo et al. 2019;Deng et al. 2021;Loomba et al. 2021;Yao et al 2021;Li et al. 2015;Zhao et al. 2017). The presence of functional HC has been also confirmed by numerous meta-analytic studies based on neuroimaging data and by applying different methods (Laird et al. 2011;Toro et al. 2008). Indeed, previous work from our group adopted a meta-analytical framework to evaluate the homotopic co-activation patterns (Mancuso et al. 2019b), confirming that homotopic FC is stronger for primary than associative regions, as observed by others (Stark et al. 2008;Aboitiz and Montiel 2003). However, that study did not address the possible influence of confounding variables such as age and sex.
The study of sex differences in the human brain is of great interest, also and especially to identify mechanisms underlying behavioral differences or to explain the prevalence of some psychiatric and neurological disorders (Salminen et al. 2021). However, data are often ambiguous and their interpretations are contradictory (Wiersch, L. and Weis, S. 2021;McCarthy, M. M. 2016), and different views on the matter are livelily debated (Wiersch and Weis 2021;McCarthy 2016;Joel et al. 2015;Del Giudice 2021;Joel et al. 2020;Grabowska 2017;Zhang et al. 2021;Joel et al. 2018;Joel 2021;Byne et al. 1988;Hirnstein and Hausmann 2021;Williams et al. 2021). More specifically, it has long been discussed if there are any anatomical differences in the corpus callosum between females and males, with most studies suggesting stronger callosal connectivity in females (Schmied et al. 2020;Ardekani et al. 2013;Suganthy et al. 2003;Habib et al. 1991;Luders et al. 2014;Bruner et al. 2012;Kanaan et al. 2012;Shino et al. 2017). Regarding functional HC, Zuo and colleagues tested the local differences between sexes, reporting greater voxel-wise resting-state HC in the posterior cingulate cortex and in the medial and lateral prefrontal cortex in females; on the contrary, males showed greater HC in the cerebellum, parahippocampal gyrus and fusiform gyrus (Zuo et al. 2010).
These findings are suggestive of different mechanisms of segregation and integration in women and men. In fact, results from Ingalhalikar and colleagues (Ingalhalikar et al. 2014) indicate that females exhibit stronger structural interhemispheric connectivity than men, unlike males in which the intra-hemispheric connections prevail compared to the other sex. However, other works suggest that such results do not hold if differences in brain volume are taken into account (Hänggi et al. 2014;Martinez et al. 2017). A related issue is that of sexual differences in hemispheric lateralization. If the characterization by Ingalhaikar and colleagues was true, regardless of any brain volume nuisance, we should expect that women had generally less lateralization and thus stronger interhemispherical co-activations than men. In fact, better microstructural organization of the corpus callosum has been associated with more symmetric processing of dichotic stimuli (Friederich et al. 2017), while regions with weaker callosal connections showed more asymmetric activations (Karolis et al. 2019). These observations suggest that stronger callosal connections may lead to less lateralization and thus more likely homotopic co-activations. Some empirical and theoretical accounts support the hypothesis of a lesser lateralization in female subjects (McGlone 1980;Hiscock et al. 1995), but the literature is inconsistent (Voyer et al. 1995;Boles 2005). For instance, there is some evidence for leftward lateralization from behavioural and FC studies in females, which would explain their greater ability in verbal tasks than males (Hjelmervik 2012;Tomasi and Volkow 2012). On the contrary, males would show lateralization to the right consistently with their better performance in visual-spatial tasks (Voyer et al. 1995;Tomasi and Volkow 2012;Wendt and Risberg 1994). However, other studies provide contradicting evidence (Agcaoglu et al. 2015;Liu et al. 2009;Nielsen et al. 2013).
In the present study, we investigated the matter of homotopic co-activations in women and men, employing the meta-analytical method previously presented by our group (Mancuso et al. 2019b). We assessed the meta-analytical homotopic connectivity (MHC) of female and male subjects on a large database of activation studies obtained from the BrainMap database. There is a large body of evidence showing the similarity between co-activations and functional connectivity (Toro et al. 2008;Smith et al. 2009;Laird et al. 2011;Laird et al. 2013;Torta et al. 2013;Cauda et al. 2021); in particular, the MHC technique has previously been shown to produce similar results to the VMHC (Mancuso et al. 2019b). This would suggest that meta-analytic homotopic co-activations could be taken both as an indicator of interhemispheric integration (vs. lateralization) and of functional connectivity. Furthermore, we chose to focus on singular cognitive domains, to investigate if the patterns of interhemispheric integration of the two sexes varie between them. If it is true that males are more lateralized than females, we should expect to find stronger homotopic co-activations in females. However, given the complex and conflicting results of sex neuroscience, we anticipated that results could be not straightforward, with complex local variations in HC between sexes and across cognitive domains.

Data collection
To obtain the activation data to model the homotopic coactivations, the BrainMap (Fox et al. 2005b;Fox and Lancaster 2002;Laird et al. 2005)  where, as Behavioural Domain, we selected each time one of the five options (i.e., Action, Emotion, Cognition, Interoception, Perception), always including all the available sub-domains for each main domain. Each query was repeated for both sexes.
Subsequently, the list of experiments was carefully screened to verify that the Sleuth output was consistent with the given query (Manuello et al. 2022), that is, that the studies did not include individuals of the opposite sex, or that all the subjects were healthy and were not taking medications. We also decided to exclude those experiments investigating a group of female subjects during their menstrual cycle, as it was indicated that it could influence the fMRI signal (Fernàndez et al. 2003;Hjelmervik et al. 2014;Donishi et al. 2017;Pritschet et al. 2020). Furthermore, hormonal fluctuation during the menstrual cycle also influences functional laterality (Hausmann et al. 2002;Hausmann 2005;Pletzer et al. 2019

Meta-analytic homotopic connectivity calculation
As a first step, the cerebral cortex was parcellated in couple of homotopic regions. To do so, the 800 volumes functional atlas by Schaefer et al. (2018) was symmetrized, substituting the right hemisphere with a flipped image of the left one. The purpose of such a passage was to obtain symmetrical MHC maps similar to their non-meta-analytical counterpart, the VMHC. The atlas was also converted to Talairach space. We worked in Talairach space to take advantage of the Talairach Daemon client (Lancaster et al. 2000;Lancaster et al. 1997) to label the areas.
Meanwhile, the 12 datasets of coordinates (i.e., 2 domaingeneral and 10 domain-specific) resulting from the queries were exported from Sleuth in Talairach standard space (Fig. 1A). As pointed out by Müller and colleagues (Müller et al 2018), multiple experiments conducted on the same group of subjects may constitute an invalidation of the requisite of statistical independence between the observations. For this reason, each dataset was scanned to identify such experiments (Fig. 1B). The issue was then addressed with an iterative sampling procedure, with 1000 repetitions. In each iteration, for every set of experiments conducted on the same group of subjects, only one was randomly chosen (Fig. 1C).
The resultant datasets were elaborated as in Mancuso and colleagues (Mancuso et al. 2019b). The activation map of every experiment was reconstructed feeding its list of foci to GingerALE (v. 3.0.2; https:// www. brain map. org/ ale/) (Eickhoff et al. 2012;Turkeltaub et al. 2012). In a GingerALE modelled activation map, each focus of activation is used as the center of a 3D Gaussian distribution of probability ( Fig. 1E), that estimates the voxel-wise probability of activation of the experiment given its foci (Eickhoff et al. 2009). Such maps are automatically converted to z-scores by the software. To avoid that foci close to the midline could generate unrealistic bilateral activations, we introduced an offset adjustment: If an activation fell in a range of 12 mm from the midline, voxels in the 12 contralateral mm were attenuated in the function of their distance from the midline, with values of more distant voxels being progressively reduced (Mancuso et al. 2019b).
Once obtained the modelled activations for every experiment of the list, each one of the volumes was determined as activated in a given experiment if at least 20% of its voxels were significantly activated, i.e. different from zero in the corresponding modelled activation map, thresholded at p = 0.05 (Fig. 1F). The 20% threshold is arbitrary, but we previously showed that different thresholds do not greatly affect the results (Mancuso et al. 2019b). Then, the probability of two homotopic areas to be found co-activated in our sample was computed by means of the Patel's κ, a Bayesian measure of co-activation, the calculation of which is explained in detail elsewhere (Cauda et al. 2020;Cauda et al. 2021;Patel et al. 2006;Manuello et al. 2018), thresholded at p = 0.05. In short, Patel's κ assesses the probability of a couple of areas of being more co-activated than expected given the observed activations of the two regions, and it is tested for significance with a Monte Carlo simulation. The resulting values were assigned to all the voxels of the two homotopic regions to produce an MHC map (Fig. 1G).
The steps above were repeated 1000 times, each time producing an MHC map that was calculated only with experiments with independent groups of subjects (Fig. 1H). The resulting set of 1000 MHC maps was finally averaged to obtain a map (Fig. 1I) representing the homotopic coactivation of that specific condition. The whole procedure was implemented separately for the 6 female datasets (i.e., 1 3 Fig. 1 Workflow of the calculation of the MHC map of a given domain in a given sex 1 3 5 behavioural domains plus the domain-general condition, see next subsection for males). To ensure the validity of such maps, the significance of the averages was also tested with SPM12 as one-sample t tests, with p = 0.05 and FWE correction for multiple comparisons.

Controlling for females/males imbalance
The BrainMap searches for the male sex returned a greater number of experiments compared to the same searches for females (see Results). The use of a Bayesian method as the Patel's κ should have prevented issues regarding statistical power, as significance is not directly influenced by sample size in such form of statistics (Costa et al. 2021). However, it may be argued that given the heterogeneity typical of neuroimaging results (Botvinik-Nezer et al. 2020) a wider set of experiments, reporting a larger number of foci, might cover a wider portion of the cortex by mere chance. Thus, the comparison between maps obtained from datasets with different sizes might be biased, with a comparatively wider array of co-activations in the map produced with the largest database, because of mere sampling artifacts. More importantly, the domain-general dataset was composed of a different proportion of studies belonging to the different cognitive domains. Thus, the differences between the domain-general maps of the two sexes might have been attributed to their different cognitive composition, with the maps of each sex likely more co-activated in those areas associated with the domains in which they are more often studied.
To overcome such issues, the MHC of males was calculated using an additional sampling passage (Fig. 1D). After the first step in which repeated subject groups were randomly purged from the experiment list of each specific domain, the pool of its experiments was sampled again, reducing its size to match the number of the corresponding female condition. Then, an MHC map was calculated as detailed above, and the procedure was repeated 1000 times. To obtain the male domain-general map, we proceeded in a similar way, sampling its database so each domain subset had the same number of experiments as in the female list. This produced a set of experiments that was as large as the female one and had the same proportion of investigated cognitive functions. These samples were finally averaged and tested for statistical significance as those of women.

Comparisons between maps
To compare the MHC maps of the two sexes, the average maps of each condition were subtracted as female-male, so that regions that showed stronger co-activations in the females compared to males had positive values, while the opposite situation resulted in negative values. Such calculation can be seen as a weighted winner-takes-all map, and it was shown in Mancuso et al. (Mancuso et al. 2019b) to be able to highlight statistically relevant divergences. As further validation, the 1000 maps resulting from the sampling procedure for each sex were treated as groups in a two-sample t-test, performed with SPM12 with p = 0.05 and FWE correction for multiple comparisons.

Search results
The domain-general search resulted in 963 experiments for the female gender and 2199 for the male gender, with 7186 and 18,031 total number of foci and 4796 and 8751 number of subjects, respectively. After randomly sampling them by groups (see "Methods") we obtained 1000 sets of 234 experiments for the females and 667 experiments for the males. The males dataset was then further sampled to reduce its size to that of the females. Details for the domain-specific searches are provided in Table 1.
Importantly, the two domain-general datasets are not balanced in their respective proportions of behavioral domains (Fig. 2). For instance, in both females and males, the most studied category is cognition, but in different proportions (29 vs 42%). Also, the cognitive domains relating to the category of emotions are investigated more in females (28%) than in males (15%). This unbalancement motivated us to implement the second stage of sampling on the males dataset, to match it to the size and proportions of the female one (see "Methods").
Similarly, there are similarities but also divergences between the sexes within the behavioral domains. With regards to the 'action' category, the analysis of executive tasks prevails in both sexes, 81 vs 71% for females and males, respectively. For cognitive tasks, areas concerning language are more investigated in males than females, 34 vs 26%, respectively. Concerning the cognitive domain of emotions, in females more negative emotions were analyzed than in males, 65 vs 58%, respectively. In the context of interoception, sexuality was the most investigated both in females and in males (74 vs 64%, respectively). Finally, the observation of responses to visual stimuli in the 'perception' task prevails in both females and males (35 vs 44%, respectively). Supplementary Figure S3 shows further details on behavioral sub-domains. Unfortunately, due to the very small number of experiments in some of these sub-domains, it became unfeasible to match the male datasets to the female ones at the sub-domain level.
We also investigated if the newer research was at least mitigating those biases. Plotting the number of experiments by year (Fig. 3) showed that after approximately year 2005 there was a substantial increase of female-only experiments in most domains. However, they still remain somewhat in the minority, and the relatively scarce number of post-2010 studies makes it hard to speculate if there is an ongoing   Figure S4) show that the sexes are still investigated preferentially in different domains, despite the increase of female-only studies.

Meta-analytic homotopic connectivity of the two sexes
The MHC maps of the two sexes show a series of obvious differences between both domains and sexes. Figure 4 shows the average MHC maps, while Fig. 5 illustrates their differences. The Action domain mainly activates the primary somatomotor cortex in both sexes. However, the peak is in its inferior portion in women, while in males the most active areas are those at the level of the parietal lobe or primary motor cortex. In addition, on the medial surface, the peak of homotopic co-activation in females is anterior to that of males. Cognition tasks co-activate especially the medial and lateral prefrontal cortex (mPFC and lPFC) in both sexes, but females seem to co-activate the ventral mPFC, the insula and the temporo-parietal junction more than males, while males show peaks in the cuneus and superior temporal gyrus (STG). The Emotion domain mostly involves the insula and mPFC, but the first is more engaged in females, while the two sexes show differential homotopic co-activations in different portions of the latter (more dorsal and ventral in females, more central in males). Interoception tasks reveal a marked difference between women and men, with the first co-activating the bilateral middle frontal gyrus (MFG) and a portion of mPFC, and the latter showing higher homotopic co-activations in particular in the posterior insula and the temporo-occipital cortex. Perception tends to co-activate the homotopic areas of a large part of the cortex but the female map is preferentially engaged in areas such as the anterior insula and the superior parietal cortex, while MHC of males peaks in the STG and medial Brodmann Area 6.
The domain-general condition shows strong homotopic coactivations in medial cortices in both females and males.
However, the two sexes display different landscapes of coactivations. While males' MHC map shows the typical patterns of homotopic FC with higher values in primary regions compared to associative areas, the same is not so obvious for females. Comparatively, women show more homotopic coactivations than men in TPJ and temporo-occipital cortices, premotor cortices and MFG. Details on the most co-activated homotopic regions in the two sexes for each domain  Table 2. The Jaccard indices calculated for the six couple of maps are presented in Table 3, and confirmed that the Action and Interoception domains show a lesser extent of overlap between significant homotopic coactivations between sexes.
A series of one-sample t tests (p = 0.05, FWE corrected) supported the results of the averaged maps, with the exception of the female Interoception map, which showed no significantly co-activated areas (Supplementary Figure S5). Similarly, a series of two-sample t tests confirmed the results of the female-male differences (Supplementary Figure S6).
In general, males seemed to show more widespread homotopic co-activations than females, although many of those areas had relatively weak values. Conversely, with the exception of the Interoception and Perception maps, females showed larger regions of strong homotopic co-activation. In fact, the number of non-zero voxels was always smaller in the female than the male maps. Also, the average of their values was higher in the female maps, with the exception of Interoception and Perception ones (Table 4).
To verify these findings, we calculated these two parameters in each of the 1000 MHC maps used to build the final average maps. So, we performed a series of two-sample, two-tailed t-tests on the difference of these parameters between sex for each domain. All tests were significant at the Bonferroni corrected threshold α = 0.008, except for the number of non-zero voxels between the Action maps, which was still very close to significance, and the mean of non-zero voxels of the Interoception domain, which was non-significant. We further tested the significance of these two parameters building a null model based on some sort of sex-neutral brains, obtained averaging maps sampled from both the female and male pools of MHC maps. Methodological details and results are provided in the Supplementary Materials (Table S1). Such tests confirmed that all ratios were highly higher or lower than 1, with the exception of that of the average value of non-zero voxels in the Interoception domain. Furthermore, we tested the hypothesis of a sex imbalance in the distribution of homotopic co-activation with a chi-square test. A 6✕2 (domain-by-sex) matrix has been built (Supplementary Table S2), using our difference maps as input. For each domain, if an atlas area had stronger female co-activations (F-M > 0), it was counted in the female column, otherwise (F-M < 0) it was assigned to men. Regions whose difference was 0 were discarded. The chi-square test was highly significant (Χ 2 = 34.22; p ≪ 0.01), indicating that the maps were unbalanced between sexes.

Discussion
Our method examined sex differences in homotopic co-activation and yielded a number of results, some of which were rather unexpected.
First, we found a sex bias in the Brainmap database, an observation that has sociological and metascientific relevance. In particular, there are more experiments that were performed only with male groups than with female groups. We also found a sex imbalance in the most studied cognitive domains. For example, females are more studied than males in the emotion domain, while males are more studied in the area of cognition. It is important to point out that our results  do not necessarily mean that males are more commonly studied than females in fMRI research, as it remains theoretically possible that mixed-sex experiments, which were not the object of our study, enroll large numbers of female subjects. However, if the BrainMap database can be considered representative of the literature, it seems that biases come into play when conducting an fMRI experiment on a homogeneous sex group of subjects. On a practical level, this fact could have negative consequences for neuroscience research, such as producing an incomplete or distorted overview of human brain functioning. A consequence of the present work was that we were forced to use a complicated method to circumvent the problem.
The results of such methods tend to confirm the hypothesis that women show more interhemispheric co-activation than men, as suggested by the structural findings by Ingalhalikar et al. (2014). For example, most homotopic coactivation maps show significantly higher homotopic coactivations in women than in men. This would suggest that the female brain is less lateralized and more integrated than that of males, as it indicates that the two hemispheres tend to respond to task demands activating in tandem more often in women than men. On the other hand, as predicted, the results also draw a complex landscape and are not attributable to a simple and binary model (Joel 2021). In fact, Interception and Perception maps appear to be exceptions to the notion that women exhibit more interhemispheric integration than men. In the case of the Perception domain, women show significantly less average co-activation than men, while in the case of Interoception the difference is nominally in favor of men but not statistically significant.
Furthermore, we observed that in all maps, males have a higher number of significantly homotopically co-activated areas, although with low κ-values. This could potentially suggest that men show less (or rarer) intense but more extensive co-activations than women. This finding should not be the result of database biases, as our sampling procedures reduced the male dataset to the same size as the female one. Also, it might be feared that the more widespread co-activations of the male maps were just the effect of the higher spatial variability that may arise from sampling from a larger   dataset. However, this seems not to be the case. In fact, the result of the one-sample t-tests should rule out that they were simply the result of a high variance since the t-test maps are virtually identical to those obtained with the average (Supplementary Figure S3). Furthermore, the higher number of voxels of each one of the 1000 male maps of the various domains compared to the female ones was confirmed by a series of two-sample t test. Nonetheless, due to the substantial differences between the female and male datasets, the effect of database bias cannot be completely ruled out.
Our results appear to confirm those of Zuo (Zuo et al. 2010) in reporting stronger HC for women in the posterior cingulate cortex and dorsolateral prefrontal cortex. Furthermore, the cluster in the posteromedial temporal lobe observed by Zuo as more correlated in males than in females falls into an area of greater male co-activation as we found in our results (Fig. 4). However, there are also a number of notable differences. In general, the areas of greatest difference in our results do not coincide with those indicated by Zuo (Zuo et al. 2010). For example, stronger female coactivation stands out in areas of social cognition such as the TPJ and the medial prefrontal cortex. Conversely, many primary perceptual and motor areas were found to be more co-activated in males. It would be difficult to explain these results on the basis of the sex imbalance of the cognitive subcategories. Although there is a higher percentage of social cognition studies in the female and male datasets, the difference seems too small to explain such a marked effect. Similarly, the greater proportion of male perception studies belonging to the 'audition' subdomain does not seem such as to justify the marked effect of male co-activation among bilateral STG. Curiously, this same area was found by Mancuso et al. (Mancuso et al. 2019b) as a region of significant difference between homotopic co-activation and resting state HC. Therefore, a possible explanation of the differences between our results and those of Zuo et al. could be given by the difference between HC at rest and co-activations during task. The level of convergence and divergence between co-activations and resting state connectivity is quite complex (Laird et al. 2013;Sepulcre et al. 2010;Cole et al. 2014;Torta et al. 2013) . For instance, Goparaju et al. (Goparaju et al. 2014) suggested a rewiring of network paths during task execution in spite of overall maintained connectivity and hubness. And while Smith et al. (Smith et al. 2009) revealed a strong similarity between resting state and coactivation networks, it is important to consider that these networks are not identical.
We also observe a number of differences between the maps representing the specific cognitive domains. Unfortunately, it was not possible to perform any sampling within the domain-specific datasets to balance their cognitive fingerprint, due to the extremely small number of experiments belonging to some subdomains. So, some doubts remain about the reliability of these results. Nonetheless, the male and female datasets did not appear particularly unbalanced in this regard, showing an overall similar proportion of experiments (Fig. 2). This reassures us that domain-specific differences between sexes should not be primarily driven by database biases.
In addition, a number of repeating differences emerge across domains, which can be also observed in the domaingeneral map. TPJ, anterior insula, inferior frontal gyrus and, to a certain extent ventral mPFC show greater female co-activations; STG and sensorimotor areas show generally stronger co-activations in males. These observations might be supported by anatomic data. In fact, Björnholm and colleagues (Björnholm et al. 2017) reported that the fractional anisotropy of the corpus callosum is higher in males, while that of the splenium is higher in females. Since the body of the CC connects the motor and auditory areas, this could clarify the greater male homotopic co-activation of these regions. Partly in agreement with this, Genc and colleagues (Genc et al. 2018) report that women have more density than men in the splenium and genu, which connect the TPJ and frontal areas such as insula and IFG (Chao et al. 2009). Interestingly, our results show that these areas are more coactivated in women (Fig. 4). However, it should be pointed out that the presence of sex differences in the CC is a longdebated issue (Ardekani et al. 2013;Luders et al. 2014;Allen et al. 1991;Prendergast et al. 2015;Kurth et al. 2018;Vannucci et al. 2017;Lee et al. 2009). In addition, it should be considered that the callosal connections probably have both excitatory and inhibitory functions (Bloom and Hynd 2005), so it is not possible to identify a clear relationship 1 3 between the CC structure and homotopic co-activations. On the other hand, Karolis and colleagues (2019) reported that cortical regions showing asymmetries in task-evoked activity have reduced connections with the opposite hemisphere, thus encouraging a comparison between our data and callosal connectivity. Similarly, our results are in part consistent with some proposed psychological characterization of sexes. In the healthy population, men and women show differences in social cognition as well as in the social brain (Proverbio 2021). Women have greater emotional processing capacity, recognizing facial expressions of basic emotions more accurately and faster than men. Furthermore, women seem to engage more emotional brain areas during social cognition tasks (Navarra-Ventura et al. 2018). The fact that we found areas associated with social cognition such as TPJ showing stronger homotopic co-activation might be related to these observations.

Limitation and future directions
A limitation of this study derives from the lack of balance in the number and type of experiments between the two sexes, which forced us to design a complex methodology to limit the bias. Despite our efforts, it cannot be excluded that this imbalance did not affect our results. In particular, the effect of larger co-activations in men could simply be due to the greater number of studies on the male population. Furthermore, while the domain-general datasets were counterbalanced in respect of their behavioral domains, the same was not possible at the subdomain level, as it would have produced too small samples of experiments. Therefore, the local differences between sex maps could be at least in part artefactual. We hope that, in the future, these disparities will be overcome and it will be possible to deal with more homogeneous databases.
A limitation of this study is that due to the meta-analytic nature of the data, it was not possible to assess handedness. Handedness has been shown to be related to functional connectivity (Tejavibulya et al. 2022), it is known to have a complex interplay with brain lateralization , in particular with that of language (Szaflarski et al. 2002). Left-handedness is also slightly more frequent in men (Papadatou-Paston et al. 2008), indicating an interaction between handedness, sex and interhemispheric connectivity. While we would expect the number of lefthanded individuals in our database to be low, future studies might focus on the relation between these three variables.
A possible future development could be to evaluate the strength of interhemispheric co-activations in different domains in a large sample of female and male subjects, possibly relating it to the underlying functional and/or structural homotopic connectivity.
This study was based on the hypothesis that women have more homotopic co-activations due to greater integration of their connectome. However, it appears that sex differences in connectivity are actually an effect of brain volume (Hänggi et al. 2014;Martínez et al. 2017). For our analysis this is not particularly relevant: given that women tend to have smaller brains than men, it follows that they will tend to have greater interhemispheric integration and, based on our results, stronger homotopic co-activations. Nonetheless, future research will investigate the question and evaluate the effect of brain volume on HC and homotopic co-activation mechanisms.

Conclusions
Our study provides empirical support for the hypothesis of greater interhemispheric integration in the female brain. However, there are considerable spatial differences in sex-specific differences in homotopic co-activation. In particular, it appears that primary motor and perceptual areas are more co-activated in males, in contrast to the general model. This argues for a multidimensional view of sex brain differences and suggests that the issue should be approached with more complex models than previously thought.

Supplementary Information
The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s00429-022-02572-0. Funding Open access funding provided by Università degli Studi di Torino within the CRUI-CARE Agreement. The authors have not disclosed any funding.

Competing interests
The authors declare no competing interests.
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 1 3 need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.