Genetic architecture of skin and eye color in an African-European admixed population

Why humans evolved blue eyes

A surprising number of humans are equipped with a subpar eye model-featuring pale, colorful irides that are nowhere as good as the original dark ones at guarding the retina from sunlight and do, in fact, raise one's risk of eye disease. Here I apply evolutionary theory to understand why. I propose that the allele for human blue eyes, which arose just once, managed to spread from one individual to millions at an astonishing speed because it is a greenbeard. "Greenbeards"-imaginary genes, or groups of genes, that produce both a green beard and a behavior that favors other bearers of a green beard-have been deemed exceedingly unlikely to show up in the real world. And yet, as individuals who prefer blue eyes are more inclined to mate with blue-eyed partners and invest in blue-eyed offspring, any blue-eye preference (whether random or arising from the bias for colorful stimuli shared by all recognition systems) becomes rapidly linked to the blue-eye trait. Thus, blue eyes gain an edge by working like a peacock's colorful tail and a nestling's colorful mouth: twice self-reinforcing, "double runaway" evolution via sexual and parental selection. The blue-eye ornament gene, by binding to a behavior that favors other bearers of the blue-eye ornament gene, is ultimately recognizing and helping copies of itself in both kin and strangers-and greatly prospering, just like theory predicts.

Genome wide association study and meta-analysis identified multiple new risk loci for freckles in 4813 Chinese individuals

Freckle is a prevalent pigmentary dermatosis with an obvious hereditary component. Dozens of freckles risk loci have been discovered through research on multiple traits or other diseases, rather than as an independent trait. To discover novel variants associated with freckles, we performed GWAS and meta-analysis in 4813 Chinese individuals. We conducted GWAS and meta-analysis of two cohorts: 197 patients and 1603 controls (Cohort I), and 336 patients and 2677 controls (Cohort II), both from China. Then we performed linkage disequilibrium (LD) analysis, eQTL study, and enrichment analysis with association results for functional implications. Finally, we discovered 59 new SNPs and 13 novel susceptibility genes associated with freckles (Pmeta <5 × 10-8), which has enriched the genetic research on freckles.

Sickle Cell Anemia Screening in Newborns and Analysis of Haplotypes in Patients from Santiago Island, Cape Verde

Sickle cell anemia (SCA) results from a mutation in the β-globin gene, leading to the production of mutant hemoglobin, known as hemoglobin S (HbS). Despite being a genetic disorder, the phenotype of SCA can be influenced by the level of fetal hemoglobin (HbF), which is associated with beta S-globin haplotypes. In this study, we conducted newborn screening (NBS) using samples collected from umbilical cord blood in two hospitals on Santiago Island, Cape Verde. In newborns, HbS was detected using high-performance liquid chromatography (HPLC) on dried blood spot, with confirmation through polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP). In addition, we assessed the hematological and clinical characteristics of a second population group consisting of patients diagnosed with SCA. Haplotype determination was performed on both newborns with HbS and patients with SCA. Beta S-globin haplotypes were determined using PCR-RFLP. Hematological values were analyzed using standard methods. Out of 346 newborns, 21 (6%) were carriers of the sickle cell trait (HbAS) while none were identified as homozygous for sickle cell disease (HbSS). Among both groups of individuals, four haplotypes were identified: Senegal, Arabi-Indian, Bantu, and Benin. The Senegal haplotype was the most prevalent, possibly reflecting the ethnic origin of the mutations observed. Hematological values did not differ significantly among haplotypes. However, higher levels of HbF were associated with better hematological values. These findings suggest a positive impact of elevated HbF levels on reducing the severity of SCA. Finally, we demonstrated how the combination of technics, HPLC and molecular analysis, provided a consistent and reproducible results that can be used for NBS for SCA.

Interpreting SNP heritability in admixed populations

SNP heritabilityh s n p 2 is defined as the proportion of phenotypic variance explained by genotyped SNPs and is believed to be a lower bound of heritability (h 2 ), being equal to it if all causal variants are known. Despite the simple intuition behindh s n p 2 , its interpretation and equivalence toh 2 is unclear, particularly in the presence of population structure and assortative mating. It is well known that population structure can lead to inflation inh ˆ s n p 2 estimates because of confounding due to linkage disequilibrium (LD) or shared environment. Here we use analytical theory and simulations to demonstrate thath s n p 2 estimates can be biased in admixed populations, even in the absence of confounding and even if all causal variants are known. This is because admixture generates LD, which contributes to the genetic variance, and therefore to heritability. Genome-wide restricted maximum likelihood (GREML) does not capture this contribution leading to under- or over-estimates ofh s n p 2 relative toh 2 , depending on the genetic architecture. In contrast, Haseman-Elston (HE) regression exaggerates the LD contribution leading to biases in the opposite direction. For the same reason, GREML and HE estimates of local ancestry heritabilityh γ 2 are also biased. We describe this bias inh ˆ s n p 2 andh ˆ γ 2 as a function of admixture history and the genetic architecture of the trait and show that it can be recovered under some conditions. We clarify the interpretation ofh ˆ s n p 2 in admixed populations and discuss its implication for genome-wide association studies and polygenic prediction.

Weakened tanning ability is an important mechanism for evolutionary skin lightening in East Asians

The evolution of light-skin pigmentation among Eurasians is considered as an adaptation to the high-latitude environments. East Asians are ideal populations for studying skin color evolution because of the complex environment of East Asia. Here, we report a strong selection signal for the pigmentation gene phenylalanine hydroxylase (PAH) in light-skinned Han Chinese individuals. The intron mutation rs10778203 in PAH is enriched in East Asians and is significantly associated with skin color of the back of the hand in Han Chinese males (P < 0.05). In vitro luciferase and transcription factor binding assays show that the ancestral allele of rs10778203 could bind to SMAD2 and has a significant enhancer activity for PAH. However, the derived T allele (the major allele in East Asians) of rs10778203 decreases the binding activity of transcription factors and enhancer activity. Meanwhile, the derived T allele of rs10778203 shows a weaker ultraviolet radiation response in A375 cells and zebrafish embryos. Furthermore, rs10778203 decreases melanin production in transgenic zebrafish embryos after ultraviolet B (UVB) treatment. Collectively, PAH is a potential pigmentation gene that regulates skin tanning ability. Natural selection has enriched the adaptive allele, resulting in weakened tanning ability in East Asians, suggesting a unique genetic mechanism for evolutionary skin lightening in East Asians.

Global and Local Ancestry and its Importance: A Review

The fastest way to significantly change the composition of a population is through admixture, an evolutionary mechanism. In animal breeding history, genetic admixture has provided both short-term and long-term advantages by utilizing the phenomenon of complementarity and heterosis in several traits and genetic diversity, respectively. The traditional method of admixture analysis by pedigree records has now been replaced greatly by genome-wide marker data that enables more precise estimations. Among these markers, SNPs have been the popular choice since they are cost-effective, not so laborious, and automation of genotyping is easy. Certain markers can suggest the possibility of a population's origin from a sample of DNA where the source individual is unknown or unwilling to disclose their lineage, which are called Ancestry-Informative Markers (AIMs). Revealing admixture level at the locus-specific level is termed as local ancestry and can be exploited to identify signs of recent selective response and can account for genetic drift. Considering the importance of genetic admixture and local ancestry, in this mini-review, both concepts are illustrated, encompassing basics, their estimation/identification methods, tools/software used and their applications.

iGWAS: Image-based genome-wide association of self-supervised deep phenotyping of retina fundus images

Existing imaging genetics studies have been mostly limited in scope by using imaging-derived phenotypes defined by human experts. Here, leveraging new breakthroughs in self-supervised deep representation learning, we propose a new approach, image-based genome-wide association study (iGWAS), for identifying genetic factors associated with phenotypes discovered from medical images using contrastive learning. Using retinal fundus photos, our model extracts a 128-dimensional vector representing features of the retina as phenotypes. After training the model on 40,000 images from the EyePACS dataset, we generated phenotypes from 130,329 images of 65,629 British White participants in the UK Biobank. We conducted GWAS on these phenotypes and identified 14 loci with genome-wide significance (p<5×10-8 and intersection of hits from left and right eyes). We also did GWAS on the retina color, the average color of the center region of the retinal fundus photos. The GWAS of retina colors identified 34 loci, 7 are overlapping with GWAS of raw image phenotype. Our results establish the feasibility of this new framework of genomic study based on self-supervised phenotyping of medical images.

Integrative functional genomic analyses identify genetic variants influencing skin pigmentation in Africans

Skin color is highly variable in Africans, yet little is known about the underlying molecular mechanism. Here we applied massively parallel reporter assays to screen 1,157 candidate variants influencing skin pigmentation in Africans and identified 165 single-nucleotide polymorphisms showing differential regulatory activities between alleles. We combine Hi-C, genome editing and melanin assays to identify regulatory elements for MFSD12, HMG20B, OCA2, MITF, LEF1, TRPS1, BLOC1S6 and CYB561A3 that impact melanin levels in vitro and modulate human skin color. We found that independent mutations in an OCA2 enhancer contribute to the evolution of human skin color diversity and detect signals of local adaptation at enhancers of MITF, LEF1 and TRPS1, which may contribute to the light skin color of Khoesan-speaking populations from Southern Africa. Additionally, we identified CYB561A3 as a novel pigmentation regulator that impacts genes involved in oxidative phosphorylation and melanogenesis. These results provide insights into the mechanisms underlying human skin color diversity and adaptive evolution.

A genome-wide genetic screen uncovers determinants of human pigmentation

Skin color, one of the most diverse human traits, is determined by the quantity, type, and distribution of melanin. In this study, we leveraged the light-scattering properties of melanin to conduct a genome-wide screen for regulators of melanogenesis. We identified 169 functionally diverse genes that converge on melanosome biogenesis, endosomal transport, and gene regulation, of which 135 represented previously unknown associations with pigmentation. In agreement with their melanin-promoting function, the majority of screen hits were up-regulated in melanocytes from darkly pigmented individuals. We further unraveled functions of KLF6 as a transcription factor that regulates melanosome maturation and pigmentation in vivo, and of the endosomal trafficking protein COMMD3 in modulating melanosomal pH. Our study reveals a plethora of melanin-promoting genes, with broad implications for human variation, cell biology, and medicine.

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

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

Localizing Post-Admixture Adaptive Variants with Object Detection on Ancestry-Painted Chromosomes

Gene flow between previously differentiated populations during the founding of an admixed or hybrid population has the potential to introduce adaptive alleles into the new population. If the adaptive allele is common in one source population, but not the other, then as the adaptive allele rises in frequency in the admixed population, genetic ancestry from the source containing the adaptive allele will increase nearby as well. Patterns of genetic ancestry have therefore been used to identify post-admixture positive selection in humans and other animals, including examples in immunity, metabolism, and animal coloration. A common method identifies regions of the genome that have local ancestry "outliers" compared with the distribution across the rest of the genome, considering each locus independently. However, we lack theoretical models for expected distributions of ancestry under various demographic scenarios, resulting in potential false positives and false negatives. Further, ancestry patterns between distant sites are often not independent. As a result, current methods tend to infer wide genomic regions containing many genes as under selection, limiting biological interpretation. Instead, we develop a deep learning object detection method applied to images generated from local ancestry-painted genomes. This approach preserves information from the surrounding genomic context and avoids potential pitfalls of user-defined summary statistics. We find the method is robust to a variety of demographic misspecifications using simulated data. Applied to human genotype data from Cabo Verde, we localize a known adaptive locus to a single narrow region compared with multiple or long windows obtained using two other ancestry-based methods.

Exploring regional aspects of 3D facial variation within European individuals

Facial ancestry can be described as variation that exists in facial features that are shared amongst members of a population due to environmental and genetic effects. Even within Europe, faces vary among subregions and may lead to confounding in genetic association studies if unaccounted for. Genetic studies use genetic principal components (PCs) to describe facial ancestry to circumvent this issue. Yet the phenotypic effect of these genetic PCs on the face has yet to be described, and phenotype-based alternatives compared. In anthropological studies, consensus faces are utilized as they depict a phenotypic, not genetic, ancestry effect. In this study, we explored the effects of regional differences on facial ancestry in 744 Europeans using genetic and anthropological approaches. Both showed similar ancestry effects between subgroups, localized mainly to the forehead, nose, and chin. Consensus faces explained the variation seen in only the first three genetic PCs, differing more in magnitude than shape change. Here we show only minor differences between the two methods and discuss a combined approach as a possible alternative for facial scan correction that is less cohort dependent, more replicable, non-linear, and can be made open access for use across research groups, enhancing future studies in this field.

Genetic architecture and evolution of color variation in American black bears

Color variation is a frequent evolutionary substrate for camouflage in small mammals, but the underlying genetics and evolutionary forces that drive color variation in natural populations of large mammals are mostly unexplained. The American black bear, Ursus americanus (U. americanus), exhibits a range of colors including the cinnamon morph, which has a similar color to the brown bear, U. arctos, and is found at high frequency in the American southwest. Reflectance and chemical melanin measurements showed little distinction between U. arctos and cinnamon U. americanus individuals. We used a genome-wide association for hair color as a quantitative trait in 151 U. americanus individuals and identified a single major locus (p < 10-13). Additional genomic and functional studies identified a missense alteration (R153C) in Tyrosinase-related protein 1 (TYRP1) that likely affects binding of the zinc cofactor, impairs protein localization, and results in decreased pigment production. Population genetic analyses and demographic modeling indicated that the R153C variant arose 9.36 kya in a southwestern population where it likely provided a selective advantage, spreading both northwards and eastwards by gene flow. A different TYRP1 allele, R114C, contributes to the characteristic brown color of U. arctos but is not fixed across the range.

Candidate pigmentation genes related to feather color variation in an indigenous chicken breed revealed by whole genome data

Chicken plumage color is an inheritable phenotype that was naturally and artificially selected for during domestication. The Baicheng You chicken is an indigenous Chinese chicken breed presenting three main feather colors, lavender, black, and yellow plumages. To explore the genetic mechanisms underlying the pigmentation in Baicheng You chickens, we re-sequenced the whole genome of Baicheng You chicken with the three plumage colors. By analyzing the divergent regions of the genome among the chickens with different feather colors, we identified some candidate genomic regions associated with the feather colors in Baicheng You chickens. We found that EGR1, MLPH, RAB17, SOX5, and GRM5 genes were the potential genes for black, lavender, and yellow feathers. MLPH, GRM5, and SOX5 genes have been found to be related to plumage colors in birds. Our results showed that EGR1 is a most plausible candidate gene for black plumage, RAB17, MLPH, and SOX5 for lavender plumage, and GRM5 for yellow plumage in Baicheng You chicken.

Evaluation of Personalized Skincare Through in-silico Gene Interactive Networks and Cellular Responses to UVR and Oxidative Stress

Purpose

Personalized approaches in dermatology are designed to match the specific requirements based on the individual genetic makeup. One major factor accounting for the differences in skin phenotypes is single nucleotide polymorphism (SNP) within several genes with diverse roles that extend beyond skin tone and pigmentation. Therefore, the cellular sensitivities to the environmental stress and damage linked to extrinsic aging could also underlie the individual characteristics of the skin and dictate the unique skin care requirements. This study aimed to identify the likely biomarkers and molecular signatures expressed in skin cells of different ethnic backgrounds, which could aid further the design of personalized skin products based on specific demands.

Methods

Using data mining and in-silico modeling, the association of SNP-affected genes with three major skin types of European, Asian and African origin was analyzed and compared within the structure-function gene interaction networks. Cultured dermal fibroblasts were subsequently subjected to ultraviolet radiation and oxidative stress and analyzed for DNA damage and senescent markers. The protective applications of two cosmetic ingredients, Resveratrol and Quercetin, were validated in both cellular and in-silico models.

Results

Each skin type was characterized by the presence of SNPs in the genes controlling facultative and constitutive pigmentation, which could also underlie the major differences in responses to photodamage, such as oxidative stress, inflammation, and barrier homeostasis. Skin-type-specific dermal fibroblasts cultured in-vitro demonstrated distinctive sensitivities to ultraviolet radiation and oxidative stress, which could be modulated further by the bioactive compounds with the predicted capacities to interact with some of the genes in the in-silico models.

Conclusion

Evaluation of the SNP-affected gene networks and likely sensitivities of skin cells, defined as low threshold levels to extrinsic stress factors, can provide a valuable tool for the design and formulation of personalized skin products that match more accurately diverse ethnic backgrounds.

Genetic adaptation of skin pigmentation in highland Tibetans

Strong ultraviolet (UV) radiation at high altitude imposes a serious selective pressure, which may induce skin pigmentation adaptation of indigenous populations. We conducted skin pigmentation phenotyping and genome-wide analysis of Tibetans in order to understand the underlying mechanism of adaptation to UV radiation. We observe that Tibetans have darker baseline skin color compared with lowland Han Chinese, as well as an improved tanning ability, suggesting a two-level adaptation to boost their melanin production. A genome-wide search for the responsible genes identifies GNPAT showing strong signals of positive selection in Tibetans. An enhancer mutation (rs75356281) located in GNPAT intron 2 is enriched in Tibetans (58%) but rare in other world populations (0 to 18%). The adaptive allele of rs75356281 is associated with darker skin in Tibetans and, under UVB treatment, it displays higher enhancer activities compared with the wild-type allele in in vitro luciferase assays. Transcriptome analyses of gene-edited cells clearly show that with UVB treatment, the adaptive variant of GNPAT promotes melanin synthesis, likely through the interactions of CAT and ACAA1 in peroxisomes with other pigmentation genes, and they act synergistically, leading to an improved tanning ability in Tibetans for UV protection.

Sex-biased admixture and assortative mating shape genetic variation and influence demographic inference in admixed Cabo Verdeans

Genetic data can provide insights into population history, but first, we must understand the patterns that complex histories leave in genomes. Here, we consider the admixed human population of Cabo Verde to understand the patterns of genetic variation left by social and demographic processes. First settled in the late 1400s, Cabo Verdeans are admixed descendants of Portuguese colonizers and enslaved West African people. We consider Cabo Verde's well-studied historical record alongside genome-wide SNP data from 563 individuals from 4 regions within the archipelago. We use genetic ancestry to test for patterns of nonrandom mating and sex-specific gene flow, and we examine the consequences of these processes for common demographic inference methods and genetic patterns. Notably, multiple population genetic tools that assume random mating underestimate the timing of admixture, but incorporating nonrandom mating produces estimates more consistent with historical records. We consider how admixture interrupts common summaries of genomic variation such as runs of homozygosity. While summaries of runs of homozygosity may be difficult to interpret in admixed populations, differentiating runs of homozygosity by length class shows that runs of homozygosity reflect historical differences between the islands in their contributions from the source populations and postadmixture population dynamics. Finally, we find higher African ancestry on the X chromosome than on the autosomes, consistent with an excess of European males and African females contributing to the gene pool. Considering these genomic insights into population history in the context of Cabo Verde's historical record, we can identify how assumptions in genetic models impact inference of population history more broadly.

Conditional GWAS of non-CG transposon methylation in Arabidopsis thaliana reveals major polymorphisms in five genes

Genome-wide association studies (GWAS) have revealed that the striking natural variation for DNA CHH-methylation (mCHH; H is A, T, or C) of transposons has oligogenic architecture involving major alleles at a handful of known methylation regulators. Here we use a conditional GWAS approach to show that CHG-methylation (mCHG) has a similar genetic architecture-once mCHH is statistically controlled for. We identify five key trans-regulators that appear to modulate mCHG levels, and show that they interact with a previously identified modifier of mCHH in regulating natural transposon mobilization.

Applicability of the IrisPlex system for eye color prediction in an admixed population from Argentina

Eye color prediction based on an individual's genetic information is of interest in the field of forensic genetics. In recent years, researchers have studied different genes and markers associated with this externally visible characteristic and have developed methods for its prediction. The IrisPlex represents a validated tool for homogeneous populations, though its applicability in populations of mixed ancestry is limited, mainly regarding the prediction of intermediate eye colors. With the aim of validating the applicability of this system in an admixed population from Argentina (n = 302), we analyzed the six single nucleotide variants used in that multiplex for eye color and four additional SNPs, and evaluated its prediction ability. We also performed a genotype-phenotype association analysis. This system proved to be useful when dealing with the extreme ends of the eye color spectrum (blue and brown) but presented difficulties in determining the intermediate phenotypes (green), which were found in a large proportion of our population. We concluded that these genetic tools should be used with caution in admixed populations and that more studies are required in order to improve the prediction of intermediate phenotypes.

Which field of research would Gregor Mendel choose in the 21st century?

Gregor Mendel's work on segregation of traits in plants established the basic methodology and rules of genetics. The interruption of Mendel's research activities in 1870 impeded the immediate recognition of the value of his work until the dawn of the 20th century. Only then were his founding laws of genetics validated, propelling the development of biological research toward the birth of molecular biology in the second half of the 20th century. While molecular plant genetics can be viewed as the spiritual heir of Mendel's research, one might wonder whether in the 21st century Gregor Mendel would prefer to practice scientific approaches other than molecular genetics such as population genetics, comparative genomics, or the emerging field of evo-chromo. In this perspective, I review aspects of these fields that might have attracted or perplexed a 21st century Mendel.

Investigating the genetic architecture of eye colour in a Canadian cohort

Eye color is highly variable in populations with European ancestry, ranging from low to high quantities of melanin in the iris. Polymorphisms in the HERC2/OCA2 locus have the largest effect on eye color in these populations, although other genomic regions also influence eye color. We performed genome-wide association studies of eye color in a Canadian cohort of European ancestry (N = 5,641) and investigated candidate causal variants. We uncovered several candidate causal signals in the HERC2/OCA2 region, whereas other loci likely harbor a single causal signal. We observed colocalization of eye color signals with the expression or methylation profiles of cultured primary melanocytes. Genetic correlations of eye and hair color suggest high genome-wide pleiotropy, but locus-level differences in the genetic architecture of both traits. Overall, we provide a better picture of the polymorphisms underpinning eye color variation, which may be a consequence of specific molecular processes in the iris melanocytes.

•

Genome-wide association studies of eye color in 5,641 participants

•

Multiple independent candidate causal variants were identified across HERC2/OCA2

•

Single candidate causal variants observed on or near IRF4, SLC24A4, TYR, and TYRP1

•

Colocalization of eye color signals with expression and methylation profiles

Genetic Connections and Convergent Evolution of Tropical Indigenous Peoples in Asia

Tropical indigenous peoples in Asia (TIA) attract much attention for their unique appearance, whereas their genetic history and adaptive evolution remain mysteries. We conducted a comprehensive study to characterize the genetic distinction and connection of broad geographical TIAs. Despite the diverse genetic makeup and large interarea genetic differentiation between the TIA groups, we identified a basal Asian ancestry (bASN) specifically shared by these populations. The bASN ancestry was relatively enriched in ancient Asian human genomes dated as early as ∼50,000 years before the present and diminished in more recent history. Notably, the bASN ancestry is unlikely to be derived from archaic hominins. Instead, we suggest it may be better modeled as a survived lineage of the initial peopling of Asia. Shared adaptations inherited from the ancient Asian ancestry were detected among the TIA groups (e.g., LIMS1 for hair morphology, and COL24A1 for bone formation), and they are enriched in neurological functions either at an identical locus (e.g., NKAIN3), or different loci in an identical gene (e.g., TENM4). The bASN ancestry could also have formed the substrate of the genetic architecture of the dark pigmentation observed in the TIA peoples. We hypothesize that phenotypic convergence of the dark pigmentation in TIAs could have resulted from parallel (e.g., DDB1/DAK) or genetic convergence driven by admixture (e.g., MTHFD1 and RAD18), new mutations (e.g., STK11), or notably purifying selection (e.g., MC1R). Our results provide new insights into the initial peopling of Asia and an advanced understanding of the phenotypic convergence of the TIA peoples.

Advances in integrative African genomics

There has been a rapid increase in human genome sequencing in the past two decades, resulting in the identification of millions of previously unknown genetic variants. However, African populations are under-represented in sequencing efforts. Additional sequencing from diverse African populations and the construction of African-specific reference genomes is needed to better characterize the full spectrum of variation in humans. However, sequencing alone is insufficient to address the molecular and cellular mechanisms underlying variable phenotypes and disease risks. Determining functional consequences of genetic variation using multi-omics approaches is a fundamental post-genomic challenge. We discuss approaches to close the knowledge gaps about African genomic diversity and review advances in African integrative genomic studies and their implications for precision medicine.

Interactions between PDA-associated polymorphisms and genetic ancestry alter ductus arteriosus gene expression

BACKGROUND

DNA polymorphisms in PTGIS and TFAP2B have been identified as risk factors for patent ductus arteriosus (PDA) in a population composed of preterm infants with European genetic ancestry but not in more genetically diverse populations.

GOAL

To determine if the effects of TFAP2B and PTGIS polymorphisms on ductus arteriosus (DA) gene expression differ based on genetic ancestry.

METHODS

DA from 273 human second trimester fetuses were genotyped for TFAP2B and PTGIS polymorphisms and for polymorphisms distributing along genetic ancestry lines. RT-PCR was used to measure the RNA expression of 49 candidate genes involved with DA closure.

RESULTS

Seventeen percent of the DA analyzed were of European ancestry. In multivariable regression analyses we found consistent associations between four PDA-related TFAP2B polymorphisms (rs2817399(A), rs987237(G), rs760900(C), and rs2817416(C)) and expression of the following genes: EPAS1, CACNB2, ECE1, KCNA2, ATP2A3, EDNRA, EDNRB, BMP9, and BMP10, and between the PTGIS haplotype rs493694(G)/rs693649(A) and PTGIS and NOS3. These changes only occurred in DA with European ancestry. No consistent positive or negative associations were found among DA samples unless an interaction between the polymorphisms and genetic ancestry was taken into account.

CONCLUSION

PTGIS and TFAP2B polymorphisms were associated with consistent changes in DA gene expression when present in fetuses with European ancestry.

IMPACT

DNA polymorphisms in PTGIS and TFAP2B have been identified as risk factors for patent ductus arteriosus (PDA) in a population composed primarily of preterm infants with European genetic ancestry but not in more genetically diverse populations. The same PTGIS and TFAP2B polymorphisms are associated with changes in ductus gene expression when present in ductus from fetuses with European genetic ancestry. No consistent associations with gene expression can be found unless an interaction between the polymorphisms and genetic ancestry is taken into account.

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.

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.

Association of pigmentation related-genes polymorphisms and geographic environmental variables in the Chinese population

Background

Human skin color is highly heritable and one of the most variable phenotypic traits. However, the genetic causes and environmental selective pressures underlying this phenotypic variation have remained largely unknown. To investigate whether the pigmentation related-genes polymorphisms are associated with the geographic environmental variables. We selected randomly 795 healthy individuals from eight ethnic groups in nine provinces in China. Six single nucleotide polymorphisms (SNPs) of SLC45A2 and TYR were genotyped using Agena MassARRAY. The Chi-square test and Spearman correlation analysis were used to compare the frequency distribution of genotypes among different ethnic groups and evaluate the relationship between SNP genetic diversity and environmental variables, respectively.

Results

The results indicated that rs28777 and rs183671 (SLC45A2) and rs1042602 (TYR) genotype frequency distributions were significantly different between the Xinjiang-Uighur and other ethnic groups (P < 0.05). Spearman correlation analysis found that rs28777-A (r = − 0.090, P = 0.011), rs183671-G (r = − 0.105, P = 0.003), rs1042602-A (r = − 0.108, P = 0.002), rs1126809-A (r = − 0.151, P < 0.001) allele frequencies were negatively correlated with the longitude; rs183671-G (r = 0.151), rs1042602-A (r = 0.157) and rs1126809-A (r = 0.138) allele frequencies were positively associated with the latitude (P < 0.001); rs183671-G (r = 0.116, P = 0.001), rs1042602-A (r = 0.105, P = 0.003) and rs1126809-A (r = 0.070, P = 0.048) allele frequencies were positively correlated with the sunshine hours; rs183671-G (r = − 0.076, P = 0.033), rs1042602-A (r = − 0.079, P = 0.027) and rs1126809-A (r = − 0.076, P = 0.031) were negatively correlated with the annual average temperature.

Conclusions

Our results confirmed the idea that environmental factors have been an important selective pressure upon pigmentation related gene polymorphisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-021-00189-7.

The evolution of human skin pigmentation involved the interactions of genetic, environmental, and cultural variables

The primary biological role of human skin pigmentation is as a mediator of penetration of ultraviolet radiation (UVR) into the deep layers of skin and the cutaneous circulation. Since the origin of Homo sapiens, dark, protective constitutive pigmentation and strong tanning abilities have been favored under conditions of high UVR and represent the baseline condition for modern humans. The evolution of partly depigmented skin and variable tanning abilities has occurred multiple times in prehistory, as populations have dispersed into environments with lower and more seasonal UVR regimes, with unique complements of genes and cultural practices. The evolution of extremes of dark pigmentation and depigmentation has been rare and occurred only under conditions of extremely high or low environmental UVR, promoted by positive selection on variant pigmentation genes followed by limited gene flow. Over time, the evolution of human skin pigmentation has been influenced by the nature and course of human dispersals and modifications of cultural practices, which have modified the nature and actions of skin pigmentation genes. Throughout most of prehistory and history, the evolution of human skin pigmentation has been a contingent and non-deterministic process.

The constraints of racialization: How classification and valuation hinder scientific research on human variation

Human biological variation has historically been studied through the lens of racialization. Despite a general shift away from the use of overt racial terminologies, the underlying racialized frameworks used to describe and understand human variation still remain. Even in relatively recent anthropological and biomedical work, we can observe clear manifestations of such racial thinking. This paper shows how classification and valuation are two specific processes which facilitate racialization and hinder attempts to move beyond such frameworks. The bias induced by classification distorts descriptions of phenotypic variation in a way that erroneously portrays European populations as more variable than others. Implicit valuation occurs in tandem with classification and produces narratives of superiority/inferiority for certain phenotypic variants without an objective biological basis. The bias of racialization is a persistent impediment stemming from the inheritance of scientific knowledge developed under explicitly racial paradigms. It is also an internalized cognitive distortion cultivated through socialization in a world where racialization is inescapable. Though undeniably challenging, this does not present an insurmountable barrier, and this bias can be mitigated through the critical evaluation of past work, the active inclusion of marginalized perspectives, and the direct confrontation of institutional structures enforcing racialized paradigms.

Analysis of Skin Pigmentation and Genetic Ancestry in Three Subpopulations from Pakistan: Punjabi, Pashtun, and Baloch

Skin pigmentation is one of the most prominent and variable phenotypes in humans. We compared the alleles of 163 SNPs and indels from the Human Pigmentation (HuPi) AmpliSeq™ Custom panel, and biogeographic ancestry with the quantitative skin pigmentation levels on the upper arm, lower arm, and forehead of 299 Pakistani individuals from three subpopulations: Baloch, Pashtun, and Punjabi. The biogeographic ancestry of each individual was estimated using the Precision ID Ancestry Panel. All individuals were mainly of mixed South-Central Asian and European ancestry. However, the Baloch individuals also had an average proportion of Sub-Saharan African ancestry of approximately 10%, whereas it was <1% in the Punjabi and Pashtun individuals. The pairwise genetic distances between the Pashtun, Punjabi, and Baloch subpopulations based on the ancestry markers were statistically significantly different. Individuals from the Pashtun subpopulation had statistically significantly lower skin pigmentation than individuals from the Punjabi and Baloch subpopulations (p < 0.05). The proportions of European and Sub-Saharan African ancestry and five SNPs (rs1042602, rs10831496, rs1426654, rs16891982, and rs12913832) were statistically significantly associated with skin pigmentation at either the upper arm, lower arm or forehead in the Pakistani population after correction for multiple testing (p < 10⁻³). A model based on four of these SNPs (rs1426654, rs1042602, rs16891982, and rs12913832) explained 33% of the upper arm skin pigmentation. The four SNPs and the proportions of European and Sub-Saharan African ancestry explained 37% of the upper arm skin pigmentation. Our results indicate that the four likely causative SNPs, rs1426654, rs1042602, rs16891982, and rs12913832 located in SLC24A5, TYR, SLC45A2, and HERC2, respectively, are essential for skin color variation in the admixed Pakistani subpopulations.

The evolution of skin pigmentation-associated variation in West Eurasia

Some of the genes responsible for the evolution of light skin pigmentation in Europeans show signals of positive selection in present-day populations. Recently, genome-wide association studies have highlighted the highly polygenic nature of skin pigmentation. It is unclear whether selection has operated on all of these genetic variants or just a subset. By studying variation in over a thousand ancient genomes from West Eurasia covering 40,000 y, we are able to study both the aggregate behavior of pigmentation-associated variants and the evolutionary history of individual variants. We find that the evolution of light skin pigmentation in Europeans was driven by frequency changes in a relatively small fraction of the genetic variants that are associated with variation in the trait today.

Skin pigmentation is a classic example of a polygenic trait that has experienced directional selection in humans. Genome-wide association studies have identified well over a hundred pigmentation-associated loci, and genomic scans in present-day and ancient populations have identified selective sweeps for a small number of light pigmentation-associated alleles in Europeans. It is unclear whether selection has operated on all of the genetic variation associated with skin pigmentation as opposed to just a small number of large-effect variants. Here, we address this question using ancient DNA from 1,158 individuals from West Eurasia covering a period of 40,000 y combined with genome-wide association summary statistics from the UK Biobank. We find a robust signal of directional selection in ancient West Eurasians on 170 skin pigmentation-associated variants ascertained in the UK Biobank. However, we also show that this signal is driven by a limited number of large-effect variants. Consistent with this observation, we find that a polygenic selection test in present-day populations fails to detect selection with the full set of variants. Our data allow us to disentangle the effects of admixture and selection. Most notably, a large-effect variant at SLC24A5 was introduced to Western Europe by migrations of Neolithic farming populations but continued to be under selection post-admixture. This study shows that the response to selection for light skin pigmentation in West Eurasia was driven by a relatively small proportion of the variants that are associated with present-day phenotypic variation.

Prediction of eye and hair pigmentation phenotypes using the HIrisPlex system in a Brazilian admixed population sample

Human pigmentation is a complex trait, probably involving more than 100 genes. Predicting phenotypes using SNPs present in those genes is important for forensic purpose. For this, the HIrisPlex tool was developed for eye and hair color prediction, with both models achieving high accuracy among Europeans. Its evaluation in admixed populations is important, since they present a higher frequency of intermediate phenotypes, and HIrisPlex has demonstrated limitations in such predictions; therefore, the performance of this tool may be impaired in such populations. Here, we evaluate the set of 24 markers from the HIrisPlex system in 328 individuals from Ribeirão Preto (SP) region, predicting eye and hair color and comparing the predictions with their real phenotypes. We used the HaloPlex Target Enrichment System and MiSeq Personal Sequencer platform for massively parallel sequencing. The prediction of eye and hair color was accomplished by the HIrisPlex online tool, using the default prediction settings. Ancestry was estimated using the SNPforID 34-plex to observe if and how an individual's ancestry background would affect predictions in this admixed sample. Our sample presented major European ancestry (70.5%), followed by African (21.1%) and Native American/East Asian (8.4%). HIrisPlex presented an overall sensitivity of 0.691 for hair color prediction, with sensitivities ranging from 0.547 to 0.782. The lowest sensitivity was observed for individuals with black hair, who present a reduced European contribution (48.4%). For eye color prediction, the overall sensitivity was 0.741, with sensitivities higher than 0.85 for blue and brown eyes, although it failed in predicting intermediate eye color. Such struggle in predicting this phenotype category is in accordance with what has been seen in previous studies involving HIrisPlex. Individuals with brown eye color are more admixed, with European ancestry decreasing to 62.6%; notwithstanding that, sensitivity for brown eyes was almost 100%. Overall sensitivity increases to 0.791 when a 0.7 threshold is set, though 12.5% of the individuals become undefined. When combining eye and hair prediction, hit rates between 51.3 and 68.9% were achieved. Despite the difficulties with intermediate phenotypes, we have shown that HIrisPlex results can be very helpful when interpreted with caution.

Skin deep: The decoupling of genetic admixture levels from phenotypes that differed between source populations

OBJECTIVES

In genetic admixture processes, source groups for an admixed population possess distinct patterns of genotype and phenotype at the onset of admixture. Particularly in the context of recent and ongoing admixture, such differences are sometimes taken to serve as markers of ancestry for individuals-that is, phenotypes initially associated with the ancestral background in one source population are assumed to continue to reflect ancestry in that population. Such phenotypes might possess ongoing significance in social categorizations of individuals, owing in part to perceived continuing correlations with ancestry. However, genotypes or phenotypes initially associated with ancestry in one specific source population have been seen to decouple from overall admixture levels, so that they no longer serve as proxies for genetic ancestry. Here, we aim to develop an understanding of the joint dynamics of admixture levels and phenotype distributions in an admixed population.

METHODS

We devise a mechanistic model, consisting of an admixture model, a quantitative trait model, and a mating model. We analyze the behavior of the mechanistic model in relation to the model parameters.

RESULTS

We find that it is possible for the decoupling of genetic ancestry and phenotype to proceed quickly, and that it occurs faster if the phenotype is driven by fewer loci. Positive assortative mating attenuates the process of dissociation relative to a scenario in which mating is random with respect to genetic admixture and with respect to phenotype.

CONCLUSIONS

The mechanistic framework suggests that in an admixed population, a trait that initially differed between source populations might serve as a reliable proxy for ancestry for only a short time, especially if the trait is determined by few loci. It follows that a social categorization based on such a trait is increasingly uninformative about genetic ancestry and about other traits that differed between source populations at the onset of admixture.

Human pegivirus 1 in Cabo Verde: prevalence and genotypic distribution among HIV-infected individuals

Human pegivirus 1 (HPgV-1) belongs to the genus Pegivirus, family Flaviviridae, and until now has been considered a non-pathogenic agent, despite being considered a risk factor for non-Hodgkin lymphoma. However, a beneficial impact of HPgV-1 on HIV disease progression has been extensively reported. Given the high prevalence of HIV in sub-Saharan Africa and the scarcity of epidemiological data for many countries of West Africa, we conducted the first study of HPgV-1 in HIV-infected individuals from Cabo Verde. To obtain new data regarding prevalence and genetic diversity of HPgV-1 in Africa, serum samples from 102 HIV-infected Cabo Verdeans were tested for the presence of viral RNA, and the circulating genotypes were identified by sequencing of the 5' untranslated region. HPgV-1 RNA was detected in 19.6% (20/102) of the samples. In 72.2% (13/18) of the samples, the virus was identified as genotype 2 (11/13 subtype 2a and 2/13 subtype 2b), and in 27.8% (5/18), it was identified as genotype 1. The estimated substitution rate of HPgV-1 genotype 2 was 5.76 × 10^-4, and Bayesian analysis indicated the existence of inner clusters within subtypes 2a and 2b. The prevalence of HPgV-1 viremia in Cabo Verde agrees with that reported previously in Africa. Genotypes 1 and 2 cocirculate, with genotype 2 being more common, and HIV/HPgV-1 coinfection was not associated with higher CD4 T cell counts in the studied population. This finding contributes for the expansion of the pegivirus research agenda in African countries.

Human genetic admixture

Throughout human history, large-scale migrations have facilitated the formation of populations with ancestry from multiple previously separated populations. This process leads to subsequent shuffling of genetic ancestry through recombination, producing variation in ancestry between populations, among individuals in a population, and along the genome within an individual. Recent methodological and empirical developments have elucidated the genomic signatures of this admixture process, bringing previously understudied admixed populations to the forefront of population and medical genetics. Under this theme, we present a collection of recent PLOS Genetics publications that exemplify recent progress in human genetic admixture studies, and we discuss potential areas for future work.

Genetic loci associated with skin pigmentation in African Americans and their effects on vitamin D deficiency

A recent genome-wide association study (GWAS) in African descent populations identified novel loci associated with skin pigmentation. However, how genomic variations affect skin pigmentation and how these skin pigmentation gene variants affect serum 25(OH) vitamin D variation has not been explored in African Americans (AAs). In order to further understand genetic factors that affect human skin pigmentation and serum 25(OH)D variation, we performed a GWAS for skin pigmentation with 395 AAs and a replication study with 681 AAs. Then, we tested if the identified variants are associated with serum 25(OH) D concentrations in a subset of AAs (n = 591). Skin pigmentation, Melanin Index (M-Index), was measured using a narrow-band reflectometer. Multiple regression analysis was performed to identify variants associated with M-Index and to assess their role in serum 25(OH)D variation adjusting for population stratification and relevant confounding variables. A variant near the SLC24A5 gene (rs2675345) showed the strongest signal of association with M-Index (P = 4.0 x 10-30 in the pooled dataset). Variants in SLC24A5, SLC45A2 and OCA2 together account for a large proportion of skin pigmentation variance (11%). The effects of these variants on M-Index was modified by sex (P for interaction = 0.009). However, West African Ancestry (WAA) also accounts for a large proportion of M-Index variance (23%). M-Index also varies among AAs with high WAA and high Genetic Score calculated from top variants associated with M-Index, suggesting that other unknown genomic factors related to WAA are likely contributing to skin pigmentation variation. M-Index was not associated with serum 25(OH)D concentrations, but the Genetic Score was significantly associated with vitamin D deficiency (serum 25(OH)D levels less than 12 ng/mL) (OR, 1.30; 95% CI, 1.04-1.64). The findings support the hypothesis suggesting that skin pigmentation evolved responding to increased demand for subcutaneous vitamin D synthesis in high latitude environments.

Rapid adaptation to malaria facilitated by admixture in the human population of Cabo Verde

Humans have undergone large migrations over the past hundreds to thousands of years, exposing ourselves to new environments and selective pressures. Yet, evidence of ongoing or recent selection in humans is difficult to detect. Many of these migrations also resulted in gene flow between previously separated populations. These recently admixed populations provide unique opportunities to study rapid evolution in humans. Developing methods based on distributions of local ancestry, we demonstrate that this sort of genetic exchange has facilitated detectable adaptation to a malaria parasite in the admixed population of Cabo Verde within the last ~20 generations. We estimate that the selection coefficient is approximately 0.08, one of the highest inferred in humans. Notably, we show that this strong selection at a single locus has likely affected patterns of ancestry genome-wide, potentially biasing demographic inference. Our study provides evidence of adaptation in a human population on historical timescales.

The color of normal: How a Eurocentric focus erases pigmentation complexity

OBJECTIVES

Skin pigmentation is both a highly variable and highly visible human phenotypic trait. Investigations into the biology and origins of this variation have been the focus of research in the fields of dermatology, anthropology, and forensic science, among others. This manuscript explores how much of what we know about the biology, genetics, and evolutionary origins of pigmentation has been strongly influenced by investigations and applications that focus on lighter skin.

METHODS

I reviewed literature from the fields of dermatology, anthropology and evolutionary genetics, and forensic science to assess how perceptions of lighter skin as the "normal" state in humans can shape the ways that knowledge is gathered and applied in these fields.

RESULTS

This normalization of lighter skin has impacted common tools used in dermatology and shaped the framework of dermatological education. A strong Eurocentric bias has shaped our understanding of the genetic architecture of pigmentary traits, which influences the ways in we understand the evolutionary processes leading to modern pigmentation diversity. Finally, I discuss how these biases in pigmentation genetics work in combination with phenotypic systems that privilege predicting lighter pigmentation variation to impede accurate prediction of intermediate phenotypes, particularly in individuals with ancestry from multiple populations. This can lead to a disproportionate targeting of already over-policed populations with darker skin.

CONCLUSIONS

Potential changes to how we conceptualize clinical and basic pigmentation research may help to reduce existing health disparities and improve understanding of pigmentation genetic architecture and how this knowledge is applied in forensic contexts.

Anthropological genetics perspectives on the transatlantic slave trade

During the Trans-Atlantic Slave Trade (TAST), around twelve million Africans were enslaved and forcibly moved from Africa to the Americas and Europe, durably influencing the genetic and cultural landscape of a large part of humanity since the 15th century. Following historians, archaeologists, and anthropologists, population geneticists have, since the 1950's mainly, extensively investigated the genetic diversity of populations on both sides of the Atlantic. These studies shed new lights into the largely unknown genetic origins of numerous enslaved-African descendant communities in the Americas, by inferring their genetic relationships with extant African, European, and Native American populations. Furthermore, exploring genome-wide data with novel statistical and bioinformatics methods, population geneticists have been increasingly able to infer the last 500 years of admixture histories of these populations. These inferences have highlighted the diversity of histories experienced by enslaved-African descendants, and the complex influences of socioeconomic, political, and historical contexts on human genetic diversity patterns during and after the slave trade. Finally, the recent advances of paleogenomics unveiled crucial aspects of the life and health of the first generation of enslaved-Africans in the Americas. Altogether, human population genetics approaches in the genomic and paleogenomic era need to be coupled with history, archaeology, anthropology, and demography in interdisciplinary research, to reconstruct the multifaceted and largely unknown history of the TAST and its influence on human biological and cultural diversities today. Here, we review anthropological genomics studies published over the past 15 years and focusing on the history of enslaved-African descendant populations in the Americas.

Ecological adaptation in Atlantic herring is associated with large shifts in allele frequencies at hundreds of loci

Atlantic herring is widespread in North Atlantic and adjacent waters and is one of the most abundant vertebrates on earth. This species is well suited to explore genetic adaptation due to minute genetic differentiation at selectively neutral loci. Here, we report hundreds of loci underlying ecological adaptation to different geographic areas and spawning conditions. Four of these represent megabase inversions confirmed by long read sequencing. The genetic architecture underlying ecological adaptation in herring deviates from expectation under a classical infinitesimal model for complex traits because of large shifts in allele frequencies at hundreds of loci under selection.

GWAS Analysis of 17,019 Korean Women Identifies the Variants Associated with Facial Pigmented Spots

Variation in skin pigmentation can be affected by both environmental factors and intrinsic factors such as age, gender, and genetic variation. Recent GWASs revealed that genetic variants of genes functionally related to a pigmentation pathway were associated with skin pigmentary traits. However, these GWASs focused on populations with European ancestry, and only a few studies have been performed on Asian populations, limiting our understanding of the genetic basis of skin pigmentary traits in Asians. To evaluate the genetic variants associated with facial pigmented spots, we conducted a GWAS analysis of objectively measured facial pigmented spots in 17,019 Korean women. This large-scale GWAS identified several genomic loci that were significantly associated with facial pigmented spots (five previously reported loci and two previously unreported loci, to our knowledge), which were detected by UV light: BNC2 at 9p22 (rs16935073; P-value = 2.11 × 10^-46), PPARGC1B at 5q32 (rs32579; P-value = 9.04 × 10^-42), 10q26 (rs11198112; P-value = 9.66 × 10^-38), MC1R at 16q24 (rs2228479; P-value = 6.62 × 10^-21), lnc01877 at 2q33 (rs12693889; P-value = 1.59 × 10^-11), CDKN2B-AS1 at 9p21 (rs643319; P-value = 7.76 × 10^-9), and MFSD12 at 19p13 (rs2240751; P-value = 9.70 × 10^-9). Further functional characterization of the candidate genes needs to be done to fully evaluate their contribution to facial pigmented spots.

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

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

Skin pigmentation and genetic variants in an admixed Brazilian population of primarily European ancestry

Although many genes have been shown to be associated with human pigmentary traits and forensic prediction assays exist (e.g. HIrisPlex-S), the genetic knowledge about skin colour remains incomplete. The highly admixed Brazilian population is an interesting study population for investigation of the complex genotype-phenotype architecture of human skin colour because of its large variation. Here, we compared variants in 22 pigmentary genes with quantitative skin pigmentation levels on the buttock, arm, and forehead areas of 266 genetically admixed Brazilian individuals. The genetic ancestry of each individual was estimated by typing 46 AIM-InDels. The mean proportion of genetic ancestry was 68.8% European, 20.8% Sub-Saharan African, and 10.4% Native American. A high correlation (adjusted R² = 0.65, p < 0.05) was observed between nine SNPs and quantitative skin pigmentation using multiple linear regression analysis. The correlations were notably smaller between skin pigmentation and biogeographic ancestry (adjusted R² = 0.45, p < 0.05), or markers in the leading forensic skin colour prediction system, the HIrisPlex-S (adjusted R² = 0.54, p < 0.05). Four of the nine SNPs, OCA2 rs1448484 (rank 2), APBA2 rs4424881 (rank 4), MFSD12 rs10424065 (rank 8), and TYRP1 1408799 (rank 9) were not investigated as part of the HIrisPlex-S selection process, and therefore not included in the HIrisPlex-S model. Our results indicate that these SNPs account for a substantial part of the skin colour variation in individuals of admixed ancestry. Hence, we suggest that these SNPs are considered when developing future skin colour prediction models.

Identification of SNPs involved in Zn and Cu accumulation in the Pacific oyster (Crassostrea gigas) by genome-wide association analysis

Oysters accumulate high concentrations of zinc (Zn) and copper (Cu), which can be transferred to human due to sea food consumption. Breeding new oyster varieties with low Zn and Cu accumulations is one important way to improve food safety. However, the genetic basis for metal accumulation in mollusks is not well understood. To address this issue, oysters collected in the field were used for genome-wide association study (GWAS) and then the identified genes were used for mRNA expressions analysis in laboratory. First, GWAS were conducted for Zn and Cu accumulation in 288 wild Pacific oysters (Crassostrea gigas) farmed in the same ocean environment. The oysters did not show obvious population structure or kinship but exhibited 8.43- and 10.0- fold changes of Zn and Cu contents respectively. GWAS have identified 11 and 12 single nucleotide polymorphisms (SNPs) associated with Zn and Cu, respectively, as well as 16 genes, which were Zn-containing proteins or participated in caveolae-dependent endocytosis. Second, the mRNA expressions of these 16 genes were observed under Zn and Cu exposure. After 9 days of Zn exposure, Zn contents increased 3.1-fold, while the mRNA expression of cell number regulator 3 increased 1.65-fold. Under 9 days of Cu exposure, Cu contents increased 1.97-fold, while the mRNA expression of caveolin-1 decreased 0.61-fold. These provide the evidence for their roles in regulating physiological levels of these two metals. The findings advance our understanding of the genetic basis of Zn and Cu accumulation in mollusks, which can be useful for breeding new, less toxic varieties of oysters.

Genome-wide Analyses Identifies Known and New Markers Responsible of Chicken Plumage Color

Simple Summary

In order to assess sources of variation related to Polverara breed plumage color (black vs. white), we carried out genome-wide analyses to identify the genomic regions involved in this trait. The present work has revealed new candidate genes involved in the phenotypic variability in local chicken populations. These results also contribute insights into the genetic basis for plumage color in poultry, and confirm the great complexity of the mechanisms that control this trait.

Abstract

Through the development of the high-throughput genotyping arrays, molecular markers and genes related to phenotypic traits have been identified in livestock species. In poultry, plumage color is an important qualitative trait that can be used as phenotypic marker for breed identification. In order to assess sources of genetic variation related to the Polverara chicken breed plumage colour (black vs. white), we carried out a genome-wide association study (GWAS) and a genome-wide fixation index (F_ST) scan to uncover the genomic regions involved. A total of 37 animals (17 white and 20 black) were genotyped with the Affymetrix 600 K Chicken single nucleotide polymorphism (SNP) Array. The combination of results from GWAS and F_ST revealed a total of 40 significant markers distributed on GGA 01, 03, 08, 12 and 21, and located within or near known genes. In addition to the well-known TYR, other candidate genes have been identified in this study, such as GRM5, RAB38 and NOTCH2. All these genes could explain the difference between the two Polverara breeds. Therefore, this study provides the basis for further investigation of the genetic mechanisms involved in plumage color in chicken.

Darwinian Positive Selection on the Pleiotropic Effects of KITLG Explain Skin Pigmentation and Winter Temperature Adaptation in Eurasians

Human skin color diversity is considered an adaptation to environmental conditions such as UV radiation. Investigations into the genetic bases of such adaptation have identified a group of pigmentation genes contributing to skin color diversity in African and non-African populations. Here, we present a population analysis of the pigmentation gene KITLG with previously reported signal of Darwinian positive selection in both European and East Asian populations. We demonstrated that there had been recurrent selective events in the upstream and the downstream regions of KITLG in Eurasian populations. More importantly, besides the expected selection on the KITLG variants favoring light skin in coping with the weak UV radiation at high latitude, we observed a KITLG variant showing adaptation to winter temperature. In particular, compared with UV radiation, winter temperature showed a much stronger correlation with the prevalence of the presumably adaptive KITLG allele in Asian populations. This observation was further supported by the in vitro functional test at low temperature. Consequently, the pleiotropic effects of KITLG, that is, pigmentation and thermogenesis were both targeted by natural selection that acted on different KITLG sequence variants, contributing to the adaptation of Eurasians to both UV radiation and winter temperature at high latitude areas.

The Evolutionary History of Human Skin Pigmentation

Skin pigmentation is a complex, conspicuous, highly variable human trait that exhibits a remarkable correlation with latitude. The evolutionary history and genetic basis of skin color variation has been the subject of intense research in the last years. This article reviews the major hypotheses explaining skin color diversity and explores the implications of recent findings about the genes associated with skin pigmentation for understanding the evolutionary forces that have shaped the current patterns of skin color variation. A major aspect of these findings is that the genetic basis of skin color is less simple than previously thought and that geographic variation in skin pigmentation was influenced by the concerted action of different types of natural selection, rather than just by selective sweeps in a few key genes.

Meta-analysis of GWA studies provides new insights on the genetic architecture of skin pigmentation in recently admixed populations

BACKGROUND

Association studies in recently admixed populations are extremely useful to identify the genetic architecture of pigmentation, due to their high genotypic and phenotypic variation. However, to date only four Genome-Wide Association Studies (GWAS) have been carried out in these populations.

RESULTS

We present a GWAS of skin pigmentation in an admixed sample from Cuba (N = 762). Additionally, we conducted a meta-analysis including the Cuban sample, and admixed samples from Cape Verde, Puerto Rico and African-Americans from San Francisco. This meta-analysis is one of the largest efforts so far to characterize the genetic basis of skin pigmentation in admixed populations (N = 2,104). We identified five genome-wide significant regions in the meta-analysis, and explored if the markers observed in these regions are associated with the expression of relevant pigmentary genes in human melanocyte cultures. In three of the regions identified in the meta-analysis (SLC24A5, SLC45A2, and GRM5/TYR), the association seems to be driven by non-synonymous variants (rs1426654, rs16891982, and rs1042602, respectively). The rs16891982 polymorphism is strongly associated with the expression of the SLC45A2 gene. In the GRM5/TYR region, in addition to the rs1042602 non-synonymous SNP located on the TYR gene, variants located in the nearby GRM5 gene have an independent effect on pigmentation, possibly through regulation of gene expression of the TYR gene. We also replicated an association recently described near the MFSD12 gene on chromosome 19 (lead variant rs112332856). Additionally, our analyses support the presence of multiple signals in the OCA2/HERC2/APBA2 region on chromosome 15. A clear causal candidate is the HERC2 intronic variant rs12913832, which has a profound influence on OCA2 expression. This variant has pleiotropic effects on eye, hair, and skin pigmentation. However, conditional and haplotype-based analyses indicate the presence of other variants with independent effects on melanin levels in OCA2 and APBA2. Finally, a follow-up of genome-wide signals identified in a recent GWAS for tanning response indicates that there is a substantial overlap in the genetic factors influencing skin pigmentation and tanning response.

CONCLUSIONS

Our meta-analysis of skin pigmentation GWAS in recently admixed populations provides new insights about the genetic architecture of this complex trait.

A Genome-Wide Association Study of Skin and Iris Pigmentation among Individuals of South Asian Ancestry

South Asia has a complex history of migrations and is characterized by substantial pigmentary and genetic diversity. For this reason, it is an ideal region to study the genetic architecture of normal pigmentation variation. Here, we present a meta-analysis of two genome-wide association studies (GWASs) of skin pigmentation using skin reflectance (M-index) as a quantitative phenotype. The meta-analysis includes a sample of individuals of South Asian descent living in Canada (N = 348), and a sample of individuals from two caste and four tribal groups from West Maharashtra, India (N = 480). We also present the first GWAS of iris color in South Asian populations. This GWAS was based on quantitative measures of iris color obtained from high-resolution iris pictures. We identified genome-wide significant associations of variants within the well-known gene SLC24A5, including the nonsynonymous rs1426654 polymorphism, with both skin pigmentation and iris color, highlighting the pleiotropic effects of this gene on pigmentation. Variants in the HERC2 gene (e.g., rs12913832) were also associated with iris color and iris heterochromia. Our study emphasizes the usefulness of quantitative methods to study iris color variation. We also identified novel genome-wide significant associations with skin pigmentation and iris color, but we could not replicate these associations due to the lack of independent samples. It will be critical to expand the number of studies in South Asian populations in order to better understand the genetic variation driving the diversity of skin pigmentation and iris color observed in this region.

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.