Over the last 500 years, admixture among Amerindians, Europeans, and Africans, principally, has come to shape the present-day gene pool of Mexicans, particularly Mestizos, who represent about 93% of the total Mexican population. In this work, we analyze the genetic data of 13 combined DNA index system-short tandem repeats (CODIS-STRs) in 1,984 unrelated Mestizos representing 10 population samples from different regions of Mexico, namely North, West, Central, and Southeast. The analysis of molecular variance (AMOVA) test demonstrated low but significant differentiation among Mestizos from different regions (F(ST) = 0.34%; P = 0.0000). Although the spatial analysis of molecular variance (SAMOVA) predicted clustering Mestizo populations into four well-delimited groups, the main differentiation was observed between Northwest when compared with Central and Southeast regions. In addition, we included analysis of individuals of Amerindian (Purepechas), European (Huelva, Spain), and African (Fang) origin. Thus, STRUCTURE analysis was performed identifying three well-differentiated ancestral populations (k = 3). STRUCTURE results and admixture estimations by means of LEADMIX software in Mestizo populations demonstrated genetic heterogeneity or asymmetric admixture throughout Mexico, displaying an increasing North-to-South gradient of Amerindian ancestry, and vice versa regarding the European component. Interestingly, this distribution of Amerindian ancestry roughly reflects pre-Hispanic Native-population density, particularly toward the Mesoamerican area. The forensic, epidemiological, and evolutionary implications of these findings are discussed herein.

It has been suggested that the higher susceptibility of Hispanics to metabolic disease is related to their Native American heritage. A frequent cholesterol transporter ABCA1 (ATP-binding cassette transporter A1) gene variant (R230C, rs9282541) apparently exclusive to Native American individuals was associated with low high-density lipoprotein cholesterol (HDL-C) levels, obesity and type 2 diabetes in Mexican Mestizos. We performed a more extensive analysis of this variant in 4405 Native Americans and 863 individuals from other ethnic groups to investigate genetic evidence of positive selection, to assess its functional effect in vitro and to explore associations with HDL-C levels and other metabolic traits. The C230 allele was found in 29 of 36 Native American groups, but not in European, Asian or African individuals. C230 was observed on a single haplotype, and C230-bearing chromosomes showed longer relative haplotype extension compared with other haplotypes in the Americas. Additionally, single-nucleotide polymorphism data from the Human Genome Diversity Panel Native American populations were enriched in significant integrated haplotype score values in the region upstream of the ABCA1 gene. Cells expressing the C230 allele showed a 27% cholesterol efflux reduction (P< 0.001), confirming this variant has a functional effect in vitro. Moreover, the C230 allele was associated with lower HDL-C levels (P = 1.77 x 10(-11)) and with higher body mass index (P = 0.0001) in the combined analysis of Native American populations. This is the first report of a common functional variant exclusive to Native American and descent populations, which is a major determinant of HDL-C levels and may have contributed to the adaptive evolution of Native American populations.

A current issue on the settlement of the Americas refers to the lack of morphological affinities between early Holocene human remains (Palaeoamericans) and modern Amerindian groups, as well as the degree of contribution of the former to the gene pool of the latter. A different origin for Palaeoamericans and Amerindians is invoked to explain such a phenomenon. Under this hypothesis, the origin of Palaeoamericans must be traced back to a common ancestor for Palaeoamericans and Australians, which departed from somewhere in southern Asia and arrived in the Australian continent and the Americas around 40,000 and 12,000 years before present, respectively. Most modern Amerindians are believed to be part of a second, morphologically differentiated migration. Here we present evidence of a modern Amerindian group from the Baja California Peninsula in Mexico, showing clearer affinities with Palaeoamerican remains than with modern Amerindians. Climatic changes during the Middle Holocene probably generated the conditions for isolation from the continent, restricting the gene flow of the original group with northern populations, which resulted in the temporal continuity of the Palaeoamerican morphological pattern to the present.

Extensive European and African admixture coupled with loss of Amerindian lineages makes the reconstruction of pre-Columbian history of Native Americans based on present-day genomes extremely challenging. Still open questions remain about the dispersals that occurred throughout the continent after the initial peopling from the Beringia, especially concerning the number and dynamics of diffusions into South America. Indeed, if environmental and historical factors contributed to shape distinct gene pools in the Andes and Amazonia, the origins of this East-West genetic structure and the extension of further interactions between populations residing along this divide are still not well understood. To this end, we generated new high-resolution genome-wide data for 229 individuals representative of one Central and ten South Amerindian ethnic groups from Mexico, Peru, Bolivia, and Argentina. Low levels of European and African admixture in the sampled individuals allowed the application of fine-scale haplotype-based methods and demographic modeling approaches. These analyses revealed highly specific Native American genetic ancestries and great intragroup homogeneity, along with limited traces of gene flow mainly from the Andes into Peruvian Amazonians. Substantial amount of genetic drift differentially experienced by the considered populations underlined distinct patterns of recent inbreeding or prolonged isolation. Overall, our results support the hypothesis that all non-Andean South Americans are compatible with descending from a common lineage, while we found low support for common Mesoamerican ancestors of both Andeans and other South American groups. These findings suggest extensive back-migrations into Central America from non-Andean sources or conceal distinct peopling events into the Southern Continent.

Y-linked markers are suitable loci to analyze genetic diversity of human populations, offering knowledge of medical, forensic, and anthropological interest. In a population sample of 206 Mestizo males from western Mexico, we analyzed two binary loci (M3 and YAP) and six Y-STRs, adding to the analysis data of Mexican Mestizos and Amerindians, and relevant worldwide populations. The paternal ancestry estimated in western Mexican-Mestizos was mainly European (60-64%), followed by Amerindian (25-21%), and African ( approximately 15%). Significant genetic heterogeneity was established between Mestizos from western (Jalisco State) and northern Mexico (Chihuahua State) compared with Mexicans from the center of the Mexican Republic (Mexico City), this attributable to higher European ancestry in western and northern than in central and southeast populations, where higher Amerindian ancestry was inferred. This genetic structure has important implications for medical and forensic purposes. Two different Pre-Hispanic evolutionary processes were evident. In Mesoamerican region, populations presented higher migration rate (N(m) = 24.76), promoting genetic homogeneity. Conversely, isolated groups from the mountains and canyons of the Western and Northern Sierra Madre (Huichols and Tarahumaras, respectively) presented a lower migration rate (N(m) = 10.27) and stronger genetic differentiation processes (founder effect and/or genetic drift), constituting a Pre-Hispanic population substructure. Additionally, Tarahumaras presented a higher frequency of Y-chromosomes without Q3 that was explained by paternal European admixture (15%) and, more interestingly, by a distinctive Native-American ancestry. In Purepechas, a special admixture process involving preferential integration of non-Purepecha women in their communities could explain contrary genetic evidences (autosomal vs. Y-chromosome) for this tribe.

An evolutionary, diachronic approach to the phenotypic craniofacial pattern arisen in a human population after high levels of admixture and gene flow was achieved by means of geometric morphometrics. Admixture has long been studied after molecular data. Nevertheless, few efforts have been made to explain the morphological outcome in human craniofacial samples. The Spanish-Amerindian contact can be considered a good scenario for such an analysis. Here we present a comparative analysis of craniofacial shape changes observed between two putative ancestor groups, Spanish and precontact Aztecs, and two diachronic admixed groups, corresponding to early and late colonial periods from the Mexico's Central Valley. Quantitative shape comparisons of Amerindian, Spanish, and admixed groups were used to test the expectations of quantitative genetics for admixture events. In its simplest form, this prediction states that an admixed group will present phenotypic values falling between those of both parental groups. Results show that, in general terms, although the human skull is a complex, integrated structure, the craniofacial morphology observed fits the theoretical expectations of quantitative genetics. Thus, it is predictive of population structure and history. In fact, results obtained after the craniofacial analysis are in accordance with previous molecular and historical interpretations, providing evidence that admixture is a main microevolutionary agent influencing modern Mexican gene pool. However, expectations are not straightforward when moderate shape changes are considered. Deviations detected at localized structures, such as the upper and lower face, highlight the evolution of a craniofacial pattern exclusively inherent to the admixed groups, indicating that quantitative characters might respond to admixture in a complicated, nondirectional way.

The study of genetic information can reveal a reconstruction of human population’s history. We sequenced the entire mtDNA control region (positions 16.024 to 576 following Cambridge Reference Sequence, CRS) of 605 individuals from seven Mesoamerican indigenous groups and one Aridoamerican from the Greater Southwest previously defined, all of them in present Mexico. Samples were collected directly from the indigenous populations, the application of an individual survey made it possible to remove related or with other origins samples. Diversity indices and demographic estimates were calculated. Also AMOVAs were calculated according to different criteria. An MDS plot, based on FST distances, was also built. We carried out the construction of individual networks for the four Amerindian haplogroups detected. Finally, barrier software was applied to detect genetic boundaries among populations. The results suggest: a common origin of the indigenous groups; a small degree of European admixture; and inter-ethnic gene flow. The process of Mesoamerica’s human settlement took place quickly influenced by the region’s orography, which development of genetic and cultural differences facilitated. We find the existence of genetic structure is related to the region’s geography, rather than to cultural parameters, such as language. The human population gradually became fragmented, though they remained relatively isolated, and differentiated due to small population sizes and different survival strategies. Genetic differences were detected between Aridoamerica and Mesoamerica, which can be subdivided into “East”, “Center”, “West” and “Southeast”. The fragmentation process occurred mainly during the Mesoamerican Pre-Classic period, with the Otomí being one of the oldest groups. With an increased number of populations studied adding previously published data, there is no change in the conclusions, although significant genetic heterogeneity can be detected in Pima and Huichol groups. This result may be explained because populations historically assigned as belonging to the same group were, in fact, different indigenous populations.

No abstract