Integrative functional genomic analyses identify novel genetic variants influencing skin pigmentation in indigenous Africans [RNA-seq]

Skin color is highly variable in Africans, yet little is known about the underlying molecular mechanism. We identified 1,157 candidate variants influencing skin pigmentation in indigenous Africans by genome-wide association studies and scans of natural selection based on differentiation in allele frequencies between lightly pigmented southern African Khoesan populations and other darkly pigmented African populations. We applied massively parallel reporter and chromosome conformation capture assays to identify novel regulatory variants and their target genes related to skin pigmentation in melanocytic cells. We identified 165 SNPs showing strong differential regulatory activities between alleles. Combining CRISPR-mediated genome editing, transcriptome profiling and melanin assays, we identified causal regulatory variants impacting pigmentation near MFSD12/HMG20B, MITF, OCA2, and DDB1/CYB561A3/TMEM138. We identified CYB561A3 as a novel gene regulating pigmentation by impacting genes involved in oxidative phosphorylation and melanogenesis. Our results broaden our understanding of the genetic basis of human skin color diversity and human adaptation. To test the role of candidate enhancers and variants in skin pigmentation, we performed CRISPR inhibition or knockout of enhancers containing the functional variants identified by MPRA in melanocytic cells. Then, we performed gene expression profiling analysis using data obtained from RNA-seq of these CRISPR-edited cells. We also performed RNA-seq using CYB561A3-koncout MNT1 cells or CYB561A3-overexpressing MNT1 cells Overall design: Comparative gene expression profiling analysis of RNA-seq data for MNT-1 cells and enhancer knockdown or knockout derivatives (CRISPR KO of enhancer E1 of MITF, CRISPR inhibition of enhancer E2 of OCA2, CRISPR KO or inhibition of enhancer E1 of CYB561A3, CYB561A3 KO, CYB561A3 overexpression).