<title>Introduction</title> <p>Understanding the extent to which sex-specific processes shape human genetic diversity has long been a matter of great interest for human population geneticists ##REF##9889112##[1]##,##REF##17067791##[2]##. To date, as detailed in ##TAB##0##Table 1##, the focus has mainly been on the analysis of uniparentally inherited markers: mitochondrial DNA (mtDNA) and the non-recombining portion of the Y chromosome (NRY). A large number of studies have found that the level of differentiation was greater for the Y chromosome than for mtDNA, both at a global ##REF##9806547##[3]## and a local scale ##REF##8751879##[4]##–##REF##17208185##[11]##, for a review see ##REF##16479583##[12]##. This result has mainly been explained by patrilocality, a widespread tendency for men to stay in their birthplace while women move to their husband's house ##UREF##0##[13]## (see ##TAB##0##Table 1## for more detailed interpretations). This hypothesis of a higher migration rate of women has been especially strengthened by the comparison of patrilocal and matrilocal populations at a local scale ##REF##11528385##[14]##–##REF##16916941##[17]##. These studies have shown that in patrilocal populations, genetic differentiation is stronger among men than among women, while the reverse is observed in matrilocal populations. It is also noteworthy that the absolute difference between male and female genetic structure is more pronounced in patrilocal than in matrilocal populations ##REF##15894624##[16]##. Interestingly, while social practices seem to consistently influence the sex-specific demography at a local scale, the robustness of a sex-specific genetic structure at a global scale is still a challenging issue (see ##TAB##0##Table 1##). A recent analysis of mtDNA and NRY variation at a global scale, which used the same panel of populations for both categories of markers (an omission that was criticized in Seielstad et al.'s ##REF##9806547##[3]## study ##REF##9806532##[18]##) showed no difference between the male and female genetic structure ##REF##15378061##[19]##. Consistent with this result, an analysis of the autosomal and X-linked microsatellite markers in the HGDP-CEPH Human Genome Diversity Cell Line Panel showed no major differences between the demographic history of men and women ##REF##15601537##[20]##. The apparent paradox between local and global trends can be resolved though, since the geographical clustering of populations with potentially different lifestyles may minimize the differences in sex-specific demography at a global scale ##REF##16617372##[21]##,##REF##11420360##[22]##. It may also be that the global structure reflects more ancient, pre-agricultural, social patterns, as patrilocality may only have increased in human societies only with the recent transition to agriculture ##REF##16479583##[12]##.</p> <p>The higher differentiation level found on the NRY as compared to mtDNA at a local scale could also be the consequence of a higher effective number of women, for example through the practice of polygyny, a tendency for men (but not for women) to have multiple mates ##REF##8751879##[4]##, ##REF##12532283##[7]##, ##REF##15190128##[15]##, ##REF##12962309##[23]##–##UREF##1##[25]##, and/or through the paternal transmission of reproductive success ##REF##17208185##[11]##. However, the influence of such processes on genetic structure has often been considered as negligible, since realistic rates of polygyny cannot create large differences in male and female genetic structure ##REF##9806547##[3]##,##REF##10364534##[5]##,##REF##11528385##[14]##. Hence, until now, the effect of local social processes on male and female effective numbers has not been investigated directly, possibly because current methods fail to unravel the relative contribution of effective number and migration rate on the differentiation level ##REF##11023812##[26]##. The consequence is that the vast majority of studies fail to show whether the observed differentiation arises from sex-specific differences in migration rate, effective numbers, or both (see ##TAB##0##Table 1##). New methods need therefore to be developed in order to appreciate the relative influence of sex-biased dispersal and differences in effective numbers on genetic structure.</p> <p>Another limitation to the use of uniparentally inherited markers stems from the fact that each of them is, in effect, a single genetic locus. For that reason, we cannot test for the robustness of the sex-specific genetic structure on these markers. We cannot either rule out the possibility that mtDNA and NRY, which contain multiple linked genes, may be shaped by selection ##REF##17912352##[27]##,##REF##16645093##[28]##. This raises the question of whether results based on uniparentally inherited markers simply reflect stochastic variation, or real differences in sex-specific demography. To answer this question, we propose a novel approach based on the joint analysis of autosomal and X-linked markers. This multi-locus analysis has the potential of providing more robust information, as these markers give an independent picture of sex-specific demography. This approach also aims to disentangle the effects of sex-biased dispersal and effective numbers on genetic structure.</p> <p>In order to recognize the impact of social organization on these differences, we investigate sex-specific genetic structure in human populations of Central Asia (##FIG##0##Figure 1##), where various ethnic groups, characterized by different languages, lifestyles and social organizations, co-exist. Although all groups share a patrilocal organization, Tajiks (sedentary agriculturalists) are bilineal, i.e. they are organized into nuclear or extended families where blood links and rights of inheritance through both male and female ancestors are of equal importance, and they preferentially establish endogamous marriages with cousins. By contrast, Kazaks, Karakalpaks, Kyrgyz and Turkmen (traditionally nomadic herders) are patrilineal, i.e. they are organized into paternal descent groups (tribes, clans, lineages), and they practice exogamous marriages, in which a man chooses a bride from a different clan.</p>
<title>Methods</title> <title>DNA Samples</title> <p>We sampled 10 populations of bilineal agriculturalists and 11 populations of patrilineal herders from West Uzbekistan to East Kyrgyzstan, representing 780 healthy adult men from 5 ethnic groups (Tajiks, Kyrgyz, Karakalpaks, Kazaks, and Turkmen) (see ##TAB##1##Table 2##). The geographic distribution of the samples and information about lifestyle is provided in ##FIG##0##Figure 1##. Also living in Central Asia, Uzbeks are traditionally patrilineal herders too, but they have recently lost their traditional social organization ##REF##17208185##[11]##, and we therefore chose not to include any sample from this ethnic group for the purpose of this study. We collected ethnologic data prior to sampling, including the recent genealogy of the participants. Using this information, we retained only those individuals that were unrelated for at least two generations back in time. All individuals gave their informed consent for participation in this study. Total genomic DNA was isolated from blood samples by a standard phenol-chloroform extraction ##UREF##7##[50]##.</p> <title>Uniparentally Inherited Markers</title> <p>The mtDNA first hypervariable segment of the mtDNA control region (HVS-I) was amplified using primers L15987 (<named-content content-type="gene">5′TCAAATGGGCCTGTCCTTGTA</named-content>) and H580 (<named-content content-type="gene">5′TTGAGGAGGTAAGCTACATA</named-content>) in 18 populations out of 21 (674 individuals, see ##TAB##1##Table 2##). The amplification products were subsequently purified with the EXOSAP standard procedure. The sequence reaction was performed using primers L15925 (<named-content content-type="gene">5′TAATACACCAGTCTTGTAAAC</named-content>) and HH23 (<named-content content-type="gene">5′AATAGGGTGATAGACCTGTG</named-content>). Sequences from positions 16 024–16 391 were obtained. Eleven Y-linked microsatellite markers (see ##TAB##2##Table 3##) were genotyped in the same individuals, following the protocol described by Parkin et al. ##REF##16289902##[51]##.</p> <title>Multi-Locus Markers</title> <p>27 autosomal and 9 X-linked microsatellite markers (see ##TAB##3##Table 4##) were genotyped in the same individuals. We used the informativeness for assignment index <italic>I</italic>
<sub>n</sub>
##REF##14631557##[52]## to select subsets of microsatellite markers on the X chromosome and the autosomes from the set of markers used in Rosenberg et al.'s worldwide study ##REF##12493913##[43]##. This statistic measures the amount of information that multiallelic markers provide about individual ancestry ##REF##14631557##[52]##. This index was calculated among a subset of 14 populations, chosen from the Rosenberg et al.'s dataset ##REF##12493913##[43]## to be genetically the closest to the Central Asian populations (Balochi, Brahui, Burusho, Hazara, Pathan, Shindi, Uygur, Han, Mongola, Yakut, Adygei, Russian, Druze and Palestinian). The rationale was to infer the information provided by individual loci about ancestry from this subset of populations, and to extrapolate the results to the populations studied here. For the X chromosome data, we pooled the ‘Screening Set10’ and ‘Screening Set52’ from the HGDP-CEPH Human Genome Diversity Cell Line Panel ##REF##11954565##[53]## analyzed by Rosenberg et al. ##REF##12493913##[43]## which represented a total of 36 microsatellites. We chose 9 markers among the 11 with the highest <italic>I</italic>
<sub>n</sub>. For autosomal data, we used the ‘Screening Set10’, which represented a total of 377 microsatellites, and chose 27 markers among the 30 with the highest <italic>I</italic>
<sub>n</sub>. All markers were chosen at a minimum of 2 cM apart from each others ##REF##10961910##[54]##. PCR amplifications were performed in a 20 µl final volume composed of 1× Eppendorf buffer, 125 µM each dNTP, 0.5U Eppendorf Taq polymerase, 125 nM of each primer, and 10 ng DNA. The reactions were performed in a Eppendorf Mastercycler with an initial denaturation step at 94°C for 5 min; followed by 36 cycles at 94°C for 30 s, 55°C for 30 s, 72°C for 20 s, and 72°C for 10 min as final extension. Forward primers were fluorescently labeled and reactions were further analyzed by capillary electrophoresis (ABI 310, Applied Biosystems). We used the software package Genemarker (SoftGenetics LLC) to obtain allele sizes from the analysis of PCR products (allele calling).</p> <title>Statistical Analyses</title> <p>We calculated the total allelic richness (<italic>AR</italic>) (over all populations), the unbiased estimate of expected heterozygosity <italic>H</italic>
<sub>e</sub>
##REF##17248844##[55]##, the total number of polymorphic sites and <italic>F</italic>
<sub>ST</sub> for mtDNA using Arlequin version 3.1. ##UREF##8##[56]##. Genetic differentiation among populations for the autosomes, the X and the Y chromosome was measured both per locus and overall loci using Weir and Cockerham's <italic>F</italic>
<sub>ST</sub> estimator ##UREF##9##[57]##, as calculated in G<sc>enepop</sc> 4.0. ##UREF##10##[58]##. The 95% confidence intervals were obtained by bootstrapping over loci ##UREF##10##[58]##, using the approximate bootstrap confidence intervals (ABC) method described by DiCiccio and Efron ##UREF##11##[59]##. Isolation by distance (i.e. the correlation between the genetic and the geographic distances) was analyzed by computing the regression of pairwise <italic>F</italic>
<sub>ST</sub>/(1−<italic>F</italic>
<sub>ST</sub>) estimates between pairs of populations to the natural logarithm of their geographical distances, and rank correlations were tested using the Mantel permutation procedure ##REF##6018555##[60]##, as implemented in G<sc>enepop</sc> 4.0. ##UREF##10##[58]##. All other statistical tests were performed using the software package R v. 2.2.1 ##UREF##12##[61]##.</p> <title>Sex-Biased Dispersal in the Island Model</title> <p>Let us consider an infinite island model of population structure ##REF##17246615##[62]##, with two classes of individuals (males and females), which describes a infinite set of populations with constant and equal sizes that are connected by gene flow. Then the expected values of <italic>F</italic>
<sub>ST</sub> for uniparentally inherited markers depend on the effective number <italic>N</italic>
<sub>m</sub> (resp. <italic>N</italic>
<sub>f</sub>) of adult males (resp. females) per population and the migration rate <italic>m</italic>
<sub>m</sub> (resp. <italic>m</italic>
<sub>f</sub>) of males (resp. females) per generation, as:  (see, e.g., ##REF##17976181##[63]##). We can therefore calculate the female-to-male ratio of the effective number of migrants per generation as: .</p> <p>In this model, we can also compute for the autosomes and the X chromosome the reproductive values for each class (sex), which are interpreted here as the probability that an ancestral gene lineage was in a given class in a distant past ##UREF##13##[64]##. From these, we can obtain the well-known expressions of effective size <italic>N</italic>
<sub>e</sub> for autosomal and X-linked genes: , respectively ##UREF##5##[45]##. Note that <italic>N</italic>
<sub>e</sub> is expressed here as a number of gene copies (i.e., twice the effective number of diploid individuals for autosomes). Likewise, the effective migration rate, i.e. the average dispersal rate of an ancestral gene lineage, is given by  for autosomal genes, and  for X-linked genes, respectively. Substituting these expressions into the well-known equation: <italic>F</italic>
<sub>ST</sub>≈1/(1+2<italic>N</italic>
<sub>e</sub>
<italic>m</italic>
<sub>e</sub>) ##UREF##13##[64]##, we get:for autosomal genes, andfor X-linked genes.</p> <title>Evaluation of the Approach through Stochastic Simulations</title> <p>We performed coalescent simulations, using an algorithm in which coalescence and migration events are considered generation-by-generation until the common ancestor of the whole sample has been reached (see ##REF##12654930##[65]##). We simulated a finite island model with 50 demes, each made of <italic>N</italic> = <italic>N</italic>
<sub>f</sub>+<italic>N</italic>
<sub>m</sub> = 500 diploid individuals, with a migration parameter <italic>m</italic> = <italic>m</italic>
<sub>f</sub>+<italic>m</italic>
<sub>m</sub> = 0.2. Using these total values for <italic>N</italic> and <italic>m</italic>, we then varied the sex-specific parameters to cover the (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) parameter space evenly. Note that the parameter <italic>m</italic> is the total migration rate, which corresponds to twice the effective migration rate for autosomal markers. Hence, for each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, the total number of individuals is 500 (although the number of females may vary from 1 to 499) and the effective migration rate for autosomal markers is . We chose these total values for <italic>N</italic> and <italic>m</italic> such that, for a ratio <italic>N</italic>
<sub>f</sub>
<italic>m</italic>
<sub>f</sub>/<italic>N</italic>
<sub>m</sub>
<italic>m</italic>
<sub>m</sub> = 21.6 (as observed for the herder populations), the distribution of <italic>F</italic>
<sub>ST</sub> estimates on uniparentally-inherited markers in the simulations were close to the observations: for mtDNA, the 95% highest posterior density interval (see ##UREF##14##[66]##, pp. 38–39) for the distribution of <italic>F</italic>
<sub>ST</sub> estimates in the simulations was [0.007; 0.033] with a mode at 0.014 (estimated value from the real dataset:  among the herders) while for the NRY, the 95% highest posterior density interval was [0.088; 0.374] with a mode at 0.187 (estimated value from the real dataset: ).</p> <p>Each simulated sample consisted in 330 sampled males from 11 populations (30 males per population), genotyped at 27 autosomal, 9 X-linked markers as well as 10 Y-linked markers and a single mtDNA locus. Each locus was assumed to follow a Generalized Stepwise Model (GSM) ##REF##12207711##[67]## with a possible range of 40 contiguous allelic states, except the mtDNA, which was assumed to follow an infinite allele model of mutation. The average mutation rate was 5.10<sup>−3</sup>, and the mean parameter of the geometric distribution of the mutation step lengths for microsatellites was set to 0.2 ##REF##12207711##[67]##,##REF##8600387##[68]##.</p>
<title>Results/Discussion</title> <title>Uniparentally-Inherited Markers</title> <p>We sampled 780 healthy adult men from 10 populations of bilineal agriculturalists and 11 populations of patrilineal herders from West Uzbekistan to East Kyrgyzstan, representing 5 ethnic groups (Tajiks, Kyrgyz, Karakalpaks, Kazaks, and Turkmen) (see ##FIG##0##Figure 1## and ##TAB##1##Table 2##). We genotyped all bilineal populations, and 8 out of 11 patrilineal populations at the HVS-I locus of mtDNA, and at 11 microsatellite markers on the NRY (for more details on the markers used, see ##TAB##2##Table 3##). The overall genetic differentiation was higher for NRY, as compared to mtDNA, both among the 10 bilineal agriculturalist populations , and among the subset of 8 patrilineal herder populations . Assuming an island model of population structure, this implies that female migration rate (<italic>m</italic>
<sub>f</sub>), and/or the effective number of females (<italic>N</italic>
<sub>f</sub>), is higher than of the corresponding parameters for males (<italic>m</italic>
<sub>m</sub> and <italic>N</italic>
<sub>m</sub>). These results also suggest that the differences in sex-specific genetic structure are much more pronounced in the patrilineal herders than in the bilineal agriculturalists. From the above <italic>F</italic>
<sub>ST</sub> estimates, we obtained the female-to-male ratio of the effective number of migrants per generation (see the <xref ref-type="sec" rid="s3">Methods</xref> section for details): <italic>N</italic>
<sub>f</sub>
<italic>m</italic>
<sub>f</sub>/<italic>N</italic>
<sub>m</sub>
<italic>m</italic>
<sub>m</sub>≈2.1 for bilineal populations and <italic>N</italic>
<sub>f</sub>
<italic>m</italic>
<sub>f</sub>/<italic>N</italic>
<sub>m</sub>
<italic>m</italic>
<sub>m</sub>≈21.6 for patrilineal populations. The ratio in patrilineal populations is thus one order of magnitude higher than in bilineal populations. However, since each of these markers is a single genetic locus, we cannot test for the robustness of the sex-specific genetic structure on these markers. We therefore examined the amount of information contained in multi-locus data on autosomal and X-linked markers, both of which average over male and female histories.</p> <title>A New Multi-Locus Approach</title> <p>In the infinite island model of population structure with two classes of individuals (males and females), we obtained the following expressions of <italic>F</italic>
<sub>ST</sub> (see the <xref ref-type="sec" rid="s3">Methods</xref> section for details):for autosomal genes, andfor X-linked genes. A special case of interest occurs when , i.e. when the differentiation of X-linked genes exactly equals that of autosomal genes. Combining eqs (1) and (2), we find that this occurs for , with <italic>N</italic> = <italic>N</italic>
<sub>f</sub>+<italic>N</italic>
<sub>m</sub> and <italic>m</italic> = <italic>m</italic>
<sub>f</sub>+<italic>m</italic>
<sub>m</sub>. Furthermore, as shown in ##FIG##1##Figure 2##, if we observe a lower genetic differentiation of autosomal markers, as compared to X-linked markers (blue zone in ##FIG##1##Figure 2##), this suggests that . This may happen, e.g., for <italic>N</italic>
<sub>f</sub> = <italic>N</italic>
<sub>m</sub> and <italic>m</italic>
<sub>f</sub> = <italic>m</italic>
<sub>m</sub>, i.e. for equal effective numbers of males and females and unbiased dispersal. But if autosomal markers are more differentiated than X-linked markers (, see the red upper-right triangle in ##FIG##1##Figure 2##), this implies that . In this case, since <italic>m</italic>
<sub>f</sub>/<italic>m</italic> and <italic>N</italic>
<sub>f</sub>/<italic>N</italic> are ratios varying between 0 and 1, the effective number of females must be higher than that of males (<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub>), and the female migration rate must be higher than half the male migration rate (<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub>/2). Hence, a prediction from this model is that when , the effective number of females is higher than that of males, whatever the pattern of sex-specific dispersal. This suggests that it is indeed possible to test for differences in effective numbers between males and females from the joint analysis of autosomal and X-linked data. We note however that when , we cannot conclude on the relative male and female effective numbers and migration rates.</p> <p>We tested the above prediction in the 10 bilineal agriculturalist populations and 11 patrilineal herder populations sampled in Central Asia by comparing the genetic structure estimated from 27 unlinked polymorphic autosomal microsatellite markers (<italic>AR</italic> = 16.2, <italic>H</italic>
<sub>e</sub> = 0.803 on average) to that from 9 unlinked polymorphic X-linked microsatellite markers (<italic>AR</italic> = 12.6, <italic>H</italic>
<sub>e</sub> = 0.752 on average) (for more details on the markers used, see ##TAB##3##Table 4##). Overall heterozygosity was not significantly different between X-linked and autosomal markers, neither in the pooled sample (two-tailed Wilcoxon sum rank test; <italic>p</italic> = 0.09), nor in the bilineal agriculturalists (<italic>p</italic> = 0.13) or the patrilineal herders (<italic>p</italic> = 0.12). The overall population structure was significantly higher for autosomal as compared to X-linked markers among patrilineal herders:  (one-tailed Wilcoxon sum rank test; ; <italic>p</italic> = 0.02). Among bilineal agriculturalists, the result was not significant:  (<italic>p</italic> = 0.36). From these results, and following our model predictions, we conclude that in patrilineal herders (where ), the effective number of females is higher than that of males. This conclusion does not hold for the bilineal agriculturalists.</p> <p>From our model, it is possible to get more precise indications on the sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values that are compatible with our data. Rearranging eqs (1–2), we get:i.e.:
</p> <p>For any given set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, we can therefore calculate from eq. (4) the expected value of  for each  estimate in the dataset. We can then test the null hypothesis  by comparing the distribution of observed and expected  values. If the hypothesis can be rejected at the <italic>α</italic> = 0.05 level, then the corresponding set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values can also be rejected. Following Ramachandran et al. ##REF##15601537##[20]##, we varied the values of the ratios <italic>N</italic>
<sub>f</sub>/<italic>N</italic> and <italic>m</italic>
<sub>f</sub>/<italic>m</italic> (respectively, the female fraction of effective number, and the female fraction of the total migration rate) from 0 to 1, with an interval of 0.01 between consecutive values. For each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, we applied the transformation in eq. (4) to each of the 27 locus-specific  values observed. Thus, for each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, we obtained 27 expected values of , given our data. These expected values of  were then compared to the 9 observed locus-specific  in our dataset, and we calculated the <italic>p</italic>-value for a two-sided Wilcoxon sum rank test between the list of 27 expected  values and the 9  observed in the dataset. The results are depicted in ##FIG##2##Figure 3##. Significant <italic>p</italic>-values (<italic>p</italic>≤0.05) correspond to a significant difference between the observed and expected values, thus to sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values that are rejected, given our data (see the blue region in ##FIG##2##Figure 3##). Conversely, non-significant <italic>p</italic>-values (<italic>p</italic>&gt;0.05) correspond to sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values that cannot be rejected (see the red region in ##FIG##2##Figure 3##).</p> <p>For the patrilineal herder populations (##FIG##2##Figures 3A–3B##), most sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values are rejected, except those corresponding to larger effective numbers for females (from ##FIG##2##Figures 3A–3B##: <italic>N</italic>
<sub>f</sub>/<italic>N</italic>&gt;0.55, i.e. <italic>N</italic>
<sub>f</sub>&gt;1.27<italic>N</italic>
<sub>m</sub>) and <italic>m</italic>
<sub>f</sub>&gt;0.67<italic>m</italic>
<sub>m</sub>. Because the multi-locus estimate of  is significantly higher than the estimate of , we expected to find such patterns of non-significant values (see ##FIG##1##Figure 2##). For the bilineal agriculturalist populations, we could not reject the hypothesis that the effective numbers and migration rates are equal across males and females or even lower in females (see ##FIG##2##Figures 3C–3D##). This is also reflected by the fact that the estimates of  were not significantly higher than the estimates of  in those populations.</p> <p>Finally, we have shown that the effective number of women is higher than that of men among patrilineal herders, but not necessarily among bilineal agriculturalists. Furthermore, a close inspection of the results depicted in ##FIG##2##Figures 3A and 3B## reveals that, among herders, we reject all the sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values for which <italic>m</italic>
<sub>f</sub>&lt;<italic>m</italic>
<sub>m</sub> at the <italic>α</italic> = 0.10 level. This is not true for agriculturalists. This suggests that the migration rates are also likely to be higher for women than for men in patrilineal populations, as compared to bilineal populations (compare ##FIG##2##Figures 3B and 3D##). Although both groups are patrilocal, such a difference in sex-specific migration patterns might be expected, since patrilineal herders are exogamous (among clans) and bilineal agriculturalists are preferentially endogamous. For example, it was observed that in patrilocal and matrilocal Indian populations, where migrations are strictly confined within endogamous groups, sex-specific patterns were not influenced by post-marital residence ##REF##16617372##[21]##.</p> <title>What Could Explain a Larger Effective Number of Females?</title> <p>While an influence of post-marital residence on the migration rate of women and men has already been widely proposed ##REF##11528385##[14]##–##REF##16916941##[17]## (see also ##TAB##0##Table 1##), the factors that may locally affect the effective number of women, relatively to that of men, are not well recognized. As seen in ##TAB##0##Table 1##, although a number of studies have compared matrilocal and patrilocal populations, few have compared contrasting groups of populations with respect to other factors as, e.g., the tendency for polygyny ##REF##15190128##[15]##. Furthermore, a number of these studies lack ethnological information a priori, concerning social organization, marriage rules, etc., which makes interpretation somewhat difficult (see ##TAB##0##Table 1##). Here, we compared two groups of patrilocal populations with contrasting social organizations, and at least five non-mutually exclusive interpretations for a larger effective number of females can be invoked:</p> <p>
<italic>Social organization</italic>, i.e. the way children are affiliated to their parents, can deeply affect sex-specific genetic variation. In Central Asia, herder populations are organized in patrilineal descent groups (tribes, clans, lineages). This implies that children are systematically affiliated with the descent groups of the father. Chaix et al. ##REF##17208185##[11]## showed that the average number of individuals carrying the same Y chromosome haplotype was much higher in patrilineal herder populations than in bilineal agriculturalist populations (where children are affiliated both to the mother and the father). These “identity cores” would be the direct consequence of the internal dynamics of their patrilineal organization. Indeed, the descent groups are not formed randomly and related men tend to cluster together, e.g. through the recurrent lineal fission of one population into new groups. This particular dynamics increases relatedness among men, and may therefore reduce the effective number of men, as compared to women.</p> <p>Indirectly, the social organization can also deflate the effective number of men through <italic>the transmission of reproductive success</italic>
##REF##15797619##[29]## if this success is culturally transmitted exclusively from fathers to sons. Because herders are patrilineal (so that inheritance is organized along paternal descent groups), social behaviors are more likely to be inherited through the paternal line of descent only. It has recently been argued that the rapid spread of Genghis Khan's patrilineal descendants throughout Central Asia was explained by this social selection phenomenon ##REF##12592608##[30]##. The correlation of fertility through the patriline has also been described in patrilineal tribes in South America ##REF##5473413##[31]##. By contrast, in bilineal societies such as the agriculturalists of Central Asia, social behaviors that influence reproductive success are more likely to be transmitted by both sexes. Furthermore, differences of cultural transmission of fitness between hunter-gatherers and agriculturalists have already been reported ##REF##16933997##[32]##. Interestingly, a slightly higher matrilineal intergenerational correlation in offspring number has been observed in the Icelandic population, which suggests that in some populations, reproductive behaviors can be maternally-inherited ##REF##12721957##[33]##.</p> <p>
<italic>Polygyny</italic>, in which the husband may have multiple wives, has often been invoked as a factor that could reduce the effective number of men ##REF##8751879##[4]##, ##REF##12532283##[7]##, ##REF##15190128##[15]##, ##REF##12962309##[23]##–##UREF##1##[25]##. While we could not find any evidence of polygyny in present-day Central Asian populations, this custom was traditionally practiced in the nomadic herder Kazak populations, although limited to the top 10 percent of men from the highest social rank ##REF##10364534##[5]##,##UREF##2##[34]##. Hence, even though we lack ethnological data to determine to what extent herders are or were practicing polygyny in a recent past, the practice of polygyny among herders in Central Asia might have influenced (at least partially) the observed differences in men and women effective numbers.</p> <p>
<italic>Recurrent bottlenecks in men</italic> due to a higher pre-reproductive mortality could also severely reduce the effective numbers of men. From the study of several groups in West Papua and Papua New Guinea ##REF##12532283##[7]##,##UREF##3##[35]##, it appears that warfare may indeed lead to the quasi-extinction of adult men in some communities, while the mass killing of adult women is far more rarely reported. However, this differential mortality could also be balanced by potentially high death rates of women during childbirth. In any case, a differential mortality is equally likely to arise in herder and agriculturalist populations. It may therefore not be relevant in explaining why we detect higher effective numbers of women (as compared to men) in patrilineal herders and not in bilineal agriculturalists.</p> <p>Since our approach implicitly assumes equal male and female generation time, the observed higher effective number of women, relatively to that of men, could result from a <italic>shorter generation time for women</italic>, due to the tendency of women to reproduce earlier in life than men and the ability of men to reproduce at a later age than women. This has indeed been described in a number of populations with different lifestyles, from complete genealogical records or mean-age-at-first-marriage databases ##REF##12721957##[33]##,##REF##10677323##[36]##,##REF##15795887##[37]##. It has even been proposed to be a nearly universal trait in humans, although its magnitude varies across regions and cultures ##REF##15795887##[37]##. Tang et al. ##REF##12019257##[38]## suggested that accounting for longer generation time in males could minimize the difference between maternal and paternal demography. However, the differences in sex-specific generation times that have been reported (e.g., 28 years for the matrilines and 31 years for the patrilines in Iceland ##REF##12721957##[33]##, 29 years for the matrilines and 35 years for the patrilines in Quebec ##REF##10677323##[36]##) are unlikely to explain the observed differences in male and female effective numbers ##REF##15317874##[24]##.</p> <title>Limits of the Approach</title> <p>There might also be non-biological explanations of our results, however, as they are based on the simplifying assumptions of Wright's infinite island model of population structure ##REF##10098262##[39]##. This model assumes (<italic>i</italic>) that there is no selection and that mutation is negligible, (<italic>ii</italic>) that each population has the same size, and sends and receives a constant fraction of its individuals to or from a common migrant pool each generation (so that geographical structure is absent), and (<italic>iii</italic>) that equilibrium is reached between migration, mutation and drift. On the first point, we did not find any evidence of selection, for any marker, based on Beaumont and Nichols' method ##UREF##4##[40]## for detecting selected markers from the analysis of the null distribution generated by a coalescent-based simulation model (data not shown). As for the second point, we tested for the significance of the correlation between the pairwise <italic>F</italic>
<sub>ST</sub>/(1−<italic>F</italic>
<sub>ST</sub>) estimates and the natural logarithm of their geographical distances ##REF##9093870##[41]##. We found no evidence for isolation by distance, either for X-linked markers (<italic>p</italic> = 0.47 for agriculturalists, <italic>p</italic> = 0.24 for herders), or for autosomal markers (<italic>p</italic> = 0.92 for agriculturalists, <italic>p</italic> = 0.45 for herders). As for the third point, the X-to-autosomes (X/A) effective size ratio can significantly deviate from the expected three-quarters (assuming equal effective numbers of men and women) following a bottleneck or an expansion ##REF##17971168##[42]##. This is because X-linked genes have a smaller effective size, and hence reach equilibrium more rapidly. After a reduction of population size, the X/A diversity ratio is lower than expected, while after an expansion, the diversity of X-linked genes recovers faster than on the autosomes, and the X/A diversity ratio is then closer to unity. In the latter case,  would be reduced and could then tend towards . However, neither reduction nor expansion should lead to , as we found in herder populations of Central Asia. Therefore, we do not expect the limits of Wright's island model to undermine our approach.</p> <title>Evaluation by Means of Stochastic Simulations</title> <p>We aimed to investigate to what extent the approach proposed here is able to detect differences in male and female effective numbers. To do this, we performed coalescent simulations in a finite island model, for a wide range of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values. The simulation parameters were set to match those of our dataset: 11 sampled demes, 30 males genotyped at 27 autosomal and 9 X-linked markers per deme (for further details concerning the simulations, see the <xref ref-type="sec" rid="s3">Methods</xref> section). We used 1421 sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, covering the whole parameter space (represented as white dots in ##FIG##3##Figure 4B##). For each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) parameter values, we simulated 100 independent datasets. For each dataset, we calculated the estimates of  at all loci, and we calculated the <italic>p</italic>-value for a one-sided Wilcoxon sum rank test for the list of  estimates . Hence, for each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) parameter values, we could calculate the proportion of significant tests at the <italic>α</italic> = 0.05 level, among the 100 independent datasets. ##FIG##3##Figure 4## shows the distribution of the percentage of significant tests in the (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) parameter space. Theory predicts that in the upper-right triangle where , we should have . One can see from ##FIG##3##Figure 4## that, given the simulation parameters used, the method is conservative: the proportion of significant tests at the <italic>α</italic> = 0.05 level is null outside of the upper-right triangle. However, we find a fairly large proportion of significant tests for large <italic>N</italic>
<sub>f</sub>/<italic>N</italic> and <italic>m</italic>
<sub>f</sub>/<italic>m</italic> ratios which indicates (<italic>i</italic>) that the method presented here has the potential to detect differences in male and female effective numbers, but (<italic>ii</italic>) that only strong differences might be detected, for similarly sized datasets as the one considered here.</p> <title>Robustness to the Sampling Scheme</title> <p>We also aimed to investigate whether the results obtained here were robust to our sampling scheme, and that our results were not biased by the inclusion of particular populations. To do this, we re-analyzed both the bilineal agriculturalists and the patrilineal herders datasets, removing one population at a time in each group. For each of these jackknifed datasets, we calculated the <italic>p</italic>-value of a one-sided Wilcoxon sum rank test , as done on the full datasets. The results are given in ##TAB##4##Table 5##. We found no significant test for any of the bilineal agriculturalist groupings (<italic>p</italic>&gt;0.109), which supports the idea that, in those populations, both the migration rate and the number of reproductive individuals can be equal for both sexes. In patrilineal herders, the tests were significant at the <italic>α</italic> = 0.05 level for 8 out of 11 population groupings. For the 3 other groupings, the <italic>p</italic>-values were 0.068, 0.078 and 0.073 (see ##TAB##4##Table 5##). Overall, the ratio of  multi-locus estimates ranged from 1.7 to 3.5 in patrilineal herders (and from 0.9 to 1.2 in bilineal agriculturalists). Although in some particular groupings of patrilineal herder populations, the difference in the distributions of  may not be strong enough to be significant, we can clearly distinguish the pattern of differentiation for autosomal and X-linked markers in patrilineal and bilineal groups. Results from coalescent simulations (see above) suggest that this lack of statistical power might be expected for  ratios close to unity. Indeed, we found that the tests were more likely to be significant for fairly large <italic>N</italic>
<sub>f</sub>/<italic>N</italic> and <italic>m</italic>
<sub>f</sub>/<italic>m</italic> ratios (the upper-right red region in ##FIG##3##Figure 4##) which would correspond to  ratios much greater than one.</p> <title>Comparison with Uniparentally-Inherited Markers</title> <p>Importantly, our results on X-linked and autosomal markers are consistent with those obtained from NRY and mtDNA (see ##FIG##2##Figures 3B–3D##): in these figures, the dashed line gives all the sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values that are compatible with the observed  estimates. These are the sets of values that satisfy  for the bilineal populations, and  for the patrilineal populations, since we inferred <italic>N</italic>
<sub>f</sub>
<italic>m</italic>
<sub>f</sub>/<italic>N</italic>
<sub>m</sub>
<italic>m</italic>
<sub>m</sub>≈2.1 and <italic>N</italic>
<sub>f</sub>
<italic>m</italic>
<sub>f</sub>/<italic>N</italic>
<sub>m</sub>
<italic>m</italic>
<sub>m</sub>≈21.6, respectively, for the two groups. For the bilineal agriculturalists (##FIG##2##Figure 3D##), the set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values inferred from the  estimates fall within the range that was not rejected, given our data on X-linked and autosomal markers. For the patrilineal herders (##FIG##2##Figure 3B##), the overlap is only partial: from the NRY and mtDNA data only, low <italic>N</italic>
<sub>f</sub>/<italic>N</italic> ratios associated with high <italic>m</italic>
<sub>f</sub>/<italic>m</italic> ratios are as likely as high <italic>N</italic>
<sub>f</sub>/<italic>N</italic> ratios associated with low <italic>m</italic>
<sub>f</sub>/<italic>m</italic> ratios. Yet, it is clear from this figure that a large set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values inferred from the single-locus estimates  can be rejected, given the observed differentiation on X-linked and autosomal markers. All genetic systems (mtDNA, NRY, X-linked and autosomal markers) converge toward the notion that patrilineal herders, in contrast to bilineal agriculturalists, have a strong sex-specific genetic structure. Yet, the information brought by X-linked and autosomal markers is substantial, since we show that this is likely due to both higher migration rates and larger effective numbers for women than for men.</p> <title>Comparison with Other Studies</title> <p>Our results, based on the X chromosome and the autosomes, also confirm previous analyses based on the mtDNA and the NRY, showing that men are genetically more structured than women in other patrilocal populations ##REF##9806547##[3]##–##REF##15996169##[10]##, ##REF##11528385##[14]##–##REF##16916941##[17]## (see also ##TAB##0##Table 1##). A handful of studies have also shown a reduced effective number of men compared to that of women, based on coalescent methods ##REF##12962309##[23]##,##REF##15317874##[24]##, but none have considered the influence of social organization on this dissimilarity (see ##TAB##0##Table 1##).</p> <p>In some respects, our results contrast with those of Wilder and Hammer ##UREF##1##[25]##, who studied sex-specific population genetic structure among the Baining of New Britain, using mtDNA, NRY, and X-linked markers. Interestingly, they found that <italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub>, but <italic>m</italic>
<sub>f</sub>&lt;<italic>m</italic>
<sub>m</sub>, and claimed that a similar result, although left unexplored by the authors, was to be found in a recent study by Hamilton et al. ##REF##15894624##[16]##. This raises the interesting point that sex-specific proportions of migrants (<italic>m</italic>) are likely to be shaped by factors that may only partially overlap with those that affect the sex-specific effective numbers (<italic>N</italic>). Further studies of human populations with contrasted social organizations, as well as further theoretical developments, are needed to appreciate this point.</p> <p>In order to ask to what extent our results generalize to other human populations, we investigated sex-specific patterns in the 51 worldwide populations represented in the HGDP-CEPH Human Genome Diversity Cell Line Panel dataset ##REF##12493913##[43]##, for which the data on the differentiation of 784 autosomal microsatellites and 36 X-linked microsatellites are available (data not shown). By doing this, we found a larger differentiation for X-linked than for autosomal markers . Therefore, we confirmed Ramachandran et al.'s ##REF##15601537##[20]## result that no major differences in demographic parameters between males and females are required to explain the X-chromosomal and autosomal results in this worldwide sample. Ramachandran et al.'s approach ##REF##15601537##[20]## is based upon a pure divergence model from a single ancestral population, which is very different from the migration-drift equilibrium model considered here. In real populations, however, genetic differentiation almost certainly arises both through divergence and limited dispersal, which places these two models at two ends of a continuum. Yet, importantly, if we apply Ramachandran et al.'s ##REF##15601537##[20]## model to the Central Asian data, our conclusions are left unchanged. In their model, the differentiation among populations is , where <italic>t</italic> is the time since divergence from an ancestral population and <italic>N</italic>
<sub>e</sub> the effective size of the populations (see, e.g., ##REF##17246175##[44]##). Hence, we get  for autosomal and X-linked markers, respectively. Therefore, our observation that  implies that , which requires that <italic>N</italic>
<sub>f</sub>&gt;7<italic>N</italic>
<sub>m</sub> since  (see, e.g., ##UREF##5##[45]##). In this case, the female fraction of effective number is larger than that of males, which is consistent with our findings in a model with migration.</p> <p>The HGDP-CEPH dataset does not provide any detailed ethnic information for the sampled groups, and we can therefore not distinguish populations with different lifestyles. However, at a more local scale in Pakistan, we were able to analyze a subset of 5 populations (Brahui, Balochi, Makrani, Sindhi and Pathan), which are presumed to be patrilineal ##UREF##6##[46]##. For this subset, we found a higher differentiation for autosomal  than for X-linked markers , although non-significantly (<italic>p</italic> = 0.12). This result seems to suggest that other patrilineal populations may behave like the Central Asian sample presented here. Therefore, because the geographical clustering of populations with potentially different lifestyles may minimize the differences in sex-specific demography at a global scale ##REF##16617372##[21]##,##REF##11420360##[22]##, and/or because the global structure may reflect ancient (pre-agricultural) marital residence patterns with less pronounced patrilocality ##REF##16479583##[12]##, we emphasize the point that large-scale studies may not be relevant to detect sex-specific patterns, which supports a claim made by many authors.</p> <title>Conclusion</title> <p>In conclusion, we have shown here that the joint analysis of autosomal and X-linked polymorphic markers provides an efficient tool to infer sex-specific demography and history in human populations, as suggested recently ##REF##16479583##[12]##,##REF##17208169##[47]##. This new multilocus approach is, to our knowledge, the first attempt to combine the information contained in mtDNA, NRY, X-linked and autosomal markers (see ##TAB##0##Table 1##), which allowed us to test for the robustness of a sex-specific genetic structure at a local scale. Unraveling the respective influence of migration and drift upon neutral genetic structure is a long-standing quest in population genetics ##REF##17402974##[48]##,##REF##11472529##[49]##. Here, our analysis allowed us to show that differences in sex-specific migration rates may not be the only cause of contrasted male and female differentiation in humans and that, contrary to the conclusion of a number of studies (see ##TAB##0##Table 1##), differences in effective numbers may also play an important role. Indeed, we have demonstrated that sex-specific differences in population structure in patrilineal herders may be the consequence of both higher female effective numbers and female effective dispersal. Our results also illustrate the importance of analyzing human populations at a local scale, rather than global or even continental scale ##REF##17067791##[2]##,##REF##15378061##[19]##,##REF##16617372##[21]##. The originality of our approach lies in the comparison of identified ethnic groups that differ in well-known social structures and lifestyles. In that respect, our study is among the very few which compare patrilineal vs. bilineal or matrilineal groups (see ##TAB##0##Table 1##), and we believe that it contributes to the growing body of evidence showing that social organization and lifestyle have a strong impact on the distribution of genetic variation in human populations. Moreover, our approach could also be applied on a wide range of animal species with contrasted social organizations. Therefore, we expect our results to stimulate research on the comparison of X-linked and autosomal data to disentangle sex-specific demography.</p>
<title>Results/Discussion</title> <title>Uniparentally-Inherited Markers</title> <p>We sampled 780 healthy adult men from 10 populations of bilineal agriculturalists and 11 populations of patrilineal herders from West Uzbekistan to East Kyrgyzstan, representing 5 ethnic groups (Tajiks, Kyrgyz, Karakalpaks, Kazaks, and Turkmen) (see ##FIG##0##Figure 1## and ##TAB##1##Table 2##). We genotyped all bilineal populations, and 8 out of 11 patrilineal populations at the HVS-I locus of mtDNA, and at 11 microsatellite markers on the NRY (for more details on the markers used, see ##TAB##2##Table 3##). The overall genetic differentiation was higher for NRY, as compared to mtDNA, both among the 10 bilineal agriculturalist populations , and among the subset of 8 patrilineal herder populations . Assuming an island model of population structure, this implies that female migration rate (<italic>m</italic>
<sub>f</sub>), and/or the effective number of females (<italic>N</italic>
<sub>f</sub>), is higher than of the corresponding parameters for males (<italic>m</italic>
<sub>m</sub> and <italic>N</italic>
<sub>m</sub>). These results also suggest that the differences in sex-specific genetic structure are much more pronounced in the patrilineal herders than in the bilineal agriculturalists. From the above <italic>F</italic>
<sub>ST</sub> estimates, we obtained the female-to-male ratio of the effective number of migrants per generation (see the <xref ref-type="sec" rid="s3">Methods</xref> section for details): <italic>N</italic>
<sub>f</sub>
<italic>m</italic>
<sub>f</sub>/<italic>N</italic>
<sub>m</sub>
<italic>m</italic>
<sub>m</sub>≈2.1 for bilineal populations and <italic>N</italic>
<sub>f</sub>
<italic>m</italic>
<sub>f</sub>/<italic>N</italic>
<sub>m</sub>
<italic>m</italic>
<sub>m</sub>≈21.6 for patrilineal populations. The ratio in patrilineal populations is thus one order of magnitude higher than in bilineal populations. However, since each of these markers is a single genetic locus, we cannot test for the robustness of the sex-specific genetic structure on these markers. We therefore examined the amount of information contained in multi-locus data on autosomal and X-linked markers, both of which average over male and female histories.</p> <title>A New Multi-Locus Approach</title> <p>In the infinite island model of population structure with two classes of individuals (males and females), we obtained the following expressions of <italic>F</italic>
<sub>ST</sub> (see the <xref ref-type="sec" rid="s3">Methods</xref> section for details):for autosomal genes, andfor X-linked genes. A special case of interest occurs when , i.e. when the differentiation of X-linked genes exactly equals that of autosomal genes. Combining eqs (1) and (2), we find that this occurs for , with <italic>N</italic> = <italic>N</italic>
<sub>f</sub>+<italic>N</italic>
<sub>m</sub> and <italic>m</italic> = <italic>m</italic>
<sub>f</sub>+<italic>m</italic>
<sub>m</sub>. Furthermore, as shown in ##FIG##1##Figure 2##, if we observe a lower genetic differentiation of autosomal markers, as compared to X-linked markers (blue zone in ##FIG##1##Figure 2##), this suggests that . This may happen, e.g., for <italic>N</italic>
<sub>f</sub> = <italic>N</italic>
<sub>m</sub> and <italic>m</italic>
<sub>f</sub> = <italic>m</italic>
<sub>m</sub>, i.e. for equal effective numbers of males and females and unbiased dispersal. But if autosomal markers are more differentiated than X-linked markers (, see the red upper-right triangle in ##FIG##1##Figure 2##), this implies that . In this case, since <italic>m</italic>
<sub>f</sub>/<italic>m</italic> and <italic>N</italic>
<sub>f</sub>/<italic>N</italic> are ratios varying between 0 and 1, the effective number of females must be higher than that of males (<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub>), and the female migration rate must be higher than half the male migration rate (<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub>/2). Hence, a prediction from this model is that when , the effective number of females is higher than that of males, whatever the pattern of sex-specific dispersal. This suggests that it is indeed possible to test for differences in effective numbers between males and females from the joint analysis of autosomal and X-linked data. We note however that when , we cannot conclude on the relative male and female effective numbers and migration rates.</p> <p>We tested the above prediction in the 10 bilineal agriculturalist populations and 11 patrilineal herder populations sampled in Central Asia by comparing the genetic structure estimated from 27 unlinked polymorphic autosomal microsatellite markers (<italic>AR</italic> = 16.2, <italic>H</italic>
<sub>e</sub> = 0.803 on average) to that from 9 unlinked polymorphic X-linked microsatellite markers (<italic>AR</italic> = 12.6, <italic>H</italic>
<sub>e</sub> = 0.752 on average) (for more details on the markers used, see ##TAB##3##Table 4##). Overall heterozygosity was not significantly different between X-linked and autosomal markers, neither in the pooled sample (two-tailed Wilcoxon sum rank test; <italic>p</italic> = 0.09), nor in the bilineal agriculturalists (<italic>p</italic> = 0.13) or the patrilineal herders (<italic>p</italic> = 0.12). The overall population structure was significantly higher for autosomal as compared to X-linked markers among patrilineal herders:  (one-tailed Wilcoxon sum rank test; ; <italic>p</italic> = 0.02). Among bilineal agriculturalists, the result was not significant:  (<italic>p</italic> = 0.36). From these results, and following our model predictions, we conclude that in patrilineal herders (where ), the effective number of females is higher than that of males. This conclusion does not hold for the bilineal agriculturalists.</p> <p>From our model, it is possible to get more precise indications on the sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values that are compatible with our data. Rearranging eqs (1–2), we get:i.e.:
</p> <p>For any given set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, we can therefore calculate from eq. (4) the expected value of  for each  estimate in the dataset. We can then test the null hypothesis  by comparing the distribution of observed and expected  values. If the hypothesis can be rejected at the <italic>α</italic> = 0.05 level, then the corresponding set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values can also be rejected. Following Ramachandran et al. ##REF##15601537##[20]##, we varied the values of the ratios <italic>N</italic>
<sub>f</sub>/<italic>N</italic> and <italic>m</italic>
<sub>f</sub>/<italic>m</italic> (respectively, the female fraction of effective number, and the female fraction of the total migration rate) from 0 to 1, with an interval of 0.01 between consecutive values. For each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, we applied the transformation in eq. (4) to each of the 27 locus-specific  values observed. Thus, for each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, we obtained 27 expected values of , given our data. These expected values of  were then compared to the 9 observed locus-specific  in our dataset, and we calculated the <italic>p</italic>-value for a two-sided Wilcoxon sum rank test between the list of 27 expected  values and the 9  observed in the dataset. The results are depicted in ##FIG##2##Figure 3##. Significant <italic>p</italic>-values (<italic>p</italic>≤0.05) correspond to a significant difference between the observed and expected values, thus to sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values that are rejected, given our data (see the blue region in ##FIG##2##Figure 3##). Conversely, non-significant <italic>p</italic>-values (<italic>p</italic>&gt;0.05) correspond to sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values that cannot be rejected (see the red region in ##FIG##2##Figure 3##).</p> <p>For the patrilineal herder populations (##FIG##2##Figures 3A–3B##), most sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values are rejected, except those corresponding to larger effective numbers for females (from ##FIG##2##Figures 3A–3B##: <italic>N</italic>
<sub>f</sub>/<italic>N</italic>&gt;0.55, i.e. <italic>N</italic>
<sub>f</sub>&gt;1.27<italic>N</italic>
<sub>m</sub>) and <italic>m</italic>
<sub>f</sub>&gt;0.67<italic>m</italic>
<sub>m</sub>. Because the multi-locus estimate of  is significantly higher than the estimate of , we expected to find such patterns of non-significant values (see ##FIG##1##Figure 2##). For the bilineal agriculturalist populations, we could not reject the hypothesis that the effective numbers and migration rates are equal across males and females or even lower in females (see ##FIG##2##Figures 3C–3D##). This is also reflected by the fact that the estimates of  were not significantly higher than the estimates of  in those populations.</p> <p>Finally, we have shown that the effective number of women is higher than that of men among patrilineal herders, but not necessarily among bilineal agriculturalists. Furthermore, a close inspection of the results depicted in ##FIG##2##Figures 3A and 3B## reveals that, among herders, we reject all the sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values for which <italic>m</italic>
<sub>f</sub>&lt;<italic>m</italic>
<sub>m</sub> at the <italic>α</italic> = 0.10 level. This is not true for agriculturalists. This suggests that the migration rates are also likely to be higher for women than for men in patrilineal populations, as compared to bilineal populations (compare ##FIG##2##Figures 3B and 3D##). Although both groups are patrilocal, such a difference in sex-specific migration patterns might be expected, since patrilineal herders are exogamous (among clans) and bilineal agriculturalists are preferentially endogamous. For example, it was observed that in patrilocal and matrilocal Indian populations, where migrations are strictly confined within endogamous groups, sex-specific patterns were not influenced by post-marital residence ##REF##16617372##[21]##.</p> <title>What Could Explain a Larger Effective Number of Females?</title> <p>While an influence of post-marital residence on the migration rate of women and men has already been widely proposed ##REF##11528385##[14]##–##REF##16916941##[17]## (see also ##TAB##0##Table 1##), the factors that may locally affect the effective number of women, relatively to that of men, are not well recognized. As seen in ##TAB##0##Table 1##, although a number of studies have compared matrilocal and patrilocal populations, few have compared contrasting groups of populations with respect to other factors as, e.g., the tendency for polygyny ##REF##15190128##[15]##. Furthermore, a number of these studies lack ethnological information a priori, concerning social organization, marriage rules, etc., which makes interpretation somewhat difficult (see ##TAB##0##Table 1##). Here, we compared two groups of patrilocal populations with contrasting social organizations, and at least five non-mutually exclusive interpretations for a larger effective number of females can be invoked:</p> <p>
<italic>Social organization</italic>, i.e. the way children are affiliated to their parents, can deeply affect sex-specific genetic variation. In Central Asia, herder populations are organized in patrilineal descent groups (tribes, clans, lineages). This implies that children are systematically affiliated with the descent groups of the father. Chaix et al. ##REF##17208185##[11]## showed that the average number of individuals carrying the same Y chromosome haplotype was much higher in patrilineal herder populations than in bilineal agriculturalist populations (where children are affiliated both to the mother and the father). These “identity cores” would be the direct consequence of the internal dynamics of their patrilineal organization. Indeed, the descent groups are not formed randomly and related men tend to cluster together, e.g. through the recurrent lineal fission of one population into new groups. This particular dynamics increases relatedness among men, and may therefore reduce the effective number of men, as compared to women.</p> <p>Indirectly, the social organization can also deflate the effective number of men through <italic>the transmission of reproductive success</italic>
##REF##15797619##[29]## if this success is culturally transmitted exclusively from fathers to sons. Because herders are patrilineal (so that inheritance is organized along paternal descent groups), social behaviors are more likely to be inherited through the paternal line of descent only. It has recently been argued that the rapid spread of Genghis Khan's patrilineal descendants throughout Central Asia was explained by this social selection phenomenon ##REF##12592608##[30]##. The correlation of fertility through the patriline has also been described in patrilineal tribes in South America ##REF##5473413##[31]##. By contrast, in bilineal societies such as the agriculturalists of Central Asia, social behaviors that influence reproductive success are more likely to be transmitted by both sexes. Furthermore, differences of cultural transmission of fitness between hunter-gatherers and agriculturalists have already been reported ##REF##16933997##[32]##. Interestingly, a slightly higher matrilineal intergenerational correlation in offspring number has been observed in the Icelandic population, which suggests that in some populations, reproductive behaviors can be maternally-inherited ##REF##12721957##[33]##.</p> <p>
<italic>Polygyny</italic>, in which the husband may have multiple wives, has often been invoked as a factor that could reduce the effective number of men ##REF##8751879##[4]##, ##REF##12532283##[7]##, ##REF##15190128##[15]##, ##REF##12962309##[23]##–##UREF##1##[25]##. While we could not find any evidence of polygyny in present-day Central Asian populations, this custom was traditionally practiced in the nomadic herder Kazak populations, although limited to the top 10 percent of men from the highest social rank ##REF##10364534##[5]##,##UREF##2##[34]##. Hence, even though we lack ethnological data to determine to what extent herders are or were practicing polygyny in a recent past, the practice of polygyny among herders in Central Asia might have influenced (at least partially) the observed differences in men and women effective numbers.</p> <p>
<italic>Recurrent bottlenecks in men</italic> due to a higher pre-reproductive mortality could also severely reduce the effective numbers of men. From the study of several groups in West Papua and Papua New Guinea ##REF##12532283##[7]##,##UREF##3##[35]##, it appears that warfare may indeed lead to the quasi-extinction of adult men in some communities, while the mass killing of adult women is far more rarely reported. However, this differential mortality could also be balanced by potentially high death rates of women during childbirth. In any case, a differential mortality is equally likely to arise in herder and agriculturalist populations. It may therefore not be relevant in explaining why we detect higher effective numbers of women (as compared to men) in patrilineal herders and not in bilineal agriculturalists.</p> <p>Since our approach implicitly assumes equal male and female generation time, the observed higher effective number of women, relatively to that of men, could result from a <italic>shorter generation time for women</italic>, due to the tendency of women to reproduce earlier in life than men and the ability of men to reproduce at a later age than women. This has indeed been described in a number of populations with different lifestyles, from complete genealogical records or mean-age-at-first-marriage databases ##REF##12721957##[33]##,##REF##10677323##[36]##,##REF##15795887##[37]##. It has even been proposed to be a nearly universal trait in humans, although its magnitude varies across regions and cultures ##REF##15795887##[37]##. Tang et al. ##REF##12019257##[38]## suggested that accounting for longer generation time in males could minimize the difference between maternal and paternal demography. However, the differences in sex-specific generation times that have been reported (e.g., 28 years for the matrilines and 31 years for the patrilines in Iceland ##REF##12721957##[33]##, 29 years for the matrilines and 35 years for the patrilines in Quebec ##REF##10677323##[36]##) are unlikely to explain the observed differences in male and female effective numbers ##REF##15317874##[24]##.</p> <title>Limits of the Approach</title> <p>There might also be non-biological explanations of our results, however, as they are based on the simplifying assumptions of Wright's infinite island model of population structure ##REF##10098262##[39]##. This model assumes (<italic>i</italic>) that there is no selection and that mutation is negligible, (<italic>ii</italic>) that each population has the same size, and sends and receives a constant fraction of its individuals to or from a common migrant pool each generation (so that geographical structure is absent), and (<italic>iii</italic>) that equilibrium is reached between migration, mutation and drift. On the first point, we did not find any evidence of selection, for any marker, based on Beaumont and Nichols' method ##UREF##4##[40]## for detecting selected markers from the analysis of the null distribution generated by a coalescent-based simulation model (data not shown). As for the second point, we tested for the significance of the correlation between the pairwise <italic>F</italic>
<sub>ST</sub>/(1−<italic>F</italic>
<sub>ST</sub>) estimates and the natural logarithm of their geographical distances ##REF##9093870##[41]##. We found no evidence for isolation by distance, either for X-linked markers (<italic>p</italic> = 0.47 for agriculturalists, <italic>p</italic> = 0.24 for herders), or for autosomal markers (<italic>p</italic> = 0.92 for agriculturalists, <italic>p</italic> = 0.45 for herders). As for the third point, the X-to-autosomes (X/A) effective size ratio can significantly deviate from the expected three-quarters (assuming equal effective numbers of men and women) following a bottleneck or an expansion ##REF##17971168##[42]##. This is because X-linked genes have a smaller effective size, and hence reach equilibrium more rapidly. After a reduction of population size, the X/A diversity ratio is lower than expected, while after an expansion, the diversity of X-linked genes recovers faster than on the autosomes, and the X/A diversity ratio is then closer to unity. In the latter case,  would be reduced and could then tend towards . However, neither reduction nor expansion should lead to , as we found in herder populations of Central Asia. Therefore, we do not expect the limits of Wright's island model to undermine our approach.</p> <title>Evaluation by Means of Stochastic Simulations</title> <p>We aimed to investigate to what extent the approach proposed here is able to detect differences in male and female effective numbers. To do this, we performed coalescent simulations in a finite island model, for a wide range of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values. The simulation parameters were set to match those of our dataset: 11 sampled demes, 30 males genotyped at 27 autosomal and 9 X-linked markers per deme (for further details concerning the simulations, see the <xref ref-type="sec" rid="s3">Methods</xref> section). We used 1421 sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, covering the whole parameter space (represented as white dots in ##FIG##3##Figure 4B##). For each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) parameter values, we simulated 100 independent datasets. For each dataset, we calculated the estimates of  at all loci, and we calculated the <italic>p</italic>-value for a one-sided Wilcoxon sum rank test for the list of  estimates . Hence, for each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) parameter values, we could calculate the proportion of significant tests at the <italic>α</italic> = 0.05 level, among the 100 independent datasets. ##FIG##3##Figure 4## shows the distribution of the percentage of significant tests in the (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) parameter space. Theory predicts that in the upper-right triangle where , we should have . One can see from ##FIG##3##Figure 4## that, given the simulation parameters used, the method is conservative: the proportion of significant tests at the <italic>α</italic> = 0.05 level is null outside of the upper-right triangle. However, we find a fairly large proportion of significant tests for large <italic>N</italic>
<sub>f</sub>/<italic>N</italic> and <italic>m</italic>
<sub>f</sub>/<italic>m</italic> ratios which indicates (<italic>i</italic>) that the method presented here has the potential to detect differences in male and female effective numbers, but (<italic>ii</italic>) that only strong differences might be detected, for similarly sized datasets as the one considered here.</p> <title>Robustness to the Sampling Scheme</title> <p>We also aimed to investigate whether the results obtained here were robust to our sampling scheme, and that our results were not biased by the inclusion of particular populations. To do this, we re-analyzed both the bilineal agriculturalists and the patrilineal herders datasets, removing one population at a time in each group. For each of these jackknifed datasets, we calculated the <italic>p</italic>-value of a one-sided Wilcoxon sum rank test , as done on the full datasets. The results are given in ##TAB##4##Table 5##. We found no significant test for any of the bilineal agriculturalist groupings (<italic>p</italic>&gt;0.109), which supports the idea that, in those populations, both the migration rate and the number of reproductive individuals can be equal for both sexes. In patrilineal herders, the tests were significant at the <italic>α</italic> = 0.05 level for 8 out of 11 population groupings. For the 3 other groupings, the <italic>p</italic>-values were 0.068, 0.078 and 0.073 (see ##TAB##4##Table 5##). Overall, the ratio of  multi-locus estimates ranged from 1.7 to 3.5 in patrilineal herders (and from 0.9 to 1.2 in bilineal agriculturalists). Although in some particular groupings of patrilineal herder populations, the difference in the distributions of  may not be strong enough to be significant, we can clearly distinguish the pattern of differentiation for autosomal and X-linked markers in patrilineal and bilineal groups. Results from coalescent simulations (see above) suggest that this lack of statistical power might be expected for  ratios close to unity. Indeed, we found that the tests were more likely to be significant for fairly large <italic>N</italic>
<sub>f</sub>/<italic>N</italic> and <italic>m</italic>
<sub>f</sub>/<italic>m</italic> ratios (the upper-right red region in ##FIG##3##Figure 4##) which would correspond to  ratios much greater than one.</p> <title>Comparison with Uniparentally-Inherited Markers</title> <p>Importantly, our results on X-linked and autosomal markers are consistent with those obtained from NRY and mtDNA (see ##FIG##2##Figures 3B–3D##): in these figures, the dashed line gives all the sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values that are compatible with the observed  estimates. These are the sets of values that satisfy  for the bilineal populations, and  for the patrilineal populations, since we inferred <italic>N</italic>
<sub>f</sub>
<italic>m</italic>
<sub>f</sub>/<italic>N</italic>
<sub>m</sub>
<italic>m</italic>
<sub>m</sub>≈2.1 and <italic>N</italic>
<sub>f</sub>
<italic>m</italic>
<sub>f</sub>/<italic>N</italic>
<sub>m</sub>
<italic>m</italic>
<sub>m</sub>≈21.6, respectively, for the two groups. For the bilineal agriculturalists (##FIG##2##Figure 3D##), the set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values inferred from the  estimates fall within the range that was not rejected, given our data on X-linked and autosomal markers. For the patrilineal herders (##FIG##2##Figure 3B##), the overlap is only partial: from the NRY and mtDNA data only, low <italic>N</italic>
<sub>f</sub>/<italic>N</italic> ratios associated with high <italic>m</italic>
<sub>f</sub>/<italic>m</italic> ratios are as likely as high <italic>N</italic>
<sub>f</sub>/<italic>N</italic> ratios associated with low <italic>m</italic>
<sub>f</sub>/<italic>m</italic> ratios. Yet, it is clear from this figure that a large set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values inferred from the single-locus estimates  can be rejected, given the observed differentiation on X-linked and autosomal markers. All genetic systems (mtDNA, NRY, X-linked and autosomal markers) converge toward the notion that patrilineal herders, in contrast to bilineal agriculturalists, have a strong sex-specific genetic structure. Yet, the information brought by X-linked and autosomal markers is substantial, since we show that this is likely due to both higher migration rates and larger effective numbers for women than for men.</p> <title>Comparison with Other Studies</title> <p>Our results, based on the X chromosome and the autosomes, also confirm previous analyses based on the mtDNA and the NRY, showing that men are genetically more structured than women in other patrilocal populations ##REF##9806547##[3]##–##REF##15996169##[10]##, ##REF##11528385##[14]##–##REF##16916941##[17]## (see also ##TAB##0##Table 1##). A handful of studies have also shown a reduced effective number of men compared to that of women, based on coalescent methods ##REF##12962309##[23]##,##REF##15317874##[24]##, but none have considered the influence of social organization on this dissimilarity (see ##TAB##0##Table 1##).</p> <p>In some respects, our results contrast with those of Wilder and Hammer ##UREF##1##[25]##, who studied sex-specific population genetic structure among the Baining of New Britain, using mtDNA, NRY, and X-linked markers. Interestingly, they found that <italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub>, but <italic>m</italic>
<sub>f</sub>&lt;<italic>m</italic>
<sub>m</sub>, and claimed that a similar result, although left unexplored by the authors, was to be found in a recent study by Hamilton et al. ##REF##15894624##[16]##. This raises the interesting point that sex-specific proportions of migrants (<italic>m</italic>) are likely to be shaped by factors that may only partially overlap with those that affect the sex-specific effective numbers (<italic>N</italic>). Further studies of human populations with contrasted social organizations, as well as further theoretical developments, are needed to appreciate this point.</p> <p>In order to ask to what extent our results generalize to other human populations, we investigated sex-specific patterns in the 51 worldwide populations represented in the HGDP-CEPH Human Genome Diversity Cell Line Panel dataset ##REF##12493913##[43]##, for which the data on the differentiation of 784 autosomal microsatellites and 36 X-linked microsatellites are available (data not shown). By doing this, we found a larger differentiation for X-linked than for autosomal markers . Therefore, we confirmed Ramachandran et al.'s ##REF##15601537##[20]## result that no major differences in demographic parameters between males and females are required to explain the X-chromosomal and autosomal results in this worldwide sample. Ramachandran et al.'s approach ##REF##15601537##[20]## is based upon a pure divergence model from a single ancestral population, which is very different from the migration-drift equilibrium model considered here. In real populations, however, genetic differentiation almost certainly arises both through divergence and limited dispersal, which places these two models at two ends of a continuum. Yet, importantly, if we apply Ramachandran et al.'s ##REF##15601537##[20]## model to the Central Asian data, our conclusions are left unchanged. In their model, the differentiation among populations is , where <italic>t</italic> is the time since divergence from an ancestral population and <italic>N</italic>
<sub>e</sub> the effective size of the populations (see, e.g., ##REF##17246175##[44]##). Hence, we get  for autosomal and X-linked markers, respectively. Therefore, our observation that  implies that , which requires that <italic>N</italic>
<sub>f</sub>&gt;7<italic>N</italic>
<sub>m</sub> since  (see, e.g., ##UREF##5##[45]##). In this case, the female fraction of effective number is larger than that of males, which is consistent with our findings in a model with migration.</p> <p>The HGDP-CEPH dataset does not provide any detailed ethnic information for the sampled groups, and we can therefore not distinguish populations with different lifestyles. However, at a more local scale in Pakistan, we were able to analyze a subset of 5 populations (Brahui, Balochi, Makrani, Sindhi and Pathan), which are presumed to be patrilineal ##UREF##6##[46]##. For this subset, we found a higher differentiation for autosomal  than for X-linked markers , although non-significantly (<italic>p</italic> = 0.12). This result seems to suggest that other patrilineal populations may behave like the Central Asian sample presented here. Therefore, because the geographical clustering of populations with potentially different lifestyles may minimize the differences in sex-specific demography at a global scale ##REF##16617372##[21]##,##REF##11420360##[22]##, and/or because the global structure may reflect ancient (pre-agricultural) marital residence patterns with less pronounced patrilocality ##REF##16479583##[12]##, we emphasize the point that large-scale studies may not be relevant to detect sex-specific patterns, which supports a claim made by many authors.</p> <title>Conclusion</title> <p>In conclusion, we have shown here that the joint analysis of autosomal and X-linked polymorphic markers provides an efficient tool to infer sex-specific demography and history in human populations, as suggested recently ##REF##16479583##[12]##,##REF##17208169##[47]##. This new multilocus approach is, to our knowledge, the first attempt to combine the information contained in mtDNA, NRY, X-linked and autosomal markers (see ##TAB##0##Table 1##), which allowed us to test for the robustness of a sex-specific genetic structure at a local scale. Unraveling the respective influence of migration and drift upon neutral genetic structure is a long-standing quest in population genetics ##REF##17402974##[48]##,##REF##11472529##[49]##. Here, our analysis allowed us to show that differences in sex-specific migration rates may not be the only cause of contrasted male and female differentiation in humans and that, contrary to the conclusion of a number of studies (see ##TAB##0##Table 1##), differences in effective numbers may also play an important role. Indeed, we have demonstrated that sex-specific differences in population structure in patrilineal herders may be the consequence of both higher female effective numbers and female effective dispersal. Our results also illustrate the importance of analyzing human populations at a local scale, rather than global or even continental scale ##REF##17067791##[2]##,##REF##15378061##[19]##,##REF##16617372##[21]##. The originality of our approach lies in the comparison of identified ethnic groups that differ in well-known social structures and lifestyles. In that respect, our study is among the very few which compare patrilineal vs. bilineal or matrilineal groups (see ##TAB##0##Table 1##), and we believe that it contributes to the growing body of evidence showing that social organization and lifestyle have a strong impact on the distribution of genetic variation in human populations. Moreover, our approach could also be applied on a wide range of animal species with contrasted social organizations. Therefore, we expect our results to stimulate research on the comparison of X-linked and autosomal data to disentangle sex-specific demography.</p>
<title>Conclusion</title> <p>In conclusion, we have shown here that the joint analysis of autosomal and X-linked polymorphic markers provides an efficient tool to infer sex-specific demography and history in human populations, as suggested recently ##REF##16479583##[12]##,##REF##17208169##[47]##. This new multilocus approach is, to our knowledge, the first attempt to combine the information contained in mtDNA, NRY, X-linked and autosomal markers (see ##TAB##0##Table 1##), which allowed us to test for the robustness of a sex-specific genetic structure at a local scale. Unraveling the respective influence of migration and drift upon neutral genetic structure is a long-standing quest in population genetics ##REF##17402974##[48]##,##REF##11472529##[49]##. Here, our analysis allowed us to show that differences in sex-specific migration rates may not be the only cause of contrasted male and female differentiation in humans and that, contrary to the conclusion of a number of studies (see ##TAB##0##Table 1##), differences in effective numbers may also play an important role. Indeed, we have demonstrated that sex-specific differences in population structure in patrilineal herders may be the consequence of both higher female effective numbers and female effective dispersal. Our results also illustrate the importance of analyzing human populations at a local scale, rather than global or even continental scale ##REF##17067791##[2]##,##REF##15378061##[19]##,##REF##16617372##[21]##. The originality of our approach lies in the comparison of identified ethnic groups that differ in well-known social structures and lifestyles. In that respect, our study is among the very few which compare patrilineal vs. bilineal or matrilineal groups (see ##TAB##0##Table 1##), and we believe that it contributes to the growing body of evidence showing that social organization and lifestyle have a strong impact on the distribution of genetic variation in human populations. Moreover, our approach could also be applied on a wide range of animal species with contrasted social organizations. Therefore, we expect our results to stimulate research on the comparison of X-linked and autosomal data to disentangle sex-specific demography.</p>
<p><bold>¤:</bold> Current address: Unidad de Biología Evolutiva, Departamento de Ciencias Experimentales y de la Salud, Universidad Pompeu Fabra, Barcelona, Spain</p> <p>Conceived and designed the experiments: EH RV. Performed the experiments: LS BMC LQM PB MG. Analyzed the data: LS RV. Contributed reagents/materials/analysis tools: BMC TH AA FN MJ EH. Wrote the paper: LS RV. Collected the samples: LS, BMC, EH.</p> <p>In the last two decades, mitochondrial DNA (mtDNA) and the non-recombining portion of the Y chromosome (NRY) have been extensively used in order to measure the maternally and paternally inherited genetic structure of human populations, and to infer sex-specific demography and history. Most studies converge towards the notion that among populations, women are genetically less structured than men. This has been mainly explained by a higher migration rate of women, due to patrilocality, a tendency for men to stay in their birthplace while women move to their husband's house. Yet, since population differentiation depends upon the product of the effective number of individuals within each deme and the migration rate among demes, differences in male and female effective numbers and sex-biased dispersal have confounding effects on the comparison of genetic structure as measured by uniparentally inherited markers. In this study, we develop a new multi-locus approach to analyze jointly autosomal and X-linked markers in order to aid the understanding of sex-specific contributions to population differentiation. We show that in patrilineal herder groups of Central Asia, in contrast to bilineal agriculturalists, the effective number of women is higher than that of men. We interpret this result, which could not be obtained by the analysis of mtDNA and NRY alone, as the consequence of the social organization of patrilineal populations, in which genetically related men (but not women) tend to cluster together. This study suggests that differences in sex-specific migration rates may not be the only cause of contrasting male and female differentiation in humans, and that differences in effective numbers do matter.</p> <title>Author Summary</title> <p>Human evolutionary history has been investigated mainly through the prism of genetic variation of the Y chromosome and mitochondrial DNA. These two uniparentally inherited markers reflect the demographic history of males and females, respectively. Their contrasting patterns of genetic differentiation reveal that women are more mobile than men among populations, which might be due to specific marriage rules. However, these two markers provide only a limited understanding of the underlying demographic processes. To obtain an independent picture of sex-specific demography, we developed a new multi-locus approach based on the analysis of markers from the autosomal and X-chromosomal compartments. We applied our method to 21 human populations sampled in Central Asia, with contrasting social organizations and lifestyles. We found that, in patrilineal populations, not only the migration rate but also the number of reproductive individuals is likely to be higher for women. This result does not hold for bilineal populations, for which both the migration rate and the number of reproductive individuals can be equal for both sexes. The social organization of patrilineal populations is the likely cause of this pattern. This study suggests that differences in sex-specific migration rates may not be the only cause of contrasting male and female differentiation in humans, and that differences in effective numbers do matter.</p>
<p>We thank all the people who volunteered to participate in this study, or who helped us in the field. We are grateful to Sylvain Théry for valuable help in handling geographic data, to Hélène Fréville and Nicolas Perrin for helpful comments on previous versions of this manuscript, as well as to three anonymous reviewers for insightful and constructive comments. We acknowledge the “Service de Systématique Moléculaire” (SSM) at the Museum National d'Histoire Naturelle (MNHN) and the Biological Resource Center of the Foundation Jean Dausset-CEPH for genotyping facilities. Part of this work was carried out by using the resources of the Computational Biology Service Unit from the Museum National d'Histoire Naturelle (MNHN) which was partially funded by Saint Gobain.</p>
<fig id="pgen-1000200-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000200.g001</object-id><label>Figure 1</label><caption><title>Geographic map of the sampled area, with the 21 populations studied.</title><p>Bilineal agriculturalist populations are in blue (Tajiks); Patrilineal herders with a semi-nomadic lifestyle are in red (Kazaks, Karakalpaks, Kyrgyz and Turkmen).</p></caption></fig> <fig id="pgen-1000200-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000200.g002</object-id><label>Figure 2</label><caption><title>Diagram representing the relative values of expected genetic differentiation for autosomal markers  and for X-linked markers .</title><p>In the red upper right triangle, the <italic>F</italic>
<sub>ST</sub> estimates for autosomal markers are higher than for X-linked markers. In this case, <italic>N</italic>
<sub>f</sub>/<italic>N</italic> is necessarily larger than 0.5. In the blue region of the figure, the <italic>F</italic>
<sub>ST</sub> estimates for autosomal markers are lower than for X-linked markers. The white plain line, at which , represents the set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values where the autosomal and X-linked <italic>F</italic>
<sub>ST</sub> estimates are equal. In this case , if <italic>N</italic>
<sub>f</sub> = <italic>N</italic>
<sub>m</sub>, then the lower effective size of X-linked markers (which would be three-quarters that of autosomal markers) can only be balanced by a complete female-bias in dispersal (<italic>m</italic>
<sub>f</sub>/<italic>m</italic> = 1). Conversely, if <italic>m</italic>
<sub>f</sub> = <italic>m</italic>
<sub>m</sub>, the large female fraction of effective numbers compensates exactly the low effective size of X-linked markers only for <italic>N</italic>
<sub>f</sub> = 7<italic>N</italic>
<sub>m</sub>. Last, if <italic>m</italic>
<sub>f</sub> = <italic>m</italic>
<sub>m</sub>/2, then the autosomal and X-linked <italic>F</italic>
<sub>ST</sub> estimates can only be equal as the number of males tends towards zero.</p></caption></fig> <fig id="pgen-1000200-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000200.g003</object-id><label>Figure 3</label><caption><title>
<italic>p</italic>-values of Wilcoxon tests plotted in the (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) parameter space.</title><p>For each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, we applied the transformation in eq. (4), and tested whether our data on autosomal and X-linked markers were consistent, given the hypothesis defined by the set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values. (A) Surface plot of the <italic>p</italic>-values, as a function of the female fraction of effective number and the female fraction of migration rate, for the herders (11 populations). The arrow indicates the line that separates the region where <italic>p</italic>≤0.05 from that where <italic>p</italic>&gt;0.05. Non-significant <italic>p</italic>-values (<italic>p</italic>&gt;0.05) correspond to the values of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) that could not be rejected, given our data. (B) Contour plots, for the same data. The dashed line indicates the range of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values inferred from the ratio of NRY and mtDNA population structure, as obtained from the relationship: . The dotted lines correspond to the cases where <italic>N</italic>
<sub>f</sub> = <italic>N</italic>
<sub>m</sub> (vertical line) and <italic>m</italic>
<sub>f</sub> = <italic>m</italic>
<sub>m</sub> (horizontal line). (C) and (D) as (A) and (B), respectively, for the agriculturalists (10 populations).</p></caption></fig> <fig id="pgen-1000200-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000200.g004</object-id><label>Figure 4</label><caption><title>Percentage of significant tests in the (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) parameter space, for simulated data.</title><p>We chose a range of 49 (<italic>N</italic>
<sub>f</sub>
<italic>m</italic>
<sub>f</sub>/<italic>N</italic>
<sub>m</sub>
<italic>m</italic>
<sub>m</sub>) ratios, varying from 0.0004 to 2401, and for each of these ratios we chose 29 sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values. By doing this, we obtained 1421 sets of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values, represented as white dots in the right-hand side panel B, covering the whole parameter space. For each set, we simulated 100 independent datasets using a coalescent-based algorithm, and taking the same number of individuals and the same number of loci for each genetic system as in the observed data. For each dataset, we calculated the <italic>p</italic>-value for a one-sided Wilcoxon sum rank test , and for each set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values we calculated the percentage of significant <italic>p</italic>-values (at the <italic>α</italic> = 0.05 level). A. Surface plot of the proportion of significant <italic>p</italic>-values (at the <italic>α</italic> = 0.05 level), as a function of the female fraction of effective number and the female fraction of migration rate. B. Contour plot, for the same data. The dotted line, at which , represents the set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values where the autosomal and X-linked <italic>F</italic>
<sub>ST</sub>'s are equal. The theory predicts that we should only find  in the upper-right triangle defined by the dotted line. Hence, the proportion of significant <italic>p</italic>-values for any set of (<italic>N</italic>
<sub>f</sub>/<italic>N</italic>, <italic>m</italic>
<sub>f</sub>/<italic>m</italic>) values in this upper right triangle gives an indication of the power of the method.</p></caption></fig>
<table-wrap id="pgen-1000200-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000200.t001</object-id><label>Table 1</label><caption><title>Human sex-specific demography inferred from genetic data.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Region</td><td align="left" rowspan="1" colspan="1">Markers</td><td align="left" rowspan="1" colspan="1">Method</td><td align="left" rowspan="1" colspan="1">Social organization<xref ref-type="table-fn" rid="nt102">a</xref>
</td><td colspan="3" align="left" rowspan="1">Differences in demographic parameters between males and females<xref ref-type="table-fn" rid="nt103">b</xref>
</td><td align="left" rowspan="1" colspan="1">References</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Sex-biased migration</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Skewed effective population size</td><td align="left" rowspan="1" colspan="1"/></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GLOBAL</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY SNPs<xref ref-type="table-fn" rid="nt104">c</xref>
</td><td align="left" rowspan="1" colspan="1">Genetic structure (AMOVA<xref ref-type="table-fn" rid="nt105">d</xref>)</td><td align="left" rowspan="1" colspan="1">NA<xref ref-type="table-fn" rid="nt106">e</xref>
</td><td colspan="3" align="left" rowspan="1">None</td><td align="left" rowspan="1" colspan="1">
##REF##15378061##[19]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">GLOBAL</td><td align="left" rowspan="1" colspan="1">Autosomal STRs<xref ref-type="table-fn" rid="nt107">f</xref>, X-linked STRs</td><td align="left" rowspan="1" colspan="1">Genetic structure (AMOVA)</td><td align="left" rowspan="1" colspan="1">NA</td><td colspan="3" align="left" rowspan="1">None</td><td align="left" rowspan="1" colspan="1">
##REF##15601537##[20]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">GLOBAL</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY SNPs<xref ref-type="table-fn" rid="nt104">c</xref>
</td><td align="left" rowspan="1" colspan="1">Coalescent-based (TMRCA<xref ref-type="table-fn" rid="nt108">g</xref> estimates)</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1">and/or</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub> (polygyny)</td><td align="left" rowspan="1" colspan="1">
##REF##15317874##[24]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">GLOBAL<xref ref-type="table-fn" rid="nt109">h</xref>
</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY STRs+SNPs, Autosomal STRs+SNPs</td><td align="left" rowspan="1" colspan="1">Genetic structure (<italic>F</italic>
<sub>ST</sub>)</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Considered as negligible<xref ref-type="table-fn" rid="nt110">i</xref>
</td><td align="left" rowspan="1" colspan="1">
##REF##9806547##[3]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">GLOBAL<xref ref-type="table-fn" rid="nt109">h</xref>
</td><td align="left" rowspan="1" colspan="1">NRY SNPs</td><td align="left" rowspan="1" colspan="1">Coalescent-based (mismatch distributions)</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">Not considered<xref ref-type="table-fn" rid="nt111">j</xref>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub> (polygyny)</td><td align="left" rowspan="1" colspan="1">
##REF##12962309##[23]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">India</td><td align="left" rowspan="1" colspan="1">mtDNA</td><td align="left" rowspan="1" colspan="1">Genetic structure (<italic>R</italic>
<sub>ST</sub>, haplotype sharing)</td><td align="left" rowspan="1" colspan="1">Endogamy, patrilocality</td><td colspan="3" align="left" rowspan="1">None</td><td align="left" rowspan="1" colspan="1">
##REF##16617372##[21]##
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NRY STRs</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Endogamy, matrilocality</td><td colspan="3" align="left" rowspan="1">None</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Sinai peninsula</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY</td><td align="left" rowspan="1" colspan="1">Genetic diversity</td><td align="left" rowspan="1" colspan="1">Endogamy and rare patrilocal exogamy, polygyny</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1">and/or</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub> (polygyny)</td><td align="left" rowspan="1" colspan="1">
##REF##8751879##[4]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">West New Guinea</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY STRs+SNPs</td><td align="left" rowspan="1" colspan="1">Genetic structure and diversity (<italic>F</italic>
<sub>ST</sub>, <italic>R</italic>
<sub>ST</sub>, haplotype diversity)</td><td align="left" rowspan="1" colspan="1">Exogamy, patrilocality, patrilineality, polygyny</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1">and/or</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub> (polygyny, warfare)</td><td align="left" rowspan="1" colspan="1">
##REF##12532283##[7]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">Sub-Saharan Africa</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY STRs+SNPs</td><td align="left" rowspan="1" colspan="1">Genetic structure (AMOVA)</td><td align="left" rowspan="1" colspan="1">FPP<xref ref-type="table-fn" rid="nt112">k</xref>: patrilocality, high polygyny</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1">and/or</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub> (polygyny)</td><td align="left" rowspan="1" colspan="1">
##REF##15190128##[15]##
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">HGP<xref ref-type="table-fn" rid="nt113">l</xref>: moderate patrilocality, low polygyny</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&lt;<italic>m</italic>
<sub>m</sub> (multilocality)</td><td align="left" rowspan="1" colspan="1">and/or</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&lt;<italic>N</italic>
<sub>m</sub>
</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Thailand</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY STRs</td><td align="left" rowspan="1" colspan="1">Coalescent-based (Approximate Bayesian Computation)</td><td align="left" rowspan="1" colspan="1">Patrilocality</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1">and/or</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1">
##REF##15894624##[16]##
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Matrilocality</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&lt;<italic>m</italic>
<sub>m</sub> (matrilocality)</td><td align="left" rowspan="1" colspan="1">and/or</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&lt;<italic>N</italic>
<sub>m</sub> (matrilocality)</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Eastern North America</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY STRs+SNPs</td><td align="left" rowspan="1" colspan="1">Genetic structure (AMOVA), coalescent-based (MIGRATE<xref ref-type="table-fn" rid="nt114">m</xref>)</td><td align="left" rowspan="1" colspan="1">Patrilocality, patrilineality</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1">and/or</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1">
##REF##16916941##[17]##
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Matrilocality, matriliny</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&lt;<italic>m</italic>
<sub>m</sub> (matrilocality)</td><td align="left" rowspan="1" colspan="1">and/or</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&lt;<italic>N</italic>
<sub>m</sub> (matrilocality)</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Central Asia (pastoral populations)</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY STRs</td><td align="left" rowspan="1" colspan="1">Genetic structure and diversity (AMOVA, <italic>R</italic>
<sub>ST</sub>)</td><td align="left" rowspan="1" colspan="1">Exogamy, patrilineality</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilineality, exogamy)</td><td align="left" rowspan="1" colspan="1">and/or</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub> (patrilineality,VRS<xref ref-type="table-fn" rid="nt115">n</xref>)</td><td align="left" rowspan="1" colspan="1">
##REF##17208185##[11]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">New Britain</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY SNPs, X-linked loci</td><td align="left" rowspan="1" colspan="1">Coalescent-based (θ<xref ref-type="table-fn" rid="nt116">o</xref> and TMRCA estimates)</td><td align="left" rowspan="1" colspan="1">No strong endogamy, ambilocality, polygyny</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&lt;<italic>m</italic>
<sub>m</sub>
</td><td align="left" rowspan="1" colspan="1">and</td><td align="left" rowspan="1" colspan="1">
<italic>N</italic>
<sub>f</sub>&gt;<italic>N</italic>
<sub>m</sub> (polygyny)</td><td align="left" rowspan="1" colspan="1">
##UREF##1##[25]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">Central Asia</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY STRs</td><td align="left" rowspan="1" colspan="1">Genetic structure (AMOVA)</td><td align="left" rowspan="1" colspan="1">Exogamy, patrilocality, polygyny</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Considered as negligible</td><td align="left" rowspan="1" colspan="1">
##REF##10364534##[5]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">Thailand</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY STRs</td><td align="left" rowspan="1" colspan="1">Genetic structure and diversity (haplotype diversity, <italic>R</italic>
<sub>ST</sub>)</td><td align="left" rowspan="1" colspan="1">Patrilocality</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Considered as negligible</td><td align="left" rowspan="1" colspan="1">
##REF##11528385##[14]##
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Matrilocality</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&lt;<italic>m</italic>
<sub>m</sub> (matrilocality)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Considered as negligible</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Sub-Saharan Africa<xref ref-type="table-fn" rid="nt109">h</xref>
</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY SNPs</td><td align="left" rowspan="1" colspan="1">Genetic structure and diversity (haplotype diversity, AMOVA)</td><td align="left" rowspan="1" colspan="1">NA<xref ref-type="table-fn" rid="nt104">c</xref>
</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&lt;<italic>m</italic>
<sub>m</sub>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Not considered</td><td align="left" rowspan="1" colspan="1">
##REF##11420360##[22]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">Continental Asia<xref ref-type="table-fn" rid="nt109">h</xref>
</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY SNPs</td><td align="left" rowspan="1" colspan="1">Genetic structure (<italic>F</italic>
<sub>ST</sub>)</td><td align="left" rowspan="1" colspan="1">NA<xref ref-type="table-fn" rid="nt104">c</xref>
</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Not considered</td><td align="left" rowspan="1" colspan="1">
##REF##12012367##[6]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">Russia</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY SNPs</td><td align="left" rowspan="1" colspan="1">Genetic structure (<italic>F</italic>
<sub>ST</sub>)</td><td align="left" rowspan="1" colspan="1">Patrilocality, patrilineality</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Not considered</td><td align="left" rowspan="1" colspan="1">
##REF##15974299##[8]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">Caucasus</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY SNPs</td><td align="left" rowspan="1" colspan="1">Genetic structure (AMOVA)</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Not considered</td><td align="left" rowspan="1" colspan="1">
##REF##15180701##[9]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">Turkey</td><td align="left" rowspan="1" colspan="1">mtDNA, NRY STRs+SNPs</td><td align="left" rowspan="1" colspan="1">Genetic structure (AMOVA)</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">
<italic>m</italic>
<sub>f</sub>&gt;<italic>m</italic>
<sub>m</sub> (patrilocality)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Not considered</td><td align="left" rowspan="1" colspan="1">
##REF##15996169##[10]##
</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000200-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000200.t002</object-id><label>Table 2</label><caption><title>Sample description.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Sampled populations (area)</td><td align="left" rowspan="1" colspan="1">Acronym</td><td align="left" rowspan="1" colspan="1">Location</td><td align="left" rowspan="1" colspan="1">Long.</td><td align="left" rowspan="1" colspan="1">Lat.</td><td align="left" rowspan="1" colspan="1">
<italic>n</italic>
<sub>X</sub>
</td><td align="left" rowspan="1" colspan="1">
<italic>n</italic>
<sub>A</sub>
</td><td align="left" rowspan="1" colspan="1">
<italic>n</italic>
<sub>Y</sub>
</td><td align="left" rowspan="1" colspan="1">
<italic>n</italic>
<sub>mt</sub>
</td></tr></thead><tbody><tr><td colspan="9" align="left" rowspan="1">
<bold>Bilineal agriculturalists</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Tajiks (Samarkand)</td><td align="left" rowspan="1" colspan="1">TJA</td><td align="left" rowspan="1" colspan="1">Uzbekistan/Tajikistan border</td><td align="left" rowspan="1" colspan="1">
<bold>39.54</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>66.89</bold>
</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">31</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">32</td></tr><tr><td align="left" rowspan="1" colspan="1">Tajiks (Samarkand)</td><td align="left" rowspan="1" colspan="1">TJU</td><td align="left" rowspan="1" colspan="1">Uzbekistan/Tajikistan border</td><td align="left" rowspan="1" colspan="1">
<bold>39.5</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>67.27</bold>
</td><td align="left" rowspan="1" colspan="1">27</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">29</td></tr><tr><td align="left" rowspan="1" colspan="1">Tajiks (Ferghana)</td><td align="left" rowspan="1" colspan="1">TJR</td><td align="left" rowspan="1" colspan="1">Tajikistan/Kyrgyzstan border</td><td align="left" rowspan="1" colspan="1">
<bold>40.36</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>71.28</bold>
</td><td align="left" rowspan="1" colspan="1">30</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">29</td></tr><tr><td align="left" rowspan="1" colspan="1">Tajiks (Ferghana)</td><td align="left" rowspan="1" colspan="1">TJK</td><td align="left" rowspan="1" colspan="1">Tajikistan/Kyrgyzstan border</td><td align="left" rowspan="1" colspan="1">
<bold>40.25</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>71.87</bold>
</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1">40</td></tr><tr><td align="left" rowspan="1" colspan="1">Tajiks (Gharm)</td><td align="left" rowspan="1" colspan="1">TJE</td><td align="left" rowspan="1" colspan="1">Northern Tajikistan</td><td align="left" rowspan="1" colspan="1">
<bold>39.12</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>70.67</bold>
</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">27</td><td align="left" rowspan="1" colspan="1">31</td></tr><tr><td align="left" rowspan="1" colspan="1">Tajiks (Gharm)</td><td align="left" rowspan="1" colspan="1">TJN</td><td align="left" rowspan="1" colspan="1">Western Tajikistan</td><td align="left" rowspan="1" colspan="1">
<bold>38.09</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>68.81</bold>
</td><td align="left" rowspan="1" colspan="1">33</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">30</td><td align="left" rowspan="1" colspan="1">35</td></tr><tr><td align="left" rowspan="1" colspan="1">Tajiks (Gharm)</td><td align="left" rowspan="1" colspan="1">TJT</td><td align="left" rowspan="1" colspan="1">Northern Tajikistan</td><td align="left" rowspan="1" colspan="1">
<bold>39.11</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>70.86</bold>
</td><td align="left" rowspan="1" colspan="1">31</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">30</td><td align="left" rowspan="1" colspan="1">32</td></tr><tr><td align="left" rowspan="1" colspan="1">Tajiks (Penjinkent)</td><td align="left" rowspan="1" colspan="1">TDS</td><td align="left" rowspan="1" colspan="1">Uzbekistan/Tajikistan border</td><td align="left" rowspan="1" colspan="1">
<bold>39.28</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>67.81</bold>
</td><td align="left" rowspan="1" colspan="1">30</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">31</td><td align="left" rowspan="1" colspan="1">31</td></tr><tr><td align="left" rowspan="1" colspan="1">Tajiks (Penjinkent)</td><td align="left" rowspan="1" colspan="1">TDU</td><td align="left" rowspan="1" colspan="1">Uzbekistan/Tajikistan border</td><td align="left" rowspan="1" colspan="1">
<bold>39.44</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>68.26</bold>
</td><td align="left" rowspan="1" colspan="1">40</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">31</td><td align="left" rowspan="1" colspan="1">40</td></tr><tr><td align="left" rowspan="1" colspan="1">Tajiks (Yagnobs from Douchambe)</td><td align="left" rowspan="1" colspan="1">TJY</td><td align="left" rowspan="1" colspan="1">Western Tajikistan</td><td align="left" rowspan="1" colspan="1">
<bold>38.57</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>68.78</bold>
</td><td align="left" rowspan="1" colspan="1">39</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">36</td><td align="left" rowspan="1" colspan="1">40</td></tr><tr><td colspan="9" align="left" rowspan="1">
<bold>Patrilineal herders with a semi-nomadic lifestyle</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Karakalpaks (Qongrat from Karakalpakia)</td><td align="left" rowspan="1" colspan="1">KKK</td><td align="left" rowspan="1" colspan="1">Western Uzbekistan</td><td align="left" rowspan="1" colspan="1">
<bold>43.77</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>59.02</bold>
</td><td align="left" rowspan="1" colspan="1">56</td><td align="left" rowspan="1" colspan="1">45</td><td align="left" rowspan="1" colspan="1">54</td><td align="left" rowspan="1" colspan="1">55</td></tr><tr><td align="left" rowspan="1" colspan="1">Karakalpaks (On Tört Uruw from Karakalpakia)</td><td align="left" rowspan="1" colspan="1">OTU</td><td align="left" rowspan="1" colspan="1">Western Uzbekistan</td><td align="left" rowspan="1" colspan="1">
<bold>42.94</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>59.78</bold>
</td><td align="left" rowspan="1" colspan="1">49</td><td align="left" rowspan="1" colspan="1">45</td><td align="left" rowspan="1" colspan="1">54</td><td align="left" rowspan="1" colspan="1">53</td></tr><tr><td align="left" rowspan="1" colspan="1">Kazaks (Karakalpakia)</td><td align="left" rowspan="1" colspan="1">KAZ</td><td align="left" rowspan="1" colspan="1">Western Uzbekistan</td><td align="left" rowspan="1" colspan="1">
<bold>43.04</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>58.84</bold>
</td><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">49</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">50</td></tr><tr><td align="left" rowspan="1" colspan="1">Kazaks (Bukara)</td><td align="left" rowspan="1" colspan="1">LKZ</td><td align="left" rowspan="1" colspan="1">Southern Uzbekistan</td><td align="left" rowspan="1" colspan="1">
<bold>40.08</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>63.56</bold>
</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">31</td></tr><tr><td align="left" rowspan="1" colspan="1">Kyrgyz (Andijan)</td><td align="left" rowspan="1" colspan="1">KRA</td><td align="left" rowspan="1" colspan="1">Tajikistan/Kyrgyzstan border</td><td align="left" rowspan="1" colspan="1">
<bold>40.77</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>72.31</bold>
</td><td align="left" rowspan="1" colspan="1">31</td><td align="left" rowspan="1" colspan="1">45</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">48</td></tr><tr><td align="left" rowspan="1" colspan="1">Kyrgyz (Narin)</td><td align="left" rowspan="1" colspan="1">KRG</td><td align="left" rowspan="1" colspan="1">Middle Kyrgyzstan</td><td align="left" rowspan="1" colspan="1">
<bold>41.6</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>75.8</bold>
</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">20</td></tr><tr><td align="left" rowspan="1" colspan="1">Kyrgyz (Narin)</td><td align="left" rowspan="1" colspan="1">KRM</td><td align="left" rowspan="1" colspan="1">Middle Kyrgyzstan</td><td align="left" rowspan="1" colspan="1">
<bold>41.45</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>76.22</bold>
</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">26</td></tr><tr><td align="left" rowspan="1" colspan="1">Kyrgyz (Narin)</td><td align="left" rowspan="1" colspan="1">KRL</td><td align="left" rowspan="1" colspan="1">Middle Kyrgyzstan</td><td align="left" rowspan="1" colspan="1">
<bold>41.36</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>75.5</bold>
</td><td align="left" rowspan="1" colspan="1">36</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Kyrgyz (Narin)</td><td align="left" rowspan="1" colspan="1">KRB</td><td align="left" rowspan="1" colspan="1">Middle Kyrgyzstan</td><td align="left" rowspan="1" colspan="1">
<bold>41.25</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>76</bold>
</td><td align="left" rowspan="1" colspan="1">31</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Kyrgyz (Issyk Kul)</td><td align="left" rowspan="1" colspan="1">KRT</td><td align="left" rowspan="1" colspan="1">Eastern Kyrgyzstan</td><td align="left" rowspan="1" colspan="1">
<bold>42.16</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>77.57</bold>
</td><td align="left" rowspan="1" colspan="1">33</td><td align="left" rowspan="1" colspan="1">37</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">Turkmen (Karakalpakia)</td><td align="left" rowspan="1" colspan="1">TUR</td><td align="left" rowspan="1" colspan="1">Western Uzbekistan</td><td align="left" rowspan="1" colspan="1">
<bold>41.55</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>60.63</bold>
</td><td align="left" rowspan="1" colspan="1">42</td><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">51</td><td align="left" rowspan="1" colspan="1">51</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000200-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000200.t003</object-id><label>Table 3</label><caption><title>Level of diversity and differentiation for NRY markers and mtDNA.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">NRY markers</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">
<italic>F</italic>
<sub>ST</sub>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Locus name</td><td align="left" rowspan="1" colspan="1">Allelic richness (AR)</td><td align="left" rowspan="1" colspan="1">
<italic>H</italic>
<sub>e</sub>
</td><td align="left" rowspan="1" colspan="1">Herders</td><td align="left" rowspan="1" colspan="1">Agriculturalists</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">DYS426</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">0.500</td><td align="left" rowspan="1" colspan="1">0.3326</td><td align="left" rowspan="1" colspan="1">0.0068</td></tr><tr><td align="left" rowspan="1" colspan="1">DYS393</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">0.492</td><td align="left" rowspan="1" colspan="1">0.1095</td><td align="left" rowspan="1" colspan="1">0.0517</td></tr><tr><td align="left" rowspan="1" colspan="1">DYS390</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">0.739</td><td align="left" rowspan="1" colspan="1">0.1229</td><td align="left" rowspan="1" colspan="1">0.1253</td></tr><tr><td align="left" rowspan="1" colspan="1">DYS385 a/b</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">0.858</td><td align="left" rowspan="1" colspan="1">0.1414</td><td align="left" rowspan="1" colspan="1">0.0278</td></tr><tr><td align="left" rowspan="1" colspan="1">DYS388</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">0.531</td><td align="left" rowspan="1" colspan="1">0.3003</td><td align="left" rowspan="1" colspan="1">0.0736</td></tr><tr><td align="left" rowspan="1" colspan="1">DYS19</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">0.743</td><td align="left" rowspan="1" colspan="1">0.1081</td><td align="left" rowspan="1" colspan="1">0.1310</td></tr><tr><td align="left" rowspan="1" colspan="1">DYS392</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">0.516</td><td align="left" rowspan="1" colspan="1">0.1345</td><td align="left" rowspan="1" colspan="1">0.0701</td></tr><tr><td align="left" rowspan="1" colspan="1">DYS391</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">0.495</td><td align="left" rowspan="1" colspan="1">0.2533</td><td align="left" rowspan="1" colspan="1">0.0686</td></tr><tr><td align="left" rowspan="1" colspan="1">DYS389I</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">0.541</td><td align="left" rowspan="1" colspan="1">0.1537</td><td align="left" rowspan="1" colspan="1">0.1395</td></tr><tr><td align="left" rowspan="1" colspan="1">DYS439</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">0.725</td><td align="left" rowspan="1" colspan="1">0.1638</td><td align="left" rowspan="1" colspan="1">0.0291</td></tr><tr><td align="left" rowspan="1" colspan="1">DYS389II</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">0.763</td><td align="left" rowspan="1" colspan="1">0.1556</td><td align="left" rowspan="1" colspan="1">0.0395</td></tr></tbody></table><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">mtDNA</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">
<italic>F</italic>
<sub>ST</sub>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Locus name</td><td align="left" rowspan="1" colspan="1">Polymorphic sites</td><td align="left" rowspan="1" colspan="1">
<italic>H</italic>
<sub>e</sub>
</td><td align="left" rowspan="1" colspan="1">Herders</td><td align="left" rowspan="1" colspan="1">Agriculturalists</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">HVS-I</td><td align="left" rowspan="1" colspan="1">121</td><td align="left" rowspan="1" colspan="1">0.0156</td><td align="left" rowspan="1" colspan="1">0.0098</td><td align="left" rowspan="1" colspan="1">0.0343</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000200-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000200.t004</object-id><label>Table 4</label><caption><title>Level of diversity and differentiation for X-linked and autosomal markers.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">
<italic>F</italic>
<sub>ST</sub>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Locus name</td><td align="left" rowspan="1" colspan="1">Allelic richness (AR)</td><td align="left" rowspan="1" colspan="1">
<italic>H</italic>
<sub>e</sub>
</td><td align="left" rowspan="1" colspan="1">Herders</td><td align="left" rowspan="1" colspan="1">Agriculturalists</td></tr></thead><tbody><tr><td colspan="5" align="left" rowspan="1">
<bold>X-linked markers</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">CTAT014</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">0.746</td><td align="left" rowspan="1" colspan="1">0.0018</td><td align="left" rowspan="1" colspan="1">0.0225</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA124E07</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">0.847</td><td align="left" rowspan="1" colspan="1">0.0024</td><td align="left" rowspan="1" colspan="1">0.0136</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA31D10</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">0.697</td><td align="left" rowspan="1" colspan="1">0.0069</td><td align="left" rowspan="1" colspan="1">0.0007</td></tr><tr><td align="left" rowspan="1" colspan="1">ATA28C05</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">0.722</td><td align="left" rowspan="1" colspan="1">0.0086</td><td align="left" rowspan="1" colspan="1">0.0179</td></tr><tr><td align="left" rowspan="1" colspan="1">AFM150xf10</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">0.832</td><td align="left" rowspan="1" colspan="1">−0.0021</td><td align="left" rowspan="1" colspan="1">0.0152</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA100G03</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">0.734</td><td align="left" rowspan="1" colspan="1">−0.0019</td><td align="left" rowspan="1" colspan="1">0.0084</td></tr><tr><td align="left" rowspan="1" colspan="1">AGAT121P</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">0.593</td><td align="left" rowspan="1" colspan="1">−0.0016</td><td align="left" rowspan="1" colspan="1">0.0048</td></tr><tr><td align="left" rowspan="1" colspan="1">ATCT003</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">0.797</td><td align="left" rowspan="1" colspan="1">0.0095</td><td align="left" rowspan="1" colspan="1">0.0261</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA31F01</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">0.804</td><td align="left" rowspan="1" colspan="1">0.0069</td><td align="left" rowspan="1" colspan="1">0.0053</td></tr><tr><td colspan="5" align="left" rowspan="1">
<bold>Autosomal markers</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">AFM249XC5</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">0.848</td><td align="left" rowspan="1" colspan="1">0.0080</td><td align="left" rowspan="1" colspan="1">0.0081</td></tr><tr><td align="left" rowspan="1" colspan="1">ATA10H11</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">0.680</td><td align="left" rowspan="1" colspan="1">0.0128</td><td align="left" rowspan="1" colspan="1">0.0193</td></tr><tr><td align="left" rowspan="1" colspan="1">AFM254VE1</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">0.837</td><td align="left" rowspan="1" colspan="1">0.0105</td><td align="left" rowspan="1" colspan="1">0.0086</td></tr><tr><td align="left" rowspan="1" colspan="1">AFMA218YB5</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">0.852</td><td align="left" rowspan="1" colspan="1">0.0030</td><td align="left" rowspan="1" colspan="1">0.0151</td></tr><tr><td align="left" rowspan="1" colspan="1">GGAA7G08</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">0.896</td><td align="left" rowspan="1" colspan="1">0.0096</td><td align="left" rowspan="1" colspan="1">0.0138</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA11H10</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">0.776</td><td align="left" rowspan="1" colspan="1">0.0017</td><td align="left" rowspan="1" colspan="1">0.0056</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA12A07</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">0.857</td><td align="left" rowspan="1" colspan="1">0.0001</td><td align="left" rowspan="1" colspan="1">0.0163</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA193A07</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">0.825</td><td align="left" rowspan="1" colspan="1">0.0064</td><td align="left" rowspan="1" colspan="1">0.0087</td></tr><tr><td align="left" rowspan="1" colspan="1">AFMB002ZF1</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">0.820</td><td align="left" rowspan="1" colspan="1">0.0028</td><td align="left" rowspan="1" colspan="1">0.0169</td></tr><tr><td align="left" rowspan="1" colspan="1">AFMB303ZG9</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">0.858</td><td align="left" rowspan="1" colspan="1">0.0090</td><td align="left" rowspan="1" colspan="1">0.0148</td></tr><tr><td align="left" rowspan="1" colspan="1">ATA34G06</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">0.675</td><td align="left" rowspan="1" colspan="1">0.0088</td><td align="left" rowspan="1" colspan="1">0.0132</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA72G09</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">0.884</td><td align="left" rowspan="1" colspan="1">−0.0023</td><td align="left" rowspan="1" colspan="1">0.0131</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA22F11</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">0.897</td><td align="left" rowspan="1" colspan="1">0.0152</td><td align="left" rowspan="1" colspan="1">0.0144</td></tr><tr><td align="left" rowspan="1" colspan="1">GGAA6D03</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">0.831</td><td align="left" rowspan="1" colspan="1">0.0048</td><td align="left" rowspan="1" colspan="1">0.0176</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA88H02</td><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">0.892</td><td align="left" rowspan="1" colspan="1">0.0063</td><td align="left" rowspan="1" colspan="1">0.0056</td></tr><tr><td align="left" rowspan="1" colspan="1">SE30</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">0.762</td><td align="left" rowspan="1" colspan="1">0.0084</td><td align="left" rowspan="1" colspan="1">0.0103</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA43C11</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">0.870</td><td align="left" rowspan="1" colspan="1">0.0028</td><td align="left" rowspan="1" colspan="1">0.0093</td></tr><tr><td align="left" rowspan="1" colspan="1">AFM203YG9</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">0.753</td><td align="left" rowspan="1" colspan="1">0.0105</td><td align="left" rowspan="1" colspan="1">0.0084</td></tr><tr><td align="left" rowspan="1" colspan="1">AFM157XG3</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">0.753</td><td align="left" rowspan="1" colspan="1">0.0147</td><td align="left" rowspan="1" colspan="1">0.0196</td></tr><tr><td align="left" rowspan="1" colspan="1">UT2095</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">0.738</td><td align="left" rowspan="1" colspan="1">0.0032</td><td align="left" rowspan="1" colspan="1">0.0112</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA28D01</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">0.896</td><td align="left" rowspan="1" colspan="1">0.0156</td><td align="left" rowspan="1" colspan="1">0.0139</td></tr><tr><td align="left" rowspan="1" colspan="1">GGAA4B09</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">0.707</td><td align="left" rowspan="1" colspan="1">0.0034</td><td align="left" rowspan="1" colspan="1">0.0208</td></tr><tr><td align="left" rowspan="1" colspan="1">ATA3A07</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">0.746</td><td align="left" rowspan="1" colspan="1">0.0078</td><td align="left" rowspan="1" colspan="1">0.0070</td></tr><tr><td align="left" rowspan="1" colspan="1">AFM193XH4</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">0.716</td><td align="left" rowspan="1" colspan="1">0.0164</td><td align="left" rowspan="1" colspan="1">0.0129</td></tr><tr><td align="left" rowspan="1" colspan="1">GATA11B12</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">0.896</td><td align="left" rowspan="1" colspan="1">0.0104</td><td align="left" rowspan="1" colspan="1">0.0265</td></tr><tr><td align="left" rowspan="1" colspan="1">AFM165XC11</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">0.785</td><td align="left" rowspan="1" colspan="1">0.0058</td><td align="left" rowspan="1" colspan="1">0.0185</td></tr><tr><td align="left" rowspan="1" colspan="1">AFM248VC5</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">0.620</td><td align="left" rowspan="1" colspan="1">0.0246</td><td align="left" rowspan="1" colspan="1">0.0145</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000200-t005" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000200.t005</object-id><label>Table 5</label><caption><title>Autosomal and X-linked differentiation on jackknifed samples.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Sample removed</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">
<italic>p</italic>-value</td><td align="left" rowspan="1" colspan="1">

</td></tr></thead><tbody><tr><td colspan="5" align="left" rowspan="1">
<bold>Patrilineal groups</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">KAZ</td><td align="left" rowspan="1" colspan="1">0.0084</td><td align="left" rowspan="1" colspan="1">0.0050</td><td align="left" rowspan="1" colspan="1">0.068</td><td align="left" rowspan="1" colspan="1">1.7</td></tr><tr><td align="left" rowspan="1" colspan="1">KKK</td><td align="left" rowspan="1" colspan="1">0.0085</td><td align="left" rowspan="1" colspan="1">0.0050</td><td align="left" rowspan="1" colspan="1">0.078</td><td align="left" rowspan="1" colspan="1">1.7</td></tr><tr><td align="left" rowspan="1" colspan="1">KRA</td><td align="left" rowspan="1" colspan="1">0.0078</td><td align="left" rowspan="1" colspan="1">0.0027</td><td align="left" rowspan="1" colspan="1">0.022</td><td align="left" rowspan="1" colspan="1">2.9</td></tr><tr><td align="left" rowspan="1" colspan="1">KRB</td><td align="left" rowspan="1" colspan="1">0.0080</td><td align="left" rowspan="1" colspan="1">0.0030</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">2.7</td></tr><tr><td align="left" rowspan="1" colspan="1">KRG</td><td align="left" rowspan="1" colspan="1">0.0078</td><td align="left" rowspan="1" colspan="1">0.0035</td><td align="left" rowspan="1" colspan="1">0.037</td><td align="left" rowspan="1" colspan="1">2.2</td></tr><tr><td align="left" rowspan="1" colspan="1">KRL</td><td align="left" rowspan="1" colspan="1">0.0086</td><td align="left" rowspan="1" colspan="1">0.0038</td><td align="left" rowspan="1" colspan="1">0.018</td><td align="left" rowspan="1" colspan="1">2.3</td></tr><tr><td align="left" rowspan="1" colspan="1">KRM</td><td align="left" rowspan="1" colspan="1">0.0069</td><td align="left" rowspan="1" colspan="1">0.0023</td><td align="left" rowspan="1" colspan="1">0.018</td><td align="left" rowspan="1" colspan="1">3.0</td></tr><tr><td align="left" rowspan="1" colspan="1">KRT</td><td align="left" rowspan="1" colspan="1">0.0081</td><td align="left" rowspan="1" colspan="1">0.0044</td><td align="left" rowspan="1" colspan="1">0.047</td><td align="left" rowspan="1" colspan="1">1.8</td></tr><tr><td align="left" rowspan="1" colspan="1">LKZ</td><td align="left" rowspan="1" colspan="1">0.0088</td><td align="left" rowspan="1" colspan="1">0.0025</td><td align="left" rowspan="1" colspan="1">0.002</td><td align="left" rowspan="1" colspan="1">3.5</td></tr><tr><td align="left" rowspan="1" colspan="1">OTU</td><td align="left" rowspan="1" colspan="1">0.0089</td><td align="left" rowspan="1" colspan="1">0.0038</td><td align="left" rowspan="1" colspan="1">0.022</td><td align="left" rowspan="1" colspan="1">2.3</td></tr><tr><td align="left" rowspan="1" colspan="1">TUR</td><td align="left" rowspan="1" colspan="1">0.0054</td><td align="left" rowspan="1" colspan="1">0.0025</td><td align="left" rowspan="1" colspan="1">0.073</td><td align="left" rowspan="1" colspan="1">2.2</td></tr><tr><td colspan="5" align="left" rowspan="1">
<bold>Bilineal groups</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">TDS</td><td align="left" rowspan="1" colspan="1">0.0125</td><td align="left" rowspan="1" colspan="1">0.0109</td><td align="left" rowspan="1" colspan="1">0.443</td><td align="left" rowspan="1" colspan="1">1.1</td></tr><tr><td align="left" rowspan="1" colspan="1">TDU</td><td align="left" rowspan="1" colspan="1">0.0132</td><td align="left" rowspan="1" colspan="1">0.0153</td><td align="left" rowspan="1" colspan="1">0.705</td><td align="left" rowspan="1" colspan="1">0.9</td></tr><tr><td align="left" rowspan="1" colspan="1">TJA</td><td align="left" rowspan="1" colspan="1">0.0144</td><td align="left" rowspan="1" colspan="1">0.0123</td><td align="left" rowspan="1" colspan="1">0.109</td><td align="left" rowspan="1" colspan="1">1.2</td></tr><tr><td align="left" rowspan="1" colspan="1">TJE</td><td align="left" rowspan="1" colspan="1">0.0140</td><td align="left" rowspan="1" colspan="1">0.0133</td><td align="left" rowspan="1" colspan="1">0.148</td><td align="left" rowspan="1" colspan="1">1.1</td></tr><tr><td align="left" rowspan="1" colspan="1">TJK</td><td align="left" rowspan="1" colspan="1">0.0134</td><td align="left" rowspan="1" colspan="1">0.0131</td><td align="left" rowspan="1" colspan="1">0.457</td><td align="left" rowspan="1" colspan="1">1.0</td></tr><tr><td align="left" rowspan="1" colspan="1">TJN</td><td align="left" rowspan="1" colspan="1">0.0148</td><td align="left" rowspan="1" colspan="1">0.0144</td><td align="left" rowspan="1" colspan="1">0.387</td><td align="left" rowspan="1" colspan="1">1.0</td></tr><tr><td align="left" rowspan="1" colspan="1">TJR</td><td align="left" rowspan="1" colspan="1">0.0140</td><td align="left" rowspan="1" colspan="1">0.0141</td><td align="left" rowspan="1" colspan="1">0.401</td><td align="left" rowspan="1" colspan="1">1.0</td></tr><tr><td align="left" rowspan="1" colspan="1">TJT</td><td align="left" rowspan="1" colspan="1">0.0139</td><td align="left" rowspan="1" colspan="1">0.0121</td><td align="left" rowspan="1" colspan="1">0.225</td><td align="left" rowspan="1" colspan="1">1.1</td></tr><tr><td align="left" rowspan="1" colspan="1">TJU</td><td align="left" rowspan="1" colspan="1">0.0139</td><td align="left" rowspan="1" colspan="1">0.0127</td><td align="left" rowspan="1" colspan="1">0.283</td><td align="left" rowspan="1" colspan="1">1.1</td></tr><tr><td align="left" rowspan="1" colspan="1">TJY</td><td align="left" rowspan="1" colspan="1">0.0139</td><td align="left" rowspan="1" colspan="1">0.0116</td><td align="left" rowspan="1" colspan="1">0.259</td><td align="left" rowspan="1" colspan="1">1.2</td></tr></tbody></table></alternatives></table-wrap>
<inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(1)</label></disp-formula> <disp-formula><label>(2)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(3)</label></disp-formula> <disp-formula><label>(4)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(5)</label></disp-formula> <disp-formula><label>(6)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula>
<table-wrap-foot><fn id="nt101"><label/><p>This table summarizes the observed patterns of sex-specific differences in demographic parameters reported in a number of recent studies. The first column lists the location of the sampled populations, or indicates whether the study is conducted at a global scale. The second column gives the markers used, and the third column indicates the statistical methods employed. The fourth column provides indications on social organization, available a priori for the populations under study. In the fifth and sixth columns, the authors' interpretations of sex-specific differences in demographic parameters are given, with respect to skewed gene flow and/or effective numbers.</p></fn><fn id="nt102"><label>a</label><p>Indications on social organization, marriage rules, etc., as provided by the authors.</p></fn><fn id="nt103"><label>b</label><p>The differences in demographic parameters between males and females, as inferred by the authors, are given in terms of sex-biased gene flow, and skewed effective numbers; the authors' interpretation to the observed pattern is given in parentheses, when available.</p></fn><fn id="nt104"><label>c</label><p>Single nucleotide polymorphisms.</p></fn><fn id="nt105"><label>d</label><p>Analysis of molecular variance ##REF##1644282##[69]##.</p></fn><fn id="nt106"><label>e</label><p>Not available (no detailed information given by the authors concerning social organization, marriage rules, etc.).</p></fn><fn id="nt107"><label>f</label><p>Short tandem repeats.</p></fn><fn id="nt108"><label>g</label><p>Time to the most recent common ancestor.</p></fn><fn id="nt109"><label>h</label><p>mtDNA and NRY were not sampled in the same individuals or populations.</p></fn><fn id="nt110"><label>i</label><p>The authors discussed a possible difference in demographic parameters between males and females, but considered it as negligible.</p></fn><fn id="nt111"><label>j</label><p>The authors did not consider this pattern.</p></fn><fn id="nt112"><label>k</label><p>Food-producer populations.</p></fn><fn id="nt113"><label>l</label><p>Hunter-gatherer populations.</p></fn><fn id="nt114"><label>m</label><p>Monte Carlo Markov chain method to estimate population sizes and migration rates ##REF##11287657##[70]##.</p></fn><fn id="nt115"><label>n</label><p>Variance in Reproductive Success.</p></fn><fn id="nt116"><label>o</label><p>population-mutation parameter.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt117"><label/><p>Long., longitude; Lat., latitude. <italic>n</italic>
<sub>X</sub>, <italic>n</italic>
<sub>A</sub>, <italic>n</italic>
<sub>Y</sub> and <italic>n</italic>
<sub>mt</sub>: sample size for X-linked, autosomal, Y-linked and mitochondrial markers, respectively.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt118"><label/><p>We calculated the total allelic richness (<italic>AR</italic>) (over all populations) and the expected heterozygosity <italic>H</italic>
<sub>e</sub>
##REF##17248844##[55]## using Arlequin version 3.1 ##UREF##8##[56]##. Genetic differentiation among populations was measured both per locus and overall loci, using Weir and Cockerham's <italic>F</italic>
<sub>ST</sub> estimator ##UREF##9##[57]##, as calculated in G<sc>enepop</sc> 4.0 ##UREF##10##[58]##. We calculated the total number of polymorphic sites, the unbiased estimate of expected heterozygosity <italic>H</italic>
<sub>e</sub>
##REF##17248844##[55]##, and <italic>F</italic>
<sub>ST</sub> using Arlequin version 3.1 ##UREF##8##[56]##.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt119"><label/><p>We calculated the allelic richness (<italic>AR</italic>) and unbiased estimates of expected heterozygosity <italic>H</italic>
<sub>e</sub>
##REF##17248844##[55]##, obtained both by locus and on average with Arlequin version 3.1 ##UREF##8##[56]##. Genetic differentiation among populations was measured both per locus and overall loci, using Weir and Cockerham's <italic>F</italic>
<sub>ST</sub> estimator ##UREF##9##[57]## as calculated in G<sc>enepop</sc> 4.0 ##UREF##10##[58]##.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt120"><label/><p>For each group, we removed one sample in turn and calculated the differentiation on autosomal and X-linked markers. The <italic>p</italic>-value gives the result of a one-sided Wilcoxon sum rank test , as performed on the full dataset.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This work was supported by the Centre National de la Recherche Scientifique (CNRS) ATIP programme (to EH), by the CNRS interdisciplinary programme “Origines de l'Homme du Langage et des Langues” (OHLL) and by the European Science Foundation (ESF) EUROCORES programme “The Origin of Man, Language and Languages” (OMLL). We also thank the “Fondation pour la Recherche Médicale” (FRM) for financial support. LS is financed by the French Ministry of Higher Education and Research. MAJ is supported by a Wellcome Trust Senior Fellowship in Basic Biomedical Science (grant number 057559).</p></fn></fn-group>
<graphic id="pgen-1000200-t001-1" xlink:href="pgen.1000200.t001"/> <graphic xlink:href="pgen.1000200.g001"/> <inline-graphic xlink:href="pgen.1000200.e001.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e002.jpg" mimetype="image"/> <graphic id="pgen-1000200-t002-2" xlink:href="pgen.1000200.t002"/> <graphic id="pgen-1000200-t003-3" xlink:href="pgen.1000200.t003"/> <graphic xlink:href="pgen.1000200.e003.jpg" mimetype="image" position="float"/> <graphic xlink:href="pgen.1000200.e004.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pgen.1000200.e005.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e006.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e007.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e008.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e009.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e010.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e011.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e012.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e013.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e014.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e015.jpg" mimetype="image"/> <graphic xlink:href="pgen.1000200.g002"/> <inline-graphic xlink:href="pgen.1000200.e016.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e017.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e018.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e019.jpg" mimetype="image"/> <graphic id="pgen-1000200-t004-4" xlink:href="pgen.1000200.t004"/> <graphic xlink:href="pgen.1000200.e020.jpg" mimetype="image" position="float"/> <graphic xlink:href="pgen.1000200.e021.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pgen.1000200.e022.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e023.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e024.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e025.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e026.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e027.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e028.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e029.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e030.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e031.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e032.jpg" mimetype="image"/> <graphic xlink:href="pgen.1000200.g003"/> <inline-graphic xlink:href="pgen.1000200.e033.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e034.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e035.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e036.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e037.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e038.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e039.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e040.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e041.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e042.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e043.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e044.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e045.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e046.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e047.jpg" mimetype="image"/> <graphic xlink:href="pgen.1000200.g004"/> <inline-graphic xlink:href="pgen.1000200.e048.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e049.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e050.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e051.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e052.jpg" mimetype="image"/> <graphic id="pgen-1000200-t005-5" xlink:href="pgen.1000200.t005"/> <inline-graphic xlink:href="pgen.1000200.e053.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e054.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e055.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e056.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e057.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e058.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e059.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e060.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e061.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e062.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e063.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e064.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e065.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e066.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e067.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e068.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e069.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e070.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e071.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e072.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e073.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e074.jpg" mimetype="image"/> <graphic xlink:href="pgen.1000200.e075.jpg" mimetype="image" position="float"/> <graphic xlink:href="pgen.1000200.e076.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pgen.1000200.e077.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e078.jpg" mimetype="image"/> <inline-graphic xlink:href="pgen.1000200.e079.jpg" mimetype="image"/>
[{"label": ["13"], "element-citation": ["\n"], "surname": ["Burton", "Moore", "Whiting", "Romney"], "given-names": ["ML", "CC", "JWM", "AK"], "year": ["1996"], "article-title": ["Regions based on social structure."], "source": ["Curr Anthro"], "volume": ["37"], "fpage": ["87"], "lpage": ["123"]}, {"label": ["25"], "element-citation": ["\n"], "surname": ["Wilder", "Hammer", "Friedlaender"], "given-names": ["JA", "MF", "JS"], "year": ["2007"], "article-title": ["Extraordinary population structure among the Baining of New Britain."], "source": ["Genes, Language, and Culture History in the Southwest Pacific"], "publisher-loc": ["Oxford, UK"], "publisher-name": ["Oxford University Press"], "fpage": ["199"], "lpage": ["207"]}, {"label": ["34"], "element-citation": ["\n"], "surname": ["White"], "given-names": ["DR"], "year": ["1988"], "article-title": ["Rethinking polygyny: co-wives, codes, and cultural systems."], "source": ["Curr Anthro"], "volume": ["29"], "fpage": ["529"], "lpage": ["558"]}, {"label": ["35"], "element-citation": ["\n"], "surname": ["Heider", "S", "S"], "given-names": ["KG", "G", "L"], "year": ["1997"], "article-title": ["Grand valley Dani: peaceful warriors."], "source": ["Case studies in cultural anthropology"], "publisher-loc": ["Forth Worth, Texas"], "publisher-name": ["Harcourt Brace College Publishers"]}, {"label": ["40"], "element-citation": ["\n"], "surname": ["Beaumont", "Nichols"], "given-names": ["M", "RA"], "year": ["1996"], "article-title": ["Evaluating loci for use in the genetic analysis of population structure."], "source": ["Proc R Soc Lond"], "volume": ["263"], "fpage": ["1619"], "lpage": ["1626"]}, {"label": ["45"], "element-citation": ["\n"], "surname": ["Wright"], "given-names": ["S"], "year": ["1939"], "source": ["Statistical genetics in relation to evolution. Actualit\u00e9s scientifiques et industrielles 802 Expos\u00e9s de Biom\u00e9trie et de Statistique Biologique XIII"], "publisher-loc": ["Paris"], "publisher-name": ["Hermann et Cie"]}, {"label": ["46"], "element-citation": ["\n"], "surname": ["Tamisier"], "given-names": ["JC"], "year": ["1998"], "source": ["Dictionnaire des peuples. Soci\u00e9t\u00e9s d'Afrique, d'Am\u00e9rique, d'Asie et d'Oc\u00e9anie"], "publisher-loc": ["Paris"], "publisher-name": ["Larousse-Bordas"]}, {"label": ["50"], "element-citation": ["\n"], "surname": ["Maniatis", "Fritsh", "S"], "given-names": ["T", "EF", "J"], "year": ["1982"], "source": ["Molecular cloning. A laboratory manual"], "publisher-loc": ["New York"], "publisher-name": ["Cold Spring Laboratory"]}, {"label": ["56"], "element-citation": ["\n"], "surname": ["Excoffier", "Laval", "Schneider"], "given-names": ["L", "LG", "S"], "year": ["2005"], "article-title": ["Arlequin ver. 3.0: An integrated software package for population genetics data analysis."], "source": ["Evol Bioinfo Online"], "volume": ["1"], "fpage": ["47"], "lpage": ["50"]}, {"label": ["57"], "element-citation": ["\n"], "surname": ["Weir", "Cockerham"], "given-names": ["BS", "CC"], "year": ["1984"], "article-title": ["Estimating "], "italic": ["F"], "source": ["Evolution"], "volume": ["38"], "fpage": ["1358"], "lpage": ["1370"]}, {"label": ["58"], "element-citation": ["\n"], "surname": ["Rousset"], "given-names": ["F"], "year": ["2008"], "article-title": ["Genepop'007: a complete re-implementation of the genepop software for Windows and Linux."], "source": ["Mol Ecol Res"], "volume": ["8"], "fpage": ["103"], "lpage": ["106"]}, {"label": ["59"], "element-citation": ["\n"], "surname": ["DiCiccio", "Efron"], "given-names": ["TJ", "B"], "year": ["1996"], "article-title": ["Bootstrap confidence intervals."], "source": ["Statistical Science"], "volume": ["11"], "fpage": ["189"], "lpage": ["228"]}, {"label": ["61"], "element-citation": ["\n"], "collab": ["R Development Core Team"], "year": ["2007"], "source": ["R: A Language and Environment for Statistical Computing"], "publisher-loc": ["Vienna"], "publisher-name": ["R Foundation for Statistical Computing"]}, {"label": ["64"], "element-citation": ["\n"], "surname": ["Rousset"], "given-names": ["F"], "year": ["2004"], "source": ["Genetic Structure and Selection in Subdivided Populations"], "publisher-loc": ["Princeton, New Jersey"], "publisher-name": ["Princeton University Press"]}, {"label": ["66"], "element-citation": ["\n"], "surname": ["Gelman", "Carlin", "Stern", "Rubin"], "given-names": ["A", "JB", "HS", "DB"], "year": ["2004"], "source": ["Bayesian Data Analysis. Second Edition"], "publisher-loc": ["New York"], "publisher-name": ["Chapman & Hall/CRC"]}]
PLoS Genet. 2008 Sep 26; 4(9):e1000200