key: cord-0048793-o2gx5s90
authors: Butovskaya, Marina; Rostovtseva, Victoria; Butovskaya, Polina; Burkova, Valentina; Dronova, Daria; Filatova, Vasilisa; Sukhodolskaya, Eugenia; Vasiliev, Vasiliy; Mesa, Tania; Rosa, Araceli; Lazebny, Oleg
title: Oxytocin receptor gene polymorphism (rs53576) and digit ratio associates with aggression: comparison in seven ethnic groups
date: 2020-08-14
journal: J Physiol Anthropol
DOI: 10.1186/s40101-020-00232-y
sha: 1f36d2b597ecdc39a151fe97d19500ba32e31246
doc_id: 48793
cord_uid: o2gx5s90

BACKGROUND: The specific role of the oxytocin receptor (OXTR) gene polymorphisms in emotional support seeking, related to social norms and culturally normative behavior, has been discussed in several studies. Evidence on the association between aggression and OXTR polymorphisms has also been reported. The goal of the current study was to analyze the effect of the OXTR rs53576 polymorphism, prenatal testosterone effect (second-to-fourth digit ratio, or 2D:4D), and culture on aggression assessed with the Buss-Perry Aggression Questionnaire (BPAQ). METHODS: The data were collected in Russia and Tanzania and included seven ethnic groups of European, Asian, and African origin. The total sample included 1705 adults (837 males, 868 females). All the subjects were evaluated with the BPAQ. As a measure of prenatal androgenization, the second and fourth digits were measured directly from hand, and the digit ratios were calculated. All the participants provided buccal samples, from which genomic DNA was extracted, and the OXTR gene rs53576 polymorphism was genotyped. Statistical analysis was performed using SPSS version 23.0; the alpha level for all analyses was set at 0.05. RESULTS: The ethnic group factor was the most significant predictor of ratings on BPAQ (medium effect size for physical aggression, anger and hostility scales, and low for verbal aggression). To study the effect of sex, the OXTR polymorphism, and prenatal androgenization, we conducted the z-score transformation for BPAQ scales and 2D:4D for each ethnic group and pooled these data into new z-score variables. According to the GLM analysis after leveling the effects of culture (z-transformation), all four scales of BPAQ demonstrated association with sex (main effects), with men scoring higher on physical and verbal aggression and women scoring higher on anger and hostility. Anger and hostility scales were also associated with OXTR polymorphism and 2D:4D of the right hand. The lowest levels of anger and hostility were observed in individuals with the AA genotype, especially in men. CONCLUSIONS: Our data suggest that both oxytocin (OXTR gene polymorphism) and fetal testosterone (2D:4D) may significantly affect emotional (anger) and cognitive (hostility) aggression in humans, given the leveling the role of culture.

As the capacity to form social bonds and parental behavior, general sociality is a cornerstone of human society, deeply rooted in human evolution [1, 2] . The way of how prosocial behavior develops through interactions between culturally varying norms, social cognition, emotions, and, potentially, genes, is at the center of attention for numerous theoretical and empirical studies [3] . Through social norms, humans are intrinsically motivated to enforce rules of social co-existence and cooperation, as well as rules of punishment in the direction of those who do not behave prosocially [4] , and do not punish violators [5] . Strong tendency to follow and enforce social norms in our species turns these norms into a powerful tool for enforcing cooperation at a large scale [6] . "Culturally constructed environments create powerfuland often autocatalyticselection pressures on genes" [7] . Currently accumulated data suggest that human culture has impacted the modern human genome and its variation [8, 9] .

It may be useful to keep in mind this perspective while discussing the gene-environmental interactions in the expression of human aggression. Human aggression is multi-dimensional and may be expressed physically, verbally, emotionally (e.g., anger), or cognitively (e.g., impulsivity, hostility) [10] . There are multiple risk factors that directly or indirectly influence aggression, including age, sex, genetics, and psychopathological and environmental factors [11, 12] . Men are exhibiting more violent (physical) same-sex aggression than women in most cultures [13] [14] [15] [16] [17] , which may be attributed to the higher impulsiveness of men and a stronger fear of physical danger in women. On the other hand, the critical finding in this respect is that cultural and religious norms affect the personal potential for aggressive behavior, regulating physical or verbal expression, reducing or reinforcing their level, and regulating emotional and cognitive potential.

The role of oxytocin in human social integration, as well as aggression, deviant, and antisocial behavior, has been at the focus of attention during recent years. Oxytocin neuropeptide, produced by the hypothalamus and secreted by the pituitary gland, has been frequently mentioned as an essential regulator of mother-child attachment, social affiliation, and social bonding [18, 19] .

Oxytocin's secretion decreases anxiety and protects against stress [18] , increases empathy and trust [20, 21] , and positively affects generosity in the context of perspective-taking [22] . Some authors also report on its role in group-serving dishonesty/deception, by promoting dishonesty when the outcome favored the group to which an individual belonged [23] . Oxytocin produces its effects through oxytocin receptor (OXTR) [24] , which is present both in the brain and other body tissues and is involved in the development of the social brain [25] . Since the discovery of the OXTR gene structure [26] , numerous studies have focused on the association between OXTR gene polymorphisms and different aspects of human physiology and behavior. The genetic variation of the OXTR affects the influence of deviant peer affiliation on antisocial behavior and is associated with proactive aggression, but not with reactive aggression [27] . Variation in the OXTR gene is commonly assessed through single nucleotide polymorphisms (SNPs). The interaction effect of the two OXTR SNPs, rs1488467 and rs4564970, with alcohol consumption on trait anger was reported for the Finnish sample of men and women [28] . The A allele at rs53576 of the OXTR gene was associated with prior suicide attempts, and neither abuse history nor attachment style moderated this relationship [29] . Data on Chinese Han adolescent highlight the compound effect of stressful life events on aggression and provided evidence of the relationship between the oxytocin system and aggressive behavior [30] . High life stress during the past 12 months was associated with high levels of physical aggression and hostility in OXTR rs53576 AA carriers, but not in G carrier boys, but this association was not significant for girls. The relevance of OXTR rs237885 to aggression was also demonstrated. Mainly, rs237885 TT carriers with a history of childhood abuse had a higher risk of aggression [31] . The specific role of OXTR gene polymorphisms in emotional support seeking, related to social norms and culturally normative behavior, has been discussed in several studies [32] [33] [34] [35] [36] [37] , with OXTR rs53576 being one of the main focuses of attention. OXTR rs53576 has been related to attachment security and marital satisfaction, and GG carriers rated higher compared A carriers [38] . Li and colleagues reported a positive association between the rs53576 polymorphism and general sociality and concluded that the G allele homozygotes had higher general sociality than the A allele carriers [39] . Kim and colleagues [32] found that emotional support seeking was more evident in distressed North Americans having the G allele at the OXTR rs53576 than in those with the AA genotype, but not in Koreans, whose emotional support seeking is not a culturally normative behavior. Later, the same authors reported that Koreans with the GG genotype were more likely to use emotional suppression than those with the AA genotype and hesitated to seek emotional support in response to stress, which is more in line with normative behavior in Korea [32] . Individuals carrying the G allele of the OXTR rs53576 appeared to be more sensitive to the cultural environment [37] . The general conclusion might be that certain people regulate emotions according to cultural norms because they are biologically susceptible to be sensitive to the socio-emotional cues in a culture [36] . However, the only previous cross-sectional study, analyzing the association between rs53576 and a wide variety of emotional traits and states in a sample of young adults of various ethnicities (European, Asian, Maori/Pacific Islander, others), did not find any significant associations [40] .

Another focus of the current study is the impact of prenatal androgen exposure on the aggressive behavior and its possible association with OXTR effects. Digit ratio (2D:4D), as a possible proxy to fetal androgenization, has been discussed in numerous studies [41] [42] [43] [44] . Sexual dimorphism of the digit ratios is found across the majority of species starting from amphibians [45, 46] , although it is not unidirectional in all species [47, 48] . In humans, 2D:4D is also known to be sexually dimorphic, with males having lower ratios than females in the majority of studied populations [15, 16, [49] [50] [51] (although see [52] [53] [54] for criticism). For today, the direct evidence on the association between digit ratios and prenatal androgen concentrations in humans is still quite scarce and somewhat contradictive [55] [56] [57] . However, a part of the mechanisms underlying the impact of prenatal androgen/estrogen exposure on the formation of differences in the 2nd and 4th finger lengths has been revealed in mice [58] . The genetic basis linking gonads' and digits' development has also been investigated during recent years [59] [60] [61] . It is well-known that in men, testosterone concentrations have two peaks during the lifespan: (1) in the period of early prenatal development (peak at~24th week of gestation) [62] and (2) in puberty. The fetal and postnatal Leydig cells, which are responsible for testosterone production during prenatal and pubertal periods respectively, most likely represent morphologically and functionally different cell generations with different origins [63, 64] . This is also supported by a well-known lack of association between the 2D:4D ratio and adult testosterone levels in humans [65] [66] [67] . The latter suggests an existence of a relatively independent path for basic masculinization through exposure to androgens in utero. The period of the prenatal testosterone peak takes place at the time of the early brain maturation, which may play an important role in the formation of the sexual dimorphism of human brain structures [68] [69] [70] [71] and concomitant behavior. Of particular interest, within the scope of the current study, are findings on associations between the digit ratio and personality traits [72] , including aggression at various stages of life history [73] [74] [75] [76] [77] , risk-taking [78] [79] [80] [81] , higher physical aggression in adult men, but not women [73, 82] , as well as adolescent boys but not girls [15, 17] . Data from various studies pointed on the mutual influence of oxytocin and testosterone systems on the development of social cognition in humans [83] . In experiments with oxytocin administration versus placebo, conducted on men, individuals with a higher digit ratio (low fetal testosterone) preferred to include low-threat, rather than highthreat targets, into their groups [84] . On the contrary, men with high testosterone exposure in the uterus, after intranasal oxytocin treatment, preferred to select highthreat targets more than controls. According to other studies [85] , the administration of oxytocin increased plasma testosterone in young men. Many studies currently are taking efforts to reveal the magnitude of environmental and genetic effects on aggression [86] [87] [88] [89] [90] [91] [92] [93] , as well as antisocial behavior [94] .

It seems clear that genes may shape psychological predispositions, but also that culture might influence how these predispositions are behaviorally manifested, and culture provides particular contexts that afford opportunities and constraints for the development of psychological tendencies by presenting cultural, social norms, and limitations [34, 75, 89, 95] . Given that genes and culture may interact to produce different outcomes, we tested the three-way interaction of genes, culture, and sex in predicting direct aggression, such as physical and verbal aggression; aggressive emotional feelings, such as anger; and cognitive component, such as hostility. To this aim, we investigated the association between culture (ethnic group, religion), sex, OXTR gene SNP rs53576, digit ratio, and self-ratings on aggression assessed with the Buss-Perry Aggression Questionnaire (BPAQ) [10] .

The following hypotheses were tested: (1) selfratings on aggression in men and women are significantly influenced by cultural norms and prescriptions (socialization); (2) the OXTR rs53576 polymorphism may cause differences in BPAQ self-ratings, especially in emotional dimension of aggressive manifestations;

(3) the OXTR rs53576 effect on BPAQ ratings is sexspecific; and (4) the digit ratio associations with selfratings on BPAQ may be different in A-carriers, compared to GG OXTR rs53576 genotypes.

The sample included 1705 adults from seven ethnic groups of European, Asian, and African origin. Four ethnic groups were from Russia (Russians [89] , Tatars [75] , Ob-Ugric [96, 97] , and Buryats [75, 98] ) with the mean age of 21 years. These samples were collected between 2010 and 2019. Three ethnic groups were sampled from Tanzania between 2006 and 2014 (Hadza [16, 74, 87-89, 97, 99-102] , Datoga [16, 88, 89, 97, [99] [100] [101] [102] [103] , and Isanzu), with a mean age of 36 years. These groups are representatives of traditional African cultures. More information about these cultures may be found elsewhere [15, 16, 74, [87] [88] [89] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] .

All the Russian participants gave written informed consents, while in Tanzania, as most of the participants were illiterate there, the oral informed consents were provided, following the Declaration of Helsinki. The study was approved by the Institute of Ethnology and Anthropology of the Russian Academy of Sciences (RAS) and the Commission for Science and Technology of Tanzania (COSTECH).

All the participants provided buccal samples. Genomic DNA was isolated using Diatom DNA Prep 200 (Isogen Lab, Moscow, Russia). DNA quality from all the samples was assessed by spectrophotometer readings (A260/280) using NanoDrop.

The polymorphism of the OXTR gene (rs53576) was genotyped using Taqman 5′ exonuclease assay (Applied Biosystems). The probes for genotyping were ordered through the TaqMan SNP genotyping assays (ID: C___ 3290335_10) Applied Biosystems assay-on-demand service. The final volume was 5 μl, which contained 5 ng of genomic DNA, 2.5 μl of Taqman Master Mix, and 0.25 μl of 40 genotyping assays. Polymerase chain reaction plates were read on an ABI PRISM 7900HT instrument, and SDS v2.3 software (Applied Biosystems) was used for the genotype analysis of data.

All population statistical data processing was carried out using GenAlEx software v6.5 [106, 107] : genotype and allele frequencies, Hardy-Weinberg equilibrium (HWE) test, test of homogeneity, linkage disequilibrium test, estimations of heterozygosity, and fixation index (FST) and their significances.

Self-reported aggression was assessed with the Buss-Perry Aggression Questionnaire (BPAQ) [10] . The BPAQ includes 29 statements, grouped into four scales-physical aggression (9 items), verbal aggression (5 items), anger (7 items), and hostility (8 items)-answered on a Likert scale anchored by 1 (extremely uncharacteristic of me) and 5 (extremely characteristic of me). The translation of the BPAQ into Russian and Swahili was done, following accepted standards (translations and back translations by four bilingual assistants [108, 109] ), and translated versions have been already used in several previous studies [15-17, 74, 75, 78, 87-89, 96] .

In Russia, all the participants filled out the forms by themselves. In Tanzania, almost all respondents were illiterate; thus, they were personally interviewed in Swahili by the first author or a trained local assistant. Consequently, for Tanzanian respondents, all questions were read aloud in one-to-one dialogs, and further explanations were provided, if necessary. In Tanzania, the local assistant provided examples of actions related to the trait common in the culture to ensure that participants understood the questions. Scores on individual scales and total BPAQ scores were calculated only for respondents who answered all items. Cronbach's alpha for the total scores was 0.77 for the Russian sample, and 0.79 for the Tanzanian, and ranged between 0.60 and 0.70 for separate scales in both samples.

The second and fourth digits of the participants were measured twice with a Vernier caliper, measuring to 0.01 mm, from the basal crease to the tip of the finger. Mean values of the two measurements were used to calculate the 2D:4D ratios. In this study, we provide the data on the right-hand digit ratio (R2D:4D), given the general information that the right-hand digit ratio usually is more sexually dimorphic compared to the left hand [110] . In the cases where there was a band of creases at the base of the digit, the most proximal crease was used [111] . Participants who reported injuries or deformities of the second or fourth digits were excluded from the statistical analysis. All measurements were done by experienced anthropologists.

The 2D:4D ratios were calculated following the procedure described earlier [112] . The repeated measures of the first and second digits for the whole sample gave an intra-class correlation of 0.94. Hence, it was assumed that the differences in between-individual measurements of the 2D:4D were significantly higher than the withinindividual measurement error.

Mean and standard deviation (SD) for continuous variables was used to describe the sample's characteristics. T test was used to estimate the sex differences in ratings on BPAQ scales in each of the studied groups, as well as for estimating differences in aggression between AA and G-allele carriers. The main and interaction effects of ethnic group, sex, and OXTR genotypes on four scales of the BPAQ were calculated using multivariate analysis of variance (MANOVA). Due to the contradictory results obtained in MAN(C)OVA, the effect of the OXTR gene rs53576 on the z-transformed self-ratings on the BPAQ scales was tested with one-way analysis of variance (ANOVA), followed by the use of the Dunnett's T3 post hoc test [113, 114] . This approach allowed to solve the problems related to unequal sample sizes and unequal SDs and to test for differences among OXTR genotype groups separately in men and women. One-way ANOVA was used for analysis of the z-transformed scales of BPAQ for male and female samples to estimate the effect of OXTR genotypes. Univariate analysis of covariance (ANCOVA) was used to assess the impact of sex, digit ratio, and OXTR rs53576 on the BPAQ aggression scales after leveling the cultural effects. Linear regression was used to test the effects of the right-hand 2D:4D and rs53576 polymorphism on the BPAQ scales. Statistical analysis was performed using SPSS version 23.0 (IBM Corp., Armonk, NY, USA). The alpha level for all analyses was set at 0.05.

The OXTR rs53576 genotype and allele frequencies are presented in Table 1 . The genotype distributions by sex are presented in Fig. 1 . As seen from Table 1 , the frequencies of the alleles A and G are gradually changing from east of Asia (Buryats) to Europe (Russians and Tatars) and Africa. In the Asian populations, allele A prevails over allele G and vice versa in the rest of populations. Accordingly, genotypes AA and AG are prevailing in Buryats and Ob-Urgic people, while AG and GG genotypes in Russians, Tatars, and African populations. As χ 2 test was applied for HWE several times, the Benjamini-Hochberg approach was used, and the corrected level of significance (p < 0.05) was obtained, q = 0.007. Examining the individual p values justified HWE for all samples, including the Buryat sample, accepting its null hypothesis rejection as accidental.

Ethnic group, sex, OXTR genotype, and 2D:4D as predictors of aggression groups and by sex, are shown in Table 2 . Sexual differences were found almost in all four aggression scales and for the 2D:4D ratio (Fig. 2) . The MANOVA two-way analysis with ratings on four aggression scales as dependent variables; the ethnic group, sex, and rs53576 OXTR polymorphism as independent variables; and the two-way interaction effects of independent variables was conducted (Table 3) .

These three independent factors, along with their twoway interactions, explained between 35.4 and 18.1% of variations of scores on four BPAQ scales. The ethnic group factor was the most significant predictor (the medium effect sizes for physical aggression, anger and hostility scales, and low for verbal aggression). Post hoc Bonferroni tests revealed highly significant differences between seven tested groups on physical aggression, anger, and hostility (p ranges from 4.217E−25 to 0.001).

The main effect of sex was significant for all scales, however, with small effect sizes. The effect of OXTR rs53576 polymorphism on aggression was not significant under such conditions. However, we found that the interaction effect of ethnic group × OXTR was significant for the anger scale, and the OXTR × sex interaction for the hostility scale. In both of these cases, the effect sizes were small ( Table 3 ). The interaction effect of ethnic group × sex was significant for all scales but verbal aggression (again, the effect sizes were small) ( Table 3) .

Given the prevailing impacts of the ethnic origin and sex on the studied aggression parameters, possible effects of the OXTR gene polymorphism on aggression could appear to be masked by such a strong population-specific and sexspecific effects. Hence, to reveal the possible impact of the OXTR on aggression independently of culture and sex, we have leveled the aggression measures for all studied ethnic groups and analyzed males and females separately. Leveling the population effects was achieved through the procedure of z-transformation of the BPAQ scales for each of the ethnic groups with subsequent pooling of these data into new z-transformed variables (z-BPAQ scales). They were later tested separately for males and females using one-way ANOVA with OXTR as an independent factor. In men, the ANOVA revealed a significant effect of OXTR genotypes on BPAQ anger (F = 3.06, df1 = 2, df2 = 834, p = 0.047) and hostility (F = 3.68, df1 = 2, df2 = 834, p = 0.026) scales. The post hoc Dunnett'sT3 test (for unequal variances) has revealed that the contrasting genotypes were AA and GG (p = 0.031 and p = 0.016, accordingly). The effects were small, but significant (anger: d = 0.26 and ɳ 2 = 0.007; hostility: d = 0.28 and ɳ 2 = 0.009).

Among women, the OXTR effect was found for hostility scale only (F = 3.096, df1 = 2, df2 = 865, p = 0.046), but the differences were not significant according to the Dunnett T3 test. To sum up, the lowest levels of anger and hostility were observed in individuals with AA genotype of the OXTR rs53576 irrespective of their ethnic identity, and such association was especially pronounced in men (Fig. 3) .

Some previous findings suggested the influence of the oxytocin and testosterone systems on shaping human social cognition during early development and mutual contributions of these systems to hypo-or hypersocio-cognitive manifestations [83] . Consequently, we tested the possible interplay between effects of the OXTR polymorphism and the right-hand 2D:4D ratio (as a proxy for prenatal androgenization) after leveling the population differences in aggression ratings and digit ratio (within-population z-score standardization) ( Table 4 ). According to the ANCOVA analysis, and corresponding to what has been previously mentioned above (Table 3) , all four scales of BPAQ demonstrated association with sex (main effects): men, in comparison with women, generally scored higher on physical and verbal aggression, whereas women scored higher on anger and hostility (sex differences are statistically significant in all cases; Fig. 4 ). Physical and verbal aggression, however, were predicted only by sex, with no significant impacts of other predictors (OXTR (rs53576), R2D:4D Z ), and their interactions. However, the main effect of the R2D:4D ratio on anger was significant only at the level of trend, which occurred due to previously mentioned significant interaction effect between R2D:4D and OXTR gene polymorphism (Table 4) . To see how these two factors interacted in determining the level of anger, we divided the general sample into three parts according to the type of the OXTR rs53576 (AA, AG, GG) and ran linear regression analysis within each type, setting anger (z-score) as a response variable, and R2D:4D (z-score) with sex as predictors.

According to the results of the regression analysis, significant association between anger and 2D:4D ratio of the right hand, when controlling for sex, occurred only in carriers of the AA genotype of the OXTR gene rs53576 (N = 282, beta = 0.190, t = 3.232, p = 0.001; model R 2 = 0.051), whereas for the G allele carriers, there were no significant associations. The partial regression plots with control for sex are presented in Fig. 5 .

To sum up, these results demonstrate that OXTR rs53576 polymorphism is associated with emotional aggressive manifestations irrespective of individual's population background. G allele carriers of the OXTR were more predisposed to express anger and hostility than individuals with AA genotype, and this tendency was especially pronounced in men (Fig. 3) . At the same time, low prenatal androgen levels (high 2D:4D ratio) also led to higher levels of anger and hostility in both sexes. This result fully corresponds to the sex differences in BPAQ subscales' ratings, which demonstrate that high levels of anger and hostility are female-specific (Fig. 4) . However, the OXTR effects on aggression even overweighed those of the prenatal androgens exposure (2D:4D) effects. We Fig. 2 Distributions of BPAQ scales and right-hand 2D:4D for men and women from seven studied populations. R2D:4D-digit ratio of the right hand. Mean scores for BPAQ scales are presented as mean total scores per each scale. Inter-scale comparisons of raw total scores are not assumed, since each scale contains different number of questions came to this conclusion, since the 2D:4D effects were evident only in those individuals who did not have a genetic predisposition to high levels of anger (AA carriers) (Fig. 5 ).

In this study, we present the data on the role of cultural (ethnic-group), sex, and genetic (OXTR rs53576 genotypes) factors in human aggression. According to our knowledge, this is the first study confirming the role of the OXTR rs53576 in aggression on the individual level in adult representatives from industrial and traditional ethnic groups.

Our findings supported the hypothesis that (1) selfratings on aggression are significantly influenced by cultural norms and prescriptions. We demonstrate that ratings on direct (physical and verbal aggression) as well as emotional (anger) and cognitive (hostility) aggression scales were highly culture-specific. To this extent, culture is a "cluster of cognitions, emotions, and practices" [115] . Gender differences in the BPAQ ratings were registered both for industrial and traditional groups. The effect of sex was in the expected direction for physical aggression [13, 15, 17, [116] [117] [118] [119] [120] [121] . In Russians, Tatars, Ob-Ugric, Buryats, and Hadza, men rated higher on physical aggression. However, in other traditional African cultures, Datoga and Isanzu, the gender differences in physical aggression were not significant. These findings are important, given the general consensus that universally, this trait is more male-specific as the outcome of sexual selection and social learning [10, 118, 119, [122] [123] [124] [125] [126] . In most ethnic groups, gender differences were not found for anger and hostility, with exceptions of anger in Isanzu (higher in women), as well as anger and hostility in Buryats (higher in women). We suggest that the difference obtained in our study reflect cultural norms and socio-cultural network structure. Particularly, relatively high level of ratings on physical aggression in Datoga women may be attributed to the polygynous marriage system and cultural beliefs about individual self-esteem, along with the absence of taboo for the use MANOVA two-way analysis with the BPAQ scales ratings as dependent variables; ethnic group, sex, and OXTR gene polymorphism (rs53576) as independent factors; and the two-way interaction effects of these factors are presented F test statistics, df Degrees of freedom, p statistical significance, partial eta 2 effect size of physical aggression in conflicts among women, which is especially evident in non-kin interactions [127] . The prevailance of physical aggression in men in most still be due to gender differences in motivational sphere. Men are more motivated to revenge physically, and contra to women, such aggression may be mainly proactive, but not reactive (emotional) [125] . Our data demonstrated that groups of European, African, and Asian origin differ in the distribution of OXTR rs53576 genotypes. Buryats and Ob-Ugric, being the most distanced from the rest of our sample, demonstrated a significantly higher prevalence of the AA genotype compared to European and African samples, which is generally in line with previous findings [97] . These differences may suggest certain specific benefits for carriers of the OXTR AA genotype in populations with Asian origin as opposed to Europeans and Africans. Earlier, Kim and colleagues hypothesized that the culturally normative behavior might be subjected to genetic influence [128, 129] . According to their findings, the emotional support seeking in distress was typical for the G allele carriers North Americans, compared to AA genotype carriers, but this was not the case for Koreans [32] . At the same time, Koreans having the GG genotype were more likely to use emotional suppression compared to the AA genotype carriers [34] . Supposedly, the G carriers may be more environmentally sensitive. The question then remained whether positive selection in the direction of the AA genotype in Asian populations makes these populations generally more resistant to social-environmental pressure.

Because of cultural-specific differences in social norms and morality, obviously visible in ratings on BPAQ, we conducted the standardization procedure separately for each study sample and used these z-transformed data in further analysis. Given the general sex differences obtained, the effect of the OXTR gene polymorphism was tested separately for men and women.

Our findings fully support the hypotheses that (2) the OXTR rs53576 polymorphism may cause differences in BPAQ ratings, especially in emotional dimension of aggressive manifestations, and that (3) such effects are sexspecific. We demonstrated that the OXTR rs53576 was a significant predictor for ratings on anger and hostility scales, which was more evident for men. The GG carriers rated higher on anger (emotional aggression) and hostility (cognitive aggression) than the AA carriers. Mind that those differences were obtained after minimizing the effect of culture, by z-transformation procedure; hence, higher anger and hostility ratings reveal genetic-associated patterns, rather than cultural-specific differences. Our findings may be of special interest in light of the results, reported on the relationships between the OXTR rs53576 polymorphism and general sociality and empathetic abilities, with GG carriers being more empathetic [39, 130] . The association between sensitivity to stressful events with the rs53576 genotypes may be gender-specific. Windle and Mrug reported that young adult female (but not males) carriers of the GG genotype, who previously experienced psychological trauma in adolescence due to parental divorce, had significantly more depressive symptoms, compared to females with AA and AG genotypes [131] . The GG OXTR rs53576 genotype carriers were slower in the recognition of emotion of fear compared to the GA genotype [132] . Our findings that GG carriers rated higher on anger and hostility (this tendency was more expressed in men) extend the previous conclusion about higher threshold on emotional sensitivity and lower stress resistance of this OXTR rs53576 genotype. Assuming GG genotype carriers are less competent in recognizing of fear in others, they may be more emotionally aggressive (anger) and more hostile to others.

According to our knowledge, this is the first study reporting the interaction effect of oxytocin and prenatal androgenization (assessed through 2D:4D ratio) on the BPAQ anger and hostility ratings. The significant positive effects of the right-hand 2D:4D on anger and hostility were found only for AA carriers and were more evident for men than women. Earlier, Weisman with co-authors suggested that women are doing better in Baron-Cohen's "Reading the Mind in the Eyes" test (RMET), compared to men due to sex differences in testosterone and oxytocin levels [83] . However, in the study, the main differences were obtained for GG carriers. Particularly, men with higher 2D:4D ratio (low fetal testosterone), performed better on the RMET, possibly because they were better at understanding the emotional state of others than individuals with lower 2D:4D [83] . Our current findings do not support the hypothesis that (4) the presence of the A allele of the OXTR plays a crucial role in determining behavioral predispositions, which in turn would distinguish the GG genotype against other variants. A number of studies by other authors have previously reported that AA homozygotes have weaker functional connectivity of the hypothalamus and increased right amygdala activation [130, 133] , which most likely affect pro-social behavior and perception. This may suggest that the effects of the A and G alleles could have more or less uniform gradients, with the AA and GG genotypes representing more or less It is important to consider these results in light of some limitations. The currently accumulated data are inconsistent, obviously demanding replication in other samples with higher power, keeping in mind the small effect sizes, usually reported. Along with testing the gene-environment interaction factor, the gene-gene interaction effect should be considered; hence, the genome-wide association studies may provide much more profound information on the genetic effect on behavior. Currently, the functionality of the OXTR rs53576 polymorphism is not clear, as it is located in an untranslated region [134] . At the moment, only associative evidence for the role of OXTR in the regulation of brain activity exists. Despite the fact that neuroimaging research has revealed specific brain differences associated with rs53576, such as reduced volumes in hypothalamic gray matter and amygdala for A allele carriers [130, 133, 135] , the particular functional role of this very locus remained to be proved in the future studies.

The ethnic group factor was the most significant predictor of ratings on all four BPAQ scales. All four scales of the BPAQ demonstrated association with sex (main effects), with men scoring higher on physical and verbal aggression and women scoring higher on anger and hostility. After leveling the effects of culture (z-transformation), all four scales of the BPAQ demonstrated association with sex, with men being higher on physical and verbal aggression and women scoring higher on anger and hostility. Anger and hostility scales were also associated with OXTR polymorphism and 2D:4D of the right hand. The lowest levels of anger and hostility were observed in individuals with the AA genotype; this effect was especially obvious in men. Both oxytocin (OXTR rs53576 polymorphism) and fetal testosterone (2D:4D) significantly affected emotional (anger) and cognitive (hostility) aggression in humans, given the leveling the role of culture. 

Culture and the evolutionary process

Proximate mechanisms regulating sociality and social monogamy, in the context of evolution

How do social norms influence prosocial development?

Altruistic punishment in humans

A different kind of animal: how culture transformed our species

The cultural evolution of human nature

Culture-gene coevolution, norm-psychology and the emergence of human prosociality

How culture shaped the human genome: bringing genetics and the human sciences together

The dopamine D4 receptor gene (DRD4) modulates cultural variation in emotional experience

The aggression questionnaire

Early predictors of adolescent aggression and adult violence

Risk factors and behaviors associated with adolescent violence and aggression

Sex differences in aggression in real-world settings: a meta-analytic review

Sex differences in covert aggression among adults

Sex differences in 2D: 4D and aggression in children and adolescents from five regions of Russia

Digit ratio (2D: 4D), aggression, and dominance in the Hadza and the Datoga of Tanzania

The association between 2D: 4D ratio and aggression in children and adolescents: cross-cultural and gender differences

Social support and oxytocin interact to suppress cortisol and subjective responses to psychosocial stress

Nonsocial functions of hypothalamic oxytocin

Oxytocin and social perception: oxytocin increases perceived facial trustworthiness and attractiveness

Oxytocin increases willingness to socially share one's emotions

Oxytocin increases generosity in humans

Oxytocin promotes group-serving dishonesty

Structure and expression of a human oxytocin receptor

Oxytocin receptor dynamics in the brain across development and species

Structural organization of the human oxytocin receptor gene

Oxytocin receptor gene (OXTR) and deviant peer affiliation: a gene-environment interaction in adolescent antisocial behavior

Associations between oxytocin receptor gene (OXTR) polymorphisms and self-reported aggressive behavior and anger: interactions with alcohol consumption

Attempted suicide and oxytocin-related gene polymorphisms

Effect of the interaction between oxytocin receptor gene polymorphism (rs53576) and stressful life events on aggression in Chinese Han adolescents

Genetic variants in oxytocin receptor gene (OXTR) and childhood physical abuse collaborate to modify the risk of aggression in Chinese adolescents

Culture, distress, and oxytocin receptor polymorphism (OXTR) interact to influence emotional support seeking

Culture, serotonin receptor polymorphism and locus of attention

Gene-culture interaction: oxytocin receptor polymorphism (OXTR) and emotion regulation

Cultural neuroscience: biology of the mind in cultural context

Nature, nurture, and their interplay: a review of cultural neuroscience

Do culture and oxytocin receptor polymorphisms interact to influence emotional expressivity? Cult Brain

Associations between spouses' oxytocin receptor gene polymorphism, attachment security, and marital satisfaction

Association of oxytocin receptor gene (OXTR) rs53576 polymorphism with sociality: a meta-analysis

The oxytocin receptor gene (OXTR) variant rs53576 is not related to emotional traits or states in young adults

Fingers as a marker of prenatal androgen exposure

Fetal development of the hand, digits and digit ratio (2D: 4D)

The use of digit ratios as markers for perinatal androgen action

Age of pubertal onset and 2nd to 4th digit ratios: preliminary findings

Digit ratio in birds

Patterns of second-to-fourth digit length ratios (2D: 4D) in two species of frogs and two species of lizards at La Selva, Costa Rica

Early-to-mid gestation fetal testosterone increases right hand 2D∶ 4D finger length ratio in polycystic ovary syndrome-like monkeys

Sex differences in rhesus monkeys' digit ratio (2D: 4D ratio) and its association with maternal social dominance rank

The 2nd:4th digit ratio, sexual dimorphism, population differences, and reproductive success

Meta-analysis of digit ratio 2D:4D shows greater sex difference in the right hand

Understanding COVID-19: Digit ratio (2D: 4D) and sex differences in national case fatality rates

Differences in the 2nd to 4th digit length ratio in humans reflect shifts along the common allometric line

A comprehensive allometric analysis of 2nd digit length to 4th digit length in humans

Avoiding misinterpretation of regression lines in allometry: is sexual dimorphism in digit ratio spurious?

2nd to 4th digit ratios, fetal testosterone and estradiol

Digit ratio (2D:4D) in newborns: influences of prenatal testosterone and maternal environment

What is the evidence for a link between digit ratio (2D:4D) and direct measures of prenatal sex hormones?

Developmental basis of sexually dimorphic digit ratios

Evans DM A variant in LIN28B is associated with 2D: 4D finger-length ratio, a putative retrospective biomarker of prenatal testosterone exposure

Genetic association suggests that SMOC1 mediates between prenatal sex hormones and digit ratio

Association of CYP19A1 single-nucleotide polymorphism with digit ratio (2D: 4D) in a sample of men and women from Ningxia (China)

Involution of human fetal Leydig cells. Animmunohistochemical, ultrastructural and quantitative study

Gene expression during development of fetal and adult Leydig cells

Key factors in the regulation of fetal and postnatal Leydig cell development

Second to fourth digit length ratio (2D:4D) and adult sex hormone levels: new data and a meta-analytic review

Revisiting the relationships of 2D: 4D with androgen receptor (AR) gene and current testosterone levels: replication study and meta-analyses

No relationship between the digit ratios (2D: 4D) and salivary testosterone change: study on men under an acute exercise

Left-right asymmetries of the temporal speech areas of the human fetus

Characterization of sexual dimorphism in the human corpus callosum

Evolving knowledge of sex differences in brain structure, function, and chemistry

The role of androgen receptors in the masculinization of brain and: what we've learned from the testicular feminization mutation

Second to fourth digit ratio (2D: 4D), tomboyism, and temperament

Finger length ratio (2D: 4D) correlates with physical aggression in men but not in women

Aggression and conflict resolution among the nomadic Hadza of Tanzania as compared with their pastoralist neighbors

Aggression self-assessment indicators and its association with sex, confession and FKBP5 rs1360780 polymorphism

Exposure to fetal testosterone, aggression, and violent behavior: a meta-analysis of the 2D: 4D digit ratio

The association between fetal testosterone and violent behavior: additional evidence using the 2D: 4D digit ratio

2D:4D, self-rated aggression, risk taking and personality traits in parachutists

Further evidence of an association between low second-to-fourth digit ratio (2d: 4d) and selection for the uniformed services: a study among police personnel in Wrocław, Poland

Experimental and self-reported measures of risk taking and digit ratio (2d: 4d): evidence from a large, systematic study

2D: 4D, optimism, and risk taking

Meta-analysis of the relationship between digit-ratio 2D: 4D and aggression

The association between 2D: 4D ratio and cognitive empathy is contingent on a common polymorphism in the oxytocin receptor gene (OXTR rs53576)

Oxytocin-motivated ally selection is moderated by fetal testosterone exposure and empathic concern

Oxytocin plasma concentrations after single intranasal oxytocin administration-a study in healthy men

Herzilya series on personality and social psychology. Human aggression and violence: Causes, manifestations, and consequences

Aggression, digit ratio, and variation in the androgen receptor, serotonin transporter, and dopamine D4 receptor genes in African foragers: the Hadza

Androgen receptor gene polymorphism, aggression, and reproduction in Tanzanian foragers and pastoralists

Serotonergic gene polymorphisms (5-HTTLPR, 5 HTR1A, 5 HTR2A), and population differences in aggression: traditional (Hadza and Datoga) and industrial (Russians) populations compared

Neural mechanisms mediating gene-behavior associations in aggression

Genetics and behaviour

An integrated analysis of genes and functional pathways for aggression in human and rodent models

Aggression based genome-wide, glutamatergic, dopaminergic and neuroendocrine polygenic risk scores predict callous-unemotional traits

The association of oxytocin receptor gene (OXTR) polymorphisms antisocial behavior: a meta-analysis

Religion and well-being: the moderating role of culture and the oxytocin receptor (OXTR) gene

2d: 4d, Big fives and aggression in young men of Caucasian, Ural and Asian origin

Polymorphisms of two loci at the oxytocin receptor gene in populations of Africa, Asia and South Europe

Sexual dimorphism in facial shape of modern Buryats of Southern Siberia

Polymorphisms of dopamine receptor genes DRD2 and DRD4 in African populations of Hadza and Datoga differing in the level of culturally permitted aggression

Comparative analysis of polymorphisms of the serotonin receptor genes HTR1A, HTR2A, and HTR1B in Hadza and Datoga males

Polymorphism of 5-HTTLPR and Stin2 loci of the serotonin transporter gene in males of African ethnic populations Hadza and Datoga

Polymorphism of the two genes encoding catecholamine degradation enzymes (COMT and MAOA) in the Hadza and Datoga African ethnic populations

Aggression and polymorphisms in AR, DAT1, DRD2, and COMT genes in Datoga pastoralists of

Waist-to-hip ratio, body-mass index, age and number of children in seven traditional societies

Approach to resource management and physical strength predict differences in helping: evidence from two small-scale societies

GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update

GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research

Back-translation for cross-cultural research

Translation procedures for the cross cultural use of measurement instruments

Digit ratio (2D: 4D), dominance, reproductive success, asymmetry, and sociosexuality in the BBC Internet Study

Second to fourth digit ratio and hand skill in Austrian children

The ratio of 2nd to 4th digit length: a predictor of sperm numbers and concentrations of testosterone, luteinizing hormone and oestrogen

Comparison of post hoc tests for unequal variance

Pairwise multiple comparisons in the unequal variance case

World religions and social evolution of the old world Oikumene civilizations: a cross-cultural perspective

Does sexual selection explain human sex differences in aggression?

Reliability and validity of Chinese version of Buss-Perry aggression questionnaire in adolescents

The generalizability of the Buss-Perry aggression questionnaire

Are there sex differences in physical aggression in the elderly? Pers Individ Differ

The Buss-Perry Aggression Questionnaire: construct validity and gender invariance among Argentinean adolescents

Aggression, conflict resolution, popularity, and attitude to school in Russian adolescents

On the dimensionality of the Buss/Perry aggression questionnaire

Gender and aggression II: personal aggressiveness

Cultural and sex differences in aggression: a comparison between Japanese and Spanish students using two different inventories

Men don't just get mad; they get even: revenge but not anger mediates gender differences in physical aggression

The reliability and validity of the Buss-Perry Aggression Questionnaire (BAQ)-Turkish version

Wife battering and traditional methods of its control in contemporary Datoga pastoralists of Tanzania

Culture and social support

Culture and social support: who seeks it and why?

A common allele in the oxytocin receptor gene (OXTR) impacts prosocial temperament and human hypothalamic-limbic structure and function

Hypothesis-driven research for G× E interactions: the relationship between oxytocin, parental divorce during adolescence, and depression in young adulthood

Oxytocin receptor (OXTR) gene polymorphisms and recognition memory for emotional and neutral faces: a pilot study

Variant in OXTR gene and functional connectivity of the hypothalamus in normal subjects

Genetic and epigenetic regulatory mechanisms of the oxytocin receptor gene (OXTR) and the (clinical) implications for social behavior

Neural mechanisms of oxytocin receptor gene mediating anxiety-related temperament

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

We thank to our local collaborator, Professor Audax Mabulla, for his help and friendship during all these years of work in Tanzania. We also thank the Russian Centre for Science and Culture in Dar-es-Salaam for constant support and valuable assistance during our stay in Tanzania, the East-Siberian State Institute of Culture (VSGIK), and International Society for Human Ethology (ISHE) for collaboration. We thank all our subjects for tolerance, friendliness, and willingness to participate in the project. Many thanks for anonymous reviewers for useful comments and suggestions on the previous version of our manuscript.