Introduction

The decisions we make in our highly dynamic social world have profound implications for individual fitness, both in terms of survival chances and opportunities to reproduce successfully. For humans, as for all primates, that social context involves not only the kinds of dyadic relationships that most evolutionary and psychological studies have focused on, but also the wider social network, or community, within which these dyadic relationships are embedded. The result is a trade off between the short term benefits to be gained at the level of individual interactions and the long term benefits to be gained through being embedded in a functional group (group augmentation selection: Kokko and Johnstone 2001; Dunbar 2012). Forming a new close romantic relationship, for example, usually results in the downgrading of two other close relationships (Burton-Chellew and Dunbar 2015), and that has consequences for the benefits that these provided. Similarly, exposure to a new social environment commonly results in the replacement of up to 40% of one’s social network, including many close friendships (Saramäki et al. 2014; Roberts and Dunbar 2015), and this involves trade offs between the costs and benefits of different relationships. It has been the capacity to manage these trade offs so as to gain the benefits of living in groups while not completely foregoing short term benefits at the individual level that, in many ways, has been the key to primate, and especially human, evolutionary success.

We have previously shown that the different domains of our social world are influenced by different neurochemicals (in particular, endorphins, oxytocin and dopamine) (Pearce et al. 2017, 2018a, b). Nonetheless, it is unlikely that this involves simple genes/behaviour relationships. It is likely that any genetic effects on behaviour are mediated through other aspects of cognition. Past work has indicated that individual variation in social behaviour is likely to be influenced by both genetic variation and dispositional profiles comprising facets such as personality, attachment style (internal cognitive-affective representation with regard to close relationship partners: Brennan et al. 1998), impulsivity (the capacity to inhibit prepotent responses) and sociocognitive skills such as the ability to identify others’ emotional expressions. However, as yet, little has been done to investigate whether these putative genetic effects on social behaviour are mediated by these dispositional factors. In this paper, we explore whether an individual’s personality, attachment style and sociocognitive skills mediate genetic influences on individual differences in two key aspects of sociability: (i) dyadic relationships in the form of romantic relationships (indexed as sociosexual orientation: Penke and Asendorpf 2008) and (ii) the wider support network the individual is embedded in (indexed as the number of individuals who we turn to for help and support in times of difficulty and distress: Dunbar and Spoors 1995; Sutcliffe et al. 2012).

That aspects of personality might mediate between the genetic underpinnings of behaviour and these behavioural outcomes is suggested by evidence for both the steps in this chain. We briefly summarise this evidence, beginning with the links between aspects of personality and our outcome behaviours. We consider not just conventional personality dimensions, but related aspects of social cognition like impulsivity, attachment style and sociocognitive skills like the ability to identify emotional signals.

Extraverts have been found to have more people in each intimacy layer of their social networks (Pollet et al. 2011), and a longitudinal study looking at the effect of personality on social relationships during the transition to university found that Extraversion related positively both to network size and whether someone became romantically involved (Asendorpf and Wilpers 1998). Moreover, Agreeableness was negatively related to having conflict with peers and Conscientiousness was positively related to the frequency of contact with family members (Asendorpf and Wilpers 1998), both of which might be expected to impact on network size and maintenance. Extraversion and Agreeableness have also been associated with various dimensions of sexuality (Schmitt and Buss 2000), and, in a cross-cultural study, short-term and promiscuous mating was related positively to Openness and Extraversion, negatively to Conscientiousness and Agreeableness, and not at all to Neuroticism in both sexes (Schmitt and Shackelford 2008). In addition, individuals with partners low in Agreeableness and Conscientiousness are less satisfied with their relationship and are more likely to consider affairs (Shackelford and Goetz 2008), suggesting that scores on these traits may affect both one’s own and one’s partner’s behaviour. Other studies have suggested that the relationship between personality and sociosexual attitudes might be mediated by other physiological mechanisms such as cortisol (Wilson et al. 2015).

Impulsivity per se does not normally feature in the list of personality dimensions, but we are particularly interested in this trait because it indexes an individual’s willingness to inhibit prepotent responses in the light of their potential consequences (e.g. the effect that a one night stand might have on an existing romantic relationship). The capacity to inhibit prepotent responses in this way has been proposed as a key requisite for maintaining complex primate societies: individuals delaying or foregoing their own gratification in order to compromise with others allows greater coordination between group members (Dunbar and Shultz 2017; Launay et al. 2015). Indeed, the capacity to inhibit prepotent responses increases with brain size across primate species (Shultz and Dunbar 2010; MacLean et al., 2014). The ability to inhibit inappropriate responses or delay gratification has been linked to better social skills in both children (Carlson et al. 2002, 2004; Casey et al. 2011; Rhoades et al. 2009; Mischel and Shoda 1976) and adults (Moffitt et al. 2001; Launay et al. 2015). Several studies have found that individuals who score high on trait Impulsivity find it more difficult to inhibit prepotent responses (Aichert et al. 2012; Kam et al. 2012; Logan et al. 1997). Individuals who are highly Impulsive may be more likely to engage impulsively in short-term sexual relationships and thus have high sociosexual orientation scores. In addition, we might expect that highly impulsive individuals would be less reliable as friends and less likely to spend the considerable time investment required to create and maintain close relationships (Miritello et al. 2013; Roberts and Dunbar 2011; Sutcliffe et al. 2012; Hall 2019). Highly impulsive people might therefore be expected to have smaller support social network layers comprising close intimates (as opposed to overall network size) than those who are more dependable and consistent in servicing their friendships.

In addition to personality, Attachment Style in close friendships is also likely to influence support network size (Pearce et al. 2017; Wei et al. 2007). Attachment style ranges across two dimensions of insecurity: Anxious Attachment measures the degree to which an individual fears rejection and abandonment and doubts their own desirability as a relationship partner, whereas Avoidant Attachment characterizes the extent to which an individual seeks to emotionally distance themselves from relationships and prefers to be self-reliant due to an underlying discomfort with interpersonal intimacy and interdependence (Brennan et al. 1998; Gillath et al. 2008). We expect that those with high scores on the Avoidant dimension are likely to have fewer close friendships, whereas individuals who score highly on the Anxious dimension, with concomitant fears of abandonment, might maintain larger support networks as an insurance policy. In terms of Sociosexual Orientation, a plausible prediction would be that those scoring high on the Avoidant dimension would lean towards shorter-term relationships, and individuals with an Anxious Attachment style might be more concerned with longer-term relationships.

As well as personality and Attachment Style influencing individuals’ social behaviour, we might expect individuals who have superior sociocognitive skills to be better able to maintain larger social networks. One important skill is the ability to identify emotional expressions on people’s faces, which is an aspect of cognitive empathy (the capacity to construct a working model of the affective states of others: Reniers et al. 2009). Someone who is more skilled at identifying others’ emotive expressions might require less cognitive load, and so might be able to manage a larger number of social partners. In support of this, higher-level mentalising abilities are positively correlated with social network size, suggesting that sociocognitive abilities support the maintenance of more social ties (Launay et al. 2015; Powell et al. 2012; Stiller and Dunbar 2007). It may be that reading a potential partner’s facial expressions also plays a role in successful mate attraction, but sociocognitive skills would not be expected to influence other aspects of Sociosexual Orientation such as attitudes to promiscuity and fantasising behaviour, and so may not be associated with Sociosexual Orientation overall.

Finally, we briefly consider the other half of the chain, namely the genetic underpinnings of behaviour. There is accumulating evidence to suggest that there are genetic influences on the conventional personality dimensions, attachment style, cognitive empathy, behaviour in sexual relationships, and engagement with wider social networks (Chen and Johnson 2012; Costa et al. 2009; Denes 2015; Lee Raby et al. 2013; Lo et al. 2017; Lucht et al. 2013; Pearce et al. 2017; Rodrigues et al. 2009; Schneiderman et al. 2014; Walum et al. 2008, 2012; Way et al. 2009; Wu et al. 2016). Moreover, personality also seems to be affected by genetic variation, having a relatively high heritability (Jang et al. 1996; Kim et al. 2013; Power and Pluess 2015; Vukasović and Bratko 2015; Yamagata et al. 2006). A number of studies have identified potential candidate genes underlying individual variation in personality (Carver et al. 2015; de Moor et al. 2012; Eisenberg et al. 2007; Hamidovic et al. 2009; Kim et al. 2013; Lo et al. 2017; Lukaszewski and Roney 2011; Varga et al. 2012). Overall, there is strong evidence of a genetic basis for personality variation, as well as other aspects of disposition (attachment style and sociocognitive capacities), sociosexual behaviour and sociality beyond dyadic interactions (see Pearce et al. 2017).

Thus, there is a prima facie case for testing the hypothesis that many of the socio-cognitive aspects of our psychology mediate the previously observed relationships between the brain’s social neuropeptides and neuroendocrines and sexual and social behaviour. Indeed, the feasibility of this is supported by Creswell et al.’s (2015) finding that social temperament (negative affectivity/Neuroticism, and social inhibition) mediates the relationship between oxytocin gene variation and social network diversity and social support.

Our previous work suggests genetic correlates of attitudes and behaviour in relation to both sexual relationships and the size of an individual’s support network (Pearce et al. 2017; Pearce et al. 2018b). Here, we test whether these relationships are mediated by personality traits, Attachment Style, or sociocognitive skills. Having already reported on genetic associations with attachment style and cognitive empathy (Pearce et al. 2017), we firstly present new analyses testing for associations between personality traits (the Big Five and Impulsivity) and variation in receptor genes linked to the six neurochemicals known to play a role in social behaviour (oxytocin, vasopressin, testosterone, β-endorphin, dopamine and serotonin). We then use backwards-stepwise regression to find best-fit models predicting Sociosexual Orientation, as a measure of attitudes and behaviour in relation to sexual relationships, and Support Network Size, as a measure of social engagement beyond the dyad, from personality, attachment style, and sociocognitive skills. Finally, we explore possible mediation effects between genes found to influence each of the dispositional effects in the ‘best-fit’ models, and Sociosexual Orientation or Support Network Size. Because the two sexes differ in their reproductive strategies and behaviour (Vigil 2007; Mehta and Strough 2009; Benenson et al. 2009, 2011; Del Giudice 2011; Cross et al. 2013; Machin and Dunbar 2013; Coates 2015; Dyble et al. 2015; David-Barrett et al. 2015; Dunbar and Machin 2014; Bhattacharya et al. 2016; Dunbar 2016, 2018; Ghosh et al. 2019; Pearce et al. 2019; Archer 2019), and this may affect the way in which cognitive differences interact with these mechanisms (Feldman Barrett et al. 2000; Hall and Matsumoto 2004; Gardner and Gabriel 2004; Bell et al. 2006; Proverbio et al. 2008; Kiesow et al. 2019), we run all analyses separately for the two sexes.

Material and Methods

Participants

Healthy adults were recruited at three science festivals and a museum in the UK. Here we focus on the 757 participants (423 female) who identified as Caucasian and did not have a history of mental illness, as previously reported (Pearce et al. 2017; Pearce et al. 2018a). All participants had ≥90% genetic coverage (i.e. each participant was missing data for <10% of the SNPs), but not all participants provided full responses to the survey questions, so sample size varies between analyses: see Table 1 for survey descriptive statistics and Table S1 for genotypic frequencies for each SNP.

Table 1 Descriptive statistics (mean, standard deviation, minimum and maximum) for female and male participants

The original data are available in full in the Online Supplementary Material.

Procedure

Participants completed a survey about their attitudes and behaviour with regard to social relationships, and provided a saliva sample for DNA extraction using OrageneDNA collection kits.

As our outcome variables, we measured general attitudes and behaviours in relation to sexual relationships using the revised Sociosexual Orientation Inventory (SOI) (Penke and Asendorpf 2008) and Support Social Network Size (following Roberts and Dunbar (2015), Stiller and Dunbar (2007), Dunbar and Spoors (1995) and Hill and Dunbar (2003), defined as the number of people participants would turn to for help and support during times of difficulty and distress).

For mediation analysis, we used four separate measures: (1) the Standard Ten Item Personality Inventory (TIPI) scale (Gosling et al. 2003), which measures Openness, Conscientiousness, Extraversion, Agreeableness and Neuroticism/Emotional Stability (the ‘Big 5’: McCrae and Costa 1992) on a 7-point scale; (2) the short-form Barratt Impulsiveness Scale as a self-report measure of Impulsivity (Spinella 2007); (3) the Reading the Mind in the Eyes task (RMET, a measure of cognitive empathy: Baron-Cohen et al. 2001); and (4) the short form Experiences of Close Relationships scale as a measure of Anxious and Avoidant Attachment Styles (Wei et al. 2007), which we modified to relate to ‘close friendships’ rather than romantic relationships.

The independent variables are the SNP loci identified by Pearce et al. (2017) as being significantly correlated with our two outcome variables. Although we originally genotyped 33 SNPs from 9 genes coding for 6 neurochemicals, pruning based on linkage disequilibrium of 7 of these SNPs, and exclusion of a further 2 SNPs due to low minor allele frequencies, left 24 target SNPs in the final analyses (Pearce et al. 2017). These comprised 10 oxytocin receptor gene (OXTR) SNPs, 2 vasopressin receptor gene (AVPR1a) SNPs, 5 β-endorphin mu-opioid receptor gene (OPRM1) SNPs, 2 serotonin receptor gene (one each for HTR1A and HTR2A) SNPs, 2 dopamine receptor gene SNPs (one each for DRD1 and DRD2), 1 SNP for ANNK1 (which is downstream of DRD2 and is closely functionally associated with it), and finally 1 SNP for the androgen receptor gene (AR), to which testosterone binds. All these 24 SNPs had at least 95% coverage (for each of the SNPs <5% of participants were missing data). As reported in Pearce et al. (2017), only rs237897, rs228485, rs265981, and rs648893 showed significant deviations from Hardy-Weinberg distributions (see Table S1 for genotypic frequencies).

Analysis

We tested for associations between the 24 SNPs and both (i) the Big 5 dimensions and (ii) Impulsivity scores using PLINK version 1.9 (Chang et al. 2014; Purcell et al. 2007). Since rs6152 is located on the X chromosome and is therefore haploid in males, only additive models were run for this SNP. For the remaining 23 SNPs, we ran genotypic models that test for (1) additive effects, add, (2) deviation from additivity, domdev, and (3) genotypic effects (the combined effect of additivity and deviation from additivity, geno-2df). All models controlled for age and sex. Since the focus of this paper is on potential mediating effects, we summarise the genetic association results in the main text and Fig. 1. For the full statistical analysis, we direct the reader to the Supplementary Material and Table S2.

Fig. 1
figure 1

Heatmap showing the percentage of SNPs from different neurochemical systems (testosterone 1 SNP, β-endorphin 5 SNPs, vasopressin 2 SNPs, oxytocin 10 SNPs, dopamine 4 SNPs, serotonin 2 SNPs) that are significantly associated with the different social domains (Pearce et al. 2017): “Disposition” here follows Pearce et al.’s (2017) terminology, relating to both attachment style and empathy measured using both RMET and the Empathy Quotient, “Dyad” refers to SOI scores and relationship quality measured by Relationship Assessment Scale, and “Network” refers both to support network size and closeness to the community measured using the Inclusion of Other in Self scale) and different dimensions of the Big 5 (self-ratings of the degree to Openness, Conscientiousness, Extravertion, Agreeableness, and Neuroticism) and with self-reported Impulsivity from the current analysis. Values for dopamine were given incorrectly in Pearce et al. (2017), and are here corrected

To explore whether individual differences in disposition can predict how individuals relate to sexual relationships and their wider social networks, we ran backwards stepwise regressions in SPSS to identify best-fit models for males and females separately. We initially entered all five Big 5 personality dimensions, Impulsivity, Anxious Attachment scores, Avoidant Attachment scores, RMET scores, and age. Sex differences in all independent variables are examined in the Supplementary Material, and correlations between all variables are given separately for males and females in Table S3.

To build on the backward regression results, we then used the dispositional (independent) variables identified in the best-fit models to conduct mediation analyses using the PROCESS plug-in for SPSS (Hayes 2013). To do so, we took the independent variables from the best-fit models from the backward regressions, and identified any SNPs that were significantly associated with them in this dataset (for personality those found here, or in Pearce et al. (2017) for Attachment Style and RMET scores). We then ran separate mediation models for each of those SNPs predicting SOI or Network Size, with each independent variable from the corresponding backward regression models becoming the mediating variable in the corresponding mediation analysis. For example, if Extraversion was found in the best-fit model predicting SOI scores, we added this as the mediator in a mediation model, with the associated SNP as the independent variable and SOI scores as the dependent variable.

All mediation models included age as a covariate. Outlying RMET scores <15 (>3 SD below the mean) were excluded. Network Size of 20 were also excluded, since 20 spaces were provided for names in the survey and this created a ceiling effect: the frequency of participants naming 20 close network members was thus artificially elevated. Where appropriate, we tested both additive and domdev mediation models, and if visual inspection of graphs suggested a dominance effect we also ran a comparison between one homozygote and the carriers of the other allele e.g. GG versus A-carrier (AA/AG). As with the backwards regression models, we ran mediation analyses for males and females separately.

Despite the very large number of subjects included in this study, the number of dependent variables and SNPs in the genetics analysis inevitably raises issues of multiple testing. How to respond to these issues has been a matter of on-going debate for some considerable time. Bonferroni correction has often been the solution of choice, but it is widely considered to be overly conservative, not least because it disproportionately increases the risk of type 2 errors (accepting the null hypothesis when it should be rejected) in a way that only impossibly large samples could overcome (Nakagawa 2004; Garamszegi 2006; Narum 2006). Even professional statisticians disagree radically about its value (Rothman 1990; Schulz and Grimes 2005; Streiner and Norman 2011). Moreover, in analysing genetics data, Bonferroni correction does not take into account dependence between SNPs that are physically or functionally associated (Rice et al. 2008). One recommended solution, with a sounder statistical basis, is simply to discount the N least significant results (where N is the number of replicates) (Streiner and Norman 2011; see also Aldrich-Blake et al. 1971). Following professional advice, we used the mperm permutation function in PLINK to control for multiple testing, recognising that, in largescale genetics studies, few results will survive this correction; however, we can at least use this informally to gauge the magnitude of the relationships identified. In exploratory studies such as the present one, this is a statistically more justifiable approach because of the type 2 error problem. This issue does not affect the multiple regression or path analyses, which incorporate their own multiple testing controls.

Results

Genetic Associations

To provide the background for the mediation analyses, the broad pattern of significant associations across neurochemical/social domain categories and the different dimensions of personality is summarized as a heatmap in Fig. 1. Those for the three levels of social domain (first three columns on left) are from Pearce et al. (2017); those for the five OCEAN personality dimensions and impulsivity (last six columns on the right) are new analyses carried out for the present paper. The statistical analysis for each SNP and personality dimension is reported in full in the Electronic Supplementary Material (ESM) and Table S2.

Backwards-Stepwise Regressions

For females, backwards-stepwise regressions indicated that the best model (adjusted R2 = 0.218) predicting SOI scores included age, Extroversion, Openness, and Impulsivity (Table 2). A similar model was obtained for males (adjusted R2 = 0.196), except that it also included Agreeableness (Table 2).

Table 2 Best-fit models produced by backward regression models predicting SOI scores and Support Network size in males and females, giving the t statistics for all independent variables with degrees of freedom in parentheses

For males, Support Network size was best predicted (adjusted R2 = 0.120) by a model including age, Neuroticism and Avoidant Attachment style (Table 2). In contrast, Network size for females was best predicted (adjusted R2 = 0.110) by age, RMET, Extraversion, Agreeableness, Anxious Attachment, and Avoidant Attachment scores, although RMET and Extraversion scores did not show significant partial effects (Table 2). Thus, both age and Avoidant Attachment were negatively related to network size in both sexes, but otherwise the best models differed between the sexes. Nonetheless, the two models explained a similar proportion of the variance in network size overall: 12% for males and 11% for females.

Mediation Analyses for SOI

We report only models where significant indirect mediation effects were found. However, a number of additional models showing trend effects are given in the Electronic Supplementary Material, along with non-significant results.

For males only, Impulsivity was found to have a significant mediating effect between variation in the dopamine DRD1 SNP rs265981 and SOI when minor allele homozygotes (AA) were compared to carriers of the major allele (AG/GG): see Fig. 2. This indirect effect survived the inclusion of additional control variables: see ESM. Carriers of the G allele were found to have significantly lower Impulsivity scores than those with AA genotypes (see Fig. 2). Impulsivity then had a significant positive relationship with SOI (Fig. 2). No direct (i.e. partial) or total effects were found between DRD1 and SOI.

Fig. 2
figure 2

Mediation analysis models between DRD1 rs265981 and SOI scores, mediated via Impulsivity scores for male participants, and mean Impulsivity scores for males for different DRD1 rs265981 genotypes. Means are shown ±2 S.E. * p < 0.05, ** p < 0.005, *** p < 0.0001

In contrast, for females, significant indirect additive effects on SOI were found via Openness for the endorphin OPRM1 rs3778151 gene, the vasopressin AVPR1a rs7294536 gene and the dopamine DRD2 rs4648317 (Fig. 3). No direct effects were found between any of these SNPs and SOI. The indirect effects for OPRM1 rs3778151 and AVPR1a rs7294536 survived the inclusion of additional control variables (see ESM).

Fig. 3
figure 3

Mediation analysis models between (A) OPRM1 rs3778151, (B) AVPR1a rs7294536 and (C) DRD2 rs4648317, and SOI scores, mediated via Openness scores for female participants. * p < 0.05, ** p < 0.005, *** p < 0.0001

Mediation Analyses for Support Network Size

There were no significant models for males. In females, a significant indirect domdev effect was found between endorphin OPRM1 rs495491 and Network Size, mediated by Extraversion (Fig. 4): heterozygotes were significantly less extraverted than homozygotes, and less extraverted individuals had smaller Support Networks. Further analysis suggested that this effect might be driven by the T-allele, since T carriers were significantly less extraverted than CC homozygotes (t391 = −2.691, p = 0.007) and this contrast also creates a significant indirect effect on Network Size (Z = -1.984, p = 0.047). Although rs495491 continued to have a significant indirect effect on Network Size independently of OPRM1 rs3778151 (which itself showed a trend effect: ESM), including additional psychological traits as covariates meant that Extraversion was no longer significantly positively linked to Network Size: the partial effect of Extraversion was masked by the opposing negative effect of Avoidant Attachment (see ESM).

Fig. 4
figure 4

Mediation analysis models between OPRM1 rs495491 and Network Size, mediated via Extraversion scores for female participants, and mean Extraversion scores for females for different OPRM1 rs495491 genotypes. Means are shown ±2 S.E. * p < 0.05, ** p < 0.005, *** p < 0.0001

Discussion

As expected, we found significant genetic associations for all of the Big 5 personality dimensions and Impulsivity (see ESM for results and discussion). Moreover, we demonstrate that Sociosexual Orientation in both men and women is best predicted by Extraversion, Openness and Impulsivity: on average, individuals scoring highly on these traits were more likely to engage in, and condone, more promiscuous behaviour and casual sexual encounters. In addition, men who rated themselves as less Agreeable were also more likely to orientate towards relationships with minimal commitment, whereas Agreeableness did not appear to influence female SOI scores as strongly. These results broadly agree with the cross-cultural work of Schmitt and Shackelford (2008), except that Conscientiousness was not included in the best-fit models for either sex, perhaps due to overlap in the variance explained with the ‘premeditation’ facet of Impulsivity. The inclusion of Impulsivity in the best-fit model supports the suggestion that successful inhibition of prepotent responses is essential for the maintenance of dyadic social relationships, at least in the sexual domain.

Furthermore, we found evidence that variation in specific genes indirectly influences how individuals relate to sexual relationships through their impact on personality traits in a sex-specific way. Whereas male Sociosexual Orientation was indirectly affected by variation in the dopamine DRD1 gene through its impact on Impulsivity, this gene was found to affect Sociosexual Orientation only directly in females (see ESM). This supplements the previous finding that interaction between the dopamine DRD1 and DRD2 genes influences the age at which individuals become sexually active (Miller et al. 1999), corroborating a previous finding that dopamine receptor genes play a role in individual differences in sexual behaviour as well as impulsivity per se (Varga et al. 2012). Our findings add to these two lines of evidence by directly demonstrating that the effects of dopamine receptor gene variation on sexual cognition and behaviour are mediated by impulsivity in men. At the same time, the fact that this relates only to men indicates that a different causal pathway exists for women.

Indeed, female Sociosexual Orientation was indirectly affected by endorphin OPRM1, vasopressin AVPR1a and dopamine DRD2 genes through their independent additive effects on Openness. Although the effect of DRD2 disappeared once the effect of the other genes were partialled out (see ESM), the mediated effect for DRD2 again indicates a role for the dopaminergic system in cognition and behaviour relating to sexual relationships, but this time through its impact on a different dimension of personality. Given the role of β-endorphin in experiencing social reward (Troisi et al. 2010) and social bonding (Dunbar 2010; Machin and Dunbar 2011; Loseth et al. 2014), it is perhaps unsurprising that it affects SOI. As well as these mediated effects, variation in OPRM1 (rs2075572 and rs648893) and DRD1 (rs265981) were found to yield direct effects on relationship quality and Sociosexual Orientation, respectively (Pearce et al. 2017), again highlighting that individual differences in these reward systems impact on sexual cognition and behaviour.

Previous work has linked behaviour relating to pair-bonding, such as relationship status and quality, to AVPR1a variation in men (Walum et al. 2008) but OXTR variation in women (Walum et al. 2012). In contrast, our findings suggest that AVPR1a also plays a role in female attitudes and behaviour regarding sexual relationships, and this is supported by the findings that AVPR1a is linked to age at first sexual intercourse in both sexes (Prichard et al. 2007). We did not find any mediated effects for OXTR, perhaps suggesting that this gene influences sexual behaviour along different causal pathways that do not involve personality, attachment styles or cognitive empathy. It is plausible that the interplay between specific gene variants and personality traits identified here, and the different patterns observed between the sexes, underlie recent suggestions of different mating strategies within both sexes as well as between them (Wlodarski and Dunbar 2015).

In comparison with Sociosexual Orientation, the factors influencing Support Network size appears to differ between the sexes. Moreover, whereas the best-fit model of SOI explained ~19-22% of the variance in each sex, the corresponding models explained only ~11-12% of Support Network size variance, suggesting that factors other than the dispositional traits measured here play an important role in explaining the latter. Identifying these further factors is a potentially productive line of future research and may involve differing cognitive capacities for processing social information, which will become increasingly complex as social networks expand.

Although individuals who reported being highly avoidant in close friendships tended to have smaller support networks in both sexes, only men’s networks were negatively impacted by Neuroticism. Highly neurotic men (that is, who self-reported low emotional stability) had fewer people to turn to for help. Although OXTR, DRD1, and OPRM1 SNPs were found to be associated with Neuroticism, which in turn correlated with Network Size, these trends towards mediated effects failed to reach significance (see ESM). These trends for OXTR broadly support Creswell et al.’s (2015) finding that Neuroticism mediates the relationship between OXTR variation and social network size and diversity. However, we failed to replicate this result specifically for rs1042778. Since Avoidant Attachment showed no genetic associations with the available SNPs, we could not test for mediated effects for this variable (Pearce et al. 2017).

For women, on the other hand, Neuroticism did not seem to play a key role in Support Network Size, but women who were highly anxious in their attachment styles specifically were more likely to have larger support networks. Neuroticism and Anxious Attachment are broadly related: in the present sample, they are significantly positively related in both men and women: p < 0.0001, R2 = 0.245 and 0.132 respectively (Table S3; see also Gillath et al. 2008). However, these two variables seem to have opposite effects on Support Network size in the two sexes.

The size of women’s support networks was also found to be larger for individuals who were more Agreeable, better able to read emotional expressions in others’ faces, and more Extraverted (although RMET and Extraversion scores failed to show significant partial relationships). The endorphin OPRM1 SNP rs495491 was found to indirectly influence women’s Network Size via its effect on their degree of Extraversion, so variation in the OPRM1 gene may thus underlie, at least partially, previous findings that Extraverts have larger social networks than Introverts (Pollet et al. 2011). The association between OPRM1 and Extraversion echoes previous work demonstrating that OPRM1 variation, albeit in relation to a different SNP (rs1799971), is linked to individual differences in social hedonia, which relates to how much pleasure an individual experiences in social interactions (Troisi et al. 2010). An individual who experiences more reward through social engagement is likely to be motivated to maintain more relationships, resulting in a larger support network.

There are, inevitably, a number of limitations to this study. Firstly, notwithstanding the large sample size in this study (>750 individuals), the issue of multiple testing means that SNP-specific results presented here should be interpreted with some caution until they have been confirmed by more samples. However, discussions over many years within the genetics community also caution us against being too hasty about throwing out the baby with the bathwater. Given that it is likely that complex phenotypes such as those underlying human social cognition and behaviour will involve multiple genes, all contributing small effects, the risk of false negatives is as problematic as the risk of false positives and a balance between these always needs to be struck (Rice et al. 2008). Statistical naivety can often create more problems than it solves. Where suitable largescale databases are available, it may be possible to exploit larger sample sizes that are less prone to type 2 errors. However, few, if any, such databases provide data on the specific social measures of interest here. The main lesson is, perhaps, that those putting such databases together need to have a wider understanding of sociality beyond the simple dyad level (Pearce et al. 2017; Dunbar 2018). That said, some of the concerns in this case are mitigated by the fact that an independent sample, despite being much smaller, yielded broadly similar results to those in Fig. 1 (Pearce et al. 2018a, b). A second issue is the fact that we relied on self-report measures of personality and sociality, and these are always at risk of a subjectivity bias. Future work could incorporate more objective measures such as how many other people in a social network list the participant as a support group member or close friend (e.g. Kwak et al. 2018). However, such a method is logistically much more onerous and certainly outside the scope of studies with the present egocentric design.

Despite these caveats, overall these findings indicate that individual differences in SOI and Support Network Size are partly driven by genetic effects on personality traits, and that these pathways differ between the sexes. Although associations between these genes and human sexual behaviour have been reported previously, we add to our understanding of the complex interplay between variables influencing human sociosexual behaviour by demonstrating that, in part, these genes exert their influence indirectly via Impulsivity in men and Openness in women. Variation in OPRM1 also influences Extraversion, which is positively associated with the size of women’s Support Networks. Together, these findings add an important dimension to our understanding of the genetic and dispositional factors underlying how individuals relate to social relationships both with sexual partners and their wider social networks. The fact that there are individual differences in how these are organised both between and within the two sexes raises interesting questions about the selection processes involved. In this respect, these results parallel those reported by Wlodarski et al. (2015; see also Wlodarski and Dunbar 2015) for psychological and anatomical indices of mating strategies. Whether they are held in dynamic equilibrium by frequency-dependent selection or are under active directional selection remains an interesting question yet to be determined.