Revisiting the peopling of Japan: an admixture perspective

The development and spread of farming, referred to as the ‘Neolithic transition,’ was one of the major demographic events of human prehistory.¹This process took place independently in different geographical areas, each one most likely associated with different demographic changes and with different domesticated animals and plants. In principle, each of these changes can be described as a process by which at least two human groups (Paleolithic hunter-gatherers (HG) and Neolithic farmers) admixed to different extents. These processes can be seen as admixture models and although they have been used to study the Neolithic transition in Europe,^2,3,4this has not been the case for Asia. Here, we focus on Eastern Asia, where the transition to agriculture has long been controversial, specifically regarding the prehistory of Japan.^5,6,7,8

Archeological data suggest that there were probably two migratory waves of incoming people, both from the Asian continent to Japan. The first migration took place c. 38000–37000 BP, before the Pleistocene land bridges were submerged,⁹and later gave rise to the Jomon culture (12000 BP).¹⁰Although they were an HG society, the Jomon were the holders of one of the oldest pottery cultures known in the world and probably also led a sedentary or semi-sedentary life, well before showing any clear evidence of having developed agriculture.^1,11A long time after this period, c. 2300 BP, a second wave of people, together with a ‘wet rice culture,’ weaving and metalwork, entered the southern Kyushu island (Figure 1), through the Korean Peninsula,¹²and then spread northeastward, starting the Yayoi period.

The transformation and the replacement models represent the two opposite extremes of the demographic models that have been proposed to explain the peopling of Japan and the contribution of both Jomon and Yayoi populations to modern Japanese. While the latter model claims that modern Japanese should be descendents of the incoming Yayoi who replaced completely the Jomon people,⁵the former entails a movement of the Yayoi culture and ideas rather than people, with consequently no genetic contribution of the Yayoi to modern Japanese.⁶然而,现实一定是那么极端和currently it is widely accepted that the modern Japanese are the result of an admixture between the two populations that produced both the Jomon and Yayoi cultures. This was suggested by Hanihara⁷and Matsumura⁸on the basis of dental and cranial characteristics, and more recently by a number of authors who used genetic data,^{13,14,15,16,17}including ancient DNA.^18,19

As one of the few points on which all studies agree is that at least two human groups admixed at some point in the past, a simplified way to explain the data is the use of an admixture approach. However, one of the limitations, of most of admixture models, is that they usually ignore genetic drift since the admixture event. This is why we used an approach²⁰that has already been applied to address the Neolithic transition in Europe^2,21,22and where drift is explicitly accounted for. We expect that the admixture process varied geographically, as the incomers (early farmers) were meeting and admixing with the local populations and their descendants were themselves mixing with other populations. Although the admixture process must have been complex, we can predict that a correlation should exist between the admixture level at a particular location, measured by the contribution of one parental population, and the geographical distance from that parental population, as has been shown in Europe.²We also expect that this relationship would not hold, if the same analysis was performed using parental populations that could not have contributed. We note here that to carry out the admixture analysis, two modern populations are chosen to approximate the haplotype frequencies of the original parental populations (Jomon and Yayoi). The choice of these parental populations is based on archeological evidence and is described in the Materials and methods section.

Thus, the aim of this study was to determine whether an admixture approach could be fruitful to studying the Neolithic transition in Japan. To do this we analyzed Y-chromosomal data from the literature, using different ‘parental’ populations, to test different hypotheses. In a first set of analyses, the parental populations were chosen among a set of Asian populations (see below for details). The data were also analyzed by using, as a negative test, populations that were unlikely to have contributed to the gene pool of modern Japanese, namely a European (Sardinia) and a geographically closer (Oceania) population, and for which comparable Y-chromosomal data were available. Altogether, we show that admixture models can indeed provide interesting insights into the peopling of Japan. In particular, our results strongly suggest that the Yayoi immigrants spread by a process similar to the demic diffusion, first proposed for Europe by Ammerman and Cavalli-Sforza.²³

Materials and methods

Populations used

The analyses presented in this study were based on published non-recombining Y-chromosome data of Japanese and other Asian populations. A total of 275 individuals, representing each of the Japanese islands (Figure 1), were analyzed: Ainu (20), Aomori (26), Shizuoka (61), Tokushima (70), Kyushu (53) and Okinawa (45). All the Japanese data were published by Hammeret al.,¹⁵except the Ainu data that were pooled with data from Tajimaet al.²⁴Mainland Asian data¹⁵were obtained for populations from Northeast (441), Southeast (683) and Central (419) Asia and also a sample from Korea²⁵(43). We also used two additional populations, Sardinia²⁶(77) and Oceania¹⁵(209),作为掺合料莫父母的人群del used (see below). Y-chromosome binary haplogroups were defined by the analysis of the binary polymorphisms described in Hammeret al.¹⁵The Y-chromosome lineages from Japan, mainland Asia, Korea, Oceania and Sardinia followed the haplogroups nomenclature of the Y Chromosome Consortium.²⁷

The admixture model

The admixture method used assumes that an ‘admixed’ or ‘hybrid’ population (H), of sizeN_h, is the result of the admixture of two independent parental populations,P₁andP₂, of sizeN₁andN₂,Tgenerations ago, with respective contributionsp₁andp₂(p₂=1−p₁). After the admixture event, the three populations are isolated and assumed to evolve independently under pure genetic drift. The advantage of this model, and of the associated inference methods, is that (i) the three populations have differentN_i(whereican be 1, 2 orh)和(2)漂移和掺合料分离。我s important to note that, by explicitly accounting for drift after the admixture event, the method allows for present-day samples from the parental populations to have drifted significantly from the original unknown parental populations. In addition, the method does not fix the original parental allele frequencies. Instead, they were allowed to vary and this uncertainty is explicitly taken into account. A Bayesian full-likelihood method based on this model was developed by Chikhiet al.,²⁰implemented in the LEA (Likelihood-based Estimation of Admixture) software freely available atwww.igc.gulbenkian.pt/static/docs/LEA.zip.²⁸LEA implements a Monte Carlo Markov Chain algorithm to jointly infer all the parameters of the admixture model, including the ancestral allelic configurations that are compatible with the present, observed allelic frequencies. For each analysis, LEA was run for 300000 steps, as it has been shown that it is enough to reach an equilibrium for the Y-chromosomal data.^2,20

Choice of parental populations

For simplicity and consistency, theP₁population was always used to represent the HG or Jomon, whereas the populationP₂was used to represent the farmers of the Yayoi period. Hence, the parameterp₁represents the ‘Jomon’ contribution, at the moment of admixture, whereasp₂would represent the ‘Yayoi’ contribution. However, like all admixture methods, it requires that these parental populations be defined. Although it is unlikely that today’s populations are direct descendents from any of the original groups, we can use current archeological and anthropological data to identify populations that are likely to be less admixed, and use them as descendents from the original parental populations. It is noteworthy that if there has been a lot of admixture in these parental populations, the general effect should be to blur the original signal, and make it less clear. Therefore, any signal observed today should be an indication that some information is still present in the data. Although the Jomon culture has almost been replaced across Japan, there are some indigenous minority ethnic groups who live in the peripheral areas of Japan, which are considered descendents of this ancient culture.^7,16,17,24Those are the Ainu people, in the northern part of the Hokkaido Island, and the Ryukyuans, in the southern Ryukyu Islands. Moreover, the Ainu lived in relative isolation until the end of the nineteenth century²⁹and show unique physical characteristics; such as hairiness, wavy hair and deep-set eyes, which are different from those of most Japanese. In contrast, the Ryukyuan kingdom had past relations with mainland Japan since medieval times, with possibly frequent gene flow,^30,31but it is thought to have nevertheless maintained genetic differentiation from mainland Japan.³²For these reasons, the admixture analyses were performed using either the Ainu or the Ryukyuans, the latter represented by the Okinawa sample, as descendents of theP₁population, in the different analyses. For the descendents of the Yayoi (considered theP₂population), different parental populations from mainland Asia were also used, namely NEA (northeast Asia), SEA (southeast Asia), CAS (central Asia) and Korea.

To determine whether our approach was robust to incorrect specification of the parental populations, we also used asP₂two populations that are unlikely to have contributed to the gene pool of the Japanese: one from Europe (Sardinia) and the other from a closer geographical area, Oceania. We expected that there should be no correlation (or at least much less) between admixture and geographical distances in these cases.

Altogether, each of the four Japanese ‘admixed’ populations (Aomori, Shizuoka, Tokushima, Kyushu) were analyzed using two populations forP₁(Ainu and Okinawa), six populations forP₂(including Sardinia and Oceania), making a total of 12 different sets of admixture analyses. In addition, for each admixture analysis, the parental populations were also considered as ‘admixed’ populations. For example, we used the Ainu asP₁, asH, against the six differentP₂populations. This kind of analysis allowed us to quantify the uncertainty around the estimation ofp₁, as the hybrid and one parental population (hereP₁) are exactly identical. Thus, thep₁后验应该总是有一个模式或相等close to one, with a variance related to both the sample size and drift since the admixture event. Of course, when the Ainu and Okinawa are used as ‘pseudo-hybrids’ the corresponding posteriors were not used in the regression analysis described below.

Calculating drift

The LEA software also allowed us to estimate genetic drift since the admixture event in the three populations, through the parameterst_i=T/N_i, whereicorresponds to1(Jomon parental population),2(Yayoi parental population) orh(Japanese hybrid population). Populations that had developed agriculture earlier would have increased in size earlier and would thus exhibit lower amounts of drift since the admixture event. Consequently, if the admixture model is consistent, thet₁values should in general be higher than thet₂values, whereast_hvalues should be more variable across populations.

Spatial variation of admixture: regression analysis

To detect, quantify and assess the significance of any geographical trend in the admixture proportions across Japan, we used a linear regression approach similar to that used by Chikhiet al.²The idea is to determine whether there is a correlation between the ‘Yayoi contribution,’ measured byp₂, and the geographical distance from the population used forP₂. For each location sampled in Japan, we computed a geographical distance from the sample used asP₂and then estimated a linear regression between this distance andp₂. To account for the uncertainty aroundp₂, we followed the resampling approach used by Chikhiet al.²For each of the Japanese samples, onep₂value was randomly sampled from the corresponding posterior distribution. This process was repeated 1000 times to obtain the empirical distribution of regression lines. This was done independently, for each set of admixture analyses performed, using a particular pair of parental populations. A similar approach was used fort_h, to determine whether drift in the admixed populations was also correlated with the geographical distance. The geographical distance was calculated as a straight line from the central point of the area corresponding to the population used asP₂(for example, close to Seoul for Korea), taking into account an entering in Japan from Korea, through Kyushu. It is worth noting that the spatial points used are necessarily non-independent, because of local gene flow, as was, for instance, noted by Sokalet al.³³As a consequence, allele or haplotype frequencies are spatially autocorrelated, and hence violate the assumption of independence necessary to calculate the significance of a linear regression, using classical approaches. This is why we did not perform such tests, and used the values to represent the relationship between the parameter of interest and geographical distance.

F_STanalysis

The genetic structure of the populations was also described usingF_STvalues. These values were computed with the equationF_ST=(_T−_S)/_T,³⁴using the Okinawans and Ainu against all other Asian populations.

Results

Admixture proportions

Figure 2shows, for the different pairs of parentals tested, the mode of thep₁posterior distributions (represented inFigure 3a), wherep₁represents the HG contribution to the Japanese populations. The modes represent the most probable values, but as thep₁distributions are wide they should be interpreted with caution. In fact, earlier simulation results suggest that any point estimate (median, mean or mode) should be interpreted with care²⁰and we should rather focus on spatial trends, if any, across the populations.²Although it was difficult to infer precise contributions, this figure shows a clear difference in the results, whether Asian or non-Asian populations are used asP₂: (i) with Sardinia or Oceania, thep₁modes are equal to one, indicating no contribution of these populations to the Japanese populations and (ii) with Asian populations, the modalp₁values are widespread and are most of the time higher than 0.2, suggesting a variable but significant contribution of the HG. There is one exception, when we use the CAS population asP₂, thep₁modes for the Ainu are very low, which is counterintuitive. A closer look at these posteriors shows that they are very flat and hence that there is little information to inferp₁.

Even though there is an uncertainty onp₁values for specific populations (seeFigure 3a), a geographical trend in the ‘Yayoi’ contribution (p₂) is found (Figure 3b). The randomization approach applied to test this trend is summarized inTable 1(补充图1 SI), correlation coefficients obtained from the linear regressions are represented. Although the values are small, this figure again shows a clear signal. First, when the Sardinian or Oceanian samples are used asP₂, which are unlikely to have contributed to the Japanese gene pool, there is a positive or no correlation, suggesting that the contribution of these two populations increases (or stays constant) with geographical distance from their current location. Second, in all the other analyses (that is, when theP₂populations are Korea, NEA, SEA and CAS), the correlation coefficients values are negative; that is,p₂logically decreases with geographical distance from the population used asP₂. In addition, with the exception of SEA asP₂, when comparing the analyses of Ainu vs Okinawa asP₁, the most negative correlation values seem to be associated with the Ainu. This trend is, however, only close to significance (P=0.058).

Drift

InFigures 4a and b, the amount of drift between the present-day samples of the populations used asP₁orP₂, and the ancestral populations, is represented throught₁andt₂respectively. Comparing thet₁andt₂distributions clearly indicates that the two estimates are extremely different. In fact, in all the analyses done, thet₁modal values, of all the Japanese populations studied, were always greater thant₂and, at the same time, were higher if we used the Ainu sample asP₁, instead of Okinawa (data not shown). Moreover, thet₁posteriors were wider and had more variable modal values thant₂, especially when using the Ainu asP₁, but nevertheless were similar to each other.

When we analyzed thet_hestimates, they appeared to be highly variable (data not shown) and to display a geographical trend. As forp₂, we applied a regression of these estimates against the geographical distance, from the population used asP₂, to all the possible combinations of parentals. The correlation coefficients obtained from these linear regressions are represented inTable 1(and Supplementary Figure SI 1b). All the correlation coefficients are positive (except with Ainu vs Sardinia), showing that thet_hvalues increase as geographical distance increases, but at the same time are much smaller when we use the Ainu asP₁, on the order of 0.034–0.067 vs the 0.486–0.666 for the Okinawans.

F_ST

Table 2(see also Supplementary Figure SI 2) shows higherF_STvalues in the pairwise comparisons involving the Ainu, compared with those involving the Okinawans. It also shows, that for both sets of pairwise comparisons, theF_STvalues increase when the geographic distance increases in a southwestern direction from the northern tip of Honshu (Aomori) toward Kyushu. This is particularly interesting as the Ainu and Okinawan are located on opposite sides of the Japanese archipelago and hence of this axis. We also note that this trend ofF_STvalues shows a clear and sudden increase when the samples are taken from the Asian continent, starting with Korea. This is particularly clear despite the fact that, when we consider only the Japanese populations, theF_STvalues involving Okinawans and Ainu are on different scales, the first set of values all being below 0.1 and the others all being above 0.1. In fact, theF_STbetween the Ainu and Okinawans (0.096) is the smallestF_STvalue among the pairs involving the Ainu (which vary from 0.096 to 0.219 in Kyushu), but it is the largest among the pairs involving the Okinawans (which vary from 0.018 in Aomori to 0.096 against the Ainu). Another consequence of this cline is that the Okinawan population seems genetically close to the northernmost populations of Honshu, but strangely not to the Ainu, the only sampled population north of Honshu. Indeed, theF_STwith the latter is higher than theF_ST冲绳人之间的值和其他日本订单开云体彩官网pulations. Thus, the Okinawans appear to be the Japanese who are genetically closest to the Ainu from the Ainu viewpoint, whereas it is exactly the opposite from the Okinawans viewpoint. It is as if the Okinawans were ‘virtually’ located in northern Honshu, and the Ainu were genetically close to them but had been submitted to significant drift.

Topof page

Discussion

Dual origins of Japanese

Although our results may only reflect the paternal history of the Japanese, they confirm the idea that a significant admixture took place and thus do not support either the replacement or the transformation models between the incoming Yayoi and the local Jomon. Indeed, in the replacement model the estimate ofp₁should be equal to zero (or at least very close, because of statistical uncertainty), whereas in the transformation modelp₁would be close or equal to 1. This is clearly not what we observe. Moreover, our results show a decreasing geographical trend in the Yayoi contribution across Japan, when populations are sampled in a southwest–northeast direction (Figures 3a and b). These results agree with a model in which the first farmers entered in Japan from Korea, through the closest island (Kyushu), and then spread across most of Japan moving to the northeast (until the geographical limits of the Honshu Island). During the expansion of farmers, it is expected that the rise in population density, because of food production, should lead to a more limited drift (as populations were larger). This can be seen in the gradient observed with thet_hestimates, which suggests that drift is higher in the northernmost populations (with a maximum observed in the Ainu), where the archeological record suggests a later arrival of agriculture.¹¹In addition, the differences encountered betweent₁andt₂estimates (t₁t₂) are consistent with a model of an expanding population that dispersed to a less populated area; that is, theP₂populations (Yayoi) increased in size earlier in time, having suffered a lower amount of drift. Interestingly, thet₁values are higher than those found in Europe,²which could be because of an earlier introduction of agriculture in Europe.

During the admixture process, it is important to note that the indigenous populations, the Jomon, who admixed with the Yayoi, were probably genetically differentiated from each other across the Japanese islands. How differentiated they were, some 2000 years ago, is difficult to say, but this pre-admixture differentiation should have some implication in the analyses, and their interpretation. The fact that allt₁posteriors were very similar to each other suggests that even if there was differentiation between HG populations, before the arrival of the Yayoi, they were not dramatically different, compared with the amount of drift that occurred since the admixture event. This indirectly shows that our model, despite its simplicity, captures important aspects of the ‘Neolithic transition’ in the Japanese archipelago. However, we note that, depending on whether we use the Ainu or Okinawans asP₁, thet₁estimates are rather different, suggesting a much higher drift when the Ainu are used. This, together with the fact that the Ainu have much higherF_STthan the Okinawans against all other populations, suggests that the Ainu have probably had a much lower effective size than the Okinawans. It could be because of a greater isolation, a later and more limited influence of agriculture or a combination of both.

The continental origin of the Yayoi farmers

Several hypotheses have been suggested regarding the geographic origin of the Asian populations, which gave rise to the Yayoi, even though it is usually accepted that they probably entered Japan through South Korea. Nevertheless, skull and teeth morphology inference⁷and classical markers^13,16support an NEA origin.

More recently, Hammeret al.¹⁵placed the Yayoi farmers as having originated in SEA. However, with the same data as Hammeret al.¹⁵(with the exception of the Korean data), the admixture model we used could not establish with so much accuracy the continental origins of these populations, but our results suggest that they entered Japan through the Korean Peninsula. It may be important to note that our approach is model-based and has been tested on simulated data,²⁰whereas the conclusion reached by Hammeret al.¹⁵was based on visual patterns of allele or haplogroup frequencies and were neither justified by any statistical test, nor by analyses of simulated data. Thus, our approach is not ‘just’ confirming established results, but rather adding more solid results to conclusions whose statistical validity was not determined. In addition, when we found that the exact location of the Yayoi cannot be ascertained with certainty, whereas Hammeret al.¹⁵assert that they arrived from SEA, one should question the strength of the latter statement. This should not be taken as a criticism of the Hammeret al.¹⁵study, which provided both new results and hypotheses to test. Rather, what our results show is that it might be possible, using an admixture approach, to test different hypotheses, something that has not been done so far. Indeed, our method was able to identify populations that clearly could not have contributed to the modern Japanese gene pool at that time (namely Sardinia and Oceania). If some of the Asian parental populations that were analyzed had generated results similar to those of Sardinia and Oceania, they could have been identified as unlikely parentals. This type of results was not observed, which suggests that these data do not contain enough information to clearly identify the most likely descendent of one of the parental populations of the modern Japanese, namely the Yayoi. This is not necessarily surprising as the Y-chromosome represents only one set of linked markers. We believe thus that until morelociare obtained this question may not be easily answered, and should remain open. If we have contributed to making this statement, we feel that a significant step will have been taken.

In summary, our results support at least one admixture event in the peopling of Japan, namely the spread of Yayoi farmers by a process of demic diffusion, similar to the one in Europe during the Neolithic age^2,23,35We suggest that when the Yayoi men entered Japan, and brought with them agriculture and new technologies, they also raised the carrying capacity of the area first colonized, leading to an increase in the size of the newly admixed populations. When this area could no longer support the increased population, their descendents expanded into new territories, repeating the admixture process. By the time the geographic limits of Japan were reached (northeastward until the Hokkaido Island and southwestward in the Ryukyu Islands), there was a gradual dilution of the Yayoi’s gene pool. However, in spite of having detected the presence of Jomon and Yayoi contributions in the Japanese populations, the method we used was not capable of locating precisely the area of origin of the ancestral populations, and different populations seemed to produce similarly consistent results. Nevertheless, the general approach appears to provide interesting and promising results, which should open new avenues for research.

Topof page

References

1. Bellwood, P.First Farmers: The Origins of Agricultural Societies(Blackwell Publishing, Oxford, 2005).
2. Chikhi, L., Nichols, R. A., Barbujani, G. & Beaumont, M. A. Y genetic data support the Neolithic demic diffusion model.Proc. Natl Acad. Sci. USA99, 11008–11013 (2002). |Article|PubMed|ChemPort|
3. Currat, M. & Excoffier, L. The effect of the Neolithic expansion on European molecular diversity.Proc. Biol. Sci.272, 679–688 (2005). |Article|PubMed|ISI|
4. Dupanloup, I., Bertorelle, G., Chikhi, L. & Barbujani, G. Estimating the impact of prehistoric admixture on the genome of Europeans.Mol. Biol. Evol.21, 1361–1372 (2004). |Article|PubMed|ChemPort|
5. Cavalli-Sforza, L. L., Mennozi, P. & Piazza, A.The History and Geography of Human Genes(Princeton University Press, Princeton, 1994).
6. Mizoguchi, Y. inPrehistoric Hunter-Gatherers in Japan(eds. Akazawa, T. and Aikens, C.M.) 107–136 (University of Tokyo Press, Tokyo, 1986).
7. Hanihara, K. Dual structure model for the population history of the Japanese.Jpn. Rev.2, 1–33 (1991).
8. Matsumura, H. Differentials of Yayoi immigration to Japan as derived from dental metrics.Homo52, 135–156 (2001). |Article|PubMed|ChemPort|
9. Pope, K. O. & Terrell, J. E. Environmental setting of human migrations in the circum-Pacific region.J. Biogeogr.35, 1–21 (2008).
10. Ono, A., Sato, H., Tsutsumi, T. & Kudo, Y. Radiocarbon Dates and Archaeology of the Late Pleistocene in the Japanese Islands.Radiocarbon44, 447–494 (2002).
11. Highman, C. inThe Human Past: World Prehistory and the Development of Human Societies(ed. Scarre, C.) 234–263 (Thames and Hudson, London, 2005).
12. Jin, H., Kwak, K., Hammer, M. F., Nakahori, Y., Shinka, T., Lee, J.et al. Y-chromosomal DNA haplogroups and their implications for the dual origins of the Koreans.Hum. Genet.114, 27–35 (2003). |Article|PubMed|ChemPort|
13. Sokal, R. R. & Thomson, B. A. Spatial genetic structure of human populations in Japan.Hum. Biol.70, 1–22 (1998). |PubMed|ChemPort|
14. Hammer, M. F. & Horai, S. Y chromosomal DNA variation and the peopling of Japan.Am. J. Hum. Genet.56, 951–962 (1995). |PubMed|ISI|ChemPort|
15. Hammer, M. F., Karafet, T. M., Park, H., Omoto, K., Harihara, S., Stoneking, M.et al. Dual origins of the Japanese: common ground for hunter-gatherer and farmer Y chromosomes.J. Hum. Genet.51, 47–58 (2006). |Article|PubMed|
16. Omoto, K. & Saitou, N. Genetic origins of the Japanese: a partial support for the dual structure hypothesis.Am. J. Phys. Anthropol.102, 437–446 (1997). |Article|PubMed|ChemPort|
17. Horai, S., Murayama, K., Hayasaka, K., Matsubayashi, S., Hattori, Y., Fucharoen, G.et al. mtDNA polymorphism in East Asian Populations, with special reference to the peopling of Japan.Am. J. Hum. Genet.59, 579–590 (1996). |PubMed|ChemPort|
18. Horai, S., Kondo, R., Murayama, K., Hayashi, S., Koike, H. & Nakai, N. Phylogenetic affiliation of ancient and contemporary humans inferred from mitochondrial DNA.Philos. Trans. R. Soc. Lond., B, Biol. Sci.333, 409–417 (1991). |Article|PubMed|ChemPort|
19. Oota, H., Saitou, N., Matsushita, T. & Ueda, S. A genetic study of 2,000-year-old human remains from Japan using mitochondrial DNA sequences.Am. J. Phys. Anthropol.98, 133–145 (1995). |Article|PubMed|ChemPort|
20. Chikhi, L., Bruford, M. W. & Beaumont, M. A. Estimation of admixture proportions: a likelihood-based approach using Markov chain Monte Carlo.Genetics158, 1347–1362 (2001). |PubMed|ISI|ChemPort|
21. 兰德里,美女,e . m . S p & Barbujani gs and evolution of the Europeans’ genome: evidence from multiple microsatellite loci.Proc. Biol. Sci.273, 1595–1602 (2006). |Article|PubMed|ISI|ChemPort|
22. Chikhi, L. inExamining the Farming/Language Dispersal Hypothesis(eds. Bellwood, P. & Renfrew, C.) 435–447 (McDonald Institute for Archaeological Research, Cambridge, 2003).
23. Ammerman, A. J. & Cavalli-Sforza, L. L.The Neolithic Transition and the Genetics of Populations in Europe(Princeton Univ. Press, Princeton, 1984).
24. Tajima, A., Hayami, M., Tokunaga, K., Juji, T., Matsuo, M., Marzuki, S.et al. Genetic origins of the Ainu inferred from combined DNA analyses of maternal and paternal lineages.J. Hum. Genet.49, 187–193 (2004). |Article|PubMed|ISI|ChemPort|
25. Xue, Y., Zerjal, T., Bao, W., Zhu, S., Shu, Q., Xu, J.et al. Male demography in East Asia: a north-south contrast in human population expansion times.Genetics172, 2431–2439 (2006). |Article|PubMed|ChemPort|
26. Semino, O., Passarino, G., Oefner, P. J., Lin, A. A., Arbuzova, S., Beckman, L. E.et al. The genetic legacy of PaleolithicHomo sapiens sapiensin extant Europeans: a Y chromosome perspective.Science290, 1155–1159 (2000). |Article|PubMed|ISI|ChemPort|
27. YCC. A nomenclature system for the tree of human Y-chromosomal binary haplogroups.Genome Res.12, 339–348 (2002). |Article|PubMed|ISI|
28. Langella, O., Chikhi, L. & Beaumont, M. LEA (likelihood-based estimation of admixture): a program to estimate simultaneously admixture and time since the admixture event.Mol. Eco. Notes1, 357–358 (2001). |Article|
29. Hudson, M. The Linguistic Prehistory of Japan: some archaeological speculations.Anthropol. Sci.102, 231–255 (1994).
30. Toma, T., Tsunehiko, H., Hajime, S., Kuniaki, H. & Hajime, I. Metric dental diversity of Ryukyu Islanders: a comparative study among Ryukyu and other Asian populations.Anthropol. Sci.115, 119–131 (2007). |Article|
31. Haneji, K., Hanihara, T., Sunakawa, H., Toma, T. & Ishida, H. Non-metric dental variation of Sakishima Islanders, Okinawa, Japan: a comparative study among Sakishima and neighboring populations.Anthropol. Sci.115, 35–45 (2007). |Article|
32. Yamaguchi-Kabata, Y., Nakazono, K., Takahashi, A., Saito, S., Hosono, N., Kubo, M.et al. Japanese population structure, based on SNP genotypes from 7003 individuals compared to other ethnic groups: effects on population-based association studies.Am. J. Hum. Genet.83, 445–456 (2008). |Article|PubMed|ChemPort|
33. Sokal, R. R., Harding, R. M. & Oden, N. L. Spatial patterns of human gene frequencies in Europe.Am. J. Phys. Anthropol.80, 267–294 (1989). |Article|PubMed|ChemPort|
34. Nei, M. F-statistics and analysis of gene diversity in subdivided populations.Ann. Hum. Genet.41, 225–233 (1977). |Article|PubMed|ISI|ChemPort|
35. Barbujani, G., Sokal, R. R. & Oden, N. L. Indo-European origins: a computer-simulation test of five hypotheses.Am. J. Phys. Anthropol.96, 109–132 (1995). |Article|PubMed|ChemPort|