Journal of Human Genetics(2009)54,349–354; doi:10.1038/jhg.2009.39; published online 8 May 2009 http://www.nature.com/jhg/journal/v54/n6/full/jhg200939a.html
- 1Instituto Gulbenkian de Ciência, Rua da Quinta Grande 6, Oeiras, Portugal
- 2UMR CNRS 5174 Evolution and Diversité Biologique, Université Paul Sabatier, 118 route de Narbonne, Bât. 4R3 b2, Toulouse cedex 4, France
Correspondence: R Rasteiro, Instituto Gulbenkian de Ciência, Rua da Quinta Grande 6, 2780-156 Oeiras, Portugal. E-mail:email@example.com
Received 8 January 2009; Revised 6 April 2009; Accepted 8 April 2009; Published online 8 May 2009.
The first inhabitants of Japan, the Jomon hunter-gatherers, had their culture significantly modified by that of the Yayoi farmers, who arrived at a later stage from mainland Asia. How this change took place is still debated, but it has been suggested that modern Japanese are the product of an admixture between these two populations. Here, we applied for the first time an admixture approach to study the Jomon–Yayoi transition, using Y-chromosomal data published earlier. Our results suggest that the Neolithic transition, in this part of the world, probably took place by a process of demic diffusion. We also show that for two populations that could not have contributed to this process, our approach is able to detect inconsistencies when they are used as parental populations. However, despite these promising results, we could not locate precisely the geographical origin of the Yayoi in mainland Asia, as different potential sources gave similarly good results. This suggests that more loci would be required for a better understanding of the peopling of Japan.
The development and spread of farming, referred to as the ‘Neolithic transition,’ was one of the major demographic events of human prehistory.1This 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,9and later gave rise to the Jomon culture (12000 BP).10Although 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,12and then spread northeastward, starting the Yayoi period.
Map of the Japanese Islands. Approximate geographical locations of the Japanese populations analyzed in this study. The other samples used as parental are not represented on the map
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,5the former entails a movement of the Yayoi culture and ideas rather than people, with consequently no genetic contribution of the Yayoi to modern Japanese.6然而,现实一定是那么极端和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 Hanihara7and Matsumura8on the basis of dental and cranial characteristics, and more recently by a number of authors who used genetic data,13,14,15,16,17including 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 approach20that has already been applied to address the Neolithic transition in Europe2,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.2We 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.23
Materials and methods
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.,15except the Ainu data that were pooled with data from Tajimaet al.24Mainland Asian data15were obtained for populations from Northeast (441), Southeast (683) and Central (419) Asia and also a sample from Korea25(43). We also used two additional populations, Sardinia26(77) and Oceania15(209),作为掺合料莫父母的人群del used (see below). Y-chromosome binary haplogroups were defined by the analysis of the binary polymorphisms described in Hammeret al.15The Y-chromosome lineages from Japan, mainland Asia, Korea, Oceania and Sardinia followed the haplogroups nomenclature of the Y Chromosome Consortium.27
The admixture model
The admixture method used assumes that an ‘admixed’ or ‘hybrid’ population (H), of sizeNh, is the result of the admixture of two independent parental populations,P1andP2, of sizeN1andN2,Tgenerations ago, with respective contributionsp1andp2(p2=1−p1). 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 differentNi(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.,20implemented in the LEA (Likelihood-based Estimation of Admixture) software freely available atwww.igc.gulbenkian.pt/static/docs/LEA.zip.28LEA 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, theP1population was always used to represent the HG or Jomon, whereas the populationP2was used to represent the farmers of the Yayoi period. Hence, the parameterp1represents the ‘Jomon’ contribution, at the moment of admixture, whereasp2would 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 century29and 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.32For these reasons, the admixture analyses were performed using either the Ainu or the Ryukyuans, the latter represented by the Okinawa sample, as descendents of theP1population, in the different analyses. For the descendents of the Yayoi (considered theP2population), 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 asP2two 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 forP1(Ainu and Okinawa), six populations forP2(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 asP1, asH, against the six differentP2populations. This kind of analysis allowed us to quantify the uncertainty around the estimation ofp1, as the hybrid and one parental population (hereP1) are exactly identical. Thus, thep1后验应该总是有一个模式或相等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.
The LEA software also allowed us to estimate genetic drift since the admixture event in the three populations, through the parametersti=T/Ni, 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, thet1values should in general be higher than thet2values, whereasthvalues 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.2The idea is to determine whether there is a correlation between the ‘Yayoi contribution,’ measured byp2, and the geographical distance from the population used forP2. For each location sampled in Japan, we computed a geographical distance from the sample used asP2and then estimated a linear regression between this distance andp2. To account for the uncertainty aroundp2, we followed the resampling approach used by Chikhiet al.2For each of the Japanese samples, onep2value 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 forth, 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 asP2(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.33As 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.
The genetic structure of the populations was also described usingFSTvalues. These values were computed with the equationFST=(T−S)/T,34using the Okinawans and Ainu against all other Asian populations.
Fig. 2 Jomon contribution, across Japan. Mode of thep1后验分布,所有日本人不爱你开云体彩官网lations analyzed, withp1representing the hunter-gatherers Jomon contribution to modern Japanese. In the x axis are represented the parental populations used (P1followed byP2). The letter codes are as follows: AINU—Ainu, OKIN—Okinawa, KOR—Korea, NEA—northeast Asia, SEA—southeast Asia, CAS—central Asia, OCEA—Oceania and SARD—Sardinia
Figure 2shows, for the different pairs of parentals tested, the mode of thep1posterior distributions (represented inFigure 3a), wherep1represents the HG contribution to the Japanese populations. The modes represent the most probable values, but as thep1distributions 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 care20and we should rather focus on spatial trends, if any, across the populations.2Although it was difficult to infer precise contributions, this figure shows a clear difference in the results, whether Asian or non-Asian populations are used asP2: (i) with Sardinia or Oceania, thep1modes are equal to one, indicating no contribution of these populations to the Japanese populations and (ii) with Asian populations, the modalp1values 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 asP2, thep1modes 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 inferp1.
Jomon and Yayoi contributions, across Japan. (a的后验分布p1for all Japanese populations (Kyushu—gray, Tokushima—dashed gray, Shizuoka—black, Aomori—dashed black and Okinawa—dotted black) and parental populations (thinner lines: Ainu—dark gray and Korea—dashed dark gray) used. Each curve corresponds to the analysis of a specific hybrid population. (b) Linear regression ofp2against geographical distance fromP2(Korea). The circles represent the mean value for each population. These analyses were carried out using the Ainu and Korean populations asP1andP2, respectively.
Even though there is an uncertainty onp1values for specific populations (seeFigure 3a), a geographical trend in the ‘Yayoi’ contribution (p2) 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 asP2, 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 theP2populations are Korea, NEA, SEA and CAS), the correlation coefficients values are negative; that is,p2logically decreases with geographical distance from the population used asP2. In addition, with the exception of SEA asP2, when comparing the analyses of Ainu vs Okinawa asP1, the most negative correlation values seem to be associated with the Ainu. This trend is, however, only close to significance (P=0.058).
InFigures 4a and b, the amount of drift between the present-day samples of the populations used asP1orP2, and the ancestral populations, is represented throught1andt2respectively. Comparing thet1andt2distributions clearly indicates that the two estimates are extremely different. In fact, in all the analyses done, thet1modal values, of all the Japanese populations studied, were always greater thant2and, at the same time, were higher if we used the Ainu sample asP1, instead of Okinawa (data not shown). Moreover, thet1posteriors were wider and had more variable modal values thant2, especially when using the Ainu asP1, but nevertheless were similar to each other.
Distributions of theti‘s for all Japanese populations. (a的后验分布t1. The different curves represent the amount of genetic drift, since the admixture event, between the present sample of Ainu and the ancestral populations of HG (Jomon) that interbred with the incoming farmers (Yayoi). (b的后验分布t2. As ina, but for the drift between the Korean and Yayoi populations instead. The color codes are as inFigure 3.
When we analyzed thethestimates, they appeared to be highly variable (data not shown) and to display a geographical trend. As forp2, we applied a regression of these estimates against the geographical distance, from the population used asP2, 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 thethvalues increase as geographical distance increases, but at the same time are much smaller when we use the Ainu asP1, on the order of 0.034–0.067 vs the 0.486–0.666 for the Okinawans.
Table 2(see also Supplementary Figure SI 2) shows higherFSTvalues in the pairwise comparisons involving the Ainu, compared with those involving the Okinawans. It also shows, that for both sets of pairwise comparisons, theFSTvalues 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 ofFSTvalues 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, theFSTvalues 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, theFSTbetween the Ainu and Okinawans (0.096) is the smallestFSTvalue 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, theFSTwith the latter is higher than theFST冲绳人之间的值和其他日本订单开云体彩官网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.
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 ofp1should be equal to zero (or at least very close, because of statistical uncertainty), whereas in the transformation modelp1would 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 thethestimates, 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.11In addition, the differences encountered betweent1andt2estimates (t1t2) are consistent with a model of an expanding population that dispersed to a less populated area; that is, theP2populations (Yayoi) increased in size earlier in time, having suffered a lower amount of drift. Interestingly, thet1values are higher than those found in Europe,2which 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 allt1posteriors 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 asP1, thet1estimates are rather different, suggesting a much higher drift when the Ainu are used. This, together with the fact that the Ainu have much higherFSTthan 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 inference7and classical markers13,16support an NEA origin.
More recently, Hammeret al.15placed the Yayoi farmers as having originated in SEA. However, with the same data as Hammeret al.15(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,20whereas the conclusion reached by Hammeret al.15was 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.15assert 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.15study, 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 age2,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.
- Bellwood, P.First Farmers: The Origins of Agricultural Societies(Blackwell Publishing, Oxford, 2005).
- 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|
- Currat, M. & Excoffier, L. The effect of the Neolithic expansion on European molecular diversity.Proc. Biol. Sci.272, 679–688 (2005). |Article|PubMed|ISI|
- 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|
- Cavalli-Sforza, L. L., Mennozi, P. & Piazza, A.The History and Geography of Human Genes(Princeton University Press, Princeton, 1994).
- Mizoguchi, Y. inPrehistoric Hunter-Gatherers in Japan(eds. Akazawa, T. and Aikens, C.M.) 107–136 (University of Tokyo Press, Tokyo, 1986).
- Hanihara, K. Dual structure model for the population history of the Japanese.Jpn. Rev.2, 1–33 (1991).
- Matsumura, H. Differentials of Yayoi immigration to Japan as derived from dental metrics.Homo52, 135–156 (2001). |Article|PubMed|ChemPort|
- Pope, K. O. & Terrell, J. E. Environmental setting of human migrations in the circum-Pacific region.J. Biogeogr.35, 1–21 (2008).
- Ono, A., Sato, H., Tsutsumi, T. & Kudo, Y. Radiocarbon Dates and Archaeology of the Late Pleistocene in the Japanese Islands.Radiocarbon44, 447–494 (2002).
- Highman, C. inThe Human Past: World Prehistory and the Development of Human Societies(ed. Scarre, C.) 234–263 (Thames and Hudson, London, 2005).
- 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|
- Sokal, R. R. & Thomson, B. A. Spatial genetic structure of human populations in Japan.Hum. Biol.70, 1–22 (1998). |PubMed|ChemPort|
- 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|
- 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|
- 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|
- 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|
- 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|
- 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|
- 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|
- 兰德里,美女,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|
- Chikhi, L. inExamining the Farming/Language Dispersal Hypothesis(eds. Bellwood, P. & Renfrew, C.) 435–447 (McDonald Institute for Archaeological Research, Cambridge, 2003).
- Ammerman, A. J. & Cavalli-Sforza, L. L.The Neolithic Transition and the Genetics of Populations in Europe(Princeton Univ. Press, Princeton, 1984).
- 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|
- 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|
- 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|
- YCC. A nomenclature system for the tree of human Y-chromosomal binary haplogroups.Genome Res.12, 339–348 (2002). |Article|PubMed|ISI|
- 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|
- Hudson, M. The Linguistic Prehistory of Japan: some archaeological speculations.Anthropol. Sci.102, 231–255 (1994).
- 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|
- 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|
- 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|
- 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|
- Nei, M. F-statistics and analysis of gene diversity in subdivided populations.Ann. Hum. Genet.41, 225–233 (1977). |Article|PubMed|ISI|ChemPort|
- 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|