Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Highly Accessed Research article

Large scale mitochondrial sequencing in Mexican Americans suggests a reappraisal of Native American origins

Satish Kumar*, Claire Bellis, Mark Zlojutro, Phillip E Melton, John Blangero and Joanne E Curran*

Author affiliations

Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas 78227, USA

For all author emails, please log on.

Citation and License

BMC Evolutionary Biology 2011, 11:293  doi:10.1186/1471-2148-11-293

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/11/293


Received:3 June 2011
Accepted:7 October 2011
Published:7 October 2011

© 2011 Kumar et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The Asian origin of Native Americans is largely accepted. However uncertainties persist regarding the source population(s) within Asia, the divergence and arrival time(s) of the founder groups, the number of expansion events, and migration routes into the New World. mtDNA data, presented over the past two decades, have been used to suggest a single-migration model for which the Beringian land mass plays an important role.

Results

In our analysis of 568 mitochondrial genomes, the coalescent age estimates of shared roots between Native American and Siberian-Asian lineages, calculated using two different mutation rates, are A4 (27.5 ± 6.8 kya/22.7 ± 7.4 kya), C1 (21.4 ± 2.7 kya/16.4 ± 1.5 kya), C4 (21.0 ± 4.6 kya/20.0 ± 6.4 kya), and D4e1 (24.1 ± 9.0 kya/17.9 ± 10.0 kya). The coalescent age estimates of pan-American haplogroups calculated using the same two mutation rates (A2:19.5 ± 1.3 kya/16.1 ± 1.5 kya, B2:20.8 ± 2.0 kya/18.1 ± 2.4 kya, C1:21.4 ± 2.7 kya/16.4 ± 1.5 kya and D1:17.2 ± 2.0 kya/14.9 ± 2.2 kya) and estimates of population expansions within America (~21-16 kya), support the pre-Clovis occupation of the New World. The phylogeography of sublineages within American haplogroups A2, B2, D1 and the C1b, C1c andC1d subhaplogroups of C1 are complex and largely specific to geographical North, Central and South America. However some sub-branches (B2b, C1b, C1c, C1d and D1f) already existed in American founder haplogroups before expansion into the America.

Conclusions

Our results suggest that Native American founders diverged from their Siberian-Asian progenitors sometime during the last glacial maximum (LGM) and expanded into America soon after the LGM peak (~20-16 kya). The phylogeography of haplogroup C1 suggest that this American founder haplogroup differentiated in Siberia-Asia. The situation is less clear for haplogroup B2, however haplogroups A2 and D1 may have differentiated soon after the Native American founders divergence. A moderate population bottle neck in American founder populations just before the expansion most plausibly resulted in few founder types in America. The similar estimates of the diversity indices and Bayesian skyline analysis in North America, Central America and South America suggest almost simultaneous (~ 2.0 ky from South to North America) colonization of these geographical regions with rapid population expansion differentiating into more or less regional branches across the pan-American haplogroups.

Background

Archeological, anatomical, linguistic, and genetic evidence have collectively shown that the original human inhabitants of the Western Hemisphere most likely arrived from Asia across an exposed land mass between northeastern Siberia and Alaska, referred to as Beringia, during the last glacial maximum (LGM) approximately 23 to 19 thousand years ago (kya) [1-6]. However, despite the general consensus of an Asian origin for Native Americans, a number of uncertainties persist regarding the dynamic nature of the peopling of the Americas, including the ancestral source population(s) within Asia, the arrival time(s) of the founder groups, the number of expansion events, and the specific migration routes into the New World [2,7-15].

Mitochondrial DNA (mtDNA) data presented over the past two decades have shown that Native American populations exhibit, almost exclusively, five mtDNA haplogroups (A-D and X) [6] classified in the autochthonous American lineages A2, B2, C1, D1, and X2a [16]. Haplogroups A - D are found throughout the New World and are frequent in Asia, supporting a northeastern Asian origin of these lineages [17,18]. The pan-American distribution of haplogroups A2, B2, C1 and D1, along with their similar levels of diversity and estimates for coalescence time, has been used to support a single-migration model [19-22]. Nested within this model, the Beringian land mass has been suggested to play an important role as a region where ancestral populations crossed and perhaps settled prior to their eventual expansion into the Americas, leading to the diversification of the New World founding lineages [23-26]. The evolutionary history of the fifth founder lineage, X2a, remains elusive and highly debated, with some scholars arguing for an additional migration event into the New World to account for the X lineage, one that was independent from the founder populations responsible for introducing the pan-American haplogroups [27,28].

The "Beringian incubation model" (BIM) [23] and its variants [21,25,26] emphasize that the Native American founder population reached greater Beringia by 30 kya, marked by the earliest evidence of human habitation in northeastern Siberia Yana Rhinoceros Horn Site [29]. According to these models, the founder populations remained in Beringian LGM refugium for about 5-15 kya, ecologically isolated to the west and physically isolated to the east by the glaciers that are believed to have effectively blocked the way to America until near the end of the LGM [24]. During this time, mtDNA lineages of American founders are likely to have differentiated from their Asian precursors. As the Laurentide and Cordilleran ice sheets retreated, the paused populations rapidly colonized the double continent around 14-16 kya, resulting in largely autochthonous patterns of variation within the continental founder haplogroups [23,25,26].

However, the timing of this peopling scenario has been questioned due to the existence of pre-Clovis archaeological sites, such as the Monte Verde site in southern South America that dates human habitation by at least 14.5 kya [30], as well as demographic models based on mtDNA variation that indicate pronounced growth of Native American founder populations starting at ~19 to18 kya [24], thus marking the effective colonization of the New World prior to the opening of the ice free corridors in North America. To accommodate these dates, an alternate early coastal route has been suggested for the initial colonization of the Americas, which would have been largely ice free at ~19 kya [24,31-33].

In northeastern Siberia, the putative origin of the Native American founder populations, there is archaeological evidence for human presence as early as 30 kya [29]. The estimated dates of the Beringian 'pause', (ranging from 5 to 15 thousand years) [23,24,26] and its potential role as a time when American founders diverged from their Siberian-Asian/Asian ancestors represent interesting avenues for microevolutionary reconstruction. Furthermore, the proposition for a Beringian pause [23,24,26] is largely based on the assumption that the American founder haplogroups differentiated during the Beringian pause/isolation. Several related issues remain to be addressed. First, despite the fact that much of interior Beringia remains in contemporary Alaska and northeast Siberia, no archaeological evidence of this population in residence has been found [34]. Second, a resident, presumably stable and possibly growing, population in Beringia seems an unlikely candidate for reduced genetic diversity in American founders as a result of founder effect. Third, while Chukchi and Siberian Eskimos share mtDNA lineages with Alaskan Eskimos, Aleuts and other Native North Americans [9,35], this sharing falls almost exclusively under haplogroups A2 and D2. The restricted genetic variation observed among these populations (only haplogroups A and D) seems inconsistent as a source area for the more extensive genetic variation observed in the rest of the Americas. Furthermore, the presence of haplogroup D2 in the Beringians is attributed to a second expansion from southern Siberia [36-39] and A2 is limited in diversity in Beringia [35,38].

In this study, we analyzed 568 mitochondrial genomes (of which 215 are newly sequenced from a Mexican American population) belonging to four pan-Native American (A2, B2, C1 and D1) lineages and A2a, A2b, A4b, C1a, C4 and D4e1 sister clades from Siberia-Asia, to better understand the underlying processes of the diversification of the Native American founders from their Asian counterparts and their expansion into the New World.

Results

Maternal legacy of Mexican Americans

Mexicans are, by and large, descendants of Native American and European (Spanish) ancestors [40]. Historical accounts also document African slavery in Mexico during the 16th-18th centuries [41-43], another source of admixture in the Mexican population. The admixture estimates compiled by Lisker et al. [44] using data derived from classical genetic systems reported in previous studies in Mexico [45-52] identified African and/or European genetic variation in all Mexican regions and groups analyzed. For mtDNA variation, some studies have measured Native American, European and African contributions to Mexican and Mexican American populations, revealing 85 to 90% of mtDNA lineages are of Native American origin [53,54], with the remainder having European (5-7%) or African ancestry (3-5%) [54]. Thus the observed frequency of Native American mtDNA in Mexican/Mexican Americans is higher than was expected on the basis of autosomal estimates of Native American admixture for these populations i.e. ~ 30-46% [53,55]. The difference is indicative of directional mating involving preferentially immigrant men and Native American women. This type of genetic asymmetry has been observed in other populations, including Brazilian individuals of African ancestry, as the analysis of sex specific and autosomal markers has revealed evidence for substantial European admixture that was mediated mostly through men [56]. In our 384 completely sequenced Mexican American mitochondrial genomes, 12 (3.1%) are of African ancestry belonging to haplogroups L0a1a'3', L2a1, L3b, L3d and U6a7; 52 (13.6%) belong to European haplogroups HV, JT, U1, U4, U5; and K and the majority (320, 83.3%) are of Native American ancestry, which is very similar to previous reports [53,54].

After removing the related individuals from our maternal lines and those of European or African descent, 215 newly sequenced unrelated mtDNA genomes of Native American ancestry were identified for founder haplogroups A2, B2, C1, D1 and D4e1c (new sister lineage of D2 found in Mexican Americans). These complete mtDNA sequences were analyzed together with 353 previously published sequences belonging to the four primary Native American founder haplogroups A2, B2, C1 and D1 and their Asian sister clades (A4, C4 and D4e; see additional file 1 for details).

Additional file 1. Details of mitochondrial genomes used in the present study. The file contains a table showing ethnographic details, GenBank accession numbers and references of the mitochondrial genome sequences used in the present study.

Format: DOC Size: 843KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Phylogenetic reconstructions and coalescent age estimates

The phylogenetic reconstruction and coalescent age estimates of the aforesaid 568 mitochondrial genomes is presented in additional file 2 (panels A-D) and summarized in Figure 1, Figure 2 and table 1.

Additional file 2. The detailed mtDNA phylogenetic tree of Native American haplogroups and their immediate Siberian-Asian sister clades (panel A-D). The file contains figures showing phylogenetic reconstruction of 568 complete mitochondrial DNA sequences (Maximum parsimony tree), belonging to Native American haplogroups and Siberian-Asian sister clades presented in four panels "A-D" (see additional file 3 for more details).

Format: XLS Size: 316KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 3. Figure Legends for additional file 2. The file contains details on phylogenetic reconstruction, sequences and mtDNA mutation rates used for coalescent age estimates in additional file 2.

Format: DOC Size: 27KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

thumbnailFigure 1. Schematic representation of mtDNA phylogenetic tree of Native American haplogroups A2 and B2 and immediate Siberian-Asian sister clades (A2a, A2b, A4a, A4b and A4c). Coalescent age calculated in thousand years (ky) as per the slow mutation rate of Mishmar et al. [58] and as per calibrated mutation rate of Soares et al. [59] are indicated in blue and red color respectively. The founder age wherever calculated are italicized. The geographical locations of the samples are identified with colors. For more details see complete phylogenetic reconstruction in additional file 2 (panels A-B) and additional file 3.

thumbnailFigure 2. Schematic representation of mtDNA phylogenetic tree of Native American haplogroups C1, C4c, D1, and D4e1c and their immediate Siberian-Asian sister clades (C1a, C4a, C4b, D2 and D4e1a). Coalescent age calculated in thousand years (ky) as per the slow mutation rate of Mishmar et al. [58] and as per calibrated mutation rate of Soares et al. [59] are indicated in blue and red color respectively. The founder age wherever calculated are italicized. The geographical locations of the samples are identified with colors. For more details see complete phylogenetic reconstruction in additional file 2 (panels C-D) and additional file 3.

Table 1. Haplogroup coalescence time estimates

Haplogroup 'A2' and its sister clades

A total of 134 American sequences (49 published earlier and 85 from this study) belonging to pan-American haplogroup A2, form 42 distinct haplotypes, one nested within the Beringian subtype A2a and the other 41 largely autochthonous to geographical North, Central and South America, except Central American which shares a few branches with both North and South Americans. 12 of these 41 haplotypes were previously labeled A2c through A2q [57]. However, based on our phylogenetic reconstruction of the haplotypes, we changed the previously reported definitions of A2e, A2f and A2n. Clade A2e is now defined by only one substitution at (13708), A2f sequences share a back mutation at np 153, and A2n has a back mutation at 16111 (but not at 153), although it should be noted that the control region mutations are recurrent. The other new haplotypes with three or more sequences have been labeled A2r through A2w.

Of the two Beringian specific branches (A2a and A2b), A2a is the only sublineage shared by the Chukchi, Siberian and North American Eskimos, Aleuts, and Na-Dene populations. The A2a and A2b sublineages might have diverged from the rest of the American subtypes sometime between 17 to 9 kya, as suggested by the founder coalescent ages of A2a and A2b and the coalescent age of haplogroup A2 calculated using the mutation rates of Mishmar et al. [58] and Soares et al. [59]. The coalescent age of haplogroup A2 calculated using these two mutation rates over America specific lineages (19.5 ± 1.3 kya/16.1 ± 1.5 kya) indicates its early expansion into America. The haplogroup A2 shares the root A4 with the Siberia-Asia specific sister clades A4a, A4b and A4c. The coalescent age estimate of A4, as per mutation rates of Mishmar et al. [58] and Soares et al. [59], are 27.5 ± 6.8 kya and 22.7 ± 7.4 kya respectively.

Haplogroup 'B2'

The phylogeny of haplogroup B2 reveals a total of 36 haplotypes, with 9 previously defined subhaplogroups named B2a through B2i [57], and two new lineages, B2j and B2k, defined here by control region substitutions at 16278 and 152, respectively. The remaining haplotypes were not classified because they consist of only two or fewer representative sequences. Similar to the other pan American haplogroups, sublineages within the B2 haplogroup are specific to geographical regions (North, Central and South America). Only a few branches harbor sequences across geographical regions. The sub branch B2b by virtue of its early differentiation (coalescent age 23.1 ± 4.8 kya/16.8 ± 4.9 kya calculated as per [58] and [59] respectively) encompasses the North, Central and South America. The coalescent age estimates of B2 calculated as per [58] and [59] mutation rates are 20.8 ± 2.0 kya and 18.1 ± 2.4 kya respectively. The haplogroup B2 lacks any known immediate sister/ancestral clades in Siberia other then the B4b1 in South East Asia.

Haplogroup 'C1' and 'C4'

The C root is represented by two lineages in America: the pan American haplogroup C1; and the recently identified subhaplogroup 'C4c', found in only two sequences, an Ijka sample from South America [23] and a Shuswap speaker from North America [60]. The C4c subhaplogroup shares root C4 with two wide spread lineages in Siberia-Asia, subhaplogroup C4a (coalescent age 25.7 ± 7.4 kya/21.6 ± 8.3 kya calculated as per [58] and [59] respectively) and C4b (coalescent age 9.3 ± 2.2 kya/6.8 ± 2.8 kya calculated as per [58] and [59] respectively). The coalescence age estimates of the C4c (20.6 ± 7.3 kya/11.8 ± 6.8 kya calculated as per [58] and [59] respectively) are based on only two known sequences and are therefore not very reliable; coalescent age of root C4 (21.0 ± 4.6 kya/20.0 ± 6.4 calculated as per [58] and [59] respectively) suggest its early differentiation.

In the case of the more abundant and widely distributed haplogroup C1, three basal branches are ubiquitous in the Americas (C1b, C1c and C1d), and spread across the three major geographic regions. The fourth branch, C1a, is found in Siberian and East Asian populations. The C1c branch exhibits an older coalescent age of 22.2 ± 3.3 kya/21.3 ± 4.3 kya (calculated as per [58] and [59] respectively), harboring sublineages (e.g., C1c3) that are spread across the three major geographical regions. The coalescent age of C1b is estimated to be 17.9 ± 2.3kya/18.4 ± 3.1kya (calculated as per [58] and [59] respectively) and for C1d, it is 20.1 ± 4.4 kya/13.5 ± 1.7 kya (calculated as per [58] and [59] respectively). The estimated coalescent age of Siberian branch C1a is 7.7 ± 3.2 kya/1.9 ± 1.9 kya (calculated as per [58] and [59] respectively) which is considerably younger. However its founder age is estimated to be 18.0 ± 7.9 kya/10.2 ± 8.1 kya which is comparable to the American branches (after allowing for standard error). The estimated coalescent ages of the roots for C4 (21.0 ± 4.6 kya/20.0 ± 6.4) and C1 (21.4 ± 2.7 kya/16.4 ± 1.5) calculated as per mutation rates of Mishmar et al. [58] and Soares et al. [59], along with their subclades C4a and C1c, respectively, indicate that these Native American clusters diverged from Siberian-Asian sister clades sometime during the LGM peak.

Haplogroup 'D1', 'D4e'and 'D4h'

The phylogenetic reconstruction of 19 new and 60 published mitochondrial genomes from America and Siberia-Asia reveals that haplogroup 'D' is represented by at least three branches in America. The D1 haplogroup is found in high frequency throughout the Americas. D4h3a is rare but also widely distributed. We identified two Mexican American sequences belonging to a new sub branch named D4e1c, defined by the substitutions 14207 and 16092. This new branch has been defined in our tree based on the only two sequences, but an exception was made due to its phylogenetic importance. The D4e1c shares the root D4e1 (substitutions 3316 and 9536) with D2, which is a diverse lineage found from southern Siberia to Beringia [9,37,38,61-63] and D4e1a found in Japan [64]. Since D4e1c has only two known sequences, coalescent age estimates are not reliable. However coalescent age estimates of D4e1 24.1 ± 9.0 kya/17.9 ± 10.0 kya (calculated as per [58] and [59] respectively) suggests early divergence of D4e1c from its sister clades Siberian/Beringian D2 and Japanese D4e1a. The phylogeny of D1 reveals 25 haplotypes, all more or less autochthonous to respective geographical regions except sub branch D1f which harbors genomes from North, Central and South America. Based on its phylogeography and coalescent age (19.3 ± 5.9 kya/13.8 ± 5.9kya), D1f appears to have differentiated early. The overall coalescent age of D1 haplogroup is 17.2 ± 2.0 kya/14.9 ± 2.2 kya (calculated as per [58] and [59] respectively).

We did not find sequences belonging to rare haplogroup X2a in our sample, whereas the other rare haplogroup D4h3 is represented by only one sequence. A recent study [28], based on 22 X2a and 46 D4h3 complete mtDNA sequences concluded that although the coding region divergence of these haplogroups shows considerable overlap, their phylogeography in America is strikingly different. The haplogroup D4h3 is found along the Pacific coast; whereas X2a is restricted to northern North America and hence envisaged two different paths of arrival. The haplogroup D4h3 followed the Pacific coastal path just like other pan American haplogroups, whereas X2a might have arrived from Beringia through a path represented by the ice free corridor between the Laurentide and Cordilleran ice sheets [28]. We concentrated our further analysis on A2, B2, C1, C4c, D4e1c and D1 American haplogroups and their Siberian/Asian sister clades.

Diversity indices and neutrality tests

Since we were interested in detecting the genetic consequences of demographic events that occurred during and after the peopling of the double continent, we computed nucleotide diversity and different estimators of the parameter θ = 2 Neμ to compare the amount of within population variability in Beringia, North America, Central America and South America. Also, in order to investigate the possibility of later expansions/migrations of Native American populations as suggested previously [28,65], we also performed the analysis by categorizing the Native American samples with respect to the four pan American haplogroups in the respective geographic regions.

The similar estimates of nucleotide diversity and θπ in North America (0.001747, 27.0), Central America (0.001735, 26.8) and South America (0.001834, 28.3) and the significantly negative value of Tajima's D and Fu's Fs (Table 2) suggest almost simultaneous (over a period of approximately 2.0 ky from South to North America) colonization of the three geographical regions, with rapid population expansion differentiating into more or less regional branches across the haplogroups (Figure 1 and 2). However the comparatively higher value of θk (1338.1 and 917.5) and significantly negative Fu's Fs (-24.0, p = 0.000 and -23.8, p = 0.004) in North and Central America, respectively, may also suggest an excess of rare alleles/haplotypes. The value of θk (413.5) in South America is comparatively low when compared to North and Central America.

Table 2. Diversity Indices and Neutrality tests for geographical region Beringia, North America, Central America and South America based on mtDNA coding region (np 577 - 16023) sequences.

The low values of diversity indices (nucleotide diversity = 0.000959; θπ = 14.8) and non-significantly negative values of both neutrality tests (Tejima's D = 1.02; Fu's Fs = 0.448) in Beringia suggest very different demographic patterns to that of Native Americans. The analysis of within haplogroups diversity in the respective geographic region (Table 2), though not significantly different when considered with standard deviation, shows two contrasting patterns of distribution i.e. nucleotide diversity and θπ levels for A2 and B2 haplogroups that are highest in North America, followed by South America and the lowest observed for Central America. In contrast, haplogroups C1 and D1 show a cline from South America to North America. As evident in the few shared branches across geographical region (Figure 1 and 2), the recent gene flow from the North and South might have resulted in the low diversity indices and higher values of θk in the Central American populations. The high values of diversity indices (nucleotide diversity = 0.000615 and θπ = 9.5) along with high θk (415.4) and highly significant negative Fu's Fs (-21.7, p = 0.00) of haplogroup A2 in North America may indicate a secondary expansion of the A2 haplogroup following the path represented by the ice free corridor between the Laurentide and Cordilleran ice sheets. A similar path for the expansion of haplogroup X2a in America has been suggested previously [28].

Bayesian skyline plot analysis

The above results suggest a complex demographic history associated with the colonization of the three geographic regions and of the Americas as a whole. We used Bayesian skyline plot analysis [66] to visually illustrate the change in effective population size in Beringia and three geographical regions of America. In this analysis, Neτ (where Ne is the population of breeding females and τ is the generation length) is unaffected by the different proposals for the actual number of founder lineages among present day populations [66]. The Bayesian skyline plots presented in Figure 3 (panels A-D) show that Beringian populations (Chukchi and Eskimos including the Aleutian Islands) suffered a gradual population decline following their divergence from Native Americans and Siberian ancestors (~15-10 kya) and are establishing moderate growth ~2.5 kya. In contrast, a pronounced population growth starting ~21 kya (95% CI: 23-19 kya) in South America and ~19 kya (95% CI: 21-17 kya) in North and Central America, continuing up ~16-15 kya in South and Central America and a little prolonged ~14-12 kya in North America has been observed. As reported earlier [24], each region in our analysis identifies a moderate population bottle neck just before the manifold rapid growth, which probably suggests a time frame of Native American founder's divergence and differentiation from their Asian ancestors sometime around LGM and then rapid expansion into America.

thumbnailFigure 3. Bayesian skyline plot showing Effective population size (Neτ) through time derived from mtDNA coding region (np 577 - 16023) sequences. The thick solid line is the median estimate and the thin lines show the 95% highest posterior density limits. In this analysis Neτ, where Ne is the population of breeding females and τ is the generation length is plotted on Y axis. Time estimates obtained using a log-normal relaxed clock with the standard substitution rate of 1.26 × 10-8 sites per year [58] are plotted on X axis and limited to < 45 kya. Panel 'A' - The Bayesian skyline plot (m = 10) for North America (n = 75). Panel 'B' - The Bayesian skyline plot (m = 10) for Central America (n = 151, randomly selected from total 215 sequences). Panel 'C' - The Bayesian skyline plot (m = 10) for South America (n = 76). Panel 'D' - The Bayesian skyline plot (m = 10) for Beringia (n = 66).

Discussion

Our analysis of a large data set of 568 mitochondrial genomes (215 this study and 353 from published sources; see additional file 1 for details) from North, Central and South America, as well as Beringia and Siberia-Asia shows that the coalescent age estimates of shared roots calculated as per the mutation rates of Mishmar et al. [58] and Soares et al. [59] between Siberian-Asian and Native American lineages i.e. A4 (27.5 ± 6.8 kya/22.7 ± 7.4 kya), C1 (21.4 ± 2.7 kya/16.4 ± 1.5 kya), C4 (21.0 ± 4.6 kya/20.0 ± 6.4 kya), and D4e1 (24.1 ± 9.0 kya/17.9 ± 10.0 kya) average around LGM when considered together with their standard deviations (Figure 4 panel A and B). These coalescent estimates along with the phylogeography of pan American haplogroups and Siberian-Asian sister clades strongly suggest that American founders diverged from their Siberian-Asian progenitors sometime during LGM. It is most likely that the hostile conditions during the LGM (centered at ~21.0 kya and extending from at least 23.0 to 19.0 kya) [67], caused by extensive glaciations, with a large ice lobe extending eastwards across northern Siberia [68,69] might have forced some populations back into southern Siberian LGM refuge, plausibly in Altai-Sayan and/or mid-lower Amur region [35,38] and the American founders to the east into Beringia. A swift expansion further south into the American continent occurred right after the LGM [24,65], most plausibly around 20-16 kya (Figure 4 panel A and B), as suggested by the coalescent age estimates of pan American haplogroups (A2:19.5 ± 1.3 kya/16.1 ± 1.5 kya, B2:20.8 ± 2.0 kya/18.1 ± 2.4 kya, C1:21.4 ± 2.7 kya/16.4 ± 1.5 kya and D1:17.2 ± 2.0 kya/14.9 ± 2.2 kya) calculated as per the mutation rates of Mishmar et al. [58] and Soares et al. [59] and estimates of populations expansions using Bayesian Skyline approach (Figure 3 panels A-C). The estimated average time difference between the divergence of Native Americans from their Siberian-Asian precursors and the expansion of founding population into America is less than 5.0 ky irrespective of the mutation rate used (Figure 4 panel A-C). The hostile conditions during LGM are also evident in the compilation of the records of fossil woody material or other plant macrofossils across northern Siberia i.e. records become progressively rarer right after 25 kya and then are absent completely for several thousand years in the period that includes the LGM [70,71]. Given that the opening of the ice free corridor is dated not earlier than ~14 kya, a coastal (Pacific) route would have been the only option for initial expansion [24,31-33]. Our data also show that the expansion of pan American haplogroups within the double continent is far more complex, contrary to the uniformity assumed by the BIM models and its variants [23-26]. The earliest expansion observed in the population reaching South America support the pre-Clovis occupation of the New World via the coastal route right after LGM. The expansion in North America occurred after the LGM (~19 kya) indicating a still prevailing LGM condition in the north. These expansions, separated by time and space and perhaps involving pan American lineages differentially, are responsible for a complex pattern of mtDNA diversity within the haplogroups A2, B2, C1, and D1 in the three geographical regions of America. Some intra-haplogroup variation, at least differentiation of sub branches B2b, C1b, C1c, C1d and D1f, already existed in American founders before expansion into the Americas. The population bottleneck, also reported previously [24], though moderate during this period (LGM) and a many fold rapid expansion right after LGM (~19 kya) perhaps better explains the few founder types in America.

thumbnailFigure 4. Schematic presentation of the average divergence of the American clades from Siberian-Asian sister clades and their dispersal into the America. Horizontal bars indicate the coalescent age estimates (at the center) ± Standard errors for each haplogroup shown on the Y axis. The average divergence time and the average time of the peopling of America across the haplogroups are also superimposed. X axis shows the time in thousand years ago. The coalescent estimates were calculated by Rho (ρ) statistics and three different mutation rates (i) panel 'A': one base substitution (i.e. one mutation other than indel) in the coding region (577 - 16023) per 5, 140 years [58], (ii) panel 'B': calibrated mutation rate of Soares et al. [59] based on all synonymous substitutions, (iii) panel 'C': one synonymous transition per 6, 764 year [76]

The phylogeography and diversity within haplogroup C1 shows a cline south to north in America. However C1 is subdivided in the three pan American braches (C1b, C1c, C1d) [23,24,65,72] and a Siberian-Asian representative C1a [38,62,64,73]. The subhaplogroup C1a, which shares root substitutions with American sub branches (C1b, C1c and C1d), is found in Japan [64] and southern Siberia [38,62,73]. Finding C1a in America or Beringia could help to explain the presence of C1a in Siberia-Asia as a result of a back migration from Beringia or America, a scenario enforced in BIM by Tamm et al. [23] to explain origin of haplogroup C1 in Beringia and its phylogeography. Yet, to date this has not been achieved even after large sampling from this area [23,24,28,37-39,58,60-62,65,72-76]. Therefore, one still cannot rule out the most parsimonious scenario that haplogroup C1 originated in Siberia-Asia [38,62]. It is also worth noting that even though C1a is found across diverse populations in Siberia-Asia, it is rare; whereas C1, which differentiated into C1b, C1c and C1d, is found in high frequency in America, and its phylogeography and coalescent age in America is very similar to the other pan American haplogroups (A2, B2 and D1) [23,24,65]. Therefore disappearing C1a from the expanding population in America appears unlikely given that it would require disappearance of a lineage from an expanding population in an unpopulated territory, a demographic situation in which loss of lineage due to genetic drift is minimized [77]. Furthermore back migration of only C1a from Beringia to Siberia-Asia as suggested in BIM by Tamm et al. [23], when all the pan American haplogroup existed in a stable, plausibly expanding population in Beringia seems less parsimonious.

In addition to C1, another founder lineage 'C4c' has been defined recently in Native Americans [23,60], which shares root 'C4' with the diverse Siberian-Asian sister clades C4a and C4b. The coalescent age estimates of C1 (21.4 ± 2.7 kya/16.4 ± 1.5 kya) and C4 (20.7 ± 4.3 kya/20.0 ± 6.4 kya) calculated as per [58] and [59] further support a recent divergence of American founders from their Siberian-Asian progenitors.

The haplogroup A2 coalescent age (17.9 ± 1.2 kya/16.1 ± 1.5 kya calculated as per [58] and [59] respectively), phylogeography and diversity distribution in America and Beringia, has been used to suggest its in-situ origin in American founders plausibly in the Beringian-Alaskan region [38]. However this diversity is restricted to only two branches (A2a and A2b) in Beringia. Therefore differentiation of its precursor in Siberia-Asia or soon after the divergence of Native American founders from their Siberian-Asian precursor cannot be ruled out. The time window for such differentiation would have been any time between the age of its shared root A4 (27.5 ± 6.8 kya/22.7 ± 7.4 kya calculated as per [58] and [59] respectively) with Siberian-Asian sister clades (A4a, A4b, A4c) and its own American coalescent age (19.5 ± 1.3 kya/16.4 ± 1.5 kya calculated as per [58] and [59] respectively), or more plausibly ~21-19 kya marked by the first expansion of modern human into America [24,65]. The diversity distribution along with demographic parameter also indicates that there could be more than one expansion of A2 in North America plausibly following two different paths, one Pacific coastal path just like the other pan American haplogroups and another path represented by the ice free corridor between the Laurentide and Cordilleran ice sheets, similar to X2a [28,65]. The environmental and paleoecological data, indicates that such a path existed and was represented by the ice-free corridor between the Laurentide and Cordilleran ice sheets, which opened approximately 15 kya [78] or possibly was never completely closed [79]. Through such a corridor, small glacial-refuge areas have been recently identified [80].

The phylogeography and within haplogroup diversity of haplogroup B2 is very similar to haplogroup A2 and supports a major expansion from North America to South America. So far haplogroup B2 lacks any known immediate sister/ancestral clades in Siberia other then the B4b1 in South East Asia, and it remains unclear whether the B2, dating to 20.8 ± 2.0 kya/18.1 ± 2.4 kya, has evolved in situ in American founders after their divergence from Siberian-Asian ancestors or Siberia-Asia itself.

The D4 clusters coalescent age (28.3 ± 4.4 kya/24.8 ± 7.9 kya calculated as per [58] and [59] respectively), represents the divergence age of its sub branches, D1 found only in America and D4e1 and D4h3 having distributed from Asia to America, suggest a more complex process of differentiation. The within haplogroup diversity of haplogroup D1 showing a cline south to north in America and having at least one sub branch 'D1f', already differentiated before the expansion and numerous later regional sub branches strongly suggest its early differentiation and a sizable population with basal D1 prior to expansion into America. Contrary to the absence of any known progenitor in Siberia-Asia, a Siberian-Asian origin of D1 cannot be ruled out [35,62]. The D4e1 is represented by Japanese branch D4e1a in Asia [64], D2 in Siberia/Beringia and by our newly discovered D4e1c branch in America. D2 constitutes two different clusters with contrasting geographic distribution. The cluster D2a is found in high frequency among Chukchi, Aleuts, and Eskimos [36-39,61-63], whereas the second cluster D2b is reported in Buryat, Khamnigan and Yakut populations [38]. The phylogeography of D2 and the coalescent age estimate (14.6 ± 6.3 kya/10.9 ± 7.1 kya calculated as per [58] and [59] respectively) may suggest its Southern Siberian origin and later expansion to Beringia. But American clade D4e1c, though rare, known by only two sequences from our current work, diverge from the root D4e1 and the coalescent age of D4e1 (24.1 ± 9.0 kya/17.9 ± 10.0 kya calculated as per [58] and [59] respectively), strongly suggest that D4e1c was the part of initial colonization of Americas.

Conclusions

Our analysis of a data set of 568 mitochondrial genomes from North, Central and South America, as well as Beringia and Siberia-Asia suggest that American founders diverged from their Siberian-Asian progenitors sometime during LGM and expanded into America soon after the LGM peak (~20-16 kya). Further, time between the Native American divergence from Siberian-Asians to the expansion into America was shorter than the previous estimates. The phylogeography of haplogroup C1 suggest that this American founder haplogroup differentiated in Siberia-Asia. Although it is not clear for the haplogroup B2, haplogroups A2 and D1 might have differentiated soon after the Native American founder's divergence. A moderate population bottle neck in American founder populations just before the expansion (Figure 3) most plausibly resulted in few founder types in America. The similar estimates of the diversity indices and Bayesian skyline analysis in North America, Central America and South America suggest almost simultaneous (over a period of approximately 2.0 ky from South to North America) colonization of these geographical regions with rapid population expansion differentiating into more or less regional branches across the pan American haplogroups A2, B2, D1 and the C1b, C1c, C1d subhaplogroups of C1. However, some sub branches (B2b, C1b, C1c, C1d and D1f) already existed in American founder haplogroups before expansion into the Americas.

Methods

There have been considerable advances in the understanding of Native American mtDNA phylogeography and the peopling processes of the American continent in recent years; however this has been based largely on mitochondrial genomes sampled from North and South American populations, with very few coming from southern North America and Central America [23,24,58,65,74-76]. Therefore, in order to address this insufficiency, we have completely sequenced mtDNA of 384 Mexican American maternal lines from our San Antonio Family Heart Study (SAFHS). Most Mexican Americans are the descendants of the Indigenous peoples of Mexico and/or Europeans [40], especially Spaniards and Mexico is considered as the southern tip of North America in general. Considering the pre-modern societies such as the Maya, who thrived in one cultural sphere that spread from El Salvador, Cozumel and Central Mexico, and into Veracruz, Honduras and Guatemala [81]; we consider our population sample has strong ethnic ties to Central America and therefore for the analysis combined with very few Central American mitochondrial genomes published previously.

Population samples

Mexican Americans account for more than 12.5% of the United States' population: 30.7 million Americans listed their ancestry as Mexican as of 2008, forming about 64% of all Hispanics and Latinos in the United States [82,83].

We have sampled 1, 236 Mexican American individuals organized into 42 extended families from the SAFHS [84]. From these extended pedigrees we selected 384 pedigree founders, representing all maternal lines in the family sample, for complete mtDNA sequencing. All participants gave informed consent and all protocols were approved by the Institutional Review Board of the University of Texas Health Science Center at San Antonio (San Antonio, TX).

As reported previously [53,54] and also observed in our study, mtDNA sequence showed some of the Mexican American individuals were of European and African ancestry due to recent admixture in historic times. Additionally, some of the founder individuals sequenced were related (from the same maternal line: based on our pedigree information). To focus our analyses, we removed all such genomes from this study and analyzed our 215 newly sequenced unrelated mitochondrial genomes belonging to Native American ancestry along with 353 from published source belonging to Native American haplogroups (A2, B2, C1 and D1) and their Asian sister clades (A4, C4 and D4e) [23,24,38,58,60-62,64,65,72-76,85]. The mtDNA sequences of Herrnstadt et al. [75] and Kivisild et al. [76] lack control region information. The 63 published source genomes [72], which were sequenced selectively for subhaplogroup C1d, were excluded from diversity and Bayesian skyline analysis because preselected sequences of single subhaplogroup would have skewed the results of these demographic parameters. See additional file 1 for more details on the mtDNA sequences used in this study.

Molecular analysis

The DNA was extracted from collected white blood cell from 10 ml blood samples using standard phenol-chloroform methods [86] with minor modifications. For complete mtDNA resequencing we used the mitoSEQr™ kit (Applied Biosystems) designed to amplify human mitochondrial DNA in 46 fragments using universal polymerase chain reaction (PCR) conditions. However we had mixed success with some of the fragments, so the resequencing of such regions of the genome was completed using alternate PCR primers and conditions of Rieder et al. [87]. Successful PCR amplification was verified on 1% ethidium bromide stained agarose gels, and purified using ExoSap-IT to remove excess primers and deoxynucleotide triphosphates (GE Healthcare). The amplicons were then used as templates in cycle sequencing reactions using the BigDye Terminator v3.1 Cycle Sequencing Ready Reaction Kit (Applied Biosystems), according to the manufacturer's instructions. The cycle sequencing products were purified using BigDye XTerminator Purification Kit (Applied Biosystems) to remove unincorporated BigDye and primer and then analyzed on an Applied Biosystems 3730 DNA Analyzer. Both strands were sequenced to promote resolution of polymorphisms. Applied Biosystems supplied DNA Sequencing Analysis Software version 5.1.1 was used for first pass base calling quality assessment. Comprehensive contig assembly and sequence alignment was performed in SeqScape Software version 2.6 (Applied Biosystems). Mutations were scored relative to the revised Cambridge Reference Sequence (rCRS) [88] with each deviation confirmed by manual checking of electropherograms. All (n = 215) mtDNA complete genome sequences used in this study have been submitted to GenBank (accession numbers HQ012049-HQ012263).

Statistical analysis

Phylogeny reconstruction and age estimation

Besides our newly sequenced 215 unrelated mitochondrial genomes of Native American ancestry, 353 additional published sequences belonging to Native American haplogroups (A2, B2, C1 and D1) and their Siberian-Asian sister clades (A4, C4 and D4e) as detailed in additional file 1 were employed for tree reconstruction. The phylogenetic trees of each of the pan American mtDNA haplogroups and their sister clades were reconstructed from median joining networks rooted to L3, using NETWORK 4.2.0.1 software [89]. The trees were checked manually to resolve homoplasies. Since we lack control region information for previously published sequences of Herrnstadt et al. [75] and Kivisild et al. [76], we used slow evolving coding region information for coalescent age estimates. The coalescent age estimates were calculated by Rho (ρ) statistics [90] and three different mutation rates: (i) one base substitution (i.e. one mutation other than indel) in the coding region (577 - 16023) per 5, 140 years [58]; (ii) one synonymous transition per 6, 764 year [76]; and (iii) calibrated mutation rate of [59] based on all synonymous substitutions. All of these mutation rates are calibrated on the basis of an assumed human-chimp split. Standard errors for coalescence estimates were calculated as per Saillard et al. [90]. The coalescent age estimates using the three aforesaid mutation rates are presented in table 1. It has been observed in our analysis and previously [23] that coalescent estimates based on the effectively faster rate of Kivisild et al. [76], yields younger age for most Native American haplogroups than the well documented archaeological date of modern human occupation at Monte Verde site in southern South America [30]. The mutation rates of Mishmar et al. [58] and Soares et al. [59] yielded higher and similar coalescent ages which fit well with the archaeological estimates of modern human occupation at Monte Verde site in southern South America [30]. Therefore, we largely based our inferences on the coalescent age estimates of these two mutation rates.

Diversity indices and neutrality tests

We computed nucleotide diversity and different estimators of the parameter θ = 2 Neμ to compare the amount of within population variability in Beringia, North America, Central America and South America: (1) θk [91], based on the number of observed alleles (k), (2) θS [92], based on the number of observed segregating sites (S), and (3) θπ [93], based on the mean number of pair-wise differences between sequences (π). We used Tajima's D [94] and Fu's Fs [95] statistics to investigate population expansion events. To investigate any event of second expansion/migration affecting or involving haplogroups differentially, we also performed the analysis by categorizing the Native American samples with respect to four pan American haplogroups in the respective geographic regions. The aforesaid diversity indices and neutrality test were performed on the slow evolving mtDNA coding region (position 577-16023) using the software package ARLEQUIN version 3.1 [96].

Bayesian skyline plot analysis

We used Bayesian skyline model of effective population size [66] to visually illustrate the demographic history of Beringian and Native American populations from the Most Recent Common Ancestor (MRCA). Effective population size is a compound population genetic parameter generally considered linearly proportional to census population size, the population of breeding females in this analysis. It is influenced by many factors, including local extinction, re-colonization and various forms of nonrandom mating [97]. The model assumes regional isolation and no phylogenetic structure priori, such as the existence of haplogroups or the number of founding haplotypes. Estimates of effective populations over time for geographic regions Beringia (n = 66), North America (n = 75), Central America (n = 151; randomly selected from 215 to overcome computation limitation in BEAST), and South America (n = 76), were derived from coding region mtDNA sequences. We considered Chukchi and Eskimos sampled from Siberia, and Aleuts sampled from Commander Island as the representative Beringian populations in this analysis. The analysis was carried out assuming HKY+G model and a log-normal relaxed clock [98] with the standard substitution rate of 1.26 × 108 sites per year [58] using Markov Chain Monte Carlo (MCMC) [99] sampling with 10 groups (m = 10). The software packages BEAST v1.4.8 [100] and Tracer v1.4.1 [101] used in this analysis are available from http://beast.bio.ed.ac.uk/ webcite. The plots were obtained using a stepwise (constant) model. The analysis was run for 30 million iterations for Beringia, North and South America and 40 million iterations for Central America with the first 10% discarded as burn-in in each case, genealogies and model parameters were sampled at every 1000 iterations thereafter.

List of Abbreviations

LGM: Last Glacial Maximum; kya: Thousand Years Ago; mtDNA: Mitochondrial DNA; BIM: Beringian Incubation Model; SAFHS: San Antonio Family Heart Study; PCR: Polymerase Chain Reaction; rCRS: revised Cambridge Reference Sequence; MRCA: Most Recent Common Ancestor; MCMC: Markov Chain Monte Carlo.

Authors' contributions

SK, CB and JEC performed the complete mtDNA sequencing and sequence alignments. SK performed the database search and phylogenetic analysis. SK and JEC drafted the manuscript. MZ and PEM provided input into the analysis and helped to improve the manuscript. JB and JEC conceived the study, participated in its design and coordination, and also helped to improve the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We are grateful to the participants in the San Antonio Family Heart Study. Data collection was supported by a grant from the US National Institute for Heart, Lung and Blood (HL045222). Additionally, this work was supported by National Institutes of Health grant MH059490, and a research grant from the San Antonio Area Foundation. This investigation was conducted in facilities constructed with support from Research Facilities Improvement Program Grant Numbers C06 RR013556 and C06 RR017515 from the National Center for Research Resources, National Institutes of Health. The AT&T Genomics Computing Center supercomputing facilities used for this work was supported in part by a gift from the AT&T Foundation.

References

  1. Wallace DC, Garrison K, Knowler WC: Dramatic founder effects in Amerindian mitochondrial DNAs.

    Am J Phys Anthropol 1985, 68:149-155. PubMed Abstract | Publisher Full Text OpenURL

  2. Greenberg JH, Turner CG, Zegura SL: The settlement of the Americas: A comparison of the linguistic, dental and genetic evidence.

    Curr Anthropol 1986, 27:477-497. Publisher Full Text OpenURL

  3. Crawford MH: The Origin of Native Americans. Cambridge, UK: Cambridge University Press; 1998. OpenURL

  4. Goebel T: Pleistocene human colonization of Siberia and peopling of the Americas: An ecological approach.

    Evol Anthropol 1999, 8:208-227. Publisher Full Text OpenURL

  5. Brace CL, Nelson AR, Seguchi N, Oe H, Sering L, Qifeng P, Yongyi L, Tumen D: Old World sources of the first New World human inhabitants: A comparative Craniofacial view.

    Proc Natl Acad Sci USA 2001, 98:10017-10022. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Schurr TG: The peopling of the New World: perspectives from molecular anthropology.

    Ann Rev Anthropol 2004, 33:551-583. Publisher Full Text OpenURL

  7. Nichols J: Linguistic diversity and the first settlement of the New World.

    Language 1990, 66:475-521. Publisher Full Text OpenURL

  8. Lell JT, Brown MD, Schurr TG, Sukernik RI, Starikovskaya YB, Torroni A, Moore LG, Troup GM, Wallace DC: Y chromosome polymorphisms in native American and Siberian populations: Identification of native American Y chromosome haplotypes.

    Hum Genet 1997, 100:536-543. PubMed Abstract | Publisher Full Text OpenURL

  9. Starikovskaya YB, Sukernik RI, Schurr TG, Kogelnik AM, Wallace DC: mtDNA diversity in Chukchi and Siberian Eskimos: Implications for the genetic history of Ancient Beringia and the peopling of the New World.

    Am J Hum Genet 1998, 63:1473-1491. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Nettle D: Linguistic diversity of the Americas can be reconciled with a recent colonization.

    Proc Natl Acad Sci USA 1999, 96:3325-3329. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Fiedel SJ: The peopling of the New World: Present evidence, new theories, and future directions.

    J Archaeol Res 2000, 8:39-103. Publisher Full Text OpenURL

  12. Lell JT, Sukernik RI, Starikovskaya YB, Su B, Jin L, Schurr TG, Underhill PA, Wallace DC: The dual origin and Siberian affinities of Native American Y chromosomes.

    Am J Hum Genet 2002, 70:192-206. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Tarazona-Santos E, Santos FR: The peopling of the Americas: A second major 'migration?

    Am J Hum Genet 2002, 70:1377-1380. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Goebel T, Waters MR, Dikova M: The archaeology of Ushki Lake, Kamchatka, and the Pleistocene peopling of the Americas.

    Science 2003, 301:501-505. PubMed Abstract | Publisher Full Text OpenURL

  15. Schurr TG, Sherry ST: Mitochondrial DNA and Y chromosome diversity and the peopling of the Americas: Evolutionary and demographic evidence.

    Am J Hum Biol 2004, 16:420-439. PubMed Abstract | Publisher Full Text OpenURL

  16. Bandelt HJ, Herrnstadt C, Yao YG, Kong QP, Kivisild T, Rengo C, Scozzari R, Richards M, Villems R, Macaulay V, Howell N, Torroni A, Zhang YP: Identification of Native American founder mtDNAs through the analysis of complete mtDNA sequences: Some caveats.

    Ann Hum Genet 2003, 67:512-524. PubMed Abstract | Publisher Full Text OpenURL

  17. Schurr TG, Ballinger SW, Gan Y, Hodge JA, Merriwether DA, Lawrence DN, Knowler WC, Weiss KM, Wallace DC: Amerindian mitochondrial DNAs have rare Asian mutations at high frequencies, suggesting they derived from four primary maternal lineages.

    Am J Hum Genet 1990, 46:613-623. PubMed Abstract | PubMed Central Full Text OpenURL

  18. Merriwether DA, Hall WW, Vahlne A, Ferrell RE: mtDNA variation indicates Mongolia may have been the source for the founding population for the New World.

    Am J Hum Genet 1996, 59:204-212. PubMed Abstract | PubMed Central Full Text OpenURL

  19. Merriwether DA, Rothhammer F, Ferrell RE: Distribution of the four-founding lineage haplotypes in Native Americans suggests a single wave of migration for the New World.

    Am J Phys Anthropol 1995, 98:411-430. PubMed Abstract | Publisher Full Text OpenURL

  20. Forster P, Harding R, Torroni A, Bandelt HJ: Origin and evolution of Native American mtDNA variation: A reappraisal.

    Am J Hum Genet 1996, 59:935-945. PubMed Abstract | PubMed Central Full Text OpenURL

  21. Bonatto SL, Salzano FM: A single and early migration for the peopling of the Americas supported by mitochondrial DNA sequence data.

    Proc Natl Acad Sci USA 1997, 94:1866-1871. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Stone AC, Stoneking M: mtDNA analysis of a prehistoric Oneota population: Implications for the peopling of the New World.

    Am J Hum Genet 1998, 62:1153-1170. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Tamm E, Kivisild T, Reidla M, Metspalu M, Glenn Smith D, Mulligan CJ, Bravi CM, Rickards O, Martinez-Labarga C, Khusnutdinova EK, Fedorova SA, Golubenko MV, Stepanov VA, Gubina MA, Zhadanov SI, Ossipova LP, Damba L, Voevoda MI, Dipierri JE, Villems R, Malhi RS: Beringian standstill and spread of Native American founders.

    PLoS ONE 2007, 2:e829. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Fagundes NJ, Kanitz R, Eckert R, Valls AC, Bogo MR, Salzano FM, Smith DG, Silva WA Jr, Zago MA, Ribeiro-dos-Santos AK, Santos SEB, Petzl-Erler ML, Bonatto SL: Mitochondrial population genomics supports a single pre-Clovis origin with a coastal route for the peopling of the Americas.

    Am J Hum Genet 2008, 82:583-592. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Kitchen A, Miyamoto MM, Mulligan CJ: A three-stage colonization model for the peopling of the Americas.

    PLoS ONE 2008, 3:e1596. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Mulligan CJ, Kitchen A, Miyamoto MM: Updated three stage model for the peopling of the Americas.

    PLoS ONE 2008, 3:e3199. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Brown MD, Hosseini SH, Torroni A, Bandelt HJ, Allen JC, Schurr TG, Scozzari R, Cruciani F, Wallace DC: 1998. MtDNA haplogroup X: An ancient link between Europe/Western Asia and North America?

    Am J Hum Genet 1998, 63:1852-1861. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Perego UA, Achilli A, Angerhofer N, Accetturo M, Pala M, Olivieri A, Kashani BH, Ritchie KH, Scozzari R, Kong QP, Myres NM, Salas A, Semino O, Bandelt HJ, Woodward SR, Torroni A: Distinctive Paleo-Indian Migration Routes from Beringia Marked by Two Rare mtDNA Haplogroups.

    Curr Biol 2009, 19:1-8. PubMed Abstract | Publisher Full Text OpenURL

  29. Pitulko VV, Nikolsky PA, Girya EY, Basilyan AE, Tumskoy VE, Koulakov SA, Astakhov SN, Pavlova EYu, Anisimov MA: The Yana RHS site: Humans in the Arctic before the last glacial maximum.

    Science 2004, 303:52-56. PubMed Abstract | Publisher Full Text OpenURL

  30. Dillehay TD: The late Pleistocene cultures of South America.

    Evol Anthropol 1999, 7:206-216. Publisher Full Text OpenURL

  31. Fladmark KR: Routes: Alternate migration corridors for early man in North America.

    Am Antiq 1979, 44:55-69. Publisher Full Text OpenURL

  32. Dixon EJ: Quest for the Origins of the First Americans. Albuquerque, NM: University of New Mexico Press; 1993. OpenURL

  33. Dalton R: The coast road.

    Nature 2003, 422:10-12. PubMed Abstract | Publisher Full Text OpenURL

  34. O'Rourke DH, Raff JA: The human genetic history of Americas: Review the final frontier.

    Current Biology 2010, 20:R202-R207. PubMed Abstract | Publisher Full Text OpenURL

  35. Volodko NV, Starikovskaya EB, Mazunin IO, Eltsov NP, Naidenko PV, Wallace DC, Sukernik RI: Mitochondrial Genome Diversity in Arctic Siberians, with Particular Reference to the Evolutionary History of Beringia and Pleistocenic Peopling of the Americas.

    Am J Hum Genet 2008, 82:1084-1100. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Rubicz R, Schurr T, Babb P, Crawford MH: Mitochondrial DNA variation and the origins of the Aleutes.

    Hum Biol 2003, 75:809-835. PubMed Abstract | Publisher Full Text OpenURL

  37. Zlojutro M, Rubicz R, Devor EJ, Spitsyn VA, Makarov S, Wilson K, Crawford MH: Genetic structure of the Aleuts and Circumpolar populations based on mitochondrial DNA sequences: a synthesis.

    Am J Phy Anthropol 2006, 129:446-464. Publisher Full Text OpenURL

  38. Derenko M, Malyarchuk B, Grzybowski T, Denisova G, Dambueva I, Perkova M, Dorzhu C, Luzina F, Lee H, Vanecek T, Villems R, Zakharov I: Phylogeographic analysis of mitochondrial DNA in Northern Asian Populations.

    Am J Hum Genet 2007, 81:1025-1041. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Gilbert MTP, Kivisild T, Grønnow B, Andersen PK, Metspalu E, Reidla M, Tamm E, Axelsson E, Götherström A, Campos PF, Rasmussen M, Metspalu M, Higham TFG, Nathan R, Hoog C-JD, Koch A, Møller LN, Andreasen C, Meldgaard M, Villems R, Bendixen C, Willerslev E: Paleo-Eskimo mtDNA Genome Reveals Matrilineal Discontinuity in Greenland.

    Science 2008, 320:1787-1789. PubMed Abstract | Publisher Full Text OpenURL

  40. Enciso A, Fernando S: Los extranjeros en Mexico, la immigracion y el gobierno: Tolerancia o intolerancia religiosa?, 1821-1830.

    Historia Mexicana XLV(3):539-566. OpenURL

  41. Beltrán AG: The slave trade in Mexico.

    Hispanic Am Hist Rev 1944, 24:412-431. Publisher Full Text OpenURL

  42. Curtin PD: The Atlantic Slave Trade. Milwaukee: University of Wisconsin Press; 1969. OpenURL

  43. Muhammad JS: Mexico and Central America. In No longer invisible: Afro-Latin Americans today. Edited by Minority Rights Group. London: Minority Rights Publications; 1995:163-180. OpenURL

  44. Lisker R, Ramírez E, Babinsky V: Genetic structure of autochthonous populations of Meso-America: Mexico.

    Hum Biol 1996, 68:395-404. PubMed Abstract OpenURL

  45. Crawford MH, Leyshon WC, Brown K, Lees F, Taylor L: Human biology in Mexico. II. A comparison of blood group, serum and red cell enzyme frequencies, and genetic distances of the Indian populations of Mexico.

    Am J Phys Anthropol 1974, 41:251-268. PubMed Abstract OpenURL

  46. Crawford MH, Lisker R, Pérez-Briceño R: Genetic microdifferentiation of two transplanted Tlaxcaltecan populations. In The Tlaxcaltecans: prehistory, demography, morphology and genetics. Volume 7. Edited by Crawford MH. University of Kansas, Lawrence: Publications in anthropology; 1976::169-175. OpenURL

  47. Tiburcio V, Romero A, De Garay AL: Gene frequencies and racial intermixture in a Mestizo population from Mexico City.

    Ann Hum Biol 1978, 5:131-138. PubMed Abstract | Publisher Full Text OpenURL

  48. Lisker R, Babinsky V: Admixture estimates in nine Mexican Indian groups and five east coast localities.

    Rev Invest Clin 1986, 38:14-149. OpenURL

  49. Lisker R, Pérez-Briceño R, Granados J, Babinsky V: Gene frequencies and admixture estimates in the state of Puebla, Mexico.

    Am J Phys Anthropol 1988, 76:331-335. PubMed Abstract | Publisher Full Text OpenURL

  50. Cerda-Flores RM, Garza-Chapa R: Variation in the gene frequencies of three generations of humans from Monterrey, Nuevo León, Mexico.

    Hum Biol 1989, 61:249-261. PubMed Abstract OpenURL

  51. Lisker R, Ramírez E, Pérez-Briceño R, Granados J, Babinsky V: Gene frequencies and admixture estimates in four Mexican urban centers.

    Hum Biol 1990, 62:791-801. PubMed Abstract OpenURL

  52. Lisker R, Rameriz E, Gonzalez-Villapando C, Stern MP: Racial admixture in a Mestizo population from Mexico City.

    Am J Hum Biol 1995, 7:213-216. Publisher Full Text OpenURL

  53. Merriwether DA, Huston S, Iyengar S, Hamman R, Norris JM, Shetterly SM, Kamboh MI, Ferrell RE: Mitochondrial versus nuclear admixture estimates demonstrate a past history of directional mating.

    Am J Phys Anthropol 1997, 102:153-159. PubMed Abstract | Publisher Full Text OpenURL

  54. Green LD, Derr JN, Knight A: mtDNA affinities of the peoples of North-Central Mexico.

    Am J Hum Genet 2000, 66:989-998. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Beuten J, Halder I, Fowler SP, Goring HH, Duggirala R, Arya R, Thompson IM, Leach RJ, Lehman DM: Wide disparity in genetic admixture among Mexican Americans from San Antonio, TX.

    Ann Hum Genet 2011, 75:529-538. PubMed Abstract | Publisher Full Text OpenURL

  56. Bortolini MC, Silva WA, Castro-de-Guerra D, Remonatto G, Mirandola RM, Hutz MH, Weimer TA, Zago MA, Salzano FM: African-derived South American populations: a history of symmetrical and asymmetrical matings according to sex revealed by bi and uniparental genetic markers.

    Am J Hum Biol 1999, 11:551-563. PubMed Abstract | Publisher Full Text OpenURL

  57. van Oven M, Kayser M: Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation.

    Hum Mutat 2009, 30:E386-E394. PubMed Abstract | Publisher Full Text OpenURL

  58. Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark A, Hosseini S, Brandon M, Easley K, Chen E, Brown M, Sukerniki RI, Olckersj A, Wallace DC: Natural selection shaped regional mtDNA variation in humans.

    Proc Natl Acad Sci USA 2003, 100:171-176. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Soares P, Ermini L, Thomson N, Mormina M, Rito T, Rohl A, Salas A, Oppenheimer S, Macaulay V, Richards MB: Correcting for purifying selection: An improved human mitochondrial molecular clock.

    Am J Hum Genet 2009, 84:740-759. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Malhi RS, Cybulski JS, Tito RY, Johnson J, Harry H, Dan C: Brief communication: Mitochondrial haplotype C4c confirmed as the founding genome in the Americas.

    Am J Phy Anthropol 2010, 141:494-497. OpenURL

  61. Derbeneva O, Starikovskaia E, Volodko N, Wallace D, Sukernik R: Mitochondrial DNA variation in Kets and Nganasans and the early peoples of Northern Eurasia.

    Genetika 2002, 38:1554-1560. PubMed Abstract OpenURL

  62. Starikovskaya E, Sukernik R, Derbeneva O, Volodko N, Ruiz-Pesini E, Torroni A, Brown M, Lott M, Hosseini S, Huoponen K, Wallace D: Mitochondrial DNA diversity in indigenous populations of the southern extent of Siberia, and the origins of Native American haplogroups.

    Ann Hum Genet 2005, 69:67-89. PubMed Abstract | Publisher Full Text OpenURL

  63. Zlojutro M, Tarskaia LA, Sorensen M, Snodgrass JJ, Leonard WR, Crawford MH: Coalescent simulations of Yakut mtDNA variation suggest small founding population.

    Am J Phy Anthopol 2009, 139:474-482. Publisher Full Text OpenURL

  64. Tanaka M, Cabrera V, Conzales A, Larruga J, Takeyasu T, Fuku N, Guo LJ, Hirose R, Fujita Y, Kurata M, Shinoda K-I, Umetsu K, Yamada Y, Oshida Y, Sato Y, Hattori N, Mizuno Y, Arai Y, Hirose N, Ohta S, Ogawa O, Tanaka Y, Kawamori R, Shamoto-Nagai M, Maruyama W, Shimokata H, Suzuki R: Mitochondrial genome variation in Eastern Asia and the Peopling of Japan.

    Genet Res 2004, 14:1832-1850. OpenURL

  65. Achilli A, Perego UA, Bravi CM, Coble MD, Kong QP, Woodward SR, Salas A, Torroni A, Bandelt HJ: The phylogeny of the four pan-American mtDNA haplogroups: Implications for evolutionary and disease studies.

    PLoS ONE 2008, 3:e1764. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  66. Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian coalescent inference of past population dynamics from molecular sequences.

    Mol Biol Evol 2005, 22:1185-1192. PubMed Abstract | Publisher Full Text OpenURL

  67. Mix AC, Bard E, Schneider R: Environmental processes of the ice age: land, oceans, glaciers (EPILOG).

    Quaternary Science Reviews 2001, 20:627-657. Publisher Full Text OpenURL

  68. Grosswald MG: Late Quaternary ice sheet of Northern Eurasia.

    Quaternary Research 1980, 13:1-32. Publisher Full Text OpenURL

  69. Denton GH, Hughes TJ: The Last Great Ice Sheets. Wiley-Interscience; 1981. OpenURL

  70. Kremenetski CV: 1995: Shifts of ranges of spruce, linden and black alder in West Siberia and Kazakhstan in Late Glacial and Holocene.

    Palynology in Russia

    Moscow. Part. 2, p 56-71. (In Russian), sited from Adams and Faure, 1997.

    OpenURL

  71. Adams JM, Faure H: QEN members. Review and Atlas of Palaeovegetation: Preliminary land ecosystem maps of the world since the Last Glacial Maximum. [http://www.esd.ornl.gov/projects/qen/adams1.html] webcite

    Oak Ridge National Laboratory, TN, USA; 1997.

  72. Perego UA, Angerhofer N, Pala M, Olivieri A, Lancioni H, Kashani BH, Carossa V, Ekins JE, Gomez-Carballa A, Huber G, Zimmermann B, Corach D, Babudri N, Panara F, Myres NM, Parson W, Semino O, Salas A, Woodward SR, Achilli A, Torroni A: The initial peopling of the Americas: Agrowing number of founding mitochondrial genomes from Beringia.

    Genome Res 2010, 20:1174-1179. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  73. Ingman M, Gyllensten U: Rate variation between mitochondrial domains and adaptive evolution in humans.

    Hum Mol Genet 2007, 16:2281-2287. PubMed Abstract | Publisher Full Text OpenURL

  74. Ingman M, Kaessmann H, Paabo S, Gyllensten U: Mitochondrial genome variation and the origin of modern humans.

    Nature 2000, 408:708-713. PubMed Abstract | Publisher Full Text OpenURL

  75. Herrnstadt C, Elson J, Fahy E, Preston G, Turnbull D, Anderson C, Ghosh S, Olefsky J, Beal M, Davis R, Howell N: Reduced-median-network analysis of complete mitochondrial DNA coding-region sequences for the major African, Asian, and European haplogroups.

    Am J Hum Genet 2002, 70:1152-1171. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  76. Kivisild T, Shen P, Wall DP, Do B, Sung R, Davis K, Passarino G, Underhill PA, Scharfe C, Torroni A, Scozzari R, Modiano D, Coppa A, Knijff PD, Feldman M, Cavalli-Sforza LL, Oefner PJ: The role of selection in the evolution of human mitochondrial genomes.

    Genetics 2006, 172:373-387. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  77. Edmonds CA, Lillie AS, Cavalli-Sforza LL: Mutations arising in the wave front of an expanding population.

    Proc Natl Acad Sci USA 2004, 101:975-979. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Goebel T, Waters MR, O'Rourke DH: The late Pleistocene dispersal of modern humans in the Americas.

    Science 2008, 319:1497-1502. PubMed Abstract | Publisher Full Text OpenURL

  79. Fagan BM: The Great Journey: The Peopling of Ancient America, Updated Edition. Gainesville, FL: University Press of Florida; 2004. OpenURL

  80. Loehr J, Worley K, Grapputo A, Carey J, Veitch A, Coltman DW: Evidence for cryptic glacial refugia from North American mountain sheep mitochondrial DNA.

    J Evol Biol 2006, 19:419-430. PubMed Abstract | Publisher Full Text OpenURL

  81. Lohse JC, Valdez FJ: Ancient Maya Commoners. Austin, USA: University of Texas Press; 2004. OpenURL

  82. Pew Hispanic Centre: Country of origin profiles: 2008 American Community Survey. [http://pewhispanic.org/data/origins/] webcite

    Pew Hispanic Centre; 2008.

    Last accessed 2010

  83. Tafoya S: "Shades of Belonging". [http://pewhispanic.org/files/reports/35.pdf] webcite

    Pew Hispanic Center; 2004.

    Last accessed 2010

  84. Mitchell BD, Kammerer CM, Blangero J, Mahaney MC, Rainwater DL, Dyke B, Hixson JE, Henkel RD, Sharp RM, Comuzzie AG, VandeBerg JL, Stern MP, MacCluer JW: Genetic and environmental contributions to cardiovascular risk factors in Mexican Americans. The San Antonio Family Heart Study.

    Circulation 1996, 94:2159-2170. PubMed Abstract | Publisher Full Text OpenURL

  85. Kong QP, Yao YG, Sun C, Bandelt HJ, Zhu CL, Zhang YP: Phylogeny of East Asian mitochondrial DNA lineages inferred from complete sequences.

    Am J Hum Genet 2003, 73:671-676. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  86. Sambrook J, Fritsch E, Maniatis T: Molecular Cloning: A Laboratory Manual. New York: Cold Spring Harbor Laboratory; 1989. OpenURL

  87. Rieder MJ, Taylor SL, Tobe VO, Nickerson DA: Automating the identification of DNA variations using quality-based fluorescence re-sequencing: analysis of the human mitochondrial genome.

    Nucleic Acids Res 1998, 26:967-973. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  88. Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N: Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA.

    Nat Genet 1999, 23:147. PubMed Abstract | Publisher Full Text OpenURL

  89. Bandelt HJ, Forster P, Rohl A: Median-joining networks for inferring intraspecific phylogenies.

    Mol Biol Evol 1999, 16:37-48. PubMed Abstract | Publisher Full Text OpenURL

  90. Saillard J, Forster P, Lynnerup N, Bandelt HJ, Norby S: mtDNA variation among Greenland Eskimos: the edge of the Beringian expansion.

    Am J Hum Genet 2000, 67:718-726. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  91. Ewens WJ: The sampling theory of selectively neutral alleles.

    Theor Popul Biol 1972, 3:87-112. PubMed Abstract | Publisher Full Text OpenURL

  92. Watterson G: On the number of segregating sites in genetical models without recombination.

    Theor Popul Biol 1975, 7:256-276. PubMed Abstract | Publisher Full Text OpenURL

  93. Tajima F: Evolutionary relationship of DNA sequences in finite populations.

    Genetics 1983, 105:437-460. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  94. Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism.

    Genetics 1989, 123:585-595. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  95. Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and backgroud selection.

    Genetics 1997, 147:915-925. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  96. Excoffier L, Laval G, Schneider S: Arlequin (ver. 3.0): An integrated software package for population genetics data analysis.

    Evolutionary Bioinformatics Online 2005, 1:47-50. OpenURL

  97. Wakeley J: The effects of subdivision on the genetic divergence of populations and species.

    Evolution 2000, 54:1092-1101. PubMed Abstract OpenURL

  98. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed Phylogenetics and Dating with Confidence.

    PLoS Biology 2006, 4:e88. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  99. Metropolis N, Rosenbluth AW, Rosenbluth MN, Teller AH, Teller E: Equations of state calculations by fast computing machines.

    J Chem Phy 1953, 21:1087-1091. Publisher Full Text OpenURL

  100. Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees.

    BMC Evol Biol 2007, 7:214. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  101. Rambaut A, Drummond AJ: Tracer v1.4. [http://beast.bio.ed.ac.uk/Tracer] webcite

    2007.