Mitochondrial Echoes of First Settlement and Genetic Continuity in El Salvador

Abstract

Background

From Paleo-Indian times to recent historical episodes, the Mesoamerican isthmus played an important role in the distribution and patterns of variability all around the double American continent. However, the amount of genetic information currently available on Central American continental populations is very scarce. In order to shed light on the role of Mesoamerica in the peopling of the New World, the present study focuses on the analysis of the mtDNA variation in a population sample from El Salvador.

Methodology/Principal Findings

We have carried out DNA sequencing of the entire control region of the mitochondrial DNA (mtDNA) genome in 90 individuals from El Salvador. We have also compiled more than 3,985 control region profiles from the public domain and the literature in order to carry out inter-population comparisons. The results reveal a predominant Native American component in this region: by far, the most prevalent mtDNA haplogroup in this country (at ∼90%) is A2, in contrast with other North, Meso- and South American populations. Haplogroup A2 shows a star-like phylogeny and is very diverse with a substantial proportion of mtDNAs (45%; sequence range 16090–16365) still unobserved in other American populations. Two different Bayesian approaches used to estimate admixture proportions in El Salvador shows that the majority of the mtDNAs observed come from North America. A preliminary founder analysis indicates that the settlement of El Salvador occurred about 13,400±5,200 Y.B.P.. The founder age of A2 in El Salvador is close to the overall age of A2 in America, which suggests that the colonization of this region occurred within a few thousand years of the initial expansion into the Americas.

Conclusions/Significance

As a whole, the results are compatible with the hypothesis that today's A2 variability in El Salvador represents to a large extent the indigenous component of the region. Concordant with this hypothesis is also the observation of a very limited contribution from European and African women (∼5%). This implies that the Atlantic slave trade had a very small demographic impact in El Salvador in contrast to its transformation of the gene pool in neighbouring populations from the Caribbean facade.

About the Authors:

Antonio Salas

* E-mail: [email protected]

Affiliation: Unidade de Xenética, Departamento de Anatomía Patolóxica e Ciencias Forenses, and Instituto de Medicina Legal, Facultade de Medicina, Universidade de Santiago de Compostela, Galicia, Spain

José Lovo-Gómez

Affiliations Unidade de Xenética, Departamento de Anatomía Patolóxica e Ciencias Forenses, and Instituto de Medicina Legal, Facultade de Medicina, Universidade de Santiago de Compostela, Galicia, Spain, Laboratorio de Genética Forense, Instituto de Medicina Legal, Dr. Roberto Masferrer, Corte Suprema de Justicia, San Salvador, El Salvador

Vanesa Álvarez-Iglesias

Affiliation: Unidade de Xenética, Departamento de Anatomía Patolóxica e Ciencias Forenses, and Instituto de Medicina Legal, Facultade de Medicina, Universidade de Santiago de Compostela, Galicia, Spain

María Cerezo

Affiliation: Unidade de Xenética, Departamento de Anatomía Patolóxica e Ciencias Forenses, and Instituto de Medicina Legal, Facultade de Medicina, Universidade de Santiago de Compostela, Galicia, Spain

María Victoria Lareu

Affiliation: Unidade de Xenética, Departamento de Anatomía Patolóxica e Ciencias Forenses, and Instituto de Medicina Legal, Facultade de Medicina, Universidade de Santiago de Compostela, Galicia, Spain

Vincent Macaulay

Affiliation: Department of Statistics, University of Glasgow, Glasgow, United Kingdom

Martin B. Richards

Affiliation: Institute of Integrative and Comparative Biology, Faculty of Biological Sciences, University of Leeds, Leeds, United Kingdom

Ángel Carracedo

Affiliation: Unidade de Xenética, Departamento de Anatomía Patolóxica e Ciencias Forenses, and Instituto de Medicina Legal, Facultade de Medicina, Universidade de Santiago de Compostela, Galicia, Spain

Introduction

El Salvador lies on the Pacific coast (without an Atlantic seaboard) and it is the smallest of the Central American countries. Most of the country rests on a fertile volcanic plateau. It is segmented by two volcanic ranges running roughly west to east, separated by broad, fertile valleys, such as that of the river Lempa. El Salvador was inhabited by Native American groups who were in part descendants of the Aztecs and Toltec of Mexico, such as the Pipil (a Nahua tribe) and the Lenca. These two Native American communities inhabited mainly the western regions, constituting about 60% of the population throughout the colonial era and into the early decades of independence [1], [2].

The development of coffee estates led to the slow but continuous dissolution of most of the communal lands of Native villages [1], [2]. Thus, the 1930 census, the last to contain the category, designated only 5.6% of the population as “Indian” – although it is not clear what criteria were used in arriving at this figure. Other independent estimates (considering religious activities, distinctive women's dress, language, and involvement in various handicrafts) placed the mid-twentieth-century Indian population at 20% (∼400,000 persons). The abandonment of Indian language and customs was hastened by political repression; most natives stopped wearing traditional dress, abandoned the Pipil language, and adopted ladino customs. By 1975 no more than ∼1% of the population wore distinctive Indian clothing or followed Indian customs. Nowadays, the official language in El Salvador is Spanish, although Nahua is still spoken among some natives.

Although the American continent has been the target of many forensic and population genetic studies, there are nevertheless many American regions, such as El Salvador, that remain genetically uncharacterized. The mtDNA molecule is commonly used in anthropological contexts because of particular features (maternal inheritance, lack of recombination and high average mutation rate) that confer great power for phylogenetic and phylogeographic inferences. Many mtDNA studies of Native Americans have, however, been limited to genotyping a handful of mtDNA coding region sites that simply distinguish the four major Native American mtDNA haplogroups, A2, B2, C1 and D1 (generally using RFLP typing); unfortunately, the information provided by these few SNPs is of limited value in forensic and population genetics.

Here we have sequenced the mtDNA control region in a sample from El Salvador in order to investigate to what extent the Native American component has survived the impact of European colonialism and the concomitant influx of African slaves to the Caribbean and Meso-America.

Materials and Methods

Sample collection and DNA extraction

A total of 90 saliva samples were collected from healthy unrelated individuals from El Salvador. DNA extraction was undertaken following standard phenol-chloroform protocol. DNA quantification was carried out using DyNA Quant 200 Fluorometer, Hoefer (APB, Uppsala, Sweden).

All the samples were collected anonymously by the Laboratorio de Genética Forense from the Instituto de Medicina Legal that belongs to the Corte Suprema de Justicia from San Salvador (El Salvador). Oral informed consent was required in all the cases. The study, including the oral informed consent protocol, was approved by the Ethical committee of the University of Santiago de Compostela, and it conforms to the Spanish Law for Biomedical Research (Law 14/2007- 3 of July).

PCR and sequence analysis

We analyzed the first and second hypervariable segments (HVS-I and HVS-II) of the mtDNA genome. We performed PCR amplifications using a 2700 Thermocycler (Applied Biosystems), using PCR and sequencing primers as reported in [3]. Cycling parameters were 95°C for 1 min, followed by 36 cycles of 95°C for 10 sec, 55°C for 30 sec and 72°C for 30 sec, and followed by 10 min at 15°C. We checked amplification products on a polyacrylamide gel visualized by silver staining and purified with Montage (Multiscreen PCR, Millipore Corporation, USA). We performed sequence reaction products on each strand by means of the ABI Prism dRhodamine Terminator cycle sequencing reaction kit (Applied Biosystems). DNA products were then purified by ethanol precipitation and sequence reaction products analyzed on the ABI Prism 3100 automatic sequencer (Applied Biosystems). We omitted population variation at the hypervariable sites (mainly related to the cytosine homopolymeric track around 310 and the CA-dinucleotide repeat around positions 522) from inter-population comparisons and phylogeographic analyses. We have used the same primers for amplification and sequencing described in [4].

Sequences were edited using the numbering system of the revised Cambridge Reference Sequence [5]. Most of the sequences could be read from np 16042–16569 and 21–550; for convenience, we will refer to these as HVS-I and HVS-II, although these sequence ranges encompass more than the canonical ranges of these control-region segments.

Quality checking

Problems with the quality of mtDNA data in forensic, clinical, and population genetic studies are unfortunately rather common; see, for instance, [6], [7], [8], [9], [10], [11]. In order to minimize the effects of potential laboratory and documentation errors, the data were read separately by two independent persons in the light of the known phylogeny. We checked phylogenetic inconsistencies by hand with special attention to private or unusual variants (e.g. rare transitions or indels). In some cases, we confirmed the sequences by repeated extraction and sequencing. In addition, to detect potential “phantom mutations” [7], we also checked the data using the computer program SPECTRA ([7], available at http://www.stats.gla.ac.uk/~vincent/fingerprint/index.html).

Statistical analysis and population comparison

Haplogroup nomenclature follows the most recently updated versions of the Native American phylogeny given in [12], [13], [14]. Diversity indices of HVS-I sequences (haplotype diversity, nucleotide diversity, average number of pairwise differences) were calculated using Arlequin 3.0 software [15]. Nucleotide and sequence diversity was computed as in the manner of Nei [16].

We estimated median-joining networks of HVS-I sequences using the Network 4.1.1.2 software [17], [18]. Coalescence times were calculated using the ρ statistic [19], [20] with an HVS-I mutation rate of one transition per 18,845 years applied for the sequence range 16090-16365 using the most recent estimates provided by Soares et al. [21].

An mtDNA database of Native American populations was compiled for population comparisons: (i) from North America: Aleuts [22], Athapaskans, Inupiaq, Yakima [23], Chukchi and Siberian Eskimos [24], Bella Coola and Haida [25], Nuu-Chah-Nulth [26], Cheyenne [27], North Native Americans [various ethnic groups; 28], [29], [30], [31], [32], [33], Apache and Navajo [34], (ii) from Meso-America: Pima [27], Maya [33], Huetar [35], Kuna [36], Ngöbe [37], Quiché [38], Emberá and Wounan [39], Mexico [40], Central Native Americans [various ethnic groups] [28], El Salvador (present study), and (iii) from South America: Native Brazilians and Araucanians or Chileans [33], Equador [41], Embera and Gavião [42], Amazonas [43], Ayoreo [44], Chilean Mapuche and Pehueche, Yaghan [45], Argentinian Mapuche [46], Cayapas [41], Xavante, Zoró, and Gavião [42], Yanomami [47], [48], South Native Americans [various ethinic groups] [29], Tuacuarembó [49], Uruguay [50], Guahibo [51], Colombia [33], [52], Yuracaré, Trinitario, Movima, and Ignaciano [53], and 105 from Arequipa, Tayacaja and San Martin in Peru [54]. We also included the data collected from several studies of ancient DNA [55], [56], [57], [58], [59], [60], [61]. In addition, other datasets were additionally used for haplotype matching comparisons [34], [62], [63], [64], [65], [66]. In total, 3,843 mtDNAs profiles (mainly HVS-I segments) were compiled for comparison with our sample from El Salvador. Those population samples consisting of less than 15 individuals were only used for haplotype matching between populations. For comparison purposes the common reading frame 16090–16365 of the HVS-I was used.

Admixture analysis

We took two different approaches to carry out an admixture analysis.

The first model was applied as described in [67] although, instead of using haplogroup frequencies as variables, we used the frequencies of the shared haplotypes (matching haplotypes) between the source populations (North and South America) and El Salvador. The number of mtDNAs within each matching haplotype in El Salvador (ni: 1≤i≤C, the number of different matching haplotypes in the sample) was assumed to be a draw from a multinomial distribution with parameters and (1≤i≤C), where R is the number of source regions in America, fji is the frequency of the ith cluster in the jth source region (assumed to be known), and αj are the admixture coefficients. This model describes samples from an urn containing C different kinds of ball, where the urn has been created by mixing together R other urns in proportions given by the admixture coefficients. We chose to analyze this model in a Bayesian framework, which meant that we had to explore the distribution of the admixture coefficients, given the data. The prior distribution of the admixture coefficients was taken to be uninformative—namely, uniform on . The posterior distribution of the {αj} was explored with the Metropolis-Hastings algorithm, using a simple proposal, and was summarized by the posterior mean of each aj and its root-mean-square deviation about the mean. To assess model fit, we examined plots of standardized residuals.

The second admixed model was applied as described [68]. The probability of origin of each of the sub-continental region was computed as where, n is the number of El Salvador sequences with matches (≥1) in the whole continental dataset; ki, the number of times the sequence i is found in El Salvador; pis, the frequency of the sequence i in the sub-continental region dataset; and pic, the frequency of the sequence i in whole continental dataset.

Founder analysis

The time to the most recent common ancestor (TMRCA) of haplogroup A2 in the phylogeny was estimated as described [19], [20].

In order to carry out a founder analysis [19], [69], we made some simplifying assumptions about the founding of El Salvador. We assumed (i) a single migration to El Salvador and (ii) that North America was the unique source population. Founder sequences were inferred as matches with samples from North America. An estimate of the time of the migration event was determined by averaging diversity over the clusters derived from each founder in El Salvador, as follows. Suppose there are r founder clusters. Let, ρi be the ρ value (average distance of the haplotypes of a clade from the respective root [19], [20]) for the ith founder cluster, σi be its estimated standard error [20] and ni be the number of sampled individuals in that cluster. Define,Then,Values of ρ and σ were converted to age using the most recent mutation rate available for the HVS-I segment of 1 transition per 18,845 years (in the sequence range 16090–16365).

Results

Summary statistics

We observed a total of 55 different HVS-I, 40 different HVS-II, and 76 different combined HVSI/II mtDNA haplotypes in El Salvador (N = 90). Some HVS-I profiles are quite common, such as C16111T–T16223C–C16290T–G16319A–T16362T, appearing in 12 mtDNAs, and its one step-mutation ‘neighbour’ haplotype (16519 on top) appearing 12 times.

As shown in Table 1, El Salvador shows haplotype and nucleotide diversity values slightly lower than those observed in the continental North, South, and other Meso American populations, which is in part due to the fact that there is virtually only one Native American haplogroup (A2) represented in El Salvador sample. Note that these comparisons have to beviewed with care because the terms “North”, “South” and “Meso-American” refer to groups of population samples of different nature; some are Native American groups that have passed through severe prehistoric bottlenecks while others are at different levels of admixture with e.g. Europeans and Africans.

[Figure omitted. See PDF.]

Table 1. Native American population mtDNA database considering the sequence range 16090–16365.

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

Phylogeography of Salvadorian Native American mtDNAs

Table 2 shows the full list of control region profiles from El Salvador and their haplogroup allocation. Frequencies of the typical Native American haplogroups A2, B2, and C1 are ∼91%, ∼2%, and ∼2%, respectively.

[Figure omitted. See PDF.]

Table 2. MtDNA haplotypes in El Salvador.

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

Figure 1 shows the frequency distribution of the main mtDNA American haplogroups in Native American populations. Although haplogroup A2 is at high frequencies in Meso America, El Salvador is particularly distinct from the other populations by its extremely high A2 haplogroup frequency. Note also that there exists substantial heterogeneity of haplogroup frequency patterns in America (even between neighbouring populations).

[Figure omitted. See PDF.]

Figure 1. Haplogroup patterns in America.

Codes for populations are as follow: North America: 1 = Chukchy, 2 = Eskimos [24]; 3 = Inuit (collected from the HvrBase database [84]; 4 = Aleuts [22]; 5 = Athapaskan [23]; 6 = Haida [25]; 7 = Apache [85], 8 = Bella Coola [25]; 9 = Navajo [85]; 10 = Sioux, 11 = Chippewa [66], 12 = Nuu-Chah-Nult [26]; 13 = Cheyenne [66]; 14 = Muskogean populations [31]; 15 = Cheyenne-Arapaho [30]; 16 = Yakima [23]; 17 = Stillwell Cherokee [30]; Meso-America: 18 = Pima [27]; 19 = Mexico [40]; 20 = Quiche [38]; 21 = Cuba [38]; 22 = El Salvador (present study); 23 = Huetar [35]; 24 = Emberá [39]; 25 = Kuna [36]; 26 = Ngöbé [37]; 27 = Wounan [37]; South America: 28 = Guahibo [51]; 29 = Yanomamo from Venezuela [48]; 30 = Gaviao [42]; 31 = Yanomamo from Venezuela and Brazil [86]; 32 = Colombia [52]; 33 = Ecuador (general population), 34 = Cayapa [41]; 35 = Xavante [42]; 36 = North Brazil [43]; 37 = Brazil [64]; 38 = Curiau [60]; 39 = Zoró [42]; 40 = Ignaciano, 41 = Yuracare [53]; 42 = Ayoreo [44]; 43 = Araucarians [33]; 44 = Pehuenche, 45 = Mapuche from Chile [45]; 46 = Coyas [4]; 47 = Tacuarembó [49]; 48 = Uruguay [50]; 49 = Mapuches from Argentina [46]; 50 = Yaghan [59].

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

The phylogeny of A2 in El Salvador is clearly star-like (Figure 2); its root is, identified by the diagnostic sites C16111T–T16223C–C16290T–G16319A–T16362C in HVS-I, and C64T–A73G–T146C–A153G –A235G–A263G–315+C in HVS-II. There are no very solid diagnostic sites in the control region that would allow us to classify A2 sub-lineages from El Salvador [12], [14]. Moreover, several control-region variants regarded as haplogroup diagnostic, such as C64T and A153G, show reversions: complete genome sequence data confirm the existence of multiple back and parallel mutations within haplogroup A2 [12], [14]. Although many of them are well-known hotspots (e.g. T146C), others such as position 64, seem to behave as hotspots only within A2 (see e.g. [12]). Other Native American lineages, like D1, D4h3 and X2a [13] are absent from our sample from El Salvador.

[Figure omitted. See PDF.]

Figure 2. Phylogenetic network of haplogroup A2 mtDNA sequences from El Salvador.

Only the variation contained in the HVS-I segment (range 16090–16365) was used. Positions are referred to the rCRS minus 16000; a transversion is specifically indicated as a suffix. The areas of the circles are proportional to the number of individuals bearing the corresponding haplotype in the El Salvador sample. The histograms indicate the frequency of the El Salvadorian mtDNA founder clusters in other continental regions (North, Meso, and South America). To compute the TMRCA and founder ages, we resolved the network to a tree (shown by the bold connections between haplotypes), using the positional mutation rates reported in [21]. Resolved parallel mutations are shown underlined. The ancestral haplotype is indicated by a hash.

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

The sub-clade of A2 carrying C16360T is particularly prevalent in Meso America, especially in the Huetar (12 matches; ∼44% of the Huetar sample) from Costa Rica [35] and the Ngöbé (three matches; ∼7% of the Ngöbe sample) from Panama [37]; in El Salvador this variant was also present in two individuals. The haplotype C16111T–C16187T–T16223C–C16290T–G16319A–T16362C is virtually only shared with the Ngöbe (19 matches that make up ∼41% of the Ngöbe sample) but was also detected in one Uruguayan [50]. El Salvador shares a higher number of haplotypes with North America (19), followed by Meso-America (10), and then South America (8); note however that the database for Meso-America (n = 395) is of a much lower sample size than the one from North (n = 2,010) and South (n = 1,596). These results roughly indicate a clear imprint of North in Meso-America and also the existence of lineages that are mainly concentrated in Meso-America (probably due to the fact that these were founders in the region and experienced posterior expansion); in some instances, some of these South mtDNAs could have been carried from Meso-America in some wave of migration towards the South (such as the one indicated above observed in Uruguay, or e.g. C16111T–T16189C–T16223C–C16290T–16311–G16319A–T16362C, which was also found in three Brazilians [64]).

We found only four Native American mtDNAs not belonging to haplogroup A2: two haplogroup B2 and two haplogroup C1 mtDNAs. We did not find any exact match amongst published data for the B2 sequence #33 that carries the distinctive variant A16269G. The haplogroup C1 sequence #41 carries C16256T; this uncommon variant within haplogroup C1 has been also observed in the Yanomama from Venezuela [47] and the Zoró from Brazil [42]. Haplotype C1 #39 was only observed in one Brazilian [43] and one Guahibo from Venezuela [51], but also in two ancient Taino samples from the Caribbean [57].

Non-Native American haplotypes in El Salvador

Signals of a European contribution to our sample from El Salvador are limited to three haplotypes (see Table 1): haplotype #96 belongs to haplogroup U2e, with exact matches in several West European locations (e.g. Northwest Spain and Portugal [70]); in Madeira [71], etc. Haplotype #12 can be assigned most plausibly to haplogroup H, while #30 probably belongs to haplogroup T2; this sequence curiously matches published sequences only observed in Portugal and Brazil [64], [72] but also a single hit in Poland [73].

We detected only one sequence belonging to a typical sub-Saharan haplogroup in El Salvador. It belongs to L0a1a, a sub-clade highly prevalent in southeast Africa [67], [74], [75], where we find exact matches in HVS-I and HVS-II. Exact matches are also found in, for example, the Atlantic African southwest coast, in Cabinda [76], and in the Tongas [77]. Although it is not possible to determine with precision the African origin of this haplotype, southeast Africa (Mozambique) is probably the best candidate population source.

Admixture analysis

The admixture analysis carried out as in [67] indicated that North America accounts for ∼92% of the lineages in El Salvador, the remaining ∼8% coming from South America. The method described in [68] indicated that North America contributed to El Salvador ∼76% of the mtDNA lineages, in contrast to the ∼24% coming from South America.

Founder analysis

We inferred seven founders in our Salvadorian sample, all present in North American populations. Some sequence matches were not considered founders because they were detected only in Mexico and not in North American populations; they are more likely the result of recent gene flow between El Salvador and neighbouring populations. Some other potential founders were also rejected because they were present mostly as singletons analyzed in North American laboratories but belonging to e.g. ‘Hispanic’ or other sample populations without information about their ethnic affiliation. All the founders are indicated in Figure 1. The founder age of haplogroup A2 in El Salvador was estimated as 13,400±5,200 Y.B.P.

Discussion

El Salvador is the smallest Latin American republic and also the most densely populated. Although historically El Salvador has been home to a culturally diverse mix of peoples, including Native Americans, Africans, and west Europeans, by the 1980s the population of the country was essentially considered to be homogeneous in terms of ethnicity and basic cultural identity. Virtually all Salvadorans speak Spanish, the official language, as their mother tongue, and the vast majority are generally characterized as “mestizos” (or “ladinos”, a term more commonly used in Central America), popularly used to refer to those persons of mosaic geographic ancestry who follow a wide variety of indigenous and “hispanic” customs and habits that over the centuries have come to constitute Spanish-American cultural patterns. In the late 1980s, the ethnic composition of the population was estimated as 89% “mestizo”, 10% Native American, and 1% “white” [78]. Therefore, in contrast to most other Central American countries, El Salvador no longer possessed an ethnically or linguistically distinct Native American population, although persons of native-like ethnicity or cultural heritage still lived in the western parts of the country. Similarly, there was no ethnically or culturally distinct African-American population as there is in neighbouring populations [79]. However, there is a general belief that much of the Salvadorian population in the 1980s had a predominantly Native American ancestry[1].

The results of the present study have shown that, in contrast to the cultural patterns observed in the today's El Salvador population, most of the mtDNA profiles found are typically Native American; haplogroup A2 account for ∼90% of the Salvadorian sample. Correspondingly, the impact of Europeans on the mtDNA pool of El Salvador is very low (∼2%). It seems that the Spanish conquerors and more recent European demographic influences did not contribute significantly to the today's genetic composition of El Salvador in the maternal side. This contrasts with the European Y-chromosome contribution to the El Salvador gene pool. According to [80] about one half in metropolitan areas and two thirds in rural populations of El Salvador belong to non-Native American haplogroups; for instance, the most common Y-chromosome haplogroup in Europe (namely, R1b) is present in El Salvador at 24% in metropolitan areas and 43% in rural regions. Concomitantly, the Native American Y-chromosome proportion in El Salvador (represented by haplogroup Q3) is about 31–49%.

Therefore, the mtDNA and Y-chromosome variation in El Salvador displays an extreme version of a pattern that was also observed in other American populations [81], [82]: the indigenous female contribution is much higher than the indigenous male contribution.

Our results show that the impact of African-American lineages on the mtDNA pool of El Salvador was very low, as indicated by the presence of only one mtDNA of sub-Saharan origin in our sample. The scarcity of the sub-Saharan component strikingly contrasts with the situation on the Caribbean coast, where (as a consequence of the Atlantic slave trade) it is clearly predominant [67], [74], [75], [79]. The Y-chromosome variation shows a similar pattern: no lineages of African ancestry have been detected in El Salvador [80].

There are no clear signals of recent genetic drift events in the general population from El Salvador, as observed in, for instance, neighbouring but isolated Native American populations such as the Ngöbé from Panamá [37] which shows extremely reduced levels of mtDNA diversity (reflecting passage through post-conquest population bottlenecks). Haplogroup A2 is at high frequency in El Salvador (∼90% of the sample) and a high percent of the lineages (45%; computed using the sequence range 16090–16365) remain unobserved in other American populations. Admixture analysis indicates that the main mtDNA influence in El Salvador can be attributed to North America. The phylogeny of A2 is rather star-like and the founder age was 12,600±4,900 years. The shape of this phylogeny points to the existence of a prehistoric demographic expansion. Considering the most recently estimated age of A2 in the American continent as a whole of 13,400±5,200 [21] (largely determined from North American samples) as a proxy for the time of the expansion into the Americas, it can be tentatively suggested that the initial settlement of El Salvador occurred rather soon after the initial colonization of the American continent, and that El Salvador largely contains the descendants of the mtDNAs in that original pool with scarce subsequent demographic influence from other American or non-American populations. Indeed, since we have genotyped samples collected in urban areas we would expect to have an even higher prevalence of the Native American component in more isolated groups from the country, as is in fact observed on the Y-chromosome side where the Native American component is higher in rural than in metropolitan areas [80].

In contrast to the high impact of the Atlantic slave trade on the Central American Caribbean coast [79], the Pacific side (at least for El Salvador) appears to have preserved its Native American mtDNA heritage intact to the present day. At the same time, this study has also shown that El Salvador harbours haplogroup frequency patterns quite different from other modern Native American communities. At the individual haplotype level, El Salvador shows numerous mtDNAs that have never been observed in other American regions, even within Central America. These features provide little support to those that assume (or claim) that “Hispanics” or Native American communities are sufficiently homogeneous to justify the portability of forensic databases from one country to another (e.g. SWGDAM; [34]); see [83] for a discussion.

Acknowledgments

We would like to thank Vilma de Aguilar (Ministerio de Salud Pública de El Salvador) for helping with the sample collection.

Author Contributions

Conceived and designed the experiments: AS. Performed the experiments: JLG VAI MC. Analyzed the data: AS VM MR. Contributed reagents/materials/analysis tools: AS MVL AC. Wrote the paper: AS VM MR.

Citation: Salas A, Lovo-Gómez J, Álvarez-Iglesias V, Cerezo M, Lareu MV, Macaulay V, et al. (2009) Mitochondrial Echoes of First Settlement and Genetic Continuity in El Salvador. PLoS ONE4(9): e6882. https://doi.org/10.1371/journal.pone.0006882

References

1. Ministerio de Educación ES (1994) Historia de El Salvador. ES Ministerio de Educación1994Historia de El Salvador

2. Lardé y Larín J (2000) J. Lardé y Larín2000 El Salvador: descubrimiento, conquista y colonización: Dirección de Publicaciones e Impresos, El Salvador. El Salvador: descubrimiento, conquista y colonización: Dirección de Publicaciones e Impresos, El Salvador.

3. Álvarez-Iglesias V, Salas A, Cerezo M, Ramos-Luis E, Jaime JC, et al. (2006) Genotyping coding region mtDNA SNPs for Asian and Native American haplogroup assignation. Int Congress Series 11: 4–6.V. Álvarez-IglesiasA. SalasM. CerezoE. Ramos-LuisJC Jaime2006Genotyping coding region mtDNA SNPs for Asian and Native American haplogroup assignation.Int Congress Series1146

4. Álvarez-Iglesias V, Jaime JC, Carracedo Á, Salas A (2007) Coding region mitochondrial DNA SNPs: targeting East Asian and Native American haplogroups. Forensic Sci Int: Genet 1: 44–55.V. Álvarez-IglesiasJC JaimeÁ. CarracedoA. Salas2007Coding region mitochondrial DNA SNPs: targeting East Asian and Native American haplogroups.Forensic Sci Int: Genet14455

5. Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, et al. (1999) Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet 23: 147.RM AndrewsI. KubackaPF ChinneryRN LightowlersDM Turnbull1999Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA.Nat Genet23147

6. Bandelt H-J, Salas A, Bravi CM (2004) Problems in FBI mtDNA database. Science 305: 1402–1404.H-J BandeltA. SalasCM Bravi2004Problems in FBI mtDNA database.Science30514021404

7. Bandelt H-J, Quintana-Murci L, Salas A, Macaulay V (2002) The fingerprint of phantom mutations in mitochondrial DNA data. Am J Hum Genet 71: 1150–1160.H-J BandeltL. Quintana-MurciA. SalasV. Macaulay2002The fingerprint of phantom mutations in mitochondrial DNA data.Am J Hum Genet7111501160

8. Salas A, Carracedo Á, Macaulay V, Richards M, Bandelt H-J (2005) A practical guide to mitochondrial DNA error prevention in clinical, forensic, and population genetics. Biochem Biophys Res Commun 335: 891–899.A. SalasÁ. CarracedoV. MacaulayM. RichardsH-J Bandelt2005A practical guide to mitochondrial DNA error prevention in clinical, forensic, and population genetics.Biochem Biophys Res Commun335891899

9. Salas A, Yao Y-G, Macaulay V, Vega A, Carracedo Á, et al. (2005) A critical reassessment of the role of mitochondria in tumorigenesis. PLoS Med 2: e296.A. SalasY-G YaoV. MacaulayA. VegaÁ. Carracedo2005A critical reassessment of the role of mitochondria in tumorigenesis.PLoS Med2e296

10. Yao Y-G, Salas A, Bravi CM, Bandelt H-J (2006) A reappraisal of complete mtDNA variation in East Asian families with hearing impairment. Hum Genet 119: 505–515.Y-G YaoA. SalasCM BraviH-J Bandelt2006A reappraisal of complete mtDNA variation in East Asian families with hearing impairment.Hum Genet119505515

11. Bandelt HJ, Salas A (2009) Contamination and sample mix-up can best explain some patterns of mtDNA instabilities in buccal cells and oral squamous cell carcinoma. BMC Cancer 16: 113.HJ BandeltA. Salas2009Contamination and sample mix-up can best explain some patterns of mtDNA instabilities in buccal cells and oral squamous cell carcinoma.BMC Cancer16113

12. Achilli A, Perego UA, Bravi CM, Coble MD, Kong QP, et al. (2008) The phylogeny of the four pan-American MtDNA haplogroups: implications for evolutionary and disease studies. PLoS ONE 3: e1764.A. AchilliUA PeregoCM BraviMD CobleQP Kong2008The phylogeny of the four pan-American MtDNA haplogroups: implications for evolutionary and disease studies.PLoS ONE3e1764

13. Perego UA, Achilli A, Angerhofer N, Accetturo M, Pala M, et al. (2009) Distinctive Paleo-Indian migration routes from Beringia marked by two rare mtDNA haplogroups. Curr Biol 19: 1–8.UA PeregoA. AchilliN. AngerhoferM. AccetturoM. Pala2009Distinctive Paleo-Indian migration routes from Beringia marked by two rare mtDNA haplogroups.Curr Biol1918

14. Tamm E, Kivisild T, Reidla M, Metspalu M, Smith DG, et al. (2007) Beringian standstill and spread of Native American founders. PLoS ONE 2: e829.E. TammT. KivisildM. ReidlaM. MetspaluDG Smith2007Beringian standstill and spread of Native American founders.PLoS ONE2e829

15. Excoffier L, Laval G, Schneider S (2005) Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinformatics Online 1: 47–50.L. ExcoffierG. LavalS. Schneider2005Arlequin ver. 3.0: An integrated software package for population genetics data analysis.Evol Bioinformatics Online14750

16. Nei N (1987) Molecular evolutionary genetics:. New York: Columbi University Press. N. Nei1987Molecular evolutionary genetics:New YorkColumbi University Press

17. Bandelt H-J, Forster P, Röhl A (1999) Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16: 37–48.H-J BandeltP. ForsterA. Röhl1999Median-joining networks for inferring intraspecific phylogenies.Mol Biol Evol163748

18. Bandelt H-J, Forster P, Sykes BC, Richards MB (1995) Mitochondrial portraits of human populations using median networks. Genetics 141: 743–753.H-J BandeltP. ForsterBC SykesMB Richards1995Mitochondrial portraits of human populations using median networks.Genetics141743753

19. Forster P, Harding R, Torroni A, Bandelt H-J (1996) Origin and evolution of Native American mtDNA variation: a reappraisal. Am J Hum Genet 59: 935–945.P. ForsterR. HardingA. TorroniH-J Bandelt1996Origin and evolution of Native American mtDNA variation: a reappraisal.Am J Hum Genet59935945

20. Saillard J, Forster P, Lynnerup N, Bandelt H-J, Nørby S (2000) mtDNA variation among Greenland Eskimos: the edge of the Beringian expansion. Am J Hum Genet 67: 718–726.J. SaillardP. ForsterN. LynnerupH-J BandeltS. Nørby2000mtDNA variation among Greenland Eskimos: the edge of the Beringian expansion.Am J Hum Genet67718726

21. Soares P, Ermini L, Thomson N, Mormina M, Rito T, et al. (2009) Correcting for purifying selection: an improved human mitochondrial molecular clock. Am J Hum Genet in press. P. SoaresL. ErminiN. ThomsonM. MorminaT. Rito2009Correcting for purifying selection: an improved human mitochondrial molecular clock.Am J Hum Genet in press

22. Rubicz R, Schurr TG, Babb PL, Crawford MH (2003) Mitochondrial DNA variation and the origins of the Aleuts. Hum Biol 75: 809–835.R. RubiczTG SchurrPL BabbMH Crawford2003Mitochondrial DNA variation and the origins of the Aleuts.Hum Biol75809835

23. Shields GF, Schmiechen AM, Frazier BL, Redd A, Voevoda MI, et al. (1993) mtDNA sequences suggest a recent evolutionary divergence for Beringian and northern North American populations. Am J Hum Genet 53: 549–562.GF ShieldsAM SchmiechenBL FrazierA. ReddMI Voevoda1993mtDNA sequences suggest a recent evolutionary divergence for Beringian and northern North American populations.Am J Hum Genet53549562

24. Starikovskaya YB, Sukernik RI, Schurr TG, Kogelnik AM, Wallace DC (1998) 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 63: 1473–1491.YB StarikovskayaRI SukernikTG SchurrAM KogelnikDC Wallace1998mtDNA diversity in Chukchi and Siberian Eskimos: implications for the genetic history of ancient Beringia and the peopling of the New World.Am J Hum Genet6314731491

25. Ward RH, Redd A, Valencia D, Frazier B, Pääbo S (1993) Genetic and linguistic differentiation in the Americas. Proc Natl Acad Sci U S A 90: 10663–10667.RH WardA. ReddD. ValenciaB. FrazierS. Pääbo1993Genetic and linguistic differentiation in the Americas.Proc Natl Acad Sci U S A901066310667

26. Ward RH, Frazier BL, Dew-Jager K, Pääbo S (1991) Extensive mitochondrial diversity within a single Amerindian tribe. Proc Natl Acad Sci U S A 88: 8720–8724.RH WardBL FrazierK. Dew-JagerS. Pääbo1991Extensive mitochondrial diversity within a single Amerindian tribe.Proc Natl Acad Sci U S A8887208724

27. Kittles RA, Bergen AW, Urbanek M, Virkkunen M, Linnoila M, et al. (1999) Autosomal, mitochondrial, and Y chromosome DNA variation in Finland: evidence for a male-specific bottleneck. Am J Phys Anthropol 108: 381–399.RA KittlesAW BergenM. UrbanekM. VirkkunenM. Linnoila1999Autosomal, mitochondrial, and Y chromosome DNA variation in Finland: evidence for a male-specific bottleneck.Am J Phys Anthropol108381399

28. Torroni A, Schurr TG, Cabell MF, Brown MD, Neel JV, et al. (1993) Asian affinities and continental radiation of the four founding Native American mtDNAs. Am J Hum Genet 53: 563–590.A. TorroniTG SchurrMF CabellMD BrownJV Neel1993Asian affinities and continental radiation of the four founding Native American mtDNAs.Am J Hum Genet53563590

29. Torroni A, Sukernik RI, Schurr TG, Starikovskaya YB, Cabell MF, et al. (1993) mtDNA variation of aboriginal Siberians reveals distinct genetic affinities with Native Americans. Am J Hum Genet 53: 591–608.A. TorroniRI SukernikTG SchurrYB StarikovskayaMF Cabell1993mtDNA variation of aboriginal Siberians reveals distinct genetic affinities with Native Americans.Am J Hum Genet53591608

30. Malhi RS, Schultz BA, Smith DG (2001) Distribution of mitochondrial DNA lineages among Native American tribes of Northeastern North America. Hum Biol 73: 17–55.RS MalhiBA SchultzDG Smith2001Distribution of mitochondrial DNA lineages among Native American tribes of Northeastern North America.Hum Biol731755

31. Bolnick DA, Smith DG (2003) Unexpected patterns of mitochondrial DNA variation among Native Americans from the southeastern United States. Am J Phys Anthropol 122: 336–354.DA BolnickDG Smith2003Unexpected patterns of mitochondrial DNA variation among Native Americans from the southeastern United States.Am J Phys Anthropol122336354

32. Lorenz JG, Smith DG (1994) Distribution of the 9-bp mitochondrial DNA region V deletion among North American Indians. Hum Biol 66: 777–788.JG LorenzDG Smith1994Distribution of the 9-bp mitochondrial DNA region V deletion among North American Indians.Hum Biol66777788

33. Horai S, Kondo R, Nakagawa-Hattori Y, Hayashi S, Sonoda S, et al. (1993) Peopling of the Americas, founded by four major lineages of mitochondrial DNA. Mol Biol Evol 10: 23–47.S. HoraiR. KondoY. Nakagawa-HattoriS. HayashiS. Sonoda1993Peopling of the Americas, founded by four major lineages of mitochondrial DNA.Mol Biol Evol102347

34. Monson KL, Miller KWP, Wilson MR, DiZinno JA, Budowle B (2002) The mtDNA Population Database: an integrated software and database resource for forensic comparison. Forensic Sci Commun 4: no 2.KL MonsonKWP MillerMR WilsonJA DiZinnoB. Budowle2002The mtDNA Population Database: an integrated software and database resource for forensic comparison.Forensic Sci Commun4no 2

35. Santos M, Ward RH, Barrantes R (1994) mtDNA variation in the Chibcha Amerindian Huetar from Costa Rica. Hum Biol 66: 963–977.M. SantosRH WardR. Barrantes1994mtDNA variation in the Chibcha Amerindian Huetar from Costa Rica.Hum Biol66963977

36. Batista O, Kolman CJ, Bermingham E (1995) Mitochondrial DNA diversity in the Kuna Amerinds of Panama. Hum Mol Genet 4: 921–929.O. BatistaCJ KolmanE. Bermingham1995Mitochondrial DNA diversity in the Kuna Amerinds of Panama.Hum Mol Genet4921929

37. Kolman CJ, Bermingham E, Cooke R, Ward RH, Arias TD, et al. (1995) Reduced mtDNA diversity in the Ngöbé Amerinds of Panamá. Genetics 140: 275–283.CJ KolmanE. BerminghamR. CookeRH WardTD Arias1995Reduced mtDNA diversity in the Ngöbé Amerinds of Panamá.Genetics140275283

38. Boles TC, Snow CC, Stover E (1995) Forensic DNA testing on skeletal remains from mass graves: a pilot project in Guatemala. J Forensic Sci 40: 349–355.TC BolesCC SnowE. Stover1995Forensic DNA testing on skeletal remains from mass graves: a pilot project in Guatemala.J Forensic Sci40349355

39. Kolman CJ, Bermingham E (1997) Mitochondrial and nuclear DNA diversity in the Choco and Chibcha Amerinds of Panama. Genetics 147: 1289–1302.CJ KolmanE. Bermingham1997Mitochondrial and nuclear DNA diversity in the Choco and Chibcha Amerinds of Panama.Genetics14712891302

40. Green LD, Derr JN, Knight A (2000) mtDNA affinities of the peoples of North-Central Mexico. Am J Hum Genet 66: 989–998.LD GreenJN DerrA. Knight2000mtDNA affinities of the peoples of North-Central Mexico.Am J Hum Genet66989998

41. Rickards O, Martinez-Labarga C, Lum JK, De Stefano GF, Cann RL (1999) mtDNA history of the Cayapa Amerinds of Ecuador: detection of additional founding lineages for the Native American populations. Am J Hum Genet 65: 519–530.O. RickardsC. Martinez-LabargaJK LumGF De StefanoRL Cann1999mtDNA history of the Cayapa Amerinds of Ecuador: detection of additional founding lineages for the Native American populations.Am J Hum Genet65519530

42. Ward RH, Salzano FM, Bonatto SL, Hutz MH, Coimbra CEA, et al. (1996) Mitochondrial DNA polymorphism in 3 Brazilian Indian tribes. Am J Hum Biol 8: 317–323.RH WardFM SalzanoSL BonattoMH HutzCEA Coimbra1996Mitochondrial DNA polymorphism in 3 Brazilian Indian tribes.Am J Hum Biol8317323

43. Santos SE, Ribeiro-Dos-Santos AK, Meyer D, Zago MA (1996) Multiple founder haplotypes of mitochondrial DNA in Amerindians revealed by RFLP and sequencing. Ann Hum Genet 60 (Pt 4): 305–319.SE SantosAK Ribeiro-Dos-SantosD. MeyerMA Zago1996Multiple founder haplotypes of mitochondrial DNA in Amerindians revealed by RFLP and sequencing.Ann Hum Genet60 (Pt 4)305319

44. Dornelles CL, Battilana J, Fagundes NJ, Freitas LB, Bonatto SL, et al. (2004) Mitochondrial DNA and Alu insertions in a genetically peculiar population: the Ayoreo Indians of Bolivia and Paraguay. Am J Hum Biol 16: 479–488.CL DornellesJ. BattilanaNJ FagundesLB FreitasSL Bonatto2004Mitochondrial DNA and Alu insertions in a genetically peculiar population: the Ayoreo Indians of Bolivia and Paraguay.Am J Hum Biol16479488

45. Moraga ML, Rocco P, Miquel JF, Nervi F, Llop E, et al. (2000) Mitochondrial DNA polymorphisms in Chilean aboriginal populations: implications for the peopling of the southern cone of the continent. Am J Phys Anthropol 113: 19–29.ML MoragaP. RoccoJF MiquelF. NerviE. Llop2000Mitochondrial DNA polymorphisms in Chilean aboriginal populations: implications for the peopling of the southern cone of the continent.Am J Phys Anthropol1131929

46. Ginther C, Corach D, Penacino GA, Rey JA, Carnese FR, et al. (1993) Genetic variation among the Mapuche Indians from the Patagonian region of Argentina: mitochondrial DNA sequence variation and allele frequencies of several nuclear genes. Exs 67: 211–219.C. GintherD. CorachGA PenacinoJA ReyFR Carnese1993Genetic variation among the Mapuche Indians from the Patagonian region of Argentina: mitochondrial DNA sequence variation and allele frequencies of several nuclear genes.Exs67211219

47. Merriwether DA, Kemp BM, Crews DE, Neel JV (2000) Gene flow and genetic variation in the Yanomama as revealed by mitochondrial DNA. In: Renfrew C, editor. America Past, America Present: Genes and Languages in the Americas and Beyond. Cambridge: McDonald Institute for Archaeological Research. pp. 89–124.DA MerriwetherBM KempDE CrewsJV Neel2000Gene flow and genetic variation in the Yanomama as revealed by mitochondrial DNA.C. RenfrewAmerica Past, America Present: Genes and Languages in the Americas and BeyondCambridgeMcDonald Institute for Archaeological Research89124

48. Williams SR, Chagnon NA, Spielman RS (2002) Nuclear and mitochondrial genetic variation in the Yanomamo: a test case for ancient DNA studies of prehistoric populations. Am J Phys Anthropol 117: 246–259.SR WilliamsNA ChagnonRS Spielman2002Nuclear and mitochondrial genetic variation in the Yanomamo: a test case for ancient DNA studies of prehistoric populations.Am J Phys Anthropol117246259

49. Bonilla C, Bertoni B, Gonzalez S, Cardoso H, Brum-Zorrilla N, et al. (2004) Substantial Native American female contribution to the population of Tacuarembo, Uruguay, reveals past episodes of sex-biased gene flow. Am J Hum Biol 16: 289–297.C. BonillaB. BertoniS. GonzalezH. CardosoN. Brum-Zorrilla2004Substantial Native American female contribution to the population of Tacuarembo, Uruguay, reveals past episodes of sex-biased gene flow.Am J Hum Biol16289297

50. Pagano S, Sans M, Pimenoff V, Cantera AM, Alvarez JC, et al. (2005) Assessment of HV1 and HV2 mtDNA variation for forensic purposes in an Uruguayan population sample. J Forensic Sci 50: 1239–1242.S. PaganoM. SansV. PimenoffAM CanteraJC Alvarez2005Assessment of HV1 and HV2 mtDNA variation for forensic purposes in an Uruguayan population sample.J Forensic Sci5012391242

51. Vona G, Falchi A, Moral P, Calo CM, Varesi L (2005) Mitochondrial sequence variation in the Guahibo Amerindian population from Venezuela. Am J Phys Anthropol 127: 361–369.G. VonaA. FalchiP. MoralCM CaloL. Varesi2005Mitochondrial sequence variation in the Guahibo Amerindian population from Venezuela.Am J Phys Anthropol127361369

52. Torres MM, Bravi CM, Bortolini MC, Duque C, Callegari-Jacques S, et al. (2006) A revertant of the major founder Native American haplogroup C common in populations from northern South America. Am J Hum Biol 18: 59–65.MM TorresCM BraviMC BortoliniC. DuqueS. Callegari-Jacques2006A revertant of the major founder Native American haplogroup C common in populations from northern South America.Am J Hum Biol185965

53. Bert F, Corella A, Gene M, Perez-Perez A, Turbon D (2004) Mitochondrial DNA diversity in the Llanos de Moxos: Moxo, Movima and Yuracare Amerindian populations from Bolivia lowlands. Ann Hum Biol 31: 9–28.F. BertA. CorellaM. GeneA. Perez-PerezD. Turbon2004Mitochondrial DNA diversity in the Llanos de Moxos: Moxo, Movima and Yuracare Amerindian populations from Bolivia lowlands.Ann Hum Biol31928

54. Vernesi C, Fuselli S, Castri L, Bertorelle G, Barbujani G (2002) Mitochondrial diversity in linguistic isolates of the Alps: a reappraisal. Hum Biol 74: 725–730.C. VernesiS. FuselliL. CastriG. BertorelleG. Barbujani2002Mitochondrial diversity in linguistic isolates of the Alps: a reappraisal.Hum Biol74725730

55. Monsalve MV, Stone AC, Lewis CM, Rempel A, Richards M, et al. (2002) Brief communication: molecular analysis of the Kwaday Dan Ts'finchi ancient remains found in a glacier in Canada. Am J Phys Anthropol 119: 288–291.MV MonsalveAC StoneCM LewisA. RempelM. Richards2002Brief communication: molecular analysis of the Kwaday Dan Ts'finchi ancient remains found in a glacier in Canada.Am J Phys Anthropol119288291

56. Stone AC, Stoneking M (1993) Ancient DNA from a pre-Columbian Amerindian population. Am J Phys Anthropol 92: 463–471.AC StoneM. Stoneking1993Ancient DNA from a pre-Columbian Amerindian population.Am J Phys Anthropol92463471

57. Lalueza-Fox C, Calderon FL, Calafell F, Morera B, Bertranpetit J (2001) MtDNA from extinct Tainos and the peopling of the Caribbean. Ann Hum Genet 65: 137–151.C. Lalueza-FoxFL CalderonF. CalafellB. MoreraJ. Bertranpetit2001MtDNA from extinct Tainos and the peopling of the Caribbean.Ann Hum Genet65137151

58. Lalueza-Fox C, Gilbert MT, Martinez-Fuentes AJ, Calafell F, Bertranpetit J (2003) Mitochondrial DNA from pre-Columbian Ciboneys from Cuba and the prehistoric colonization of the Caribbean. Am J Phys Anthropol 121: 97–108.C. Lalueza-FoxMT GilbertAJ Martinez-FuentesF. CalafellJ. Bertranpetit2003Mitochondrial DNA from pre-Columbian Ciboneys from Cuba and the prehistoric colonization of the Caribbean.Am J Phys Anthropol12197108

59. Moraga M, Santoro CM, Standen VG, Carvallo P, Rothhammer F (2005) Microevolution in prehistoric Andean populations: chronologic mtDNA variation in the desert valleys of northern Chile. Am J Phys Anthropol 127: 170–181.M. MoragaCM SantoroVG StandenP. CarvalloF. Rothhammer2005Microevolution in prehistoric Andean populations: chronologic mtDNA variation in the desert valleys of northern Chile.Am J Phys Anthropol127170181

60. Ribeiro-Dos-Santos AK, Santos SE, Machado AL, Guapindaia V, Zago MA (1996) Heterogeneity of mitochondrial DNA haplotypes in Pre-Columbian Natives of the Amazon region. Am J Phys Anthropol 101: 29–37.AK Ribeiro-Dos-SantosSE SantosAL MachadoV. GuapindaiaMA Zago1996Heterogeneity of mitochondrial DNA haplotypes in Pre-Columbian Natives of the Amazon region.Am J Phys Anthropol1012937

61. García-Bour J, Pérez-Pérez A, Álvarez S, Fernández E, López-Parra AM, et al. (2004) Early population differentiation in extinct aborigines from Tierra del Fuego-Patagonia: ancient mtDNA sequences and Y-chromosome STR characterization. Am J Phys Anthropol 123: 361–370.J. García-BourA. Pérez-PérezS. ÁlvarezE. FernándezAM López-Parra2004Early population differentiation in extinct aborigines from Tierra del Fuego-Patagonia: ancient mtDNA sequences and Y-chromosome STR characterization.Am J Phys Anthropol123361370

62. Brown MD, Hosseini SH, Torroni A, Bandelt H-J, Allen JC, et al. (1998) mtDNA haplogroup X: an ancient link between Europe/Western Asia and North America? Am J Hum Genet 63: 1852–1861.MD BrownSH HosseiniA. TorroniH-J BandeltJC Allen1998mtDNA haplogroup X: an ancient link between Europe/Western Asia and North America?Am J Hum Genet6318521861

63. Dornelles CL, Bonatto SL, De Freitas LB, Salzano FM (2005) Is haplogroup X present in extant South American Indians? Am J Phys Anthropol 127: 439–448.CL DornellesSL BonattoLB De FreitasFM Salzano2005Is haplogroup X present in extant South American Indians?Am J Phys Anthropol127439448

64. Alves-Silva J, da Silva Santos M, Guimaraes PE, Ferreira AC, Bandelt H-J, et al. (2000) The ancestry of Brazilian mtDNA lineages. Am J Hum Genet 67: 444–461.J. Alves-SilvaM. da Silva SantosPE GuimaraesAC FerreiraH-J Bandelt2000The ancestry of Brazilian mtDNA lineages.Am J Hum Genet67444461

65. Ribeiro-Dos-Santos AK, Carvalho BM, Feio-Dos-Santos AC, Santos SE (2006) Nucleotide variability of HV-I in Afro-descendents populations of the Brazilian Amazon Region. Forensic Sci Int. AK Ribeiro-Dos-SantosBM CarvalhoAC Feio-Dos-SantosSE Santos2006Nucleotide variability of HV-I in Afro-descendents populations of the Brazilian Amazon Region.Forensic Sci Int

66. Smith DG, Malhi RS, Eshleman J, Lorenz JG, Kaestle FA (1999) Distribution of mtDNA haplogroup X among Native North Americans. Am J Phys Anthropol 110: 271–284.DG SmithRS MalhiJ. EshlemanJG LorenzFA Kaestle1999Distribution of mtDNA haplogroup X among Native North Americans.Am J Phys Anthropol110271284

67. Salas A, Richards M, Lareu MV, Scozzari R, Coppa A, et al. (2004) The African diaspora: mitochondrial DNA and the Atlantic slave trade. Am J Hum Genet 74: 454–465.A. SalasM. RichardsMV LareuR. ScozzariA. Coppa2004The African diaspora: mitochondrial DNA and the Atlantic slave trade.Am J Hum Genet74454465

68. Mendizabal I, Sandoval K, Berniell-Lee G, Calafell F, Salas A, et al. (2008) Genetic origin, admixture, and asymmetry in maternal and paternal human lineages in Cuba. BMC Evol Biol 8: 213.I. MendizabalK. SandovalG. Berniell-LeeF. CalafellA. Salas2008Genetic origin, admixture, and asymmetry in maternal and paternal human lineages in Cuba.BMC Evol Biol8213

69. Richards M, Macaulay V, Hickey E, Vega E, Sykes B, et al. (2000) Tracing European founder lineages in the Near Eastern mtDNA pool. Am J Hum Genet 67: 1251–1276.M. RichardsV. MacaulayE. HickeyE. VegaB. Sykes2000Tracing European founder lineages in the Near Eastern mtDNA pool.Am J Hum Genet6712511276

70. González AM, Brehm A, Pérez JA, Maca-Meyer N, Flores C, et al. (2003) Mitochondrial DNA affinities at the Atlantic fringe of Europe. Am J Phys Anthropol 120: 391–404.AM GonzálezA. BrehmJA PérezN. Maca-MeyerC. Flores2003Mitochondrial DNA affinities at the Atlantic fringe of Europe.Am J Phys Anthropol120391404

71. Brehm A, Pereira L, Kivisild T, Amorim A (2003) Mitochondrial portraits of the Madeira and Acores archipelagos witness different genetic pools of its settlers. Hum Genet 114: 77–86.A. BrehmL. PereiraT. KivisildA. Amorim2003Mitochondrial portraits of the Madeira and Acores archipelagos witness different genetic pools of its settlers.Hum Genet1147786

72. Pereira L, Prata MJ, Amorim A (2000) Diversity of mtDNA lineages in Portugal: not a genetic edge of European variation. Ann Hum Genet 64: 491–506.L. PereiraMJ PrataA. Amorim2000Diversity of mtDNA lineages in Portugal: not a genetic edge of European variation.Ann Hum Genet64491506

73. Malyarchuk BA, Grzybowski T, Derenko MV, Czarny J, Wozniak M, et al. (2002) Mitochondrial DNA variability in Poles and Russians. Ann Hum Genet 66: 261–283.BA MalyarchukT. GrzybowskiMV DerenkoJ. CzarnyM. Wozniak2002Mitochondrial DNA variability in Poles and Russians.Ann Hum Genet66261283

74. Salas A, Richards M, De la Fé T, Lareu MV, Sobrino B, et al. (2002) The making of the African mtDNA landscape. Am J Hum Genet 71: 1082–1111.A. SalasM. RichardsT. De la FéMV LareuB. Sobrino2002The making of the African mtDNA landscape.Am J Hum Genet7110821111

75. Salas A, Carracedo Á, Richards M, Macaulay V (2005) Charting the Ancestry of African Americans. Am J Hum Genet 77: 676–680.A. SalasÁ. CarracedoM. RichardsV. Macaulay2005Charting the Ancestry of African Americans.Am J Hum Genet77676680

76. Beleza S, Gusmão L, Amorim A, Carracedo Á, Salas A (2005) The genetic legacy of western Bantu migrations. Hum Genet 117: 366–375.S. BelezaL. GusmãoA. AmorimÁ. CarracedoA. Salas2005The genetic legacy of western Bantu migrations.Hum Genet117366375

77. Trovoada MJ, Pereira L, Gusmao L, Abade A, Amorim A, et al. (2004) Pattern of mtDNA variation in three populations from Sao Tome e Principe. Ann Hum Genet 68: 40–54.MJ TrovoadaL. PereiraL. GusmaoA. AbadeA. Amorim2004Pattern of mtDNA variation in three populations from Sao Tome e Principe.Ann Hum Genet684054

78. CONCULTURA (2000) Biblioteca de Historia Salvadoreña, El Salvador, Historia de sus pueblos, villas y ciudades. Consejo Nacional para la Cultura y el Arte. CONCULTURA2000Biblioteca de Historia Salvadoreña, El Salvador, Historia de sus pueblos, villas y ciudades.Consejo Nacional para la Cultura y el Arte

79. Salas A, Richards M, Lareu MV, Sobrino B, Silva S, et al. (2005) Shipwrecks and founder effects: Divergent demographic histories reflected in Caribbean mtDNA. Am J Phys Anthropol 128: 855–860.A. SalasM. RichardsMV LareuB. SobrinoS. Silva2005Shipwrecks and founder effects: Divergent demographic histories reflected in Caribbean mtDNA.Am J Phys Anthropol128855860

80. Lovo-Gómez J, Blanco-Verea A, Lareu MV, Brion M, Carracedo Á (2007) The genetic male legacy from El Salvador. Forensic Sci Int 171: 198–203.J. Lovo-GómezA. Blanco-VereaMV LareuM. BrionÁ. Carracedo2007The genetic male legacy from El Salvador.Forensic Sci Int171198203

81. Dipierri JE, Alfaro E, Martínez-Marignac VL, Bailliet G, Bravi CM, et al. (1998) Paternal directional mating in two Amerindian subpopulations located at different altitudes in northwestern Argentina. Hum Biol 70: 1001–1010.JE DipierriE. AlfaroVL Martínez-MarignacG. BaillietCM Bravi1998Paternal directional mating in two Amerindian subpopulations located at different altitudes in northwestern Argentina.Hum Biol7010011010

82. Salas A, Jaime JC, Álvarez-Iglesias V, Carracedo Á (2008) Gender bias in the multi-ethnic genetic composition of Central Argentina. J Hum Genet 53: 662–674.A. SalasJC JaimeV. Álvarez-IglesiasÁ. Carracedo2008Gender bias in the multi-ethnic genetic composition of Central Argentina.J Hum Genet53662674

83. Salas A, Bandelt H-J, Macaulay V, Richards MB (2007) Phylogeographic investigations: The role of trees in forensic genetics. Forensic Sci Int 168: 1–13.A. SalasH-J BandeltV. MacaulayMB Richards2007Phylogeographic investigations: The role of trees in forensic genetics.Forensic Sci Int168113

84. Burckhardt F, von Haeseler A, Meyer S (1999) HvrBase: compilation of mtDNA control region sequences from primates. Nucleic Acids Res 27: 138–142.F. BurckhardtA. von HaeselerS. Meyer1999HvrBase: compilation of mtDNA control region sequences from primates.Nucleic Acids Res27138142

85. Budowle B, Allard MW, Fisher CL, Isenberg AR, Monson KL, et al. (2002) HVI and HVII mitochondrial DNA data in Apaches and Navajos. Int J Legal Med 116: 212–215.B. BudowleMW AllardCL FisherAR IsenbergKL Monson2002HVI and HVII mitochondrial DNA data in Apaches and Navajos.Int J Legal Med116212215

86. Merriwether DA, Kemp BM, Crews DE, Neel JV (2002) Gene flow and genetic variation in the Yanomama as revealed by mitochondrial DNA.Cambridge: McDonald Institute for Archaeological Research. 89–124.DA MerriwetherBM KempDE CrewsJV Neel2002Gene flow and genetic variation in the Yanomama as revealed by mitochondrial DNA.CambridgeMcDonald Institute for Archaeological Research89124