Biogeographic Perspectives on Human Genetic Diversification

Modern humans originated in Africa 300,000 yr ago, and before leaving their continent of origin, they underwent a process of intense diversification involving complex demographic dynamics. Upon exiting Africa, different populations emerged on the four other inhabited continents, shaped by the interplay of various evolutionary processes, such as migrations, founder effects, and natural selection. Within each region, continental populations, in turn, diversified and evolved almost independently for millennia. As a backdrop to this diversification, introgressions from archaic species contributed to establishing different patterns of genetic diversity in different geographic regions, reshaping our understanding of our species' variability. With the increasing availability of genomic data, it has become possible to delineate the subcontinental human population structure precisely. However, the bias toward the genomic research focused on populations from the global North has limited our understanding of the real diversity of our species and the processes and events that guided different human groups throughout their evolutionary history. This perspective is part of a series of articles celebrating 40 yr since our journal, Molecular Biology and Evolution, was founded (Russo et al. 2024). The perspective is accompanied by virtual issues, a selection of papers on human diversification published by Genome Biology and Evolution and Molecular Biology and Evolution.

Human genetic determinants of COVID-19 in Brazil: challenges and future plans

COVID-19 pandemic represented a worldwide major challenge in different areas, and efforts undertaken by the scientific community led to the understanding of some of the genetic determinants that influence the different COVID-19 outcomes. In this paper, we review the studies about the role of human genetics in COVID-19 severity and how Brazilian studies also contributed to those findings. Rare variants in genes related to Inborn Errors of Immunity (IEI) in the type I interferons pathway, and its phenocopies, have been described as being causative of severe outcomes. IEI and its phenocopies are present in Brazil, not only in COVID-19 patients, but also in autoimmune conditions and severe reactions to yellow fever vaccine. In addition, studies focusing on common variants and GWAS studies encompassing worldwide patients have found several loci associated with COVID-19 severity. A GWAS study including only Brazilian COVID-19 patients identified a new locus 1q32.1 associated with COVID-19 severity. Thus, more comprehensive studies considering the Brazilian genomic diversity should be performed, since they can help to reveal not only what are the genetic determinants that contribute to the different outcomes for COVID-19 in the Brazilian population, but in the understanding of human genetics in different health conditions.

A multidisciplinary overview on the Tupi-speaking people expansion

The cultural and biological diversity of South American indigenous groups represent extremes of human variability, exhibiting one of the highest linguistic diversities alongside a remarkably low within-population genetic variation and an extremely high inter-population genetic differentiation. On top of that, this region has seen some of the most dramatic demographic events in human history unleashed by the European colonization of the Americas. As a result of this process, the distribution of indigenous populations has been radically changed. In this review we focus on the Tupi, the largest and most widespread linguistic family in eastern South America. Tupi are believed to have originated in southwestern Amazon, from where some of its subfamilies expanded into other parts of the Amazon and, in the case of the Tupi-Guarani, beyond its borders. Recent evidence from archaeology, linguistics, and genetics aligns with José Brochado's Tupi Expansion model. He proposed that the gradual development of agricultural systems within the Amazon resulted in population growth and, eventually, territorial expansion. This model also supports separate Tupi Expansion branches: Tupinambá (Atlantic coast) and Guarani (south, midwest Brazil). Although being the most populous group on Brazil's Atlantic coast, which was the most affected by European colonization, the Tupi still account for roughly 20% of the country's overall indigenous population. Finally, despite its importance and more than a century of research on the Tupi and their expansion history, many key questions remain unanswered, which we attempt to summarize and explore here.

Neanderthal introgression in SCN9A impacts mechanical pain sensitivity

The Nav1.7 voltage-gated sodium channel plays a key role in nociception. Three functional variants in the SCN9A gene (encoding M932L, V991L, and D1908G in Nav1.7), have recently been identified as stemming from Neanderthal introgression and to associate with pain symptomatology in UK BioBank data. In 1000 genomes data, these variants are absent in Europeans but common in Latin Americans. Analysing high-density genotype data from 7594 Latin Americans, we characterized Neanderthal introgression in SCN9A. We find that tracts of introgression occur on a Native American genomic background, have an average length of ~123 kb and overlap the M932L, V991L, and D1908G coding positions. Furthermore, we measured experimentally six pain thresholds in 1623 healthy Colombians. We found that Neanderthal ancestry in SCN9A is significantly associated with a lower mechanical pain threshold after sensitization with mustard oil and evidence of additivity of effects across Nav1.7 variants. Our findings support the reported association of Neanderthal Nav1.7 variants with clinical pain, define a specific sensory modality affected by archaic introgression in SCN9A and are consistent with independent effects of the Neanderthal variants on Nav1.7 function.

Genomic history of coastal societies from eastern South America

Sambaqui (shellmound) societies are among the most intriguing archaeological phenomena in pre-colonial South America, extending from approximately 8,000 to 1,000 years before present (yr BP) across 3,000 km on the Atlantic coast. However, little is known about their connection to early Holocene hunter-gatherers, how this may have contributed to different historical pathways and the processes through which late Holocene ceramists came to rule the coast shortly before European contact. To contribute to our understanding of the population history of indigenous societies on the eastern coast of South America, we produced genome-wide data from 34 ancient individuals as early as 10,000 yr BP from four different regions in Brazil. Early Holocene hunter-gatherers were found to lack shared genetic drift among themselves and with later populations from eastern South America, suggesting that they derived from a common radiation and did not contribute substantially to later coastal groups. Our analyses show genetic heterogeneity among contemporaneous Sambaqui groups from the southeastern and southern Brazilian coast, contrary to the similarity expressed in the archaeological record. The complex history of intercultural contact between inland horticulturists and coastal populations becomes genetically evident during the final horizon of Sambaqui societies, from around 2,200 yr BP, corroborating evidence of cultural change.

Automatic landmarking identifies new loci associated with face morphology and implicates Neanderthal introgression in human nasal shape

We report a genome-wide association study of facial features in >6000 Latin Americans based on automatic landmarking of 2D portraits and testing for association with inter-landmark distances. We detected significant associations (P-value <5 × 10^-8) at 42 genome regions, nine of which have been previously reported. In follow-up analyses, 26 of the 33 novel regions replicate in East Asians, Europeans, or Africans, and one mouse homologous region influences craniofacial morphology in mice. The novel region in 1q32.3 shows introgression from Neanderthals and we find that the introgressed tract increases nasal height (consistent with the differentiation between Neanderthals and modern humans). Novel regions include candidate genes and genome regulatory elements previously implicated in craniofacial development, and show preferential transcription in cranial neural crest cells. The automated approach used here should simplify the collection of large study samples from across the world, facilitating a cosmopolitan characterization of the genetics of facial features.

Population-based change-point detection for the identification of homozygosity islands

MOTIVATION

This work is motivated by the problem of identifying homozygosity islands on the genome of individuals in a population. Our method directly tackles the issue of identification of the homozygosity islands at the population level, without the need of analysing single individuals and then combine the results, as is made nowadays in state-of-the-art approaches.

RESULTS

We propose regularised offline change-point methods to detect changes in the parameters of a multidimensional distribution when we have several aligned, independent samples of fixed resolution. We present a penalised maximum likelihood approach that can be efficiently computed by a dynamic programming algorithm or approximated by a fast binary segmentation algorithm. Both estimators are shown to converge almost surely to the set of change-points without the need of specifying a priori the number of change-points. In simulation we observed similar performances from the exact and greedy estimators. Moreover, we provide a new methodology for the selection of the regularisation constant which has the advantage of being automatic, consistent and less prone to subjective analysis.

AVAILABILITY

The data used in the application is from the Human Genome Diversity Project (HGDP) and is publicly available. Algorithms were implemented using the R software R Core Team (2020) in the R package blockcpd, found at https://github.com/Lucas-Prates/blockcpd.

SUPPLEMENTARY INFORMATION

Supplementary material is available online at Bioinformatics.

Indigenous people from Amazon show genetic signatures of pathogen-driven selection

Ecological conditions in the Amazon rainforests are historically favorable for the transmission of numerous tropical diseases, especially vector-borne diseases. The high diversity of pathogens likely contributes to the strong selective pressures for human survival and reproduction in this region. However, the genetic basis of human adaptation to this complex ecosystem remains unclear. This study investigates the possible footprints of genetic adaptation to the Amazon rainforest environment by analyzing the genomic data of 19 native populations. The results based on genomic and functional analysis showed an intense signal of natural selection in a set of genes related to Trypanosoma cruzi infection, which is the pathogen responsible for Chagas disease, a neglected tropical parasitic disease native to the Americas that is currently spreading worldwide.

Precision medicine implementation challenges for APOL1 testing in chronic kidney disease in admixed populations

Chronic Kidney Disease (CKD) is a public health problem that presents genetic and environmental risk factors. Two alleles in the Apolipoprotein L1 (APOL1) gene were associated with chronic kidney disease; these alleles are common in individuals of African ancestry but rare in European descendants. Genomic studies on Afro-Americans have indicated a higher prevalence and severity of chronic kidney disease in people of African ancestry when compared to other ethnic groups. However, estimates in low- and middle-income countries are still limited. Precision medicine approaches could improve clinical outcomes in carriers of risk alleles in the Apolipoprotein L1 gene through early diagnosis and specific therapies. Nevertheless, to enhance the definition of studies on these variants, it would be necessary to include individuals with different ancestry profiles in the sample, such as Latinos, African Americans, and Indigenous peoples. There is evidence that measuring genetic ancestry improves clinical care for admixed people. For chronic kidney disease, this knowledge could help establish public health strategies for monitoring patients and understanding the impact of the Apolipoprotein L1 genetic variants in admixed populations. Therefore, researchers need to develop resources, methodologies, and incentives for vulnerable and disadvantaged communities, to develop and implement precision medicine strategies and contribute to consolidating diversity in science and precision medicine in clinical practice.

A genomic perspective on South American human history

It has generally been accepted that the current indigenous peoples of the Americas are derived from ancestors from northeastern Asia. The latter were believed to have spread into the American continent by the end of the Last Glacial Maximum. In this sense, a joint and in-depth study of the earliest settlement of East Asia and the Americas is required to elucidate these events accurately. The first Americans underwent an adaptation process to the Americas' vast environmental diversity, mediated by biological and cultural evolution and niche construction, resulting in enormous cultural diversity, a wealth of domesticated species, and extensive landscape modifications. Afterward, in the Late Holocene, the advent of intensive agricultural food production systems, sedentism, and climate change significantly reshaped genetic and cultural diversity across the continent, particularly in the Andes and Amazonia. Furthermore, starting around the end of the 15th century, European colonization resulted in massive extermination of indigenous peoples and extensive admixture. Thus, the present review aims to create a comprehensive picture of the main events involved in the formation of contemporary South American indigenous populations and the dynamics responsible for shaping their genetic diversity by integrating current genetic data with evidence from archeology, linguistics and other disciplines.

Disentangling Signatures of Selection Before and After European Colonization in Latin Americans

Throughout human evolutionary history, large-scale migrations have led to intermixing (i.e., admixture) between previously separated human groups. Although classical and recent work have shown that studying admixture can yield novel historical insights, the extent to which this process contributed to adaptation remains underexplored. Here, we introduce a novel statistical model, specific to admixed populations, that identifies loci under selection while determining whether the selection likely occurred post-admixture or prior to admixture in one of the ancestral source populations. Through extensive simulations, we show that this method is able to detect selection, even in recently formed admixed populations, and to accurately differentiate between selection occurring in the ancestral or admixed population. We apply this method to genome-wide SNP data of ∼4,000 individuals in five admixed Latin American cohorts from Brazil, Chile, Colombia, Mexico, and Peru. Our approach replicates previous reports of selection in the human leukocyte antigen region that are consistent with selection post-admixture. We also report novel signals of selection in genomic regions spanning 47 genes, reinforcing many of these signals with an alternative, commonly used local-ancestry-inference approach. These signals include several genes involved in immunity, which may reflect responses to endemic pathogens of the Americas and to the challenge of infectious disease brought by European contact. In addition, some of the strongest signals inferred to be under selection in the Native American ancestral groups of modern Latin Americans overlap with genes implicated in energy metabolism phenotypes, plausibly reflecting adaptations to novel dietary sources available in the Americas.

Population Histories and Genomic Diversity of South American Natives

South America is home to one of the most culturally diverse present-day native populations. However, the dispersion pattern, genetic substructure, and demographic complexity within South America are still poorly understood. Based on genome-wide data of 58 native populations, we provide a comprehensive scenario of South American indigenous groups considering the genomic, environmental, and linguistic data. Clear patterns of genetic structure were inferred among the South American natives, presenting at least four primary genetic clusters in the Amazonian and savanna regions and three clusters in the Andes and Pacific coast. We detected a cline of genetic variation along a west-east axis, contradicting a hard Andes-Amazon divide. This longitudinal genetic variation seemed to have been shaped by both serial population bottlenecks and isolation by distance. Results indicated that present-day South American substructures recapitulate ancient macroregional ancestries and western Amazonia groups show genetic evidence of cultural exchanges that led to language replacement in precontact times. Finally, demographic inferences pointed to a higher resilience of the western South American groups regarding population collapses caused by the European invasion and indicated precontact population reductions and demic expansions in South America.

A genome-wide association study identifies novel gene associations with facial skin wrinkling and mole count in Latin Americans

BACKGROUND

Genome-wide association studies (GWASs) have identified genes influencing skin ageing and mole count in Europeans, but little is known about the relevance of these (or other genes) in non-Europeans.

OBJECTIVES

To conduct a GWAS for facial skin ageing and mole count in adults < 40 years old, of mixed European, Native American and African ancestry, recruited in Latin America.

METHODS

Skin ageing and mole count scores were obtained from facial photographs of over 6000 individuals. After quality control checks, three wrinkling traits and mole count were retained for genetic analyses. DNA samples were genotyped with Illumina's HumanOmniExpress chip. Association testing was performed on around 8 703 729 single-nucleotide polymorphisms (SNPs) across the autosomal genome.

RESULTS

Genome-wide significant association was observed at four genome regions: two were associated with wrinkling (in 1p13·3 and 21q21·2), one with mole count (in 1q32·3) and one with both wrinkling and mole count (in 5p13·2). Associated SNPs in 5p13·2 and in 1p13·3 are intronic within SLC45A2 and VAV3, respectively, while SNPs in 1q32·3 are near the SLC30A1 gene, and those in 21q21·2 occur in a gene desert. Analyses of SNPs in IRF4 and MC1R are consistent with a role of these genes in skin ageing.

CONCLUSIONS

We replicate the association of wrinkling with variants in SLC45A2, IRF4 and MC1R reported in Europeans. We identify VAV3 and SLC30A1 as two novel candidate genes impacting on wrinkling and mole count, respectively. We provide the first evidence that SLC45A2 influences mole count, in addition to variants in this gene affecting melanoma risk in Europeans.

Prediction of eye, hair and skin colour in Latin Americans

Here we evaluate the accuracy of prediction for eye, hair and skin pigmentation in a dataset of > 6500 individuals from Mexico, Colombia, Peru, Chile and Brazil (including genome-wide SNP data and quantitative/categorical pigmentation phenotypes - the CANDELA dataset CAN). We evaluated accuracy in relation to different analytical methods and various phenotypic predictors. As expected from statistical principles, we observe that quantitative traits are more sensitive to changes in the prediction models than categorical traits. We find that Random Forest or Linear Regression are generally the best performing methods. We also compare the prediction accuracy of SNP sets defined in the CAN dataset (including 56, 101 and 120 SNPs for eye, hair and skin colour prediction, respectively) to the well-established HIrisPlex-S SNP set (including 6, 22 and 36 SNPs for eye, hair and skin colour prediction respectively). When training prediction models on the CAN data, we observe remarkably similar performances for HIrisPlex-S and the larger CAN SNP sets for the prediction of hair (categorical) and eye (both categorical and quantitative), while the CAN sets outperform HIrisPlex-S for quantitative, but not for categorical skin pigmentation prediction. The performance of HIrisPlex-S, when models are trained in a world-wide sample (although consisting of 80% Europeans, https://hirisplex.erasmusmc.nl), is lower relative to training in the CAN data (particularly for hair and skin colour). Altogether, our observations are consistent with common variation of eye and hair colour having a relatively simple genetic architecture, which is well captured by HIrisPlex-S, even in admixed Latin Americans (with partial European ancestry). By contrast, since skin pigmentation is a more polygenic trait, accuracy is more sensitive to prediction SNP set size, although here this effect was only apparent for a quantitative measure of skin pigmentation. Our results support the use of HIrisPlex-S in the prediction of categorical pigmentation traits for forensic purposes in Latin America, while illustrating the impact of training datasets on its accuracy.

A GWAS in Latin Americans identifies novel face shape loci, implicating VPS13B and a Denisovan introgressed region in facial variation

To characterize the genetic basis of facial features in Latin Americans, we performed a genome-wide association study (GWAS) of more than 6000 individuals using 59 landmark-based measurements from two-dimensional profile photographs and ~9,000,000 genotyped or imputed single-nucleotide polymorphisms. We detected significant association of 32 traits with at least 1 (and up to 6) of 32 different genomic regions, more than doubling the number of robustly associated face morphology loci reported until now (from 11 to 23). These GWAS hits are strongly enriched in regulatory sequences active specifically during craniofacial development. The associated region in 1p12 includes a tract of archaic adaptive introgression, with a Denisovan haplotype common in Native Americans affecting particularly lip thickness. Among the nine previously unidentified face morphology loci we identified is the VPS13B gene region, and we show that variants in this region also affect midfacial morphology in mice.

Adaptation and co-adaptation of skin pigmentation and vitamin D genes in native Americans

We carried out an exhaustive review regarding human skin color variation and how much it may be related to vitamin D metabolism and other photosensitive molecules. We discuss evolutionary contexts that modulate this variability and hypotheses postulated to explain them; for example, a small amount of melanin in the skin facilitates vitamin D production, making it advantageous to have fair skin in an environment with little radiation incidence. In contrast, more melanin protects folate from degradation in an environment with a high incidence of radiation. Some Native American populations have a skin color at odds with what would be expected for the amount of radiation in the environment in which they live, a finding challenging the so-called "vitamin D-folate hypothesis." Since food is also a source of vitamin D, dietary habits should also be considered. Here we argue that a gene network approach provides tools to explain this phenomenon since it indicates potential alleles co-evolving in a compensatory way. We identified alleles of the vitamin D metabolism and pigmentation pathways segregated together, but in different proportions, in agriculturalists and hunter-gatherers. Finally, we highlight how an evolutionary approach can be useful to understand current topics of medical interest.

Demographic history and selection at HLA loci in Native Americans

The American continent was the last to be occupied by modern humans, and native populations bear the marks of recent expansions, bottlenecks, natural selection, and population substructure. Here we investigate how this demographic history has shaped genetic variation at the strongly selected HLA loci. In order to disentangle the relative contributions of selection and demography process, we assembled a dataset with genome-wide microsatellites and HLA-A, -B, -C, and -DRB1 typing data for a set of 424 Native American individuals. We find that demographic history explains a sizeable fraction of HLA variation, both within and among populations. A striking feature of HLA variation in the Americas is the existence of alleles which are present in the continent but either absent or very rare elsewhere in the world. We show that this feature is consistent with demographic history (i.e., the combination of changes in population size associated with bottlenecks and subsequent population expansions). However, signatures of selection at HLA loci are still visible, with significant evidence selection at deeper timescales for most loci and populations, as well as population differentiation at HLA loci exceeding that seen at neutral markers.

Roadmap for Establishing Large-Scale Genomic Medicine Initiatives in Low- and Middle-Income Countries

In the post-genomic era, genomic medicine interventions as a key component of personalized medicine and tailored-made health care are greatly anticipated following recent scientific and technological advances. Indeed, large-scale sequencing efforts that explore human genomic variation have been initiated in several, mostly developed, countries across the globe, such as the United States, the United Kingdom, and a few others. Here, we highlight the successful implementation of large-scale national genomic initiatives, namely the Genome of Greece (GoGreece) and the DNA do Brasil (DNABr), aiming to emphasize the importance of implementing such initiatives in developing countries. Based on this experience, we also provide a roadmap for replicating these projects in other low-resource settings, thereby bringing genomic medicine in these countries closer to clinical fruition.

Insights on hair, skin and eye color of ancient and contemporary Native Americans

Over the past few years, tools capable of predicting pigmentation phenotypes have been developed aiming to contribute for criminal and anthropological investigations. In this study, we used eight genetic systems to infer eye, hair, and skin color of ancient and contemporary Native Americans. To achieve this goal, we retrieved 61 SNPs from 42 samples available in free online repositories of DNA sequences. We performed pigmentation predictions using two freely available tools, HIrisPlex-S and Snipper, in addition to two other published models. This workflow made possible to predict all three phenotypes with at least one tool for 29 out of the 42 samples. Considering these 29 individuals, predictions for eye, hair, and skin color were obtained with HIrisPlex-S for 27, 28 and 27 individuals, respectively, while 24, 25 and 25 individuals had such predictions with Snipper. In general, ancient and contemporary Native Americans were predicted to have intermediate/brown eyes, black hair, and intermediate/darker skin pigmentation.

Reconstructed Lost Native American Populations from Eastern Brazil Are Shaped by Differential Jê/Tupi Ancestry

After the colonization of the Americas by Europeans and the consequent Trans-Atlantic Slave Trade, most Native American populations in eastern Brazil disappeared or went through an admixture process that configured a population composed of three main genetic components: the European, the sub-Saharan African, and the Native American. The study of the Native American genetic history is challenged by the lack of availability of genome-wide samples from Native American populations, the technical difficulties to develop ancient DNA studies, and the low proportions of the Native American component in the admixed Brazilian populations (on average 7%). We analyzed genome-wide data of 5,825 individuals from three locations of eastern Brazil: Salvador (North-East), Bambui (South-East), and Pelotas (South) and we reconstructed populations that emulate the Native American groups that were living in the 16th century around the sampling locations. This genetic reconstruction was performed after local ancestry analysis of the admixed Brazilian populations, through the rearrangement of the Native American haplotypes into reconstructed individuals with full Native American ancestry (51 reconstructed individuals in Salvador, 45 in Bambui, and 197 in Pelotas). We compared the reconstructed populations with nonadmixed Native American populations from other regions of Brazil through haplotype-based methods. Our results reveal a population structure shaped by the dichotomy of Tupi-/Jê-speaking ancestry related groups. We also show evidence of a decrease of the diversity of nonadmixed Native American groups after the European contact, in contrast with the reconstructed populations, suggesting a reservoir of the Native American genetic diversity within the admixed Brazilian population.

Latin Americans show wide-spread Converso ancestry and imprint of local Native ancestry on physical appearance

Historical records and genetic analyses indicate that Latin Americans trace their ancestry mainly to the intermixing (admixture) of Native Americans, Europeans and Sub-Saharan Africans. Using novel haplotype-based methods, here we infer sub-continental ancestry in over 6,500 Latin Americans and evaluate the impact of regional ancestry variation on physical appearance. We find that Native American ancestry components in Latin Americans correspond geographically to the present-day genetic structure of Native groups, and that sources of non-Native ancestry, and admixture timings, match documented migratory flows. We also detect South/East Mediterranean ancestry across Latin America, probably stemming mostly from the clandestine colonial migration of Christian converts of non-European origin (Conversos). Furthermore, we find that ancestry related to highland (Central Andean) versus lowland (Mapuche) Natives is associated with variation in facial features, particularly nose morphology, and detect significant differences in allele frequencies between these groups at loci previously associated with nose morphology in this sample.

Reconstructing the Deep Population History of Central and South America

We report genome-wide ancient DNA from 49 individuals forming four parallel time transects in Belize, Brazil, the Central Andes, and the Southern Cone, each dating to at least ∼9,000 years ago. The common ancestral population radiated rapidly from just one of the two early branches that contributed to Native Americans today. We document two previously unappreciated streams of gene flow between North and South America. One affected the Central Andes by ∼4,200 years ago, while the other explains an affinity between the oldest North American genome associated with the Clovis culture and the oldest Central and South Americans from Chile, Brazil, and Belize. However, this was not the primary source for later South Americans, as the other ancient individuals derive from lineages without specific affinity to the Clovis-associated genome, suggesting a population replacement that began at least 9,000 years ago and was followed by substantial population continuity in multiple regions.

Selection scan reveals three new loci related to high altitude adaptation in Native Andeans

The Andean Altiplano has been occupied continuously since the late Pleistocene, ~12,000 years ago, which places the Andean natives as one of the most ancient populations living at high altitudes. In the present study, we analyzed genomic data from Native Americans living a long-time at Andean high altitude and at Amazonia and Mesoamerica lowland areas. We have identified three new candidate genes - SP100, DUOX2 and CLC - with evidence of positive selection for altitude adaptation in Andeans. These genes are involved in the TP53 pathway and are related to physiological routes important for high-altitude hypoxia response, such as those linked to increased angiogenesis, skeletal muscle adaptations, and immune functions at the fetus-maternal interface. Our results, combined with other studies, showed that Andeans have adapted to the Altiplano in different ways and using distinct molecular strategies as compared to those of other natives living at high altitudes.

Contrasting patterns of RUNX2 repeat variations are associated with palate shape in phyllostomid bats and New World primates

Establishing the genetic basis that underlies craniofacial variability in natural populations is one of the main topics of evolutionary and developmental studies. One of the genes associated with mammal craniofacial variability is RUNX2, and in the present study we investigated the association between craniofacial length and width and RUNX2 across New World bats (Phyllostomidae) and primates (Catarrhini and Platyrrhini). Our results showed contrasting patterns of association between the glutamate/alanine ratios (Q/A ratio) and palate shape in these highly diverse groups. In phyllostomid bats, we found an association between shorter/broader faces and increase of the Q/A ratio. In New World monkeys (NWM) there was a positive correlation of increasing Q/A ratios to more elongated faces. Our findings reinforced the role of the Q/A ratio as a flexible genetic mechanism that would rapidly change the time of skull ossification throughout development. However, we propose a scenario in which the influence of this genetic adjustment system is indirect. The Q/A ratio would not lead to a specific phenotype, but throughout the history of a lineage, would act along with evolutionary constraints, as well as other genes, as a facilitator for adaptive morphological changes.

Developmental pathways inferred from modularity, morphological integration and fluctuating asymmetry patterns in the human face

Facial asymmetries are usually measured and interpreted as proxies to developmental noise. However, analyses focused on its developmental and genetic architecture are scarce. To advance on this topic, studies based on a comprehensive and simultaneous analysis of modularity, morphological integration and facial asymmetries including both phenotypic and genomic information are needed. Here we explore several modularity hypotheses on a sample of Latin American mestizos, in order to test if modularity and integration patterns differ across several genomic ancestry backgrounds. To do so, 4104 individuals were analyzed using 3D photogrammetry reconstructions and a set of 34 facial landmarks placed on each individual. We found a pattern of modularity and integration that is conserved across sub-samples differing in their genomic ancestry background. Specifically, a signal of modularity based on functional demands and organization of the face is regularly observed across the whole sample. Our results shed more light on previous evidence obtained from Genome Wide Association Studies performed on the same samples, indicating the action of different genomic regions contributing to the expression of the nose and mouth facial phenotypes. Our results also indicate that large samples including phenotypic and genomic metadata enable a better understanding of the developmental and genetic architecture of craniofacial phenotypes.

Relación entre tratamiento hormonal, cirugía-ortodoncia maxilofacial, traumatismos y malformaciones craneofaciales y la asimetría fluctuante

En este trabajo se evalúa la relación entre la asimetría fluctuante facial (AFF) y los tratamientos hormonales, cirugías maxilofaciales, ortodoncia, traumatismos y malformaciones. En el marco del proyecto CANDELA, se tomaron cinco fotografías faciales de 3162 voluntarios entre los 18 y 85 años. Por fotogrametría se colocaron 34 landmarks o puntos en 3D y mediante el método Procrustes ANOVA se obtuvieron valores individuales de asimetría fluctuante facial. Se realizó una prueba de ANOVA de una vía y la prueba de Welch y Levene para conocer las diferencias entre media y varianza de los valores de asimetría facial y las variables respuesta. También, se caracterizó la variación morfológica del componente asimétrico de la forma facial mediante técnicas multivariadas sobre los grupos que resultaran diferentes significativamente. Las mujeres que reportaron haber recibido algún tipo de tratamiento hormonal mostraron mayores valores de asimetría fluctuante facial respecto al grupo sin tratamiento. Esta asociación se mantuvo una vez removido el efecto de la ancestría genética y sin interactuar con el resto de variables incluidas en el análisis. Los cambios morfológicos asociados a este factor se concentran en el mentón, maxilar labio inferior, región perifrontal, región nasal y orejas. Algunos trabajos anteriores dieron cuenta de la posible relación entre la asimetría facial y los niveles de hormonas, pero no hay estudios que sustenten la relación causal o directa entre la asociación aquí planteada. El presente trabajo es una evidencia más de la asociación entre el consumo de hormonas y modificaciones de caracteres faciales en para poblaciones urbanas mestizas latinoamericanas.

MRPL53, a New Candidate Gene for Orofacial Clefting, Identified Using an eQTL Approach

A valuable approach to understand how individual and population genetic differences can predispose to disease is to assess the impact of genetic variants on cellular functions (e.g., gene expression) of cell and tissue types related to pathological states. To understand the genetic basis of nonsyndromic cleft lip with or without cleft palate (NSCL/P) susceptibility, a complex and highly prevalent congenital malformation, we searched for genetic variants with a regulatory role in a disease-related tissue, the lip muscle (orbicularis oris muscle [OOM]), of affected individuals. From 46 OOM samples, which are frequently discarded during routine corrective surgeries on patients with orofacial clefts, we derived mesenchymal stem cells and correlated the individual genetic variants with gene expression from these cultured cells. Through this strategy, we detected significant cis-eQTLs (i.e., DNA variants affecting gene expression) and selected a few candidates to conduct an association study in a large Brazilian cohort (624 patients and 668 controls). This resulted in the discovery of a novel susceptibility locus for NSCL/P, rs1063588, the best eQTL for the MRPL53 gene, where evidence for association was mostly driven by the Native American ancestry component of our Brazilian sample. MRPL53 (2p13.1) encodes a 39S protein subunit of mitochondrial ribosomes and interacts with MYC, a transcription factor required for normal facial morphogenesis. Our study illustrates not only the importance of sampling admixed populations but also the relevance of measuring the functional effects of genetic variants over gene expression to dissect the complexity of disease phenotypes.

Haplotype Study in SCA10 Families Provides Further Evidence for a Common Ancestral Origin of the Mutation

Spinocerebellar ataxia type 10 (SCA10) is an autosomal dominant neurodegenerative disorder characterized by progressive cerebellar ataxia and epilepsy. The disease is caused by a pentanucleotide ATTCT expansion in intron 9 of the ATXN10 gene on chromosome 22q13.3. SCA10 has shown a geographical distribution throughout America with a likely degree of Amerindian ancestry from different countries so far. Currently available data suggest that SCA10 mutation might have spread out early during the peopling of the Americas. However, the ancestral origin of SCA10 mutation remains under speculation. Samples of SCA10 patients from two Latin American countries were analysed, being 16 families from Brazil (29 patients) and 21 families from Peru (27 patients) as well as 49 healthy individuals from Indigenous Quechua population and 51 healthy Brazilian individuals. Four polymorphic markers spanning a region of 5.2 cM harbouring the ATTCT expansion were used to define the haplotypes, which were genotyped by different approaches. Our data have shown that 19-CGGC-14 shared haplotype was found in 47% of Brazilian and in 63% of Peruvian families. Frequencies from both groups are not statistically different from Quechua controls (57%), but they are statistically different from Brazilian controls (12%) (p < 0.001). The most frequent expanded haplotype in Quechuas, 19-15-CGGC-14-10, is found in 50% of Brazilian and in 65% of Peruvian patients with SCA10. These findings bring valuable evidence that ATTCT expansion may have arisen in a Native American chromosome.

A tale of agriculturalists and hunter-gatherers: Exploring the thrifty genotype hypothesis in native South Americans

OBJECTIVES

To determine genetic differences between agriculturalist and hunter-gatherer southern Native American populations for selected metabolism-related markers and to test whether Neel's thrifty genotype hypothesis (TGH) could explain the genetic patterns observed in these populations.

MATERIALS AND METHODS

375 Native South American individuals from 17 populations were genotyped using six markers (APOE rs429358 and rs7412; APOA2 rs5082; CD36 rs3211883; TCF7L2 rs11196205; and IGF2BP2 rs11705701). Additionally, APOE genotypes from 39 individuals were obtained from the literature. AMOVA, main effects, and gene-gene interaction tests were performed.

RESULTS

We observed differences in allele distribution patterns between agriculturalists and hunter-gatherers for some markers. For instance, between-groups component of genetic variance (F_CT ) for APOE rs429358 showed strong differences in allelic distributions between hunter-gatherers and agriculturalists (p = 0.00196). Gene-gene interaction analysis indicated that the APOE E4/CD36 TT and APOE E4/IGF2BP2 A carrier combinations occur at a higher frequency in hunter-gatherers, but this combination is not replicated in archaic (Neanderthal and Denisovan) and ancient (Anzick, Saqqaq, Ust-Ishim, Mal'ta) hunter-gatherer individuals.

DISCUSSION

A complex scenario explains the observed frequencies of the tested markers in hunter-gatherers. Different factors, such as pleotropic alleles, rainforest selective pressures, and population dynamics, may be collectively shaping the observed genetic patterns. We conclude that although TGH seems a plausible hypothesis to explain part of the data, other factors may be important in our tested populations.

Correction: Socioeconomic Status Is Not Related with Facial Fluctuating Asymmetry: Evidence from Latin-American Populations

[This corrects the article DOI: 10.1371/journal.pone.0169287.].

Socioeconomic Status Is Not Related with Facial Fluctuating Asymmetry: Evidence from Latin-American Populations

The expression of facial asymmetries has been recurrently related with poverty and/or disadvantaged socioeconomic status. Departing from the developmental instability theory, previous approaches attempted to test the statistical relationship between the stress experienced by individuals grown in poor conditions and an increase in facial and corporal asymmetry. Here we aim to further evaluate such hypothesis on a large sample of admixed Latin Americans individuals by exploring if low socioeconomic status individuals tend to exhibit greater facial fluctuating asymmetry values. To do so, we implement Procrustes analysis of variance and Hierarchical Linear Modelling (HLM) to estimate potential associations between facial fluctuating asymmetry values and socioeconomic status. We report significant relationships between facial fluctuating asymmetry values and age, sex, and genetic ancestry, while socioeconomic status failed to exhibit any strong statistical relationship with facial asymmetry. These results are persistent after the effect of heterozygosity (a proxy for genetic ancestry) is controlled in the model. Our results indicate that, at least on the studied sample, there is no relationship between socioeconomic stress (as intended as low socioeconomic status) and facial asymmetries.

Genetic evidence for two founding populations of the Americas

Genetic studies have been consistent with a single common origin of Native American groups from Central and South America^1-4. However, some morphological studies have suggested a more complex picture, whereby the northeast Asian affinities of present-day Native Americans contrast with a distinctive morphology seen in some of the earliest American skeletons, which share traits with present-day Australasians (indigenous groups in Australia, Melanesia, and island southeast Asia)^5-8. Here we analyze genome-wide data to show that some Amazonian Native Americans descend partly from a Native American founding population that carried ancestry more closely related to indigenous Australians, New Guineans and Andaman Islanders than to any present-day Eurasians or Native Americans. This signature is not present to the same extent or at all in present-day Northern and Central Americans or a ~12,600 year old Clovis genome, suggesting a more diverse set of founding populations of the Americas than previously accepted.

Ancient remains and the first peopling of the Americas: Reassessing the Hoyo Negro skull

OBJECTIVE

A noticeably well-preserved ∼12.500 years-old skeleton from the Hoyo Negro cave, Yucatán, México, was recently reported, along with its archaeological, genetic and skeletal characteristics. Based exclusively on an anatomical description of the skull (HN5/48), Chatters and colleagues stated that this specimen can be assigned to a set of ancient remains that differ from modern Native Americans, the so called "Paleoamericans". Here, we aim to further explore the morphological affinities of this specimen with a set of comparative cranial samples covering ancient and modern periods from Asia and the Americas.

METHODS

Images published in the original article were analyzed using geometric morphometrics methods. Shape variables were used to perform Principal Component and Discriminant analysis against the reference samples.

RESULTS

Even thought the Principal Component Analysis suggests that the Hoyo Negro skull falls in a subregion of the morphospace occupied by both "Paleoamericans" and some modern Native Americans, the Discriminant analyses suggest greater affinity with a modern Native American sample.

DISCUSSION

These results reinforce the idea that the original population that first occupied the New World carried high levels of within-group variation, which we have suggested previously on a synthetic model for the settlement of the Americas. Our results also highlight the importance of developing formal classificatory test before deriving settlement hypothesis purely based on macroscopic descriptions.

Self-assessment of color categories and its relationship with HLA profiling in Brazilian bone marrow donors

The Brazil Ministry of Health maintains a Registry of Bone Marrow Donors that corresponds to approximately 12% of the Bone Marrow Donors Worldwide registry. This registry contains information on ethnicity (by self-assessment of color) and HLA-A, -B, and -DRB1 type. The self-assessment of color tool has been extensively used for admixed population characterization. In this context, Brazil represents a highly admixed population, resulting from 5 centuries of colonization and interbreeding, mainly, but not exclusively, among Native Americans, Europeans, and Africans. Here we evaluated self-assessed skin color and HLA genetic information from 71,291 bone marrow donors of southern Brazil to verify how likely is the HLA profiling correspondence within and between self-assessed color groups. We found that HLA itself was a better ancestry indicator than was self-assessed color. Therefore, self-assessment of color in highly admixed populations, such as that of Brazil, is not indicative of higher correspondence in the HLA profiles within skin color groups.

Implications of the admixture process in skin color molecular assessment

The understanding of the complex genotype-phenotype architecture of human pigmentation has clear implications for the evolutionary history of humans, as well as for medical and forensic practices. Although dozens of genes have previously been associated with human skin color, knowledge about this trait remains incomplete. In particular, studies focusing on populations outside the European-North American axis are rare, and, until now, admixed populations have seldom been considered. The present study was designed to help fill this gap. Our objective was to evaluate possible associations of 18 single nucleotide polymorphisms (SNPs), located within nine genes, and one pseudogene with the Melanin Index (MI) in two admixed Brazilian populations (Gaucho, N = 352; Baiano, N = 148) with different histories of geographic and ethnic colonization. Of the total sample, four markers were found to be significantly associated with skin color, but only two (SLC24A5 rs1426654, and SLC45A2 rs16891982) were consistently associated with MI in both samples (Gaucho and Baiano). Therefore, only these 2 SNPs should be preliminarily considered to have forensic significance because they consistently showed the association independently of the admixture level of the populations studied. We do not discard that the other two markers (HERC2 rs1129038 and TYR rs1126809) might be also relevant to admixed samples, but additional studies are necessary to confirm the real importance of these markers for skin pigmentation. Finally, our study shows associations of some SNPs with MI in a modern Brazilian admixed sample, with possible applications in forensic genetics. Some classical genetic markers in Euro-North American populations are not associated with MI in our sample. Our results point out the relevance of considering population differences in selecting an appropriate set of SNPs as phenotype predictors in forensic practice.

Ocular and craniofacial phenotypes in a large Brazilian family with congenital aniridia

Congenital aniridia is a rare genetic disorder characterized by varying degrees of iris hypoplasia that are associated with additional ocular abnormalities. More than 90% of the causal mutations identified are found in the PAX6 gene, a transcription factor of critical importance in the process of neurogenesis and ocular development. Here, we investigate clinical, molecular, and craniofacial features of a large Brazilian family with congenital aniridia. Among the 56 eyes evaluated, phenotype variation encompassed bilateral total aniridia to mild iris defects with extensive variation between eyes of the same individual. PAX6 molecular screening indicated a heterozygous splice mutation (c.141 + 1G>A). Thus, we hypothesize that this splicing event may cause variation in the expression of the wild-type transcript, which may lead to the observed variation in phenotype. Affected individuals were more brachycephalic, even though their face height and cephalic circumference were not significantly different when compared to those of non-affected relatives. From this, we infer that the head shape of affected subjects may also be a result of the PAX6 splice-site mutation. Our data summarize the clinical variability associated with the ocular phenotype in a large family with aniridia, and help shed light on the role of PAX6 in neurocranial development.

The cognitive ability of extinct hominins: bringing down the hierarchy using genomic evidences

BACKGROUND

The availability of the full genomes of Homo sapiens, Homo neanderthalensis, and Denisovans, as well as modern bioinformatic tools, are opening new possibilities for the understanding of the differences and similarities present in these taxa.

METHODS

We searched for cognitive genes, examined their status in the genomes of these three entities. All substitutions present among them were retrieved.

RESULTS

We found 93 nonsynonymous substitutions in 51 cognitive genes, in which the derived allele was present in archaic and modern humans and the ancestral allele in other nonhuman primates.

CONCLUSIONS

The general picture obtained is of similarity in cognitive genes between extinct and extant humans.

Brief communication: 5-HTTLPR genetic diversity and mode of subsistence in Native Americans

The relationship between the "individualism-collectivism" and the serotonin transporter functional polymorphism (5-HTTLPR), suggested in the previous reports, was tested in Native South Amerindian populations. A total of 170 individuals from 21 populations were genotyped for the 5-HTTLPR alleles. For comparative purposes, these populations were classified as individualistic (recent history of hunter-gathering) or collectivistic (agriculturalists). These two groups showed an almost identical S allele frequency (75 and 76%, respectively). The analysis of molecular variance showed no structural differences between them. Behavioral typologies like those suggested by JY Chiao and KD Blizinsky (Proc R Soc B 277 () 529-537) are always a simplification of complex phenomena and should be regarded with caution. In addition, classification of a whole nation in the individualist/collectivist dichotomy is controversial. The focus on modes of subsistence in preindustrial societies, as was tested here, may be a good alternative although the postulated association between the 5-HTTLPR S allele and the collectivist societies was not confirmed.

A bayesian approach to genome/linguistic relationships in native South Americans

The relationship between the evolution of genes and languages has been studied for over three decades. These studies rely on the assumption that languages, as many other cultural traits, evolve in a gene-like manner, accumulating heritable diversity through time and being subjected to evolutionary mechanisms of change. In the present work we used genetic data to evaluate South American linguistic classifications. We compared discordant models of language classifications to the current Native American genome-wide variation using realistic demographic models analyzed under an Approximate Bayesian Computation (ABC) framework. Data on 381 STRs spread along the autosomes were gathered from the literature for populations representing the five main South Amerindian linguistic groups: Andean, Arawakan, Chibchan-Paezan, Macro-Jê, and Tupí. The results indicated a higher posterior probability for the classification proposed by J.H. Greenberg in 1987, although L. Campbell's 1997 classification cannot be ruled out. Based on Greenberg's classification, it was possible to date the time of Tupí-Arawakan divergence (2.8 kya), and the time of emergence of the structure between present day major language groups in South America (3.1 kya).

Lack of support for the association between facial shape and aggression: a reappraisal based on a worldwide population genetics perspective

Antisocial and criminal behaviors are multifactorial traits whose interpretation relies on multiple disciplines. Since these interpretations may have social, moral and legal implications, a constant review of the evidence is necessary before any scientific claim is considered as truth. A recent study proposed that men with wider faces relative to facial height (fWHR) are more likely to develop unethical behaviour mediated by a psychological sense of power. This research was based on reports suggesting that sexual dimorphism and selection would be responsible for a correlation between fWHR and aggression. Here we show that 4,960 individuals from 94 modern human populations belonging to a vast array of genetic and cultural contexts do not display significant amounts of fWHR sexual dimorphism. Further analyses using populations with associated ethnographical records as well as samples of male prisoners of the Mexico City Federal Penitentiary condemned by crimes of variable level of inter-personal aggression (homicide, robbery, and minor faults) did not show significant evidence, suggesting that populations/individuals with higher levels of bellicosity, aggressive behaviour, or power-mediated behaviour display greater fWHR. Finally, a regression analysis of fWHR on individual's fitness showed no significant correlation between this facial trait and reproductive success. Overall, our results suggest that facial attributes are poor predictors of aggressive behaviour, or at least, that sexual selection was weak enough to leave a signal on patterns of between- and within-sex and population facial variation.

Homo sapiens, Homo neanderthalensis and the Denisova specimen: New insights on their evolutionary histories using whole-genome comparisons

After a brief review of the most recent findings in the study of human evolution, an extensive comparison of the complete genomes of our nearest relative, the chimpanzee (Pan troglodytes), of extant Homo sapiens, archaic Homo neanderthalensis and the Denisova specimen were made. The focus was on non-synonymous mutations, which consequently had an impact on protein levels and these changes were classified according to degree of effect. A total of 10,447 non-synonymous substitutions were found in which the derived allele is fixed or nearly fixed in humans as compared to chimpanzee. Their most frequent location was on chromosome 21. Their presence was then searched in the two archaic genomes. Mutations in 381 genes would imply radical amino acid changes, with a fraction of these related to olfaction and other important physiological processes. Eight new alleles were identified in the Neanderthal and/or Denisova genetic pools. Four others, possibly affecting cognition, occured both in the sapiens and two other archaic genomes. The selective sweep that gave rise to Homo sapiens could, therefore, have initiated before the modern/archaic human divergence.

Fibroblast growth factor receptor 1 (FGFR1) variants and craniofacial variation in Amerindians and related populations

OBJECTIVES

The polymorphic site rs4647905 of the FGFR1 gene was previously associated with a decrease in cephalic index (CI). Here, we evaluate the relationships between genotypes and cephalometric measurements and indices in one Mexican Native and two mestizo Mexican populations using two haplotype-tag SNPs (rs4647905 and rs3213849) that represent >85% of the FGFR1 variability, plus three other SNPs (rs2293971, rs2304000, and rs930828) situated nearby. In addition, we genotyped five South American natives, two European, one African, and one Siberian populations to evaluate their intra and intercontinental population diversity.

METHODS

The five SNPs were tested and the craniofacial measurements and indices were collected using standardized procedures. Principal Component Analysis was used to verify individual/population comparisons. Associations were performed through the generalized linear model (GLM), coefficient of determination R(2) and linear regression tests.

RESULTS

We found a tendency for a decrease in CI in individuals homozygous for allele rs4647905C, regardless of the population to which they belong, though the effect is more pronounced in mestizo. When the GLM analyses were performed using the absolute/linear cephalometric measurements, a statistically significant association was found between four SNPs and head length in the mestizo population.

CONCLUSIONS

FGFR1 polymorphisms, especially rs4647905, can have an important role in the normal human skull variation, primarily due to their influence in head length, which would affect other cephalometric absolute/linear measures as well as indices like CI as a result of the pervasive nature of the morphological integration that characterizes the human skull.

Evolutionary responses to a constructed niche: ancient Mesoamericans as a model of gene-culture coevolution

Culture and genetics rely on two distinct but not isolated transmission systems. Cultural processes may change the human selective environment and thereby affect which individuals survive and reproduce. Here, we evaluated whether the modes of subsistence in Native American populations and the frequencies of the ABCA1*Arg230Cys polymorphism were correlated. Further, we examined whether the evolutionary consequences of the agriculturally constructed niche in Mesoamerica could be considered as a gene-culture coevolution model. For this purpose, we genotyped 229 individuals affiliated with 19 Native American populations and added data for 41 other Native American groups (n = 1905) to the analysis. In combination with the SNP cluster of a neutral region, this dataset was then used to unravel the scenario involved in 230Cys evolutionary history. The estimated age of 230Cys is compatible with its origin occurring in the American continent. The correlation of its frequencies with the archeological data on Zea pollen in Mesoamerica/Central America, the neutral coalescent simulations, and the F_ST-based natural selection analysis suggest that maize domestication was the driving force in the increase in the frequencies of 230Cys in this region. These results may represent the first example of a gene-culture coevolution involving an autochthonous American allele.

Predicting Homo pigmentation phenotype through genomic data: from Neanderthal to James Watson

BACKGROUND

Human pigmentation is regulated by several genes acting at different stages of melanin formation. Functional and association studies have elucidated the role of several of these genes in pigmentation phenotypes. Forensic and evolutionary studies can benefit from this knowledge.

OBJECTIVES

To evaluate the reliability of the prediction of pigmentation phenotypes using a large database of genetic markers in individuals with known phenotypes; and from this try to predict the pigmentation phenotype of prehistoric Homo specimens and of contemporary individuals whose visible phenotypes are not known.

METHODS

We compared predicted and observed phenotypic data through an analysis of 124 single nucleotide polymorphisms in 33 genic and seven intergenic regions of 30 subjects, five of them prehistoric, whose complete nuclear genomes are available in UCSC and PSU UCSC public databases.

RESULTS

For the molecular predicted versus observed phenotypes, the percentage of agreement was as follows: freckles: 91; skin: 64; hair: 44; eyes: 36; total: 59; while the molecular predicted versus probable (no visible observation available; inferences based on ethnic population characteristics) it was, respectively, 83, 60, 42, 67, and 63. The difference between two sets is statistically nonsignificant (P = 0.75).

CONCLUSION

To our knowledge, this is the first article to examine the effect of a large number of genetics markers for phenotype prediction. The approach could be useful for forensic applications, as well as for the determination of possible phenotypes of extinct prehistoric individuals.

A functional ABCA1 gene variant is associated with low HDL-cholesterol levels and shows evidence of positive selection in Native Americans

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.