A genomewide association study of skin pigmentation in a South Asian population

Assessment of IrisPlex-based multiplex for eye and skin color prediction with application to a Portuguese population

DNA phenotyping research is one of the most emergent areas of forensic genetics. Predictions of externally visible characteristics are possible through analysis of single nucleotide polymorphisms. These tools can provide police with "intelligence" in cases where there are no obvious suspects and unknown biological samples found at the crime scene do not result in any criminal DNA database hits. IrisPlex, an eye color prediction assay, revealed high prediction rates for blue and brown eye color in European populations. However, this is less predictive in some non-European populations, probably due to admixing. When compared to other European countries, Portugal has a relatively admixed population, resulting from a genetic influx derived from its proximity to and historical relations with numerous African territories. The aim of this work was to evaluate the utility of IrisPlex in the Portuguese population. Furthermore, the possibility of supplementing this multiplex with additional markers to also achieve skin color prediction within this population was evaluated. For that, IrisPlex was augmented with additional SNP loci. Eye and skin color prediction was estimated using the multinomial logistic regression and binomial logistic regression models, respectively. The results demonstrated eye color prediction accuracies of the IrisPlex system of 90 and 60% for brown and blue eye color, respectively, and 77% for intermediate eye color, after allele frequency adjustment. With regard to skin color, it was possible to achieve a prediction accuracy of 93%. In the future, phenotypic determination multiplexes must include additional loci to permit skin color prediction as presented in this study as this can be an advantageous tool for forensic investigation.

Exploration of SNP variants affecting hair colour prediction in Europeans

DNA profiling is a key tool for forensic analysis; however, current methods identify a suspect either by direct comparison or from DNA database searches. In cases with unidentified suspects, prediction of visible physical traits e.g. pigmentation or hair distribution of the DNA donors can provide important probative information. This study aimed to explore single nucleotide polymorphism (SNP) variants for their effect on hair colour prediction. A discovery panel of 63 SNPs consisting of already established hair colour markers from the HIrisPlex hair colour phenotyping assay as well as additional markers for which associations to human pigmentation traits were previously identified was used to develop multiplex assays based on SNaPshot single-base extension technology. A genotyping study was performed on a range of European populations (n = 605). Hair colour phenotyping was accomplished by matching donor's hair to a graded colour category system of reference shades and photography. Since multiple SNPs in combination contribute in varying degrees to hair colour predictability in Europeans, we aimed to compile a compact marker set that could provide a reliable hair colour inference from the fewest SNPs. The predictive approach developed uses a naïve Bayes classifier to provide hair colour assignment probabilities for the SNP profiles of the key SNPs and was embedded into the Snipper online SNP classifier ( http://mathgene.usc.es/snipper/ ). Results indicate that red, blond, brown and black hair colours are predictable with informative probabilities in a high proportion of cases. Our study resulted in the identification of 12 most strongly associated SNPs to hair pigmentation variation in six genes.

SLC transporters as therapeutic targets: emerging opportunities

Solute carrier (SLC) transporters - a family of more than 300 membrane-bound proteins that facilitate the transport of a wide array of substrates across biological membranes - have important roles in physiological processes ranging from the cellular uptake of nutrients to the absorption of drugs and other xenobiotics. Several classes of marketed drugs target well-known SLC transporters, such as neurotransmitter transporters, and human genetic studies have provided powerful insight into the roles of more-recently characterized SLC transporters in both rare and common diseases, indicating a wealth of new therapeutic opportunities. This Review summarizes knowledge on the roles of SLC transporters in human disease, describes strategies to target such transporters, and highlights current and investigational drugs that modulate SLC transporters, as well as promising drug targets.

Exploratory investigation of genetic associations with basal cell carcinoma risk: genome-wide association study in Jeju Island, Korea

AIM

Little is known about the genetic associations with Basal cell carcinoma (BCC) risk in non-Caucasian populations, in which BCC is rare, as in Korea. We here conducted a pilot genome-wide association study (GWAS) in 12 patients and 48 standard controls.

METHOD

A total of 263,511 SNPs were analyzed with the Illumina HumanOmni1 Quad v1.0 DNA Analysis BeadChip for cases and Korean HapMap 570K for controls.

RESULTS

SNP-based analyses, based on the allele genetic model with adjustment for sex and age showed suggestive associations with BCC risk for 6 SNPs with a P-value (P < 0.0005). However, these associations were not statistically significant after Bonferroni correction: rs1040503, rs2216491, rs13407683, rs4751072, rs9891263, and rs1368474. In addition, results from gene-based analyses showed suggestive associations with BCC risk for 33 candidate genes with a P-value (P <0.0005). Consistent with previous GWAS and replication studies in Caucasian populations, PADI6, RHOU and SLC45A2 were identified as having null associations with BCC (P > 0.05), likely due to the smaller sample size.

CONCLUSIONS

Although this was a small-scale negative study, to our knowledge, we have conducted the first GWAS for BCC risk in an Asian population. Further large studies in non-Caucasian populations are required to achieve statistical significance and confirm these findings.

Two novel tyrosinase (TYR) gene mutations with pathogenic impact on oculocutaneous albinism type 1 (OCA1)

Oculocutaneous albinism (OCA) is a heterogeneous group of autosomal recessive disorders resulting from mutations of the tyrosinase (TYR) gene and presents with either complete or partial absence of pigment in the skin, hair and eyes due to a defect in an enzyme involved in the production of melanin. In this study, mutations in the TYR gene of 30 unrelated Iranian OCA1 patients and 100 healthy individuals were examined using PCR-sequencing. Additionally, in order to predict the possible effects of new mutations on the structure and function of tyrosinase, these mutations were analyzed by SIFT, PolyPhen and I-Mutant 2 software. Here, two new pathogenic p.C89S and p.H180R mutations were detected in two OCA1 patients. Moreover, the R402Q and S192Y variants, which are common non-pathogenic polymorphisms, were detected in 17.5% and 35% of the patients, respectively. The outcome of this study has extended the genotypic spectrum of OCA1 patients, which paves the way for more efficient carrier detection and genetic counseling.

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

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

Characterizing the genetic differences between two distinct migrant groups from Indo-European and Dravidian speaking populations in India

Background

India is home to many ethnically and linguistically diverse populations. It is hypothesized that history of invasions by people from Persia and Central Asia, who are referred as Aryans in Hindu Holy Scriptures, had a defining role in shaping the Indian population canvas. A shift in spoken languages from Dravidian languages to Indo-European languages around 1500 B.C. is central to the Aryan Invasion Theory. Here we investigate the genetic differences between two sub-populations of India consisting of: (1) The Indo-European language speaking Gujarati Indians with genome-wide data from the International HapMap Project; and (2) the Dravidian language speaking Tamil Indians with genome-wide data from the Singapore Genome Variation Project.

Results

We implemented three population genetics measures to identify genomic regions that are significantly differentiated between the two Indian populations originating from the north and south of India. These measures singled out genomic regions with: (i) SNPs exhibiting significant variation in allele frequencies in the two Indian populations; and (ii) differential signals of positive natural selection as quantified by the integrated haplotype score (iHS) and cross-population extended haplotype homozygosity (XP-EHH). One of the regions that emerged spans the SLC24A5 gene that has been functionally shown to affect skin pigmentation, with a higher degree of genetic sharing between Gujarati Indians and Europeans.

Conclusions

Our finding points to a gene-flow from Europe to north India that provides an explanation for the lighter skin tones present in North Indians in comparison to South Indians.

Skin pigmentation, sun exposure and vitamin D levels in children of the Avon Longitudinal Study of Parents and Children

Background

It has been hypothesised that light skin pigmentation has arisen to ensure adequate levels of vitamin D as human populations moved out of Africa and into higher latitudes. Vitamin D, which is primarily obtained through exposure to sunlight (specifically ultraviolet radiation B (UVR-B)), has been inversely associated with several complex diseases. Greater sun exposure, on the other hand, is a well-known cause of skin cancer. The potential of UVR to be beneficial for some health outcomes but detrimental for others has prompted a public health debate on how to balance the positive and negative consequences of sun exposure. In this study we aimed to determine the validity of the evolutionary hypothesis linking lighter skin with higher vitamin D concentrations in a European population. Additionally, we aimed to examine the influence of pigmentation on personal behaviour towards sunlight exposure and the effects of this behaviour on vitamin D.

Methods

We combined genetic variants strongly associated with skin colour, tanning or freckling to create genetic scores for each of these phenotypes. We examined the association of the scores with pigmentary traits, sun exposure and serum 25-hydroxyvitamin D (25(OH)D) levels among children of the Avon Longitudinal Study of Parents and Children (ALSPAC, N = 661 to 5649).

Results

We found that fairer-skinned children, i.e. those with higher pigmentation score values, had higher levels of 25(OH)D (0.6 nmol/l; 95% CI 0.2, 1.0; per unit increase in skin colour score; N = 5649). These children also used more protection against the damaging effects of UVR.

Conclusions

In this population taking protective measures against sunburn and skin cancer does not seem to remove the positive effect that having a less pigmented skin has on vitamin D production. Our findings require further replication as skin pigmentation showed only a small effect on circulating 25(OH)D.

Neanderthal origin of the haplotypes carrying the functional variant Val92Met in the MC1R in modern humans

Skin color is one of the most visible and important phenotypes of modern humans. Melanocyte-stimulating hormone and its receptor played an important role in regulating skin color. In this article, we present evidence of Neanderthal introgression encompassing the melanocyte-stimulating hormone receptor gene MC1R. The haplotypes from Neanderthal introgression diverged with the Altai Neanderthal 103.3 ka, which postdates the anatomically modern human-Neanderthal divergence. We further discovered that all of the putative Neanderthal introgressive haplotypes carry the Val92Met variant, a loss-of-function variant in MC1R that is associated with multiple dermatological traits including skin color and photoaging. Frequency of this Neanderthal introgression is low in Europeans (∼5%), moderate in continental East Asians (∼30%), and high in Taiwanese aborigines (60-70%). As the putative Neanderthal introgressive haplotypes carry a loss-of-function variant that could alter the function of MC1R and is associated with multiple traits related to skin color, we speculate that the Neanderthal introgression may have played an important role in the local adaptation of Eurasians to sunlight intensity.

Molecular analysis of common polymorphisms within the human Tyrosinase locus and genetic association with pigmentation traits

We have compared the melanogenic activities of cultured melanocytes carrying two common TYR alleles as homozygous 192S-402R wild-type, heterozygous and homozygous variant. This includes assays of TYR protein, DOPAoxidase activity, glycosylation and temperature sensitivity of protein and DOPAoxidase levels. Homozygous wild-type strains on average had higher levels of TYR protein and enzyme activity than other genotypes. Homozygous 402Q/Q melanocytes produced significantly less TYR protein, displayed altered trafficking and glycosylation, with reduced DOPAoxidase. However, near wild-type TYR activity levels could be recovered at lower growth temperature. In a sample population from Southeast Queensland, these two polymorphisms were present on four TYR haplotypes, designated as WT 192S-402R, 192Y-402R and 192S-402Q with a double-variant 192Y-402Q of low frequency at 1.9%. Based on cell culture findings and haplotype associations, we have used an additive model to assess the penetrance of the ten possible TYR genotypes derived from the combination of these haplotypes.

Implications of the admixture process in skin color molecular assessment

The understanding of the complex genotype-phenotype architecture of human pigmentation has clear implications for the evolutionary history of humans, as well as for medical and forensic practices. Although dozens of genes have previously been associated with human skin color, knowledge about this trait remains incomplete. In particular, studies focusing on populations outside the European-North American axis are rare, and, until now, admixed populations have seldom been considered. The present study was designed to help fill this gap. Our objective was to evaluate possible associations of 18 single nucleotide polymorphisms (SNPs), located within nine genes, and one pseudogene with the Melanin Index (MI) in two admixed Brazilian populations (Gaucho, N = 352; Baiano, N = 148) with different histories of geographic and ethnic colonization. Of the total sample, four markers were found to be significantly associated with skin color, but only two (SLC24A5 rs1426654, and SLC45A2 rs16891982) were consistently associated with MI in both samples (Gaucho and Baiano). Therefore, only these 2 SNPs should be preliminarily considered to have forensic significance because they consistently showed the association independently of the admixture level of the populations studied. We do not discard that the other two markers (HERC2 rs1129038 and TYR rs1126809) might be also relevant to admixed samples, but additional studies are necessary to confirm the real importance of these markers for skin pigmentation. Finally, our study shows associations of some SNPs with MI in a modern Brazilian admixed sample, with possible applications in forensic genetics. Some classical genetic markers in Euro-North American populations are not associated with MI in our sample. Our results point out the relevance of considering population differences in selecting an appropriate set of SNPs as phenotype predictors in forensic practice.

Current status of the use of single-nucleotide polymorphisms in forensic practices

Forensic geneticists often use short tandem repeats (STRs) to solve cases. However, STRs can be insufficient when DNA samples are degraded due to environmental exposure and mass disasters, alleged and real relatives are genetically related in paternity or kinship analyses, or a suspect is lacking. In such cases, single-nucleotide polymorphisms (SNPs) can provide valuable information and thus should be seriously considered as a tool to help resolve challenging cases. In this review, the current status of SNP analyses in forensic applications and the comparative advantages and disadvantages of SNPs with other biomarkers are discussed.

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

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

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

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

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

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

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

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

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

Global patterns of diversity and selection in human tyrosinase gene

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

The SLC24 gene family of Na⁺/Ca²⁺-K⁺ exchangers: from sight and smell to memory consolidation and skin pigmentation

Members of the SLC24 gene family encode K(+)-dependent Na(+)/Ca(2+) exchangers (NCKX) that utilize both the inward Na(+) and outward K(+) gradients to extrude Ca(2+) from cells. There are five human SLC24 genes that play a role in biological process as diverse as vision in retinal rod and cone photoreceptors, olfaction, skin pigmentation and at least three of the five genes are also widely expressed in the brain. Here I review the functional, physiological and structural features of NCKX proteins that have emerged in the past few years.

Improved eye- and skin-color prediction based on 8 SNPs

Aim

To improve the 7-plex system to predict eye and skin color by increasing precision and detailed phenotypic descriptions.

Methods

Analysis of an eighth single nucleotide polymorphism (SNP), rs12896399 (SLC24A4), showed a statistically significant association with human eye color (P = 0.007) but a rather poor strength of agreement (κ = 0.063). This SNP was added to the 7-plex system (rs12913832 at HERC2, rs1545397 at OCA2, rs16891982 at SLC45A2, rs1426654 at SLC24A5, rs885479 at MC1R, rs6119471 at ASIP, and rs12203592 at IRF4). Further, the instruction guidelines on the interpretation of genotypes were changed to create a new 8-plex system. This was based on the analysis of an 803-sample training set of various populations. The newly developed 8-plex system can predict the eye colors brown, green, and blue, and skin colors light, not dark, and not light. It is superior to the 7-plex system with its additional ability to predict blue eye and light skin color.

Results

The 8-plex system was tested on an additional 212 samples, the test set. Analysis showed that the number of positive descriptions for eye colors as being brown, green, or blue increased significantly (P = 6.98e-15, z-score: -7.786). The error rate for eye-color prediction was low, at approximately 5%, while the skin color prediction showed no error in the test set (1% in training set).

Conclusions

We can conclude that the new 8-plex system for the prediction of eye and skin color substantially enhances its former version.

Why it is hard to find genes associated with social science traits: theoretical and empirical considerations

OBJECTIVES

We explain why traits of interest to behavioral scientists may have a genetic architecture featuring hundreds or thousands of loci with tiny individual effects rather than a few with large effects and why such an architecture makes it difficult to find robust associations between traits and genes.

METHODS

We conducted a genome-wide association study at 2 sites, Harvard University and Union College, measuring more than 100 physical and behavioral traits with a sample size typical of candidate gene studies. We evaluated predictions that alleles with large effect sizes would be rare and most traits of interest to social science are likely characterized by a lack of strong directional selection. We also carried out a theoretical analysis of the genetic architecture of traits based on R.A. Fisher's geometric model of natural selection and empirical analyses of the effects of selection bias and phenotype measurement stability on the results of genetic association studies.

RESULTS

Although we replicated several known genetic associations with physical traits, we found only 2 associations with behavioral traits that met the nominal genome-wide significance threshold, indicating that physical and behavioral traits are mainly affected by numerous genes with small effects.

CONCLUSIONS

The challenge for social science genomics is the likelihood that genes are connected to behavioral variation by lengthy, nonlinear, interactive causal chains, and unraveling these chains requires allying with personal genomics to take advantage of the potential for large sample sizes as well as continuing with traditional epidemiological studies.

Allelic frequencies of 20 visible phenotype variants in the korean population

The prediction of externally visible characteristics from DNA has been studied for forensic genetics over the last few years. Externally visible characteristics include hair, skin, and eye color, height, and facial morphology, which have high heritability. Recent studies using genome-wide association analysis have identified genes and variations that correlate with human visible phenotypes and developed phenotype prediction programs. However, most prediction models were constructed and validated based on genotype and phenotype information on Europeans. Therefore, we need to validate prediction models in diverse ethnic populations. In this study, we selected potentially useful variations for forensic science that are associated with hair and eye color, iris pattern, and facial morphology, based on previous studies, and analyzed their frequencies in 1,920 Koreans. Among 20 single nucleotide polymorphisms (SNPs), 10 SNPs were polymorphic, 6 SNPs were very rare (minor allele frequency < 0.005), and 4 SNPs were monomorphic in the Korean population. Even though the usability of these SNPs should be verified by an association study in Koreans, this study provides 10 potential SNP markers for forensic science for externally visible characteristics in the Korean population.

Determination of population origin: a comparison of autosomal SNPs, Y-chromosomal and mtDNA haplogroups using a Malagasy population as example

Y-chromosomal and mitochondrial DNA (mtDNA) polymorphisms have been used for population studies for a long time. However, there is another possibility to define the origin of a population: autosomal single-nucleotide polymorphisms (SNPs) whose allele frequencies differ considerably in different populations. In an attempt to compare the usefulness of these approaches we studied a population from Madagascar using all the three mentioned approaches. Former investigations of Malagasy maternal (mtDNA) and paternal (Y chromosome) lineages have led to the assumption that the Malagasy are an admixed population with an African and Asian-Indonesian heritage. Our additional study demonstrated that more than two-third of the Malagasy investigated showed clearly a West African genotype regarding only the autosomal SNPs despite the fact that 64% had an Asian mtDNA and more than 70% demonstrated an Asian-Indonesian heritage in either mtDNA or Y-chromosomal haplogroup or both. Nonetheless, the admixture of the Malagasy could be confirmed. A clear African or Asian-Indonesian heritage according to all the three DNA approaches investigated was only found in 14% and 1% of male samples, respectively. Not even the European or Northern African influences, detected in 9% of males (Y-chromosomal analysis) and 11% of samples (autosomal SNPs) were consistent. No Malagasy in our samples showed a European or Northern African origin in both categories. So, the analysis of autosomal SNPs could confirm the admixed character of the Malagasy population, even if it pointed to a greater African influence as detectable by Y-chromosomal or mtDNA analysis.

Colorful DNA polymorphisms in humans

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

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

Variation in human skin and eye color is substantial and especially apparent in admixed populations, yet the underlying genetic architecture is poorly understood because most genome-wide studies are based on individuals of European ancestry. We study pigmentary variation in 699 individuals from Cape Verde, where extensive West African/European admixture has given rise to a broad range in trait values and genomic ancestry proportions. We develop and apply a new approach for measuring eye color, and identify two major loci (HERC2[OCA2] P = 2.3×10⁻⁶², SLC24A5 P = 9.6×10⁻⁹) that account for both blue versus brown eye color and varying intensities of brown eye color. We identify four major loci (SLC24A5 P = 5.4×10⁻²⁷, TYR P = 1.1×10⁻⁹, APBA2[OCA2] P = 1.5×10⁻⁸, SLC45A2 P = 6×10⁻⁹) for skin color that together account for 35% of the total variance, but the genetic component with the largest effect (∼44%) is average genomic ancestry. Our results suggest that adjacent cis-acting regulatory loci for OCA2 explain the relationship between skin and eye color, and point to an underlying genetic architecture in which several genes of moderate effect act together with many genes of small effect to explain ∼70% of the estimated heritability.

Differences in skin and eye color are some of the most obvious traits that underlie human diversity, yet most of our knowledge regarding the genetic basis for these traits is based on the limited range of variation represented by individuals of European ancestry. We have studied a unique population in Cape Verde, an archipelago located off the West African coast, in which extensive mixing between individuals of Portuguese and West African ancestry has given rise to a broad range of phenotypes and ancestral genome proportions. Our results help to explain how genes work together to control the full range of pigmentary phenotypic diversity, provide new insight into the evolution of these traits, and provide a model for understanding other types of quantitative variation in admixed populations.

A genome-wide association study for primary open angle glaucoma and macular degeneration reveals novel Loci

Glaucoma and age-related macular degeneration (AMD) are the two leading causes of visual loss in the United States. We utilized a novel study design to perform a genome-wide association for both primary open angle glaucoma (POAG) and AMD. This study design utilized a two-stage process for hypothesis generation and validation, in which each disease cohort was utilized as a control for the other. A total of 400 POAG patients and 400 AMD patients were ascertained and genotyped at 500,000 loci. This study identified a novel association of complement component 7 (C7) to POAG. Additionally, an association of central corneal thickness, a known risk factor for POAG, was found to be associated with ribophorin II (RPN2). Linked monogenic loci for POAG and AMD were also evaluated for evidence of association, none of which were found to be significantly associated. However, several yielded putative associations requiring validation. Our data suggest that POAG is more genetically complex than AMD, with no common risk alleles of large effect.

Exome sequencing identifies SLC24A5 as a candidate gene for nonsyndromic oculocutaneous albinism

Oculocutaneous albinism (OCA) is a heterogeneous and autosomal recessive disorder with hypopigmentation in the eye, hair, and skin color. Four genes, TYR, OCA2, TYRP1, and SLC45A2, have been identified as causative genes for nonsyndromic OCA1-4, respectively. The genetic identity of OCA5 locus on 4q24 is unknown. Additional unknown OCA genes may exist as at least 5% of OCA patients have not been characterized during mutational screening in several populations. We used exome sequencing with a family-based recessive mutation model to determine that SLC24A5 is a previously unreported candidate gene for nonsyndromic OCA, which we designate as OCA6. Two deleterious mutations in this patient, c.591G>A and c.1361insT, were identified. We found apparent increase of immature melanosomes and less mature melanosomes in the patient's skin melanocytes. However, no defects in the platelet dense granules were observed, excluding typical Hermansky-Pudlak syndrome (HPS), a well-known syndromic OCA. Moreover, the SLC24A5 protein was reduced in steady-state levels in mouse HPS mutants with deficiencies in BLOC-1 and BLOC-2. Our results suggest that SLC24A5 is a previously unreported nonsyndromic OCA candidate gene and that the SLC24A5 transporter is transported into mature melanosomes by HPS protein complexes.

The impact of natural selection on health and disease: uses of the population genetics approach in humans

Investigations of the legacy of natural selection in the human genome have proved particularly informative, pinpointing functionally important regions that have participated in our genetic adaptation to the environment. Furthermore, genetic dissection of the intensity and type of selection acting on human genes can be used to predict involvement in different forms and severities of human diseases. We review here the progress made in population genetics studies toward understanding the effects of selection, in its different forms and intensities, on human genome diversity. We discuss some outstanding, robust examples of genes and biological functions subject to strong dietary, climatic and pathogen selection pressures. We also explore the possible relationship between cancer and natural selection, a topic that has been largely neglected because cancer is generally seen as a late-onset disease. Finally, we discuss how the present-day incidence of some diseases of modern societies may represent a by-product of past adaptation to other selective forces and changes in lifestyle. This perspective thus illustrates the value of adopting a population genetics approach in delineating the biological mechanisms that have played a major evolutionary role in the way humans have genetically adapted to different environments and lifestyles over time.

Cost-effective genome-wide estimation of allele frequencies from pooled DNA in Atlantic salmon (Salmo salar L.)

Background

New sequencing technologies have tremendously increased the number of known molecular markers (single nucleotide polymorphisms; SNPs) in a variety of species. Concurrently, improvements to genotyping technology have now made it possible to efficiently genotype large numbers of genome-wide distributed SNPs enabling genome wide association studies (GWAS). However, genotyping significant numbers of individuals with large number of SNPs remains prohibitively expensive for many research groups. A possible solution to this problem is to determine allele frequencies from pooled DNA samples, such ‘allelotyping’ has been presented as a cost-effective alternative to individual genotyping and has become popular in human GWAS. In this article we have tested the effectiveness of DNA pooling to obtain accurate allele frequency estimates for Atlantic salmon (Salmo salar L.) populations using an Illumina SNP-chip.

Results

In total, 56 Atlantic salmon DNA pools from 14 populations were analyzed on an Atlantic salmon SNP-chip containing probes for 5568 SNP markers, 3928 of which were bi-allelic. We developed an efficient quality control filter which enables exclusion of loci showing high error rate and minor allele frequency (MAF) close to zero. After applying multiple quality control filters we obtained allele frequency estimates for 3631 bi-allelic loci. We observed high concordance (r > 0.99) between allele frequency estimates derived from individual genotyping and DNA pools. Our results also indicate that even relatively small DNA pools (35 individuals) can provide accurate allele frequency estimates for a given sample.

Conclusions

Despite of higher level of variation associated with array replicates compared to pool construction, we suggest that both sources of variation should be taken into account. This study demonstrates that DNA pooling allows fast and high-throughput determination of allele frequencies in Atlantic salmon enabling cost-efficient identification of informative markers for discrimination of populations at various geographical scales, as well as identification of loci controlling ecologically and economically important traits.

Insights into the genetic architecture of early stage age-related macular degeneration: a genome-wide association study meta-analysis

Genetic factors explain a majority of risk variance for age-related macular degeneration (AMD). While genome-wide association studies (GWAS) for late AMD implicate genes in complement, inflammatory and lipid pathways, the genetic architecture of early AMD has been relatively under studied. We conducted a GWAS meta-analysis of early AMD, including 4,089 individuals with prevalent signs of early AMD (soft drusen and/or retinal pigment epithelial changes) and 20,453 individuals without these signs. For various published late AMD risk loci, we also compared effect sizes between early and late AMD using an additional 484 individuals with prevalent late AMD. GWAS meta-analysis confirmed previously reported association of variants at the complement factor H (CFH) (peak P = 1.5×10⁻³¹) and age-related maculopathy susceptibility 2 (ARMS2) (P = 4.3×10⁻²⁴) loci, and suggested Apolipoprotein E (ApoE) polymorphisms (rs2075650; P = 1.1×10⁻⁶) associated with early AMD. Other possible loci that did not reach GWAS significance included variants in the zinc finger protein gene GLI3 (rs2049622; P = 8.9×10⁻⁶) and upstream of GLI2 (rs6721654; P = 6.5×10⁻⁶), encoding retinal Sonic hedgehog signalling regulators, and in the tyrosinase (TYR) gene (rs621313; P = 3.5×10⁻⁶), involved in melanin biosynthesis. For a range of published, late AMD risk loci, estimated effect sizes were significantly lower for early than late AMD. This study confirms the involvement of multiple established AMD risk variants in early AMD, but suggests weaker genetic effects on the risk of early AMD relative to late AMD. Several biological processes were suggested to be potentially specific for early AMD, including pathways regulating RPE cell melanin content and signalling pathways potentially involved in retinal regeneration, generating hypotheses for further investigation.

NCKX5, a natural regulator of human skin colour variation, regulates the expression of key pigment genes MC1R and alpha-MSH and alters cholesterol homeostasis in normal human melanocytes

Natural human skin colour is determined both by environmental exposure to ultraviolet light and through inherited genetic variation in a very limited number of genes. Variation of a non-synonymous single-nucleotide polymorphism (nsSNP; rs1426654) in the gene (SLC24A5) encoding the NCKX5 protein is associated with differences in constitutive skin colour in South Asians. The nsSNP encodes the substitution of alanine for threonine at residue 111 (A111T) near a transmembrane region required for exchanger activity, a region which is highly conserved across different species and between NCKX family members. We have shown that NCKX5 is located at the trans-Golgi network of melanocytes and functions as a potassium-dependent sodium-calcium exchanger. When heterologously expressed, the 111T variant of NCKX5 shows significantly lower exchanger activity than the A111 variant. We have postulated that lower exchanger activity causes the reduced melanogenesis and lighter skin in Thr111-positive individuals. We used gene expression microarrays with qPCR replication and validation to assess the impact of siRNA-mediated knockdown of SLC24A5 on the transcriptome of cultured normal human melanocytes (NHM). Very few genes associated with melanogenesis were altered at the transcript level except for MC1R, suggesting that SLC24A5 interacts with at least one well-characterized melanogenic signalling pathway. More surprisingly, the expression of a number of cholesterol homeostatic genes was altered after SLC24A5 knockdown, and the total cholesterol content of NHM was increased. Cholesterol has previously been identified as a potential melanogenic regulator, and our data imply that NCKX5 exchanger function influences natural variation in skin pigmentation via a novel, unknown mechanism affecting cellular sterol levels.

Exome sequencing of only seven Qataris identifies potentially deleterious variants in the Qatari population

The Qatari population, located at the Arabian migration crossroads of African and Eurasia, is comprised of Bedouin, Persian and African genetic subgroups. By deep exome sequencing of only 7 Qataris, including individuals in each subgroup, we identified 2,750 nonsynonymous SNPs predicted to be deleterious, many of which are linked to human health, or are in genes linked to human health. Many of these SNPs were at significantly elevated deleterious allele frequency in Qataris compared to other populations worldwide. Despite the small sample size, SNP allele frequency was highly correlated with a larger Qatari sample. Together, the data demonstrate that exome sequencing of only a small number of individuals can reveal genetic variations with potential health consequences in understudied populations.

Genome-wide association studies of quantitatively measured skin, hair, and eye pigmentation in four European populations

Pigmentation of the skin, hair, and eyes varies both within and between human populations. Identifying the genes and alleles underlying this variation has been the goal of many candidate gene and several genome-wide association studies (GWAS). Most GWAS for pigmentary traits to date have been based on subjective phenotypes using categorical scales. But skin, hair, and eye pigmentation vary continuously. Here, we seek to characterize quantitative variation in these traits objectively and accurately and to determine their genetic basis. Objective and quantitative measures of skin, hair, and eye color were made using reflectance or digital spectroscopy in Europeans from Ireland, Poland, Italy, and Portugal. A GWAS was conducted for the three quantitative pigmentation phenotypes in 176 women across 313,763 SNP loci, and replication of the most significant associations was attempted in a sample of 294 European men and women from the same countries. We find that the pigmentation phenotypes are highly stratified along axes of European genetic differentiation. The country of sampling explains approximately 35% of the variation in skin pigmentation, 31% of the variation in hair pigmentation, and 40% of the variation in eye pigmentation. All three quantitative phenotypes are correlated with each other. In our two-stage association study, we reproduce the association of rs1667394 at the OCA2/HERC2 locus with eye color but we do not identify new genetic determinants of skin and hair pigmentation supporting the lack of major genes affecting skin and hair color variation within Europe and suggesting that not only careful phenotyping but also larger cohorts are required to understand the genetic architecture of these complex quantitative traits. Interestingly, we also see that in each of these four populations, men are more lightly pigmented in the unexposed skin of the inner arm than women, a fact that is underappreciated and may vary across the world.

Comprehensive candidate gene study highlights UGT1A and BNC2 as new genes determining continuous skin color variation in Europeans

Natural variation in human skin pigmentation is primarily due to genetic causes rooted in recent evolutionary history. Genetic variants associated with human skin pigmentation confer risk of skin cancer and may provide useful information in forensic investigations. Almost all previous gene-mapping studies of human skin pigmentation were based on categorical skin color information known to oversimplify the continuous nature of human skin coloration. We digitally quantified skin color into hue and saturation dimensions for 5,860 Dutch Europeans based on high-resolution skin photographs. We then tested an extensive list of 14,185 single nucleotide polymorphisms in 281 candidate genes potentially involved in human skin pigmentation for association with quantitative skin color phenotypes. Confirmatory association was revealed for several known skin color genes including HERC2, MC1R, IRF4, TYR, OCA2, and ASIP. We identified two new skin color genes: genetic variants in UGT1A were significantly associated with hue and variants in BNC2 were significantly associated with saturation. Overall, digital quantification of human skin color allowed detecting new skin color genes. The variants identified in this study may also contribute to the risk of skin cancer. Our findings are also important for predicting skin color in forensic investigations.

Functional assessment of human coding mutations affecting skin pigmentation using zebrafish

A major challenge in personalized medicine is the lack of a standard way to define the functional significance of the numerous nonsynonymous, single nucleotide coding variants that are present in each human individual. To begin to address this problem, we have used pigmentation as a model polygenic trait, three common human polymorphisms thought to influence pigmentation, and the zebrafish as a model system. The approach is based on the rescue of embryonic zebrafish mutant phenotypes by “humanized” zebrafish orthologous mRNA. Two hypomorphic polymorphisms, L374F in SLC45A2, and A111T in SLC24A5, have been linked to lighter skin color in Europeans. The phenotypic effect of a second coding polymorphism in SLC45A2, E272K, is unclear. None of these polymorphisms had been tested in the context of a model organism. We have confirmed that zebrafish albino fish are mutant in slc45a2; wild-type slc45a2 mRNA rescued the albino mutant phenotype. Introduction of the L374F polymorphism into albino or the A111T polymorphism into slc24a5 (golden) abolished mRNA rescue of the respective mutant phenotypes, consistent with their known contributions to European skin color. In contrast, the E272K polymorphism had no effect on phenotypic rescue. The experimental conclusion that E272K is unlikely to affect pigmentation is consistent with a lack of correlation between this polymorphism and quantitatively measured skin color in 59 East Asian humans. A survey of mutations causing human oculocutaneous albinism yielded 257 missense mutations, 82% of which are theoretically testable in zebrafish. The developed approach may be extended to other model systems and may potentially contribute to our understanding the functional relationships between DNA sequence variation, human biology, and disease.

Human pigmentation genes under environmental selection

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

Molecular basis of albinism in India: evaluation of seven potential candidate genes and some new findings

Albinism represents a group of genetic disorders with a broad spectrum of hypopigmentary phenotypes dependent on the genetic background of the patients. Oculocutaneous albinism (OCA) patients have little or no pigment in their eyes, skin and hair, whereas ocular albinism (OA) primarily presents the ocular symptoms, and the skin and hair color may vary from near normal to very fair. Mutations in genes directly or indirectly regulating melanin production are responsible for different forms of albinism with overlapping clinical features. In this study, 27 albinistic individuals from 24 families were screened for causal variants by a PCR-sequencing based approach. TYR, OCA2, TYRP1, SLC45A2, SLC24A5, TYRP2 and SILV were selected as candidate genes. We identified 5 TYR and 3 OCA2 mutations, majority in homozygous state, in 8 unrelated patients including a case of autosomal recessive ocular albinism (AROA). A homozygous 4-nucleotide novel insertion in SLC24A5 was detected in a person showing with extreme cutaneous hypopigmentation. A potential causal variant was identified in the TYRP2 gene in a single patient. Haplotype analyses in the patients carrying homozygous mutations in the classical OCA genes suggested founder effect. This is the first report of an Indian AROA patient harboring a mutation in OCA2. Our results also reveal for the first time that mutations in SLC24A5 could contribute to extreme hypopigmentation in humans.

The timing of pigmentation lightening in Europeans

The inverse correlation between skin pigmentation and latitude observed in human populations is thought to have been shaped by selective pressures favoring lighter skin to facilitate vitamin D synthesis in regions far from the equator. Several candidate genes for skin pigmentation have been shown to exhibit patterns of polymorphism that overlap the geospatial variation in skin color. However, little work has focused on estimating the time frame over which skin pigmentation has changed and on the intensity of selection acting on different pigmentation genes. To provide a temporal framework for the evolution of lighter pigmentation, we used forward Monte Carlo simulations coupled with a rejection sampling algorithm to estimate the time of onset of selective sweeps and selection coefficients at four genes associated with this trait in Europeans: KITLG, TYRP1, SLC24A5, and SLC45A2. Using compound haplotype systems consisting of rapidly evolving microsatellites linked to one single-nucleotide polymorphism in each gene, we estimate that the onset of the sweep shared by Europeans and East Asians at KITLG occurred approximately 30,000 years ago, after the out-of-Africa migration, whereas the selective sweeps for the European-specific alleles at TYRP1, SLC24A5, and SLC45A2 started much later, within the last 11,000-19,000 years, well after the first migrations of modern humans into Europe. We suggest that these patterns were influenced by recent increases in size of human populations, which favored the accumulation of advantageous variants at different loci.

Prediction of people's origin from degraded DNA--presentation of SNP assays and calculation of probability

The characterization of externally visible traits by DNA analysis is already an important tool when investigating ancient skeletal remains and may gain similar importance in future forensic DNA analysis. This, however, depends on the different legal regulations in the different countries. Besides eye or hair color, the population origin can provide crucial information in criminal prosecution. In this study, we present the analysis of 16 single-nucleotide polymorphisms (SNPs) combined to two robust SNaPshot assays with a detection threshold of 25-pg DNA. This assay was applied to 891 people from seven different populations (West Africa, North Africa, Turkey, Near East, Balkan states, North Europe, and Japan) with a thorough statistical evaluation. The prediction model was validated by an additional 125 individuals predominantly with an ancestry from those same regions. The specificity of these SNPs for the prediction of all population origins is very high (>90 %), but the sensitivity varied greatly (more than 90 % for West Africa, but only 25 % for the Near East). We could identify West Africans with a certainty of 100 %, and people from North Africa, the Balkan states, or North Europe nearly with the same reliability while determination of Turks or people from the Near East was rather difficult. In conclusion, the two SNaPshot assays are a powerful and reliable tool for the identification of people with an ancestry in one of the above listed populations, even from degraded DNA.

Gene mapping study for constitutive skin color in an isolated Mongolian population

To elucidate the genes responsible for constitutive human skin color, we measured the extent of skin pigmentation in the buttock, representative of lifelong non-sun-exposed skin, and conducted a gene mapping study on skin color in an isolated Mongolian population composed of 344 individuals from 59 families who lived in Dashbalbar, Mongolia. The heritability of constitutive skin color was 0.82, indicating significant genetic association on this trait. Through the linkage analysis using 1,039 short tandem repeat (STR) microsatellite markers, we identified a novel genomic region regulating constitutive skin color on 11q24.2 with an logarithm of odds (LOD) score of 3.39. In addition, we also found other candidate regions on 17q23.2, 6q25.1, and 13q33.2 (LOD ≥ 2). Family-based association tests on these regions with suggestive linkage peaks revealed ten and two significant single nucleotide polymorphisms (SNPs) on the linkage regions of chromosome 11 and 17, respectively. We were able to discover four possible candidate genes that would be implicated to regulate human skin color: ETS1, UBASH3B, ASAM, and CLTC.

Ethiopian genetic diversity reveals linguistic stratification and complex influences on the Ethiopian gene pool

Humans and their ancestors have traversed the Ethiopian landscape for millions of years, and present-day Ethiopians show great cultural, linguistic, and historical diversity, which makes them essential for understanding African variability and human origins. We genotyped 235 individuals from ten Ethiopian and two neighboring (South Sudanese and Somali) populations on an Illumina Omni 1M chip. Genotypes were compared with published data from several African and non-African populations. Principal-component and STRUCTURE-like analyses confirmed substantial genetic diversity both within and between populations, and revealed a match between genetic data and linguistic affiliation. Using comparisons with African and non-African reference samples in 40-SNP genomic windows, we identified "African" and "non-African" haplotypic components for each Ethiopian individual. The non-African component, which includes the SLC24A5 allele associated with light skin pigmentation in Europeans, may represent gene flow into Africa, which we estimate to have occurred ~3 thousand years ago (kya). The non-African component was found to be more similar to populations inhabiting the Levant rather than the Arabian Peninsula, but the principal route for the expansion out of Africa ~60 kya remains unresolved. Linkage-disequilibrium decay with genomic distance was less rapid in both the whole genome and the African component than in southern African samples, suggesting a less ancient history for Ethiopian populations.

Molecular genetic studies and delineation of the oculocutaneous albinism phenotype in the Pakistani population

Background

Oculocutaneous albinism (OCA) is caused by a group of genetically heterogeneous inherited defects that result in the loss of pigmentation in the eyes, skin and hair. Mutations in the TYR, OCA2, TYRP1 and SLC45A2 genes have been shown to cause isolated OCA. No comprehensive analysis has been conducted to study the spectrum of OCA alleles prevailing in Pakistani albino populations.

Methods

We enrolled 40 large Pakistani families and screened them for OCA genes and a candidate gene, SLC24A5. Protein function effects were evaluated using in silico prediction algorithms and ex vivo studies in human melanocytes. The effects of splice-site mutations were determined using an exon-trapping assay.

Results

Screening of the TYR gene revealed four known (p.Arg299His, p.Pro406Leu, p.Gly419Arg, p.Arg278*) and three novel mutations (p.Pro21Leu, p.Cys35Arg, p.Tyr411His) in ten families. Ex vivo studies revealed the retention of an EGFP-tagged mutant (p.Pro21Leu, p.Cys35Arg or p.Tyr411His) tyrosinase in the endoplasmic reticulum (ER) at 37°C, but a significant fraction of p.Cys35Arg and p.Tyr411His left the ER in cells grown at a permissive temperature (31°C). Three novel (p.Asp486Tyr, p.Leu527Arg, c.1045-15 T > G) and two known mutations (p.Pro743Leu, p.Ala787Thr) of OCA2 were found in fourteen families. Exon-trapping assays with a construct containing a novel c.1045-15 T > G mutation revealed an error in splicing. No mutation in TYRP1, SLC45A2, and SLC24A5 was found in the remaining 16 families. Clinical evaluation of the families segregating either TYR or OCA2 mutations showed nystagmus, photophobia, and loss of pigmentation in the skin or hair follicles. Most of the affected individuals had grayish-blue colored eyes.

Conclusions

Our results show that ten and fourteen families harbored mutations in the TYR and OCA2 genes, respectively. Our findings, along with the results of previous studies, indicate that the p.Cys35Arg, p.Arg278* and p.Gly419Arg alleles of TYR and the p.Asp486Tyr and c.1045-15 T > G alleles of OCA2 are the most common causes of OCA in Pakistani families. To the best of our knowledge, this study represents the first documentation of OCA2 alleles in the Pakistani population. A significant proportion of our cohort did not have mutations in known OCA genes. Overall, our study contributes to the development of genetic testing protocols and genetic counseling for OCA in Pakistani families.

Further development of forensic eye color predictive tests

In forensic analysis predictive tests for external visible characteristics (or EVCs), including inference of iris color, represent a potentially useful tool to guide criminal investigations. Two recent studies, both focused on forensic testing, have analyzed single nucleotide polymorphism (SNP) genotypes underlying common eye color variation (Mengel-From et al., Forensic Sci. Int. Genet. 4:323 and Walsh et al., Forensic Sci. Int. Genet. 5:170). Each study arrived at different recommendations for eye color predictive tests aiming to type the most closely associated SNPs, although both confirmed rs12913832 in HERC2 as the key predictor, widely recognized as the most strongly associated marker with blue and brown iris colors. Differences between these two studies in identification of other eye color predictors may partly arise from varying approaches to assigning phenotypes, notably those not unequivocally blue or dark brown and therefore occupying an intermediate iris color continuum. We have developed two single base extension assays typing 37 SNPs in pigmentation-associated genes to study SNP-genotype based prediction of eye, skin, and hair color variation. These assays were used to test the performance of different sets of eye color predictors in 416 subjects from six populations of north and south Europe. The presence of a complex and continuous range of intermediate phenotypes distinct from blue and brown eye colors was confirmed by establishing eye color populations compared to genetic clusters defined using Structure software. Our study explored the effect of an expanded SNP combination beyond six markers has on the ability to predict eye color in a forensic test without extending the SNP assay excessively - thus maintaining a balance between the test's predictive value and an ability to reliably type challenging DNA with a multiplex of manageable size. Our evaluation used AUC analysis (area under the receiver operating characteristic curves) and naïve Bayesian likelihood-based classification approaches. To provide flexibility in SNP-based eye color predictive tests in forensic applications we modified an online Bayesian classifier, originally developed for genetic ancestry analysis, to provide a straightforward system to assign eye color likelihoods from a SNP profile combining additional informative markers from the predictors analyzed by our study plus those of Walsh and Mengel-From. Two advantages of the online classifier is the ability to submit incomplete SNP profiles, a common occurrence when typing challenging DNA, and the ability to handle physically linked SNPs showing independent effect, by allowing the user to input frequencies from SNP pairs or larger combinations. This system was used to include the submission of frequency data for the SNP pair rs12913832 and rs1129038: indicated by our study to be the two SNPs most closely associated to eye color.