Association study confirms the role of twoOCA2polymorphisms in normal skin pigmentation variation in East Asian populations

Novel compound heterozygous mutations in OCA2 gene were identified in a Chinese family with oculocutaneous albinism

BACKGROUND

Oculocutaneous albinism (OCA) is a group of rare autosomal recessive disorders characterized by clinical genetic heterogeneity. OCA type II (OMIM: 203200) is the most common subtype among African and African Americans, primarily caused by pathogenic variants in the OCA2 (HGNC ID: 8101) gene. In this study, we presented a Chinese family with OCA and reported two novel variants in the OCA2 gene.

METHODS

Whole-exome sequencing (WES) was performed to identify pathogenic variants in the proband. The candidate variants were subsequently validated using Sanger sequencing and QPCR assay. Additionally, bioinformatics analyses were employed to predict the deleteriousness and conservation of the identified mutations.

RESULTS

In the 16-year-old male proband, two novel compound heterozygous OCA2 variants, NM_000275.3: c.1640T>G (NP_000266.2: p.L547R) and an exons 10-19 deletion variant, were identified. Meanwhile, a reported heterozygous variant c.1441G>A/p.A481T (NM_000275.3, NP_000266.2) in the OCA2 gene was also found in the proband. Sanger sequencing confirmed that the two variants c.1441G>A/p.A481T and c.1640T>G/p.L547R were inherited from his father. Moreover, qPCR assay revealed that the exons 10-19 deletion was inherited from the mother, his sister also carried this variant. Fortunately, the variant was not detected in the amniotic fluid of the proband's sister. Multiple online bioinformatics tools predicted the variant c.1640T>G to be damaging, leading to the replacement of a highly conserved leucine with an arginine. The gross exon 10-19 deletion in the OCA2 gene resulted in a truncated, non-functional protein losing the 3-9 transmembrane α-helices domains. According to the American College of Medical Genetics and Genomics classification, these three variants in the OCA2 gene were evaluated as likely pathogenic.

CONCLUSION

This study has identified two novel compound variants in the OCA2 gene and a previously reported variant in a Chinese family with OCA. By expanding the mutation spectrum of the OCA2 gene, our findings contribute to a better understanding of the genetic basis of OCA.

Human genetic history on the Tibetan Plateau in the past 5100 years

Using genome-wide data of 89 ancient individuals dated to 5100 to 100 years before the present (B.P.) from 29 sites across the Tibetan Plateau, we found plateau-specific ancestry across plateau populations, with substantial genetic structure indicating high differentiation before 2500 B.P. Northeastern plateau populations rapidly showed admixture associated with millet farmers by 4700 B.P. in the Gonghe Basin. High genetic similarity on the southern and southwestern plateau showed population expansion along the Yarlung Tsangpo River since 3400 years ago. Central and southeastern plateau populations revealed extensive genetic admixture within the plateau historically, with substantial ancestry related to that found in southern and southwestern plateau populations. Over the past ~700 years, substantial gene flow from lowland East Asia further shaped the genetic landscape of present-day plateau populations. The high-altitude adaptive EPAS1 allele was found in plateau populations as early as in a 5100-year-old individual and showed a sharp increase over the past 2800 years.

Association between Variants in the OCA2-HERC2 Region and Blue Eye Colour in HERC2 rs12913832 AA and AG Individuals

The OCA2-HERC2 region is strongly associated with human pigmentation, especially eye colour. The HERC2 SNP rs12913832 is currently the best-known predictor for blue and brown eye colour. However, in a previous study we found that 43 of 166 Norwegians with the brown eye colour genotype rs12913832:AA or AG, did not have the expected brown eye colour. In this study, we carried out massively parallel sequencing of a ~500 kbp HERC2-OCA2 region in 94 rs12913832:AA and AG Norwegians (43 blue-eyed and 51 brown-eyed) to search for novel blue eye colour variants. The new candidate variants were subsequently typed in a Norwegian biobank population (total n = 519) for population specific association analysis. We identified five new variants, rs74409036:A, rs78544415:T, rs72714116:T, rs191109490:C and rs551217952:C, to be the most promising candidates for explaining blue eye colour in individuals with the rs12913832:AA and AG genotype. Additionally, we confirmed the association of the missense variants rs74653330:T and rs121918166:T with blue eye colour, and observed lighter skin colour in rs74653330:T individuals. In total, 37 (86%) of the 43 blue-eyed rs12913832:AA and AG Norwegians could potentially be explained by these seven variants, and we suggest including them in future prediction models.

Genetic Variants in Telomerase Reverse Transcriptase Contribute to Solar Lentigines

Solar lentigines (SLs) are a hallmark of human skin aging. They result from chronic exposure to sunlight and other environmental stressors. Recent studies also imply genetic factors, but findings are partially conflicting and lack of replication. Through a multi-trait based analysis strategy, we discovered that genetic variants in telomerase reverse transcriptase were significantly associated with non-facial SL in two East Asian (Taizhou longitudinal cohort, n = 2,964 and National Survey of Physical Traits, n = 2,954) and one Caucasian population (SALIA, n = 462), top SNP rs2853672 (P-value for Taizhou longitudinal cohort = 1.32 × 10^‒28 and P-value for National Survey of Physical Traits = 3.66 × 10^‒17 and P-value for SALIA = 0.0007 and P_meta = 4.93 × 10^‒44). The same variants were nominally associated with facial SL but not with other skin aging or skin pigmentation traits. The SL-enhanced allele/haplotype upregulated the transcription of the telomerase reverse transcriptase gene. Of note, well-known telomerase reverse transcriptase‒related aging markers such as leukocyte telomere length and intrinsic epigenetic age acceleration were not associated with SL. Our results indicate a previously unrecognized role of telomerase reverse transcriptase in skin aging‒related lentigines formation.

Genomic diversity and post-admixture adaptation in the Uyghurs

Population admixture results in genome-wide combinations of genetic variants derived from different ancestral populations of distinct ancestry, thus providing a unique opportunity for understanding the genetic determinants of phenotypic variation in humans. Here, we used whole-genome sequencing of 92 individuals with high coverage (30–60×) to systematically investigate genomic diversity in the Uyghurs living in Xinjiang, China (XJU), an admixed population of both European-like and East-Asian-like ancestry. The XJU population shows greater genetic diversity, especially a higher proportion of rare variants, compared with their ancestral source populations, corresponding to greater phenotypic diversity of XJU. Admixture-induced functional variants in EDAR were associated with the diversity of facial morphology in XJU. Interestingly, the interaction of functional variants between SLC24A5 and OCA2 likely influences the diversity of skin pigmentation. Notably, selection has seemingly been relaxed or canceled in several genes with significantly biased ancestry, such as HERC2–OCA2. Moreover, signatures of post-admixture adaptation in XJU were identified, including genes related to metabolism (e.g. CYP2D6), digestion (e.g. COL11A1), olfactory perception (e.g. ANO2) and immunity (e.g. HLA). Our results demonstrated population admixture as a driving force, locally or globally, in shaping human genetic and phenotypic diversity as well as in adaptive evolution.

Further insight into the global variability of the OCA2-HERC2 locus for human pigmentation from multiallelic markers

The OCA2-HERC2 locus is responsible for the greatest proportion of eye color variation in humans. Numerous studies extensively described both functional SNPs and associated patterns of variation over this region. The goal of our study is to examine how these haplotype structures and allelic associations vary when highly variable markers such as microsatellites are used. Eleven microsatellites spanning 357 Kb of OCA2-HERC2 genes are analyzed in 3029 individuals from worldwide populations. We found that several markers display large differences in allele frequency (10% to 35% difference) among Europeans, East Asians and Africans. In Europe, the alleles showing increased frequency can also discriminate individuals with (IrisPlex) predicted blue and brown eyes. Distinct haplotypes are identified around the variants C and T of the functional SNP rs12913832 (associated to blue eyes), with linkage disequilibrium r² values significant up to 237 Kb. The haplotype carrying the allele rs12913832 C has high frequency (76%) in blue eye predicted individuals (30% in brown eye predicted individuals), while the haplotype associated to the allele rs12913832 T is restricted to brown eye predicted individuals. Finally, homozygosity values reach levels of 91% near rs12913832. Odds ratios show values of 4.2, 7.4 and 10.4 for four markers around rs12913832 and 7.1 for their core haplotype. Hence, this study provides an example on the informativeness of multiallelic markers that, despite their current limited potential contribution to forensic eye color prediction, supports the use of microsatellites for identifying causing variants showing similar genetic features and history.

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.

Hypomorphic alleles pose challenges in rare disease genomic variant interpretation

Exon skipping associated with an ATP7B intronic variant in a patient with Wilson's disease. (A) Sashimi plot visualization of aligned RNA sequencing data from proband liver tissue at ATP7B exons 14-13-12. The red track shows traditional RNA-seq data; the blue track shows RNA-seq enriched with exon capture (cDNA-cap) which achieves higher depth of protein-coding transcripts. The histogram indicates overall sequencing depth while arcs tabulate the number of junction-spanning reads supporting exon pairs. (B) The domain structure (top) and exon structure (bottom) of ATP7B. Loss of exon 13 (dashed box) would remove a transmembrane domain and disrupt the first phosphorylation domain.

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.

The distinctive geographic patterns of common pigmentation variants at the OCA2 gene

Oculocutaneous Albinism type 2 (OCA2) is a gene of great interest because of genetic variation affecting normal pigmentation variation in humans. The diverse geographic patterns for variant frequencies at OCA2 have been evident but have not been systematically investigated, especially outside of Europe. Here we examine population genetic variation in and near the OCA2 gene from a worldwide perspective. The very different patterns of genetic variation found across world regions suggest strong selection effects may have been at work over time. For example, analyses involving the variants that affect pigmentation of the iris argue that the derived allele of the rs1800407 single nucleotide polymorphism, which produces a hypomorphic protein, may have contributed to the previously demonstrated positive selection in Europe for the enhancer variant responsible for light eye color. More study is needed on the relationships of the genetic variation at OCA2 to variation in pigmentation in areas beyond Europe.

Studying human and nonhuman primate evolutionary biology with powerful in vitro and in vivo functional genomics tools

In recent years, tools for functional genomic studies have become increasingly feasible for use by evolutionary anthropologists. In this review, we provide brief overviews of several exciting in vitro techniques that can be paired with "-omics" approaches (e.g., genomics, epigenomics, transcriptomics, proteomics, and metabolomics) for potentially powerful evolutionary insights. These in vitro techniques include ancestral protein resurrection, cell line experiments using primary, immortalized, and induced pluripotent stem cells, and CRISPR-Cas9 genetic manipulation. We also discuss how several of these methods can be used in vivo, for transgenic organism studies of human and nonhuman primate evolution. Throughout this review, we highlight example studies in which these approaches have already been used to inform our understanding of the evolutionary biology of modern and archaic humans and other primates while simultaneously identifying future opportunities for anthropologists to use this toolkit to help answer additional outstanding questions in evolutionary anthropology.

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.

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.

The Genetics of Human Skin and Hair Pigmentation

Human skin and hair color are visible traits that can vary dramatically within and across ethnic populations. The genetic makeup of these traits-including polymorphisms in the enzymes and signaling proteins involved in melanogenesis, and the vital role of ion transport mechanisms operating during the maturation and distribution of the melanosome-has provided new insights into the regulation of pigmentation. A large number of novel loci involved in the process have been recently discovered through four large-scale genome-wide association studies in Europeans, two large genetic studies of skin color in Africans, one study in Latin Americans, and functional testing in animal models. The responsible polymorphisms within these pigmentation genes appear at different population frequencies, can be used as ancestry-informative markers, and provide insight into the evolutionary selective forces that have acted to create this human diversity.

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.

Off-target phenotypes in forensic DNA phenotyping and biogeographic ancestry inference: A resource

With recent advances in DNA sequencing technologies it has become feasible and cost effective to genotype larger marker sets for forensic purposes. Two technologies that make use of the larger marker sets have come into focus in forensic research and applications; inference of biogeographic ancestry (BGA) and forensic DNA phenotyping (FDP). These methods hold the promise to reveal information about a yet unknown perpetrator from a DNA sample. In contrast, DNA-profiling, that is a standard practice in case work, relies on matching DNA-profiles between crime scene material and suspects on a database of DNA-profiles. Markers for DNA-profiling were developed under the premise to reveal as little additional information about the human source of the profile as possible, the rationale being that personal privacy rights have to be balanced against the public interest in solving a crime. The same argument holds for markers used in BGA and FDP; these markers might also reveal information on off-target phenotypes (OTPs), that go beyond BGA and the phenotypes targeted in FDP. In particular, health related OTPs might shift the balance between privacy protection and public interest. However, to our knowledge, there is currently no convenient resource available to incorporate knowledge on OTPs in BGA and FDP assay design and application. In order to provide such a resource, we performed a systematic search for OTPs associated with a comprehensive set of markers (1766 SNPs) used or suggested to be used for BGA inference and FDP. In this set, we identified a relatively small number of 27 SNPs (1.53%) that convey information on diverse health related OTPs such as cancer risk, induced asthma, or risk of alcoholism. Some of these SNPs are commonly used for FDP and BGA across different marker sets. We conclude that the effects of SNP markers used in FDP and BGA on OTPs are currently limited, with few exceptions that should be considered in a balanced decision on assay design and application.

Clinical and Biological Characterization of Skin Pigmentation Diversity and Its Consequences on UV Impact

Skin color diversity is the most variable and noticeable phenotypic trait in humans resulting from constitutive pigmentation variability. This paper will review the characterization of skin pigmentation diversity with a focus on the most recent data on the genetic basis of skin pigmentation, and the various methodologies for skin color assessment. Then, melanocyte activity and amount, type and distribution of melanins, which are the main drivers for skin pigmentation, are described. Paracrine regulators of melanocyte microenvironment are also discussed. Skin response to sun exposure is also highly dependent on color diversity. Thus, sensitivity to solar wavelengths is examined in terms of acute effects such as sunburn/erythema or induced-pigmentation but also long-term consequences such as skin cancers, photoageing and pigmentary disorders. More pronounced sun-sensitivity in lighter or darker skin types depending on the detrimental effects and involved wavelengths is reviewed.

Estimating Time to the Common Ancestor for a Beneficial Allele

The haplotypes of a beneficial allele carry information about its history that can shed light on its age and the putative cause for its increase in frequency. Specifically, the signature of an allele's age is contained in the pattern of variation that mutation and recombination impose on its haplotypic background. We provide a method to exploit this pattern and infer the time to the common ancestor of a positively selected allele following a rapid increase in frequency. We do so using a hidden Markov model which leverages the length distribution of the shared ancestral haplotype, the accumulation of derived mutations on the ancestral background, and the surrounding background haplotype diversity. Using simulations, we demonstrate how the inclusion of information from both mutation and recombination events increases accuracy relative to approaches that only consider a single type of event. We also show the behavior of the estimator in cases where data do not conform to model assumptions, and provide some diagnostics for assessing and improving inference. Using the method, we analyze population-specific patterns in the 1000 Genomes Project data to estimate the timing of adaptation for several variants which show evidence of recent selection and functional relevance to diet, skin pigmentation, and morphology in humans.

FineMAV: prioritizing candidate genetic variants driving local adaptations in human populations

We present a new method, Fine-Mapping of Adaptive Variation (FineMAV), which combines population differentiation, derived allele frequency, and molecular functionality to prioritize positively selected candidate variants for functional follow-up. We calibrate and test FineMAV using eight experimentally validated "gold standard" positively selected variants and simulations. FineMAV has good sensitivity and a low false discovery rate. Applying FineMAV to the 1000 Genomes Project Phase 3 SNP dataset, we report many novel selected variants, including ones in TGM3 and PRSS53 associated with hair phenotypes that we validate using available independent data. FineMAV is widely applicable to sequence data from both human and other species.

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.

A Genetic Mechanism for Convergent Skin Lightening during Recent Human Evolution

Skin lightening among Eurasians is thought to have been a convergence occurring independently in Europe and East Asia as an adaptation to high latitude environments. Among Europeans, several genes responsible for such lightening have been found, but the information available for East Asians is much more limited. Here, a genome-wide comparison between dark-skinned Africans and Austro-Asiatic speaking aborigines and light-skinned northern Han Chinese identified the pigmentation gene OCA2, showing unusually deep allelic divergence between these groups. An amino acid substitution (His615Arg) of OCA2 prevalent in most East Asian populations—but absent in Africans and Europeans—was significantly associated with skin lightening among northern Han Chinese. Further transgenic and targeted gene modification analyses of zebrafish and mouse both exhibited the phenotypic effect of the OCA2 variant manifesting decreased melanin production. These results indicate that OCA2 plays an important role in the convergent skin lightening of East Asians during recent human evolution.

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.

Molecular genetic response of Xiphophorus maculatus-X. couchianus interspecies hybrid skin to UVB exposure

The phenotypic and genetic similarities between Xiphophorus and human melanoma render Xiphophorus a useful animal model for studying the genetic basis of melanoma etiology. In the Xiphophorus model, melanoma has been shown to be inducible by ultraviolet light (UVB) exposure among interspecies hybrids, but not in parental line fish similarly treated. This leads to questions of what genes are responsive to UVB exposure in the skin of the interspecies hybrids, as well as how parental alleles in hybrids may be differentially regulated and the potential roles they may play in induced melanomagenesis. To address these questions, we produced X. maculatus Jp 163 B×X. couchianus (Sp-Couch) F1 hybrid fish, exposed both hybrid and parental fish to UVB, and performed gene expression profiling of the skin using RNA-Seq methodology. We characterized a group of unique UVB-responsive genes in Sp-Couch hybrid including dct, pmela, tyr, tyrp1a, slc2a11b, rab38a, rab27, tspan10, slc45a2, oca2, slc24a5, ptn and mitfa. These genes are associated with melanin production and melanocyte proliferation. They were also up-regulated in Sp-Couch hybrid, indicating that their UVB response is hybridization initiated. In the hybrid, several melanin production and pigmentation related genes, including slc45a2, tspan10, dct, slc2a11b and ptn showed either X. couchianus or X. maculatus allele specific expression. The finding that these genes exhibit allele specific expression regulatory mechanisms in Sp-Couch hybrids, but do not exhibit a corresponding UVB response in either one of the parental fishes, may suggest UVB targets and imply mechanisms regarding the susceptibility of Sp-Couch to induced melanomagenesis.