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

Whole genome resequencing reveals genomic regions related to red plumage in ducks

Plumage color is a characteristic trait of ducks that originates as a result of natural and artificial selection. As a conspicuous phenotypic feature, it is a breed characteristic. Previous studies have identified some genes associated with the formation of black and white plumage in ducks. However, studies on the genetic basis underlying the red plumage phenotype in ducks are limited. Here, genome-wide association analysis (GWAS) and selection signal detection (Fst, θπ ratio, and cross-population composite likelihood ratio [XP-CLR]) were conducted to identify candidate regions and genes underlying duck plumage color phenotype. Selection signal detection revealed 29 overlapping genes (including ENPP1 and ULK1) significantly associated with red plumage color in Ji'an Red ducks. ENSAPLG00000012679, ESRRG, and SPATA5 were identified as candidate genes associated with red plumage using GWAS. Selection signal detection revealed that 19 overlapping genes (including GMDS, PDIA6, and ODC1) significantly correlated with light brown plumage in Brown Tsaiya ducks. GWAS to narrow down the significant regions further revealed nine candidate genes (AKT1, ATP6V1C2, GMDS, LRP4, MAML3, PDIA6, PLD5, TMEM63B, and TSPAN8). Notably, in Brown Tsaiya ducks, GMDS, ODC1, and PDIA6 exhibit significantly differentiated allele frequencies among other feather-colored ducks, while in Ji'an Red ducks, ENSAPLG00000012679 has different allele frequency distributions compared with that in other feather-colored ducks. This study offers new insights into the variation and selection of the red plumage phenotype using GWAS and selective signals.

Uncovering genes underlying coat color variation in indigenous cattle breeds through genome-wide positive selection

The identification of candidate genes related to pigmentation and under selective sweep provides insights into the genetic basis of pigmentation and the evolutionary forces that have shaped this variation. The selective sweep events in the genes responsible for normal coat color in Indian cattle groups are still unknown. To find coat color genes displaying signs of selective sweeps in the indigenous cattle, we compiled a list of candidate genes previously investigated for their association with coat color and pigmentation. After that, we performed a genome-wide scan of positive selection signatures using the BovineSNP50K Bead Chip in 187 individuals of seven indigenous breeds. We applied a wide range of methods to find evidence of selection, such as Tajima's D, CLR, iHS, varLD, ROH, and FST. We found a total of sixteen genes under selective sweep, that were involved in coat color and pigmentation physiology. These genes are CRIM1 in Gir, MC1R in Sahiwal, MYO5A, PMEL and POMC in Tharparkar, TYRP1, ERBB2, and ASIP in Red Sindhi, MITF, LOC789175, PAX3 and TYR in Ongole, and IRF2, SDR165 and, KIT in Nelore, ADAMTS19 in Hariana. These genes are related to melanin synthesis, the biology of melanocytes and melanosomes, and the migration and survival of melanocytes during development.

Epigenetics insights from perceived facial aging

Facial aging is the most visible manifestation of aging. People desire to look younger than others of the same chronological age. Hence, perceived age is often used as a visible marker of aging, while biological age, often estimated by methylation markers, is used as an objective measure of age. Multiple epigenetics-based clocks have been developed for accurate estimation of general biological age and the age of specific organs, including the skin. However, it is not clear whether the epigenetic biomarkers (CpGs) used in these clocks are drivers of aging processes or consequences of aging. In this proof-of-concept study, we integrate data from GWAS on perceived facial aging and EWAS on CpGs measured in blood. By running EW Mendelian randomization, we identify hundreds of putative CpGs that are potentially causal to perceived facial aging with similar numbers of damaging markers that causally drive or accelerate facial aging and protective methylation markers that causally slow down or protect from aging. We further demonstrate that while candidate causal CpGs have little overlap with known epigenetics-based clocks, they affect genes or proteins with known functions in skin aging, such as skin pigmentation, elastin, and collagen levels. Overall, our results suggest that blood methylation markers reflect facial aging processes, and thus can be used to quantify skin aging and develop anti-aging solutions that target the root causes of aging.

Tensor decomposition based feature extraction and classification to detect natural selection from genomic data

Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under non-convex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data while preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx , which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.

MaLAdapt Reveals Novel Targets of Adaptive Introgression From Neanderthals and Denisovans in Worldwide Human Populations

Adaptive introgression (AI) facilitates local adaptation in a wide range of species. Many state-of-the-art methods detect AI with ad-hoc approaches that identify summary statistic outliers or intersect scans for positive selection with scans for introgressed genomic regions. Although widely used, approaches intersecting outliers are vulnerable to a high false-negative rate as the power of different methods varies, especially for complex introgression events. Moreover, population genetic processes unrelated to AI, such as background selection or heterosis, may create similar genomic signals to AI, compromising the reliability of methods that rely on neutral null distributions. In recent years, machine learning (ML) methods have been increasingly applied to population genetic questions. Here, we present a ML-based method called MaLAdapt for identifying AI loci from genome-wide sequencing data. Using an Extra-Trees Classifier algorithm, our method combines information from a large number of biologically meaningful summary statistics to capture a powerful composite signature of AI across the genome. In contrast to existing methods, MaLAdapt is especially well-powered to detect AI with mild beneficial effects, including selection on standing archaic variation, and is robust to non-AI selective sweeps, heterosis from deleterious mutations, and demographic misspecification. Furthermore, MaLAdapt outperforms existing methods for detecting AI based on the analysis of simulated data and the validation of empirical signals through visual inspection of haplotype patterns. We apply MaLAdapt to the 1000 Genomes Project human genomic data and discover novel AI candidate regions in non-African populations, including genes that are enriched in functionally important biological pathways regulating metabolism and immune responses.

Genomic Uniqueness of Local Sheep Breeds From Morocco

Sheep farming is a major source of meat in Morocco and plays a key role in the country's agriculture. This study aims at characterizing the whole-genome diversity and demographic history of the main Moroccan sheep breeds, as well as to identify selection signatures within and between breeds. Whole genome data from 87 individuals representing the five predominant local breeds were used to estimate their level of neutral genetic diversity and to infer the variation of their effective population size over time. In addition, we used two methods to detect selection signatures: either for detecting selective sweeps within each breed separately or by detecting differentially selected regions by contrasting different breeds. We identified hundreds of genomic regions putatively under selection, which related to several biological terms involved in local adaptation or the expression of zootechnical performances such as Growth, UV protection, Cell maturation or Feeding behavior. The results of this study revealed selection signatures in genes that have an important role in traits of interest and increased our understanding of how genetic diversity is distributed in these local breeds. Thus, Moroccan local sheep breeds exhibit both a high genetic diversity and a large set of adaptive variations, and therefore, represent a valuable genetic resource for the conservation of sheep in the context of climate change.

A large Canadian cohort provides insights into the genetic architecture of human hair colour

Hair colour is a polygenic phenotype that results from differences in the amount and ratio of melanins located in the hair bulb. Genome-wide association studies (GWAS) have identified many loci involved in the pigmentation pathway affecting hair colour. However, most of the associated loci overlap non-protein coding regions and many of the molecular mechanisms underlying pigmentation variation are still not understood. Here, we conduct GWAS meta-analyses of hair colour in a Canadian cohort of 12,741 individuals of European ancestry. By performing fine-mapping analyses we identify candidate causal variants in pigmentation loci associated with blonde, red and brown hair colour. Additionally, we observe colocalization of several GWAS hits with expression and methylation quantitative trait loci (QTLs) of cultured melanocytes. Finally, transcriptome-wide association studies (TWAS) further nominate the expression of EDNRB and CDK10 as significantly associated with hair colour. Our results provide insights on the mechanisms regulating pigmentation biology in humans.

GWAS Identifies Multiple Genetic Loci for Skin Color in Korean Women

Human skin color is largely determined by genetic factors. Recent GWASs have reported several genetic variants associated with skin color, mostly in European and African populations. In this study, we performed GWAS in 17,019 Korean women to identify genetic variants associated with facial skin color, quantitatively measured as CIELAB color index. We identified variants in three, one, and six genomic loci associated with facial skin color index L∗, a∗, and b∗ values, respectively, and replicated the associations (combined analysis P-value < 5.0 × 10^-8). The significant loci included variants in known genes (OCA2 rs74653330, BNC2 rs16935073, rs72620727 near KITLG, and SLC6A17 rs6689641) and to our knowledge previously unreported genes (SCARB1 rs10846744, SYN2 rs12629034, and LINC00486 rs6543678). This is GWAS to elucidate genetic variants of facial skin color in a Korean female population. Further functional characterizations of the investigated genes are warranted to elucidate their contribution to skin pigmentation-related traits.

Influence of Ethnicities and Skin Color Variations in Different Populations: A Review

BACKGROUND

In the world scientific tradition, skin color is the primary physical characteristic used to divide humans into groups. Human skin has a wide range of tones and colors, which can be seen in a wide range of demographic populations. Many factors influence the color of people's skin, but the pigment melanin is by far the most important. Melanin is produced by cells called melanocytes in the skin and is the primary determinant of skin color in people with darker skin. Indeed, >150 genes have now been identified as having a direct or indirect effect on skin color. Vitamin D has recently been discovered to regulate cellular proliferation and differentiation in a variety of tissues, including the skin. The mechanisms through which the active vitamin D metabolite 1,25 dihydroxyvitamin D3 (or calcitriol) affects keratinocyte development are numerous and overlap with the mechanisms by which calcium influences keratinocyte differentiation. Ultraviolet (UV) is the most major modifiable risk factor for skin cancer and many other environmental-influenced skin disorders when it is abundant in the environment. Although the UV component of sunlight is known to cause skin damage, few researches have looked at the impact of non-UV solar radiation on skin physiology in terms of inflammation, and there is less information on the role of visible light in pigmentation.

SUMMARY

The quantity and quality of melanin are regulating by the expression of genes. The enzyme tyrosinase is primarily responsible for the genetic mechanism that controls human skin color. Genetics determines constitutive skin color, which is reinforced by facultative melanogenesis and tanning reactions. High quantities of melanin and melanogenic substances are typically accepted in darker skin to protect against UV radiation-induced molecular damage. Previous research has proposed that skin color variation is caused by a dynamic genetic mechanism, contributing to our understanding of how population demographic history and natural selection shape human genetic and phenotypic diversity. However, the most significant ethnic skin color difference is determined by melanin content. This current review aimed to assess the influence of skin color variations in skin structure and functions as well as difference in dermatological disease patterns. Also, this article reviewed several cases of skin color adaptation in different populations. Key Messages: Skin color impacts the composition and activity. Therefore, the contrast of dermatological ailments between distinct race-related categories is remarkable. Skin color adaptation is a challenging procedure. Refinement of skin color is an age-old craving of humans with ever-evolving drifts.

Quantitative Human Paleogenetics: What can Ancient DNA Tell us About Complex Trait Evolution?

Genetic association data from national biobanks and large-scale association studies have provided new prospects for understanding the genetic evolution of complex traits and diseases in humans. In turn, genomes from ancient human archaeological remains are now easier than ever to obtain, and provide a direct window into changes in frequencies of trait-associated alleles in the past. This has generated a new wave of studies aiming to analyse the genetic component of traits in historic and prehistoric times using ancient DNA, and to determine whether any such traits were subject to natural selection. In humans, however, issues about the portability and robustness of complex trait inference across different populations are particularly concerning when predictions are extended to individuals that died thousands of years ago, and for which little, if any, phenotypic validation is possible. In this review, we discuss the advantages of incorporating ancient genomes into studies of trait-associated variants, the need for models that can better accommodate ancient genomes into quantitative genetic frameworks, and the existing limits to inferences about complex trait evolution, particularly with respect to past populations.

Evolutionary genetics of skin pigmentation in African populations

Skin color is a highly heritable human trait, and global variation in skin pigmentation has been shaped by natural selection, migration and admixture. Ethnically diverse African populations harbor extremely high levels of genetic and phenotypic diversity, and skin pigmentation varies widely across Africa. Recent genome-wide genetic studies of skin pigmentation in African populations have advanced our understanding of pigmentation biology and human evolutionary history. For example, novel roles in skin pigmentation for loci near MFSD12 and DDB1 have recently been identified in African populations. However, due to an underrepresentation of Africans in human genetic studies, there is still much to learn about the evolutionary genetics of skin pigmentation. Here, we summarize recent progress in skin pigmentation genetics in Africans and discuss the importance of including more ethnically diverse African populations in future genetic studies. In addition, we discuss methods for functional validation of adaptive variants related to skin pigmentation.

The genetic variation of lactase persistence alleles in Sudan and South Sudan

Lactase persistence (LP) is a well-studied example of a Mendelian trait under selection in some human groups due to gene-culture coevolution. We investigated the frequencies of genetic variants linked to LP in Sudanese and South Sudanese populations. These populations have diverse subsistence patterns, and some are dependent on milk to various extents, not only from cows but also from other livestock such as camels and goats. We sequenced a 316-bp region involved in regulating the expression of the LCT gene on chromosome 2, which encompasses five polymorphisms that have been associated with LP. Pastoralist populations showed a higher frequency of LP-associated alleles compared with nonpastoralist groups, hinting at positive selection also among northeast African pastoralists. Among the LP variants, the -14009:G variant occurs at the highest frequency among the investigated populations, followed by the -13915:G variant, which is likely of Middle Eastern origin, consistent with Middle Eastern gene flow to the Sudanese populations. There was no incidence of the “East African” LP allele (-14010:C) in the Sudanese and South Sudanese groups, and only one heterozygous individual for the “European” LP allele (-13910:T), suggesting limited recent admixture from these geographic regions. The Beja population of the Beni Amer show three different LP variants at substantial and similar levels, resulting in one of the greatest aggregation of LP variants among all populations across the world.

Population Genomics Reveals Incipient Speciation, Introgression, and Adaptation in the African Mona Monkey (Cercopithecus mona)

Guenons (tribe Cercopithecini) are the most widely distributed nonhuman primate in the tropical forest belt of Africa and show considerable phenotypic, taxonomic, and ecological diversity. However, genomic information for most species within this group is still lacking. Here, we present a high-quality de novo genome (total 2.90 Gb, contig N50 equal to 22.7 Mb) of the mona monkey (Cercopithecus mona), together with genome resequencing data of 13 individuals sampled across Nigeria. Our results showed differentiation between populations from East and West of the Niger River ∼84 ka and potential ancient introgression in the East population from other mona group species. The PTPRK, FRAS1, BNC2, and EDN3 genes related to pigmentation displayed signals of introgression in the East population. Genomic scans suggest that immunity genes such as AKT3 and IL13 (possibly involved in simian immunodeficiency virus defense), and G6PD, a gene involved in malaria resistance, are under positive natural selection. Our study gives insights into differentiation, natural selection, and introgression in guenons.

Dissecting dynamics and differences of selective pressures in the evolution of human pigmentation

Human pigmentation is a highly diverse and complex trait among populations and has drawn particular attention from both academic and non-academic investigators for thousands of years. Previous studies detected selection signals in several human pigmentation genes, but few studies have integrated contribution from multiple genes to the evolution of human pigmentation. Moreover, none has quantified selective pressures on human pigmentation over epochs and between populations. Here, we dissect dynamics and differences of selective pressures during different periods and between distinct populations with new approaches. We use genotype data of 19 genes associated with human pigmentation from 17 publicly available datasets and obtain data for 2346 individuals of six representative population groups from across the world. Our results quantify the strength of natural selection on light pigmentation not only in modern Europeans (0.0259/generation) but also in proto-Eurasians (0.00650/generation). Our results also suggest that several derived alleles associated with human dark pigmentation may be under positive directional selection in some African populations. Our study provides the first attempt to quantitatively investigate the dynamics of selective pressures during different time periods in the evolution of human pigmentation.This article has an associated First Person interview with the first author of the article.

Probabilistic Estimation of Identity by Descent Segment Endpoints and Detection of Recent Selection

Most methods for fast detection of identity by descent (IBD) segments report identity by state segments without any quantification of the uncertainty in the endpoints and lengths of the IBD segments. We present a method for determining the posterior probability distribution of IBD segment endpoints. Our approach accounts for genotype errors, recent mutations, and gene conversions which disrupt DNA sequence identity within IBD segments, and it can be applied to large cohorts with whole-genome sequence or SNP array data. We find that our method's estimates of uncertainty are well calibrated for homogeneous samples. We quantify endpoint uncertainty for 77.7 billion IBD segments from 408,883 individuals of white British ancestry in the UK Biobank, and we use these IBD segments to find regions showing evidence of recent natural selection. We show that many spurious selection signals are eliminated by the use of unbiased estimates of IBD segment endpoints and a pedigree-based genetic map. Eleven of the twelve regions with the greatest evidence for recent selection in our scan have been identified as selected in previous analyses using different approaches. Our computationally efficient method for quantifying IBD segment endpoint uncertainty is implemented in the open source ibd-ends software package.

Functional Analysis and Fine Mapping of the 9p22.2 Ovarian Cancer Susceptibility Locus

Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. SIGNIFICANCE: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene.See related commentary by Choi and Brown, p. 439.

A GWAS in Latin Americans highlights the convergent evolution of lighter skin pigmentation in Eurasia

We report a genome-wide association scan in >6,000 Latin Americans for pigmentation of skin and eyes. We found eighteen signals of association at twelve genomic regions. These include one novel locus for skin pigmentation (in 10q26) and three novel loci for eye pigmentation (in 1q32, 20q13 and 22q12). We demonstrate the presence of multiple independent signals of association in the 11q14 and 15q13 regions (comprising the GRM5/TYR and HERC2/OCA2 genes, respectively) and several epistatic interactions among independently associated alleles. Strongest association with skin pigmentation at 19p13 was observed for an Y182H missense variant (common only in East Asians and Native Americans) in MFSD12, a gene recently associated with skin pigmentation in Africans. We show that the frequency of the derived allele at Y182H is significantly correlated with lower solar radiation intensity in East Asia and infer that MFSD12 was under selection in East Asians, probably after their split from Europeans.

Pigmentation variation in humans is influenced by complex genetic architecture in different populations. Here, the authors perform a genome-wide association analysis involving > 6,000 Latin Americans for pigmentation of skin and eyes, and identify known and novel genetic associations.

Shades of complexity: New perspectives on the evolution and genetic architecture of human skin

Like many highly variable human traits, more than a dozen genes are known to contribute to the full range of skin color. However, the historical bias in favor of genetic studies in European and European-derived populations has blinded us to the magnitude of pigmentation's complexity. As deliberate efforts are being made to better characterize diverse global populations and new sequencing technologies, better measurement tools, functional assessments, predictive modeling, and ancient DNA analyses become more widely accessible, we are beginning to appreciate how limited our understanding of the genetic bases of human skin color have been. Novel variants in genes not previously linked to pigmentation have been identified and evidence is mounting that there are hundreds more variants yet to be found. Even for genes that have been exhaustively characterized in European populations like MC1R, OCA2, and SLC24A5, research in previously understudied groups is leading to a new appreciation of the degree to which genetic diversity, epistatic interactions, pleiotropy, admixture, global and local adaptation, and cultural practices operate in population-specific ways to shape the genetic architecture of skin color. Furthermore, we are coming to terms with how factors like tanning response and barrier function may also have influenced selection on skin throughout human history. By examining how our knowledge of pigmentation genetics has shifted in the last decade, we can better appreciate how far we have come in understanding human diversity and the still long road ahead for understanding many complex human traits.

Genetic Hitchhiking and Population Bottlenecks Contribute to Prostate Cancer Disparities in Men of African Descent

Prostate cancer incidence and mortality rates in African and African American men are greatly elevated compared with other ethnicities. This disparity is likely explained by a combination of social, environmental, and genetic factors. A large number of susceptibility loci have been reported by genome-wide association studies (GWAS), but the contribution of these loci to prostate cancer disparities is unclear. Here, we investigated the population structure of 68 previously reported GWAS loci and calculated genetic disparity contribution statistics to identify SNPs that contribute the most to differences in prostate cancer risk across populations. By integrating GWAS results with allele frequency data, we generated genetic risk scores for 45 African and 19 non-African populations. Tests of natural selection were used to assess why some SNPs have large allele frequency differences across populations. We report that genetic predictions of prostate cancer risks are highest for West African men and lowest for East Asian men. These differences may be explained by the out-of-Africa bottleneck and natural selection. A small number of loci appear to drive elevated prostate cancer risks in men of African descent, including rs9623117, rs6983267, rs10896449, rs10993994, and rs817826. Although most prostate cancer-associated loci are evolving neutrally, there are multiple instances where alleles have hitchhiked to high frequencies with linked adaptive alleles. For example, a protective allele at 2q37 appears to have risen to high frequency in Europe due to selection acting on pigmentation. Our results suggest that evolutionary history contributes to the high rates of prostate cancer in African and African American men.Significance: A small number of genetic variants cause an elevated risk of prostate cancer in men of West African descent. Cancer Res; 78(9); 2432-43. ©2018 AACR.

Localization of adaptive variants in human genomes using averaged one-dependence estimation

Statistical methods for identifying adaptive mutations from population genetic data face several obstacles: assessing the significance of genomic outliers, integrating correlated measures of selection into one analytic framework, and distinguishing adaptive variants from hitchhiking neutral variants. Here, we introduce SWIF(r), a probabilistic method that detects selective sweeps by learning the distributions of multiple selection statistics under different evolutionary scenarios and calculating the posterior probability of a sweep at each genomic site. SWIF(r) is trained using simulations from a user-specified demographic model and explicitly models the joint distributions of selection statistics, thereby increasing its power to both identify regions undergoing sweeps and localize adaptive mutations. Using array and exome data from 45 ‡Khomani San hunter-gatherers of southern Africa, we identify an enrichment of adaptive signals in genes associated with metabolism and obesity. SWIF(r) provides a transparent probabilistic framework for localizing beneficial mutations that is extensible to a variety of evolutionary scenarios.

Genome-wide association study of pigmentary traits (skin and iris color) in individuals of East Asian ancestry

Background

Currently, there is limited knowledge about the genetics underlying pigmentary traits in East Asian populations. Here, we report the results of the first genome-wide association study of pigmentary traits (skin and iris color) in individuals of East Asian ancestry.

Methods

We obtained quantitative skin pigmentation measures (M-index) in the inner upper arm of the participants using a portable reflectometer (N = 305). Quantitative measures of iris color (expressed as L*, a* and b* CIELab coordinates) were extracted from high-resolution iris pictures (N = 342). We also measured the color differences between the pupillary and ciliary regions of the iris (e.g., iris heterochromia). DNA samples were genotyped with Illumina’s Infinium Multi-Ethnic Global Array (MEGA) and imputed using the 1000 Genomes Phase 3 samples as reference haplotypes.

Results

For skin pigmentation, we did not observe any genome-wide significant signal. We followed-up in three independent Chinese samples the lead SNPs of five regions showing multiple common markers (minor allele frequency ≥ 5%) with good imputation scores and suggestive evidence of association (p-values < 10⁻⁵). One of these markers, rs2373391, which is located in an intron of the ZNF804B gene on chromosome 7, was replicated in one of the Chinese samples (p = 0.003). For iris color, we observed genome-wide signals in the OCA2 region on chromosome 15. This signal is driven by the non-synonymous rs1800414 variant, which explains 11.9%, 10.4% and 6% of the variation observed in the b*, a* and L* coordinates in our sample, respectively. However, the OCA2 region was not associated with iris heterochromia.

Discussion

Additional genome-wide association studies in East Asian samples will be necessary to further disentangle the genetic architecture of pigmentary traits in East Asian populations.

Adaptation of human skin color in various populations

BACKGROUND

Skin color is a well-recognized adaptive trait and has been studied extensively in humans. Understanding the genetic basis of adaptation of skin color in various populations has many implications in human evolution and medicine.

DISCUSSION

Impressive progress has been made recently to identify genes associated with skin color variation in a wide range of geographical and temporal populations. In this review, we discuss what is currently known about the genetics of skin color variation. We enumerated several cases of skin color adaptation in global modern humans and archaic hominins, and illustrated why, when, and how skin color adaptation occurred in different populations. Finally, we provided a summary of the candidate loci associated with pigmentation, which could be a valuable reference for further evolutionary and medical studies.

CONCLUSION

Previous studies generally indicated a complex genetic mechanism underlying the skin color variation, expanding our understanding of the role of population demographic history and natural selection in shaping genetic and phenotypic diversity in humans. Future work is needed to dissect the genetic architecture of skin color adaptation in numerous ethnic minority groups around the world, which remains relatively obscure compared with that of major continental groups, and to unravel the exact genetic basis of skin color adaptation.

Evaluation of MC1R high-throughput nucleotide sequencing data generated by the 1000 Genomes Project

The advent of next-generation sequencing allows simultaneous processing of several genomic regions/individuals, increasing the availability and accuracy of whole-genome data. However, these new approaches may present some errors and bias due to alignment, genotype calling, and imputation methods. Despite these flaws, data obtained by next-generation sequencing can be valuable for population and evolutionary studies of specific genes, such as genes related to how pigmentation evolved among populations, one of the main topics in human evolutionary biology. Melanocortin-1 receptor (MC1R) is one of the most studied genes involved in pigmentation variation. As MC1R has already been suggested to affect melanogenesis and increase risk of developing melanoma, it constitutes one of the best models to understand how natural selection acts on pigmentation. Here we employed a locally developed pipeline to obtain genotype and haplotype data for MC1R from the raw sequencing data provided by the 1000 Genomes FTP site. We also compared such genotype data to Phase 3 VCF to evaluate its quality and discover any polymorphic sites that may have been overlooked. In conclusion, either the VCF file or one of the presently described pipelines could be used to obtain reliable and accurate genotype calling from the 1000 Genomes Phase 3 data.

Identifying signatures of positive selection in pigmentation genes in two South Asian populations

OBJECTIVES

Skin pigmentation is a polygenic trait showing wide phenotypic variations among global populations. While numerous pigmentation genes have been identified to be under positive selection among European and East populations, genes contributing to phenotypic variation in skin pigmentation within and among South Asian populations are still poorly understood. The present study uses data from the Phase 3 of the 1000 genomes project focusing on two South Asian populations-GIH (Gujarati Indian from Houston, Texas) and ITU (Indian Telugu from UK), so as to decode the genetic architecture involved in adaptation to ultraviolet radiation in South Asian populations.

METHODS

Statistical tests included were (1) tests 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).

RESULTS

Twenty-two pigmentation genes fall in the top 1% for at least one statistic in the GIH population, 5 of which (LYST, OCA2, SLC24A5, SLC45A2, and TYR) have been previously associated with normal variation in skin, hair, or eye color. In comparison, 17 genes fall in the top 1% for at least one statistic in the ITU population. Twelve loci which are identified as outliers in the ITU scan were also identified in the GIH population.

CONCLUSIONS

These results suggest that selection may have affected these loci broadly across the region.

Genome-wide data from two early Neolithic East Asian individuals dating to 7700 years ago

Ancient genomes have revolutionized our understanding of Holocene prehistory and, particularly, the Neolithic transition in western Eurasia. In contrast, East Asia has so far received little attention, despite representing a core region at which the Neolithic transition took place independently ~3 millennia after its onset in the Near East. We report genome-wide data from two hunter-gatherers from Devil's Gate, an early Neolithic cave site (dated to ~7.7 thousand years ago) located in East Asia, on the border between Russia and Korea. Both of these individuals are genetically most similar to geographically close modern populations from the Amur Basin, all speaking Tungusic languages, and, in particular, to the Ulchi. The similarity to nearby modern populations and the low levels of additional genetic material in the Ulchi imply a high level of genetic continuity in this region during the Holocene, a pattern that markedly contrasts with that reported for Europe.

Potential molecular characteristics in situ in response to repetitive UVB irradiation

Background

To identify molecular characteristics in situ in response to repetitive UVB (ultraviolet-B) irradiation.

Methods

Microarray data from the Gene Expression Omnibus were re-analyzed to identify DEGs (differentially expressed genes) between UVB-irradiated and non-irradiated skin biopsies. Enrichment and annotation analyses were performed respectively using DAVID, and TSGene and TAG databases. PPIs (protein-protein interactions) were analyzed using STRING, and miRNAs (microRNAs) and TFs (transcription factors) were predicted separately by miRNA-related databases and ENCODE. Accordingly, the PPI network and regulatory networks were visualized using Cytoscape, and they were merged together to obtain an integrated network for mining densely connected modules.

Results

Altogether, 151 up- and 64 down-regulated genes were identified between UVB-irradiated and non-irradiated skin biopsies, among which down-regulated DNAJB4 and SLIT2 were annotated as tumor-suppressors and up-regulated KIT was annotated as an oncogene. The up-regulated DEGs were significantly enriched in biological processes related to pigmentation (DCT, SOX10, TYRP1, TYR, MLPH, KIT and GPR143), while the down-regulated DEGs were dramatically related to haemopoiesis and the immune system (GPR183, INHBA, PTPRC, PLEK, CD8A and IKZF1). Furthermore, many miRNAs were screened for the DEGs, including miR-206 and miR-496 targeting KIT, miR-184 targeting DCT, and highly significant miR-337-5p, miR-21 and miR-16. Additionally, TFs were identified for the DEGs, among which PAX5 and HNF4A targeted MLPH and GPR143, respectively, while BATF, SPI1 and EP300 jointly target GPR183, PTPRC and PLEK.

Conclusions

The pigmentation and immune system implicated by DEGs, miRNAs and TFs might be important molecular mechanisms in response to UVB irradiation.

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.

Genes Regulated by Vitamin D in Bone Cells Are Positively Selected in East Asians

Vitamin D and folate are activated and degraded by sunlight, respectively, and the physiological processes they control are likely to have been targets of selection as humans expanded from Africa into Eurasia. We investigated signals of positive selection in gene sets involved in the metabolism, regulation and action of these two vitamins in worldwide populations sequenced by Phase I of the 1000 Genomes Project. Comparing allele frequency-spectrum-based summary statistics between these gene sets and matched control genes, we observed a selection signal specific to East Asians for a gene set associated with vitamin D action in bones. The selection signal was mainly driven by three genes CXXC finger protein 1 (CXXC1), low density lipoprotein receptor-related protein 5 (LRP5) and runt-related transcription factor 2 (RUNX2). Examination of population differentiation and haplotypes allowed us to identify several candidate causal regulatory variants in each gene. Four of these candidate variants (one each in CXXC1 and RUNX2 and two in LRP5) had a >70% derived allele frequency in East Asians, but were present at lower (20-60%) frequency in Europeans as well, suggesting that the adaptation might have been part of a common response to climatic and dietary changes as humans expanded out of Africa, with implications for their role in vitamin D-dependent bone mineralization and osteoporosis insurgence. We also observed haplotype sharing between East Asians, Finns and an extinct archaic human (Denisovan) sample at the CXXC1 locus, which is best explained by incomplete lineage sorting.

Distribution of two OCA2 polymorphisms associated with pigmentation in East-Asian populations

Two OCA2 polymorphisms (rs1800414 and rs74653330) have been associated with pigmentation in East Asians. We explored the distribution of these markers in a panel of samples from populations around the world. The derived allele of rs1800414 has high frequencies in a broad East-Asian region, whereas the derived allele of rs74653330 is primarily restricted to northern East Asia. Our data suggest that these polymorphisms may have been selected independently in different regions of East Asia.

Genetics of autoimmune diseases: insights from population genetics

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Autoimmune diseases (ADs) are a family of complex heterogeneous disorders with similar underlying mechanisms characterized by immune responses against self. Collectively, ADs are common, exhibit gender and ethnic disparities, and increasing incidence. As natural selection is an important influence on human genetic variation, and immune function genes are enriched for signals of positive selection, it is thought that the prevalence of AD risk alleles seen in different population is partially the result of differing selective pressures (for example, due to pathogens). With the advent of high-throughput technologies, new analytical methodologies and large-scale projects, evidence for the role of natural selection in contributing to the heritable component of ADs keeps growing. This review summarizes the genetic regions associated with susceptibility to different ADs and concomitant evidence for selection, including known agents of selection exerting selective pressure in these regions. Examples of specific adaptive variants with phenotypic effects are included as an evidence of natural selection increasing AD susceptibility. Many of the complexities of gene effects in different ADs can be explained by population genetics phenomena. Integrating AD susceptibility studies with population genetics to investigate how natural selection has contributed to genetic variation that influences disease risk will help to identify functional variants and elucidate biological mechanisms. As such, the study of population genetics in human population holds untapped potential for elucidating the genetic causes of human disease and more rapidly focusing to personalized medicine.

Selective Strolls: Fixation and Extinction in Diploids Are Slower for Weakly Selected Mutations Than for Neutral Ones

In finite populations, an allele disappears or reaches fixation due to two main forces, selection and drift. Selection is generally thought to accelerate the process: a selected mutation will reach fixation faster than a neutral one, and a disadvantageous one will quickly disappear from the population. We show that even in simple diploid populations, this is often not true. Dominance and recessivity unexpectedly slow down the evolutionary process for weakly selected alleles. In particular, slightly advantageous dominant and mildly deleterious recessive mutations reach fixation slightly more slowly than neutral ones (at most 5%). This phenomenon determines genetic signatures opposite to those expected under strong selection, such as increased instead of decreased genetic diversity around the selected site. Furthermore, we characterize a new phenomenon: mildly deleterious recessive alleles, thought to represent a wide fraction of newly arising mutations, on average survive in a population slightly longer than neutral ones, before getting lost. Consequently, these mutations are on average slightly older than neutral ones, in contrast with previous expectations. Furthermore, they slightly increase the amount of weakly deleterious polymorphisms, as a consequence of the longer unconditional sojourn times compared to neutral mutations.

Deep History of East Asian Populations Revealed Through Genetic Analysis of the Ainu

Despite recent advances in population genomics, much remains to be elucidated with regard to East Asian population history. The Ainu, a hunter-gatherer population of northern Japan and Sakhalin island of Russia, are thought to be key to elucidating the prehistory of Japan and the peopling of East Asia. Here, we study the genetic relationship of the Ainu with other East Asian and Siberian populations outside the Japanese archipelago using genome-wide genotyping data. We find that the Ainu represent a deep branch of East Asian diversity more basal than all present-day East Asian farmers. However, we did not find a genetic connection between the Ainu and populations of the Tibetan plateau, rejecting their long-held hypothetical connection based on Y chromosome data. Unlike all other East Asian populations investigated, the Ainu have a closer genetic relationship with northeast Siberians than with central Siberians, suggesting ancient connections among populations around the Sea of Okhotsk. We also detect a recent genetic contribution of the Ainu to nearby populations, but no evidence for reciprocal recent gene flow is observed. Whole genome sequencing of contemporary and ancient Ainu individuals will be helpful to understand the details of the deep history of East Asians.

MC1R diversity in Northern Island Melanesia has not been constrained by strong purifying selection and cannot explain pigmentation phenotype variation in the region

BACKGROUND

Variation in human skin pigmentation evolved in response to the selective pressure of ultra-violet radiation (UVR). Selection to maintain darker skin in high UVR environments is expected to constrain pigmentation phenotype and variation in pigmentation loci. Consistent with this hypothesis, the gene MC1R exhibits reduced diversity in African populations from high UVR regions compared to low-UVR non-African populations. However, MC1R diversity in non-African populations that have evolved under high-UVR conditions is not well characterized.

METHODS

In order to test the hypothesis that MC1R variation has been constrained in Melanesians the coding region of the MC1R gene was sequenced in 188 individuals from Northern Island Melanesia. The role of purifying selection was assessed using a modified McDonald Kreitman's test. Pairwise FST was calculated between Melanesian populations and populations from the 1000 Genomes Project. The SNP rs2228479 was genotyped in a larger sample (n = 635) of Melanesians and tested for associations with skin and hair pigmentation.

RESULTS

We observe three nonsynonymous and two synonymous mutations. A modified McDonald Kreitman's test failed to detect a significant signal of purifying selection. Pairwise FST values calculated between the four islands sampled here indicate little regional substructure in MC1R. When compared to African, European, East and South Asian populations, Melanesians do not exhibit reduced population divergence (measured as FST) or a high proportion of haplotype sharing with Africans, as one might expect if ancestral haplotypes were conserved across high UVR populations in and out of Africa. The only common nonsynonymous polymorphism observed, rs2228479, is not significantly associated with skin or hair pigmentation in a larger sample of Melanesians.

CONCLUSIONS

The pattern of sequence diversity here does not support a model of strong selective constraint on MC1R in Northern Island Melanesia This absence of strong constraint, as well as the recent population history of the region, may explain the observed frequencies of the derived rs2228479 allele. These results emphasize the complex genetic architecture of pigmentation phenotypes, which are controlled by multiple, possibly interacting loci. They also highlight the role that population history can play in influencing phenotypic diversity in the absence of strong natural selection.

Adaptive Evolution Favoring KLK4 Downregulation in East Asians

The human kallikrein (KLK) cluster, located at chromosome 19q13.3-13.4, encodes 15 serine proteases, including neighboring genes (KLK3, KLK2, KLK4, and KLK5) with key roles in the cascades of semen liquefaction, tooth enamel maturation, and skin desquamation. KLK2 and KLK3 were previously identified as targets of adaptive evolution in primates through different mechanisms linked to reproductive biology and, in humans, genome-wide scans of positive selection captured, a yet unexplored, evidence for KLK neutrality departure in East Asians. We perform a detailed evaluation of KLK3-KLK5 variability in the 1000 Genomes samples from East Asia, Europe, and Africa, which was sustained by our own sequencing. In East Asians, we singled out a 70-kb region surrounding KLK4 that combined unusual low levels of diversity, high frequency variants with significant levels of population differentiation (FST > 0.5) and fairly homogenous haplotypes given the large local recombination rates. Among these variants, rs1654556_G, rs198968_T, and rs17800874_A stand out for their location on putative regulatory regions and predicted functional effects, namely the introduction of several microRNA binding sites and a repressor motif. Our functional assays carried out in different cellular models showed that rs198968_T and rs17800874_A operate synergistically to reduce KLK4 expression and could be further assisted by rs1654556_G. Considering the previous findings that KLK4 inactivation causes enamel malformations in humans and mice, and that this gene is coexpressed in epidermal layers along with several substrates involved in either cell adhesion or keratinocyte differentiation, we propose KLK4 as another target of selection in East Asians correlated to tooth and epidermal morphological traits.

Fifteen years of genomewide scans for selection: trends, lessons and unaddressed genetic sources of complication

Genomewide scans for natural selection (GWSS) have become increasingly common over the last 15 years due to increased availability of genome-scale genetic data. Here, we report a representative survey of GWSS from 1999 to present and find that (i) between 1999 and 2009, 35 of 49 (71%) GWSS focused on human, while from 2010 to present, only 38 of 83 (46%) of GWSS focused on human, indicating increased focus on nonmodel organisms; (ii) the large majority of GWSS incorporate interpopulation or interspecific comparisons using, for example F(ST), cross-population extended haplotype homozygosity or the ratio of nonsynonymous to synonymous substitutions; (iii) most GWSS focus on detection of directional selection rather than other modes such as balancing selection; and (iv) in human GWSS, there is a clear shift after 2004 from microsatellite markers to dense SNP data. A survey of GWSS meant to identify loci positively selected in response to severe hypoxic conditions support an approach to GWSS in which a list of a priori candidate genes based on potential selective pressures are used to filter the list of significant hits a posteriori. We also discuss four frequently ignored determinants of genomic heterogeneity that complicate GWSS: mutation, recombination, selection and the genetic architecture of adaptive traits. We recommend that GWSS methodology should better incorporate aspects of genomewide heterogeneity using empirical estimates of relevant parameters and/or realistic, whole-chromosome simulations to improve interpretation of GWSS results. Finally, we argue that knowledge of potential selective agents improves interpretation of GWSS results and that new methods focused on correlations between environmental variables and genetic variation can help automate this approach.

Comparison of the Transcriptional Profiles of Melanocytes from Dark and Light Skinned Individuals under Basal Conditions and Following Ultraviolet-B Irradiation

We analysed the whole-genome transcriptional profile of 6 cell lines of dark melanocytes (DM) and 6 of light melanocytes (LM) at basal conditions and after ultraviolet-B (UVB) radiation at different time points to investigate the mechanisms by which melanocytes protect human skin from the damaging effects of UVB. Further, we assessed the effect of different keratinocyte-conditioned media (KCM+ and KCM-) on melanocytes. Our results suggest that an interaction between ribosomal proteins and the P53 signaling pathway may occur in response to UVB in both DM and LM. We also observed that DM and LM show differentially expressed genes after irradiation, in particular at the first 6h after UVB. These are mainly associated with inflammatory reactions, cell survival or melanoma. Furthermore, the culture with KCM+ compared with KCM- had a noticeable effect on LM. This effect includes the activation of various signaling pathways such as the mTOR pathway, involved in the regulation of cell metabolism, growth, proliferation and survival. Finally, the comparison of the transcriptional profiles between LM and DM under basal conditions, and the application of natural selection tests in human populations allowed us to support the significant evolutionary role of MIF and ATP6V0B in the pigmentary phenotype.

Defining window-boundaries for genomic analyses using smoothing spline techniques

Background

High-density genomic data is often analyzed by combining information over windows of adjacent markers. Interpretation of data grouped in windows versus at individual locations may increase statistical power, simplify computation, reduce sampling noise, and reduce the total number of tests performed. However, use of adjacent marker information can result in over- or under-smoothing, undesirable window boundary specifications, or highly correlated test statistics. We introduce a method for defining windows based on statistically guided breakpoints in the data, as a foundation for the analysis of multiple adjacent data points. This method involves first fitting a cubic smoothing spline to the data and then identifying the inflection points of the fitted spline, which serve as the boundaries of adjacent windows. This technique does not require prior knowledge of linkage disequilibrium, and therefore can be applied to data collected from individual or pooled sequencing experiments. Moreover, in contrast to existing methods, an arbitrary choice of window size is not necessary, since these are determined empirically and allowed to vary along the genome.

Results

Simulations applying this method were performed to identify selection signatures from pooled sequencing F_ST data, for which allele frequencies were estimated from a pool of individuals. The relative ratio of true to false positives was twice that generated by existing techniques. A comparison of the approach to a previous study that involved pooled sequencing F_ST data from maize suggested that outlying windows were more clearly separated from their neighbors than when using a standard sliding window approach.

Conclusions

We have developed a novel technique to identify window boundaries for subsequent analysis protocols. When applied to selection studies based on F_ST data, this method provides a high discovery rate and minimizes false positives. The method is implemented in the R package GenWin, which is publicly available from CRAN.

The Two-pore channel (TPC) interactome unmasks isoform-specific roles for TPCs in endolysosomal morphology and cell pigmentation

The two-pore channels (TPC1 and TPC2) belong to an ancient family of intracellular ion channels expressed in the endolysosomal system. Little is known about how regulatory inputs converge to modulate TPC activity, and proposed activation mechanisms are controversial. Here, we compiled a proteomic characterization of the human TPC interactome, which revealed that TPCs complex with many proteins involved in Ca(2+) homeostasis, trafficking, and membrane organization. Among these interactors, TPCs were resolved to scaffold Rab GTPases and regulate endomembrane dynamics in an isoform-specific manner. TPC2, but not TPC1, caused a proliferation of endolysosomal structures, dysregulating intracellular trafficking, and cellular pigmentation. These outcomes required both TPC2 and Rab activity, as well as their interactivity, because TPC2 mutants that were inactive, or rerouted away from their endogenous expression locale, or deficient in Rab binding, failed to replicate these outcomes. Nicotinic acid adenine dinucleotide phosphate (NAADP)-evoked Ca(2+) release was also impaired using either a Rab binding-defective TPC2 mutant or a Rab inhibitor. These data suggest a fundamental role for the ancient TPC complex in trafficking that holds relevance for lysosomal proliferative scenarios observed in disease.

Human skin color is influenced by an intergenic DNA polymorphism regulating transcription of the nearby BNC2 pigmentation gene

Single nucleotide polymorphisms (SNPs) found to be statistically significant when associated with human diseases, and other phenotypes are most often located in non-coding regions of the genome. One example is rs10765819 located in the first intron of the BNC2 gene previously associated with (saturation of) human skin color. Here, we demonstrate that a nearby intergenic SNP (rs12350739) in high linkage disequilibrium with rs10756819 is likely the causal DNA variant for the observed BNC2 skin color association. The highly conserved region surrounding rs12350739 functions as an enhancer element regulating BNC2 transcription in human melanocytes, while the activity of this enhancer element depends on the allelic status of rs12350739. When the rs12350739-AA allele is present, the chromatin at the region surrounding rs12350739 is inaccessible and the enhancer element is only slightly active, resulting in low expression of BNC2, corresponding with light skin pigmentation. When the rs12350739-GG allele is present however, the chromatin at the region surrounding rs12350739 is more accessible and the enhancer is active, resulting in a higher expression of BNC2, corresponding with dark skin pigmentation. Overall, we demonstrate the identification of the functional DNA variant that explains the BNC2 skin color association signal, providing another important step towards further understanding human pigmentation genetics beyond statistical association. We thus deliver a clear example of how an intergenic non-coding DNA variant modulates the regulatory potential of the enhancer element it is located within, which in turn results in allele-dependent differential gene expression affecting variation in common human traits.

Patterns of coding variation in the complete exomes of three Neandertals

We present the DNA sequence of 17,367 protein-coding genes in two Neandertals from Spain and Croatia and analyze them together with the genome sequence recently determined from a Neandertal from southern Siberia. Comparisons with present-day humans from Africa, Europe, and Asia reveal that genetic diversity among Neandertals was remarkably low, and that they carried a higher proportion of amino acid-changing (nonsynonymous) alleles inferred to alter protein structure or function than present-day humans. Thus, Neandertals across Eurasia had a smaller long-term effective population than present-day humans. We also identify amino acid substitutions in Neandertals and present-day humans that may underlie phenotypic differences between the two groups. We find that genes involved in skeletal morphology have changed more in the lineage leading to Neandertals than in the ancestral lineage common to archaic and modern humans, whereas genes involved in behavior and pigmentation have changed more on the modern human lineage.