Novel quantitative pigmentation phenotyping enhances genetic association, epistasis, and prediction of human eye colour

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

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

A new approach to broaden the range of eye colour identifiable by IrisPlex in DNA phenotyping

IrisPlex system represents the most popular model for eye colour prediction. Based on six polymorphisms this model provides very accurate predictions that strongly depend on the definition of eye colour phenotypes. The aim of the present study was to introduce a new approach to improve eye colour prediction using the well-validated IrisPlex system. A sample of 238 individuals from a Southern Italian population was collected and for each of them a high-resolution image of eye was obtained. By quantifying eye colour variation into CIELAB space several clustering algorithms were applied for eye colour classification. Predictions with the IrisPlex model were obtained using eye colour categories defined by both visual inspection and clustering algorithms. IrisPlex system predicted blue and brown eye colour with high accuracy while it was inefficient in the prediction of intermediate eye colour. Clustering-based eye colour resulted in a significantly increased accuracy of the model especially for brown eyes. Our results confirm the validity of the IrisPlex system for forensic purposes. Although the quantitative approach here proposed for eye colour definition slightly improves its prediction accuracy, further research is still required to improve the model particularly for the intermediate eye colour prediction.

A genome-wide epistatic network underlies the molecular architecture of continuous color variation of body extremities

Deciphering the molecular architecture of coat coloration for a better understanding of the biological mechanisms underlying pigmentation still remains a challenge. We took advantage of a rabbit French experimental population in which both a pattern and a gradient of coloration from white to brown segregated within the himalayan phenotype. The whole experimental design was genotyped using the high density Affymetrix® AxiomOrcun™ SNP Array and phenotyped into 6 different groups ordered from the lighter to the darker. Genome-wide association analyses pinpointed an oligogenic determinism, under recessive and additive inheritance, involving genes already known in melanogenesis (ASIP, KIT, MC1R, TYR), and likely processed pseudogenes linked to ribosomal function, RPS20 and RPS14. We also identified (i) gene-gene interactions through ASIP:MC1R affecting light cream/beige phenotypes while KIT:RPS responsible of dark chocolate/brown colors and (ii) a genome-wide epistatic network involving several others coloration genes such as POT1 or HPS5. Finally, we determined the recessive inheritance of the English spotting phenotype likely involving a copy number variation affecting at least the end of the coding sequence of the KIT gene. Our analyses of coloration as a continuous trait allowed us to go beyond much of the established knowledge through the detection of additional genes and gene-gene interactions that may contribute to the molecular architecture of the coloration phenotype.

Prediction of eye, hair and skin colour in Latin Americans

Here we evaluate the accuracy of prediction for eye, hair and skin pigmentation in a dataset of > 6500 individuals from Mexico, Colombia, Peru, Chile and Brazil (including genome-wide SNP data and quantitative/categorical pigmentation phenotypes - the CANDELA dataset CAN). We evaluated accuracy in relation to different analytical methods and various phenotypic predictors. As expected from statistical principles, we observe that quantitative traits are more sensitive to changes in the prediction models than categorical traits. We find that Random Forest or Linear Regression are generally the best performing methods. We also compare the prediction accuracy of SNP sets defined in the CAN dataset (including 56, 101 and 120 SNPs for eye, hair and skin colour prediction, respectively) to the well-established HIrisPlex-S SNP set (including 6, 22 and 36 SNPs for eye, hair and skin colour prediction respectively). When training prediction models on the CAN data, we observe remarkably similar performances for HIrisPlex-S and the larger CAN SNP sets for the prediction of hair (categorical) and eye (both categorical and quantitative), while the CAN sets outperform HIrisPlex-S for quantitative, but not for categorical skin pigmentation prediction. The performance of HIrisPlex-S, when models are trained in a world-wide sample (although consisting of 80% Europeans, https://hirisplex.erasmusmc.nl), is lower relative to training in the CAN data (particularly for hair and skin colour). Altogether, our observations are consistent with common variation of eye and hair colour having a relatively simple genetic architecture, which is well captured by HIrisPlex-S, even in admixed Latin Americans (with partial European ancestry). By contrast, since skin pigmentation is a more polygenic trait, accuracy is more sensitive to prediction SNP set size, although here this effect was only apparent for a quantitative measure of skin pigmentation. Our results support the use of HIrisPlex-S in the prediction of categorical pigmentation traits for forensic purposes in Latin America, while illustrating the impact of training datasets on its accuracy.

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.

The challenge of predicting human pigmentation traits in degraded bone samples with the MPS-based HIrisPlex-S system

Identification of human remains is an important part of human DNA analysis studies. STR and mitochondrial DNA markers are well suited for the analysis of degraded biological samples including bone material. However, these DNA markers may be useless when reference material is not available. In these cases, predictive DNA analysis can support the process of human identification by providing investigative leads. Forensic DNA phenotyping has progressed significantly by offering new methods based on massively parallel sequencing technology, but the frequent degradation processes observed in skeletal remains can make analysis of such samples challenging. In this study, we demonstrate the usefulness of a recently established Ion AmpliSeq^TM HIrisPlex-S panel using Ion Torrent technology for analyzing bone samples that show different levels of DNA degradation. In total, 63 bone samples at post-mortem intervals up to almost 80 years were genotyped and eye, hair and skin colour predictions were performed using the HIrisPlex-S models. Following the recommended coverage thresholds, it was possible to establish full DNA profiles comprising of 41 DNA variants for 35 samples (55.6%). For 5 samples (7.9%) no DNA profiles were generated. The remaining 23 samples (36.5%) produced partial profiles and showed a clear underperformance of 3 HIrisPlex-S SNPs - rs1545397 (OCA2), rs1470608 (OCA2) and rs10756819 (BNC2), all used for skin colour prediction only. None of the 23 samples gave complete genotypes needed for skin colour prediction was obtained, and in 7 of them (25.9%) the 3 underperformed SNPs were the cause. At the same time, the prediction of eye and hair colour using complete IrisPlex and HIrisPlex profiles could be made for these 23 samples in 20 (87.0%) and 12 cases (52.2%), respectively. Complete HIrisPlex-S profiles were generated from as little as 49 pg of template DNA. Five samples for which the HIrisPlex-S analysis failed, consistently failed in standard STR analysis. Importantly, the 3 underperforming SNPs produced significantly lower number of reads in good quality samples. Nonetheless, the AUC loss resulting from missing data for these 3 SNPs is not considered large (≤0.004) and the prediction of pigmentation from partial profiles is also available in the current HPS tool. The study shows that DNA degradation and the resulting loss of data are the most serious challenge to DNA phenotyping of skeletal remains. Although the newly developed HIrisPlex-S panel has been successfully validated in the current research, primer redesign for the 3 underperforming SNPs in the MPS design should be considered in the future.

Genome-wide association studies for iris pigmentation and heterochromia patterns in Large White pigs

Eye colour genetics have been extensively studied in humans since the rediscovery of Mendel's laws. This trait was first interpreted using simplistic genetic models but soon it was realised that it is more complex. In this study, we analysed eye colour variability in a Large White pig population (n = 897) and report the results of GWASs based on several comparisons including pigs having four main eye colour categories (three with both pigmented eyes of different brown grades: pale, 17.9%; medium, 14.8%; and dark, 54.3%; another one with both eyes completely depigmented, 3.8%) and heterochromia patterns (heterochromia iridis - depigmented iris sectors in pigmented irises, 3.2%; heterochromia iridum - one whole eye iris of depigmented phenotype and the other eye with the iris completely pigmented, 5.9%). Pigs were genotyped with the Illumina PorcineSNP60 BeadChip and GEMMA was used for the association analyses. The results indicated that SLC45A2 (on chromosome 16, SSC16), EDNRB (SSC11) and KITLG (SSC5) affect the different grades of brown pigmentation of the eyes, the bilateral eye depigmentation defect and the heterochromia iridis defect recorded in this white pig population respectively. These genes are involved in several mechanisms affecting pigmentation. Significant associations for the eye depigmented patterns were also identified for SNPs on two SSC4 regions (including two candidate genes: NOTCH2 and PREX2) and on SSC6, SSC8 and SSC14 (including COL17A1 as candidate gene). This study provided useful information to understand eye pigmentation mechanisms, further valuing the pig as animal model to study complex phenotypes in humans.

Clinical and genetic variability in children with partial albinism

Individuals who have ocular features of albinism and skin pigmentation in keeping with their familial background present a considerable diagnostic challenge. Timely diagnosis through genomic testing can help avert diagnostic odysseys and facilitates accurate genetic counselling and tailored specialist management. Here, we report the clinical and gene panel testing findings in 12 children with presumed ocular albinism. A definitive molecular diagnosis was made in 8/12 probands (67%) and a possible molecular diagnosis was identified in a further 3/12 probands (25%). TYR was the most commonly mutated gene in this cohort (75% of patients, 9/12). A disease-causing TYR haplotype comprised of two common, functional polymorphisms, TYR c.[575 C > A;1205 G > A] p.[(Ser192Tyr);(Arg402Gln)], was found to be particularly prevalent. One participant had GPR143-associated X-linked ocular albinism and another proband had biallelic variants in SLC38A8, a glutamine transporter gene associated with foveal hypoplasia and optic nerve misrouting without pigmentation defects. Intriguingly, 2/12 individuals had a single, rare, likely pathogenic variant in each of TYR and OCA2 - a significant enrichment compared to a control cohort of 4046 individuals from the 100,000 genomes project pilot dataset. Overall, our findings highlight that panel-based genetic testing is a clinically useful test with a high diagnostic yield in children with partial/ocular albinism.

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.

Genetic modifiers as relevant biological variables of eye disorders

From early in the study of mammalian genetics, it was clear that modifiers can have a striking influence on phenotypes. Today, several modifiers have now been studied in enough detail to allow a glimpse of how they function and influence our perspective of disease. With respect to diseases of the eye, some modifiers are an important source of phenotypic variation that can elucidate how genes function in networks to collectively shape ocular anatomy and physiology, thus influencing our understanding of basic biology. Other modifiers represent an opportunity for new therapeutic targets, whose manipulation could be used to mitigate ophthalmic disease. Here, we review progress in the study of genetic modifiers of eye disorders, with examples from mice and humans that together illustrate the ubiquitous nature of genetic modifiers and why they are relevant biological variables in experimental design. Special emphasis is given to ophthalmic modifiers in mice, especially those relevant to selection of genetic background and those that might inadvertently be a source of experimental variability. These modifiers are capable of influencing interpretations of many experiments using targeted genome manipulations such as knockouts or transgenics. Whereas there are fewer examples of modifiers of eye disorders in humans with a molecular identification, there is ample evidence that they exist and should be considered as a relevant biological variable in human genetic studies as well.

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.