Genomic Signatures of Positive Selection in Human Populations of the OXT, OXTR, AVP, AVPR1A and AVR1B Gene Variants Related to the Regulation of Psychoemotional Response

The neurobiological systems of maintenance and control of behavioral responses result from natural selection. We have analyzed the selection signatures for single nucleotide variants (SNV) of the genes of oxytocin (OXT, OXTR) and vasopressin (AVP, AVPR1A, AVPR1B) systems, which are associated with the regulation of social and emotional behavior in distinct populations. The analysis was performed using original WGS (whole genome sequencing) data on Eastern Slavs (SlEast), as well as publicly available data from the 1000 Genomes Project on GBR, FIN, IBR, PUR, BEB, CHB, and ACB populations (the latter were taken as reference). To identify selection signatures, we rated the integrated haplotype scores (iHS), the numbers of segregating sites by length (nSl), and the integrated haplotype homozygosity pooled (iHH12) measures; the fixation index Fst was implemented to assess genetic differentiation between populations. We revealed that the strongest genetic differentiation of populations was found with respect to the AVPR1B gene, with the greatest differentiation observed in GRB (Fst = 0.316) and CHB (Fst = 0.325) in comparison to ACB. Also, high Fst values were found for SNVs of the AVPR1B gene rs28499431, rs33940624, rs28477649, rs3883899, and rs28452187 in most of the populations. Selection signatures have also been identified in the AVP, AVPR1A, OXT, and OXTR genes. Our analysis shows that the OXT, OXTR, AVP, AVPR1A, and AVPR1B genes were subject to positive selection in a population-specific process, which was likely contributing to the diversity of adaptive emotional response types and social function realizations.


Introduction
Identification of positive selection genomic signals is a traditional population genetics approach based on the analysis of results of whole genome sequencing (WGS) from the perspective of comparison between species and populations.The signals indicate genetic variants corresponding to different adaptation models and, therefore, designate loci of phenotypic variability, thus appearing as objects of particular interest when mapping complex traits [1][2][3].
The driving forces and direction of selection during the evolution process are not always clear.Nevertheless, the development and human populations' growth upon the settlement of Homo sapiens in various landscape and climatic zones [4] came to a realization through the accumulation of the variants that are most suitable for new environments and contributing to the reinforcement of specific patterns of cooperative behavior and offspring care [5].Many of the selected variants can affect complex endophenotypes (characteristic sets of behavioral signs with a strong genetic component), which represent constellations of normal or altered emotional, social, and cognitive functions [6] and which, among others, can be determined by gene pleiotropy, hysteresis, and other phenomena that define the complex trait's variability.
At present, there is growing interest in understanding the genetic basis for emotional and sociocognitive functions in s [7][8][9][10].With regard to social and emotional behavior, the neuropeptides oxytocin (OXT), vasopressin (AVP) and their receptors (AVPR1A, AVPR1B, and OXTR) gain particular attention [11,12] because the role of OXT and AVP in the mechanisms of establishing and maintaining social interaction, social memory, empathy, prosociality, and anxiety is well recognized [11].Mechanisms of psychoemotional adaptation are gaining attention due to the standing problem of acute and chronic stress and the aging of the general population.Therefore, both genetic factors and molecular-biological mechanisms of stress resilience and susceptibility to environmental changes evoke persistent interest.A review of the molecular mechanism of oxytocin and vasopressin action can be found in [13,14].
It is of note that the OXT and AVP genes show almost no variability in their sequences, while the genes encoding their receptors are characterized by a high degree of variability (the high number of SNVs associated with phenotypic peculiarities).Therefore, the study on variability of the AVPR1A, AVPR1B, and OXTR genes makes a great contribution to understanding the genetic architecture of behavioral traits associated with social behavior and emotions [15].It has already been shown that genetic variants in the human AVPR1A, AVPR1B, and OXTR genes are associated with the diversity of perception of social signals (empathy, emotional intelligence, eye contact, face "reading," autism), prosociality (trust, altruism, generosity, loyalty), social tolerance (aggressiveness, competition), and anxiety [7,8,16,17] that contributes to the variability of personality types appearing in human population.In turn, evaluation of the positively selected variants is of particular interest because it allows one to identify functionally significant SNVs contributing to the formation of psychopathological phenotypes.
The signaling mediated by oxytocin and vasopressin is essential not only in the regulation of physiological functions, among which labor in females and reabsorption of water in the kidney, respectively, are most prominent.Polymorphisms in genes for oxytocin and vasopressin and their specific receptors are associated with diversity in traits related to the behavioral functions of positive and negative valence and general arousal.They are significant in setting the cumulative risk for autism spectrum disorders (ASD), schizophrenia, stress susceptibility and panic disorder, major depression, and aggression [18][19][20][21][22][23][24].The high biomedical significance of both signaling systems mediating psychoemotional functions and social interactions supports a persistent interest in the genetics of these two systems.
Identification of positive selection genomic signals in distinct populations would be helpful in revealing genetic mechanisms of behavioral traits.The study was based on the assumption that diversity in the climate/geographical zones is a factor of difference in the scenarios of positive gene selection during Homo sapiens migration.Therefore, the genetic diversity within the genes of the oxytocin and vasopressin systems across different human populations was our primary interest.In this study, we identified selection signatures in the OXT, OXTR, AVP, AVPR1A, and AVPR1B genes in different human populations with population genetic models.

Characteristics of Sample
An original sample of 50 individuals from the East Slavic subethnic group (SlEast) was stratified based on the criteria, which were set upon results of the principal component analysis (PCA) of the 1000 Genomes consortium data (Figure S1) [25].In addition, we performed the cluster analysis by admixture to show the ancestry proportion of SIEast versus the 1000 Genome samples (Figure S2).
WGS data for Afro-Caribbeans in Barbados (ACB), British from England and Scotland (GBR), Finns in Finland (FIN), Iberian population in Spain (IBR), Puerto Ricans in Puerto Rico (PUR), Bengalis in Bangladesh (BEB), and Han Chinese in Beijing (CHB) (50 individuals in each ethnic group) were from the 1000 Genomes Consortium public database.Fifty individuals from each ethnic group were selected in a random order from the general pool of data.

DNA Isolation, Preparation of Sequencing Libraries and WGS of Whole Blood Samples
Isolation of genomic DNA (gDNA) from whole blood samples was performed manually using a DNA Blood Mini Kit (Qiagen, Hilden, Germany) in accordance with the manufacturer's protocol.The yield and purity of the isolated gDNA were manually determined using Qubit 4.0 fluorimeter (Thermo Fisher Scientific, Waltham, MA, USA) and NanoDropTM One C Microvolume UV-Vis (Thermo Fisher Scientific, MA, Waltham, USA), respectively.Only gDNA samples with absorbance ratios A260/280 of 1.7-1.9 and A230/260 of 1.8-2.2, in accordance, were selected for further analysis.
The amount of 150-500 ng of gDNA was used to prepare Next-Generation Sequencing (NGS) libraries.Libraries were prepared using the Illumina DNA Prep kit (Illumina, Inc., San Diego, CA, USA) according to the manufacturer's recommendations using the Tecan Freedom EVO robotic station (Tecan, Männedorf, Switzerland).The gDNA concentrations in library samples were measured using the Infinite F Nano Plus tablet reader (Tecan, Männedorf, Switzerland).The size of the resulting libraries was determined with an Agilent D1000 reagent kit using an Agilent 4200 TapeStation (Agilent Technologies, Inc., Santa Clara, CA, USA).Pooling was performed automatically using a Tecan Freedom EVO robotic station (Tecan, Männedorf, Switzerland).Each of the 50 samples of the library pool was diluted to the final gDNA concentration of 1.5 nM prior to sequencing.Pool quality control was performed with an Agilent High Sensitivity D1000 Screen Tape reagent kit using an Agilent 4200 TapeStation (Agilent Technologies, Inc., Santa Clara, CA, USA).

Population Genetic Analysis
The choice of analytical approaches was based on their specificity to spot various characteristics of the selection signal.In general, selection signals are subdivided into two types: "hard" selection, in which de novo mutations are fixed and the frequency of specific haplotypes is significantly increased in a population, therefore eliminating genetic diversity in the vicinity of the adaptive locus [26]; and "soft" selection, when selection acts toward the previously neutral alleles or allele frequency increases, resulting in a higher frequency of a specific set of haplotypes [1,25].Respectively, we used the following statistics to look for regions of recent and/or ongoing positive selection.(1) iHS (integrated haplotype score) identified loci in which a strong selection led to the fixation of new alleles and to the decrease in haplotype homozygosity (alleles that may be on their way to fixation or become a balanced polymorphism) [1].(2) nSl (number of segregating sites by length), which generally serves the same purposes as iHS, allowing one to detect hard and soft selective sweeps, appears more stable, especially with regard to the difference in the rate of recombination and the influence of population factors such as population splits, population bottlenecks (when a sharp decrease in the size of the population results in a significant drop in the genetic variability), and population growth without requiring a genetic map [8].
(3) iHH12 (integrated haplotype homozygosity pooled), which was developed based on H12 statistics [27], is able to scan soft selective sweeps by adapting H12 into an integrated haplotype homozygosity structure [28].Finally, (4) the fixation index Fst was used to assess the genetic differentiation of populations [4].Because a local adaptation of a population can appear even at conditions of low genomic differentiation, a comparison of genetic variability between two populations increases the power of detection for positive selection, in contrast to what can be revealed at the level of a single population.
Quality control, cohort pooling, and removal of duplicate SNVs were performed in accordance [29].Data for all 400 (eight populations, each represented by 50 individuals) samples were combined using the plink1.9program [30].bcftools software [31] was used to remove duplicate SNVs with identical positions.Duplicates were removed to provide effective processing on the next steps of analysis.After variant calling, variants where more than 5% or 10% of the genotype calls were not available, were removed with geno0.05 and mind0.1 tools, respectively.Hard-filtering of rare minor alleles was performed with maf0.05.Within-group data pooling was performed after removing multiallelic SNVs and correcting allele flipping.All data were normalized to the GRCh38 reference.The start and end positions for the OXT, OXTR, AVP, AVPR1A, and AVPR1B gene regions (GRCh 38 assembly) were obtained from the NCBI database (accessed on 15 August 2022) [32].
Genetic variability was assessed in every population to weight intrapopulation variability.The iHS, nSl, and iHH12 approaches for the detection of selection signals were implemented using the selscan program [33].SNVs with a minor allele frequency (MAF) > 5% were kept for analysis, with each SNV sequentially treated as a major SNV.Non-standardized scores were normalized using the "norm" script of the selscan program.Windows with crit = 1 values were considered as candidate regions for selection.
We have compared a number of different populations with the ACB population aiming at the evaluation of changes in genomes, which were accompanying the migration of Homo sapiens from Africa and its accustoming to new regions of settling.The genetic differentiation of the ACB and SlEast, ACB and GBR, A and CB and FIN, ACB and IBR, ACB and PUR, ACB and BEB, and ACB and CHB populations for the OXT, AVP, AVPR1A, AVPR1B, OXTR genes was calculated based on the pairwise Fst method proposed by Weir & Cockerham, which was implemented in the VCFTools program [34].A difference was considered significant at Fst ≥ 0.3, which corresponds to the quantile level of 0.99 [35].

Overall Estimation of Population Differentiation Based on the Analysis of Genes of Oxytocin and Vasopressin Systems
The greatest genetic population differentiation was observed with respect to the AVPR1B gene.The average Fst value for all analyzed populations was of 0.277 with the differentiation of GRB to ACB and CHB to ACB of Fst = 0.316 and Fst = 0.325, respectively.Selection signatures have also been identified in the AVP, AVPR1A, OXT, and OXTR genes.The Fst values of population differentiation averaged among all variants of the OXT, OXTR, AVP, AVPR1A, and AVPR1B genes for each individual locus are presented in Table 1.The strongest genetic population differentiation was found in the case of the AVPR1B gene.For the GRB to ACB and CHB to ACB populations, the mean Fst for all AVPR1B gene variants was 0.321.At the same time, rs28499431, rs33940624, rs28477649, rs3883899, and rs28452187 showed the highest Fst values in all populations under analysis (Table 2).As the frequency of the preferred allele rises rapidly, the classic strong directional selection signal tends to be localized to an unusually long invariable haplotype.According to the results of the iHS analysis, all SNVs of the AVPR1B gene were subject to positive selection in the GBR populations (Table 2).A similar phenomenon was observed in the FIN, IBS, and PUR populations, in which almost all SNVs of the AVPR1B gene showed a positive selection signature.In contrast, four positively selected SNVs were identified in the SlEast population, and only one SNV was identified in the BEB, CHB, and ACB populations.
Based on the results of nSl analysis, the selection was confirmed for all SNVs of the AVPR1B gene (with the exception of rs28588803 in GBR, which showed a positive signal detected by the iHS method).Additional SNVs have been identified in all populations except ACB.
Based on an assessment of iHH12, selection signals in the AVPR1B gene were identified only in the GBR population.In other populations, positive selection signals were not detected by the iHH12 method (Table 2).
Variants in the AVPR1A gene also showed significant multipopulational divergence relative to the ACB population.Two variants, rs10784339 and rs10747983, have been identified in the SlEast, GBR, IBS, and PUR populations; two variants, rs11829406 and rs11829452, in the SlEast, IBS, PUR, BEB, and CHB populations and in the SlEast, PUR, BEB, and CHB populations, respectively; and the rs10047514 variant in the GBR and FIN populations (Table 3).Selection traits in the AVPR1A gene were found only in the GBR population; two variants (rs11829406 and rs11829452) were identified by the iHS and nSL methods, and four variants (rs10747983, rs10784339, rs10877968, rs1565878685) by the nSL method (Table 3).
The analysis of Fst indicators revealed the absence of significant differentiation signals between ACB and other populations in relation to the AVP gene (Table 4).At the same time, a selection signal was found in the ACB population in the AVP gene (variant rs3787482) (Table 4).

Population Differentiation with Respect to OXT and OXTR Genes
Analysis of SNVs of the OXTR gene generally did not reveal differentiation; however, some variants showed rather high Fst values.Specifically, the levels of differentiation between ACB and CHB in relation to rs2324728l, rs237884, rs1042778, and rs237895, Fst were above 0.5, and between ACB and BEB in relation to rs2324728l and rs237884, Fst reached 0.5 (Table 5).In addition, rs59190448 has been identified in all populations except PUR.
The selection signals identified in the OXTR gene were found in all populations except CHB.The largest number of positively selected SNVs (13 SNVs) was found in the PUR population using the iHS method (Table 5).Three SNVs were identified in the SlEast and ACB populations using the iHS method.At the same time, the rs34992398 variant showed positive selection in the PUR, SlEast, and ACB populations.The rs918316 variant was detected in the GBR population using the iHS method.This variant was also present in the IBS population; however, here, it was discovered using the nSl method.Two SNVs were identified in the BEB population based on the nSl approach.Positive selection of the rs237888 variant was detected in the FIN population by the iHS and nSl methods.Analysis of Fst scores revealed the absence of significant differentiation signals between ACB and other populations in relation to the OXT gene (Table 4).Of note, a selection signal was also found in the OXT gene; the rs111869749 variant was found in the GBR population (Table 4).

Discussion
The selection of specific genetic structures in human populations occurs under the pressure of various factors.In general, it aims at the fixation of adaptive changes and, in particular, at the consolidation of complex behavioral programs, which conduce the survival of the individual and the viability of an entire population.Genetic studies make a prodigious contribution not only to understanding the evolution of behavioral traits but also to identifying functional gene variants that influence behavior such as caregiving, tribal care, and cooperation, as well as competition, aggression, and stress response.
In our study, we showed that genes of the oxytocin and vasopressin systems are objects of positive selection in the population process.The selection evolves individually in each population and provides a variety of behavioral reactions.The significance of neuropeptides OXT and AVP and their receptors OXTR, AVPR1A, and AVPR1B is determined by the role they play in the regulation of social interaction and stress response, as well as in maintaining water homeostasis and carbohydrate metabolism.Respectively, the selection in the genes of the oxytocin and vasopressin systems can be forced by a combined driving force.
We found that variants of the OXT, OXTR, AVP, AVPR1A, and AVPR1B genes showed positive selection signals in populations belonging to Indo-European ethnic groups (SlEast, GBR, FIN, IBR, BEB), in the Han subpopulation (CHB), and in Puerto Ricans of Puerto Rico (PUR), who are of a mixed ethnic origin.Overall, the genetic differentiation of the scrutinized populations was stable, as evidenced by the high Fst fixation scores for the analyzed genes in all populations.It is noteworthy that the selection of SNVs of the AVPR1B gene was the most pronounced.SNVs of other genes showed selection signs that were characteristic only for selected populations.
The observed high differentiation of populations in relation to the AVPR1B gene confirms the results of the recent study [36], which indicated a positive selection for the AVPR1B gene both in the phylogenetic lineage of human primates and among populations of White Caucasians, African Americans, Yoruba Africans, and East Asians.Taking into account the role of the AVPR1B receptor in insulin secretion [37,38], the presence of significant positive selection for this gene in representatives of GBR, FIN, IBS, and PUR may be related to the peculiarities of the regulatory hormonal mechanisms of glucose utilization in African Americans and European Americans of European descent.The most notable one is the high secretion of insulin in African Americans [38][39][40][41].
The AVPR1B receptor expression has been shown both in peripheral tissues [42] and in the central nervous system [37,43].AVPR1B receptor is mainly expressed in the pituitary gland in the brain.It mediates the effect of vasopressin on the corticoliberin-stimulated release of adrenocorticotropic hormone (ACTH) there [44] and, therefore, affects the activity of the hypothalamic-pituitary-adrenal (HPA) axis [45].One of the hypotheses of the etiopathogenesis of stress-associated diseases, which considers activation of neuroinflammation and inhibition of neuroplasticity, is based on the fact of the HPA axis dysregulation and insufficient control over the basal and stimulated cortisol levels in patients [46,47].Therefore, it is not surprising that SNVs of the AVPR1B gene are increasingly singled out in GWAS studies as candidates for assessing the risk of the development of major depression and other stress-related psychopathologies.In turn, a fixation of peculiarities of the HPA axis activity occurs at a population level.As had been shown, both basal and stress-stimulated cortisol levels were related to ethnic specificity.The flatter diurnal cortisol slopes were found among African Americans and Hispanics in comparison to Caucasians in the US [48]; in addition, stress-stimulated cortisol secretion during public speaking as well as circadian secretion pattern were reduced in African Americans compared to Caucasians [45].Since the activity of AVPR1B receptors affects the level of ACTH release, differentiation of populations based on the selection of stress responses of different severity can be mediated by the positive selection of AVPR1B gene variants.Indeed, we found that the rs33933482 variant of the AVPR1B gene, which has been shown to be associated with increased cortisol secretion during public speaking stress in Europeans [45], was positively selected in Indo-European and CHB populations.Other positively selected variants of the AVPR1B gene that we found are known to be associated with side effects of vasopressin [49], panic disorders [18,19], aggression in children [20], and autism spectrum disorders (ASD) [21].Thus, AVPR1B appears to be one of the genes, the variability of which is underlying the recent adaptation of the stress response.
AVPR1A receptor-mediated signaling is deeply involved in the regulation of physiological functions.It is determined by the effects of vasopressin on glycogenolysis in the liver, insulin secretion, contraction of vascular smooth muscle cells, water reabsorption in the renal tubules (antidiuretic activity), and platelet aggregation [50,51].The importance of the AVPR1A gene in the control of social interaction, altruistic behavior, and pairing in humans has also been documented [36,52,53].Variants of the AVPR1A gene, which were identified as positively selected in our study, are associated with heroin and general drug addiction (rs10784339) [19], mental and behavioral disorders associated with substance abuse (rs10747983) [54], diabetic status, increased glucose concentration and triglycerides in the blood, and increased body mass index [43,55].Effects of AVPR1A gene variants related to the function of the nervous system can be interpreted as those which give preference to behaviors that lead to quick results and allow a more liberal risk assessment at any activity [56].
Our analysis of the selection signatures of the AVP, AVPR1A, and AVPR1B genes shows that population differentiation is accompanied by a selection for variants, which supports variability in stress reactivity among individuals and thus determines the diversity of stressadaptation strategies.The evolutionary strategy of a positive modulation of the stress response, including positive selection in the genes of the vasopressin system, has certain benefits.The intensive response to stress has advantages from a short-term perspective since it supports a quick and balanced response to changing environmental conditions at any physiological level, including the metabolic and immune systems [36].The directed selection of vasopressin system genes may result from the relatively recent spread of human populations with settlements in strictly different ecosystems.
Our analysis showed that, in contrast to the vasopressin system, the OXT and OXTR genes are main subjects of a neutral evolution, which occurs due to genetic drift.The relative stability of the OXT and OXTR genes in the population process may be due to the exclusive role of this neuropeptide hormone in childbirth (coordinating contraction of myometrial cells), lactation, and maternal behavior.However, we have shown a high level of selection in the OXTR gene in the PUR population; partially selective signals are also found in other populations but not in CHB.Selected OXTR gene variants have been shown to be associated with changes in sensitivity to exogenous oxytocin and contractility of myometrium, as well as (rs4686302) with preterm birth [57][58][59].However, a significant interaction between the effects of rs4686302 and obesity and gestational diabetes cannot be overruled.Our analysis showed that the rs4686302 variant was identified only in the PUR population.This population is characterized by a preterm birth rate of 11.6%, which is higher than averages of 8.1% for other countries of Latin America and the Caribbean, but close to the average rate of 11.9% for African populations [60].It should be noted that the actual phenotypic significance of this variant in the PUR population is difficult to ascertain due to the combination of the multifactorial nature of preterm birth and the presence of a large percentage of the West African ethnic component in the PUR population.
The neuropeptide OXT is also involved in the regulation of behavior.It determines the degree of affiliation in interpersonal interactions by increasing the efficiency and functional significance of neural networks involved in the control of social interactions [61].We have identified several variants of the OXT and OXTR genes that were under directed selection.Among them, the rs53576 OXT gene variant showed an association with the level of empathy; OXT gene variants rs6133010 and rs53576 were shown to be associated with the level of emotional intelligence in Western Slavs [54].Carriers of the G/G genotype in the OXTR rs1042778 gene variant are less sensitive to perceived family support received in order to interfere with the effect of low maternal emotional warmth and acceptance during early life [62].The OXTR gene variant rs6770632 is associated with depressive symptoms [22].In humans, the genetic variability in the OXTR rs2254298 variant determines the volume of the amygdala, a brain structure involved in the control of fear and anxiety.In turn, the OXTR rs53576 gene variant is associated with an increased level of stress-stimulated cortisol secretion [63].According to the results of our analysis, other positively selected variants of the OXTR gene are associated with aggressive behavior and depressive states [22,23,63], alcohol abuse and associated aggressive behavior [24,64], severe schizophrenia, and autism spectrum disorders [23,24,65].Although the selection for SNVs of the OXTR gene was presumably influenced by the same factors that determined the effective differentiation of populations with respect to the AVPR1A and AVPR1B genes, it should be noted that the exceptional importance of the oxytocin system in the physiological regulation of labor set apart the corresponding genes to a conservative pool that had determined the low the degree of populations' differentiation in this gene.
Table 6 summarizes the publications on the functional significance of AVPR1A, AVPR1B, and OXTR gene SNVs.Table 6.Evidence for the functional significance of SNVs of genes of the oxytocin and vasopressin system.Abbreviations: rs, reference number of a single nucleotide variant (SNV); ACB-Afro-Caribbeans in Barbados, SlEast-Eastern Slavs, GBR-British from England and Scotland, FIN-Finns in Finland, IBR-Iberian population in Spain, PUR-Puerto Ricans in Puerto Rico, BEB-Bengalis in Bangladesh, CHB-Han Chinese in Beijing; nSl is the number of segregation sites along the length, iHS is the integrated index of haplotypes, iHH12 is the integrated index of homozygosity of haplotypes, and Fst is the index of fixation.

Conclusions
An analysis of the selection signatures of isolated genomic regions of individuals, which was carried out using general population genetic models, allowed us to show positive selection signals in the AVP, AVPR1A, and AVPR1B genes, as well as, to a lesser degree, in the OXT and OXTR genes.Positive selection in these genes was actualized with the individual scenario in each case of population differentiation.High population differentiation was revealed with respect to genes of the vasopressin system and, most pronouncedly, to the AVPR1B gene, with Fst varied from 0.168 to 0.325.The significance of the OXT and OXTR genes in the control of labor (a complex of endocrine, paracrine, and immunological mechanisms providing parturition) and the maintenance of maternal behavior was, possibly, the factor determined their greater stability in population formation, which was shown in our study.In general, our results indicate that during the formation of populations, there was a positive selection towards variants that determine the diversity of emotional stress reactions.Despite the elevated probability of accumulation of individuals prone to increased stress-reactivity upon realization of the effects of these SVNs, such selection may be an advantage since it supports immediate behavioral and physiological changes, which result from a balanced reaction of the neuronal, metabolic, and immune systems upon acute environmental stimuli.

Table 2 .
Result of selection mapping in the AVPR1B gene.Abbreviations: chr-chromosome, pos-position in the given chromosome, rs-reference number of a single nucleotide variant (SNV); ACB-Afro-Caribbeans in Barbados, SlEast-East Slavs, GBR-British from England and Scotland, FIN-Finns in Finland, IBR-Iberian population in Spain, PUR-Puerto Ricans in Puerto Rico, BEB-Bengalis in Bangladesh, CHB-Han Chinese in Beijing; color codes of cells indicate the cases of significance criterion fulfillment; nSl (orange color code)-the number of segregation sites along the length, iHS (blue color code)-the integral index of haplotypes, iHH12 (yellow color code)-the integral index of homozygosity of haplotypes, Fst (green color code)-fixation index; boxes with gray fill indicate cases of exclusion from the analysis (MAF < 0.05).

Table 1 .
Fst estimates between ACB and other global populations.Each population group is represented by 50 individual data.

Table 3 .
Result of selection mapping in the AVPR1A gene.Abbreviations: see legend to Table2.

Table 4 .
The result of selection mapping in the AVP and OXT genes.Abbreviations: see legend to Table2.

Table 5 .
The result of selection mapping in the OXTR gene.Abbreviations: see legend to Table2.