HERC2 rs12913832 modulates human pigmentation by attenuating chromatin-loop formation between a long-range enhancer and the OCA2 promoter

The fork remodeler helicase-like transcription factor in cancer development: all at once

The Helicase-like Transcription Factor (HLTF) is a member of the SNF2-family of fork remodelers, primarily studied for its capacity to provide DNA Damage Tolerance (DDT) and to induce replication fork reversal (RFR). HLTF is recruited at stalled forks where both its ATPase motor and HIP116 Rad5p N-terminal (HIRAN) domains are necessary for regulating its interaction with DNA. HIRAN bestows specificity to ssDNA 3'-end and imparts branch migration as well as DNA remodeling capabilities facilitating damage repair. Both expression regulation and mutation rate affect HLTF activity. Gene hypermethylation induces loss of HLTF function, in particular in colorectal cancer (CRC), implying a tumour suppressor role. Surprisingly, a correlation between hypermethylation and HLTF mRNA upregulation has also been observed, even within the same cancer type. In many cancers, both complex mutation patterns and the presence of gene Copy Number Variations (CNVs) have been reported. These conditions affect the amount of functional HLTF and question the physiological role of this fork remodeler. This review offers a systematic collection of the presently strewed information regarding HLTF, its structural and functional characteristics, the multiple roles in DDT and the regulation in cancer progression highlighting new research perspectives.

Single-cell profiling of MC1R-inhibited melanocytes

The human red hair color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to define a MC1R-inhibited Gene Signature (MiGS) comprising a large set of previously unidentified genes which may be implicated in melanogenesis and oncogenic transformation. We show that one of the candidate MiGS genes, TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into further mechanisms by which melanocytes with reduced MC1R signaling may regulate pigmentation and offer new candidates of study toward understanding how individuals with the RHC phenotype are predisposed to melanoma.

Identification of Candidate Genes for Red-Eyed (Albinism) Domestic Guppies Using Genomic and Transcriptomic Analyses

Guppies are small tropical fish with brightly colored bodies and variable tail shapes. There are two phenotypes of domestic guppy eye color: red and black. The wild type is black-eyed. The main object of this study was to identify candidate genes for the red-eyed phenotype in domestic guppies. We hope to provide molecular genetic information for the development of new domestic guppy strains. Additionally, the results also contribute to basic research concerning guppies. In this study, 121 domestic guppies were used for genomic analysis (GWAS), and 44 genes were identified. Furthermore, 21 domestic guppies were used for transcriptomic analysis, and 874 differentially expressed genes (DEGs) were identified, including 357 upregulated and 517 downregulated genes. Through GO and KEGG enrichment, we identified some important terms or pathways mainly related to melanin biosynthesis and ion transport. qRT-PCR was also performed to verify the differential expression levels of four important candidate genes (TYR, OCA2, SLC45A2, and SLC24A5) between red-eyed and black-eyed guppies. Based on the results of genomic and transcriptomic analyses, we propose that OCA2 is the most important candidate gene for the red-eyed phenotype in guppies.

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

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

Whole Exome Sequencing as a First-Line Molecular Genetic Test in Developmental and Epileptic Encephalopathies

Developmental and epileptic encephalopathies (DEE) are severe neurodevelopmental disorders characterized by recurrent, usually early-onset, epileptic seizures accompanied by developmental impairment often related to both underlying genetic etiology and abnormal epileptiform activity. Today, next-generation sequencing technologies (NGS) allow us to sequence large portions of DNA quickly and with low costs. The aim of this study is to evaluate the use of whole-exome sequencing (WES) as a first-line molecular genetic test in a sample of subjects with DEEs characterized by early-onset drug-resistant epilepsies, associated with global developmental delay and/or intellectual disability (ID). We performed 82 WESs, identifying 35 pathogenic variants with a detection rate of 43%. The identified variants were highlighted on 29 different genes including, 3 new candidate genes (KCNC2, STXBP6, DHRS9) for DEEs never identified before. In total, 23 out of 35 (66%) de novo variants were identified. The most frequently identified type of inheritance was autosomal dominant de novo (60%) followed by autosomal recessive in homozygosity (17%) and heterozygosity (11%), autosomal dominant inherited from parental mosaicism (6%) and X-linked dominant de novo (6%). The most frequent mutations identified were missense (75%) followed by frameshift deletions (16%), frameshift duplications (5%), and splicing mutations (3%). Considering the results obtained in the present study we support the use of WES as a form of first-line molecular genetic testing in DEEs.

Genetic landscape of forensic DNA phenotyping markers among Mediterranean populations

Forensic DNA Phenotyping can reveal the appearance of an unknown individual by predicting the External Visible Characteristics (EVC) from DNA obtained at the crime scene. Our aim is to characterize the genetic landscape of Human identification markers responsible for EVC among Mediterranean populations compared to other worldwide groups. We conducted an exhaustive search for genes involved in EVC variation. Then, variants located on these genes were extracted from public genotypic data of Mediterranean, American, African and East Asiatic populations. The genetic landscape of these Human identification markers, their allelic distribution and admixture analyses, were determined using plink, R and ADMIXTURE softwares. Our results showed that the Mediterranean populations appear close to the Mexican populations and distinguished from sub Saharan African populations living in the USA and from East Asiatic populations. We highlighted a total of 103454 common variants shared between the studied populations and among them, 25 common variants associated with EVC. Interestingly, genotype frequencies results showed that the rs17646946, rs13016869, rs977588, rs1805008 and rs2240751 variants located respectively in the TCHH, PRKCE, OCA2, MC1R and MFSD12 genes are significantly different between the Mediterranean and Asiatic populations. The genotype frequencies of the variants rs977589 and rs7179994 located in the OCA2 gene, and of rs12913832 and rs2240751 located respectively in HERC2 and MFSD12 genes are significantly different between the Mediterranean and American populations. Our work generates a large number of EVC variants that could be a valuable resource for future studies in the forensic field.

Whole Genome Resequencing Reveals Selection Signals Related to Wool Color in Sheep

Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

Specific Learning Disorders: Variation Analysis of 15 Candidate Genes in 9 Multiplex Families

Background and Objectives: Specific Learning Disorder (SLD) is a complex neurobiological disorder characterized by a persistent difficult in reading (dyslexia), written expression (dysgraphia), and mathematics (dyscalculia). The hereditary and genetic component is one of the underlying causes of SLD, but the relationship between genes and the environment should be considered. Several genetic studies were performed in different populations to identify causative genes. Materials and Methods: Here, we show the analysis of 9 multiplex families with at least 2 individuals diagnosed with SLD per family, with a total of 37 persons, 21 of whom are young subjects with SLD, by means of Next-Generation Sequencing (NGS) to identify possible causative mutations in a panel of 15 candidate genes: CCPG1, CYP19A1, DCDC2, DGKI, DIP2A, DYM, GCFC2, KIAA0319, MC5R, MRPL19, NEDD4L, PCNT, PRMT2, ROBO1, and S100B. Results: We detected, in eight families out nine, SNP variants in the DGKI, DIP2A, KIAA0319, and PCNT genes, even if in silico analysis did not show any causative effect on this behavioral condition. In all cases, the mutation was transmitted by one of the two parents, thus excluding the case of de novo mutation. Moreover, the parent carrying the allelic variant transmitted to the children, in six out of seven families, reports language difficulties. Conclusions: Although the present results cannot be considered conclusive due to the limited sample size, the identification of genetic variants in the above genes can provide input for further research on the same, as well as on other genes/mutations, to better understand the genetic basis of this disorder, and from this perspective, to better understand also the neuropsychological and social aspects connected to this disorder, which affects an increasing number of young people.

Forensic DNA Phenotyping: Genes and Genetic Variants for Eye Color Prediction

In recent decades, the use of genetic polymorphisms related to specific phenotypes, such as eye color, has greatly contributed to the development of the research field called forensic DNA phenotyping (FDP), enabling the investigators of crime cases to reduce the number of suspects, making their work faster and more precise. Eye color is a polygenic phenotype, and many genetic variants have been highlighted, with the major contributor being the HERC2-OCA2 locus, where many single nucleotide variations (SNPs) were identified. Interestingly, the HERC2-OCA2 locus, containing the intronic SNP rs12913832, the major eye color determinant, shows a high level of evolutionary conservation across many species of vertebrates. Currently, there are some genetic panels to predict eye color by genomic DNA analysis, even if the exact role of the SNP variants in the formation of eye color is still poorly understood, with a low level of predictivity in the so-called intermediate eye color. Many variants in OCA2, HERC2, and other genes lie in introns or correspond to synonymous variants, highlighting greater complexity in the mechanism of action of such genes than a simple missense variation. Here, we show the main genes involved in oculocutaneous pigmentation and their structural and functional features, as well as which genetic variants show the highest level of eye color predictivity in currently used FDP assays. Despite the great recent advances and impact of FDP in criminal cases, it is necessary to enhance scientific research to better understand the mechanism of action behind each genetic variant involved in eye color, with the goal of obtaining higher levels of prediction.

Ethnicity is not biology: retinal pigment score to evaluate biological variability from ophthalmic imaging using machine learning

Background

Few metrics exist to describe phenotypic diversity within ophthalmic imaging datasets, with researchers often using ethnicity as an inappropriate marker for biological variability.

Methods

We derived a continuous, measured metric, the retinal pigment score (RPS), that quantifies the degree of pigmentation from a colour fundus photograph of the eye. RPS was validated using two large epidemiological studies with demographic and genetic data (UK Biobank and EPIC-Norfolk Study).

Findings

A genome-wide association study (GWAS) of RPS from UK Biobank identified 20 loci with known associations with skin, iris and hair pigmentation, of which 8 were replicated in the EPIC-Norfolk cohort. There was a strong association between RPS and ethnicity, however, there was substantial overlap between each ethnicity and the respective distributions of RPS scores.

Interpretation

RPS serves to decouple traditional demographic variables, such as ethnicity, from clinical imaging characteristics. RPS may serve as a useful metric to quantify the diversity of the training, validation, and testing datasets used in the development of AI algorithms to ensure adequate inclusion and explainability of the model performance, critical in evaluating all currently deployed AI models. The code to derive RPS is publicly available at: https://github.com/uw-biomedical-ml/retinal-pigmentation-score.

Funding

The authors did not receive support from any organisation for the submitted work.

Loss of MC1R signaling implicates TBX3 in pheomelanogenesis and melanoma predisposition

The human Red Hair Color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA-sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to reveal a Pheomelanin Gene Signature (PGS) comprising genes implicated in melanogenesis and oncogenic transformation. We show that TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, is part of the PGS and binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into mechanisms by which MC1R signaling regulates pigmentation and how individuals with the RHC phenotype are predisposed to melanoma.

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.

Massively parallel reporter assays and variant scoring identified functional variants and target genes for melanoma loci and highlighted cell-type specificity

The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRAs) by using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRAs, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (eQTLs or meQTLs, respectively) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis-regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription-factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRAs along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional.

Predicting Eye and Hair Color in a Turkish Population Using the HIrisPlex System

Forensic DNA Phenotyping (FDP) can reveal the appearance of an unknown individual by predicting the ancestry, phenotype (i.e., hair, eye, skin color), and age from DNA obtained at the crime scene. The HIrisPlex system has been developed to simultaneously predict eye and hair color. However, the prediction accuracy of the system needs to be assessed for the tested population before implementing FDP in casework. In this study, we evaluated the performance of the HIrisPlex system on 149 individuals from the Turkish population. We applied the single-based extension (SNaPshot chemistry) method and used the HIrisPlex online tool to test the prediction of the eye and hair colors. The accuracy of the HIrisPlex system was assessed through the calculation of the area under the receiver characteristic operating curves (AUC), sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV). The results showed that the proposed method successfully predicted the eye and hair color, especially for blue (100%) and brown (95.60%) eye and black (95.23) and brown (98.94) hair colors. As observed in previous studies, the system failed to predict intermediate eye color, representing 25% in our cohort. The majority of incorrect predictions were observed for blond hair color (40.7%). Previous HIrisPlex studies have also noted difficulties with these phenotypes. Our study shows that the HIrisPlex system can be applied to forensic casework in Turkey with careful interpretation of the data, particularly intermediate eye color and blond hair color.

Uncovering the extensive trade-off between adaptive evolution and disease susceptibility

Favored mutations in the human genome may make the carriers adapt to changing environments and lifestyles but also susceptible to specific diseases. The scale and details of the trade-off between adaptive evolution and disease susceptibility are unclear because most favored mutations in different populations remain unidentified. As no statistical test can discriminate favored mutations from nearby hitchhiking neutral ones, we report a deep-learning network (DeepFavored) to integrate multiple statistical tests and divide identifying favored mutations into two subtasks. We identify favored mutations in three human populations and analyzed the correlation between favored/hitchhiking mutations and genome-wide association study (GWAS) sites. Both favored and hitchhiking neutral mutations are enriched in GWAS sites with population-specific features, and the enrichment and population specificity are prominent in genes in specific Gene Ontology (GO) terms. These provide evidence for extensive and population-specific trade-offs between adaptive evolution and disease susceptibility. The unveiled scale helps understand and investigate differences and diseases of humans.

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.

QTL Mapping for Age-Related Eye Pigmentation in the Pink-Eyed Dilution Castaneus Mutant Mouse

Pink-eyed dilution castaneus (Oca2p-cas) is a mutant gene on mouse chromosome 7 that arose spontaneously in wild Mus musculus castaneus. Homozygotes for Oca2p-cas exhibit pink eyes and a light gray coat throughout life. In an ordinary mutant strain carrying Oca2p-cas, we previously discovered a novel spontaneous mutation that gradually increases melanin pigmentation in the eyes and coat with aging, and we developed a novel mutant strain that was fixed for the novel phenotype. The purpose of this study was to map major quantitative trait loci (QTLs) for the novel pigmentation phenotype and for expression levels of four important melanogenesis genes, microphthalmia-associated transcription factor (Mitf), tyrosinase (Tyr), tyrosinase-related protein-1 (Tyrp1) and dopachrome tautomerase (Dct). We developed 69 DNA markers and created 303 F2 mice from two reciprocal crosses between novel and ordinary mutant strains. The QTL analysis using a selective genotyping strategy revealed a significant QTL for eye pigmentation between 34 and 64 Mb on chromosome 13. This QTL explained approximately 20% of the phenotypic variance. The QTL allele derived from the novel strain increased pigmentation. Although eye pigmentation was positively correlated with Dct expression, no expression QTLs were found, suggesting that the pigmentation QTL on chromosome 13 may not be directly in the pathway of any of the four melanogenesis genes. This study is the first step toward identifying a causal gene for the novel spontaneous phenotype in mice and is expected to discover a new regulatory mechanism for complex melanin biosynthesis during aging.

Investigating the genetic architecture of eye colour in a Canadian cohort

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

•

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

•

Multiple independent candidate causal variants were identified across HERC2/OCA2

•

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

•

Colocalization of eye color signals with expression and methylation profiles

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

Zebrafish Syndromic Albinism Models as Tools for Understanding and Treating Pigment Cell Disease in Humans

Simple Summary

Zebrafish (Danio rerio) is an emerging model for studying many diseases, including disorders originating in black pigment cells, melanocytes. In this review of the melanocyte literature, we discuss the current knowledge of melanocyte biology relevant to understanding different forms of albinism and the potential of the zebrafish model system for finding novel mechanisms and treatments.

Abstract

Melanin is the pigment that protects DNA from ultraviolet (UV) damage by absorbing excess energy. Melanin is produced in a process called melanogenesis. When melanogenesis is altered, diseases such as albinism result. Albinism can result in an increased skin cancer risk. Conversely, black pigment cell (melanocyte) development pathways can be misregulated, causing excessive melanocyte growth that leads to melanoma (cancer of melanocytes). Zebrafish is an emerging model organism used to study pigment disorders due to their high fecundity, visible melanin development in melanophores (melanocytes in mammals) from 24 h post-fertilization, and conserved melanogenesis pathways. Here, we reviewed the conserved developmental pathways in zebrafish melanophores and mammalian melanocytes. Additionally, we summarized the progress made in understanding pigment cell disease and evidence supporting the strong potential for using zebrafish to find novel treatment options for albinism.

Candidate Genes in Bull Semen Production Traits: An Information Approach Review

Semen quality plays a crucial role in the successful implementation of breeding programs, especially where artificial insemination (AI) is practiced. Bulls with good semen traits have good fertility and can produce a volume of high semen per ejaculation. The aim of this review is to use an information approach to highlight candidate genes and their relation to bull semen production traits. The use of genome-wide association studies (GWAS) has been demonstrated to be successful in identifying genomic regions and individual variations associated with production traits. Studies have reported over 40 genes associated with semen traits using Illumina BeadChip single-nucleotide polymorphism (SNPs).

Epigenetics of alcohol-related liver diseases

Alcohol-related liver disease (ARLD) is a primary cause of chronic liver disease in the United States. Despite advances in the diagnosis and management of ARLD, it remains a major public health problem associated with significant morbidity and mortality, emphasising the need to adopt novel approaches to the study of ARLD and its complications. Epigenetic changes are increasingly being recognised as contributing to the pathogenesis of multiple disease states. Harnessing the power of innovative technologies for the study of epigenetics (e.g., next-generation sequencing, DNA methylation assays, histone modification profiling and computational techniques like machine learning) has resulted in a seismic shift in our understanding of the pathophysiology of ARLD. Knowledge of these techniques and advances is of paramount importance for the practicing hepatologist and researchers alike. Accordingly, in this review article we will summarise the current knowledge about alcohol-induced epigenetic alterations in the context of ARLD, including but not limited to, DNA hyper/hypo methylation, histone modifications, changes in non-coding RNA, 3D chromatin architecture and enhancer-promoter interactions. Additionally, we will discuss the state-of-the-art techniques used in the study of ARLD (e.g. single-cell sequencing). We will also highlight the epigenetic regulation of chemokines and their proinflammatory role in the context of ARLD. Lastly, we will examine the clinical applications of epigenetics in the diagnosis and management of ARLD.

Different Pigmentation Risk Loci for High-Risk Monosomy 3 and Low-Risk Disomy 3 Uveal Melanomas

BACKGROUND

Uveal melanoma (UM), a rare malignant tumor of the eye, is predominantly observed in populations of European ancestry. UMs carrying a monosomy 3 (M3) frequently relapse mainly in the liver, whereas UMs with disomy 3 (D3) are associated with more favorable outcome. Here, we explored the UM genetic predisposition factors in a large genome-wide association study (GWAS) of 1142 European UM patients and 882 healthy controls .

METHODS

We combined 2 independent datasets (Global Screening Array) with the dataset described in a previously published GWAS in UM (Omni5 array), which were imputed separately and subsequently merged. Patients were stratified according to their chromosome 3 status, and identified UM risk loci were tested for differential association with M3 or D3 subgroups. All statistical tests were 2-sided.

RESULTS

We recapitulated the previously identified risk locus on chromosome 5 on CLPTM1L (rs421284: odds ratio [OR] =1.58, 95% confidence interval [CI] = 1.35 to 1.86; P = 1.98 × 10-8) and identified 2 additional risk loci involved in eye pigmentation: IRF4 locus on chromosome 6 (rs12203592: OR = 1.76, 95% CI = 1.44 to 2.16; P = 3.55 × 10-8) and HERC2 locus on chromosome 15 (rs12913832: OR= 0.57, 95% CI = 0.48 to 0.67; P = 1.88 × 10-11). The IRF4 rs12203592 single-nucleotide polymorphism was found to be exclusively associated with risk for the D3 UM subtype (ORD3 = 2.73, 95% CI = 1.87 to 3.97; P = 1.78 × 10-7), and the HERC2 rs12913832 single-nucleotide polymorphism was exclusively associated with risk for the M3 UM subtype (ORM3 = 2.43, 95% CI = 1.79 to 3.29; P = 1.13 × 10-8). However, the CLPTM1L risk locus was equally statistically significant in both subgroups.

CONCLUSIONS

This work identified 2 additional UM risk loci known for their role in pigmentation. Importantly, we demonstrate that UM tumor biology and metastatic potential are influenced by patients' genetic backgrounds.

Study on the Correlation Between Iris Characteristics and Schizophrenia

BACKGROUND

Recently, researchers have conducted many studies on the potential contribution of the retina and other eye structures on schizophrenia. This study aimed to evaluate differences in iris characteristics between patients with schizophrenia and healthy individuals so as to find more easily accessible and easily measurable biomarkers with a view to improving clinical assessments and furthering our understanding of the disease.

METHODS

Overall, 80 patients with schizophrenia and 52 healthy individuals were included in the case group and the control group, respectively. Iris images were collected from all subjects to compare differences in the structure and color of the iris. The Positive and Negative Symptom Scale (PANSS) and the Modified Overt Aggression Scale (MOAS) were used to evaluate the clinical symptoms and characteristics of 45 first-episode untreated schizophrenics, and analyzed correlations between iris characteristics and schizophrenia symptoms.

RESULTS

There were significant differences in iris crypts (P<0.05) and pigment spots (P<0.01) between the case and control group, but no significant difference was found in iris wrinkles (P<0.05). The logistic regression analysis demonstrated that the total iris crypts [odds ratio (OR) 1.166, 95% confidence interval (CI) 1.022-1.330] and total iris pigment spots (OR 1.815, 95% CI 1.186-2.775) increased the risk of suffering from schizophrenia. Furthermore, it was demonstrated that the number of iris crypts was positively associated with the MOAS score (r=0.474, P<0.01). Moreover, the number of the iris pigment spots (r=0.395, P<0.01) and wrinkles (r=0.309, P<0.05) were positively correlated with the subjects' negative symptom scores, respectively.

CONCLUSION

Iris crypts and pigment spots were identified as potential biomarkers for detecting schizophrenia. In patients with first-episode untreated schizophrenia, iris characteristics may help psychiatrists to identify the illness and its severity, and to detect characteristic clinical symptoms.

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.

Predicting eye and hair colour in a Norwegian population using Verogen's ForenSeq™ DNA signature prep kit

Prediction of eye and hair colour from DNA can be an important investigative tool in forensic cases if conventional DNA profiling fails to match DNA from any known suspects or cannot obtain a hit in a DNA database. The HIrisPlex model for simultaneous eye and hair colour predictions was developed for forensic usage. To genotype a DNA sample, massively parallel sequencing (MPS) has brought new possibilities to the analysis of forensic DNA samples. As part of an in-house validation, this study presents the genotyping and predictive performance of the HIrisPlex SNPs in a Norwegian study population, using Verogen's ForenSeq™ DNA Signature Prep Kit on the MiSeq FGx system and the HIrisPlex webtool. DNA-profiles were successfully typed with DNA input down to 125 pg. In samples with DNA input < 125 pg, false homozygotes were observed with as many as 92 reads. Prediction accuracies in terms of AUC were high for red (0.97) and black (0.93) hair colours, as well as blue (0.85) and brown (0.94) eye colours. The AUCs for blond (0.72) and brown (0.70) hair colour were considerably lower. None of the individuals was predicted to have intermediate eye colour. Therefore, the error rates of the overall eye colour predictions were 37% with no predictive probability threshold (pmax) and 26% with a probability threshold of 0.7. We also observed that more than half of the incorrect predictions were for individuals carrying the rs12913832 GG genotype. For hair colour, 65% of the individuals were correctly predicted when using the highest probability category approach. The main error was observed for individuals with brown hair colour that were predicted to have blond hair. Utilising the prediction guide approach increased the correct predictions to 75%. Assessment of phenotype-genotype associations of eye colours using a quantitative eye colour score (PIE-score), revealed that rs12913832 AA individuals of Norwegian descent had statistically significantly higher PIE-score (less brown eye colour) than individuals of non-northern European descent. To our knowledge, this has not been reported in other studies. Our study suggests that careful assessment of the target population prior to the implementation of forensic DNA phenotyping to case work is beneficial.

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

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

GenNet framework: interpretable deep learning for predicting phenotypes from genetic data

Applying deep learning in population genomics is challenging because of computational issues and lack of interpretable models. Here, we propose GenNet, a novel open-source deep learning framework for predicting phenotypes from genetic variants. In this framework, interpretable and memory-efficient neural network architectures are constructed by embedding biologically knowledge from public databases, resulting in neural networks that contain only biologically plausible connections. We applied the framework to seventeen phenotypes and found well-replicated genes such as HERC2 and OCA2 for hair and eye color, and novel genes such as ZNF773 and PCNT for schizophrenia. Additionally, the framework identified ubiquitin mediated proteolysis, endocrine system and viral infectious diseases as most predictive biological pathways for schizophrenia. GenNet is a freely available, end-to-end deep learning framework that allows researchers to develop and use interpretable neural networks to obtain novel insights into the genetic architecture of complex traits and diseases.

van Hilten and colleagues present GenNet, a deep-learning framework for predicting phenotype from genetic data. This framework generates interpretable neural networks that provide insight into the genetic basis of complex traits and diseases.

Cell-type-specific meQTLs extend melanoma GWAS annotation beyond eQTLs and inform melanocyte gene-regulatory mechanisms

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

Biophysical evidence to support and extend the vitamin D-folate hypothesis as a paradigm for the evolution of human skin pigmentation

OBJECTIVE

To test the "vitamin D-folate hypothesis for the evolution of human skin pigmentation."

METHODS

Total ozone mapping spectrometer (TOMS) satellite data were used to examine surface UV-irradiance in a large (n = 649) Australian cross-sectional study population. Genetic analysis was used to score vitamin D- and folate-related gene polymorphisms (n = 22), along with two pigmentation gene variants (IRF4-rs12203592/HERC2-rs12913832). Red cell folate and vitamin D₃ were measured by immunoassay and HPLC, respectively.

RESULTS

Ultraviolet radiation (UVR) and pigmentation genes interact to modify blood vitamin levels; Light skin IRF4-TT genotype has greatest folate loss while light skin HERC2-GG genotype has greatest vitamin D₃ synthesis (reflected in both TOMS and seasonal data). UV-wavelength exhibits a dose-response relationship in folate loss within light skin IRF4-TT genotype (305 > 310 > 324 > 380 nm). Significant vitamin D₃ photosynthesis only occurs within light skin HERC2-GG genotype, and is maximal at 305 nm. Three dietary antioxidants (vitamins C, E, and β-carotene) interact with UVR and pigmentation genes preventing oxidative loss of labile reduced folate vitamers, with greatest benefit in light skin IRF4-TT subjects. The putative photosensitiser, riboflavin, did not sensitize red cell folate to UVR and actually afforded protection. Four genes (5xSNPs) influenced blood vitamin levels when stratified by pigmentation genotype; MTHFR-rs1801133/rs1801131, TS-rs34489327, CYP24A-rs17216707, and VDR-ApaI-rs7975232. Lightest IRF4-TT/darkest HERC2-AA genotype combination (greatest folate loss/lowest vitamin D₃ synthesis) has 0% occurrence. The opposing, commonest (39%) compound genotype (darkest IRF4-CC/lightest HERC2-GG) permits least folate loss and greatest synthesis of vitamin D₃ .

CONCLUSION

New biophysical evidence supports the vitamin D-folate hypothesis for evolution of skin pigmentation.

Search for genetic loci involved in the constitution and skin type of a Japanese women using a genome-wide association study

The constitution and skin type of individuals are influenced by various factors. Recently, the influence of genetic predispositions on these has been emphasized. To date, genome-wide association studies (GWAS) have shown several single nucleotide polymorphisms (SNPs) that affect individual's constitution and skin type. However, these studies have mainly focused on the Caucasian population, and only a few association analyses with the constitution and skin type of individuals involving a Japanese population have been conducted. In this study, we conducted a GWAS analysis of 9 phenotypes regarding the constitution or skin type of 1108 Japanese women based on a questionnaire. As a result, in addition to SNPs known to be involved in phenotypes in the past, we discovered new SNPs and genetic regions related to darkness of pigmented spots, skin flushing, frequency of rough skin and responsiveness to cosmetics.

Introgression contributes to distribution of structural variations in cattle

Structural variations (SVs) are an important source of phenotypic diversity in cattle. Here, 72 whole genome sequences representing taurine and zebu cattle were used to identify SVs. Applying multiple approaches, 16,738 SVs were identified. A comparison against the Database of Genomic Variants archives revealed that 1575 SVs were novel in our data. A novel duplication covering the entire GALNT15 gene, was observed only in N'Dama. A duplication, which was previously reported only in zebu and associated with navel length, was also observed in N'Dama. Investigation of a novel deletion located upstream of CAST13 gene and identified only in Italian cattle and zebu, revealed its introgressed origin in the former. Overall, our data highlights how the SVs distribution in cattle is also shaped by forces such as demographical differences and gene flow. The cattle SVs of this study and its meta-data can be visualized on an interactive genome browser at https://tinyurl.com/svCowArs.

Prediction of Eye Colour in Scandinavians Using the EyeColour 11 (EC11) SNP Set

Description of a perpetrator's eye colour can be an important investigative lead in a forensic case with no apparent suspects. Herein, we present 11 SNPs (Eye Colour 11-EC11) that are important for eye colour prediction and eye colour prediction models for a two-category reporting system (blue and brown) and a three-category system (blue, intermediate, and brown). The EC11 SNPs were carefully selected from 44 pigmentary variants in seven genes previously found to be associated with eye colours in 757 Europeans (Danes, Swedes, and Italians). Mathematical models using three different reporting systems: a quantitative system (PIE-score), a two-category system (blue and brown), and a three-category system (blue, intermediate, brown) were used to rank the variants. SNPs with a sufficient mean variable importance (above 0.3%) were selected for EC11. Eye colour prediction models using the EC11 SNPs were developed using leave-one-out cross-validation (LOOCV) in an independent data set of 523 Norwegian individuals. Performance of the EC11 models for the two- and three-category system was compared with models based on the IrisPlex SNPs and the most important eye colour locus, rs12913832. We also compared model performances with the IrisPlex online tool (IrisPlex Web). The EC11 eye colour prediction models performed slightly better than the IrisPlex and rs12913832 models in all reporting systems and better than the IrisPlex Web in the three-category system. Three important points to consider prior to the implementation of eye colour prediction in a forensic genetic setting are discussed: (1) the reference population, (2) the SNP set, and (3) the reporting strategy.

Validation of the Verogen ForenSeq™ DNA Signature Prep kit/Primer Mix B for phenotypic and biogeographical ancestry predictions using the Micro MiSeq® Flow Cells

In anticipation of offering phenotypic and biogeographical ancestry predictions to help resolve cases, the Verogen ForenSeq™ DNA Signature Prep kit/Primer Mix B was evaluated in the context of Micro MiSeq® Flow Cells. These flow cells were determined as the best format for a quick turnaround time response and cost effective approach compared to standard flow cells. The phenotype informative SNPs (piSNPs) and ancestry informative SNPs (aiSNPs) were thoroughly examined through sensitivity, reproducibility and repeatability, concordance, robustness (mock casework) and low level DNA mixture studies purposely selecting individuals with different phenotypes (hair and eye color) when possible and different biogeographical ancestry. SNP locus-specific interpretation thresholds were established for the Universal Analysis Software (UAS) based on surviving alleles and SNP predictor rank to minimize false homozygous genotypes and maximize the information that can be derived from an unknown sample. Dropin alleles' intensity determined an appropriate threshold to minimize false heterozygous SNP genotypes. The selection of inappropriate interpretation thresholds was shown to have major consequences on phenotypic predictions. A 3.2% and 4.8% minor DNA component contribution to a DNA mixture had no impact on ancestry predictions whereas a 9.1% contribution did. The multi-locus SNP genotypes generated using the ForenSeq™ DNA Signature Prep kit/Primer Mix B were shown to be reliable, reproducible, concordant and resulted in predictions that were also reliable, reproducible and concordant based on the limited number of donors (N = 19) used in this study.

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

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

Transcriptional co-activator regulates melanocyte differentiation and oncogenesis by integrating cAMP and MAPK/ERK pathways

The cyclic AMP pathway promotes melanocyte differentiation by activating CREB and the cAMP-regulated transcription co-activators 1-3 (CRTC1-3). Differentiation is dysregulated in melanomas, although the contributions of CRTC proteins is unclear. We report a selective differentiation impairment in CRTC3 KO melanocytes and melanoma cells, due to downregulation of oculo-cutaneous albinism II (OCA2) and block of melanosome maturation. CRTC3 stimulates OCA2 expression by binding to CREB on a conserved enhancer, a regulatory site for pigmentation and melanoma risk. CRTC3 is uniquely activated by ERK1/2-mediated phosphorylation at Ser391 and by low levels of cAMP. Phosphorylation at Ser391 is constitutively elevated in human melanoma cells with hyperactivated ERK1/2 signaling; knockout of CRTC3 in this setting impairs anchorage-independent growth, migration, and invasiveness, whereas CRTC3 overexpression supports cell survival in response to the mitogen-activated protein kinase (MAPK) inhibitor vemurafenib. As melanomas expressing gain-of-function mutations in CRTC3 are associated with reduced survival, our results suggest that CRTC3 inhibition may provide therapeutic benefit in this setting.

Evolutionary genetics of skin pigmentation in African populations

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

Genetic Variants in Transcription Factor Binding Sites in Humans: Triggered by Natural Selection and Triggers of Diseases

Variants of transcription factor binding sites (TFBSs) constitute an important part of the human genome. Current evidence demonstrates close links between nucleotides within TFBSs and gene expression. There are multiple pathways through which genomic sequences located in TFBSs regulate gene expression, and recent genome-wide association studies have shown the biological significance of TFBS variation in human phenotypes. However, numerous challenges remain in the study of TFBS polymorphisms. This article aims to cover the current state of understanding as regards the genomic features of TFBSs and TFBS variants; the mechanisms through which TFBS variants regulate gene expression; the approaches to studying the effects of nucleotide changes that create or disrupt TFBSs; the challenges faced in studies of TFBS sequence variations; the effects of natural selection on collections of TFBSs; in addition to the insights gained from the study of TFBS alleles related to gout, its associated comorbidities (increased body mass index, chronic kidney disease, diabetes, dyslipidemia, coronary artery disease, ischemic heart disease, hypertension, hyperuricemia, osteoporosis, and prostate cancer), and the treatment responses of patients.

Rapid speciation via the evolution of pre-mating isolation in the Iberá Seedeater

Behavioral isolation can catalyze speciation and permit the slow accumulation of additional reproductive barriers between co-occurring organisms. We illustrate how this process occurs by examining the genomic and behavioral bases of pre-mating isolation between two bird species (Sporophila hypoxantha and the recently discovered S. iberaensis) that belong to the southern capuchino seedeaters, a recent, rapid radiation characterized by variation in male plumage coloration and song. Although these two species co-occur without obvious ecological barriers to reproduction, we document behaviors indicating species recognition by song and plumage traits and strong assortative mating associated with genomic regions underlying male plumage patterning. Plumage differentiation likely originated through the reassembly of standing genetic variation, indicating how novel sexual signals may quickly arise and maintain species boundaries.

Genome-wide association study in almost 195,000 individuals identifies 50 previously unidentified genetic loci for eye color

Human eye color is highly heritable, but its genetic architecture is not yet fully understood. We report the results of the largest genome-wide association study for eye color to date, involving up to 192,986 European participants from 10 populations. We identify 124 independent associations arising from 61 discrete genomic regions, including 50 previously unidentified. We find evidence for genes involved in melanin pigmentation, but we also find associations with genes involved in iris morphology and structure. Further analyses in 1636 Asian participants from two populations suggest that iris pigmentation variation in Asians is genetically similar to Europeans, albeit with smaller effect sizes. Our findings collectively explain 53.2% (95% confidence interval, 45.4 to 61.0%) of eye color variation using common single-nucleotide polymorphisms. Overall, our study outcomes demonstrate that the genetic complexity of human eye color considerably exceeds previous knowledge and expectations, highlighting eye color as a genetically highly complex human trait.

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

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

Induced pluripotent stem cell reprogramming-associated methylation at the GABRA2 promoter and chr4p12 GABAA subunit gene expression in the context of alcohol use disorder

Twin studies indicate that there is a significant genetic contribution to the risk of developing alcohol use disorder (AUD). With the exception of coding variants in ADH1B and ALDH2, little is known about the molecular effects of AUD-associated loci. We previously reported that the AUD-associated synonymous polymorphism rs279858 within the GABAA α2 receptor subunit gene, GABRA2, was associated with gene expression of the chr4p12 GABAA subunit gene cluster in induced pluripotent stem cell (iPSC)-derived neural cultures. Based on this and other studies that showed changes in GABRA2 DNA methylation associated with schizophrenia and aging, we examined methylation in GABRA2. Specifically, using 69 iPSC lines and neural cultures derived from 47 of them, we examined whether GABRA2 rs279858 genotype predicted methylation levels and whether methylation was related to GABAA receptor subunit gene expression. We found that the GABRA2 CpG island undergoes random stochastic methylation during reprogramming and that methylation is associated with decreased GABRA2 gene expression, an effect that extends to the GABRB1 gene over 600 kb distal to GABRA2. Further, we identified additive effects of GABRA2 CpG methylation and GABRA2 rs279858 genotype on expression of the GABRB1 subunit gene in iPSC-derived neural cultures.

Inbreeding coefficients and runs of homozygosity islands in Brazilian water buffalo

Characterization of autozygosity is relevant to monitor genetic diversity and manage inbreeding levels in breeding programs. Identification of autozygosity hotspots can unravel genomic regions targeted by selection for economically important traits and can help identify candidate genes for selection. In this study, we estimated the inbreeding levels of a Brazilian population of Murrah buffalo undergoing selection for milk production traits, particularly milk yield. We also studied the distribution of runs of homozygosity (ROH) islands and identified putative genes and quantitative trait loci (QTL) under selection. We genotyped 422 Murrah buffalo for 51,611 SNP; 350 of these had ROH longer than 10 Mb, indicating the occurrence of inbreeding in the last 5 generations. The mean length of the ROH per animal was 4.28 ± 1.85 Mb. Inbreeding coefficients were calculated from the genomic relationship matrix, the pedigree, and the ROH, with estimates varying between 0.242 and 0.035. Inbreeding estimates from the pedigree had a low correlation with the genomic estimates, and estimates from the genomic relationship matrix were much higher than those from the pedigree or the ROH. Signatures of selection were identified in 6 genomic regions, located on chromosomes 1, 2, 3, 5, 16, and 18, encompassing a total of 190 genes and 174 QTL. Many of the genes (e.g., APRT and ACSF3) and QTL identified are related to milk production traits, such as milk yield, milk fat yield and percentage, and milk protein yield and percentage. Other genes are associated with reproduction and immune response traits as well as morphological aspects of the buffalo species. Inbreeding levels in this population are still low but are increasing due to selection and should be managed to avoid future losses due to inbreeding depression. The proximity of genes linked to milk production traits with genes associated with reproduction and immune system traits suggests the need to include these latter genes in the breeding program to avoid negatively affecting them due to selection for production traits.

eQTL Colocalization Analyses Identify NTN4 as a Candidate Breast Cancer Risk Gene

Breast cancer genome-wide association studies (GWASs) have identified 150 genomic risk regions containing more than 13,000 credible causal variants (CCVs). The CCVs are predominantly noncoding and enriched in regulatory elements. However, the genes underlying breast cancer risk associations are largely unknown. Here, we used genetic colocalization analysis to identify loci at which gene expression could potentially explain breast cancer risk phenotypes. Using data from the Breast Cancer Association Consortium (BCAC) and quantitative trait loci (QTL) from the Genotype-Tissue Expression (GTEx) project and The Cancer Genome Project (TCGA), we identify shared genetic relationships and reveal novel associations between cancer phenotypes and effector genes. Seventeen genes, including NTN4, were identified as potential mediators of breast cancer risk. For NTN4, we showed the rs61938093 CCV at this region was located within an enhancer element that physically interacts with the NTN4 promoter, and the risk allele reduced NTN4 promoter activity. Furthermore, knockdown of NTN4 in breast cells increased cell proliferation in vitro and tumor growth in vivo. These data provide evidence linking risk-associated variation to genes that may contribute to breast cancer predisposition.

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.

Association between brown eye colour in rs12913832:GG individuals and SNPs in TYR, TYRP1, and SLC24A4

The genotype of a single SNP, rs12913832, is the primary predictor of blue and brown eye colours. The genotypes rs12913832:AA and rs12913832:GA are most often observed in individuals with brown eye colours, whereas rs12913832:GG is most often observed in individuals with blue eye colours. However, approximately 3% of Europeans with the rs12913832:GG genotype have brown eye colours. The purpose of the study presented here was to identify variants that explain brown eye colour formation in individuals with the rs12913832:GG genotype. Genes and regulatory regions surrounding SLC24A4, TYRP1, SLC24A5, IRF4, TYR, and SLC45A2, as well as the upstream region of OCA2 within the HERC2 gene were sequenced in a study comprising 40 individuals with the rs12913832:GG genotype. Of these, 24 individuals were considered to have blue eye colours and 16 individuals were considered to have brown eye colours. We identified 211 variants within the SLC24A4, TYRP1, IRF4, and TYR target regions associated with eye colour. Based on in silico analyses of predicted variant effects we recognized four variants, TYRP1 rs35866166:C, TYRP1 rs62538956:C, SLC24A4 rs1289469:C, and TYR rs1126809:G, to be the most promising candidates for explanation of brown eye colour in individuals with the rs12913832:GG genotype. Of the 16 individuals with brown eye colours, 14 individuals had four alleles, whereas the alleles were rare in the blue eyed individuals. rs35866166, rs62538956, and rs1289469 were for the first time found to be associated with pigmentary traits, whilst rs1126809 was previously found to be associated with pigmentary variation. To improve prediction of eye colours we suggest that future eye colour prediction models should include rs35866166, rs62538956, rs1289469, and rs1126809.

Eye color prediction using the IrisPlex system: a limited pilot study in the Iraqi population

Background

Forensic DNA phenotyping has gained momentum in the recent past due to the prediction of externally visible characters (EVCs) from the biological sample. The most common phenotypes like eye, hair, and skin color are predicted from the biological samples using a web-based system called IrisPlex. Based on six genetic SNPs, the IrisPlex system is developed and validated for its prediction accuracy in diverse ethnic groups worldwide. In previous studies, this system proved to have significant prediction accuracy. The EVCs vary substantially based on different geographical locations. Hence, the objective of this study was to validate the accuracy of the IrisPlex system in predicting the eye colors in the Iraqi population.

Methods

Six genetic single-nucleotide polymorphisms SNPs (HERC2-rs12913832, OCA2- rs1800407, SLC24A4-rs12896399, SLC45A2- rs16891982, TYR-rs1393350, and IRF4- rs12203592) in 58 Iraqi subjects were performed using Sequenom MassARRAY Genotyping. According to Liu et al., a predicted probability of 0.7 was considered as the threshold.

Results

Participants in this study of brown color were observed in 44.83%, intermediate in 43.1%, and blue in 12.07%. Completely predictive accuracy is obtained in 1; we observed the AUC at threshold 0.7 was 0.91 for brown, 0.79 for blue, and 0.60 for intermediate eye color. The sensitivity was 42.85% for blue, 0% for intermediate eye color, and 100% for brown-colored eye. Specificity was 100% for blue, 100% for intermediate, and 78.13% for brown eye color.

Conclusion

Hence, it was concluded that the prediction accuracy of the IrisPlex system for blue and brown color eye in the Iraqi population is significant in the studied population size. However, a pivotal study with larger sample size is required to represent the prediction accuracy of the IrisPlex system in the whole Iraqi population.