Gene Flow between the Korean Peninsula and Its Neighboring Countries

Abstract

SNP markers provide the primary data for population structure analysis. In this study, we employed whole-genome autosomal SNPs as a marker set (54,836 SNP markers) and tested their possible effects on genetic ancestry using 320 subjects covering 24 regional groups including Northern ( = 16) and Southern ( = 3) Asians, Amerindians ( = 1), and four HapMap populations (YRI, CEU, JPT, and CHB). Additionally, we evaluated the effectiveness and robustness of 50K autosomal SNPs with various clustering methods, along with their dependencies on recombination hotspots (RH), linkage disequilibrium (LD), missing calls and regional specific markers. The RH- and LD-free multi-dimensional scaling (MDS) method showed a broad picture of human migration from Africa to North-East Asia on our genome map, supporting results from previous haploid DNA studies. Of the Asian groups, the East Asian group showed greater differentiation than the Northern and Southern Asian groups with respect to Fst statistics. By extension, the analysis of monomorphic markers implied that nine out of ten historical regions in South Korea, and Tokyo in Japan, showed signs of genetic drift caused by the later settlement of East Asia (South Korea, Japan and China), while Gyeongju in South East Korea showed signs of the earliest settlement in East Asia. In the genome map, the gene flow to the Korean Peninsula from its neighboring countries indicated that some genetic signals from Northern populations such as the Siberians and Mongolians still remain in the South East and West regions, while few signals remain from the early Southern lineages.

About the Authors:

Jongsun Jung

Contributed equally to this work with: Jongsun Jung, Hoyoung Kang

Affiliations Syntekabio Inc., Seoul, Korea, DNA link Inc., Seoul, Korea, Division of Structural and Functional Genomics, Center for Genome Science, Korea National Institute of Health, Seoul, Korea

Hoyoung Kang

Contributed equally to this work with: Jongsun Jung, Hoyoung Kang

Affiliations Korean Bioinformation Center, KRIBB, Daejeon, Korea, University of Science and Technology, Daejeon, Korea

Yoon Shin Cho

Affiliation: Division of Structural and Functional Genomics, Center for Genome Science, Korea National Institute of Health, Seoul, Korea

Ji Hee Oh

Affiliation: Division of Structural and Functional Genomics, Center for Genome Science, Korea National Institute of Health, Seoul, Korea

Min Hyung Ryu

Affiliation: Division of Structural and Functional Genomics, Center for Genome Science, Korea National Institute of Health, Seoul, Korea

Hye Won Chung

Affiliation: Department of Obstetrics and Gynecology, College of Medicine, Ewha Womans University, Seoul, Korea

Jeong-Sun Seo

Affiliation: Department of Biochemistry and Molecular Biology, Seoul National University of Medicine, Seoul, Korea

Jong-Eun Lee

Affiliation: DNA link Inc., Seoul, Korea

Bermseok Oh

Affiliation: Department of Biomedical Engineering, Kyung Hee University College of Medicine, Seoul, Korea

Jong Bhak

Affiliation: Korean Bioinformation Center, KRIBB, Daejeon, Korea

Hyung-Lae Kim

* E-mail: [email protected]

Affiliation: Division of Structural and Functional Genomics, Center for Genome Science, Korea National Institute of Health, Seoul, Korea

Introduction

The X and Y chromosomes and mitochondrial DNA (mtDNA) recombine to a much lesser extent than do autosomes [1]. This limited, or non-existent, recombination in haploid SNPs (single nucleotide polymorphisms) has been used to trace the evolutionary lineages of individuals and populations showing either maternal or paternal transmission. However, obstacles encountered with the use of lesser (or non-) recombining SNPs, such as genetic drift, low polymorphic nature, natural selection, and the small size of these non-recombining DNAs, limit the scope of genetic studies. Recently, an advanced technology, termed whole genome sampling analysis, has enabled the genotyping of nearly one megabyte of diploid autosomal SNPs at a time [2]. These autosomal SNPs are a reservoir for human evolutionary history, nearly 100 times larger than that of lesser (non-) recombining SNPs, and include the entire range from highly conserved to highly polymorphic SNPs. However, when individuals and populations of different races and regions are compared using the whole genome approach, a large proportion of the redundant and recombining regions exacerbate the difficulties encountered in exploring evolutionary relationships, common traits and disease susceptibilities between populations.

The origin of Amerindians, based on Y chromosome and mtDNA studies, revealed that they migrated from Siberia to Alaska [3]. Recently, over 600 autosomal short tandem repeat (STR) markers were used to delineate the migration routes of most of the Amerindians [4], while analysis of molecular variance, inferred from metric distances between DNA haplotypes of human mitochondrial DNA, was used for migration studies of Southern and Northern Europe [5]. Some diabetes-related rare markers were also found in Amerindian and Finnish families [6]. Y chromosomal DNA haplogroups disclosed the dual origins of the Korean population [7], while the origins of East Asians have been studied using certain autosomal loci [8]. Investigation into male demography in East Asia using Y chromosome haplogroups in a simulation analysis provided information on the settlement times of Northern and Southern populations [9]. Recently, genome-wide analysis of haplotypes in Asia was reported that more than 90% of East Asian haplotypes could be found in either Southeast Asian or Central-South Asian populations and showed clinal structure with haplotype diversity decreasing from south to north [10].

Given a symmetric distance (or a covariance) matrix between individuals, the dimension-reduction methodologies, such as multidimensional scaling (MDS) and principal components analysis (PCA), can display individual interacting patterns effectively by means of two, or more, major eigenvector components [11]. One of the relationship matrices for the comparison of individual whole genomes is based on identity-by-state (IBS) distances that show an opposite relationship to allele-sharing distances (ASD) [12] when applied to a neighbor-joining (NJ)-based phylogenetic tree [13]. Compared to the dimension-reduction methodologies described above, the Markov Chain Monte Carlo (MCMC)-based STRUCTURE program is also one of the most popular methods for understanding genetic structures and admixtures [14]. Recently, Salmela et al. [15] showed that 250K genome-wide autosomal SNPs could provide powerful resources for the study of population structure and differentiation in Northern Europe, including Northern Germany, Great Britain, Finland and Sweden.

Geographically, the Korean Peninsula is a strategic location in East Asia surrounded by China, Japan and Russia. Therefore, the objectives of our gene flow analyses in Northern, Eastern and Southern Asia were to determine the extent of the gene flow from these neighboring countries and the levels of genetic signals that have remained in the Korean Peninsula, using subjects collected from the ten most historical regions in South Korea. In addition, samples from Mongolia, Jinlin in China and Amerindian populations, representing Northern people; Cambodians and Vietnamese, representing Southern people; and Chinese and Japanese were used to measure their influences on the gene flow. For this study, we propose three different gene flow models based on historical events, geographical location and anthropology; Model I: SW (South West) Korea, which includes the JJ, NJ and GJ regions (Fig. 1) has some lineages in common with the BaekJae Empire (BC18-AD660) through a migration event by the Northern people of the Goguryeo Empire which ruled most of Northern Manchuria and Korean peninsula (BC37-AD568), and are also closely related to Mongolians. They ruled most of the Eastern part of China and most of the South Western part of the Korean Peninsula. Model II: SE (South East) Korea, which includes the GU, US and GR populations (Fig. 1), is the region where the Northern people from Siberia (based on their grave patterns and ancient cultural traditions) and the loyal Southern families were settled in the Shilla(BC57-AD918) and Kaya(AD42-532) Empires, respectively. Model III: MW (Middle West) Korea, which includes the YC, JC, CA and PC populations (Fig. 1), is the human melting pot of the Korean Peninsula where Western Chinese, including those from the SanDung Peninsula, crossed the Yellow Sea between China and Korea, living and trading in both regions. In addition, it is known that historically the Northern and Southern Koreans often clashed with each other in this region. The geographical locations of the models used in this study are shown in Fig. 1 and Table 1.

[Figure omitted. See PDF.]

Figure 1. Map of Asia (A) and the Korean Peninsula (B).

Historically and geographically, the Korean Peninsula can be subdivided into two regions, the West and East, based on the continuous ridge-line of mountains from Janggun Peak of Mt. Baekdu on the Northern border of the Korean Peninsula down to Mt. Jiri near its Southern end. The Western part of South Korea is further subdivided into Middle and South regions. The abbreviations for populations, their geographical regions and sample sizes are summarized in Table 1.

https://doi.org/10.1371/journal.pone.0011855.g001

[Figure omitted. See PDF.]

Table 1. Geographic origins and sizes of the samples.

https://doi.org/10.1371/journal.pone.0011855.t001

Marker selection is the most important step in data handling. Because sex chromosomes and mtDNA recombine to a much lesser extent than do autosomes during meiosis, they do not record gene flow event on themselves. Only when complete enumeration or meaningful sampling survey was used, considerable inferences can be deducted. On the other hand, because autosomes accumulate the record of the gene flow event by way of the recombination, the markers on autosomes can be a better object than the ones on sex chromosomes and mtDNA; however, it is also necessary to deal with the effects of recombination hotspots (RH) and linkage disequilibrium (LD). Therefore, in this study, we have evaluated several variables including RH, LD and genotype call rates using known clustering methods such as IBS-based Multi-Dimensional Scaling (MDS), Principal Components Analysis (PCA), the ASD-based Neighbor Joining (NJ) method, and STRUCTURE using 320 subjects from Amerindian populations, from four populations from the HapMap Project (Yoruba, Caucasian, Japanese, and Chinese), and from Vietnamese, Cambodians, North-Eastern Chinese, Mongolians, and Koreans to reconstruct a human migration route that had been partly analyzed using mtDNA and the X and Y chromosomes.

Results

Genome map of human migration

The genome map was drawn using MDS with an RH-free dataset (Supporting Methods and Discussion S1). The genome map in Fig. 2A showed that the Yoruban, Caucasian, and Amerindian populations lined up first, and then Northern and Southern Asian populations were merged in East Asia in Fig. 2B. In the extended genome map for Northern and Southern Asia, the merged pattern became clearer than that of the whole populations where the MH (Mongolian) population lined up to JPT (Tokyo) and JJ (JeJu) populations, while the CB (Cambodian) population lined up to the VN (Vietnamese), VK (Vietnamese-Korean), CHB (Peking), and JL (Jilin) populations.

[Figure omitted. See PDF.]

Figure 2. Genome Map of World (A) and Asia (B).

https://doi.org/10.1371/journal.pone.0011855.g002

In our analysis of genetics structure by STRUCTURE program [16], four populations displayed perfect population membership probabilities, the Africans, Caucasians, Amerindians, and Northern Asians, and three populations displayed partial membership, the Mongolians, Vietnamese, and Cambodians (Fig. 3). Substantial admixture patterns between populations were also observed for Mongolians and Caucasians, and Amerindians with East Asians and Caucasians.

[Figure omitted. See PDF.]

Figure 3. Genetic Structure.

A)World, B)Asia, C)Korea, Japan and China. Note that the plot with a recombination hotspot (RH) of K = 5 on the left is shown to compare with the same plot without the RH effect on the right where the purple bands were spread over all Asian populations.

https://doi.org/10.1371/journal.pone.0011855.g003

Gene flow between the Korean Peninsula and neighboring countries

Genome map and an NJ tree [18], [19] for Korea-China-Japan are shown in Fig. 4, where the genome map is a further extension of the World and Asian map (Fig. 2). In the figure, dots representing the Chinese are located at the bottom and those for Japanese are at the top, while some outliers in SE Korea and Kobe in Japan are shown on the right hand side. Of the ten historical regions in South Korea, some in SW Korea overlap with those of Japan, while most of the MW Korean regions are located at the center of the genome map. Similarly, in the NJ tree, nodes for SW Korea are close to those in Japan, MW Korea is close to China, and SE Korea is located at the right hand side of the tree.

[Figure omitted. See PDF.]

Figure 4. MDS and NJ Tree of Korean, Chinese and Japanese.

A)MDS(multi-dimensional scaling) [17], B)NJ Tree. In the NJ tree, samples for Japan, Kobe-Japan, SW Korea, SE Korea, MW Korea, China, Jinlin, and Korean-Japanese are blue, cyan, green, purple, black, brown and dark green, respectively.

https://doi.org/10.1371/journal.pone.0011855.g004

The gene flow events of the three selected models for SW, MW and SE Korea can be assessed using the genome map. The populations in Model I (SW Korea) are closer to Mongolians than are the other two regions in the genome map (Fig. 2B). Historically, some of the loyal families and their subjects in the Goguryeo Empire moved to this region and formed the BaekJae Empire in BC18-22. This region also showed connections with populations in Tokyo (JPT), as illustrated in Fig. 4. Certain outliers in Model II (SE Korea) display some similarity to the people of Kobe, a port city near Osaka, indicating that there may have been links between the two regions. In addition, considering that the SE Korea region has some connections with Siberian lineages, with respect to grave patterns and culture, it is possible that the outliers in the GU and Kobe (KB) populations could be of Siberian lineage. On the other hand, the GR and US populations showed average signals in the Korean Peninsula. Historically, the Kaya Empire, with its southern lineages, was formed in the GR region and then the Shilla and Kaya Empires became united around AD532. Very recently, the US region became one of the rapidly developing regions, and people from other provinces moved to this region. This might explain why it shows an average signal in South Korea. Model III (MW Korea): the Middle West area formed a melting pot in the Korean Peninsula because populations moving from South to North, North to South, and from Eastern China, including the SanDung peninsula, to the Middle West in Korea all came together in this region. In the genome map, the signals for MW Korea are also close to those for Peking (CHB) in China. The overall result for the Korea-Japan-China genome map indicates that some signals for Mongolia and Siberia remain in SW Korea and SE Korea, respectively, while MW Korea displays an average signal for South Korea.

Genetic drift in the Korean Peninsula

Fst statistics assume that alleles arrive in the population via migration rather than mutation. Therefore, the mean Fst of two populations displays the differentiation caused by a migration event. Table 2 indicates that the population showing the highest level of differentiation in the world (in our dataset) was the Amerindian group, while, in Asia, the East Asian group revealed the highest level of differentiation compared to Northern and Southern Asians, whereas there was no difference among East Asians, including Koreans, Japanese and Chinese. In addition, Salmela et al. [15] suggested that Eastern Finns also showed signs of genetic drift by analyzing the relationships between monomorphic markers. Using the same measurements for East Asian groups (Korea, Japan and China), nine out of ten historical regions in South Korea showed similar signs of genetic drift that might be due to the later settlement of East Asia. Of the nine historical regions, Gimje (GJ) showed the highest genetic drift. On the other hand, Gyeongju (GU) in South East Korea and Jinlin (JL) in North-Eastern China displayed the lowest drift, suggesting that the populations in the GU and JL regions could be earlier settlers of the Korean Peninsula (Table S1 in supporting information).

[Figure omitted. See PDF.]

Table 2. Mean Fst distribution of World, Asia and East Asia.

https://doi.org/10.1371/journal.pone.0011855.t002

Discussion

Current populations in the Korean Peninsula were formed through interactions between populations from Southern and Northern Asia. Therefore, the goal of this study was to determine whether there are still genetic signals remaining in the Korean Peninsula using autosomal genome-wide SNPs, which are 100 times more frequent than those in lesser recombining DNA, such as mtDNA, X and Y chromosomes. Interestingly, some of the historical regions and the regionally separated island, Jeju (JJ) and Gyeongju (GU) in Korea still harbor some old signals from Northern populations, while few signals remain from Southern populations.

In this study, we reconstructed the genome map using data from 320 individuals in 24 groups after correcting and re-evaluating all bias parameters (Supporting Methods and Discussion S1, Fig. S1, S2, S3 and S4). Interestingly, the genome map was similar to the physical map of Europe and Asia, a relationship that was reported previously in a correlation experiment between physical distances and genetic variations by Mountain and Cavalli-Sforza [12]. In addition, genetic structure analysis using the STRUCTURE method revealed the existence of five major populations, African, Caucasian, Amerindian, North-East Asian, and Southern Asian. Therefore, in between, there were significant admixtures such as Mongolian with Caucasian, Vietnamese (or Cambodian) with unknown Southern original settlers, and Amerindians with both North-East Asians and Caucasians (Fig. 2). Thus, the genome map could be used to estimate the relative evolution times between populations and genome interactions (or gene flows) (Fig. 4) such as Korean-Vietnamese, between Vietnam and Korea; Korean-Japanese, between Korea and Japan; Mongolian between Amerindians and North-East Asians; Chinese between Koreans and Vietnamese; and Vietnamese between Cambodians and Chinese. In particular, through both the genome map and STRUCTURE results, the evolutionary relationship between Koreans, Japanese, and Chinese has become clearer since it can be seen that the North-East Chinese (JL) and the Peking Chinese (CHB) are relatively closely related, and the Japanese (JPT) and Korea-JeJu (JJ) are also fairly close (Fig. 2).

In this study, we didn't include North Korean samples because of political reasons. However, historical migration event of BaekJae from Goguryeo Empire (BC37-AD568) in Northern Korea imply that Northern lineages remain in South Korea. Nevertheless, considering the fact that recent migration events between provinces in North Korea are much less than those in South Korea, the gene flow study of North Korean may provide additional clues in the future study.

Materials and Methods

Ethics Statement

DNA samples of Southern and Northern Asian groups for SNP genotyping were obtained from the Biobank of the National Genome Research Institute (BNGRI), and from the Northern East Asian Nation Functional Genome Project. The samples used for this study all provided written informed consent.

Recombination hotspots

HapMap Phase II [21] provides highly dense SNPs compared to most of the 50K Affymatrix chip SNPs used in this study. Therefore, given a marker position in the 50K Affymatrix genotype data, all lower and upper 2.5K SNPs around each of the 50K SNP markers were extracted from HapMap and named the 5K_HapMap dataset. Using the 5K_HapMap dataset, the pairwise R2 coefficients and their means in each of the four populations, CHB, JPT, CEU, and YRI, were calculated at each marker position using the EM algorithm [22]. Therefore, each of the 50K SNP markers has four different mean R2 coefficients, and if at least one of the four populations (CHB, JPT, CEU, and YRI) had an R2 value less than 0.7, the marker position was considered as a possible recombination hotspot (RH). Overall, 5K of 54K SNPs were assigned as recombination hotspots. The extraction of 5K spanning SNPs from HapMap and exhaustive pairwise LD (linkage disequilibrium) calculations for the selected 5K_HapMap dataset were performed on an IBM supercomputer (Model JS20) and a cluster computer (DualPowerPC970) with 100 nodes, and took four working days with most of the nodes.

Generation of different GCR sample sets

The average genotype call rate (GCR) of the current data is 98.7% and named as “data0”. To generate several sample sets having a series of different successful call rates, no call genotypes were inserted randomly in each sample:(1)where i is a sample index, Ngt, Nng, and Nig are the numbers of genotypes, natural no call genotypes, and randomly inserted no call genotypes, respectively, and the sum of Ngt, Nng, and Nig is the total number of SNPs used in this study. In this way, three extra genotype datasets, data1, data2, and data3, were generated as different call rate datasets, GCR (95%), GCR (90%), and GCR (85%), respectively. The marker call rate cutoff of data0 was set to 95% and those of data1, data2 and data3 were set to 80%.

Regional specific markers

The selection of regional specific markers (RSM) was performed using variance analysis. The variance (vari) for RSM can be computed as:(2)(3)(4)where i is a population, n11, n12, and n22 are the number of genotypes in a population, N11, N12 and N22 are the number of genotypes in the whole populations, and p is the allele frequency. The allele frequency in a population and the mean allele frequency for all populations are shown in equation (3). At each marker position, the variances of allele frequency in each population were computed using equations (2) to (4), and then one or more variances that are relatively different from the rest of the populations could be calculated using the z-scores of the variances:(5)where μ and σ are the mean and standard deviation of the variances over all populations at each marker position, respectively. Among the z-scores at a marker position, the marker with at least one population (z-scores >4.0 or <−4.0) was selected. Here, the z-score with 4.0> or <−4.0 indicates that the marker variance of a specific population is significantly greater, or smaller, than the variance mean. Of the selected markers, the first one of two, or more, consecutive markers in 500K, or less distance, was selected and used. In a visual inspection, most of the markers selected with LD (linkage disequilibrium) blocks were located within 200K SNP intervals or less. Therefore, these RSMs are considered as LD-free marker datasets.

Subjects and Genotyping

A total of 320 subjects from 24 regional groups were analyzed in this study. They include Yoruba (YRI), European (CEU), Japanese (JPT), Chinese (CHB), Amerindians (AI), and several population groups from Southern and Northern Asia comprising Chinese from the Jilin area (JL), Vietnamese (VN), Cambodians (CB), Mongolians (MH), and Koreans from ten cities in South Korea (Table 1).

Autosomal SNP genotyping data for the YRI, CEU, JPT, and CHB regional groups were obtained from the International HapMap database (www.hapmap.org) and Amerindian data (AI) were obtained from Affymetrix. Of these, the genotype call rates were sorted, and the best 16 samples from each population of 45 unrelated samples for the JPT and CHB groups (or 60 for the YRI and CEU groups) were selected for this study. For the Southern and Northern Asian groups, genotyping was performed with the Affymetrix GeneChip Mapping 50K_Xba array that comprises 58,960 SNPs, according to the manufacturer's protocol (Affymetrix, Santa Clara, CA, USA). The BRLMM algorithm (Bayesian Robust Linear Model with Mahalanobis distance classifier) [23], [24] was applied to genotype calling for one batch of all samples after initial genotyping with the Dynamic Model algorithm [25]. For the selection of samples in our study, several sample folds with their ages and residential histories were collected to determine whether they were native residents.

Supporting Information

[Figure omitted. See PDF.]

Supporting Methods and Discussion S1.

https://doi.org/10.1371/journal.pone.0011855.s001

(0.07 MB DOC)

Table S1.

Fst, IBS, and Monomorphic marker.

https://doi.org/10.1371/journal.pone.0011855.s002

(0.07 MB XLS)

Figure S1.

MDS without recombination hotspot (RH) filtering. MDS plots with SNPs after QC, QC+500KB spacing between SNPs, and using only regional specific markers (RSM) are shown in columns. Different genotype call rates (98% to 85%) are shown in rows.

https://doi.org/10.1371/journal.pone.0011855.s003

(0.06 MB TIF)

Figure S2.

MDS with recombination hotspot (RH) filtering.

https://doi.org/10.1371/journal.pone.0011855.s004

(0.06 MB TIF)

Figure S3.

PCA without recombination hotspot (RH) filtering. MDS plots with SNPs after QC, QC+500KB spacing between SNPs, and only regional specific markers (RSM) are shown in columns. Different genotype call rates (98% to 85%) are shown in rows.

https://doi.org/10.1371/journal.pone.0011855.s005

(0.06 MB TIF)

Figure S4.

MDS with recombination hotspot (RH) filtering.

https://doi.org/10.1371/journal.pone.0011855.s006

(0.06 MB TIF)

Acknowledgments

We would like to express our gratitude to Affymetrix Inc. for the Amerind chip data and to Sangho Oh for the Regional maps.

Author Contributions

Conceived and designed the experiments: JJ HLK. Performed the experiments: YSC JHO MHR. Analyzed the data: JJ HK. Contributed reagents/materials/analysis tools: HWC JSS JEL BO JB. Wrote the paper: JJ HK.

Citation: Jung J, Kang H, Cho YS, Oh JH, Ryu MH, Chung HW, et al. (2010) Gene Flow between the Korean Peninsula and Its Neighboring Countries. PLoS ONE5(7): e11855. https://doi.org/10.1371/journal.pone.0011855

References

1. Ross MT, Grafham DV, Coffey AJ, Scherer S, McLay K, et al. (2005) The DNA sequence of the human X chromosome. Nature 434: 325–37.MT RossDV GrafhamAJ CoffeyS. SchererK. McLay2005The DNA sequence of the human X chromosome.Nature43432537

2. Kennedy GC, Matsuzaki H, Dong S, Liu WM, Huang J, et al. (2003) Large-scale genotyping of complex DNA. Nat Biotechnol 21: 1233–7.GC KennedyH. MatsuzakiS. DongWM LiuJ. Huang2003Large-scale genotyping of complex DNA.Nat Biotechnol2112337

3. Lell JT, Sukernik RI, Starikovskaya YB, Su B, Jin L, et al. (2002) The dual origin and Siberian affinities of Native American Y chromosomes. Am J Hum Genet 70: 1377–80.JT LellRI SukernikYB StarikovskayaB. SuL. Jin2002The dual origin and Siberian affinities of Native American Y chromosomes.Am J Hum Genet70137780

4. Wang S, Lewis CM, Jakobsson M, Ramachandran S, Ray N, et al. (2007) Genetic Variation and Population Structure in Native Americans. PLoS Genet 3: e185.S. WangCM LewisM. JakobssonS. RamachandranN. Ray2007Genetic Variation and Population Structure in Native Americans.PLoS Genet3e185

5. Excoffier L, Smouse PE, Quattro JM (1992) Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131: 479–91.L. ExcoffierPE SmouseJM Quattro1992Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data.Genetics13147991

6. Watanabe RM, Valle T, Hauser ER, Ghosh S, Eriksson J, et al. (1999) Familiality of quantitative metabolic traits in Finnish families with non-insulin-dependent diabetes mellitus. Hum Hered 49: 159–68.RM WatanabeT. ValleER HauserS. GhoshJ. Eriksson1999Familiality of quantitative metabolic traits in Finnish families with non-insulin-dependent diabetes mellitus.Hum Hered4915968

7. Jin HJ, Kwak KD, Hammer MF, Nakahori Y, Shinka T, et al. (2003) Y-chromosomal DNA haplogroups and their implications for the dual origins of the Koreans. Hum Genet 114: 27–35.HJ JinKD KwakMF HammerY. NakahoriT. Shinka2003Y-chromosomal DNA haplogroups and their implications for the dual origins of the Koreans.Hum Genet1142735

8. Kim JJ, Verdu P, Pakstis AJ, Speed WC, Kidd JR, et al. (2005) Use of autosomal loci for clustering individuals and populations of East Asian origin. Hum Genet 117: 511–9.JJ KimP. VerduAJ PakstisWC SpeedJR Kidd2005Use of autosomal loci for clustering individuals and populations of East Asian origin.Hum Genet1175119

9. Xue Y, Zerjal T, Bao W, Zhu S, Shu Q, et al. (2006) Male demography in East Asia: a north-south contrast in human population expansion times. Genetics 172: 2431–9.Y. XueT. ZerjalW. BaoS. ZhuQ. Shu2006Male demography in East Asia: a north-south contrast in human population expansion times.Genetics17224319

10. The HUGO Pan-Asian SNP Consortium (2009) Mapping Human Genetic Diversity in Asia. Science 326: 1541–5.The HUGO Pan-Asian SNP Consortium2009Mapping Human Genetic Diversity in Asia.Science32615415

11. Patterson N, Price AL, Reich D (2006) Population structure and eigenanalysis. PLoS Genet 2: e190.N. PattersonAL PriceD. Reich2006Population structure and eigenanalysis.PLoS Genet2e190

12. Mountain J, Cavalli-Sforza L (1997) Multilocus genotypes, a tree of individuals, and human evolutionary history. Am J Hum Genet 61: 705–718.J. MountainL. Cavalli-Sforza1997Multilocus genotypes, a tree of individuals, and human evolutionary history.Am J Hum Genet61705718

13. Tamura K, Nei M, Kumar S (2004) Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci USA 101: 11030–11035.K. TamuraM. NeiS. Kumar2004Prospects for inferring very large phylogenies by using the neighbor-joining method.Proc Natl Acad Sci USA1011103011035

14. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959.JK PritchardM. StephensP. Donnelly2000Inference of population structure using multilocus genotype data.Genetics155945959

15. Salmela E, Lappalainen T, Fransson I, Andersen PM, Dahlman-Wright K, et al. (2008) Genome-wide analysis of single nucleotide polymorphisms uncovers population structure in Northern Europe. PLoS One 3: e3519.E. SalmelaT. LappalainenI. FranssonPM AndersenK. Dahlman-Wright2008Genome-wide analysis of single nucleotide polymorphisms uncovers population structure in Northern Europe.PLoS One3e3519

16. Rosenberg NA (2005) DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes 4: 137–138.NA Rosenberg2005DISTRUCT: a program for the graphical display of population structure.Mol Ecol Notes4137138

17. Press WH, Flannery BP, Teukolsky SA, Vetterling WT (2007) Numerical Recipes in C: the art of scientific computing. Cambridge University Press. WH PressBP FlannerySA TeukolskyWT Vetterling2007Numerical Recipes in C: the art of scientific computingCambridge University Press

18. Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 24: 1596–9.K. TamuraJ. DudleyM. NeiS. Kumar2007MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0.Mol Biol Evol2415969

19. Bowcock AM, Ruiz-Linares A, Tomfohrde J, Minch E, Kidd JR, et al. (1994) High resolution of human evolutionary trees with polymorphic microsatellites. Nature 368: 455–457.AM BowcockA. Ruiz-LinaresJ. TomfohrdeE. MinchJR Kidd1994High resolution of human evolutionary trees with polymorphic microsatellites.Nature368455457

20. Bingham J, Sudarsanam S (2000) Visualizing large hierarchical clusters in hyperbolic space. Bioinformatics 16: 660–661.J. BinghamS. Sudarsanam2000Visualizing large hierarchical clusters in hyperbolic space.Bioinformatics16660661

21. The International HapMap Consortium (2005) A haplotype map of the human genome. Nature 437: 1299–1320.The International HapMap Consortium2005A haplotype map of the human genome.Nature43712991320

22. Terwilliger JD, Ott J (1993) A novel polylocus method for linkage analysis using the lod-score or affected sib-pair method. Genet Epidemiol 10: 477–82.JD TerwilligerJ. Ott1993A novel polylocus method for linkage analysis using the lod-score or affected sib-pair method.Genet Epidemiol1047782

23. Affymetrix in Technical Report, 2006. Affymetrix in Technical Report, 2006.

24. Rabbee N, Speed TA (2006) A genotype calling algorithm for Affymetrix SNP arrays. Bioinformatics 22: 7–12.N. RabbeeTA Speed2006A genotype calling algorithm for Affymetrix SNP arrays.Bioinformatics22712

25. Di X, Matsuzaki H, Webster TA, Hubbell E, Liu G, et al. (2005) Dynamic model based algorithms for screening and genotyping over 100 K SNPs on oligonucleotide microarrays. Bioinformatics 21: 1958–1963.X. DiH. MatsuzakiTA WebsterE. HubbellG. Liu2005Dynamic model based algorithms for screening and genotyping over 100 K SNPs on oligonucleotide microarrays.Bioinformatics2119581963