The Origin of Amerindians and the Peopling of the Americas According to HLA Genes: Admixture with Asian and Pacific People.

The classical three-waves theory of American peopling through Beringia was based on a mixed anthropological and linguistic methodology. The use of mtDNA, Y chromosome and other DNA markers offers different results according to the different markers and methodologies chosen by different authors. At present, the peopling of Americas remains uncertain, regarding: time of population, number of peopling waves and place of peopling entrance among other related issues. In the present review, we have gathered most available HLA data already obtained about First Native American populations, which raise some doubts about the classical three waves of American peopling hypothesis. In summary, our conclusions are: 1) North West Canadian Athabaskans have had gene flow with: a) close neighboring populations, b) Amerindians, c) Pacific Islanders including East Australians and d) Siberians; 2) Beringia was probably not the only entrance of people to America: Pacific Ocean boat trips may have contributed to the HLA genetic American profile (or the opposite could also be true); 3) Amerindians entrance to America may have been different to that of Athabaskans and Eskimos and Amerindians may have been in their lands long before Athabaskans and Eskimos because they present and altogether different set of HLA-DRB1 allele frequencies; 4) Amerindians show very few "particular alleles", almost all are shared with other Amerindians, Athabaskans and Pacific Islanders, including East Australians and Siberians; 5) Our results do not support the three waves model of American peopling, but another model where the people entrance is not only Beringia, but also Pacific Coast. Reverse migration (America to Asia) is not discarded and different movements of people in either direction in different times are supported by the Athabaskan population admixture with Asian-Pacific population and with Amerindians, 6) HLA variability is more common than allele veriability in Amerindians. Finally, it is shown that gene genealogy analises should be completed with allele frequency analyses in population relatednes and migrations studies.

INTRODUCTION

The First Amerindian Natives are postulated to have come from Asia through the Bering land bridge between 30,000–12,000 years before the present (BP). These conclusions have been based on cultural, morphological and genetic similarities between American and Asian populations. Both Siberia [1] and Mongolia [2,3] have been put forward as the most likely places of origin in Asia.

Greenberg first postulated the triple migration theory Fig. () for explaining the peopling of the Americas [4]: Amerindians (most North and South American Indians; 12,000 years BP), Na-Dene (Athabascans, Navajo, Apache; 8,000 years BP) and Eskimo-Aleuts (6,000 years BP). Research carried out before the widespread use of Y Chromosome (Y Chr) and other nuclear DNA markers including mtDNA [5] for the study of populations [6,7] supported the three-wave model. However, other mtDNA studies have not [8,9]; other authors postulate only one wave coming from Mongolia / North China as giving rise to the First Native American ancestors [2,3]. The study of Y Chromosome DNA markers seemed to suggest the existence of a single major paternal haplotype in both North and South American Native populations [10,11]. However, other studies on Y Chromosome show that more than one paternal founder haplotypes arrived in America during different migrations [12], probably from Siberia [13]. See also Fig. () [14-17].

More recently, new mtDNA analysis has suggested that all mtDNA lineages must have been isolated in Asia before entering the New World by at least 7-15 thousand years. They even suggest that this place must have been Beringia [18]. Also, a dispersal of Amerindians coming from Asia has been put forward through Coastal Pacific line [19] based on all available archaeological, anthropological and mtDNA and genetic data.

All these calculations are done by using paternal (Y Chr) or maternal (mtDNA) lineages may be biased when populations displacements are concerned, as in the putative Amerindians displacement from Asia to the Americas. In addition, other authors [20] using nuclear histocompatibility (HLA) markers do not regard as important and possible to establish the number and timing of migration waves. The important issue is whether immigrants (Amerindians) were already differentiated (in Asia) into such ethnic groups whose descendants are still to be found in Asia. If they were differentiated then the question of how and when they crossed the Bering Land Bridge is a secondary one [20].

Y Chr and mtDNA studies seem unable to resolve this question, mainly because the studies have been gene- rather than population- based. Their emphasis has been on variant genealogies rather than on population frequencies studies. In this regard HLA data may be more informative [20] because maternal and paternal lineages and both frequencies (i.e.: genetic distances, dendrograms and correspondence analyses) and genealogies (quasi-specific HLA alleles and haplotypes) may be studied for comparing populations.

Alu-insertion investigations have also been carried out to ascertain the origin of First Americans [21]. The results are not concordant with the multiple-wave migration hypothesis; a surprisingly short genetic distance between Chinese and Native Americans was found and explained by a recent gene flow from Asia [21].

A Trans-Pacific route of American peopling from Asia or Polynesia has been suggested because HTLV-1 virus strains shared identical sequences in Japan and in the northern coast of South America [22] and some HLA alleles may have been introduced by the same Trans-Pacific route [15]. Finally, both genetic [17] and archaeological [16] evidence suggests that a two-way Trans-Atlantic traffic occurred before Columbus discovered America Fig. (); archaeologists in New Mexico have recently found tools used 20,000 years ago in Spanish Solutrean culture [16].

All these discrepancies and uncertainties about Amerindian origins may be due to methodological (sampling) differences and also to the different genealogical history of each genetic marker and/or to the phylogenetic usefulness of different DNA markers. For instance, functional molecules — cytochrome (cyt) b mtDNA— are used against an admixture of intronic and exonic DNA markers, as in the Alu or STR studies: the obtained molecular genetics history could not be the same one. In addition, population movements should be studied like a “group of genes” movements, i.e.: with gene frequencies (genetic distances, dendrograms and correspondence analyses), which better reflect a population displacement and other populations (Asian / Amerindian) relatedness, and afterwards completed with genealogies (quasi-specific HLA alleles, HLA haplotypes, mtDNA and Y Chr markers).

Thus, in the present work, we have studied the North, Meso and South American Amerindians’ HLA gene frequencies and compared it with those of other North American Indians and worldwide populations, particularly with Asian and Pacific populations. Also, we have studied the following Amerindian ethnic groups: Seri, Mixe, Mixtecans, Zapotecans, Guaranis [23], Lakota Sioux [24], Mazatecans [25], Teeneks [26], Mayans [27], Kogi, Arsario, Arhuacs, Wayu [28], Cayapa [29], Lamas [30], Aymaras [31], Quechuans [32], Terena Indians [33], Xavantes, Mayos [34], Uros [35], Nahuas [36], Tarahumaras [37], Toba Pilaga, Mataco Wichi, Eastern Toba [15], Mexican Mestizos and Jaidukama (unpublished results) and also Aleuts [14].

Our aims are: 1) To determine the HLA class I (A and B) and class II (DRB1 and DQB1) quasi-specific Amerindian allelic lineages (hereafter ‘‘alleles’’ for simplicity) or specific HLA haplotypes by using DNA sequencing and serology; in other words, the most frequent HLA alleles and haplotypes in Amerindians which do not exist or exist in very low frequency in other populations, i.e.: genealogy comparisons and 2) To compare the Amerindians HLA allele frequencies with those of other First American Natives (Na-Dene, Eskimo and Aleuts) and also those of other worldwide populations in order to clarify the still unclear peopling of the Americas and the origins of Amerindians, i.e.: groups of genes comparisons by using genetic distances, Neighbor Joining (NJ) dendrograms and correspondence analyses.

RESULTS AND DISCUSSION

DRB1 Alleles and HLA-Extended Haplotypes: North-American and Meso / South-American Populations

The low number of class I alleles found may be artificial, since many of them may not have been yet detected. In fact, many more HLA Class I alleles are defined at present day (www.anthonynolan.org.uk).

We have chosen DRB1 alleles because many populations are typed for DRB1 high resolution alleles and very few for HLA class I or other class II loci. No completely specific DRB1 alleles are found in North or South American populations: some of the alleles are found in other populations in a very low frequency Fig. (), footnote, [38]. At the moment, the only exceptions are DRB1*0411 and DRB1*0417 alleles, which are only found in all studied Meso and South Amerindians (Table ), Fig. (), bold red color (www.antonynolan.org.uk), compared to previous times.

Notwithstanding, some Meso and South American DRB1 Amerindian alleles tend to be quasi-specific Table , Fig. (), red color, but not North American alleles which are clearly shared with other non-Amerindian populations Fig. (), blue color. This is concordant with the existence of gene flow between Amerindians and Pacific or Siberian people, but not necessarily with a migration of Amerindians from Asia or Pacific Areas, although there are signs of cultural or genetic contacts with Asia [19] or even with Iberians [16,17], Fig. (). DRB1*0802 and DRB1*0407 are present in all the most frequent Meso-American haplotypes Fig. (). DRB1*0802 is also present in Siberian Eskimos and Japanese Ainu: otherwise it is present in almost all Amerindians. DRB1*0407 is present in almost all Amerindian populations and absent or in a non-significant frequency in other populations. DRB1*0403 is present in one South American most frequent haplotype Fig. () but also is found in high frequency throughout Pacific Islands (Samoa, Papua New Guinea, New Zealand Maories, Taiwan, Tonga, Cook Islands) [38]. A Pacific gene flow in either direction may not be discarded by this genealogy approach also; HLA frequency data (elaborated in dendrograms and correspondence analyses) separate more Amerindians from other populations Figs. (, ), see below). DRB1*0407 is an Amerindian allele in one of the most frequent South American Amerindian (Table ) [39-57 and previous references], Fig. (), bold red colour [58].

Both genealogy (extended haplotypes, HLA-A,-B,-DRB1,-DQB1) and allele frequency in population analyses (Correspondence and NJ multidimensional populations relatedness) have been carried out.

North-Americans

The relatively low number of class I alleles found some years ago in Amerindians [57] compared to other worldwide populations may be due to the fact that techniques were not by then detecting new class I alleles and many of the alleles had not yet been described.

Fig. () shows that the most frequent extended haplotypes in North Americans are specific for North American populations, Yupik (Eskimos), Aleuts and one of the five most frequent haplotypes are also found in Taiwan and Japanese populations (A*24-B*40-DRB1*1401-DQB1*0503). This shows that a low degree of North American haplotypes sharing is found between North American and Asian-Pacific populations. However, there is a clear genetic HLA relatedness between isolated populations close to Beringia: Eskimos, Udegeys, Nivkhs (North East coast of Siberia) and Koryaks and Chukchi from extreme North East Siberia Figs. (, , ), and North West American populations: Athabaskan, Alaskan Eskimos (Yupik) and Tlingit.

These results in which class II high resolution alleles and also specific class I-class II extended haplotypes are used suggest that admixture occurred between extreme North East Siberian groups and North American Na-Dene (including Tlingit) and Eskimo (Yupik) people. However, results do not indicate anything about direction of admixture or whether migrations in both directions occurred.

On the other hand, Asian populations which are geographically not close to Beringia (Japanese, Ainu, Manchu, Singapore Chinese, Buyi) do not cluster with North Americans neither in NJ dendrogram Fig. () or correspondence analysis Fig. ().

Finally, Lakota-Sioux Amerindians which have inhabited in North United States, are not related with Asians and West Siberians Figs. (, , ) but with Meso and South Americans.

Meso-Americans

Most frequent haplotypes Fig. (), relatedness dendrograms Fig. () and correspondence Fig. () do not relate these Amerindians with any Asian population, including North East Siberians. Haplotypes of Meso-Americans are shared with other Amerindians and one of them with Alaskan Eskimo (Yupik): A*02-B*35-*DRB1*0802-DQB1 *0402.

South-Americans

These Amerindian speaking groups are related to other South-American Amerindians and to Meso-American Amerindians Figs. (, ). Most frequent haplotypes are shared with other American Amerindians, but not with Asians Fig. ().

In summary, the general view after analyzing the most frequent extended haplotypes (genealogy) is that Amerindians have little relatedness with Asians; this is also confirmed by allele frequencies in populations and the derived analyses, Figs. (, ). Genealogy studies are less suitable for comparing and relating groups of people [20]. Also, comparing HLA four loci most frequent extended haplotypes of North-Americans only share one haplotype (A*24-B*40-DRB1*1401-DQB1*0503) with Taiwanese and Japanese in low frequencies see Fig. () footnote.

Specific Extended Haplotypes for Amerindian Ethnic Groups and Aleuts

Some new extended four loci haplotypes have been found only in Amerindian and Aleut specific groups and in no other either Amerindians or World populations (Table ). It is striking that in small groups of people apparently specific HLA four loci recombinations occur and are fixed. The evolutive forces to achieve an appropriate extended haplotype may be advantageous for a population to deal with its specific environmental pathogens [59]. In this case, evolutive forces to drive de novo HLA haplotype appearance must include pathogens and not low frequency genes driven selection. However, both of these types of selection for inducing HLA (or other genes) variability and allele fixation in populations are mathematically indistinguishable [60].

In addition, this high haplotype variability in relatively small ethnic groups is concordant with the fact that HLA haplotype frequencies show greater variation among racial groups than individual alleles [61]. Thus, new allele appearance may be a relatively rare event (usually by a gene conversion mechanism [62]), compared with a new haplotype appearance. In fact, evolution for variability in the MHC may be more frequent in haplotype and not in allele diversification.

See references: [25-27,30-32,34-36,14].

New HLA alleles are continuously being described in populations (Anthony Nolan database: http://www.anthonynolan.org.uk, [38]) but this does not necessarily mean that are continuously being produced; they may have been fixed as low frequency alleles in populations for a long time and only described at present times according to technology advances.

In conclusion, new specific haplotypes are found in North and South American Amerindians, while specific alleles for a particular Amerindian population are rarely found or not found; new alleles are being newly described in a particular population but later found in others. Selection for variability within North and South American Natives is acting upon haplotypes more than upon alleles Fig. (). This finding may be universal for all World populations.

Genes and Languages

It was postulated that genes correlated with languages [63]; however, from Fig. () (NJ dendrogram) it may be seen that Na-Dene / Eskimo / Siberian / group is genetically very close as measured by HLA-DRB1 frequencies and speak distant languages [64]; this is confirmed by HLA-DRB1 and HLA-DQB1 correspondence analysis Fig. (). Both NJ and correspondence analysis correlates quite well with geography but does not correlate with languages. Some authors find correlation between genes and languages when selected ethic groups and selected languages are used but only at macro-geographical level.

mtDNA and Y Chr Markers

Specialists have studied genealogies of haplotypes and / or other markers [5,11,19,18]. However, for studying populations relatedness and migration genetics it is more suitable to study gene frequencies (or allele frequencies) as done in Figs. (, ). It is remarkable that Northern USA Lakota-Sioux is genetically close to other Amerindians from South America as expected. However, Canadian Athabaskans go together with Northern Americans First Inhabitants and Siberians Figs. (, ).

In general, mtDNA and Y Chr markers have studied the origins of Amerindians [11,19,18], postulated time and place of entrance of Amerindians, Athabaskans and Eskimo [19], and whether the genetic findings fit with Greenberg linguistically based theory [4] (three waves of Americas peopling). Also, archeological findings have been contrasted with genetic findings [1,16,19]. In the end, conclusions are diverse and no consensus exists about Amerindian origin and relatedness [2,3,5,8-13,19].

HLA Amerindian Genetic Anthropology

HLA has been ignored by anthropologists to study genetics probably because a lack of HLA understanding and the claim that HLA is shaped by selection because of disease linkage; mtDNA and Y Chr (OMIM, http://ncbi.nlm.nih.gov/omim) are also linked to diseases. In addition, they only give a paternal or maternal view of markers genealogy which may have been divergent in small primitive colonizing populations. However, HLA is a nuclear marker giving an even genealogy and genetic history for both sexes. The best test showing that HLA is a good genetic marker for studying population relatedness is that it usually correlates with geography.

HLA dendrograms or correspondence analyses based on HLA frequencies show that Amerindians (in the sense of Greenberg definition, [4]) seem to be separated from other World wide populations, including northern Canadian Athabaskans and Eskimos. The latter cluster in Fig. () with Siberians. This means that HLA-DRB1 and HLA-DQB1 allele frequencies are completely different in Amerindian compared to from other First American Natives or other World populations.

In addition, the studied populations show particular HLA four loci haplotypes for each specific studied population (Table ). This is not the case for HLA-DRB1 alleles: except for two DRB1 alleles —DRB1*0411 and DRB1*0417, Fig. () Meso and South American Amerindians alleles are shared with the following populations: 1) Siberians, 2) other First American Inhabitants including Athabaskans and Eskimo, but not Aleuts [14], probably because Aleuts come from a Baikal Lake ancient stock that are related to both Aleuts and European Lapps (Saami), 3) Asian Pacific Coast populations (Ainu, Japanese, Taiwan) and to a lesser extent with Indochina people and, 4) East Australian Aborigins and Pacific Islands, like Papua New Guinea or Samoa groups.

In this context and because Canadian Athabaskans have been placed in the postulated entrance for American peopling [4], an haplotype study was done by computing the most frequent Athabaskan two loci haplotypes (DRB1-DQB1) [55] and looking where this part of Chromosome 6 was found around the World (i.e.: a genealogy study, which completes our population frequency study). Results are shown in Table (see footnote). Athabaskans DRB1-DQB1 genes are shared with: 1) Neighbors, including Alaskan Eskimo (Yupik), 2) Amerindians from North and South America, 3) Siberians, 4) Pacific Islands inhabitants, including those of Samoa, Papua New Guinea, Cook Islands and Japanese-Ainu and even Eastern Australia Aborigins. Only one haplotype (f), in Table is exclusively shared among Amerindians. This suggests that Athabaskans are composed of a genetic HLA admixture and that gene flow has occured between Athabaskans and all the other above mentioned Pacific-Asian populations. The direction of the hypothetical HLA gene flow is not known and may have occured in different directions in different times. Thus, there is no point to conclude about one or more waves of Americas peopling from our data. However, this also shows that not only Beringia was an active pass of primitive Amerindians, but also Pacific navigation was.

Results on North American and South American population HLA alleles also support this view Fig. (). Why nowadays North and South American Amerindians are altogether different populations to the rest of the World regarding to HLA frequencies is only a matter of speculation:

It has been calculated that about 80 million First American Natives died after 1492 AD within the following 100 years from Alaska to Patagonia [65]. It was mainly due to a lack of appropriate immune response to European-borne diseases, mainly measles, influenza and plague [65]. This may have shaped the First American Natives HLA profile by increasing rare HLA alleles able to present new pathogens to T cells [60,61,66]. However, this is not likely since First North American Natives (non-Amerindians) also suffered many epidemics [65] and do not have as different HLA profile from Asians, as Amerindians do.

First American Natives, including Alaskan Eskimo (Yupik) and Athabaskans must have been in America long before calculated (20,000 years ago), because both North and South American Fist Natives were similarly susceptible to European-borne diseases [65], indicating the existence of long population isolation.

North West Canadian Athabaskans have had gene flow with a) close neighboring populations, b) Amerindians, c) Pacific Islanders including East Australians and d) Siberians.

Beringia was not probably the only entrance of people to Americas: Pacific Ocean boat trips may have contributed to the HLA genetic American profile (or the opposite could also be true).

Amerindians entrance to America may have been different to that of Athabaskans and Eskimos and Amerindians may have been in their lands long before Athabaskans and Eskimos because they present and altogether different sets of HLA-DRB1 frequencies.

Amerindians show very few “particular alleles”, almost all are shared with other Amerindians, Athabaskans, Pacific Islanders, including East Australians and Siberians. However, specific Amerindian haplotypes are found in isolates.

Genes and languages do not correlate.

Our results do not support the three-wave model of American peopling, but another model where the people entrance is not only Beringia, but also Pacific Coast. Reverse migration (America to Asia) is not discarded and different movements of people in either direction in different times are supported by the Athabaskan population admixture with Asian-Pacific population and with Amerindians.

Table 1

Populations Studied in the Present Work. A Total of 14,698 Chromosomes were Analyzed

Population	N	Reference	Population	N	Reference
Seri	100	[23]	Manchu	50	[39]
Mixe	55	[23]	Koreans	100	[39]
Mixtecans	103	[23]	Japanese	493	[39]
Zapotecans	75	[23]	Ainu	50	[40]
Lakota Sioux	302	[24]	Khalk Mongolians	100	[41]
Mazatecans	90	[25]	Tuvinians	190	[42]
Teeneks	44	[26]	Khoton Mongolians	85	[41]
Mexican Mestizos	99	Unpublished results	Sardinians	91	[39]
Mayans	132	[27]	Italians	284	[39]
Wayu	88	[28]	French	179	[39]
Arhuaco	107	[28]	Spaniards	176	[43]
Kogi	42	[28]	Spanish Basques	82	[43]
Arsario	18	[28]	Algerians	106	[44]
Jaidukama	39	Unpublished results	Berbers (Souss)	98	[45]
Cayapa	100	[29]	Moroccans	96	[46]
Lamas	83	[30]	Macedonians	172	[47]
Aymaras	102	[31]	Cretans	135	[48]
Quechuas	80	[32]	Finns	157	[14]
Xavantes	74	[15]	Saami (Lapps)	81	[14]
Terena Indians	60	[33]	Aleuts	85	[14]
Guarani	32	[23]	Pomors	73	[14]
Toba-Pilaga	19	[15]	Danish	124	[39]
Mataco-Wichi	49	[15]	West Russians	200	[49]
Eastern Toba	135	[15]	Chuvashians	82	[50]
Mayos	60	[34]	Fidji Islands	57	[51]
Tarahumaras	44	[37]	Papua New Guinean	65	[51]
Uros	105	[35]	Central Desert	152	[52]
Nahuas	85	[36]	Yuendumu	119	[52]
Eskimos	80	[53]	Kimberley	82	[54]
Athabaskans	124	[55]	Western Samoa	51	[56]
Tlingit	53	[39]	Madang	65	[51]
Nivkhs	32	[53]	Rabaul	60	[51]
Udegeys	25	[53]	New Caledonia	65	[51]
Koryaks	92	[53]	Cape York	80	[54]
Chukchi	59	[53]	South American Blacks	59	[39]
Kets	22	[53]	North American Blacks	447	[39]
Evenks	35	[53]	Hottentots	65	[39]
Singapore Chinese	71	[39]	Bushmen	103	[39]
Buyi	70	[39]

Table 2

Extended Haplotypes Only Found in Some Amerindian Populations [14,25-27,30-32,34-36]

HAPLOTYPE	(Freq. %)	POPULATION
A*02-B*39-DRB1*1602-DQB1*0301	(3.3%)	Mazatecans
A*02-B*62-DRB1*1602-DQB1*0301	(3.3%)	Mazatecans
A*02-B*15-DRB1*0404-DQB1*0302	(1.5%)	Mayans
A*02-B*39-DRB1*0802-DQB1*0402	(3.4%)	Aymaras
A*02-B*39-DRB1*0901-DQB1*0303	(3.4%)	Aymaras
A*02-B*48-DRB1*0403-DQB1*0302	(7.8%)	Lamas
A*02-B*48-DRB1*0804-DQB1*0402	(7.8%)	Lamas
A*02-B*39-DRB1*1402-DQB1*0301	(3.6%)	Lamas
A*66-B*41-DRB1*1303-DQB1*0301	(3.6%)	Lamas
A*02-B*48-DRB1*0411-DQB1*0302	(2.4%)	Lamas
A*24-B*15-DRB1*0901-DQB1*0303	(1.8%)	Lamas
A*33-B*38-DRB1*1104-DQB1*0301	(1.8%)	Lamas
A*68-B*35-DRB1*0802-DQB1*0402	(3.6%)	Quechua
A*02-B*48-DRB1*1402-DQB1*0301	(2.9%)	Quechua
A*02-B*48-DRB1*0802-DQB1*0402	(2.2%)	Quechua
A*02 B*52 DRB1*0411 DQB1*0302	(3.7%)	Teeneks
A*68 B*35 DRB1*1402 DQB1*0301	(2.8%)	Teeneks
A*68 B*40 DRB1*1602 DQB1*0301	(2.6%)	Teeneks
A*68 B*35 DRB1*1406 DQB1*0301	(2.6%)	Teeneks
A*02-B*35-DRB1*1406-DQB1*0301	(4.2%)	Mayos
A*02-B*48-DRB1*0404-DQB1*0302	(3.3%)	Mayos
A*24-B*51-DRB1*0407-DQB1*0302	(3.3%)	Mayos
A*02-B*08-DRB1*0407-DQB1*0302	(2.5%)	Mayos
A*30-B*49-DRB1*1001-DQB1*0501	(7.5%)	Nahuas
A*02-B*52-DRB1*1402-DQB1*0301	(2.7%)	Nahuas
A*68-B*61-DRB1*1602-DQB1*0303	(2.0%)	Nahuas
A*24-B*15-DRB1*1402-DQB1*0301	(3.2%)	Uros
A*68-B*35-DRB1*0403-DQB1*0302	(3.2%)	Uros
A*24-B*48-DRB1*0403-DQB1*0302	(2.2%)	Uros
A*02-B*40-DRB1*0101-DQB1*0501	(5.6%)	Aleuts
A*24-B*37-DRB1*0801-DQB1*0402	(4.2%)	Aleuts
A*24-B*39-DRB1*0404-DQB1*0302	(4.2%)	Aleuts
A*24-B*39-DRB1*1201-DQB1*0301	(4.2%)	Aleuts
A*02-B*15-DRB1*0401-DQB1*0301	(2.8%)	Aleuts
A*02-B*51-DRB1*1501-DQB1*0602	(2.8%)	Aleuts
A*26-B*40-DRB1*1401-DQB1*0503	(2.8%)	Aleuts
A*32-B*44-DRB1*0701-DQB1*02	(2.8%)	Aleuts
A*68-B*40-DRB1*0404-DQB1*0302	(2.8%)	Aleuts
A*68-B*40-DRB1*0802-DQB1*0402	(2.8%)	Aleuts
A*68-B*39-DRB1*1201-DQB1*0301	(2.8%)	Aleuts

Table 3

DRB1-DQB1 Haplotypes also Found in Athabaskans [55] and also in other Populations

DRB1	DQB1	FREQ. (%)
1402	0301	34.7a
1401	0503	16.9b
0901	0303	10.5c
0403	0302	9.6d
1201	0301	8.9e
0802	0402	4.8f
0410	0402	4.0g

aXavantes (25.7%); Mataco Wichi (22.4%); Alaska Yupik Natives (22.0%); Tarahumara (11.9%); Russia Siberia Eskimos Chukotka Paninsula (11.3%).

Ainu (20.0%); Spain Malaga Gipsy (9.5%); Slovenia (7.9%); Papua New Guinea New Britain Tolai (7.6%); Aleuts (6.9%); Alaska Yupik Natives (6.7%).

Japan (29.5%); Russia Siberian Khabarovsk Evenki (26.0%); Samoa (25.9%); Russia Siberia Nganasan Dudinka (25.0%); Russia Siberia Koryaks North East Kamchatka (22.3%); Russia Siberia Negidal (20.0%); Russia Siberia Udegeys (19.0%); Russia Siberia Ulchi (15.8%); Russia Siberia Chukchi (14.7%).

Samoa (17.2%); Russia Siberia Nganasan (12.5%); Taiwan Aboriginal (9.5%); Yucpa (9.3%); Russian Siberia Kushun (8.0%); Lamas (7.8%); Russia Siberia Khabarovsk Evenki (6.0%); Russia Siberia Negidal (5.7%); Philippines (5.2%); Zapotecans (5.2%).

Russia Siberia Ket Lower Yenisey (17.7%); Russia Siberia Irkutsk Tofalar (17.4%); Russia Tuva (13.7%); Russia Siberia Kets Sulamai (13.6%); Russia Siberia Ulchi (13.0%); Russia Siberia Nganasan Dudinka (12.5%); Russia Siberia Kets Sulamai (11.4%); Russia Siberia Khanty Mansi (8.8%).

Mixe (28.0%); Tarahumara (27.9%); Xavantes (23.0%); Mixtecans (21.6%); Zapotecans (21.5%); Eastern Toba (18.9%); Jalisco Mestizos (14.9%); Highlands Mestizos (13.8%); Alaska Yupik Natives (13.3%); Toba Pilaga (10.5%); Aymaras (10.4%); Ainu (10.0%).

Japan (3.8%); Alaska Yupik Natives (3.6%); Alaska Yupik Natives (3.5%); Australia New South Wales Aborigines (2.4%); Papua New Guinea (2.3%); Ainu (2.0%), China (1.8%); Japan (1.7%); Japan Central (1.6%); Mixtecans (1.5%).

These footnote frequencies were taken from reference [38] and from our own publications (Table 1).