The timing of pigmentation lightening in Europeans

Detection of human adaptation during the past 2000 years

Detection of recent natural selection is a challenging problem in population genetics. Here we introduce the singleton density score (SDS), a method to infer very recent changes in allele frequencies from contemporary genome sequences. Applied to data from the UK10K Project, SDS reflects allele frequency changes in the ancestors of modern Britons during the past ~2000 to 3000 years. We see strong signals of selection at lactase and the major histocompatibility complex, and in favor of blond hair and blue eyes. For polygenic adaptation, we find that recent selection for increased height has driven allele frequency shifts across most of the genome. Moreover, we identify shifts associated with other complex traits, suggesting that polygenic adaptation has played a pervasive role in shaping genotypic and phenotypic variation in modern humans.

Understanding rare and common diseases in the context of human evolution

The wealth of available genetic information is allowing the reconstruction of human demographic and adaptive history. Demography and purifying selection affect the purge of rare, deleterious mutations from the human population, whereas positive and balancing selection can increase the frequency of advantageous variants, improving survival and reproduction in specific environmental conditions. In this review, I discuss how theoretical and empirical population genetics studies, using both modern and ancient DNA data, are a powerful tool for obtaining new insight into the genetic basis of severe disorders and complex disease phenotypes, rare and common, focusing particularly on infectious disease risk.

Local Adaptation of Sun-Exposure-Dependent Gene Expression Regulation in Human Skin

Sun-exposure is a key environmental variable in the study of human evolution. Several skin-pigmentation genes serve as classical examples of positive selection, suggesting that sun-exposure has significantly shaped worldwide genomic variation. Here we investigate the interaction between genetic variation and sun-exposure, and how this impacts gene expression regulation. Using RNA-Seq data from 607 human skin samples, we identified thousands of transcripts that are differentially expressed between sun-exposed skin and non-sun-exposed skin. We then tested whether genetic variants may influence each individual’s gene expression response to sun-exposure. Our analysis revealed 10 sun-exposure-dependent gene expression quantitative trait loci (se-eQTLs), including genes involved in skin pigmentation (SLC45A2) and epidermal differentiation (RASSF9). The allele frequencies of the RASSF9 se-eQTL across diverse populations correlate with the magnitude of solar radiation experienced by these populations, suggesting local adaptation to varying levels of sunlight. These results provide the first examples of sun-exposure-dependent regulatory variation and suggest that this variation has contributed to recent human adaptation.

Varying levels of sun-exposure across the world have significantly shaped human evolution. Previous analyses have found several skin pigmentation genes with evidence of strong evolutionary pressures throughout human history, manifesting as large differences in the frequency of genomic variants across populations. But even within populations, individuals respond differently to sun-exposure, suggesting variation in addition to the major differences in skin pigmentation across populations. Here we investigated whether genetic variants associate with response to sun-exposure within Europeans. To measure the response we analyzed gene expression in sun-exposed and non-sun-exposed skin, and identified ten genetic variants that associated with the sun-exposure response of nearby genes. One of these genetic variants, which associated with the sun-exposure response of the gene RASSF9, showed evidence of adaptation in humans in response to solar radiation. Together this evidence suggests that the regulation of gene expression is influenced by sun-exposure and that the sun-exposure dependent effect on RASSF9 expression may have had an effect on human fitness. To our knowledge, this is the first example of an environment-dependent regulatory variant with evidence of adaptation in humans.

The African diaspora: history, adaptation and health

The trans-Atlantic slave trade brought millions of Africans to the New World. Advances in genomics are providing novel insights into the history and health of Africans and the diasporan populations. Recent examples reviewed here include the unraveling of substantial hunter-gatherer and 'Eurasian' admixtures across sub-Saharan Africa, expanding our understanding of ancestral African genetics; the global ubiquity of mixed ancestry; the revealing of African ancestry in Latin Americans that likely derived from the slave trade; and understanding of the ancestral backgrounds of APOL1 and LPL found to influence kidney disease and lipid levels, respectively, providing specific insights into disease etiology and health disparities.

Food and pathogen adaptations in the Angolan Namib desert: Tracing the spread of lactase persistence and human African trypanosomiasis resistance into southwestern Africa

OBJECTIVES

We investigated the frequency distribution and haplotype diversity of human African trypanosomiasis (HAT) resistance and lactase persistence (LP) variants in populations from the Angolan Namib to trace the spread of these genetic adaptations into southwestern Africa.

MATERIALS AND METHODS

We resequenced two fragments of the LCT enhancer and the APOL1 gene and genotyped flanking short tandem repeat loci in six groups with different subsistence traditions living in the Angolan Namib, and in a comparative dataset including other populations from Africa and Europe. LP in the Angolan Namib is represented by the -14010*C allele, which is associated with a predominant haplotype that is shared with other southern and eastern African populations. While LP was found to be more frequent in foragers than in pastoralists, the frequencies of the two APOL1 variants associated with HAT-resistance (G1 and G2) did not differ between the two groups. The G1 allele is mostly associated with a single widespread haplotype. The G2 allele is linked to several haplotypes that are molecularly related to haplotypes found in other African Bantu-speaking populations. The putatively archaic G3 variant displayed more intra-allelic diversity in Africa than in Europe.

DISCUSSION

The LP adaptation was carried to southern Africa by non-Bantu speaking pastoralists from eastern Africa, but an obvious link between its presence in southern Angola and groups speaking languages of the Khoe-Kwadi family, as previously found in other areas, could not be confirmed. The presence of APOL1 variants G1 and G2 is linked to the Bantu expansions. Our results suggest that the G3 variant was retained in modern humans by incomplete lineage sorting.

Basis for the gain and subsequent dilution of epidermal pigmentation during human evolution: The barrier and metabolic conservation hypotheses revisited

The evolution of human skin pigmentation must address both the initial evolution of intense epidermal pigmentation in hominins, and its subsequent dilution in modern humans. While many authorities believe that epidermal pigmentation evolved to protect against either ultraviolet B (UV-B) irradiation-induced mutagenesis or folic acid photolysis, we hypothesize that pigmentation augmented the epidermal barriers by shifting the UV-B dose-response curve from toxic to beneficial. Whereas erythemogenic UV-B doses produce apoptosis and cell death, suberythemogenic doses benefit permeability and antimicrobial function. Heavily melanized melanocytes acidify the outer epidermis and emit paracrine signals that augment barrier competence. Modern humans, residing in the cooler, wetter climes of south-central Europe and Asia, initially retained substantial pigmentation. While their outdoor lifestyles still permitted sufficient cutaneous vitamin D3 (VD3) synthesis, their marginal nutritional status, coupled with cold-induced caloric needs, selected for moderate pigment reductions that diverted limited nutritional resources towards more urgent priorities (=metabolic conservation). The further pigment-dilution that evolved as humans reached north-central Europe (i.e., northern France, Germany), likely facilitated cutaneous VD3 synthesis, while also supporting ongoing, nutritional requirements. But at still higher European latitudes where little UV-B breaches the atmosphere (i.e., present-day UK, Scandinavia, Baltic States), pigment dilution alone could not suffice. There, other nonpigment-related mutations evolved to facilitate VD3 production; for example, in the epidermal protein, filaggrin, resulting in reduced levels of its distal metabolite, trans-urocanic acid, a potent UV-B chromophore. Thus, changes in human pigmentation reflect a complex interplay between latitude, climate, diet, lifestyle, and shifting metabolic priorities.

Latitudinal Clines of the Human Vitamin D Receptor and Skin Color Genes

The well-documented latitudinal clines of genes affecting human skin color presumably arise from the need for protection from intense ultraviolet radiation (UVR) vs. the need to use UVR for vitamin D synthesis. Sampling 751 subjects from a broad range of latitudes and skin colors, we investigated possible multilocus correlated adaptation of skin color genes with the vitamin D receptor gene (VDR), using a vector correlation metric and network method called BlocBuster. We discovered two multilocus networks involving VDR promoter and skin color genes that display strong latitudinal clines as multilocus networks, even though many of their single gene components do not. Considered one by one, the VDR components of these networks show diverse patterns: no cline, a weak declining latitudinal cline outside of Africa, and a strong in- vs. out-of-Africa frequency pattern. We confirmed these results with independent data from HapMap. Standard linkage disequilibrium analyses did not detect these networks. We applied BlocBuster across the entire genome, showing that our networks are significant outliers for interchromosomal disequilibrium that overlap with environmental variation relevant to the genes’ functions. These results suggest that these multilocus correlations most likely arose from a combination of parallel selective responses to a common environmental variable and coadaptation, given the known Mendelian epistasis among VDR and the skin color genes.

Association of genetic variants with skin pigmentation phenotype among populations of west Maharashtra, India

OBJECTIVES

South Asians exhibit extensive variation in skin melanin index (MI) which is observed across the broader region of South Asia as well as within restricted geographic regions. However, the genetic variants associated with variation in the skin pigmentation phenotype are poorly understood in these populations. The present study examines the association between MI measures and genetic variants from 5 candidate pigmentation genes among 533 individuals representing 6 populations of West Maharashtra.

METHODS

Associations between five single nucleotide polymorphisms (SNPs) known to play a role in pigmentation (rs1426654-SLC24A5, rs1042602-TYR, rs16891982-SLC45A2, rs6058017-ASIP, and rs642742-KITLG) and MI measures were tested using standard one-way analysis of variance (ANOVA) within each population. Multiple linear regression was used to test the effects of these SNPs in the full West Maharashtra sample using sex, age, and population or social group as covariates.

RESULTS

rs1426654 showed significant association with MI in all six study populations (P < 0.01). Association tests using sex, age, and population as covariates showed rs1426654 and rs1042602 to be significantly (P < 0.01) associated with lighter skin pigmentation in West Maharashtra as a whole. By contrast, when social group was added as a covariate instead of population, rs1426654, rs1042602, and rs16891982 were significantly (P < 0.01) associated with lighter skin pigmentation.

CONCLUSIONS

Only rs1426654 is significantly associated with MI in each individual population; however, rs1426654, rs1042602, and rs16891982 are significantly associated with pigmentation in the broader West Maharashtra region after controlling for population and social group, with rs1426654 (SLC24A5) explaining the majority of the observed variation. Am. J. Hum. Biol. 28:610-618, 2016. © 2016 Wiley Periodicals, Inc.

Ancient DNA and the rewriting of human history: be sparing with Occam's razor

Ancient DNA research is revealing a human history far more complex than that inferred from parsimonious models based on modern DNA. Here, we review some of the key events in the peopling of the world in the light of the findings of work on ancient DNA.

Estimating the Ages of Selection Signals from Different Epochs in Human History

Genetic variation harbors signatures of natural selection driven by selective pressures that are often unknown. Estimating the ages of selection signals may allow reconstructing the history of environmental changes that shaped human phenotypes and diseases. We have developed an approximate Bayesian computation (ABC) approach to estimate allele ages under a model of selection on new mutations and under demographic models appropriate for human populations. We have applied it to two resequencing data sets: An ultra-high depth data set from a relatively small sample of unrelated individuals and a lower depth data set in a larger sample with transmission information. In addition to evaluating the accuracy of our method based on simulations, for each SNP, we assessed the consistency between the posterior probabilities estimated by the ABC approach and the ancient DNA record, finding good agreement between the two types of data and methods. Applying this ABC approach to data for eight single nucleotide polymorphisms (SNPs), we were able to rule out an onset of selection prior to the dispersal out-of-Africa for three of them and more recent than the spread of agriculture for an additional three SNPs.

Quantitative assessment of skin, hair, and iris variation in a diverse sample of individuals and associated genetic variation

OBJECTIVES

The main goals of this study are to 1) quantitatively measure skin, hair, and iris pigmentation in a diverse sample of individuals, 2) describe variation within and between these samples, and 3) demonstrate how quantitative measures can facilitate genotype-phenotype association tests.

MATERIALS AND METHODS

We quantitatively characterize skin, hair, and iris pigmentation using the Melanin (M) Index (skin) and CIELab values (hair) in 1,450 individuals who self-identify as African American, East Asian, European, Hispanic, or South Asian. We also quantify iris pigmentation in a subset of these individuals using CIELab values from high-resolution iris photographs. We compare mean skin M index and hair and iris CIELab values among populations using ANOVA and MANOVA respectively and test for genotype-phenotype associations in the European sample.

RESULTS

All five populations are significantly different for skin (P <2 × 10(-16) ) and hair color (P <2 × 10(-16) ). Our quantitative analysis of iris and hair pigmentation reinforces the continuous, rather than discrete, nature of these traits. We confirm the association of three loci (rs16891982, rs12203592, and rs12913832) with skin pigmentation and four loci (rs12913832, rs12203592, rs12896399, and rs16891982) with hair pigmentation. Interestingly, the derived rs12203592 T allele located within the IRF4 gene is associated with lighter skin but darker hair color.

DISCUSSION

The quantitative methods used here provide a fine-scale assessment of pigmentation phenotype and facilitate genotype-phenotype associations, even with relatively small sample sizes. This represents an important expansion of current investigations into pigmentation phenotype and associated genetic variation by including non-European and admixed populations. Am J Phys Anthropol 160:570-581, 2016. © 2015 Wiley Periodicals, Inc.

A Common Genetic Origin for Early Farmers from Mediterranean Cardial and Central European LBK Cultures

The spread of farming out of the Balkans and into the rest of Europe followed two distinct routes: An initial expansion represented by the Impressa and Cardial traditions, which followed the Northern Mediterranean coastline; and another expansion represented by the LBK (Linearbandkeramik) tradition, which followed the Danube River into Central Europe. Although genomic data now exist from samples representing the second migration, such data have yet to be successfully generated from the initial Mediterranean migration. To address this, we generated the complete genome of a 7,400-year-old Cardial individual (CB13) from Cova Bonica in Vallirana (Barcelona), as well as partial nuclear data from five others excavated from different sites in Spain and Portugal. CB13 clusters with all previously sequenced early European farmers and modern-day Sardinians. Furthermore, our analyses suggest that both Cardial and LBK peoples derived from a common ancient population located in or around the Balkan Peninsula. The Iberian Cardial genome also carries a discernible hunter–gatherer genetic signature that likely was not acquired by admixture with local Iberian foragers. Our results indicate that retrieving ancient genomes from similarly warm Mediterranean environments such as the Near East is technically feasible.

Skin pigmentation variation among populations of West Maharashtra, India

OBJECTIVES

Global patterns of skin pigmentation have evolved as an adaptation to local ultraviolet radiation (UVR). Indian populations exposed to intense UVR show great variation in skin pigmentation. The UVR-based selection model cannot satisfactorily address the high prevalence of light skin among these populations. Thus, the present study examines pigmentation variation among populations of West Maharashtra and the Indian subcontinent within the context of population structure and social hierarchy.

METHODS

Melanin index (MI) was measured from 555 individuals representing six endogamous populations of West Maharashtra. Skin pigmentation was assessed in terms of variation between populations and differences between and among castes and tribes. A linear regression analysis was run to assess the relationship among MI, UVR, and social hierarchy using published MI data from 13 Indian endogamous populations.

RESULTS

Skin pigmentation differed significantly among populations of West Maharashtra. Significant pigmentation variation exists between castes and tribes of West Maharashtra as well as across the Indian subcontinent. We observe a significant negative relationship between social hierarchy and skin pigmentation, whereas the relationship between UVR and MI is weak.

CONCLUSIONS

Our results suggest that various factors may have contributed to pigmentation diversity across the Indian subcontinent. The lack of correlation between UVR and MI suggests that natural selection may not have played a significant role in shaping pigmentation variation across the subcontinent. We discuss other possible explanations, including metabolic conservation and cultural factors such as traditional social hierarchies and strict endogamy that have led to the development of population structure.

Evidence for recent, population-specific evolution of the human mutation rate

As humans dispersed out of Africa they adapted to new environmental challenges, including changes in exposure to mutagenic solar radiation. Humans in temperate latitudes have acquired light skin that is relatively transparent to UV light, and some evidence suggests that their DNA damage response pathways have also experienced local adaptation. This raises the possibility that different populations have experienced different selective pressures affecting genome integrity. Here, I present evidence that the rate of a particular mutation type has recently increased in the European population, rising in frequency by 50% during the 40,000-80,000 y since Europeans began diverging from Asians. A comparison of SNPs private to Africa, Asia, and Europe in the 1000 Genomes data reveals that private European variation is enriched for the transition 5'-TCC-3' → 5'-TTC-3'. Although it is not clear whether UV played a causal role in changing the European mutational spectrum, 5'-TCC-3' → 5'-TTC-3' is known to be the most common somatic mutation present in melanoma skin cancers, as well as the mutation most frequently induced in vitro by UV. Regardless of its causality, this change indicates that DNA replication fidelity has not remained stable even since the origin of modern humans and might have changed numerous times during our recent evolutionary history.

A hidden Markov model for investigating recent positive selection through haplotype structure

Recent positive selection can increase the frequency of an advantageous mutant rapidly enough that a relatively long ancestral haplotype will be remained intact around it. We present a hidden Markov model (HMM) to identify such haplotype structures. With HMM identified haplotype structures, a population genetic model for the extent of ancestral haplotypes is then adopted for parameter inference of the selection intensity and the allele age. Simulations show that this method can detect selection under a wide range of conditions and has higher power than the existing frequency spectrum-based method. In addition, it provides good estimate of the selection coefficients and allele ages for strong selection. The method analyzes large data sets in a reasonable amount of running time. This method is applied to HapMap III data for a genome scan, and identifies a list of candidate regions putatively under recent positive selection. It is also applied to several genes known to be under recent positive selection, including the LCT, KITLG and TYRP1 genes in Northern Europeans, and OCA2 in East Asians, to estimate their allele ages and selection coefficients.

Novel genomic signals of recent selection in an Ethiopian population

The recent feasibility of genome-wide studies of adaptation in human populations has provided novel insights into biological pathways that have been affected by adaptive pressures. However, only a few African populations have been investigated using these genome-wide approaches. Here, we performed a genome-wide analysis for evidence of recent positive selection in a sample of 120 individuals of Wolaita ethnicity belonging to Omotic-speaking people who have inhabited the mid- and high-land areas of southern Ethiopia for millennia. Using the 11 HapMap populations as the comparison group, we found Wolaita-specific signals of recent positive selection in several human leukocyte antigen (HLA) loci. Notably, the selected loci overlapped with HLA regions that we previously reported to be associated with podoconiosis-a geochemical lymphedema of the lower legs common in the Wolaita area. We found selection signals in PPARA, a gene involved in energy metabolism during prolonged food deficiency. This finding is consistent with the dietary use of enset, a crop with high-carbohydrate and low-fat and -protein contents domesticated in Ethiopia subsequent to food deprivation 10 000 years ago, and with metabolic adaptation to high-altitude hypoxia. We observed novel selection signals in CDKAL1 and NEGR1, well-known diabetes and obesity susceptibility genes. Finally, the SLC24A5 gene locus known to be associated with skin pigmentation was in the top selection signals in the Wolaita, and the alleles of single-nucleotide polymorphisms rs1426654 and rs1834640 (SLC24A5) associated with light skin pigmentation in Eurasian populations were of high frequency (47.9%) in this Omotic-speaking indigenous Ethiopian population.

Genome flux and stasis in a five millennium transect of European prehistory

The Great Hungarian Plain was a crossroads of cultural transformations that have shaped European prehistory. Here we analyse a 5,000-year transect of human genomes, sampled from petrous bones giving consistently excellent endogenous DNA yields, from 13 Hungarian Neolithic, Copper, Bronze and Iron Age burials including two to high (~22 × ) and seven to ~1 × coverage, to investigate the impact of these on Europe’s genetic landscape. These data suggest genomic shifts with the advent of the Neolithic, Bronze and Iron Ages, with interleaved periods of genome stability. The earliest Neolithic context genome shows a European hunter-gatherer genetic signature and a restricted ancestral population size, suggesting direct contact between cultures after the arrival of the first farmers into Europe. The latest, Iron Age, sample reveals an eastern genomic influence concordant with introduced Steppe burial rites. We observe transition towards lighter pigmentation and surprisingly, no Neolithic presence of lactase persistence.

Recent advances in high-throughput sequencing techniques have enabled the analysis of ancient human genomes. Here the authors sequence ancient human genomes that span a period of 5,000 years, to understand the ancestral influence on Europe's genetic landscape.

Adaptations to local environments in modern human populations

After leaving sub-Saharan Africa around 50000-100000 years ago, anatomically modern humans have quickly occupied extremely diverse environments. Human populations were exposed to further environmental changes resulting from cultural innovations, such as the spread of farming, which gave rise to new selective pressures related to pathogen exposures and dietary shifts. In addition to changing the frequency of individual adaptive alleles, natural selection may also shape the overall genetic architecture of adaptive traits. Here, we review recent advances in understanding the genetic architecture of adaptive human phenotypes based on insights from the studies of lactase persistence, skin pigmentation and high-altitude adaptation. These adaptations evolved in parallel in multiple human populations, providing a chance to investigate independent realizations of the evolutionary process. We suggest that the outcome of adaptive evolution is often highly variable even under similar selective pressures. Finally, we highlight a growing need for detecting adaptations that did not follow the classical sweep model and for incorporating new sources of genetic evidence such as information from ancient DNA.

The interplay between natural selection and susceptibility to melanoma on allele 374F of SLC45A2 gene in a South European population

We aimed to study the selective pressures interacting on SLC45A2 to investigate the interplay between selection and susceptibility to disease. Thus, we enrolled 500 volunteers from a geographically limited population (Basques from the North of Spain) and by resequencing the whole coding region and intron 5 of the 34 most and the 34 least pigmented individuals according to the reflectance distribution, we observed that the polymorphism Leu374Phe (L374F, rs16891982) was statistically associated with skin color variability within this sample. In particular, allele 374F was significantly more frequent among the individuals with lighter skin. Further genotyping an independent set of 558 individuals of a geographically wider population with known ancestry in the Spanish population also revealed that the frequency of L374F was significantly correlated with the incident UV radiation intensity. Selection tests suggest that allele 374F is being positively selected in South Europeans, thus indicating that depigmentation is an adaptive process. Interestingly, by genotyping 119 melanoma samples, we show that this variant is also associated with an increased susceptibility to melanoma in our populations. The ultimate driving force for this adaptation is unknown, but it is compatible with the vitamin D hypothesis. This shows that molecular evolution analysis can be used as a useful technology to predict phenotypic and biomedical consequences in humans.

Cuba: exploring the history of admixture and the genetic basis of pigmentation using autosomal and uniparental markers

We carried out an admixture analysis of a sample comprising 1,019 individuals from all the provinces of Cuba. We used a panel of 128 autosomal Ancestry Informative Markers (AIMs) to estimate the admixture proportions. We also characterized a number of haplogroup diagnostic markers in the mtDNA and Y-chromosome in order to evaluate admixture using uniparental markers. Finally, we analyzed the association of 16 single nucleotide polymorphisms (SNPs) with quantitative estimates of skin pigmentation. In the total sample, the average European, African and Native American contributions as estimated from autosomal AIMs were 72%, 20% and 8%, respectively. The Eastern provinces of Cuba showed relatively higher African and Native American contributions than the Western provinces. In particular, the highest proportion of African ancestry was observed in the provinces of Guantánamo (40%) and Santiago de Cuba (39%), and the highest proportion of Native American ancestry in Granma (15%), Holguín (12%) and Las Tunas (12%). We found evidence of substantial population stratification in the current Cuban population, emphasizing the need to control for the effects of population stratification in association studies including individuals from Cuba. The results of the analyses of uniparental markers were concordant with those observed in the autosomes. These geographic patterns in admixture proportions are fully consistent with historical and archaeological information. Additionally, we identified a sex-biased pattern in the process of gene flow, with a substantially higher European contribution from the paternal side, and higher Native American and African contributions from the maternal side. This sex-biased contribution was particularly evident for Native American ancestry. Finally, we observed that SNPs located in the genes SLC24A5 and SLC45A2 are strongly associated with melanin levels in the sample.

Geographical ancestry is a key determinant of epidermal morphology and dermal composition

BACKGROUND

Geographical ancestry plays a key role in determining the susceptibility of human skin to external insults and dermatological disease. Despite this, studies of skin from individuals of diverse geographical ancestry focus primarily on epidermal pigmentation. Few reports characterize the gross morphology and composition of the dermis and dermal-epidermal junction (DEJ).

OBJECTIVES

To characterize epidermal morphology and dermal composition in skin from individuals of diverse geographical ancestry.

METHODS

Immunohistochemical techniques were used to assess epidermal morphology and protein composition of the DEJ and dermal extracellular matrix in photoprotected skin from young African, Eurasian and Far East Asian individuals (n = 7 per group; age 18-30 years).

RESULTS

The epidermis of African skin was thicker, with deeper rete ridges and a more convoluted DEJ than Eurasian and Far East Asian skin. Compared with Eurasians, protein composition of the DEJ was collagen VII poor in African and Far East Asian skin (P < 0·001 and P < 0·01, respectively); the dermis of African skin was enriched in fibrillar collagens (P < 0·05), but was relatively elastin poor (P < 0·05). African dermis was abundant in fibrillin-rich microfibrils and fibulin-5 (P < 0·001 and P < 0·001, respectively) compared with Eurasian and Far East Asian skin.

CONCLUSIONS

We demonstrate that fundamental differences exist in skin structure and composition in individuals of diverse geographical ancestry. Disparate environmental pressures encountered by ancestral human populations living at different latitudes may have driven adaptations in skin structure and composition. Further research into the functional significance and clinical consequences of these differences is warranted.

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.

Sun exposure, sunbeds and sunscreens and melanoma. What are the controversies?

The association between various measures of sun exposure and melanoma risk is quite complex to dissect as many case-control studies of melanoma included different subtypes of melanomas which are likely to be biologically different, so interpretation of the data is difficult. Screening bias in countries with high levels of sun exposure is also an issue. Now that progress is being made in the genetic subclassification of melanoma tumours, it is apparent that melanomas have different somatic changes according to body sites/histological subtypes and that UV exposure may be relevant for some but not all types of melanomas. Melanoma behaviour also points to non-sun-related risk factors, and complex gene-environment interactions are likely. As UV exposure is the only environmental factor ever linked to melanoma, it is still prudent to avoid excessive sun exposure and sunburn especially in poor tanners. However, the impact of strict sun avoidance, which should not be recommended, may take years to be apparent as vitamin D deficiency is a now a common health issue in Caucasian populations, with a significant impact on health in general.

Direct evidence for positive selection of skin, hair, and eye pigmentation in Europeans during the last 5,000 y

Pigmentation is a polygenic trait encompassing some of the most visible phenotypic variation observed in humans. Here we present direct estimates of selection acting on functional alleles in three key genes known to be involved in human pigmentation pathways--HERC2, SLC45A2, and TYR--using allele frequency estimates from Eneolithic, Bronze Age, and modern Eastern European samples and forward simulations. Neutrality was overwhelmingly rejected for all alleles studied, with point estimates of selection ranging from around 2-10% per generation. Our results provide direct evidence that strong selection favoring lighter skin, hair, and eye pigmentation has been operating in European populations over the last 5,000 y.

Information compression exploits patterns of genome composition to discriminate populations and highlight regions of evolutionary interest

Background

Genomic information allows population relatedness to be inferred and selected genes to be identified. Single nucleotide polymorphism microarray (SNP-chip) data, a proxy for genome composition, contains patterns in allele order and proportion. These patterns can be quantified by compression efficiency (CE). In principle, the composition of an entire genome can be represented by a CE number quantifying allele representation and order.

Results

We applied a compression algorithm (DEFLATE) to genome-wide high-density SNP data from 4,155 human, 1,800 cattle, 1,222 sheep, 81 dogs and 49 mice samples. All human ethnic groups can be clustered by CE and the clusters recover phylogeography based on traditional fixation index (F_ST) analyses. CE analysis of other mammals results in segregation by breed or species, and is sensitive to admixture and past effective population size. This clustering is a consequence of individual patterns such as runs of homozygosity. Intriguingly, a related approach can also be used to identify genomic loci that show population-specific CE segregation. A high resolution CE ‘sliding window’ scan across the human genome, organised at the population level, revealed genes known to be under evolutionary pressure. These include SLC24A5 (European and Gujarati Indian skin pigmentation), HERC2 (European eye color), LCT (European and Maasai milk digestion) and EDAR (Asian hair thickness). We also identified a set of previously unidentified loci with high population-specific CE scores including the chromatin remodeler SCMH1 in Africans and EDA2R in Asians. Closer inspection reveals that these prioritised genomic regions do not correspond to simple runs of homozygosity but rather compositionally complex regions that are shared by many individuals of a given population. Unlike F_ST, CE analyses do not require ab initio population comparisons and are amenable to the hemizygous X chromosome.

Conclusions

We conclude with a discussion of the implications of CE for a complex systems science view of genome evolution. CE allows one to clearly visualise the evolution of individual genomes and populations through a formal, mathematically-rigorous information space. Overall, CE makes a set of biological predictions, some of which are unique and await functional validation.

Extreme population differences in the human zinc transporter ZIP4 (SLC39A4) are explained by positive selection in Sub-Saharan Africa

Extreme differences in allele frequency between West Africans and Eurasians were observed for a leucine-to-valine substitution (Leu372Val) in the human intestinal zinc uptake transporter, ZIP4, yet no further evidence was found for a selective sweep around the ZIP4 gene (SLC39A4). By interrogating allele frequencies in more than 100 diverse human populations and resequencing Neanderthal DNA, we confirmed the ancestral state of this locus and found a strong geographical gradient for the derived allele (Val372), with near fixation in West Africa. In extensive coalescent simulations, we show that the extreme differences in allele frequency, yet absence of a classical sweep signature, can be explained by the effect of a local recombination hotspot, together with directional selection favoring the Val372 allele in Sub-Saharan Africans. The possible functional effect of the Leu372Val substitution, together with two pathological mutations at the same codon (Leu372Pro and Leu372Arg) that cause acrodermatitis enteropathica (a disease phenotype characterized by extreme zinc deficiency), was investigated by transient overexpression of human ZIP4 protein in HeLa cells. Both acrodermatitis mutations cause absence of the ZIP4 transporter cell surface expression and nearly absent zinc uptake, while the Val372 variant displayed significantly reduced surface protein expression, reduced basal levels of intracellular zinc, and reduced zinc uptake in comparison with the Leu372 variant. We speculate that reduced zinc uptake by the ZIP4-derived Val372 isoform may act by starving certain pathogens of zinc, and hence may have been advantageous in Sub-Saharan Africa. Moreover, these functional results may indicate differences in zinc homeostasis among modern human populations with possible relevance for disease risk.

Association of genetic variants with self-assessed color categories in Brazilians

The Brazilian population was formed by extensive admixture of three different ancestral roots: Amerindians, Europeans and Africans. Our previous work has shown that at an individual level, ancestry, as estimated using molecular markers, was a poor predictor of color in Brazilians. We now investigate if SNPs known to be associated with human skin pigmentation can be used to predict color in Brazilians. For that, we studied the association of fifteen SNPs, previously known to be linked with skin color, in 243 unrelated Brazilian individuals self-identified as White, Browns or Blacks from Rio de Janeiro and 212 unrelated Brazilian individuals self-identified as White or Blacks from São Paulo. The significance of association of SNP genotypes with self-assessed color was evaluated using partial regression analysis. After controlling for ancestry estimates as covariates, only four SNPs remained significantly associated with skin pigmentation: rs1426654 and rs2555364 within SLC24A5, rs16891982 at SLC45A2 and rs1042602 at TYR. These loci are known to be involved in melanin synthesis or transport of melanosomes. We found that neither genotypes of these SNPs, nor their combination with biogeographical ancestry in principal component analysis, could predict self-assessed color in Brazilians at an individual level. However, significant correlations did emerge at group level, demonstrating that even though elements other than skin, eye and hair pigmentation do influence self-assessed color in Brazilians, the sociological act of self-classification is still substantially dependent of genotype at these four SNPs.

Vitamin D: calcium and bone homeostasis during evolution

Vitamin D3 is already found early in the evolution of life but essentially as inactive end products of the photochemical reaction of 7-dehydrocholestol with ultraviolet light B. A full vitamin D (refers to vitamin D2 and D3) endocrine system, characterized by a specific VDR (vitamin D receptor, member of the nuclear receptor family), specific vitamin D metabolizing CYP450 enzymes regulated by calciotropic hormones and a dedicated plasma transport-protein is only found in vertebrates. In the earliest vertebrates (lamprey), vitamin D metabolism and VDR may well have originated from a duplication of a common PRX/VDR ancestor gene as part of a xenobiotic detoxification pathway. The vitamin D endocrine system, however, subsequently became an important regulator of calcium supply for an extensive calcified skeleton. Vitamin D is essential for normal calcium and bone homeostasis as shown by rickets in vitamin D-deficient growing amphibians, reptiles, birds and mammals. From amphibians onward, bone is gradually more dynamic with regulated bone resorption, mainly by combined action of PTH and 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3) on the generation and function of multinucleated osteoclasts. Therefore, bone functions as a large internal calcium reservoir, under the control of osteoclasts. Osteocytes also display a remarkable spectrum of activities, including mechanical sensing and regulating mineral homeostasis, but also have an important role in global nutritional and energy homeostasis. Mineralization from reptiles onward is under the control of well-regulated SIBLING proteins and associated enzymes, nearly all under the control of 1,25(OH)2D3. The vitamin D story thus started as inert molecule but gained an essential role for calcium and bone homeostasis in terrestrial animals to cope with the challenge of higher gravity and calcium-poor environment.

Low frequency variants, collapsed based on biological knowledge, uncover complexity of population stratification in 1000 genomes project data

Analyses investigating low frequency variants have the potential for explaining additional genetic heritability of many complex human traits. However, the natural frequencies of rare variation between human populations strongly confound genetic analyses. We have applied a novel collapsing method to identify biological features with low frequency variant burden differences in thirteen populations sequenced by the 1000 Genomes Project. Our flexible collapsing tool utilizes expert biological knowledge from multiple publicly available database sources to direct feature selection. Variants were collapsed according to genetically driven features, such as evolutionary conserved regions, regulatory regions genes, and pathways. We have conducted an extensive comparison of low frequency variant burden differences (MAF<0.03) between populations from 1000 Genomes Project Phase I data. We found that on average 26.87% of gene bins, 35.47% of intergenic bins, 42.85% of pathway bins, 14.86% of ORegAnno regulatory bins, and 5.97% of evolutionary conserved regions show statistically significant differences in low frequency variant burden across populations from the 1000 Genomes Project. The proportion of bins with significant differences in low frequency burden depends on the ancestral similarity of the two populations compared and types of features tested. Even closely related populations had notable differences in low frequency burden, but fewer differences than populations from different continents. Furthermore, conserved or functionally relevant regions had fewer significant differences in low frequency burden than regions under less evolutionary constraint. This degree of low frequency variant differentiation across diverse populations and feature elements highlights the critical importance of considering population stratification in the new era of DNA sequencing and low frequency variant genomic analyses.

Low frequency variants are likely to play an important role in uncovering complex trait heritability; however, they are often continent or population specific. This specificity complicates genetic analyses investigating low frequency variants for two reasons: low frequency variant signals in an association test are often difficult to generalize beyond a single population or continental group, and there is an increase in false positive results in association analyses due to underlying population stratification. In order to reveal the magnitude of low frequency population stratification, we performed pairwise population comparisons using the 1000 Genomes Project Phase I data to investigate differences in low frequency variant burden across multiple biological features. We found that low frequency variant confounding is much more prevalent than one might expect, even within continental groups. The proportion of significant differences in low frequency variant burden was also dependent on the region of interest; for example, annotated regulatory regions showed fewer low frequency burden differences between populations than intergenic regions. Knowledge of population structure and the genomic landscape in a region of interest are important factors in determining the extent of confounding due to population stratification in a low frequency genomic analysis.

Population Genomics of Human Adaptation

Recent advances in genotyping technologies have facilitated genome-wide scans for natural selection. Identification of targets of natural selection will shed light on processes of human adaptation and evolution and could be important for identifying variation that influences both normal human phenotypic variation as well as disease susceptibility. Here we focus on studies of natural selection in modern humans who originated ~200,000 years go in Africa and migrated across the globe ~50,000 - 100,000 years ago. Movement into new environments, as well as changes in culture and technology including plant and animal domestication, resulted in local adaptation to diverse environments. We summarize statistical approaches for detecting targets of natural selection and for distinguishing the effects of demographic history from natural selection. On a genome-wide scale, immune-related genes appear to be major targets of positive selection. Genes associated with reproduction and fertility also appear to be fast evolving. Additional examples of recent human adaptation include genes associated with lactase persistence, eccrine glands, and response to hypoxia. Lastly, we emphasize the need to supplement scans of selection with functional studies to demonstrate the physiologic impact of candidate loci.

The light skin allele of SLC24A5 in South Asians and Europeans shares identity by descent

Skin pigmentation is one of the most variable phenotypic traits in humans. A non-synonymous substitution (rs1426654) in the third exon of SLC24A5 accounts for lighter skin in Europeans but not in East Asians. A previous genome-wide association study carried out in a heterogeneous sample of UK immigrants of South Asian descent suggested that this gene also contributes significantly to skin pigmentation variation among South Asians. In the present study, we have quantitatively assessed skin pigmentation for a largely homogeneous cohort of 1228 individuals from the Southern region of the Indian subcontinent. Our data confirm significant association of rs1426654 SNP with skin pigmentation, explaining about 27% of total phenotypic variation in the cohort studied. Our extensive survey of the polymorphism in 1573 individuals from 54 ethnic populations across the Indian subcontinent reveals wide presence of the derived-A allele, although the frequencies vary substantially among populations. We also show that the geospatial pattern of this allele is complex, but most importantly, reflects strong influence of language, geography and demographic history of the populations. Sequencing 11.74 kb of SLC24A5 in 95 individuals worldwide reveals that the rs1426654-A alleles in South Asian and West Eurasian populations are monophyletic and occur on the background of a common haplotype that is characterized by low genetic diversity. We date the coalescence of the light skin associated allele at 22–28 KYA. Both our sequence and genome-wide genotype data confirm that this gene has been a target for positive selection among Europeans. However, the latter also shows additional evidence of selection in populations of the Middle East, Central Asia, Pakistan and North India but not in South India.

Human skin color is one of the most visible aspects of human diversity. The genetic basis of pigmentation in Europeans has been understood to some extent, but our knowledge about South Asians has been restricted to a handful of studies. It has been suggested that a single nucleotide difference in SLC24A5 accounts for 25–38% European-African pigmentation differences and correlates with lighter skin. This genetic variant has also been associated with skin color variation among South Asians living in the UK. Here, we report a study based on a homogenous cohort of South India. Our results confirm that SLC24A5 plays a key role in pigmentation diversity of South Asians. Country-wide screening of the variant reveals that the light skin associated allele is widespread in the Indian subcontinent and its complex patterning is shaped by a combination of processes involving selection and demographic history of the populations. By studying the variation of SLC24A5 sequences among a diverse set of individuals, we show that the light skin associated allele in South Asians is identical by descent to that found in Europeans. Our study also provides new insights into positive selection acting on the gene and the evolutionary history of light skin in humans.

Molecular phylogeography of a human autosomal skin color locus under natural selection

Divergent natural selection caused by differences in solar exposure has resulted in distinctive variations in skin color between human populations. The derived light skin color allele of the SLC24A5 gene, A111T, predominates in populations of Western Eurasian ancestry. To gain insight into when and where this mutation arose, we defined common haplotypes in the genomic region around SLC24A5 across diverse human populations and deduced phylogenetic relationships between them. Virtually all chromosomes carrying the A111T allele share a single 78-kb haplotype that we call C11, indicating that all instances of this mutation in human populations share a common origin. The C11 haplotype was most likely created by a crossover between two haplotypes, followed by the A111T mutation. The two parental precursor haplotypes are found from East Asia to the Americas but are nearly absent in Africa. The distributions of C11 and its parental haplotypes make it most likely that these two last steps occurred between the Middle East and the Indian subcontinent, with the A111T mutation occurring after the split between the ancestors of Europeans and East Asians.

Global patterns of diversity and selection in human tyrosinase gene

Global variation in skin pigmentation is one of the most striking examples of environmental adaptation in humans. More than two hundred loci have been identified as candidate genes in model organisms and a few tens of these have been found to be significantly associated with human skin pigmentation in genome-wide association studies. However, the evolutionary history of different pigmentation genes is rather complex: some loci have been subjected to strong positive selection, while others evolved under the relaxation of functional constraints in low UV environment. Here we report the results of a global study of the human tyrosinase gene, which is one of the key enzymes in melanin production, to assess the role of its variation in the evolution of skin pigmentation differences among human populations. We observe a higher rate of non-synonymous polymorphisms in the European sample consistent with the relaxation of selective constraints. A similar pattern was previously observed in the MC1R gene and concurs with UV radiation-driven model of skin color evolution by which mutations leading to lower melanin levels and decreased photoprotection are subject to purifying selection at low latitudes while being tolerated or even favored at higher latitudes because they facilitate UV-dependent vitamin D production. Our coalescent date estimates suggest that the non-synonymous variants, which are frequent in Europe and North Africa, are recent and have emerged after the separation of East and West Eurasian populations.

Exploring signatures of positive selection in pigmentation candidate genes in populations of East Asian ancestry

Background

Currently, there is very limited knowledge about the genes involved in normal pigmentation variation in East Asian populations. We carried out a genome-wide scan of signatures of positive selection using the 1000 Genomes Phase I dataset, in order to identify pigmentation genes showing putative signatures of selective sweeps in East Asia. We applied a broad range of methods to detect signatures of selection including: 1) Tests designed to identify deviations of the Site Frequency Spectrum (SFS) from neutral expectations (Tajima’s D, Fay and Wu’s H and Fu and Li’s D* and F*), 2) Tests focused on the identification of high-frequency haplotypes with extended linkage disequilibrium (iHS and Rsb) and 3) Tests based on genetic differentiation between populations (LSBL). Based on the results obtained from a genome wide analysis of 25 kb windows, we constructed an empirical distribution for each statistic across all windows, and identified pigmentation genes that are outliers in the distribution.

Results

Our tests identified twenty genes that are relevant for pigmentation biology. Of these, eight genes (ATRN, EDAR, KLHL7, MITF, OCA2, TH, TMEM33 and TRPM1,) were extreme outliers (top 0.1% of the empirical distribution) for at least one statistic, and twelve genes (ADAM17, BNC2, CTSD, DCT, EGFR, LYST, MC1R, MLPH, OPRM1, PDIA6, PMEL (SILV) and TYRP1) were in the top 1% of the empirical distribution for at least one statistic. Additionally, eight of these genes (BNC2, EGFR, LYST, MC1R, OCA2, OPRM1, PMEL (SILV) and TYRP1) have been associated with pigmentary traits in association studies.

Conclusions

We identified a number of putative pigmentation genes showing extremely unusual patterns of genetic variation in East Asia. Most of these genes are outliers for different tests and/or different populations, and have already been described in previous scans for positive selection, providing strong support to the hypothesis that recent selective sweeps left a signature in these regions. However, it will be necessary to carry out association and functional studies to demonstrate the implication of these genes in normal pigmentation variation.

Colorful DNA polymorphisms in humans

In this review article we summarize current knowledge on how variation on the DNA level influences human pigmentation including color variation of iris, hair, and skin. We review recent progress in the field of human pigmentation genetics by focusing on the genes and DNA polymorphisms discovered to be involved in determining human pigmentation traits, their association with diseases particularly skin cancers, and their power to predict human eye, hair, and skin colors with potential utilization in forensic investigations.

Human pigmentation genes under environmental selection

Genome-wide association studies and comparative genomics have established major loci and specific polymorphisms affecting human skin, hair and eye color. Environmental changes have had an impact on selected pigmentation genes as populations have expanded into different regions of the globe.