NNT mediates redox-dependent pigmentation via a UVB- and MITF-independent mechanism.

Ultraviolet (UV) light and incompletely understood genetic and epigenetic variations determine skin color. Here we describe an UV- and microphthalmia-associated transcription factor (MITF)-independent mechanism of skin pigmentation. Targeting the mitochondrial redox-regulating enzyme nicotinamide nucleotide transhydrogenase (NNT) resulted in cellular redox changes that affect tyrosinase degradation. These changes regulate melanosome maturation and, consequently, eumelanin levels and pigmentation. Topical application of small-molecule inhibitors yielded skin darkening in human skin, and mice with decreased NNT function displayed increased pigmentation. Additionally, genetic modification of NNT in zebrafish alters melanocytic pigmentation. Analysis of four diverse human cohorts revealed significant associations of skin color, tanning, and sun protection use with various single-nucleotide polymorphisms within NNT. NNT levels were independent of UVB irradiation and redox modulation. Individuals with postinflammatory hyperpigmentation or lentigines displayed decreased skin NNT levels, suggesting an NNT-driven, redox-dependent pigmentation mechanism that can be targeted with NNT-modifying topical drugs for medical and cosmetic purposes.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci.

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

Clotting factor genes are associated with preeclampsia in high-altitude pregnant women in the Peruvian Andes.

Preeclampsia is a multi-organ complication of pregnancy characterized by sudden hypertension and proteinuria that is among the leading causes of preterm delivery and maternal morbidity and mortality worldwide. The heterogeneity of preeclampsia poses a challenge for understanding its etiology and molecular basis. Intriguingly, risk for the condition increases in high-altitude regions such as the Peruvian Andes. To investigate the genetic basis of preeclampsia in a population living at high altitude, we characterized genome-wide variation in a cohort of preeclamptic and healthy Andean families (n = 883) from Puno, Peru, a city located above 3,800 meters of altitude. Our study collected genomic DNA and medical records from case-control trios and duos in local hospital settings. We generated genotype data for 439,314 SNPs, determined global ancestry patterns, and mapped associations between genetic variants and preeclampsia phenotypes. A transmission disequilibrium test (TDT) revealed variants near genes of biological importance for placental and blood vessel function. The top candidate region was found on chromosome 13 of the fetal genome and contains clotting factor genes PROZ, F7, and F10. These findings provide supporting evidence that common genetic variants within coagulation genes play an important role in preeclampsia. A selection scan revealed a potential adaptive signal around the ADAM12 locus on chromosome 10, implicated in pregnancy disorders. Our discovery of an association in a functional pathway relevant to pregnancy physiology in an understudied population of Native American origin demonstrates the increased power of family-based study design and underscores the importance of conducting genetic research in diverse populations.

Inbreeding, Native American ancestry and child mortality: linking human selection and paediatric medicine.

The children of related parents show increased risk of early mortality. The Native American genome typically exhibits long stretches of homozygosity, and Latin Americans are highly heterogeneous regarding the individual burden of homozygosity, the proportion and the type of Native American ancestry. We analysed nationwide mortality and genome-wide genotype data from admixed Chileans to investigate the relationship between common causes of child mortality, homozygosity and Native American ancestry. Results from two-stage linear-Poisson regression revealed a strong association between the sum length of runs of homozygosity (SROH) above 1.5 Megabases (Mb) in each genome and mortality due to intracranial non-traumatic haemorrhage of foetus and newborn (5% increased risk of death per Mb in SROH, P = 1 × 10-3) and disorders related to short gestation and low birth weight (P = 3 × 10-4). The major indigenous populations in Chile are Aymara-Quechua in the north of the country and the Mapuche-Huilliche in the south. The individual proportion of Aymara-Quechua ancestry was associated with an increased risk of death due to anencephaly and similar malformations (P = 4 × 10-5), and the risk of death due to Edwards and Patau trisomy syndromes decreased 4% per 1% Aymara-Quechua ancestry proportion (P = 4 × 10-4) and 5% per 1% Mapuche-Huilliche ancestry proportion (P = 2 × 10-3). The present results suggest that short gestation, low birth weight and intracranial non-traumatic haemorrhage mediate the negative effect of inbreeding on human selection. Independent validation of the identified associations between common causes of child death, homozygosity and fine-scale ancestry proportions may inform paediatric medicine.

Deep genetic affinity between coastal Pacific and Amazonian natives evidenced by Australasian ancestry.

Different models have been proposed to elucidate the origins of the founding populations of America, along with the number of migratory waves and routes used by these first explorers. Settlements, both along the Pacific coast and on land, have been evidenced in genetic and archeological studies. However, the number of migratory waves and the origin of immigrants are still controversial topics. Here, we show the Australasian genetic signal is present in the Pacific coast region, indicating a more widespread signal distribution within South America and implicating an ancient contact between Pacific and Amazonian dwellers. We demonstrate that the Australasian population contribution was introduced in South America through the Pacific coastal route before the formation of the Amazonian branch, likely in the ancient coastal Pacific/Amazonian population. In addition, we detected a significant amount of interpopulation and intrapopulation variation in this genetic signal in South America. This study elucidates the genetic relationships of different ancestral components in the initial settlement of South America and proposes that the migratory route used by migrants who carried the Australasian ancestry led to the absence of this signal in the populations of Central and North America.

Adaptive strategies of aquatic mammals: Exploring the role of the HIF pathway and hypoxia tolerance.

Aquatic mammals (marine and freshwater species) share significant and similar adaptations, enabling them to tolerate hypoxia during regular breath-hold diving. Despite the established importance of HIF1A, a master regulator in the molecular mechanism of hypoxia response, and other associated genes, their role in the evolutionary adaptation of aquatic mammals is not fully understood. In this study, we investigated this topic by employing a candidate gene approach to analyze 11 critical genes involved in the HIF1A signaling pathway in aquatic mammals. Our gene analyses included evaluating positive and negative selection, relaxation or constriction of selection, and molecular convergence compared to other terrestrial mammals, including subterranean mammals. Evidence of selection suggested a significant role of negative selection, as well as relaxation of the selective regime in cetaceans for most of these genes. We found that the glutamine 68 variant in the HIF3α protein is unique to cetaceans and initial evaluations indicated a destabilizing effect on protein structure. However, further analyses are necessary to evaluate its functional impact and adaptive relevance in this taxon.

Oral microbiota, co-evolution, and implications for health and disease: The case of indigenous peoples.

Evidence indicates that oral microbiota plays a crucial role in human health and disease. For instance, diseases with multifactorial etiology, such as periodontitis and caries, which cause a detrimental impact on human well-being and health, can be caused by alterations in the host-microbiota interactions, where non-pathogenic bacteria give way to pathogenic orange/red-complex bacterial species (a change from a eubiotic to dysbiotic state). In this scenario, where thousands of oral microorganisms, including fungi, archaea, and phage species, and their host are co-evolving, a set of phenomena, such as the arms race and Red or Black Queen dynamics, are expected to operate. We review concepts on the subject and revisit the nature of bacterial complexes linked to oral health and diseases, as well as the problem of the bacterial resistome in the face of the use of antibiotics and what is the impact of this on the evolutionary trajectory of the members of this symbiotic ecosystem. We constructed a 16SrRNA tree to show that adaptive consortia of oral bacterial complexes do not necessarily rescue phylogenetic relationships. Finally, we remember that oral health is not exempt from health disparity trends in some populations, such as Native Americans, when compared with non-Indigenous people.

Development and internal validation of a multifactorial risk prediction model for gallbladder cancer in a high-incidence country.

Since 2006, Chile has been implementing a gallbladder cancer (GBC) prevention program based on prophylactic cholecystectomy for gallstone patients aged 35 to 49 years. The effectiveness of this prevention program has not yet been comprehensively evaluated. We conducted a retrospective study of 473 Chilean GBC patients and 2137 population-based controls to develop and internally validate three GBC risk prediction models. The Baseline Model accounted for gallstones while adjusting for sex and birth year. Enhanced Model I also included the non-genetic risk factors: body mass index, educational level, Mapuche surnames, number of children and family history of GBC. Enhanced Model II further included Mapuche ancestry and the genotype for rs17209837. Multiple Cox regression was applied to assess the predictive performance, quantified by the area under the precision-recall curve (AUC-PRC) and the number of cholecystectomies needed (NCN) to prevent one case of GBC at age 70 years. The AUC-PRC for the Baseline Model (0.44%, 95%CI 0.42-0.46) increased by 0.22 (95%CI 0.15-0.29) when non-genetic factors were included, and by 0.25 (95%CI 0.20-0.30) when incorporating non-genetic and genetic factors. The overall NCN for Chileans with gallstones (115, 95%CI 104-131) decreased to 92 (95%CI 60-128) for Chileans with a higher risk than the median according to Enhanced Model I, and to 80 (95%CI 59-110) according to Enhanced Model II. In conclusion, age, sex and gallstones are strong risk factors for GBC, but consideration of other non-genetic factors and individual genotype data improves risk prediction and may optimize allocation of financial resources and surgical capacity.

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.

Genomic insight into the origins and dispersal of the Brazilian coastal natives.

In the 15th century, ∼900,000 Native Americans, mostly Tupí speakers, lived on the Brazilian coast. By the end of the 18th century, the coastal native populations were declared extinct. The Tupí arrived on the east coast after leaving the Amazonian basin ∼2,000 y before present; however, there is no consensus on how this migration occurred: toward the northern Amazon and then directly to the Atlantic coast, or heading south into the continent and then migrating to the coast. Here we leveraged genomic data from one of the last remaining putative representatives of the Tupí coastal branch, a small, admixed, self-reported Tupiniquim community, as well as data of a Guaraní Mbyá native population from Southern Brazil and of three other native populations from the Amazonian region. We demonstrated that the Tupiniquim Native American ancestry is not related to any extant Brazilian Native American population already studied, and thus they could be considered the only living representatives of the extinct Tupí branch that used to settle the Atlantic Coast of Brazil. Furthermore, these data show evidence of a direct migration from Amazon to the Northeast Coast in pre-Columbian time, giving rise to the Tupí Coastal populations, and a single distinct migration southward that originated the Guaraní people from Brazil and Paraguay. This study elucidates the population dynamics and diversification of the Brazilian natives at a genomic level, which was made possible by recovering data from the Brazilian coastal population through the genomes of mestizo individuals.

Exploring the diversity of AVPR2 in Primates and its evolutionary implications.

The current study focuses on the investigation of AVPR2 (VTR2C) protein-coupled receptor variants specific to different primate taxa. AVPR2 is activated by the neurohormone AVP, which modulates physiological processes, including water homeostasis. Our findings reveal positive selection at three AVPR2 sites at positions 190, 250, and 346. Variation at position 250 is associated with human Congenital Nephrogenic Diabetes Insipidus (cNDI), a condition characterized by excessive water loss. Other 13 functional sites with potential adaptive relevance include positions 185, 202, 204, and 252 associated with cNDI. We identified SH3-binding motifs in AVPR2's ICL3 and N-terminus domains, with some losses observed in clades of Cercopithecidae, Callitrichinae, and Atelidae. SH3-binding motifs are crucial in regulating cellular physiology, indicating that the differences may be adaptive. Co-evolution was found between AVPR2 residues and those in the AVP signal peptide/Neurophysin-2 and AQP2, other molecules in the same signaling cascade. No significant correlation was found between these Primates' taxon-specific variants and the bioclimatic variables of the areas where they live. Distinct co-evolving amino acid sequences in functional sites were found in Platyrrhini and Catarrhini, which may have adaptive implications involving glucocorticoid hormones, suggesting varied selective pressures. Further studies are required to confirm these results.

Insights into the evolutionary history of the most skilled tool-handling platyrrhini monkey: Sapajus libidinosus from the Serra da Capivara National Park.

Sapajus libidinosus members of the Pedra Furada group, living in the Serra da Capivara National Park, use stone tools in a wider variety of behaviors than any other living animal, except humans. To rescue the evolutionary history of the Caatinga S. libidinosus and identify factors that may have contributed to the emergence and maintenance of their tool-use culture, we conducted fieldwork seasons to obtain biological samples of these capuchin monkeys. UsingCYTBsequences, we show a discrete but constant population growth from the beginning of the Holocene to the present, overlapping the emergence of the Caatinga biome. Our habitat suitability reconstruction reports the presence of plants whose hard fruits, seeds, or roots are processed by capuchins using tools. TheS. libidinosusindividuals in the Caatinga were capable of dynamically developing and maintaining their autochthonous culture thanks to: a) cognitive capacity to generate and execute innovation under selective pressure; b) tolerance favoring learning and cultural inheritance; c) an unknown genetic repertoire that underpins the adaptive traits; d) a high degree of terrestriality; e) presence and abundance of natural resources, which makes some places "hot spots" for innovation, and cultural diversification within a relatively short time.

Ionotropic Receptors as a Driving Force behind Human Synapse Establishment.

The origin of nervous systems is a main theme in biology and its mechanisms are largely underlied by synaptic neurotransmission. One problem to explain synapse establishment is that synaptic orthologs are present in multiple aneural organisms. We questioned how the interactions among these elements evolved and to what extent it relates to our understanding of the nervous systems complexity. We identified the human neurotransmission gene network based on genes present in GABAergic, glutamatergic, serotonergic, dopaminergic, and cholinergic systems. The network comprises 321 human genes, 83 of which act exclusively in the nervous system. We reconstructed the evolutionary scenario of synapse emergence by looking for synaptic orthologs in 476 eukaryotes. The Human-Cnidaria common ancestor displayed a massive emergence of neuroexclusive genes, mainly ionotropic receptors, which might have been crucial to the evolution of synapses. Very few synaptic genes had their origin after the Human-Cnidaria common ancestor. We also identified a higher abundance of synaptic proteins in vertebrates, which suggests an increase in the synaptic network complexity of those organisms.

Evolutionary history of the SARS-CoV-2 Gamma variant of concern (P.1): a perfect storm.

Our goal was to describe in more detail the evolutionary history of Gamma and two derived lineages (P.1.1 and P.1.2), which are part of the arms race that SARS-CoV-2 wages with its host. A total of 4,977 sequences of the Gamma strain of SARS-CoV-2 from Brazil were analyzed. We detected 194 sites under positive selection in 12 genes/ORFs: Spike, N, M, E, ORF1a, ORF1b, ORF3, ORF6, ORF7a, ORF7b, ORF8, and ORF10. Some diagnostic sites for Gamma lacked a signature of positive selection in our study, but these were not fixed, apparently escaping the action of purifying selection. Our network analyses revealed branches leading to expanding haplotypes with sites under selection only detected when P.1.1 and P.1.2 were considered. The P.1.2 exclusive haplotype H_5 originated from a non-synonymous mutational step (H3509Y) in H_1 of ORF1a. The selected allele, 3509Y, represents an adaptive novelty involving ORF1a of P.1. Finally, we discuss how phenomena such as epistasis and antagonistic pleiotropy could limit the emergence of new alleles (and combinations thereof) in SARS-COV-2 lineages, maintaining infectivity in humans, while providing rapid response capabilities to face the arms race triggered by host immuneresponses.

Origins, Admixture Dynamics, and Homogenization of the African Gene Pool in the Americas.

The Transatlantic Slave Trade transported more than 9 million Africans to the Americas between the early 16th and the mid-19th centuries. We performed a genome-wide analysis using 6,267 individuals from 25 populations to infer how different African groups contributed to North-, South-American, and Caribbean populations, in the context of geographic and geopolitical factors, and compared genetic data with demographic history records of the Transatlantic Slave Trade. We observed that West-Central Africa and Western Africa-associated ancestry clusters are more prevalent in northern latitudes of the Americas, whereas the South/East Africa-associated ancestry cluster is more prevalent in southern latitudes of the Americas. This pattern results from geographic and geopolitical factors leading to population differentiation. However, there is a substantial decrease in the between-population differentiation of the African gene pool within the Americas, when compared with the regions of origin from Africa, underscoring the importance of historical factors favoring admixture between individuals with different African origins in the New World. This between-population homogenization in the Americas is consistent with the excess of West-Central Africa ancestry (the most prevalent in the Americas) in the United States and Southeast-Brazil, with respect to historical-demography expectations. We also inferred that in most of the Americas, intercontinental admixture intensification occurred between 1750 and 1850, which correlates strongly with the peak of arrivals from Africa. This study contributes with a population genetics perspective to the ongoing social, cultural, and political debate regarding ancestry, admixture, and the mestizaje process in the Americas.

Genetic evidence for two founding populations of the Americas.

Genetic studies have consistently indicated a single common origin of Native American groups from Central and South America. 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). Here we analyse 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 in a ∼12,600-year-old Clovis-associated genome, suggesting a more diverse set of founding populations of the Americas than previously accepted.

Molecular evolutionary insights from PRLR in mammals.

Prolactin (PRL) is a pleiotropic neurohormone secreted by the mammalian pituitary gland into the blood, thus reaching many tissues and organs beyond the brain. PRL binds to its receptor, PRLR, eliciting a molecular signaling cascade. This system modulates essential mammalian behaviors and promotes notable modifications in the reproductive female tissues and organs. Here, we explore how the intracellular domain of PRLR (PRLR-ICD) modulates the expression of the PRLR gene. Despite differences in the reproductive strategies between eutherian and metatherian mammals, there is no clear distinction between PRLR-ICD functional motifs. However, we found selection signatures that showed differences between groups, with many conserved functional elements strongly maintained through purifying selection across the class Mammalia. We observed a few residues under relaxed selection, the levels of which were more pronounced in Eutheria and particularly striking in primates (Simiiformes), which could represent a pre-adaptive genetic element protected from purifying selection. Alternative, new motifs, such as YLDP (318-321) and others with residues Y283 and Y290, may already be functional. These motifs would have been co-opted in primates as part of a complex genetic repertoire related to some derived adaptive phenotypes, but these changes would have no impact on the primordial functions that characterize the mammals as a whole and that are related to the PRL-PRLR system.

Genome-wide association studies and CRISPR/Cas9-mediated gene editing identify regulatory variants influencing eyebrow thickness in humans.

Hair plays an important role in primates and is clearly subject to adaptive selection. While humans have lost most facial hair, eyebrows are a notable exception. Eyebrow thickness is heritable and widely believed to be subject to sexual selection. Nevertheless, few genomic studies have explored its genetic basis. Here, we performed a genome-wide scan for eyebrow thickness in 2961 Han Chinese. We identified two new loci of genome-wide significance, at 3q26.33 near SOX2 (rs1345417: P = 6.51×10(-10)) and at 5q13.2 near FOXD1 (rs12651896: P = 1.73×10(-8)). We further replicated our findings in the Uyghurs, a population from China characterized by East Asian-European admixture (N = 721), the CANDELA cohort from five Latin American countries (N = 2301), and the Rotterdam Study cohort of Dutch Europeans (N = 4411). A meta-analysis combining the full GWAS results from the three cohorts of full or partial Asian descent (Han Chinese, Uyghur and Latin Americans, N = 5983) highlighted a third signal of genome-wide significance at 2q12.3 (rs1866188: P = 5.81×10(-11)) near EDAR. We performed fine-mapping and prioritized four variants for further experimental verification. CRISPR/Cas9-mediated gene editing provided evidence that rs1345417 and rs12651896 affect the transcriptional activity of the nearby SOX2 and FOXD1 genes, which are both involved in hair development. Finally, suitable statistical analyses revealed that none of the associated variants showed clear signals of selection in any of the populations tested. Contrary to popular speculation, we found no evidence that eyebrow thickness is subject to strong selective pressure.

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.

Meta-analysis of genome-wide association studies identifies 8 novel loci involved in shape variation of human head hair.

Shape variation of human head hair shows striking variation within and between human populations, while its genetic basis is far from being understood. We performed a series of genome-wide association studies (GWASs) and replication studies in a total of 28 964 subjects from 9 cohorts from multiple geographic origins. A meta-analysis of three European GWASs identified 8 novel loci (1p36.23 ERRFI1/SLC45A1, 1p36.22 PEX14, 1p36.13 PADI3, 2p13.3 TGFA, 11p14.1 LGR4, 12q13.13 HOXC13, 17q21.2 KRTAP, and 20q13.33 PTK6), and confirmed 4 previously known ones (1q21.3 TCHH/TCHHL1/LCE3E, 2q35 WNT10A, 4q21.21 FRAS1, and 10p14 LINC00708/GATA3), all showing genome-wide significant association with hair shape (P < 5e-8). All except one (1p36.22 PEX14) were replicated with nominal significance in at least one of the 6 additional cohorts of European, Native American and East Asian origins. Three additional previously known genes (EDAR, OFCC1, and PRSS53) were confirmed at the nominal significance level. A multivariable regression model revealed that 14 SNPs from different genes significantly and independently contribute to hair shape variation, reaching a cross-validated AUC value of 0.66 (95% CI: 0.62-0.70) and an AUC value of 0.64 in an independent validation cohort, providing an improved accuracy compared with a previous model. Prediction outcomes of 2504 individuals from a multiethnic sample were largely consistent with general knowledge on the global distribution of hair shape variation. Our study thus delivers target genes and DNA variants for future functional studies to further evaluate the molecular basis of hair shape in humans.